JED 2026年6月|卫星观测几何中的角度与距离

2026-09-26

卫星观测几何中的角度与距离

JED 2026年6月EW 101深读与算例核验

从卫星向地面看,一条视线可以对应几个不同的角度:在地心量到的是地心角,在卫星上量到的是离天底角,在地面量到的是仰角。它们描述同一个空间关系,却有不同的顶点和零方向。如果把这些角度混在一起,即使公式外形熟悉,也容易得到貌似合理的错误结果。

Dave Adamy在JED 2026年6月的EW 101专栏中,讨论卫星地平线边界和定向天线的观察方向。本文围绕其中的球形地球教学模型,解释几何关系,并核对原刊两页中的算式。原刊有若干符号、代入和算术结果不一致,尤其是离天底角与地面仰角;下文将它们与按同一输入复算的结果分开列出。[1]

一、先分清在哪里量角度

图1 依据原刊几何关系自绘。地心角、离天底角和地面仰角分别在不同顶点量取;示意图夸大卫星高度,不按算例比例绘制。

图1 依据原刊几何关系自绘。地心角、离天底角和地面仰角分别在不同顶点量取;示意图夸大卫星高度,不按算例比例绘制。

先把地球当作半径固定的球。在这个模型中,从卫星沿地心方向向下,交到球面的点就是星下点。卫星、星下点和地心在同一条直线上。原刊把星下点记为SVP,并以它的经纬度说明卫星在地球表面上方的位置。这里暂不讨论地球椭球形状、地形高度或传播折射,这样才能把注意力集中在基本角度关系上。

离天底角从卫星指向地心的方向量起,描述视线偏离正下方多少;地面仰角从地面点的局部水平面量起,描述卫星高出当地地平线多少。两条局部参考方向位于不同位置,因此这两个角一般不是简单的互余关系。还必须考虑地心角,也就是星下点和地面点对应的两条地心半径之间的夹角。

原文把从nadir量起的角度称为elevation。中文如果只译成“仰角”,容易让读者误以为它以局部水平面为零点。本文统一称它为离天底角K,把地面观察者从局部水平面向上量取的角度称为地面仰角ε。名称的作用是保留参考方向,而不是争论哪一种英文说法更常见。

方位角也需要参考系。NASA喷气推进实验室的SPICE文档将方位、仰角与观察者的局部参考坐标系联系起来,并明确零方向、正方向及单位。它说明了一个基本原则:单独给出“多少度”还不完整,必须同时说明从哪里向哪个方向量。本文只借用这项定义原则,没有使用SPICE进行原刊算例计算。[2]

二、地平线是切线条件给出的几何边界

原刊先求卫星看到的地平线。连接卫星、地心与地平线点,可以形成一个直角三角形:在理想球面上,刚好擦过地球的视线与该点半径垂直。这个直角是计算的依据,并不是凭图形外观作出的近似判断。若视线进入球体,地球会遮挡;若视线没有接触球体,它也就没有交到地面。

设地球半径为R,卫星高度为H,卫星到地心的距离为R加H。由三角形得到星下点至地平线边界的地心角C。它是一个角度,不是地面距离,也不是卫星离天底方向的角度。原文采用R为6371千米、H为300千米,算得C约17.2度,这一数量级与复算相符。

同一球面模型里,以星下点为中心、与它保持这一地心角的点构成地平线轮廓。将轮廓写成经纬度表,是把空间中的几何边界转换成地图可表达的位置。这个轮廓在球面上具有简单关系,但投到平面地图之后的形状取决于投影,不能凭屏幕上的圆或椭圆反过来判断空间距离。

这里的“可见”仅表示球形地球没有挡住几何视线,并没有证明接收设备能够获得足够清晰的信号,也没有说明天线能覆盖整个区域。地平线轮廓与有效观测范围不是同一项能力。原刊在本篇主要建立前者所需的几何关系,后续性能问题仍需其他信息。

C=\arccos\!\left(\frac{R}{R+H}\right)

R与H使用相同长度单位,C为星下点至几何地平线的地心角。采用原刊6371千米与300千米,复算C约17.2482度;公式基于球面、直线传播和无遮挡假设。

三、从地心角走到经纬度 要看球面三角形

为描述边界上的某个点,原刊构造了北极、星下点和该边界点之间的球面三角形。这里的三条边不是平面上的直线长度,而是对应大圆弧的角长度。星下点余纬等于90度减去它的纬度,边界点余纬也以同样方式定义;另外一条边则是上一节求得的地心角C。

