关于#matlab#的问题：第一次使用matlab

Question

想解决以上问题。关于公式的、坐标点的输入和答案于图案的结果不正确。希望能得到各位的帮助，谢谢！

Answer 1

类似这样的用ode45求解

clc; clear
mu = 4902.8^(3/2);
eqs = @(t, x)[
    x(2)
    - mu/(x(1)^2+x(3)^2+x(5)^2)^(3/2)*x(1);
    x(4);
    - mu/(x(1)^2+x(3)^2+x(5)^2)^(3/2)*x(3);
    x(6);
    - mu/(x(1)^2+x(3)^2+x(5)^2)^(3/2)*x(5);
    ];
x0 = [1440.4, -0.77897, -510.34, 0.27599, 854.12, 1.4787];
[t, x] = ode45(eqs, [0, 7000], x0);
plot3(x(:,1), x(:,3), x(:,5), 'r-')
ds = diff(x(:,1:2:6));
S = sum(sqrt(sum(ds.^2,2))) % 路程