如何用Matlab數(shù)值模擬磁場?
網(wǎng)友解答: 我在此提供兩種磁鐵/導(dǎo)線/電偶極子磁場在各處分布的Matlab源程序:單個磁鐵/直導(dǎo)線磁場分布源程序: x=-20:0.5:20;y=-20:0.5:20;[X,Y]=mes
我在此提供兩種磁鐵/導(dǎo)線/電偶極子磁場在各處分布的Matlab源程序:
單個磁鐵/直導(dǎo)線磁場分布源程序:
x=-20:0.5:20;y=-20:0.5:20;[X,Y]=meshgrid(x,y);BX=-Y./((X).^2+Y.^2);BY=X./((X).^2+Y.^2);quiver(X,Y,BX,BY,0.8)title('單個磁鐵磁場分布'),hold ona=0;b=0;plot(a,b,'ro',a,b,'r+') ;圖1 單個磁鐵/直導(dǎo)線磁場分布圖
圖2 單個磁鐵/直導(dǎo)線磁場分布中心放大圖
電偶極子電磁場分布clear;clf;q=2e-6;k=9e9;a=2.0;b=0;x=-10:0.6:10;y=x;[X,Y]=meshgrid(x,y);rp=sqrt((X-a).^2+(Y-b).^2);rm=sqrt((X+a).^2+(Y+b).^2);V=q*k*(1./rp-1./rm);[Ex,Ey]=gradient(-V);AE=sqrt(Ex.^2+Ey.^2);Ex=Ex./AE;Ey=Ey./AE;cv=linspace(min(min(V)),max(max(V)),51);contour(X,Y,V,cv,'r-')%axis('square')title('fontname{ 宋體 }fontsize{11} 電偶極子的電場線與等勢線 '),hold onquiver(X,Y,Ex,Ey,0.6,'b')plot(a,b,'bo',a,b,'b+')plot(-a,-b,'bo',-a,-b,'w-')xlabel('x');ylabel('y'),hold off圖3 電偶極子電場線與等勢線
圖4 電偶極子電場線與等勢線中心放大圖