JED 2025年12月 EW 101 深读|Dave Adamy
一个方向角,如果没有说明从哪里量起,即使小数点后写了很多位,也可能被读成另一个量。Dave Adamy在JED 2025年12月的EW 101栏目中,用地面两点和一颗卫星组成的几何关系,解释方位、角度和距离。[1] 这篇短文的难点不是公式数量,而是几种相似表述之间的区别。
本文保留原刊的教学算例,重点说明三角形、角度基准、长度与舍入之间的关系,并复核纸面数字。这里讨论的是球形地球模型中的几何阅读,不把一次算例的结果解释成具体设备的指向精度、测量精度或任务性能。
原文使用SVP,即Sub-Vehicle Point,表示子星下点。在本文采用的球形地球模型中,在地心与卫星的连线上,取卫星下方的地表交点,就得到这个点。它把卫星在空间的位置与地表上的一个参照位置联系起来,但它本身不是卫星,也不是另一个地面观测点。
原刊把所考察的地面发射机或接收机统一称为“威胁”。这属于该栏目的应用语境;参与几何关系的仍是一个给定地面位置。为便于理解,以下在解释图形时称它为地面点。换一种称呼不改变经纬度、地心角或三角形的定义。
作者选用半径6371 km的球形地球,并在圆轨道假设下给出卫星到地心的距离10560 km。两者相减为4189 km,这是算例的卫星高度。6371、10560和4189三个数分别代表不同起点的距离,不能仅因单位相同就相互代入。
圆轨道假设在这里很关键。原刊把半长轴与卫星到地心的距离直接相等,是因为已经限定为圆。若只记住“半长轴就是当前位置半径”而忘记这一条件,就会把算例中的特殊关系误写成一般轨道定义。本文沿用其条件,不扩展到轨道预测。

图1 编辑依据原刊图2、图3整理的符号示意。球面三角形的边以角度表示,平面三角形的边以长度表示;两图均不按比例绘制。
第一组三角形由北极、子星下点和地面点构成。按原刊符号,三个顶点分别是A、B、C,大写字母表示顶点角,小写a、b、c对应其对边。这里最容易产生误会的是“小写就是长度”的习惯:在这一球面三角形里,小写边长用相应的地心角表示。
例如,a描述子星下点与地面点之间的大圆弧所对应的地心角。图画把这条边沿地表画出,并不意味着a已经变成千米。b与c分别是两点纬度的余角;因为一个顶点在北极,这样就能把经纬度资料与球面三角形的量联系起来。
原刊的例子取子星下点北纬30°、东经100°,地面点北纬45°、东经120°。因而A为20°,b为45°,c为60°。这些是同一教学实例的输入。A是北极处两个子午面之间的角,而地面两点之间的角距离是a,两者虽同为角,却不是同一个量。
图1把两种三角形的记号并列。左侧的连接线仅用于示意字母对应,不能拿普通平面三角形的内角和去解释它。读原刊图2时,应同时记住它画的是球面关系;将其抽象成纸面上的三个顶点,并不会使地球曲率自动消失。
原刊先使用球面边余弦定律,从b、c和A求a。公式中的所有三角函数必须采用同一角度约定。本例给出的30°、45°和20°均为度数,计算时若误把数字当作弧度,错误会先发生在函数输入阶段,后面增加有效位数也无法补救。
纸面上列出的取整乘加为0.707×0.5+0.707×0.866×0.940。复算得到0.92902628,而原刊写成0.930。即便只保留三位小数,前一结果通常也应写作0.929,因此这一处不能简单解释为正常的末位舍入。
若直接从原刊给出的角度重新计算,并保留中间精度,余弦值约为0.9289952496,相应的a约为21.7213°。原刊后续用0.930得到21.57°,与前面的输入不再完全一致。本文把这两条计算路径区分开,避免默默更换数字却不说明原因。

原刊球面边余弦定律。A为北极处的角;a、b、c均为球面边对应的地心角。本算例以度为角度单位;a不是千米距离。
得到a后,作者用球面正弦定律求B。B位于子星下点,表示与指向北极的大圆方向之间的角。在原刊给定的这个例子里,它与所讨论的方位角对应;按一致的中间精度复算,得到约40.8039°,原刊印出的结果为41.08°。
这里还有一个比末位数字更根本的问题:正弦值本身不能唯一识别所有方向。一个锐角和它的补角具有相同正弦。计算器给出的反正弦主值只是约定范围内的一个结果,不能用“按键有答案”代替几何位置判断。原刊算例的相对位置支持所取的锐角,但该步骤不能不加条件地搬到任意两点。
方位还涉及从哪个北向开始、向哪个方向计角,以及取值范围如何定义。原刊图示足以帮助理解本例,却不等于一套适用于全球所有位置的完整方向约定。本文强调这点,是为了保留算例成立的前提,而不是另行给出天线控制或定位算法。
在概念上,方位角是在相应当地水平参考中描述方向,天底夹角则描述视线偏离径向向下方向的程度。两者处理的是不同维度的关系。原文有关方位平面的语句比较紧凑,阅读时应与图1的水平面标注一起理解,不宜仅从一句话抽出孤立定义。

