clear;

a=30

w=2*180/25.4/24/3600

Vsw=3.4*10^5

r0=7*10^8

r=linspace(r0,10^12,102)

phi=linspace(0,360,361)

[r,phi]=meshgrid(r,phi)

theta=90+asind(sind(a)*sind(phi+w*(r-r0)/Vsw))

x=r.*sind(theta).*cosd(phi);

y=r.*sind(theta).*sind(phi);

z=r.*cosd(theta);

surf(x,y,z)

title('行星际电流片三维成图');

xlabel('x');

ylabel('y');

zlabel('z')

hold on

编程代码2

clear;

a=pi/12

w=2.7*10^(-6)

Vsw=340000

r0=7*10^8

r=linspace(r0,10^12,100)

phi=linspace(0,2*pi,100)

[r,phi]=meshgrid(r,phi)

theta=pi/2+asin(sin(a)*sin(phi+w*(r-r0)/Vsw))

[x,y,z]=sph2cart(phi,pi/2-theta,r)

surf(x,y,z)

hold on