JED 2026年3月|卫星飞得很快,为什么频移不能直接按速度算

2026-09-26

卫星飞得很快,为什么频移不能直接按速度算

JED 2026年3月|从多普勒算例读懂视线、单位与近似

一颗卫星以每秒数千米的速度运行,地面收到的信号会偏移多少?直觉很容易把轨道速度直接代入公式,但多普勒效应首先关心的是两端距离怎样变化。速度大小只是信息的一部分,方向和参考系同样重要。

JED 2026年3月的EW 101专栏中,Dave Adamy以三小时周期的圆轨道为例,讨论星地链路的多普勒频移。[1] 本文沿着这条思路解释基础物理,并对原刊算例中的单位、几何近似和最终公式逐项核验。原刊计算、公开资料与编辑推导分别标明。

一 先看两端的距离,别急着看轨道速度

图1 编辑示意。正负号按“接收频率减发射频率”定义,横向情形只表示该瞬间的几何一阶项为零。

图1 编辑示意。正负号按“接收频率减发射频率”定义,横向情形只表示该瞬间的几何一阶项为零。

设想一颗卫星正在高速运行,地面天线却始终指向天空中的同一位置。两件事并不矛盾:相对地心惯性空间的运动,与相对地面接收站的距离变化,是两个问题。把某个坐标系中的速度大小直接放进多普勒公式,可能在第一步就选错了物理量。

地球静止轨道提供了一个很好的反例。ESA的说明指出,理想地球静止卫星在赤道上空沿地球自转方向运动,周期与地球自转相同,因此相对地面保持固定位置。[2] 编辑据此判断,在这一理想模型中,固定地面站与卫星之间的距离不变,运动引起的一阶几何多普勒项为零。卫星仍有很大的轨道速度,二者并不冲突。

这不表示真实链路的全部频率偏差都消失,也不表示所有地球同步轨道都具有同样条件。同步只约束周期;要相对地面静止,还需要圆形、赤道面及同向运动等条件。实际接收频率还可能受到频率基准等因素影响。原刊以GEO的高速运动引入主题时,必须补上这个区分。

NASA/JPL的教学材料把关键说得很明确:频移对应视线方向的运动;接近时接收频率上升,远离时下降。[3] 火车汽笛可以帮助建立直觉,但声音类比并不替代电磁波公式。尤其在两端都运动时,应先问距离在增加还是减少,再讨论变化得有多快。

二 径向速度要同时说清方向与正负号

原刊将径向速度解释为发射机与接收机之间的距离变化率,又用它的幅值写频移公式。理解趋势时问题不大,真正代数运算时却容易混淆:如果把“距离增大”定义为正,就不能再把正号自动解释为频率升高。本文采用一套明确约定,把接收频率减去发射频率定义为频移。

下面是编辑为统一符号给出的低速、单程、几何近似。ρ表示两端距离,ρ的时间导数为正意味着远离。单位视线向量u从接收端指向发射端,v_s和v_g分别为卫星与地面站在同一参考系、同一几何时刻的速度。先求速度差,再向视线投影,才得到距离变化率。

这个写法同时解释了三个容易误读的场景。物体横向掠过时,可以速度很快而瞬时径向分量很小;两端以相同速度平移时,相对距离不因共同平移而改变;卫星过顶的圆轨道理想瞬间,轨道速度和地面站自转速度都与当地竖直视线垂直,几何频移可经过零值。

还要分清距离和距离变化率。两个发射端可以离接收端远近不同,却在某一瞬间具有相同的径向速度;相反,距离相同也不意味着正在以同样的速度接近。低速几何式给出的频移只约束相对速度在一条视线上的投影,无法凭一个瞬时数值恢复完整的三维速度。若把多普勒当作距离本身,或者把径向速度当作三维速度大小,都会向观测量索取它没有包含的信息。

近似符号也有含义:这里不展开光行时修正、相对论高阶项以及传播介质和钟差。它用于理解原刊的基础几何,而不是完整精密测轨模型。公式中的负号来自本文的定义,并非在宣称原刊采用另一种符号就必然错误;需要避免的是在同一推导中悄悄更换定义。

\dot{\rho}=(\mathbf v_s-\mathbf v_g)\cdot\mathbf u,\qquad \Delta f=f_r-f_t\simeq-\frac{f_t}{c}\dot{\rho}

ρ为距离,u为从地面站指向卫星的单位向量;速度与光速c使用相同单位,频率均用Hz。此式为低速单程几何近似,距离增大对应负频移。

三 地面站也在运动,但不是把两个速率直接相减

Adamy接下来分别计算地面站与卫星的速度,这个拆解很有教学价值。固定在地表并不等于在地心惯性系中静止。地球自转使地面站向东运动;纬度越高,站点到自转轴的距离越小,所以线速度随纬度变化。原文图1表达的正是这一点,而不是站点到地心的半径随纬度缩短。

