MATLAB航天器方向追踪：远距飞行时的精度损失问题

针对航天器远离太阳后，计算地心位置矢量r_A = r - r_Earth时因大数相减导致的精度损失问题，可采用以下几种实用方案：

• 改用相对位置递推计算
  由于航天器无外力作用沿直线运动，其日心速度v_sc恒定；若地球的日心速度v_Earth在验证时段内可近似为恒定，可通过初始时刻的地心位置r_A0 = r0 - r_Earth0，直接递推后续时刻的地心位置：

  r_A(t) = r_A0 + (v_sc - v_Earth) * (t - t0)
  

  这种方式完全避免了每次计算大数r与r_Earth的减法操作，从根源上消除精度损失。

• 使用高精度数值类型
  若必须保留r - r_Earth的计算逻辑，可将原有数值类型替换为更高精度的格式：

  • 编程语言中，C++可使用long double，Python可使用decimal模块（需设置足够高的精度位数）；
  • 注意：高精度计算会带来一定的性能开销，若仅用于方向验证，需权衡精度与效率的平衡。

• 直接从方向角度推导验证逻辑
  无需计算完整的r_A矢量，通过数学变换直接验证方向一致性：
  已知目标地心单位矢量dir，需验证(r - r_Earth)/|r - r_Earth| ≈ dir。利用矢量运算变换：

  1. 计算点积项：dot_r_dir = r · dir，dot_earth_dir = r_Earth · dir，则(r - r_Earth) · dir = dot_r_dir - dot_earth_dir；
  2. 计算r - r_Earth的模长近似值（当|r| >> |r_Earth|时）：

     |r - r_Earth| ≈ |r| * sqrt(1 - 2*(r·r_Earth)/|r|² + |r_Earth|²/|r|²)
     

  3. 方向误差的余弦值为(dot_r_dir - dot_earth_dir) / |r - r_Earth|，通过该值与1的偏差即可计算角度误差，全程避免大数相减操作。

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