类似这样的用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))) % 路程
最简单的是 你可以先用初始值解微分方程,然后再。。。。
请问数学方法怎么求解啊q
您好,我是有问必答小助手,您的问题已经有小伙伴帮您解答,感谢您对有问必答的支持与关注!