gobetter的个人博客分享 http://blog.sciencenet.cn/u/gobetter

博文

日球层电流片matlab成图

已有 3889 次阅读 2015-5-9 19:18 |个人分类:编程|系统分类:科研笔记| 日球层电流片成图

                            假设夹角α为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









https://wap.sciencenet.cn/blog-2292236-888823.html

上一篇:parker太阳风解-MATLAB实现
收藏 IP: 124.16.100.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-11-1 07:26

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部