matlab绘制磁场图,Matlab的三维绘图，磁场的空间分布

%线圈轴线与Z轴重合，正方形边长2a,P(x,y,z)

%将X，Y方向的直接定义正方向

clear all;

u0= 4* pi*1e-7 ; %初始化, 给定半径、电流

n=28;             %匝数

I0=3.0;

a=0.75;   %边长2a ，a=750

l=0.5445*a ; %最佳间距

C0=u0 /( 4* pi)*I0;

Nx=201; Ny=201; Nz=201;           % 设定观测点网格数

x=linspace(-1,1,Nx);

y=linspace(-1,1,Ny);

z=linspace(-1,1,Nz);

%AB边设为1

d1=sqrt((a-x).^2+(l-z).^2); %单位m,点P到导线的距离

s11=sqrt((a-x).^2+(a+y).^2+(l-z).^2);%(a+y)

s12=sqrt((a-x).^2+(a-y).^2+(l-z).^2);%(a-y)

%AB边设为1产生的磁场

B1=n*C0*(d1.^-1).*((s11.^-1).*(a+y)+(s12.^-1).*(a-y));

%AB边设为1在XYZ的分量

B1x=B1.*(d1.^-1).*(l-z);    %定义B1x为X正向

B1y=0;

B1z=B1.*(d1.^-1).*(a-x);

%BC边设为2

d2=sqrt((a-y).^2+(