卫星几何中的角距离斜距与视距

JED 2025年8月 Dave Adamy卫星几何文章深读

已知卫星正下方的位置、卫星高度和一个地面点，能否算出两者相距多远？Dave Adamy在2025年8月JED的空间电子战连载第三篇中，用球面三角形和平面三角形连接起这个问题，并进一步计算自由空间损耗。[1] 这是一个适合复习坐标、几何与单位的例子。

但距离算出来，传播条件还没有自动成立。按原刊给定的高度和经纬度复算，卫星方向实际低于地面点的几何地平线，连接两点的直线受到地球遮挡。本文保留原来的假设输入，先说明计算怎样展开，再补上原文未检查的视距条件。这个区别，决定了最后的损耗数字应该怎样理解。

01 星下点先把空间位置投到球面

Adamy把地球近似为球体，从地心向卫星作一条射线，与地表相交的位置就是星下点，原文缩写为SVP。本文用O表示地心、S表示卫星、V表示星下点、P表示地面点。O、V、S在同一直线上，卫星高度H就是从V向外到S的径向距离。

因此，卫星位置在这个模型中可以拆成两个部分：星下点的经纬度给出方向，高度给出离球面的距离。只知道星下点，并不能确定卫星距地面点有多远；还需要知道它沿这条地心射线处在什么位置。

原文明确使用地心纬度，即地心到该点的方向与赤道平面的夹角。日常坐标资料也可能给出以参考椭球法线定义的大地纬度。NASA/JPL的SPICE文档对后者有清楚定义。[3] 对球体，两种方向重合；改用椭球时，不能把“纬度”两个字当成可以省略定义的通行证。

同样，原文关于纬线不是大圆的说明需要保留赤道这一例外。球面上的经线是过地心平面形成的大圆的一部分；除赤道外的纬线是小圆。后面使用的球面三角形边，沿的是大圆弧，不能直接拿沿纬线走过的弧段代替。

02 同一张图上有四种不同的量

图1 编辑自绘的几何量示意。为区分各条线，图中高度作了放大，未按原刊算例比例绘制；H、s、d、a分别表示径向高度、地表弧长、斜距与地心角。

图1把容易混淆的量放在一起：H是卫星的径向高度，s是V与P之间的地表大圆弧长，d是S与P之间的直线距离，a是OV与OP之间的地心角。它们都参与同一个空间关系，但端点、路径和单位并不相同。

高度回答卫星离星下点多远，斜距回答它离所讨论的地面点多远。当地面点偏离星下点时，两者一般不相等。地表弧长则只连接球面上的两个位置，并没有到达卫星。把地表距离直接放进卫星传播公式，会在选择几何对象这一步就发生错位。

地心角又是另一种量。球面三角形常用角度表示边长，是因为一段大圆弧可以由地心张角来描述。这并不表示角度可以直接按千米使用；只有选定球半径并完成弧度换算，角距离才变成球面上的长度。

03 经度差与两个纬度确定地心角

原刊选取北极、星下点V和地面点P作为球面三角形的三个顶点。由北极到V的大圆弧对应余纬90°减去星下点纬度，由北极到P的弧对应另一余纬；两条经线在北极形成的夹角，则由经度差确定。

这样，就可以通过球面余弦定理求出VP所对应的地心角a。把原文的余纬表达式换成纬度表达式，得到下面的等价形式。φs和φp分别代表星下点与地面点的地心纬度，Δλ是采用同一约定的经度差。

原刊给出星下点经度200°、纬度45°，地面点经度230°、纬度20°。这里按同一正方向、0至360°范围的经度记法理解，两者经度差为30°。这是一组教学输入，原文没有将其对应到某一颗真实卫星或具体地面站。

经度改写到另一常用范围时，要对两点采用一致的规则。这里需要的是相对角关系，而不是孤立比较200与230两个数字有多大。计算器若处在度模式，可以直接使用角度；程序库通常接收弧度，输入前必须明确转换。

球面模型下的地心角关系。所有三角函数的角度单位必须一致；求得余弦值后，反余弦给出0至180°范围内的较小地心角。

04 保留精度能解释距离差异

