res3body.m

来自「大气层内航天器飞行动力学matlab仿真模型」· M 代码 · 共 11 行

% Copyright Ashish Tewari (c) 2006
function Xdot=res3body(t,X)
m=0.00095369; 
x=X(1); 
y=X(2); 
r1=sqrt((x-m)^2+y^2); 
r2=sqrt((x+1-m)^2+y^2); 
Xdot(1,1)=X(3); 
Xdot(2,1)=X(4);
Xdot(3,1)=x+2*Xdot(2,1)-(1-m)*(x-m)/r1^3-m*(x+1-m)/r2^3;
Xdot(4,1)=y-2*Xdot(1,1)-(1-m)*y/r1^3-m*y/r2^3;