物聯谢鑫实验报告2.docVIP

物联一班 谢鑫 20100740123 实验六：使用m语言对电磁场的仿真 实验目的与要求 掌握m文件调试方法； 掌握利用m语言仿真分析电磁场分布。 二、实验原理及说明 半径为a的环形载流回路周围空间的磁场分布 设载流圆环中流过的电流为I，则圆环在空间任意一点P(x,y,z)产生的磁感应强度矢量为 由于r2=x2+ y2+ z2 则R2=a2+ r2-2×a×r×cosβ = a2+ r2-2×a×ρ = a2+ r2-2×a×(x×cosα+y×sinα) = a2+ x2+ y2+ z2-2a×x×cosα-2a×y×sinα =(x - a×cosα)2+(y - a×sinα)2+ z2 得 由可求出磁感应强度矢量在x,y,z方向的分量分别为 实验内容和步骤 （一）应用MATLAB对半径为a的环形载流回路周围空间的磁场分布进行仿真分析 令B的三个分量中的x=0，仅考虑圆环电流在yoz平面上产生的磁场的分布，而不必考虑Bx分量，则可编写出下面的描绘二维磁场分布的程序： clear; figure(1) a=0.3;y=-1:0.04:1;the=0:pi/20:2*pi; I=1;u0=4*pi*1e-7;K0=I*u0/4/pi; [Y,Z,T]=meshgrid(y,y,the); r=sqrt((a*cos(T)).^2+Z.^2+(Y-a*sin(T)).^2); r3=r.^3; dby=a*Z.*sin(T)./r3; by=K0*trapz(dby,3); dbz=a*(a-Y.*sin(T))./r3; bz=K0*trapz(dbz,3); subplot(121); [bSY,bSZ]=meshgrid([0:0.05:0.2],0); h1=streamline(Y(:,:,1),Z(:,:,1),by,bz,bSY,bSZ,[0.1,1000]); h2=copyobj(h1,gca); rotate(h2,[1,0,0],180,[0,0,0]); h3=copyobj(allchild(gca),gca); rotate(h3,[0,1,0],180,[0,0,0]); title(磁场的二维图,fontsize,15); for kk=1:4 [bSY,bSZ]=meshgrid(0.2+kk*0.02,0); streamline(Y(:,:,1),Z(:,:,1),by,bz,bSY,bSZ,[0.02/(kk+1),4500]); streamline(-Y(:,:,1),Z(:,:,1),-by,bz,-bSY,bSZ,[0.02/(kk+1),4500]); end （三）要求设计编写程序，如果B的三个分量均考虑，编写程序绘出电流环的三维磁力线图。 源代码如下： clear all; figure(1) a=0.3; y=-1:0.04:1; the=0:pi/20:2*pi; I=1; u0=4*pi*1e-7; K0=I*u0/4/pi; [Y,Z,T]=meshgrid(y,y,the); r=sqrt((a*cos(T)).^2+Z.^2+(Y-a*sin(T)).^2); r3=r.^3; dby=a*Z.*sin(T)./r3; by=K0*trapz(dby,3); dbz=a*(a-Y.*sin(T))./r3; bz=K0*trapz(dbz,3); subplot(121); [bSY,bSZ]=meshgrid([0:0.05:0.2],0); h1=streamline(Y(:,:,1),Z(:,:,1),by,bz,bSY,bSZ,[0.1,1000]); h2=copyobj(h1,gca); rotate(h2,[1,0,0],180,[0,0,0]); h3=copyobj(allchild(gca),gca); rotate(h3,[0,1,0],180,[0,0,0]); title(磁场的二维图,fontsize,15); for kk=1:4 [bSY,bSZ]=meshgrid(0.2+kk*0.02,0); streamline(Y(:,:,1),Z(:,:,1),by,bz,bSY,bSZ,[0.02/(kk+1),4500]); streamline(-Y(:,:,1),Z(:,:,1),-by,bz,-bSY,bSZ,[0.02/(kk+1),4500]); end [X,Y,Z]=meshgrid(-0.5:0.04:0.5); r2=X.^2+Y.^2+Z.^