原刊把三角函数值先保留三位小数，得到cos a约为0.818，再求出a约为35.1°。如果从给定的45°、20°、30°出发，在中间步骤保留计算精度，则cos a为0.8172866216，a约为35.1859°。两种结果的差别首先来自舍入路径。

这一点值得说明，因为把复算值和原刊值直接并排，却不交代中间精度，会让读者误以为余弦定理本身发生了问题。按0.818取反余弦确实约为35.1149°，原刊写35.1°在这一舍入前提下是相符的。

更稳妥的计算方式，是让内部数值保持足够精度，最后按表达需要统一舍入。本文展示若干小数位，是为了追踪算例差异，并不意味着给定的地球球体模型、卫星高度和位置具有同样多位数的测量准确度。

从角距离求地表弧长时，用球半径乘以弧度值即可。本例得到约3912.49千米。这个量有助于理解地表两点的分离，但下一步求卫星斜距时，需要回到地心、卫星和地面点组成的平面三角形。

s为球面上V到P的较短大圆弧长，R为球半径，a必须以弧度代入。R取6371千米时，本例s约为3912.49千米。

05 斜距由平面三角形给出

球面计算确定方向关系后，问题转到由O、S、P组成的平面。OP长为R，OS长为R加H，两者夹角为a。使用平面余弦定理，便得到卫星与地面点之间的直线距离d。这里的半径与高度必须采用相同的长度单位。

原刊采用R为6371千米、H为1000千米，因此卫星到地心的距离为7371千米。若沿用已舍入的cos a等于0.818，算出的d约为4253.67千米，取整数就是原刊的4254千米。这个结果在其自身的数值路径中可以复现。

若始终使用未提前舍入的余弦值，d则约为4261.53千米，差约7.87千米。表1把两种计算口径放在一起。差异并不是因为两种方法描述了不同卫星，而是同一教学输入在中间精度上的处理不同。

余弦定理只计算两点的欧氏距离，无论连接两点的线段是否穿过球体，公式都能给出正的结果。也就是说，得到一个合理大小的斜距，并不能证明这条直线可以成为无遮挡的传播路径。可见性需要另外检查。

平面三角形OSP中的斜距公式。R、H和d使用同一长度单位；本式计算直线距离，本身不检查地球遮挡。

表1 同一原刊输入在不同舍入口径下的复算

06 地平线条件改变算例的解释

图2 按原刊原始数值复算的局部剖面，编辑自绘，横纵比例相同。卫星在画面外；红线表示从P指向卫星的直线方向，浅色区域是球体内部。

现在补上原刊没有计算的条件。仍沿用球体、地面点高度为零、传播路径按直线理解的假设。在地面点P处，垂直于OP的切线代表当地几何水平线。只有卫星位于这条水平线之上或恰在其上，连接方向才不会立即进入地球内部。

将卫星位置投影到OP方向，投影长度为(R加H)乘以cos a。它至少要达到R，才满足上述条件。因此，给定高度下，P与星下点之间允许的地心角存在一个几何边界。边界处OSP构成以P为直角的三角形。

把本例R为6371千米、H为1000千米代入，视距边界的地心角约为30.1933°。而原刊输入所对应的a约为35.1859°，已经超出边界。即便采用原文35.1°的舍入数，判断也不会改变；这不是补几位小数就能消除的差异。

图2显示P附近的局部剖面。指向卫星的直线低于当地水平线，并在起始部分进入球体内部。这个几何事实与斜距计算同时成立：距离是可以算出的，但原文随后默认的无遮挡视距条件没有得到满足。

编辑补充的几何视距条件。适用于球面上高度为零的地面点和直线传播模型；a_h是切线边界，不包含地形或传播环境的附加影响。

07 用地面仰角交叉核对几何判断

还可以直接从地面点观察卫星的方向。在OSP所在平面内，卫星相对P的水平分量为(R加H)乘以sin a，竖直分量为(R加H)乘以cos a再减去R。两者共同确定地面仰角。

下面采用atan2形式，第一项为竖直分量，第二项为朝向卫星的水平分量。它保留符号信息。本例得到仰角约为负4.6676°，表示卫星方向位于几何水平线以下，与前面的地心角边界判断一致。

