MATLAB中sin(pi/2-A)与cos(A)精度差异及大圆距离代码疑问

一、为什么用sin(pi/2-A)和cos(pi/2-A)替代？ #

这本质是三角函数诱导公式的应用：sin(A) = cos(pi/2 - A)，cos(A) = sin(pi/2 - A)。代码里把纬度lat转换成了余纬（90°-lat，即从北极点到该点的球面角度），再代入计算，和标准公式完全等价。

这么写核心是为了全纬度范围的数值精度统一：

• 当lat接近90°（北极）或-90°（南极）时，余纬接近0°，此时cos(余纬)对应原公式的sin(lat)，而cos(x)在x接近0时的数值精度极高；反之如果直接计算sin(lat)（lat接近pi/2弧度），虽然精度也不错，但用余纬的写法能避免三角函数在参数接近pi/2时的微小误差波动。
• 对于你关注的lat接近0°的场景，两种写法的精度差异几乎可以忽略，这种写法是为了兼顾高纬度的精度表现而做的统一处理。

二、引入pi会增加误差吗？ #

不会。MATLAB中的pi是IEEE双精度浮点数中最接近真实圆周率的数值，其误差约为1.2246e-16，这个量级的误差远小于绝大多数工程测量数据的误差，完全可以忽略，不会成为计算结果的主要误差来源。

三、MATLAB三角函数的工作原理与精度特性 #

MATLAB的sin、cos、acos等三角函数基于IEEE 754双精度浮点数实现，底层调用经过高度优化的数学库（如Intel MKL），针对不同参数区间选择最优近似算法（泰勒级数、多项式拟合、CORDIC算法等），在精度和效率之间取得平衡。

参数接近0或pi/2时的误差情况： #

• 参数接近0时：
  • sin(x)的精度极高，泰勒展开首项就是x，误差为x³量级，双精度下x=1e-15时，误差约为1e-45，远小于双精度的最小误差阈值（~1e-16）。
  • cos(x)接近1，此时若计算1 - cos(x)会出现浮点抵消误差（两个接近的数相减丢失有效数字），但直接计算cos(x)本身的精度依然很高。
• 参数接近pi/2时：
  • cos(x)接近0，此时cos(x)的误差会略大，因为x的微小误差会直接传递到结果中；而sin(x)接近1，精度极高，类似cos(x)在x接近0时的表现。
  • 用cos(pi/2 - x)代替sin(x)（或反之），本质是把参数转换到精度更高的区间计算，避免在pi/2附近的精度波动。

四、总结 #

你同事的代码写法是一种兼顾全纬度范围精度的优化技巧，通过诱导公式将纬度转换为余纬计算，既保持了和标准公式的等价性，又能在高纬度区域获得更稳定的数值精度。对于你关注的低纬度（lat接近0）场景，两种写法的精度差异极小，而引入的pi带来的误差完全可以忽略。

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