用matlab画极坐标方程,matlab极坐标

把polar(b,p)移到end上面，然后再在end上面加一句hold on试试

更多追问追答

追问

b=[];

p=[];

L=180;

R=50;

for i=1:1000

y=integral(@(t)sqrt(1+(pi*pi/25)*cos(pi*t/1000)^2),0,i);

b=0.0014379531169*y;

r=-((1+(pi*pi/25)*cos(pi*i/1000)^2)^(1.5))/((pi*pi*pi*i/12500)*sin(pi*i/1000));

a=asin(L/r);

s=b*tan(a);

p=R+s;

polar(b,p)

hold on

end

好像还是不行，我对MATLAB完全不懂，只是需要画一个图

追答

b=[];

p=[];

L=180;

R=50;

for i=1:1000

y=integral(@(t)sqrt(1+(pi*pi/25)*cos(pi*t/1000).^2),0,i);

b(i)=0.0014379531169*y;

r=-((1+(pi*pi/25)*cos(pi*i/1000)^2)^(1.5))/((pi*pi*pi*i/12500)*sin(pi*i/1000));

a=asin(L/r);

s=b(i)*tan(a);

p(i)=R+s;

end

polar(b,p)

b=[];

p=[];

L=180;

R=50;

for i=1:1000

y=integral(@(t)sqrt(1+(pi*pi/25)*cos(pi*t/1000).^2),0,i);

b(i)=0.0014379531169*y;

r=-((1+(pi*pi/25)*cos(pi*i/1000)^2)^(1.5))/((pi*pi*pi*i/12500)*sin(pi*i/1000));

a=asin(L/r);

s=b(i)*tan(a);

p(i)=R+s;

end

polar(b,p)