MATLAB在电磁场教学中的应用.doc

MATLAB在电磁场课程中的应用摘要：电磁场课程理论性强，概念抽象，需要较强的多维空间想象能力和逻辑思维能力，不能直观的进行观察和研究，难以很好地掌握。文中简要介绍了MATLAB语言的基本计算功能和画图功能，并对电磁场课程中的具体实例进行了理论计算及可视化仿真，这样不仅提高了计算速度，而且也进一步加深了对电磁场空间物理现象的理解。关键词：电磁场；MATLAB；可视化1 引言电磁场理论是分析各种电磁现象的基本规律、应用原理与应用方法的技术基础课，是培养合格的电气信息类专业本科生所应具备的知识结构的重要组成部分。公共基础课（数学、物理等）侧重于抽象问题的分析与计算，而专业课又侧重于工程实际中的应用，电磁场则起到了承前启后的作用，使学生们初步认识各种电磁现象及电磁过程的物理本质。掌握运用多种数学工具解决电磁问题的方法和技巧，为学生顺利进入专业课的学习打下坚实的基础1。电磁场涉及内容较广，概念抽象，是空间与时间综合性最强的课程之一。应用的很多内容在数学的教学中往往不是重点内容，可在电磁场的教学中，这些内容又是分析电磁现象的重要数学工具。可见，对数学基础薄弱的学生来说，“教”与“学”都感到非常困难。针对这种情况传统的教学模式已经逐渐不能适应时代的发展的要求，因此在教学中积极采用现代化设备，通过高科技手段使学生能够直接获取知识，成为自身学习及各个高校教学的热点。而MATLAB具有强大的计算及绘图能力，在电磁场教学中应用非常广泛。2 MATLAB特点及应用MATLAB是由美国MathWorks公司推出的一款优秀的程序仿真开发软件。经过多年的逐步发展与不断完善，已经成为国际公认的最优科学计算与数学应用软件之一。其内容涉及矩阵代数、微积分、应用数学、计算机图形学、物理等很多方面。集计算、绘图及声音处理于一体，主要特点如下2,3：（1）计算功能强大。能够实现数值与符号计算、计算结果与编程可视化、数字与文字的统一处理、离线与在线计算等，针对不同领域提供了丰富的工具箱，用户还可以根据自己的需要任意扩充函数工具库。（2）强大的绘图功能。能够实现二维、三维图形的绘制，可以从图形直观的衡量程序的效果。（3）界面友好。效率高，编程简洁，MATLAB以矩阵为基本单元的可视化程序设计语言，指令表达和标准教材的数学表达式相近。（4）简单易学，特别适用于初学者，用户可以在短时间内掌握。正是由于MATLAB强大的功能和广泛的适用性，才得到了用户的普遍认可，在自动控制、神经网络、信号处理等诸多方面，都有广泛的应用。3 应用MATLAB实现电磁场计算电磁场涉及数学知识很多，如：积分变换、矢量分析、场论等，也涉及到泛函分析、变分法、微分几何、积分方程等方面的基础知识，在函数分析中变量是三维空间，甚至是在四维空间中讨论电磁场的变化，变化量既有标量又有矢量。这是电磁场课程不容易掌握的一个主要原因。而MATLAB几乎可以解决科学计算的任何问题。应用举例一：设单芯电缆有两层绝缘体,分界面亦是同轴圆柱面,电缆上电荷体密度=0.6c/,内层绝缘体介电常数为2，外层绝缘体介电常数为3.8，内导体绝缘体半径为1cm，内层绝缘半径为3cm,外层绝缘体半径为7cm,求内导体与外壳导体之间的电压U为多少？解：在绝缘体中取任意点P，设P至O点的距离为p。过P点作同轴圆柱面，高为，该面再加上下两底面作为“高斯面S”。由于对称，显然D在上下底面上没有法相分量，在同轴圆柱面上D是均匀的且沿半径向外取向。应用高斯定律得： (1)于是各层绝缘体中电场强度分别为 和 (2)而电压为 (3)程序如下：m1=2;m2=3.8;t=0.6;p1=1;p2=3;p3=7;% m为介电常数，t为线密度，p为半径E1=(p)t./(2*pi*m1*p); %求内层绝缘体场强E2=(p)t./(2*pi*m2*p); %求外层绝缘体场强U3=quad(E1,p1,p2); %求内层绝缘体电压U4=quad(E2,p2,p3); %求内层绝缘体电压U=U3+U4; %求内外导体之间电压程序运行结果：U = 0.07374 应用MATLAB实现电磁场可视化在电磁场场量分析中抽象思维程度很高，电磁场不同于一般物质的五态，没有固定形态、没有静止质量、没有颜色，甚至没有明确的大小边界，很不容易直接感知，这也是电磁场课程不容易掌握的另一个主要原因。但如果采用 MATLAB 计算并绘图，将电力线、等位线等用二维或三维图形清晰展现出来，学生的理解会更加直观4。应用举例二：平行双线传输线可看做两根单位带电量分别是+和-的无限长细圆柱或直线，试画出其电位分布。 解：已知线密度为均匀分布的无限长线电荷周围的电场为由于线电荷无限长，零参考电位点不能取在无穷远点，一般可任意指定某一位置为零参考点，因此，单根线电荷的电位场为：平行双线的电位场是两根单线的场的叠加：以上求解过程往往容易理解，但是由唯一性定理，若要使两平行线电荷在两圆柱导体外部空间引起的电场与两圆柱导体之间原来的电场完全相同，则要找到两个与两圆柱导体表面圆周相重合的圆周来，换句话来说，圆柱导体表面是等位面，若电轴产生的电位使原来圆柱导体所在位置表面电位相同即可。这一层次学生往往难以理解，现在通过MATLAB编程画出平行双线的等位线图，就可以清楚地看到等位面所在的位置。程序如下：%平行双线的等位线图如图1所示x,y=meshgrid(-3:.01:3,-3:.01:3);f=log(sqrt(x+1).2+y.2+eps)-log(sqrt(x-1).2+y.2+eps);v=-17,-1.5,-1,-.5,-.2,0,.2,.5,1,1.5,17;C,h=contour(x,y,f,v,m);clabel(C)xlabel(x)ylabel(y)图1 平行双线的等位线图%平行双线的电位和归一化电场分布如图2所示x,y=meshgrid(-3:.25:3,-3:.25:3);f=log(sqrt(x+1).2+y.2+eps)-log(sqrt(x-1).2+y.2+eps);v=-17,-1.5,-1,-.5,-.2,0,.2,.5,1,1.5,17;C,h=contour(x,y,f,v,b);hold ondx,dy=gradient(-f,.25,.25);D=sqrt(dx.2+dy.2);dx=dx./D;dy=dy./D;quiver(x,y,dx,dy,.7);xlabel(x)ylabel(y) 图2 平行双线的电位和归一化电场分布%平行双线的电位三维立体图如图3所示syms x yV=log(sqrt(x+1).2+y.2)-log(sqrt(x-1).2+y.2);xMax=8;NGrid=40;xPlot=linspace(-xMax,xMax,NGrid);x,y=meshgrid(xPlot);VPlot=eval(V);ExPlot,EyPlot=gradient(-VPlot);clf;subplot(1,2,1),meshc(VPlot);xlabel(x);ylabel(y);zlabel(电位);图3 平行双线的电位三维立体图通过引入MATLAB强大的绘图功能，可以将数据以多种图形形式表现出来，实现了电磁场可视化，使电磁场中的概念更加直观、清晰，易于接受，使学生能够进一步深入分析、理解电磁场的各种性能。4 结语在电磁场课程教学的过程中，利用MATLAB软件进行技算、模拟、实现结果的可视化，大大提高了学生的解题速度，有效地提高了学生学习的兴趣，使学生能够进一步理解电磁场的空间物理现象，同时也丰富了教师教学的方法和手段，为电磁场理论的可视化提供了一个新的平台。参考文献1 冯慈璋,马西奎.工程电磁场导论.高等教育出版社.20002 周立鹏等. MATLAB在电磁场教学中的应用Z. 科技信息, 2009, 35：516-5173 刘美丽.MATLAB语言及应用.国防工业出版社.20124 杰荣,蔡新华,胡惟文.基于MATLAB的空间电磁场分布可视化研究.中国科技论文统计源期刊.2005 The Application of MATLAB in Electromagnetic FieldSHI Lei, HAO Jing(Northeast Dianli University, Jilin Jilin 132012 )Abstract: Electromagnetic Field theory is hard and the concepts are nonrepresentational, requiring us to have strong imaginary abilities of multidimensional space and logical thinking abilities as well. In this paper the basic calculation and painting functions of MATLAB language are introduced, and particular examples of the Electromagnetic Field are calculated and visualization processed， thu