采用原刊的球形地球半径6371 km、恒星日86164 s和北纬30°,站点到自转轴的距离为地球半径乘纬度余弦,算得向东速度约402.3 m/s。原刊第40页将分母后面写成“86164 km/sec”,但上下文已经说明它是时间;这里按86164 s计算,并将千米换成米。

卫星例子选的是周期3小时的圆轨道。原表对应高度4189 km、轨道半径10560 km,二者相差地球半径6371 km。圆周长除以10800 s,得到约6143.6 m/s,按原刊舍入为6144 m/s。这里a从地心量起,不能用离地高度替代;表头的分钟也必须先换成秒。

两个计算得到的都是速度大小,尚未得到径向速度。把“6144减402.3”当作结果,就预设了两者共线且方向合适;再加上一个角度,只能处理该角度实际代表的投影。地面自转方向、卫星瞬时运动方向和两端连线,必须在同一几何图里辨认。

v_g=\frac{2\pi R_E\cos\varphi}{T_{\rm sid}},\qquad v_s=\frac{2\pi a}{P}

R_E为球形地球半径,φ为地面站纬度,T_sid为恒星日;a为圆轨道半径,P为轨道周期。长度用m、时间用s时,两式给出的速率为m/s,仍须向视线投影。

四 方位角算对,不代表空间投影已经完成

图2 编辑重新绘制的几何检查图。地心O、站点G和卫星S构成直角三角形;地平线方向与卫星轨道切向一般不平行。图仅解释投影,不给出完整过境轨迹。

图2 编辑重新绘制的几何检查图。地心O、站点G和卫星S构成直角三角形;地平线方向与卫星轨道切向一般不平行。图仅解释投影,不给出完整过境轨迹。

原刊再用球面三角形讨论过顶轨道的方向。对于纬度30°、倾角45°的例子,作者得到角B约54.7°,并把对应方向描述为正东向北偏35.3°。中文表达应为“东偏北35.3°”;写成“北偏东35.3°”会交换基准方向,尽管数字完全没变。

这一计算帮助描述轨道在站点附近的方向,却不能自动回答卫星位于地平线时,轨道切线与视线是否重合。过顶时的局部方向、地平面内的方位关系、远处卫星所在位置的切线,是不同层面的几何对象。原刊图3与图4应带着这个区别一起读。

更直观的检查可以从图2做起。地面站看见卫星恰在几何地平线上,意味着站点到卫星的连线与地球表面相切。卫星的圆轨道速度则与“地心到卫星”的半径垂直。两处对应的半径方向不一样,两条切线通常也不会相同。

因此,“卫星出现在地平线”不能成为“其全部轨道速度都沿着视线”的理由。原文把地平线区域视为地表的一小部分,再将来向与速度方向对齐,这是需要检验尺度的近似。一个看上去局部的地面图,未必能代表跨越数千千米的空间连线。

五 一个简单比例,就能检查近似是否足够小

为单独检验卫星项,暂时忽略地面站自转,只看球形地球和同一平面内的圆轨道。地平线位置的三角形满足OG垂直GS,OG长度为地球半径R_E,OS长度为轨道半径a。由这一直角关系,轨道切向速度在视线方向上的分量幅值,是轨道速率乘R_E/a。

这不是把原文的完整问题全部求解,而是一个局部一致性检查:原文若把该分量近似成全部速率,就等价于将R_E/a近似为1。对原刊采用的6371 km与10560 km,这个比例约为0.603。它明显没有接近1到可以不加说明地忽略。

按原刊圆轨道速率6144 m/s估算,单独的卫星视线分量在这一检查模型下约为3.71 km/s,而不是6.14 km/s。这里没有加入地面站自转,也没有重建30°纬度、45°倾角轨道的完整过境,因此不能把3.71 km/s再随意减去原刊另一个时刻或近似下的地面分量,拼成“修正版答案”。

这个检查的价值,是定位问题出在哪一步:前面的周长除周期可以算得准确,后面的视线投影仍可能不充分。算术、单位和几何分别通过检查,结果才有意义。仅仅增加小数位数,不能修复将一个方向上的速度当成另一个方向速度的错误。

\left|v_{s,\parallel}\right|=v_s\frac{R_E}{a},\qquad \frac{R_E}{a}=\frac{6371}{10560}\simeq0.603

编辑几何核验:球形地球、圆轨道、站点与轨道共面的地平线剖面;暂不计地面自转。左式仅为卫星速度对视线的投影幅值,不是完整的相对径向速度。

六 96.9 kHz能算出来,仍要问它代表什么

原刊先计算6144−402.3×cos35.3°,得到约5816 m/s,再用5 GHz乘以这一速度与光速之比,得出96.9 kHz。若接受前一步给定的径向速度,这个频率换算在算术上可以复现。问题在于,径向速度的几何来源需要成立,算术一致并不能独立证明模型正确。

