Текст программы для расчета составляющих магнитной индукции в MatLab

figure

R = 0.01; L = 0.01; Ms = 868000; mu0=4*pi*1e-7;

for z = 0:0.001:0.05

Bz0 = dblquad (@(Z,f) (mu0/(4*pi)).*(Ms*R^2)./((z-Z).^2+R^2).^(3/2),-L,0,0,2*pi);

hold on

plot(z,Bz0,'k.')

grid on

title('Bz(z) for magnet')

xlabel('z,m')

ylabel('Bz(z), Tl')

end

figure

z1 = 0;

for r=0:0.001:0.05

Bz1 = dblquad(@(Z,f) (mu0/(4*pi)).*(Ms.*R.*(r-R))./((z1-Z).^2+R^2+r.^2-...

2.*R.*r.*cos(f)).^(3/2),-L,0,0,2*pi);

hold on

plot(r,Bz1,'k.')

grid on

title('Bz(z,r) for magnet')

xlabel('r,m')

ylabel('Bz(r), Tl')

end

figure

for z = 0:0.001:0.3

Br0 = dblquad(@(Z,f) (mu0/(4*pi)).*(Ms.*R.*Z)./((z-Z).^2+R^2+r.^2-...

2.*R.*r.*cos(f)).^(3/2),-L,0,0,2*pi);

hold on

plot(z,Br0,'k.')

grid on

title('Br(z) for magnet')

xlabel('z,m')

ylabel('Br(z), Tl')

end

figure

for r = 0:0.001:0.05

Br1 = dblquad(@(Z,f) (mu0/(4*pi)).*(Ms.*R.*(z1-Z))./((z1-Z).^2+R^2+r.^2-...

2.*R.*r.*cos(f)).^(3/2),-L,0,0,2*pi);

hold on

plot(r,Br1,'k.')

grid on

title('Br(z,r) for magnet')

xlabel('r,m')

ylabel('Br(r), Tl')

end