原刊球面正弦定律的此例形式。分子和分母均为无量纲的正弦值;由反正弦求B时,还须结合给定几何关系判断分支,不能把主值当作任意位置的完整方位答案。
接下来,作者把视角转入包含卫星、地面点和地心的一个平面。三个顶点分别是E、F、G。现在小写e、f、g表示普通长度:e是6371 km的地球半径,f是10560 km的地心距离,g是卫星与地面点之间的直线距离,也就是本例的斜距。
连接两组图形的是地心处的角G。它与前面球面三角形的边a表示同一个地心角。一个量从第一组图里的“边”变成第二组图里的“角”,数值却没有改变。这正是作者连续使用两类三角形的理由,也是在抄录公式时最需要保留的对应关系。
原刊采用平面余弦定律。将图形与公式一起核对可见,余弦项前是负号,原刊后面的根式也使用负号。用原刊的0.930计算,斜距约为5192.94 km,与其写出的5193 km一致;若采用从原始角度得到的余弦值,则斜距约为5205.94 km。
差异约13 km,来源在于上游角度计算路径,不能据此宣称某种实际测距设备有13 km误差。本文只复核一组纸面算例。这个区别也说明,末端结果看似与公式吻合,并不能证明整条计算链已经自洽,仍需追溯它使用了哪个中间值。

原刊平面余弦定律。e为地球半径,f为卫星至地心距离,g为两点直线斜距,三者统一用km;G是地心夹角,与球面图中的a相同。公式两侧均为长度的平方。

图2 编辑据原刊“from nadir”的定义绘制。E在卫星处从天底方向量起;地心角与地面仰角的顶点不同,不能只凭角名互换。
原刊图1和图3的图注都强调“from nadir”,即从天底方向起算。图3中E位于卫星,它的两条边分别通向地心和地面点,因此E描述视线与卫星径向向下方向之间的夹角。本文用“天底夹角”表述,避免只写“俯仰角”而丢掉零点。
地面观测者常说的仰角,则以地面点为顶点,并相对于那里的当地水平面起算。它和卫星处的E既不在同一顶点,也不共享同一基准线。不能看到两处都出现elevation一词,就把两个数当成彼此可直接替换的量。
在同一个顶点,水平参考与竖直参考之间可以建立几何关系;但把角的顶点换到地面之后,还必须考虑整个三角形。图2只标出原刊涉及的基准差别,不给不同地点的角随意套一个固定换算。角名必须连同“谁在看、从哪里量”一起阅读。
原刊再次使用平面正弦定律,按其取整路径得到26.8°;用前面保持中间精度的a与g复算,E约为26.9309°。这是同一球形地球例子的纸面差异,不是对真实观测方向所作的精度承诺。