因此,本文把96.9 kHz保留为“原刊近似算例的结果”,不把它写成这个三小时轨道在指定站点的已核实最大频移,更不写成某个卫星系统的实测指标。原文是教学讨论,没有提供实测轨迹、频率时间序列或误差统计。

第41页末尾还给出一个合并公式,其中分母含1.8×10的7次方,括号内又含27.9这一常数,同时说明轨道半径单位为米、周期单位为秒。按这些标注与原例参数直接代入,所得约1.70 MHz,无法复现前面的96.9 kHz。这里是按印刷式进行的复算,不是另一个预测值。

固定常数有可能隐含单位换算,但这不能由读者自行猜定。尤其在长度、时间和角度组合的公式中,一个数字从哪里来、带什么单位、对应哪一步投影,都应能追溯。本文没有获得作者勘误,不把某一种猜测改写成作者本意,也不把这条合并式作为可直接套用的通式。

“最大”还需要说明搜索范围。原刊针对圆轨道过顶并看到几何地平线的简化情形;真实站点可能只在高于某个仰角时接收,也可能经历不过顶的过境。可见区间内的极值,应由该区间的距离变化率决定。原文的一张局部示意图不能替代所有轨道、站点和时间的范围定义。

七 接收机看到的频差,还包含几何以外的因素

把几何频移解释清楚之后,还要避免另一种混淆:接收机测得一个频率偏差,就把它全归于卫星运动。ESA Navipedia的GNSS基本观测量说明指出,多普勒主要受天线之间相对运动影响,同时还包含与接收机钟频误差有关的公共偏置。[4] 这说明理想几何量与实际测量值不是同一个层次。

对一般读者,可以把频率基准想象成两把刻度略有差异的尺。即使被测对象的几何关系没有变化,用不同刻度比较,也可能得到偏差。这只是解释频率基准的编辑类比,并不代表所有偏差都能归结为一个固定常数,更不替代具体仪器的校准模型。

时间尺度同样值得留意。公式中的距离导数描述一个瞬间,而若在一段时间内统计载波相位变化,得到的是这段时间对应的平均变化。接近与远离的过程若跨过零值,平均频移可能很小,但这不表示期间每个时刻的频移都很小。这里是在解释导数与区间平均的数学区别;将某条实测记录与模型比较时,还要知道记录对应什么时刻和时段,不能只比较表格里两个孤立的频率数字。

同样,不能把“频移大”和“信号一定更难接收”直接画等号。频率偏离多少、变化有多快,以及接收设备在什么条件下保持同步,是不同问题。本文不由轨道速率推导某台接收机的性能;没有相应设备资料和验证条件,就没有足够证据给出这种判断。

原刊讨论的是星地之间单程信号的基础关系。若换成往返信号或带转发过程的测量,传播链路和频率关系已经改变,不能把单程式原样搬用。阅读公式时先确定“哪一端发、哪一端收、比较哪一个频率”,比看到多普勒三个字就套同一式子更可靠。

八 技术导读最值得留下的是一条可检查的推理

这篇两页专栏提供了一个清楚的教学起点:把地面自转、圆轨道速度、轨道方向和频移联系起来。深读的意义,是继续追问每一步究竟保留了什么条件。地球半径与轨道高度不能混用,分钟与秒不能混用,轨道速率与视线速度也不能混用。

读者若只记住一个判断,可以记住“距离怎样变化”。它既能解释高速运动的理想地球静止卫星为何没有相应的一阶几何频移,也能解释过顶瞬间为何可能经过零值,还能帮助识别地平线近似遗漏的投影。一个定义比一串失去条件的数字更耐用。

本文保留了原刊可复现的基础速率计算,同时指出符号、方位译法、几何投影和合并式单位的边界。新增推导与核验均为编辑说明,不是原作者的勘误声明。读懂一篇技术文章,不仅是把公式算到最后,更是能够说清结果在哪个模型里成立。

资料来源与说明

[1] Dave Adamy. Space EW (Part 10): Doppler Shift in Satellite Link. JED, March 2026, pp.40–41。原刊图1—4、表1与全文均已回核。

[2] ESA. Geostationary orbit. 2020-03-02。用于核对理想地球静止轨道与地球自转的关系。

[3] NASA/JPL Education. Exploring the Doppler Effect With NASA. 更新于2025-09-25。用于核对视线速度与接近、远离时的频移方向。

[4] J. Sanz Subirana, J. M. Juan Zornoza, M. Hernández-Pajares. GNSS Basic Observables. ESA Navipedia, 2011。用于区分多普勒观测中的相对运动和接收机钟频误差。

两幅配图均为编辑绘制;本文几何检查和数值复算不是原刊实测结果,也不代表作者发布了正式勘误。

阅读5
分享
写评论...