万有引力太阳系行星轨迹模拟

一、运行结果

二、MATLAB程序 

%太阳及水金地火四行星的运行轨迹模拟
%Author: Lei Zhen
%Date: 2019-04-26

%万有引力公式F=GMm/r^2
%物体运动是速度与加速度产生的位置变化
clc
clear all
close all
set(0,'defaultfigurecolor','w');

%T  质量        x坐标       y坐标  z坐标   x速度  y速度      z速度 x加速度 y加速度 z加速度 星球直径
%                                        （合速度，不是分量速度）
T=[
    1.9891e30   0           0       0       0   0           0   0       0       0        140e7           %太阳1
    3.30e23     56.6720e9   0       0       0   47.88e3     0   0       0       0        4880000         %水星2 
    4.869e24    108.2065e9  0       0       0   35.02e3     0   0       0       0        12103600        %金星3 
    5.98e24     149.57