49 磁场计算.docVIP

练习49 磁场计算 知识背景 磁场是一个很基本的电磁场现象。如同电场计算一样，磁场计算在科学研究和工程实际问题中有着广泛的应用。我们在这个练习中着手解决磁场问题。并试图用图形将数据可视化，从而使我们清楚地把握磁场特征。主要内容【本练习讲述知识点】 本练习考查读者综合使用编程、绘图、逻辑验证等来解决实际磁学问题的能力。我们将利用linspace语句、for循环语句、subplot和mesh绘图命令及逻辑运算符。练习中涉及到较为复杂的程序，希望读者仔细体会。练习过程 （1）电流环产生的磁场 我们来结合实际例子看一下如何解决这类问题： 我们来看看如何用毕奥-萨伐定律计算电流环产生的磁场。磁学知识告诉我们，载流导线产生的磁场规律为： 任一电流元在空间任一点P处产生的磁感应强度为：其中，为电流元到P点的矢径，l为导线圆的长度矢量。则P点的总磁场可沿载流导体全长积分各段产生的磁场来求得。我们在命令区里输入：R=2.5;I0=4;s=4*pi*1e-7;C0=I0*s/(4*pi);x=linspace(-3,3,20);y=x;N=20;t0=linspace(0,2*pi,N+1);t1=t0(1:N);y1=R*cos(t1);z1=R*sin(t1);t2=t0(2:N+1);y2=R*cos(t2);z2=R*sin(t2);dlx=0;dly=y2-y1;dlz=z2-z1;xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;for i=1:20 for j=1:20;rx=x(j)-xc;ry=y(i)-yc;rz=0-zc;r3=sqrt(rx.^2+ry.^2+rz.^2).^3;dlXrx=dly.*rz-dlz.*ry;dlXry=dlz.*rx-dlx.*rz;Bx(i,j)=sum(C0*dlXrx./r3);By(i,j)=sum(C0*dlXry./r3);endendclf;quiver(x,y,Bx,By)得到的结果如图49-1所示。从图中，我们清楚地看出整个电流环的分布情况。如果我们改变参数，比如将R改为5，看看磁场分布图，如图49-2所示。图49-1 电流环产生的磁场分布图49-2 增大R后的磁场分布从图中可以看到，增大R后，磁场变化很大。（2）赫姆霍兹线圈的验证所谓赫姆霍兹线圈，就是间距正好等于线圈半径的一对相同的共轴载流线圈。赫姆霍兹线圈沿轴线附近的磁场的大小十分均匀，而且都沿x轴方向。我们来对这个结论进行验证。我们在命令区里输入相应的程序（见下页），运行后，得到图49-3。R=1;I0=5;s=4*pi*1e-7;C0=I0*s/(4*pi);NGx=21;NGy=21;x=linspace(-R,R,NGx);y=linspace(-R,R,NGy);N=20;t0=linspace(0,2*pi,N+1);t1=t0(1:N);y1=R*cos(t1);z1=R*sin(t1);t2=t0(2:N+1);y2=R*cos(t2);z2=R*sin(t2);dlx=0;dly=y2-y1;dlz=z2-z1;xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;for i=1:NGy for j=1:NGx;rx=x(j)-xc;ry=y(i)-yc;rz=0-zc;r3=sqrt(rx.^2+ry.^2+rz.^2).^3;dlXrx=dly.*rz-dlz.*ry;dlXry=dlz.*rx-dlx.*rz;Bx(i,j)=sum(C0*dlXrx./r3);By(i,j)=sum(C0*dlXry./r3);endendBax=Bx(:,11:21)+Bx(:,1:11);Bay=By(:,11:21)+By(:,1:11);subplot(1,2,1),mesh(x(11:21),y,Bax);xlabel(x);ylabel(y);subplot(1,2,2);plot(y,Bax),grid,xlabel(y);ylabel(Bx)本例的计算模型与第一个例子差不多，但是将观测范围取在两线圈之间。B生成的线圈左边的磁场等于A线圈的左边磁场，所以A、B两线圈在中间部分的合成磁场等于A线圈的右磁场与左磁场平移R后的和。因此，我们需要观测A线圈的左右区间x=[-R,R]内的磁场。由图49-3我们可以得到这样的结论：在赫姆霍兹曲线的两个线圈之间的轴线附近，有一个很大的区域，x方向的磁场强度Bx比较均匀。类似地，该区域内的y方向的磁场强度By近似为零。再来看一下参数设置，如果我们将参数略作改动，即取R=5，I0=1。运行程序，我们将得到图49-4。图49-3 赫