ecneangle.m

来自「这是用于计算天文中的开普勒轨道元素的程序」· M 代码 · 共 12 行

function E=EcNeAngle(e,M)
Nomin_E(1)=M+e*sin(M)+0.5*e^2*sin(2*M);
Nomin_M(1)=Nomin_E(1)-e*sin(Nomin_E(1));
i=1;
while abs(M-Nomin_M(i))>=0.000000001
    Nomin_E(i+1)=Nomin_E(i)+(M-Nomin_M(i))/(1-e*cos(Nomin_E(i)));
    Nomin_M(i+1)=Nomin_E(i+1)-e*sin(Nomin_E(i+1)); 
    i=i+1;
end
E=Nomin_E(i);
end