TLAB及其在理工课程中的应用指南 第三版 十二五陈怀琛 第5－9章 第6章.pptVIP

　　◆ MATLAB程序 　clear all; 　%初始化(给定环半径、电流、图形) 　mu0=4*pi*1e-7; %真空导磁率(T*m/A) 　I0=5.0; Rh=1; %这两个常数不影响结果 　C0=mu0/(4*pi) * I0; 　 %归并常数 　%下面三行输入语句与上题不同, 观测范围x取［-Rh, Rh］, 即线圈的左右都取, %因为以后要把A线圈右边的磁场和B线圈左边的磁场相加 NGx=1; NGy = 21; %设定观测点网格数 x=linspace(-Rh, Rh, NGx); %设定观测点范围及数组 y=linspace(-Rh, Rh, NGy); %y也取［-Rh, Rh］ 主程序段同例6-5-1(从Nh=20到最后一个end) Bax=Bx(∶,11∶21)+Bx(∶, 1∶11); %x0 区域的磁场平移叠加到x0 区域 Bay=By(∶, 11∶21)+By(∶, 1∶11); subplot(1, 2, 1), mesh(x(11∶21),y,Bax);xlabel(′x′); ylabel(′y′); %画出其Bx分布三维图 subplot(1, 2, 2), plot(y, Bax), grid, xlabel(′y′); ylabel(′Bx′); 　　◆ 运行结果 　　运行结果如图 6-12 所示。可以看出, 在亥姆霍兹线圈的两个线圈之间的轴线附近相当大的一个区域内, x方向的磁场强度Bx是非常均匀的。用相仿的语句也不难得知, 该区域内的y方向的磁场强度By近似为零，这留给读者自行证明。为了看出这个区域的大小和形状, 可以用多种方法, 这里用最简单的一种非图形方法，键入 图 6-12 亥姆霍兹线圈轴线附近Bx按x、y的网格曲面和沿y向的分布图 (a) Bx按x、y的网格曲面; (b) Bx沿y向的分布图 ［q］=abs((Bax-Bax(11, 6))/Bax(11, 6))0.02 %磁场相对误差＜2%的点 q = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 0 0 1 0 0 0 1 0 0 0 0 0 0 1 1 1 1 1 0 0 0 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 0 0 0 1 1 1 1 1 1 1 0 0 →x轴线 0 1 1 1 1 1 1 1 1 1 0 0 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 1 1 1 1 1 0 0 0 0 0 0 1 0 0 0