极坐标系下画图

[云龙九现]

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

[wwpwwp1985]

十分感谢，正在寻找。。

[280448499]

xiexie!!!!!!!!!!

[sandylamb]

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

[jingjie11]

谢谢楼主:):):)

[jwping1986]

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

[jwping1986]

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

[jwping1986]

请高手帮忙，如何使我的图能画出半圆环的效果 http://eg5a39.mail.163.com/js3/main.jsp?sid=qAXIqkDDrKbEAPoAyWDDUENnWLjygsxD
我之前的程序如下：
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

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

[jwping1986]

谢谢你啊，总算有人对我的问题回答了，而且让我看到梦寐以求的圆环。可以加我QQ吗287377075，帮我看看这程序，非常感谢:)，我的是网格划分，求温度场T的。

[jwping1986]

谢谢，总算有人回答，你可以加我QQ吗，287377075，帮我看看这个程序，谢谢你啊，谢谢:)

[jwping1986]

谁能帮我啊，这个极坐标图能画吗

[未注册]

谢谢，先下下来看看。

[CodeAnt]

很基础的说明，在此谢过了。
顺便对作者那句 axis squire 表示敬意。

[horsehh]

非常感谢，很有用，希望有问题能相互讨论