3D axis rotation with meshgrid

2 次查看(过去 30 天)
Laura Mugica
Laura Mugica 2018-9-17
Hello, I would like to rotate B around the z axis and plot the contour plot. I tried the following but I don't think it is working since the contour plot after the rotation should not be symmetric.
Hope anyone can help,
thet=[0:0.01*pi:pi];
plane=[-15:0.02:15,15:0.02:15];
z = 10;
[theta,plane,z]=meshgrid(thet,plane,z);
[x,y,z]=pol2cart(theta,plane,z);
r = sqrt(x.^2+y.^2+z.^2)*1e-9;
B = const*((300e-18)./r.^5-1./r.^3)*1e6;
phi = 53;
%rotation around Z axis
X = x.*cos(phi) - y.*sin(phi);
Y = x.*sin(phi) + y.*cos(phi);
contourf(x,y,B,25,'LineStyle','none')

回答(1 个)

Naman Chaturvedi
Naman Chaturvedi 2018-9-20
编辑:Naman Chaturvedi 2018-9-20
You aren't plotting the updated values X and Y. Last line should be:-
>> contourf(X,Y,B,25,'LineStyle','none')
  1 个评论
Laura Mugica
Laura Mugica 2018-9-20
I get exactly the same plot when I use contourf(X,Y,B,25) that's why I think it's not working.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Surface and Mesh Plots 的更多信息

产品


版本

R2017a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by