MATLAB在电磁学求解方面的应用mmm.docVIP

MATLAB在求解电磁学中方程和实验数据处理方面的 应用 王敏娣 （渭南师范学院物理与电气工程学院 物理学 2011级专升本班） 摘要：本文针对求解电磁学中一些计算比较复杂或利用常规方法无法给出解析解的问题，而通过Matlab编程给出其数值解。该程序简短,结果精确,可避免烦琐、复杂的数学运算。 这些电磁学问题包括: 微分方程及超越方程、实验数据处理等。 关键字：Matlab；电磁学； 微分方程；超越方程; 实验数据 1引言 MATLAB是一种具有强大的数值分析、矩阵计算、高科技计算程序。其功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。在教学中引进计算机数值方法，有助于将一些高深的物理知识深入浅出、生动形象地讲清楚，而且有利于激发学生的学习兴趣，从而提高学习积极性，使问题变的直观、清楚、简洁明了。 本文主要通过对几个具体电磁学问题的处理来体现Matlab在处理计算方面所表达出来的优越性及方便简洁性，从而让大家对Matlab在高等物理计算机解题方面的应用更加清晰、明了。以便在以后类似的解题中更好的运用Matlab这个强大的数学计算工具帮助我们解决那些复杂、繁琐、深奥的物理问题和物理现象。 2 MATLAB求解电磁学中的微积分问题、超越方程的举例与实验数据处理分析 2.1求解微分方程问题 在一B=0.2T，方向垂直于轨道面向下的均匀磁场中，有一长L=1M，质量M=0.1KG的金属杆，沿一倾a=30金属滑道由静止下滑，如图2.1所示（1）若滑道与自感系数L=0.5H的线圈相连，试编写一计算机程序，考察该金属杆的运动速度及线圈电流随时间的变化关系；如果改变自感L的大小，金属感的运动速度及线圈内的电流将如何变化？ 在上个问题中如果将与滑轨相连的自感线圈改为一个的电阻，其他条件不 变，试编写一个计算机程序，考察该金属杆的运动速度及线圈内电流随时间的变化关系；如果改变电阻的R的大小，金属杆的运动速度及线圈内电流的关系将如何变化？ 解（1）：当金属杆下滑时，它在磁场中运动而产生感应电动势并有感应电流通过，其运动方程和回路方程分别是 (1-1) （1-2） （ 图2.1） 利用MATLAB 求解 , 程序2.1（1a）如下: （其中， clear global K1 K2 v=dsolve(D2v+K1*v+K2=0,v(0)=0,Dv(0)=0); %求解微分方程 v=simplify(v) %简化结果表达式 运行程序得： v =K2*(cos(K1^(1/2)*t)-1)/K1 因此，金属杆的运行速度为 又由（1—2）式得： 执行程序2.1（1a），可以得到金属杆的运动速度与时间的关系，如图【1a】所示 图【1a】金属杆的运动速度与时间的曲线图 线圈中电流随时间的变化程序2.1(1b)如下： clear K1=0.01; K2=-0.5; dt=0.1; t=[0:dt:200]; I=BlK2/LK1*[(sin(K1^(1/2)*t)/K1^(1/2))-t] plot(I,t); xlabel(I/A);ylabel(t/s) 执行程序2.1（1b)后，可以得到线圈中电流随时间的变化关系曲线图【1b】 【1b】线圈内电流随时间的变化关系 小结：由上述两条曲线不难看出，无论金属杆的下滑速度还是线圈内的电流均随时间做振荡。事实上，如果将式（1-1）对时间t再一次求导，并代入式（1-2），有 （1-3） 很显然，这是速度的谐振动方程，谐振动的周期为 （1-4） 由此可见，振荡周期随L的增大而有所增长。 解（2）：当金属杆下滑时，它在磁场中运动而产生感应电动势并有感应电流通过，其运动方程和电回路方程分别是