求助:有谁做了磁偶极子的磁场的matlab模拟,学习一下

给你贴个电偶极袜毁子场的MATLAB模拟吧！clearclfq=2e-6k=9e9a=1.5b=-1.5x=-6:0.6:6y=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./AEEy=Ey./AE% 场强归一化，使箭头等长cv=linspace(min(V(:)),max(V(:)),49)% 产生 49 个电孝春位值contourf(X,Y,V,cv,'k-') % 用黑实线画告慎备填色等位线图title('\fontname{隶书}偶极子的场','fontsize',20),hold onquiver(X,Y,Ex,Ey,0.7) % 第五输入宗量 0.7 使场强箭头长短适中。plot(a,b,'wo',a,b,'w+') % 用白线画正电荷位置plot(-a,-b,'wo',-a,-b,'w-') % 用白线画负电荷位置xlabel('x')ylabel('y'),hold off

x=-2:0.5:2y=-2:0.5:2

[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 on

a=0b=0

plot(a,b,'ro'或拿蠢昌,a,b,'r+')

---