原刊选取星下点东经100度、北纬30度,并选择相对当地北方向向东45度的边界点。因此星下点余纬e为60度,方向角F为45度,球面角边d等于C。由球面余弦关系先求边界点余纬f,再换回纬度。按原始输入保留精度复算,得到纬度约北纬41.2305度,与原刊北纬41.2度的表达一致。

接下来还需要经度差。原刊通过球面正弦关系求得约16.1度,再与星下点东经100度相加。按未经提前舍入的输入,本例经度差约16.1878度,即边界点约东经116.1878度。小的差异应结合原文逐步舍入理解,不应与后面明显不满足算式的结果混为一谈。

不过,原文说把F从0度转到360度,就能重复计算整个边界,这一步不能只机械重复反正弦主值。正弦函数不能独自区分全部象限,跨经度边界也需要一致约定。本例处于便于判别的东北方向,读者容易忽略这一点;换到其他方向,若没有分支判断,可能得到经度镜像或不连续的轮廓。

这一提醒不要求读者记住更多公式,而是认识到“某个点算对”与“整个轮廓都能正确生成”之间还有一步。把可用的单例公式推广到全方向时,需要确认角度范围、符号和参考方向是否仍然成立。原文的教学例子提供了关系,但不应被误当成已经完成所有边界处理的通用程序。

\cos f=\cos e\cos d+\sin e\sin d\cos F

e与f分别为星下点和边界点的余纬,d为两点间的地心角,F为星下点处的球面夹角。所有三角函数使用同一角度单位;边界点纬度为90度减f。

四、地面点的角距离不能用经纬差直接相加

原刊随后另选一个位于东经102度、北纬32度的地面点,讨论它相对星下点的方向。它与星下点的纬度差和经度差都是2度,但两者不能直接相加,也不能当成同一平面上的两条等长边。纬度、经度描述球面位置,所需的是两个位置向量之间的地心角。

沿用原刊记号,星下点与地面点的余纬分别为h和j,经度差为G,地心角为g。本文为避免与卫星高度H混淆,仅在这一球面关系中保留小写h。代入60度、58度和2度后,cos g约为0.9989434,得到g约2.6341度。它远小于本例地平线地心角,所以该点位于理想球面的几何可见范围内。

原刊在中间步骤写成反余弦0.738等于2.6度,这是不成立的:反余弦0.738约为42.4度。它既不等于前面三项乘积相加的结果,也不可能给出所写的小角度。这里应作为原文数值笔误或计算记录错误注明,而不是悄悄复制成一个看似有推导过程的结果。

方位角J在原刊例子中约为41.3度,但这个值采用了较粗的中间舍入。按完整输入复算,本例约为40.0913度。两者的差异说明,当角度通过多个三角函数关联时,过早舍入会在后续计算中放大。对读者而言,与其逐行沿用三位小数,不如保留计算精度,只在最终展示时进行统一舍入。

\cos g=\cos h\cos j+\sin h\sin j\cos G

h与j为两点余纬,G为经度差,g为星下点与所选地面点之间的地心角,均为角量。本例h=60度、j=58度、G=2度。

\sin J=\frac{\sin G\sin j}{\sin g}

J为星下点处相对当地北方向的球面方位角。本式在本例已知东北象限条件下取对应解;正弦关系单独不足以确定任意位置的全象限方向。

五、斜距与离天底角 必须使用同一个三角形

求得地心角之后,可以回到卫星、地心和地面点组成的平面三角形。两条半径边分别为m与k:m是卫星到地心的距离6671千米,k是地球半径6371千米;另一条边n是卫星至地面点的直线斜距。夹在m与k之间的角,就是刚才算出的g。

平面余弦定理的交叉项应该是2mk乘以cos g。原刊符号式却写成2nk,等于把要求解的斜距放进了错误的位置;数值代入行的首项又出现6617,而此前定义为6671。它们不是不同地球模型的选择,而是同一页内部的符号与数字不一致。先让定义和算式对应,再谈计算结果,才不会把排印错误传播到后续推导。

按完整输入复算,斜距约424.0405千米。原文写418千米;如果仅把g提前舍入成2.6度,但其他数字保持正确,斜距约421.3千米,也不是418千米。因此,应说明本文采用的完整输入与统一计算方式,不能把几个来自不同中间步骤的数字混在一张结果表里。

