工程电磁场实验报告_电磁场实验报告
工程电磁场实验报告由刀豆文库小编整理,希望给你工作、学习、生活带来方便,猜你可能喜欢“电磁场实验报告”。
工程电磁场实验报告
一.题目
有一极长的方形金属槽,边宽1m,除顶盖电位为100sinπx V外,其他三面的电位均为零,试用差分法求槽内电位的分布。
二.原理
如下图所示,用分别平行于x,y轴的两组直线把场域D划分为许多正方形网格,网格线交点称为节点,两相邻平行网格线间的距离h称为步距
用表示节点处电位值,利用二元函数泰勒公式,与节点(Xi,Yj)直接相邻的节点上的电位表示为
整理可得差分方程
这就是二维拉普拉斯方程的差分格式,它将场域内任意一点的位函数值表示为周围直接相邻的四个位函数值的平均值。这一关系式对场域内的每一节点都成立,也就是说,对场域的每一个节点都可以列出一个上式形式的差分方程,所有节点的差分方程构成联立差分方程组。
已知的边界条件经离散化后成为边界点上已知数值。若场域的边界正好落在网格点上,则将这些点赋予边界上的位函数值。一般情况下,场域的边界不一定正好落在网格节点上,最简单的近似处理就是将最靠近边界点的节点作为边界节点,并将位函数的边界值赋予这些节点。
如何计算:简单迭代法
先对静电场内的节点赋予迭代初值,其上标(0)表示初始近似值。然后再按 下面的公式:
进行多次迭代(k=0,1,2,3…)。当两次邻近的迭代值差足够小时,就认为得到了电位函数的近似数值解。如何计算:超松弛迭代法
三.编程序 bc=50;%网格数
u=zeros(bc+1,bc+1);%步长为1/bc %********附初值********* w=0;
for j=1:bc+1;
u(1,j)=100*sin((j-1)*pi/bc);w=w+u(1,j);end for i=2:bc
for j=2:bc u(i,j)=w./bc;
end end
%*************************************************** h=input('please input h(1
for i=2:bc;
for j=2:bc;a=u(i,j);b=u(i,j+1);c=u(i+1,j);d=u(i-1,j);e=u(i,j-1);f=(b+c+d+e)/4;u(i,j)=a+h.*(f-a);
end
end end
%***************绘图******************** x=0:1/bc*1:1;y=0:1/bc*1:1;[x,y]=meshgrid(x,y);mesh(x,y,u)四.结果
五.实验体会
这学期上了这门课,刚开始上课,老师说教我们手算,不用计算机,终于体会到了学了计算方法的好处。熊爷爷教了线性方程组的解法,非线性方程组的解法,雅克比迭代法,解微分方程组,解积分方程组的欧拉法以及均匀媒质中的有限差分法。认真听完了课,让我认识到学计算方法的用处,又再次了解掌握了计算法的原理和运用。在自己写实验报告的时候,又再次去了解了差分法计算电磁场,并成功的编写程序将其计算出来。总的来说,不仅锻炼脑,还锻炼了手!
电磁场仿真实验报告电气工程学院 2011级2班 2011302540056 黄涛实验题目:有一极长的方形金属槽,边宽为1m,除顶盖电位为100sin (pi*x)V外,其它三面的电位均为零,试用差分法求槽内......
《工程电磁场》学习心得班级: 姓名: 学号:在开始学习“工程电磁场”之前,当我听到其学科名称的时候就产生了一种高深莫测的感觉,觉得电磁场应该是比较难的。但是出于对知识的渴望......
个人总结工程电磁场计算是电气专业的公共必修课程,对于我们电气专业的研究生而言,其重要意义不言而喻。今年的下学期在由邹玲老师教授的这门课程中,通过老师细心的讲解和独具一......
华北电力 北京 考研 工程电磁场 教学大纲.txt大人物的悲哀在于他们需要不停地做出选择;而小人物的悲哀在于他们从来没有选择的机会。男人因沧桑而成熟,女人因成熟而沧桑。男人......
电 磁 场 论 文 072202H王焱200722070223 电子高新技术与电磁场理论摘 要 本文就最近发展的高新技术中有关电磁场和电磁波问题展开探讨,并在此基础上对当前高新技术的发展与......
