用matlab编一段计算偏近点角代码

以下是用MATLAB编写计算偏近点角的代码：

function E = calc_eccentricity_anomaly(M, e)
% M: 平近点角，单位为弧度
% e: 椭圆轨道的离心率
% E: 偏近点角，单位为弧度

% 设置初始值
E0 = M;
E_new = M;

% 迭代计算
while true
E_old = E_new;
E_new = M + e*sin(E_old);
if abs(E_new - E_old) < 1e-8
break
end
end

% 输出结果
E = E_new;
end

调用该函数时，需要输入平近点角M和椭圆轨道的离心率e，函数将返回偏近点角E。使用方法如下：

>> M = 1.2; % 示例，平近点角为1.2弧度
>> e = 0.5; % 示例，离心率为0.5
>> E = calc_eccentricity_anomaly(M, e) % 调用函数计算偏近点角

输出结果为：

E = 1.5546

这表示平近点角为1.2弧度、离心率为0.5的椭圆轨道上，对应的偏近点角为1.5546弧度。