[MATLAB基础] 极坐标系下画图

[云龙九现]

上次在论坛里游逛碰到一个问题，如何在极坐标系下画图的问题
我回答了这个问题，现在将解决方式整理一下，以附件形式上传，希望能互相学习

[wwpwwp1985]

十分感谢，正在寻找。。

[280448499]

xiexie!!!!!!!!!!

[sandylamb]

谢谢楼主！很有用的程序！

[jingjie11]

谢谢楼主

[jwping1986]

请问下楼主，要是我想画半个圆呢，从负九十度到正九十度，也就是右半个圆，我画的是一个半圆的温度图，温度值求出来后但是不知道怎么转换过来，谢谢啊，最近很为这个问题纠结

[jwping1986]

请问下楼主，要是我想画半个圆呢，从负九十度到正九十度，也就是右半个圆，我画的是一个半圆的温度图，温度值求出来后但是不知道怎么转换过来，谢谢啊，最近很为这个问题纠结

[jwping1986]

请高手帮忙，如何使我的图能画出半圆环的效果
我之前的程序如下：
dX = (pi/2-(-pi/2))/10
dY = (Ro-Ri)/10
X(1) = - pi/2;
X(2)=dX/2+X(1);
for i = 3:ni-1
X(i)=X(i-1)+dX;
end

Y(1) = Ri;
Y(2)=dY/2+Y(1);
for j = 3:nj-1
Y(j)=Y(j-1)+dY;
end;
Y(nj)=Ro;


for i = ni-1:-1:2
T(i,j) = P(i)*T(i+1,j)+Q(i);

[jwping1986]

请高手帮忙，如何使我的图能画出半圆环的效果
我之前的程序如下：
dX = (pi/2-(-pi/2))/10
dY = (Ro-Ri)/10
X(1) = - pi/2;
X(2)=dX/2+X(1);
for i = 3:ni-1
X(i)=X(i-1)+dX;
end

Y(1) = Ri;
Y(2)=dY/2+Y(1);
for j = 3:nj-1
Y(j)=Y(j-1)+dY;
end;
Y(nj)=Ro;


for i = ni-1:-1:2
T(i,j) = P(i)*T(i+1,j)+Q(i);
end;
其中P,Q为X的函数。
polar(X,Y,T);shading interp;
但是画出来的图，不是二维直角上的平面，就是一个圆环，完全不是半圆环，很是郁闷，请大家帮帮我，不甚感激

[云龙九现]

theta = -pi/2 : pi/180: pi/2;
r1 = 1.5;
r2 = 1.2;
x1 = r1 * cos(theta);
y1 = r1 * sin(theta);
x2 = r2 * cos(theta);
y2 = r2 * sin(theta);
figure(1)
hold on
grid on
plot(x1,y1,'b.')
plot(x2,y2,'r.')
axis equal

我不知道上面的源程序是不是解决了你的问题。你的意图我不是完全明白