球面余弦定理 Spherical Law of Cosines
深入解析球面余弦定理：从几何原理到现代导航的核心算法
什么是球面余弦定理？
在平面几何中，我们知道三角形三边与夹角的关系由余弦定理描述：。然而，当我们将三角形绘制在球面上时（即球面三角形），由于曲率的存在，平面几何的公式不再直接适用。球面余弦定理（Spherical Law of Cosines）正是描述球面上球面三角形的边（以角度度量）与角之间关系的核心公式。
它是解决球面三角学问题的基础工具之一，广泛应用于地理信息系统（GIS）、航空航天、航海导航以及天文学等领域。对于任何试图理解全球定位系统（GPS）如何计算距离的用户来说，掌握球面余弦定理是理解底层逻辑的关键一步。
核心概念：球面三角形
由三个大圆（Great Circles）的交点构成的三角形。其边长不是直线距离，而是圆心角的角度值（弧度或度）。
关键变量
a, b, c: 球面三角形的边长（角度）。
A, B, C: 对应的内角。
适用范围
适用于任何半径为 R 的球体。在地球近似为球体的假设下，它是计算地表两点间最短路径（大圆距离）的理论基础。
公式详解与推导逻辑
对于边长为 a, b, c 和对应角为 A, B, C 的球面三角形，球面余弦定理有两种主要形式：用于计算边和用于计算角。
1. 边的余弦定理（最常用）
这是导航计算中最重要的形式，用于已知两边及其夹角，求第三边。
| 公式名称 | 数学表达式 | 物理意义 |
|---|---|---|
| 边余弦定理 | cos(c) = cos(a)cos(b) + sin(a)sin(b)cos(C) |
已知两边 a, b 及夹角 C，求第三边 c |
| 角余弦定理 | cos(C) = -cos(A)cos(B) + sin(A)sin(B)cos(c) |
已知两角 A, B 及夹边 c，求第三角 C |
2. 在地球导航中的具体应用形式
在地理导航中，我们通常已知两个地点的纬度（Lat）和经度（Lon）。我们可以将地球视为一个单位球，极点为 P。设地点1为 A，地点2为 B。
- 边 a = 90° - Lat1 （余纬）
- 边 b = 90° - Lat2 （余纬）
- 夹角 C = |Lon1 - Lon2| （经度差）
- 边 c = 中心角 θ （两点间的角距离）
代入公式后，经过三角恒等变换（利用 sin(90°-x)=cos(x), cos(90°-x)=sin(x)），我们得到导航中常用的球面余弦定理变体：
cos(θ) = sin(Lat1)sin(Lat2) + cos(Lat1)cos(Lat2)cos(ΔLon)
其中，θ 是两点相对于地心的夹角（弧度），ΔLon 是经度之差。一旦求出 θ，距离 D 即为：D = R × θ（R 为地球半径）。
天文学中的延伸应用
除了地球导航，球面余弦定理在天文学中同样扮演着不可或缺的角色。天文学家将星空视为一个巨大的“天球”，所有天体都投影在这个球面上。
1. 计算天体角距离
假设我们要计算两颗恒星之间的角距离。已知它们的赤经（RA）和赤纬（Dec），这与地球上的经度和纬度类似。直接应用球面余弦定理即可得出两星在天球上的夹角。
2. 地平坐标系转换
在观测天体时，我们需要将赤道坐标系（RA/Dec）转换为地平坐标系（方位角 Azimuth / 高度角 Altitude）。这涉及到观测者所在地的纬度、天体的赤纬以及本地时角（LHA）。这些转换公式本质上都是球面三角学的变体，其中余弦定理是核心。
太阳高度角计算
用于太阳能板的角度调整。通过公式：sin(h) = sin(φ)sin(δ) + cos(φ)cos(δ)cos(H)
其中 h 为高度角，φ 为纬度，δ 为太阳赤纬，H 为时角。这正是边余弦定理的变形。
月相与位置预测
精确预测月食或日食时，需要计算太阳、地球、月球三者在天球上的相对位置，球面三角学提供了必要的数学框架。
与平面余弦定理及Haversine公式的对比
用户常问：球面余弦定理与其他距离计算公式有何不同？
平面余弦定理
公式： c² = a² + b² - 2ab cos(C)
适用场景： 小范围区域（如一个城市内的几公里），此时地球曲率可忽略不计，表面近似平面。
缺点： 在全球尺度下误差极大，完全不可用。
Haversine公式
公式： a = sin²(Δlat/2) + cos(lat1)cos(lat2)sin²(Δlon/2)
适用场景： 通用球面距离计算，特别是短距离。
优点： 数值稳定性好，避免了 arccos(1) 附近的精度损失问题。
缺点： 计算步骤稍多，涉及更多三角函数运算。
球面余弦定理
公式： cos(c) = cos(a)cos(b) + sin(a)sin(b)cos(C)
适用场景： 长距离计算，几何意义直观。
缺点： 当两点非常接近时，cos(c) 接近 1，arccos 的微小误差会被放大，导致距离计算不准确。
编程实现示例
以下是使用 Python 实现球面余弦定理计算距离的代码。代码中包含了详细的注释，便于开发者理解。
import math
def spherical_law_of_cosines(lat1, lon1, lat2, lon2):
"""
使用球面余弦定理计算两点间的大圆距离
:param lat1, lon1: 起点纬度、经度（十进制度）
:param lat2, lon2: 终点纬度、经度（十进制度）
:return: 距离（公里）
"""
# 地球平均半径 (km)
R = 6371.0
# 转换为弧度
lat1_rad = math.radians(lat1)
lon1_rad = math.radians(lon1)
lat2_rad = math.radians(lat2)
lon2_rad = math.radians(lon2)
# 经度差
dlon = lon2_rad - lon1_rad
# 应用球面余弦定理
# cos(c) = sin(lat1)sin(lat2) + cos(lat1)cos(lat2)cos(dlon)
cos_c = (math.sin(lat1_rad) math.sin(lat2_rad) +
math.cos(lat1_rad) math.cos(lat2_rad) math.cos(dlon))
# 防止浮点数误差导致 cos_c 超出 [-1, 1] 范围
cos_c = max(-1, min(1, cos_c))
# 计算中心角 c
c = math.acos(cos_c)
# 计算距离
distance = R c
return distance
示例：北京到纽约
dist = spherical_law_of_cosines(39.9042, 116.4074, 40.7128, -74.0060)
print(f"距离: {dist:.2f} km")常见问题解答 (FAQ)
球面余弦定理在短距离计算时可能因浮点数精度问题产生误差，而Haversine公式在短距离和全球范围内都更加稳定。但对于长距离（如跨半球）计算，球面余弦定理依然有效且直观。
地球是一个近似球体（椭球体）。在平面地图上，直线并不等于最短路径。球面几何中的大圆航线（Great Circle）才是两点间的最短路径，必须通过球面三角学公式计算。
在天文学中，它用于计算天体在天球坐标系中的角距离，例如计算两颗恒星之间的角距离，或者太阳、月亮在地平坐标系中的位置变化。
可以使用在线的大圆距离计算器（如Movable Type Studios提供的工具）输入相同坐标进行比对。或者使用Haversine公式作为交叉验证，两者结果应在允许误差范围内一致。