星星

[复制链接]

字体大小: 正常放大

李秋锐

14 主题	14 听众	76 积分

升级 74.74%

TA的每日心情

	奋斗 2015-3-29 18:44

签到天数: 17 天

[LV.4]偶尔看看III

电梯直达

1^#

发表于 2015-2-1 18:10 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

shg;R0=1;
a=12*R0;b=9*R0;
T0=2*pi;
T=5*T0;dt=pi/100;t=[0:dt:T]';
f=sqrt(a^2-b^2);
th=12.5*pi/180;
E=exp(-t/20);
x=E.*(a*cos(t)-f);y=E.*(b*cos(th)*sin(t));z=E.*(b*sin(th)*sin(t));
plot3(x,y,z,'g')
[X,Y,Z]=sphere(30);X=R0*X;Y=R0*Y;Z=R0*Z;
grid on,hold on
surf(X,Y,Z)
shading interp
x1=-18*R0;x2=6*R0;y1=-12*R0;y2=12*R0;z1=-6*R0;z2=6*R0;
axis([x1 x2 y1 y2 z1 z2])
view([133 65])
comet3(x,y,z,0.02)
hold off

zan