公司建设网站首页/怎样做推广是免费的

%线圈轴线与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+(l-z).^2); %单位m,点P到导线的距离

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

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

%BC边设为2产生的磁场

B2=n*C0*(d2.^-1).*((s21.^-1).*(a+x)+(s22.^-1).*(a-x));

%BC边设为2在XYZ的分量

B2x=0;

B2y=B2.*(d2.^-1).*(l-z);   %定义B2y为Y正向

B2z=B2.*(d2.^-1).*(a-y);

%CD边设为3

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

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

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

%CD边设为3产生的磁场

B3=n*C0*(d3.^-1).*((s31.^-1).*(a+y)+(s32.^-1).*(a-y));

%CD边设为3在XYZ的分量

B3x=B3.*(d3.^-1).*(l-z)*(-1);  %B3x取反向

B3y=0;

B3z=B3.*(d3.^-1).*(a+x);

%DA边设为4

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

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

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

%DA边设为4产生的磁场

B4=n*C0*(d4.^-1).*((s41.^-1).*(a+x)+(s42.^-1).*(a-x));

%AB边设为1在XYZ的分量

B4x=0;

B4y=B4.*(d4.^-1).*(l-z)*(-1);  %B4y取反向

B4z=B4.*(d4.^-1).*(a+y);

%AB边设为1

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

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

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

%AB边设为1产生的磁场

B5=n*C0*(d5.^-1).*((s51.^-1).*(a+y)+(s52.^-1).*(a-y));

%AB边设为1在XYZ的分量

B5x=B5.*(d5.^-1).*(l+z)*(-1);   %B5x取反向

B5y=0;

B5z=B5.*(d5.^-1).*(a-x);

%BC边设为2

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

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

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

%BC边设为2产生的磁场

B6=n*C0*(d6.^-1).*((s61.^-1).*(a+x)+(s62.^-1).*(a-x));

%BC边设为2在XYZ的分量

B6x=0;

B6y=B6.*(d6.^-1).*(l+z)*(-1);  %B6y取反向

B6z=B6.*(d6.^-1).*(a-y);

%CD边设为3

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

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

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

%CD边设为3产生的磁场

B7=n*C0*(d7.^-1).*((s71.^-1).*(a+y)+(s72.^-1).*(a-y));

%CD边设为3在XYZ的分量

B7x=B7.*(d7.^-1).*(l+z);

B7y=0;

B7z=B7.*(d7.^-1).*(a+x);

%DA边设为4

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

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

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

%DA边设为4产生的磁场

B8=n*C0*(d8.^-1).*((s81.^-1).*(a+x)+(s82.^-1).*(a-x));

%AB边设为1在XYZ的分量

B8x=0;

B8y=B8.*(d8.^-1).*(l+z);

B8z=B8.*(d8.^-1).*(a+y);

Bx=B1x+B3x+B5x+B7x;

By=B2y+B4y+B6y+B8y;

Bz1=B1z+B2z+B3z+B4z;

Bz2=B5z+B6z+B7z+B8z;

Bz=B1z+B2z+B3z+B4z+B5z+B6z+B7z+B8z;

B=B1+B2+B3+B4+B5+B6+B7+B8;

subplot(3,1,1);

plot(z,Bz1,'*')

axis([-1.00 1.00 -0.5e-4 1.5e-4])

%axis([-0.200 0.200 0 1e-4])

hold on;

subplot(3,1,1);

plot(z,Bz2,'*')

axis([-1.00 1.00 -0.5e-4 1.5e-4])

%axis([-0.200 0.200 0 1e-4])

subplot(3,1,2);

plot(z,Bz,'.')

axis([-1.00 1.00 0 1.5e-4])

%axis([-0.200 0.200 0 1e-4])

subplot(3,2,5);

plot(x,Bx,'*')

axis([-1.00 1.00 -0.5e-4 0.5e-4])

%axis([-0.200 0.200 0 1e-4])

subplot(3,2,6);

plot(y,By,'*')

axis([-1.00 1.00 -0.5e-4 0.5e-4])

%axis([-0.200 0.200 0 1e-4])

这是一个磁场分布的程序，有做过类似的可以提些建议

我想把它修改成三维图形显示的，如B,x(或y),z的三维图，希望高手给修改下

非常感谢