matlab模拟三体运动_matlab三体运动

闲的蛋疼，学学matlab动画，三体运动还是有点问题(提示line没xdate属性，脑瓜也疼了)，不多废话上代码

clc;clear;close;

figure('name','三体运动');

axis equal

grid on

hold on

colordef black

view(3)

xlabel('x')

ylabel('y')

zlabel('z')

hold on

%初始化参数设置(质量、初始位置、初始速度、万有引力常数、时间步长)

m1=rand()*10^23;

m2=rand()*10^23;

m3=rand()*10^23;

pos1=[rand()*10^7,rand()*10^7,rand()*10^7];

pos2=[rand()*10^7,rand()*10^7,rand()*10^7];

pos3=[rand()*10^7,rand()*10^7,rand()*10^7];

v1=[rand()*1000,rand()*1000,rand()*1000];

v2=[rand()*1000,rand()*1000,rand()*1000];

v3=[rand()*1000,rand()*1000,rand()*1000];

G=6.67*10^(-11);

dt=0.005;

%行星和轨迹(定义行星位置和轨迹变量)

planet1=plot3(pos1(1),pos1(2),pos1(3),'b:.','markersize',20);

planet2=plot3(pos2(1),pos2(2),pos2(3),'r:.','markersize',20);

planet3=plot3(po