Matlab动画模拟太阳系行星运动.pdf

Matlab动画模拟太阳系行星运动

figure(name,星系演示);%设置标题名字

pausetime=.02;%设置暂停时间

set(gca,xlim,[-5050],ylim,[-5030],zlim,[-5050]);

set(gcf,doublebuffer,on)%消除抖动

xlabel(x轴),ylabel(y轴),zlabel(z轴);

axisequal;

gridon;

view([352]);

holdon

a=[8.512.5203050608010090];b=[81218264555709030];

omga=[41.2510.50.10.050.250.1251];r=[0.350.80.80.532.51.

51.50.35];%长轴，短轴，角速度，球体半径

c=sqrt(a.^2-b.^2);h=pi/18;h1=pi/10;f=pi/9;g=pi/8;

aby=[hh0;h1h0;hh0;hh0;hh0;hh0;hg0;hhh;g0g];%每个轨道平面倾

斜角度，偏移设置

%colo={y,m,b,m,r,c,b,b};

[X,Y,Z]=sphere(40);

surf(5*X,5*Y,5*Z);colormap(autumn)%设置太阳

light(position,[102],style,infinite)

lightingphong

materialshiny

t=0:0.01*pi:50*pi;

-来源网络，仅供个人学习参考

t;

num=length(a);

forn=1:num

x(:,n)=a(n)*cos(omga(n)*t)+c(n);

y(:,n)=b(n)*sin(omga(n)*t);

z(:,n)=0*t;%计算未经轨道平面角度倾斜的轨道位置

xuanz(:,:)=[100;0cos(aby(n,1))-sin(aby(n,1));0sin(aby(n,

1))cos(aby(n,1))]*[cos(aby(n,2))0sin(aby(n,2));010;-sin(

aby(n,2))0cos(aby(n,2))]*[cos(aby(n,3))-sin(aby(n,3))0;s

in(aby(n,3))cos(aby(n,3))0;001];%每个轨道平面倾斜计算

xyz(:,:)=[x(:,n)y(:,n)z(:,n)]*xuanz(:,:);

x(:,n)=xyz(:,1);

y(:,n)=xyz(:,2);

z(:,n)=xyz(:,3);%计算轨道平面倾斜后的轨道位置

p(n)=surf(r(n)*X+x(1,n),r(n)*Y+y(1,n),r(n)*Z+z(1,n));sha

dinginterp%画出每个行星

plot3(x(:,n),y(:,n),z(:,n),-k);%画出所有轨迹线

end

set(p(1),facecolor,y);

set(p(2),facecolor,m);set(p(3),facecolor,b);set(

p(4),facecolor,m);

set(p(5),facecolor,r);set(p(6),facecolor,c);set(

p(7),facecolor,b);set(p(8),facecolor,b);set(p(9)

-来源网络，仅供个人学习参考

,facecolor,r);%设置所有行星的颜色

form=1:5000%旋转计算

forn=1:length(a)

set(p(n),xdata,r(n)*X+x(m,n),ydata,r(n)*Y+y(m,n),zd

ata,