MATLAB的答案

让一个矢量图在极坐标下的速度场

70(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;
[θ,ρ]= 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中央,发现社区如何帮助你!

开始狩猎!