apolloeq.m

来自「科学运算语言 MATLAB 的程序设计有关的函数代码」· M 代码 · 共 12 行

function dx=apolloeq(t,x)
%APOLLOEQ 函数用来描述 Apollo 卫星的数学模型，它可以直接用于
%  常微分方程的数值求解。

%Designed by Prof D Xue (c) 2000
mu=1/82.45; mu1=1-mu; 
r1=sqrt((x(1)+mu)^2+x(3)^2); r2=sqrt((x(1)-mu1)^2+x(3)^2);
dx=[x(2); 
    2*x(4)+x(1)-mu1*(x(1)+mu)/r1^3-mu*(x(1)-mu1)/r2^3;
    x(4); 
    -2*x(2)+x(3)-mu1*x(3)/r1^3-mu*x(3)/r2^3];