|||
假设夹角α为30°,电流片成图如下:
编程代码1 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 |
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-1 07:26
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社