zhen_pinng_pian_jindiaojiao.m

上一个的继续,连续有九个压缩文件,大家可以打开解压后自己慢慢享用,这个都是本人调试通过了的,希望共同进步!

%真近点角的计算Fs：
%Ms：平近点角    Es：偏近点角
DToR=2*pi/360;  %转换常数从角度到弧度
Ms=0:1:359; %平近点序列（角度）
es=0.01;    %地球偏心率
Ms=DToR*Ms;
%用迭代法求偏近点角
%Es=Ms+es*sin(Es);  开普勒方程，迭代计算公式
Es=zeros(1,360);    %生成偏近点角序列的存储变量
for i=1:360
    dEs=1;          %赋初值
    TempMs=Ms(i);
    TempEs=TempMs;
    while (abs(dEs)>0.0001)     %计算精度控制循环的循环次数
        TempEsNext=TempMs+es*sin(TempEs);
        dEs=TempEs-TempEsNext;
        TempEs=TempEsNext;
    end
    Es(i)=TempEs;   %存储计算所得的偏近点角
end
%Es=Es/DTor;
%Ms=Ms/DTor;
Fs=zeros(1,360);    %生成真近点角序列的存储变量，已知偏近点角求真近点，
v1=sqrt((1+es)/(1-es)); %公式中使用的常数，注意避免将此计算过程放入循环！
for i=1:360
    v2=v1*tan(Es(i)/2);
    Fs(i)=mod((2*atan(v2)+2*pi),2*pi);
end
Es=Es/DToR; %将平近点转换为角度
Ms=Ms/DToR; %将偏近点转换为角度
Fs=Fs/DToR; %将真近点转换为角度