MATLAB的答案

在极坐标下画出速度场的矢量图

40次浏览(过去30天)
玛丽娜Markaki
玛丽娜Markaki 2021年5月7日
评论道: 玛丽娜Markaki2021年5月8日
在画出压力的等高线(我很容易就画好了)之后，我被要求画出速度场的矢量图。压力为 P = (c * (R1、R2) - (1 - c)) *日志((x ^ 2 + y ^ 2) ^(1/2)) /日志(R2 / R1) + 1 - c 而且 P = (c * (R1、R2) - (1 - c)) *日志(r) /日志(R2 / R1) + 1 - c 在极坐标中，c R1 R2是常数。我可以用这个代码画出直角坐标部分的向量图但我不知道如何在极坐标下做。谢谢你！
c = 0.1;
R1 = 1;
R2 = 10;
X = 1:10;
Y = 1:10;
[X,Y] = meshgrid(X,Y);
P = (c * (R1、R2) - (1 - c)) *日志((x ^ 2 + y ^ 2) ^(1/2)) /日志(R2 / R1) + 1 - c;
p_x = (c * (R1、R2) - (1 - c)) * (X / (X ^ 2 + y ^ 2)) /日志(R2 / R1);
p_y = (c * (R1、R2) - (1 - c)) * (Y / (x ^ 2 + Y ^ 2)) /日志(R2 / R1);
图;
箭袋(X, Y, p_x p_y)
标题(“速度场图”）

登录回答这个问题。

答案(1)

乍得格林
乍得格林 2021年5月7日
你能用 cart2pol 把坐标和矢量分量转换成极坐标?
1评论
玛丽娜Markaki
玛丽娜Markaki 2021年5月8日
这是我写的代码但它给出了U和V必须是相同大小的错误。
c = 0.1;
R1 = 1;
R2 = 10;
X = -10:10;
Y = -10:10;
[theta,rho] = cart2pol(x,y);
P = (c * (R1、R2) - (1 - c)) *日志(ρ)/日志(R2 / R1) + 1 - c;
p_r = (c * (R1、R2) - (1 - c))。/(ρ*日志(R2 / R1));
p_theta = 0;
图;
箭袋(θ,ρ,p_r p_theta)
标题(“速度场图”）

登录评论。

类别

了解更多极地的情节在帮助中心而且文件交换

标签

社区寻宝

在MATLAB Central中找到宝藏，并发现社区如何帮助您!

开始狩猎!