keplerhyp.m

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

% Copyright Ashish Tewari (c) 2006
function H=keplerhyp(e,M);
H=M;
fH=e*sinh(H)-H-M;
fpH=e*cosh(H)-1;
dH=-fH/fpH;
H=H+dH;
eps=1e-10;
i=0;
while abs(fH)>eps
i=i+1  
fH=e*sinh(H)-H-M
fpH=e*cosh(H)-1;
dH=-fH/fpH;
H=H+dH
end