螺绕环截面绘制

在做电磁场的模拟随便弄的

clc,clear,close all
% n=input('螺绕环的匝数n:');
% i=input('电流大小I:');
% r1=input('内径r1:');
% r2=input('外径r2:');
n=20;
i=1;
r1=2;
r2=3;
r=(r1+r2)/2;
m=linspace(r1,r2);
k=4e-7;
b=k*i./(2*m);
t=linspace(0,2*pi,200);
x1=r1*cos(t);x2=r2*cos(t);
y1=r1*sin(t);y2=r2*sin(t);
%螺绕环绘制
plot(x1,y1,x2,y2) 
axis equal
hold on
for a=1:n
    x=r1*cos(2*pi/n*a);
    y=r1*sin(2*pi/n*a);
    plot(x,y,'o',MarkerSize=20,MarkerEdgeColor=[0.5,0.5,0.5])
    plot(x,y,'.',MarkerSize=15,LineWidth=2,MarkerEdgeColor=[0.9,0.5,0.5])
    x=r2*cos(2*pi/n*a);
    y=r2*sin(2*pi/n*a);
    plot(x,y,'o',MarkerSize=20,MarkerEdgeColor=[0.5,0.5,0.5])
    plot(x,y,'x',MarkerSize=15,LineWidth=2,MarkerEdgeColor=[0.5,0.5,0.9])
end
ylim([-r2-n/10,r2+n/10])
plot(r*cos(t),r*sin(t),LineWidth=4)
plot(r,0,'^',MarkerFaceColor='k',MarkerSize=16)
title('螺绕环演示',FontSize=16)
hold off
% 磁场变化
figure
plot(m,b)
title('环内磁场变化')
xlabel('r(r1 to r2)')
ylabel('磁感应强度B')
grid on