离天底角K位于卫星顶点,它对面的边是地球半径k。由正弦关系可求K约43.6679度。原刊把反正弦计算结果写成68.6度,与本例的边长和夹角不相容;即使直接采用它给出的简化代入6371乘0.045再除以418,反正弦也只有约43.3度。这个差异已经不是保留一位小数可以解释的。

n=\sqrt{m^{2}+k^{2}-2mk\cos g}

m、k、n均为长度;m=R+H,k=R。g是m与k的夹角。采用完整原始输入,n约424.0405千米。

\sin K=\frac{k\sin g}{n}

K是卫星处视线偏离天底的角度,本例可见几何中取相应锐角,约43.6679度。公式中的n必须与同一组g、m和k一致。

六、三个角度相加正确 仍不代表算例已经正确

图2 对原刊算例的核验次序。定义、公式、分支与数值回代各自检查不同问题,不能由单一等式代替全部核验。

图2 对原刊算例的核验次序。定义、公式、分支与数值回代各自检查不同问题,不能由单一等式代替全部核验。

在卫星、地心和地面点组成的平面三角形中,地面点处的内角M,加上K与g,应等于180度。地面点指向地心的方向与局部水平面相差90度,因此卫星在地面看来高出水平面的仰角为M减90度。把这两个关系合并,可以看出地心角也参与了离天底角与地面仰角之间的转换。

采用一致复算值,M约133.6980度,地面仰角约43.6980度。原刊由68.6度的K继续得到M为108.8度,再算出18.8度的地面仰角。最后两步在所给数字之间是连贯的,但上游K已经不符合边长关系,所以连续计算并没有把错误纠正回来。

这正是学习算例时容易遇到的陷阱:让三角形内角和等于180度,只能验证一条关系。若其中一个角来自错误代入,另外一个角用减法补出来,内角和当然仍然正确。还需要检查正弦关系、余弦关系以及角度所处的几何范围,才能确认整组结果一致。

下面的表格将原刊结果与同一输入下的复算并列。复算列保持一套精度路径,最终统一展示四位小数;它表示理想球面模型中的数学结果,不是对真实卫星位置或姿态精度的声明。真实数据的误差不能因为计算器多显示几位小数就自动消失。

M=180^\circ-g-K,\qquad \varepsilon=M-90^\circ=90^\circ-g-K

M是地面点处指向地心与指向卫星两条射线的夹角;ε为以局部水平面为零点的地面仰角。所有角按度表示,这里仍采用球形地球与直线视线模型。

表1 同一原刊输入下的数值核验

几何量原刊列值统一精度复算
地平线地心角C17.2度17.2482度
边界点纬度北纬41.2度北纬41.2305度
边界点经度东经116.1度东经116.1878度
所选点地心角g2.6度2.6341度
所选点方位角J41.3度40.0913度
斜距n418千米424.0405千米
离天底角K68.6度43.6679度
地面仰角ε18.8度43.6980度

七、把几何可见与系统能力放在不同层次理解

这篇文章虽然以定向天线观测为背景,但所完成的主要工作是把位置关系转成角度和距离。一个地面点位于地平线以内,并不能单独说明设备将它纳入了观测;一条视线方向已经定义,也不能单独说明天线在该方向有足够的响应。球面几何解决了关系中的一部分,不能替代全部系统判断。

星下点经纬度只描述某一时刻的位置条件。实际卫星会运动,姿态与参考系也需说明。本篇没有具体轨道、姿态时间序列或天线安装关系,因此教学算例的方位角不能直接等同于某台机构的控制指令。

从学习角度看,读者可以用几个直观情形检查理解:地面点趋向星下点时,视线趋向天底方向,地面仰角趋向直上方;地面点趋向理想地平线时,地面仰角趋向零。这些极限情形帮助判断角度名称有没有混淆,却也仍需与相应模型条件一致。

原刊预告后续将把角度与距离用于进一步的射频计算,但本篇并未给出完整链路性能。读完这一期,更重要的收获是能够分辨三个角度、看懂球面与平面三角形怎样衔接,并知道如何识别公式记录中的不一致。几何关系先保持自洽,后面的讨论才有可靠的起点。

资料来源与说明

[1] Dave Adamy:Intercept from Space using a Directional Antenna,Space EW Part 13,JED,2026年6月,第49—50页。

[2] NASA/JPL NAIF:SPICE azlcpo_c文档,Detailed Output、Particulars,方位、仰角及参考系定义。

阅读7
分享
写评论...