为核对代数与单位，本文还将两点的球面经纬度分别换成三维直角坐标，再直接计算位置差的长度和相对当地竖直方向的分量。独立构造的向量结果与上述斜距、仰角一致。这种交叉核验针对的是算例数学，没有引入真实轨道或实测传播数据。

原刊的图3用于说明平面三角形关系，不能代替对给定数字的可见性判断。阅读示意图时，线段画在地球外侧不代表所有代入值都能满足同样的几何条件。图形帮助理解变量，数值仍需接受条件检查。

ε为地面点的几何仰角，以当地水平线为零。第一参数为竖直分量，第二参数为水平分量；本例ε约为−4.6676°。

08 损耗数值依赖传播前提与单位

原文随后采用2 GHz频率，代入自由空间损耗表达式，得出约171 dB。ITU-R P.525-5将自由空间传播作为一种理想参考条件，并给出点对点基本传输损耗公式。[2] 使用这个表达式，需要把所计算的空间距离与公式适用的物理条件联系起来。

为复现原刊，下式保留其32.44常数。ITU文件以兆赫和千米表示时将常数取为32.4；这约0.04 dB的取整差别不影响本例的主要判断。真正重要的是频率要写成2000 MHz、距离以千米代入，两个对数都以10为底。

如果只作形式代入，使用4261.53千米会得到约171.05 dB，使用原刊4254千米则约为171.04 dB，两者取整数都是171 dB。然而，前面已经确认原算例受地球遮挡，所以不能把这个接近的数值当成有效视距链路的证明。

还要注意原文“to and from satellites”讨论的是往卫星或从卫星进行传输，并不是同一信号往返一次的雷达过程。这里的表达式是单程点对点自由空间损耗，不能因为中文出现“往返”就将其解释为两程损耗。

自由空间损耗本身也不等于完整链路预算。本文没有加入发射功率、天线增益、接收条件或实际环境资料，因此不据此判断某一系统能否通信、截获或干扰。对这篇基础几何文章，先把传播条件和单位说清楚，已经能够纠正最关键的阅读误差。

L_fs以dB表示，d为单程点对点距离，F为频率。分母明确输入单位。本例的约171.05 dB仅为公式代入结果，未满足无遮挡视距前提。

09 把教学模型的边界保留到结论

Adamy这篇文章展示了一条清楚的数学连接：用球面关系获得地心角，再用平面关系获得空间斜距。补充视距条件后，这条连接才能被正确用于解释传播。地球取球体、地面高度取零、角度和距离采用统一单位，都是结论的一部分。

这组输入只描述某一瞬间的位置关系。原文没有提供轨道要素、时间序列或真实卫星身份，因此不能由一个斜距继续推出过境时间、轨迹变化或覆盖时长。文章预告下一期再谈地面轨迹，也说明静态几何和随时间变化的轨道问题各有范围。

此外，球体模型中的几何地平线并不等于某一现场的全部可用条件。地面高度、地形遮挡以及传播环境，都可能需要另外描述。本文补充的检查只是原有假设内的一项必要核验，没有将它包装成完整的实际链路评估。

读者可以据此形成一种更稳健的阅读习惯：先确认每个量对应的端点和方向，再检查单位、舍入与适用条件，最后解释数值。即使一连串公式都算得通，只要最后借用的物理前提没有成立，结论的措辞就必须停在数学结果能够支持的位置。

让几何计算与传播条件对应

本例的收获不只是得到约4262千米的斜距。地心角、弧长和空间距离各有定义，舍入决定复算细节，而地平线检查决定视距解释能否成立。将这些条件一起保留下来，才能读懂原刊的教学意图，也看见算例需要补足的地方。

资料来源与说明

[1] Dave Adamy. Space EW (Part 3) – Satellite Geometry. JED，2025年8月，第39—40页。

[2] ITU-R. Recommendation P.525-5，Calculation of free-space attenuation，2024年11月，附件1第1节与第2.3节。

[3] NASA/JPL NAIF. SPICE Toolkit，recgeo_c，Rectangular to geodetic coordinates，参数定义与Particulars说明。