基于Jean Meeus天文算法的角度归算结果不符问题求助

问题根源 #

你当前的角度归算逻辑错误处理了负数角度：当角度为负时，你先取绝对值再计算模360的余数，得到的是正数角度的余数，但Jean Meeus书中要求的是将负数角度归算到0-360度区间，需要计算的是负数本身的模运算结果（或等价补值）。

比如：

• 对于-8234.2625，你得到Abs(-8234.2625) mod 360 = 314.2625，但正确结果应为360 - 314.2625 = 45.7375
• 对于-108984.6278，你得到Abs(-108984.6278) mod 360 = 264.6278，正确结果应为360 - 264.6278 = 95.3722

正确的归算逻辑 #

要将任意正/负角度归算到0-360度区间，核心是基于**向下取整（Floor）**的计算，公式为：

reducedDeg := deg - 360 * Floor(deg / 360);

验证示例： #

1. 计算-8234.2625：

   • deg / 360 = -8234.2625 / 360 ≈ -22.87295
   • Floor(-22.87295) = -23（向下取整，取比当前值小的最大整数）
   • reducedDeg = -8234.2625 - 360*(-23) = -8234.2625 + 8280 = 45.7375

2. 计算-108984.6278：

   • deg / 360 = -108984.6278 / 360 ≈ -302.73508
   • Floor(-302.73508) = -303
   • reducedDeg = -108984.6278 - 360*(-303) = -108984.6278 + 109080 = 95.3722

修正你的代码实现 #

替换当前的Abs(val mod x)逻辑，改用基于Floor的通用函数：

// 将角度归算到0-360度区间的函数
function ReduceAngle(deg: Double): Double;
var
  quotient: Double;
begin
  quotient := deg / 360.0;
  quotient := Floor(quotient);
  Result := deg - 360.0 * quotient;
  // 处理浮点数精度误差，确保结果严格在0-360之间
  if Result < 0 then Result := Result + 360.0;
  if Result >= 360 then Result := Result - 360.0;
end;

使用示例： #

// M的计算代码
deg := -8234.2625;
reducedDeg := ReduceAngle(deg); // 得到45.7375
convertedDeg := reducedDeg * Pi / 180; // 得到0.7982弧度

// M1的计算代码
deg := -108984.6278;
reducedDeg := ReduceAngle(deg); // 得到95.3722
convertedDeg := reducedDeg * Pi / 180; // 得到1.6646弧度

内容的提问来源于stack exchange，提问作者Akira Takeshi