2025年太阳地球月亮运动轨迹MATLAB仿真程序.pdfVIP

博观而约取，厚积而薄发。——苏轼

太阳地球月亮运动轨迹MATLAB仿真程序

太阳、地球、月亮运动轨迹MATLAB仿真程序

该代码包含了地球绕太阳运动，月亮绕地球运动的MATLAB轨迹仿真程序，实

时显示地球、月亮的运动轨迹。有兴趣的朋友可以购买，以供交流。程序从第二页

开始。

M文件1:draw_ball.m

%画三维球体的函数

%(x0,y0,z0)为球心

%r为球半径

functiondraw_ball(x0,y0,z0,r)[x1,y1,z1]=sphere;

x=x1*r+x0;

y=y1*r+y0;

z=z1*r+z0;

surf(x,y,z);

M文件2:draw_circle.m

%画二维圆形

%(x0,y0)为圆心

%r为圆半径

functiondraw_circle(x0,y0,r)theta=0:pi/100:2*pi;

x=r*cos(theta)+x0;

y=r*sin(theta)+y0;

plot(x,y,-r);

M文件3:RungeKutta_EarthSun.m%四阶Runge-kutta法解地日微分方程

百学须先立志。——朱熹

functionyh=RungeKutta_EarthSun(w,h)

oldy=w;

k1=dery(oldy);

midy=oldy+k1*h/2;

k2=dery(midy);

midy=oldy+k2*h/2;

k3=dery(midy);

midy=oldy+k3*h;

k4=dery(midy);

yh=oldy+(k1+2*k2+2*k3+k4)*h/6;

%地日微分方程组

functiondy=dery(w)

G=6.674e-11;Ms=1.989e30;miu_s=G*Ms;r2=w(4)^2+w(5)^2;

dy(1)=1;

dy(2)=(-1)*miu_s*w(4)/r2^1.5;dy(3)=(-1)*miu_s*w(5)/r2^1.5;dy(4)

=w(2);

dy(5)=w(3);

M文件4:RungeKutta_MoonEarth.m%四阶Runge-kutta法解地月微分方程

functionyh=RungeKutta_MoonEarth(w,h)

oldy=w;

k1=dery(oldy);

midy=oldy+k1*h/2;

k2=dery(midy);

midy=oldy+k2*h/2;

好学近乎知，力行近乎仁，知耻近乎勇。——《中庸》

k3=dery(midy);

midy=oldy+k3*h;