原刊平面正弦定律的重排。E为卫星处的天底夹角,G为地心角,e与g使用相同长度单位;右侧无量纲。文中数值仅复核原刊给定实例。
公式可以先放到一个容易想象的位置关系中检查,而不必一开始就相信长串小数。若地面点恰好位于子星下点,球面上的角距离为零;地心、地面点和卫星在同一直线上。此时斜距应等于卫星高度,而不应等于卫星到地心的距离。这个几何判断可以帮助辨认e、f、g是否被误放。
同样,在这一位置关系下,视线与天底方向一致,天底夹角应为零。如果有人把图中的E解释为相对于卫星当地水平面的角,就会得到另一种起算值。用这个特例检验定义,比单独争论“俯仰角”该怎样翻译更直接:先让文字和图形描述同一条线。
但方位角在这个特例中没有唯一方向。地面点没有相对于子星下点的水平偏移,不能要求图中仍保留一个唯一的“朝向哪边”。因此,当某个三角函数表达式出现分母趋近于零时,应先检查图形是否退化,而不是任意补一个零值。本文将这一点作为理解公式适用条件的编辑补充。
距离也可以作类似区分。地表大圆弧描述沿球面走过的路程,斜距描述空间中两点的直线间隔,卫星高度描述卫星与其正下方地表之间的径向距离。原刊的a、g和高度分别参与这些不同描述;不能把“从卫星看向某处”理解成沿地表量一段距离。
这些自检并不增加新的测量资料,也不会消除输入误差。它们的作用是找出定义和公式之间的矛盾,例如把地心距离当高度、把地面仰角当天底夹角,或在退化图形上强行保留方向值。先处理这类问题,再比较末位数字,复算的目的才是澄清概念,而不只是让计算器输出更多位数。
表1把原刊的纸面数值和从同一输入重新计算的数值并列。两列具有不同作用:左列帮助定位原刊内容,右列检查数学自洽。右列多保留几位,是为了使读者看清差异如何传递,并不意味着输入资料真的具有同样多位的准确度。
原刊选择6371 km作为球形地球半径,已经作了模型层面的简化。即便对这组固定数字进行非常精细的算术运算,也不能把结果称为真实地球表面的高精度位置解。计算精细程度与模型逼近现实的程度,应当分别讨论。
还有一种容易忽略的混用,是一处保留完整计算值,另一处把三位小数的显示值重新输入。此时同一份文字可能在不同步骤使用了不同精度。整理技术算例时,比较清楚的做法是说明哪些数只用于展示,哪些数实际进入后续计算,再统一给出最后的展示位数。
因此,复核不能只问结果“差得大不大”,还应问它能否由所列输入和步骤重现。本文采用原始角度保持中间精度的路径,并同时保留原刊数值以便对照。发现差异后的解释,应准确落在具体算术步骤上,不能将其泛化为整篇几何方法失效。
表1 同一原刊算例的数值回核
| 量 | 原刊印值 | 保持中间精度复算 |
|---|---|---|
| cos a | 0.930 | 0.92899525 |
| 地心角 a | 21.57° | 21.7213° |
| 本例方位角 B | 41.08° | 40.8039° |
| 斜距 g | 5193 km | 5205.94 km |
| 天底夹角 E | 26.8° | 26.9309° |
原刊把3小时周期、4189 km高度的圆轨道称为低地球轨道。ESA在2022年有关导航任务的说明中,以低于2000 km描述低地球轨道。[2] 按这一常见口径,原刊该例不宜继续称作低地球轨道。本文保留其高度和地心距离作数学输入,并明确指出名称与该口径不符。
这一名称问题不会自动改变三角形的边长,却会影响读者对案例背景的理解。要修正的是分类表述,不能为了让名称看起来一致,擅自把4189 km改成另一高度。否则,就已经换了一道题,随后得到的结果也不再是对原刊算例的核查。
纬度同样有模型含义。原文说明的是地心纬度。在参考椭球模型下,大地纬度与地心纬度的定义并不相同;JPL的SPICE文档以参考椭球表面外法线相对于赤道平面的角定义相应的大地纬度。[3] 球面模型让这一差别消失,但现实资料并不会仅因写着“纬度”就自动采用同一约定。
这并不是要求在入门文章里加入所有复杂修正,而是要求尊重简化的边界。球形地球、圆轨道、指定经纬度与给定时刻的位置,足以构成本文的教学问题。若读者拿到另一份带有参考椭球、时刻或坐标框架的数据,应先确认定义是否兼容,再讨论能否沿用相同表达。
这篇EW 101的价值,在于展示两个三角形如何通过同一个地心角联系起来:球面关系处理地表两点的位置,平面关系处理卫星、地面点与地心之间的长度和夹角。它们的字母习惯相似,量的性质却不同,理解这点比背下最后几个数字更有用。
回读原刊时,可以先在图上找到每个角的顶点和基准,再确认每条边是角度还是长度,最后检查取整值怎样进入下一步计算。这样既能保留原文的教学线索,也能指出具体的算术与术语问题,使一组看似熟悉的三角公式成为可解释、可复核的几何知识。
[1] Dave Adamy. Space EW (Part 7) – Location of an EW Threat. JED,2025年12月,第67—68页,图1—3及纸面算例。
[2] ESA. New navigation missions for enhanced satnav and Earth mapping,2022年7月6日。仅核对低地球轨道低于2000 km的表述。
[3] NASA JPL/NAIF. SPICE PGRREC程序文档,Detailed Input及Particulars中的参考椭球与纬度定义。
[4] 本文两图为编辑据原刊定义重绘的示意;数值表为原刊输入的算术复核。模型简化、反三角函数分支及精度说明为编辑解读。