翻过《雷达原理》第三版的人多少都有同感:公式量很大,但真正绕不开的,翻来覆去就是雷达方程、热噪声、脉冲重复频率和多普勒频率这么几条线。我前段时间帮朋友做一套雷达总体参数预算,又把整本教材的公式从头顺了一遍。最大的感触是,这些公式不是孤立考点,而是一条物理链条:距离方程决定链路能不能看到目标,噪声决定灵敏度下限,脉冲重复频率决定距离和速度的模糊情况,多普勒频率决定能不能测速和分辨运动目标。这篇就把我认为最核心的那几组公式拆开讲清楚,同时结合工程里常见的换算和误用场景,给正在啃教材,或者需要做雷达链路预算、信号处理算法验证的朋友做个参考。
1. 雷达距离方程:先把“距离四次方”这条主线焊死
1.1 从发射到接收:三个乘积项怎么拼出雷达方程
雷达距离方程最容易背错,是因为它不像通信链路那样只有单向衰减。电磁波从发射天线出去,先走一段距离 (R) 到达目标,被目标反射后还要再走一段距离 (R) 回到接收天线,总共是双程传播。所以接收功率密度不是按 (1/R^2) 衰减,而是按 (1/R^4) 衰减。这在物理上可以拆成三步看:
- 第一步,发射机辐射功率 (P_t) 乘天线增益 (G_t),在距离 (R) 处形成的照射功率密度是:
[ \frac{P_t G_t}{4\pi R^2} ]
第二步,目标相当于一个“再辐射天线”,它截获一部分功率后向空间散射。用雷达散射截面 (\sigma) 来衡量目标反射能力,等效为一个各向同性的辐射面积,单位是平方米。
第三步,散射回来的功率再经过一次 (4\pi R^2) 扩散,被接收天线的有效面积 (A_e) 捕获。
把三个因素乘在一起,就得到教材里最常见的接收功率公式:
[ P_r = \frac{P_t G_t \sigma A_e}{(4\pi)^2 R^4} ]
对于单基地雷达,收发共用天线,通常直接把增益写成一个 (G)。再利用天线增益和有效面积的关系 (G = \frac{4\pi A_e}{\lambda^2}),把 (A_e) 代换掉,就得到更常用的单基地雷达方程:
[ P_r = \frac{P_t G^2 \lambda^2 \sigma}{(4\pi)^3 R^4} ]
这里每个符号的单位都必须是约定好的:功率用瓦,增益用倍数而不是分贝,波长和距离用米,雷达散射截面用平方米。很多初学者喜欢上来就背这个最终形式,结果漏掉中间那个 (4\pi) 的量级,算出来差得很离谱。我的建议是,先记住“照射、散射、接收”这三次转换,再根据需要代入 (A_e) 和 (\lambda) 的关系,这样即使忘了最终形式,现场也能推出来。
1.2 天线增益与有效口径的换算关系
雷达方程里最容易被低估的变量其实是天线增益。教材里给出的基本关系是:
[ G = \frac{4\pi A_e}{\lambda^2} ]
这里的 (A_e) 是有效口径,不是天线物理面积。实际天线存在口径效率 (\eta),一般取 0.5 到 0.7 之间,因此有:
[ A_e = \eta A ]
举个具体例子:X 波段波长 (\lambda = 0.03\text{ m}),如果天线设计增益是 40 dB,也就是线性增益 (10^4),那么有效口径大约是:
[ A_e = \frac{G\lambda^2}{4\pi} \approx \frac{10^4 \times 0.0009}{12.57} \approx 0.72\text{ m}^2 ]
如果口径效率取 0.6,物理口径大约要 (1.2\text{ m}^2)。这个估算在相控阵天线布阵时非常有用:知道增益指标,就能反推单元数量和阵列面积。反过来,如果你先知道天线物理尺寸,也可以用这个公式估算增益,再去代入雷达方程。
1.3 最大作用距离的公式形态
雷达方程最核心的应用是求最大作用距离。做法很简单:当接收功率刚好等于最小可检测信号 (S_{\min}) 时,目标就在最大距离上。于是:
[ R_{\max} = \left[ \frac{P_t G^2 \lambda^2 \sigma}{(4\pi)^3 S_{\min}} \right]^{1/4} ]
这里会出现一个反直觉的结论:由于作用距离与接收功率是四次方根关系,接收灵敏度改善 1 dB,作用距离只改善 0.25 dB,也就是大约 1.06 倍。这说明“往距离上砸钱”是很不划算的。真正立竿见影的是加大天线孔径,因为增益 (G) 在公式里出现了两次,作用距离只对 (G) 开四次方,(2/4=0.5),也就是说增益提升 1 dB,距离提升约 0.5 dB,大约是 12%。相比提高灵敏度,提高增益和孔径的效率要高得多。
当然,(S_{\min}) 在教材后面会被进一步展开成噪声、带宽、信噪比和损耗的表达式,这正是第二节要聊的内容。
2. 噪声决定灵敏度:最小可检测信号那条公式
2.1 热噪声功率:先记住那条物理下限
雷达接收机不管做得再好,热噪声是躲不掉的。接收机输入端的热噪声功率可以写成:
[ N = k T_0 B ]
其中 (k = 1.38 \times 10^{-23}\text{ J/K}) 是玻尔兹曼常数,(T_0) 是参考温度,通常取 290 K,(B) 是接收机带宽,单位是 Hz。常温下,每赫兹带宽的热噪声功率 (kT_0) 大约是 (-174\text{ dBm/Hz})。很多射频工程师口里的“底噪 -174”,指的就是这个值。
如果带宽是 1 MHz,那么理想接收机的噪声底就是:
[ -174 + 10\log_{10}(10^6) = -114\text{ dBm} ]
这个数在链路预算里非常好用。实际接收机还有额外噪声,工程上会用噪声系数 (F) 来统一表示。
2.2 噪声系数与等效噪声温度
噪声系数定义为:
[ F = \frac{S_i / N_i}{S_o / N_o} ]
理想情况下 (F = 1),也就是 0 dB,实际接收机一般大于 1。噪声系数也可以换算成等效噪声温度:
[ T_e = T_0 (F - 1) ]
在相控阵雷达里,系统噪声温度经常包含天线温度、接收机噪声温度和馈线损耗几部分。如果接收机前级有 2 dB 损耗,即便后面的低噪声放大器再好,这 2 dB 损耗也会直接抬高系统噪声温度,作用距离就会掉。这是一条特别容易被忽略的工程经验:馈线损耗和接收机噪声系数同样重要,不能只看低噪声放大器的指标。
2.3 把信噪比写进雷达方程
最小可检测信号不是绝对功率,它取决于后续检测需要的最低信噪比。接收到的目标回波功率必须大于噪声乘以最低信噪比,才能可靠检测。所以把 (S_{\min}) 展开成 (S_{\min} = k T_0 F L B (S/N)_{\min}),代入雷达方程后得到:
[ R_{\max}^4 = \frac{P_t G^2 \lambda^2 \sigma}{(4\pi)^3 k T_0 B F L (S/N)_{\min}} ]
这里的 (L) 是总损耗,包括波导损耗、天线罩损耗、雨衰、扫描损耗等,通常用倍数形式进公式。如果这些都以 dB 给出,要先把 dB 换回线性倍数,或者先用 dB 形式逐项加,最后再统一换算。
还有一个常用的变形:如果接收机带宽与脉冲宽度匹配,也就是 (B \approx 1/\tau),那么公式里的 (P_t / B) 约等于 (P_t \tau),也就是一个脉冲的能量 (E)。这样最大距离公式可以写成脉冲能量形式:
[ R_{\max}^4 = \frac{E G^2 \lambda^2 \sigma}{(4\pi)^3 k T_0 F L (S/N)_{\min}} ]
这就是为什么很多人说“雷达测距靠能量,不靠瞬时功率”的原因。峰值功率再高,如果脉宽很窄、能量不够,照样看不远。相反,用脉冲压缩把大时间带宽积信号的能量压到场后再检测,效果等同于提高接收信噪比。
做一次手算校核会很有感觉。比如某雷达峰值功率 1 MW,天线增益 40 dB,波长 0.03 m,目标散射截面 1 平方米,带宽 1 MHz,噪声系数 3 dB,总损耗 3 dB,最低信噪比要求 13 dB。代入公式后估算出来的最大作用距离大约在 100 公里量级。如果目标散射截面降到 0.1 平方米,作用距离会掉到原来的 (10^{-1/4} \approx 0.56),只剩 50 多公里。这说明雷达作用距离对目标截面的变化非常敏感。
3. 脉冲体制的尺子:距离分辨率、不模糊距离和盲区
3.1 距离分辨率的两种写法
脉冲雷达最基本的一个公式是距离分辨率和脉冲宽度之间的关系:
[ \Delta R = \frac{c \tau}{2} ]
因为电磁波是双程传播,所以时间要除以 2。比如一个 1 微秒的脉冲,对应的距离单元宽度大约是 150 米。想要 3 米的分辨率,脉冲宽度就得压到 20 纳秒。但这么短的脉冲,能量又不够,于是脉冲压缩技术就登场了。
现代雷达普遍使用调频或编码信号,距离分辨率不再由脉冲宽度决定,而是由信号带宽决定:
[ \Delta R = \frac{c}{2B} ]
这里 (B) 是调制带宽。比如带宽 1 MHz 时距离分辨率是 150 米,带宽 500 MHz 时距离分辨率可以做到 0.3 米。脉冲压缩还额外带来一个处理增益:
[ G_{\text{pc}} = B \tau ]
也就是时间带宽积。简单脉冲的时间带宽积约等于 1,所以没有额外增益;线性调频信号的时间带宽积可以做到几百甚至几千,匹配滤波后信噪比收益非常可观。
3.2 最大不模糊距离与脉冲重复频率
脉冲雷达另一个绕不开的公式是最大不模糊距离:
[ R_u = \frac{c T_r}{2} = \frac{c}{2 f_r} ]
这里 (T_r) 是脉冲重复周期,(f_r) 是脉冲重复频率。如果目标距离超过 (R_u),来自第 (n) 个发射脉冲的回波会落到第 (n+1) 个接收周期里,看起来像近处目标,这就是距离混叠或距离模糊。
举个例子,脉冲重复频率 1000 Hz 时,最大不模糊距离是 150 公里。对于地波雷达这种需要看超远程目标的体制,重复频率往往压得很低,这时候距离模糊几乎是必然的,需要用多重脉冲重复频率去解模糊。教材里解模糊的数学本质相当于中国剩余定理,但你能在工程里一直用的基础公式,始终是这一个 (R_u = c/(2f_r))。
3.3 盲区和占空比:距离可能被“挖掉”
雷达发射脉冲时,接收机为了保护自己,通常会被关闭一段时间。这段时间在地面上对应的盲区距离是:
[ R_{\text{blind}} = \frac{c \tau}{2} ]
注意这个公式和距离分辨率公式长得一样,但含义完全不同。距离分辨率说的是两个目标能不能分开,盲区说的是这段时间内回波根本收不到。脉宽越宽,盲区越大,但发射能量越多,测程越远。这就是每部雷达都要面对的矛盾。
引入脉冲重复频率后,还可以定义占空比:
[ D = \frac{\tau}{T_r} = \tau f_r ]
平均发射功率等于峰值功率乘占空比:
[ P_{\text{av}} = P_t D = P_t \tau f_r ]
设计雷达时要关注这个量,因为发射机的散热、电源容量和平均功率都受它限制。大功率雷达的占空比通常只有百分之几,盲目增加脉宽去补能量,可能先撞上盲区变大和平均功率受限这两堵墙。
4. 多普勒公式:从径向速度到频移
4.1 基础表达式:多普勒频率等于两倍径向速度除以波长
雷达测速的起点是经典多普勒公式。目标径向速度为 (v_r),电磁波波长为 (\lambda),基带多普勒频率为:
[ f_d = \frac{2 v_r}{\lambda} ]
也可以写成:
[ f_d = \frac{2 v_r}{c} f_0 ]
其中 (f_0) 是发射载频。这里“径向速度”是指目标相对雷达连线方向的速度分量,不是目标的绝对运动速度。目标靠近雷达时多普勒频率为正,远离时为负,这个正负符号在具体处理中要看系统定义,但幅度计算都是一样的。
教材里还会给一个更严格的表达式,考虑相对论效应时接收频率为:
[ f_r = \frac{1 + v/c}{1 - v/c} f_0 ]
在 (v \ll c) 时近似展开,就得到上面那个线性公式。工程上等不到相对论修正起作用,直接用 (f_d = 2v_r/\lambda) 完全足够。
4.2 慢时间采样里的相位变化:脉冲多普勒的核心
如果用脉冲雷达去测多普勒,每个脉冲之间的相位差是关键。同一距离单元,目标在每个脉冲重复周期内移动一小段距离,回波相位随之变化。相邻两个脉冲的相位差为:
[ \Delta \phi = 2\pi f_d T_r ]
其中 (T_r) 是脉冲重复周期。沿“慢时间”维度看,同一距离门的回波相位在做等速旋转,旋转频率就是多普勒频率。对慢时间信号做 FFT,就能把多普勒频率分离出来。这就是脉冲多普勒雷达和 MTD 处理的物理基础。
MTI 对消器同样建立在相位变化基础上。最简单的两脉冲对消器传递函数可以写成:
[ H(z) = 1 - z^{-1} ]
它的频率响应在直流处为零,在 (f_d = n f_r) 处也出现凹口,所以固定目标能被抑制,但速度恰好多普勒频率等于重复频率整数倍时,运动目标也会被对消掉,这就引出了盲速问题。
4.3 多普勒模糊和盲速:低重复频率躲不过的坑
脉冲雷达对慢时间信号采样,采样率就是脉冲重复频率 (f_r)。因此多普勒频率能无模糊测量的范围只有:
[ |f_d| \le \frac{f_r}{2} ]
换算成速度,单边最大不模糊径向速度:
[ |v_{\max}| \le \frac{\lambda f_r}{4} ]
这里容易把人绕晕的是“盲速”公式。盲速是多普勒频率恰好落在脉冲重复频率整数倍时的速度:
[ v_{\text{blind}, n} = \frac{n \lambda f_r}{2} ]
举个例子,波长 0.03 米,脉冲重复频率 1000 Hz,那么 (n=1) 时的盲速是 15 米每秒。也就是说,以 15 米每秒径向速度运动的目标,在慢时间采样后看起来和静止目标没有区别,会被 MTI 对消器一起滤掉。所以工程上经常用参差重复频率或高重复频率波形来改变盲速位置,让目标不会同时落在所有波形的凹口里。
这里还要注意最大不模糊测速和最大不模糊测距是矛盾的。要提高测速范围,就要提高重复频率,但重复频率一高,最大不模糊距离就变短。雷达总体设计经常在这种矛盾里反复权衡,而教材里的多频工作方式就是用来同时解距离模糊和速度模糊的。
5. 公式落到预算和仿真时,我常踩的几个坑
5.1 单位不统一是最大的错误来源
雷达方程里如果用米、瓦、倍数是比较省心的。可实际工程里大家习惯说公里、兆瓦、dB,于是换算就成了第一道坎。
功率要统一:(1\text{ MW} = 10^6\text{ W}),(1\text{ kW} = 10^3\text{ W})。如果混用,方程里差出 3 个数量级,作用距离会差接近 10 倍。
距离要统一:(R) 必须用米,不能拿着 150 km 直接代进去。带宽要用 Hz,不能用 MHz。波长要用米,不能用毫米。雷达散射截面单位是平方米,可很多数据表里会给 dBsm,这时先换算回线性值再代入,不然整个链路预算都是错的。
5.2 把 dB 和倍数混着算
公式里的 (G)、(F)、(L) 在教科书里默认是线性倍数。但工程图纸上写的是 40 dB、3 dB、5 dB 这类值。如果你把 40 当作 40 代进去,而不是换成 (10^4),最后结果会差很多个数量级。
反过来,如果做链路预算,习惯上把整个方程都写成 dB 形式,用 10lg 和 20lg 逐步累加。这本身没问题,但只要中间混了一步“直接把 dB 当倍数乘”,后面全错。我的习惯是:要么全部用 dB 形式从头加到尾,要么全部换回线性值一次算完,不要来回切换。
5.3 漏掉损耗项,或者把损耗重复计入
雷达方程里那个 (L) 是一个非常“现实”的参数。天线罩损耗、波导损耗、扫描损耗、雨衰、信号处理损耗,每一笔都得记。很多初算结果偏乐观,原因不是公式用错,而是把 (L) 整个漏掉了。
但画蛇添足的也常见:如果 (F) 已经包含了接收链路的损耗贡献,又在 (L) 里重复加一遍馈线损耗,灵敏度就会被算低,预算又偏悲观。正确做法是把损耗按位置分清楚:天线罩和扫描损耗放在发射/接收链路里,接收机馈线损耗放在噪声系数换算里,还是单独放在系统损耗里,必须在预算表里写清楚口径,前后保持一致。
5.4 仿真验证时别忘了信噪比门限
做 matlab 仿真时,很多人直接把公式里最小的可检测信号换成某个任意门限,结果“发现作用距离特别远”。其实最小可检测信号应该对应一组检测概率和虚警概率,不是拍拍脑袋给个 dB 数。
教材里对检测概率、虚警概率和脉冲积累的分析,本质上就是在告诉你那个 ((S/N)_{\min}) 是怎么来的。我在实际仿真中会先用单脉冲信噪比做检测门限粗算,再加入相参积累增益。相参积累 (n) 个脉冲,理想情况下信噪比可以提高 (n) 倍:
[ \left(\frac{S}{N}\right){\text{积累后}} = n \cdot \left(\frac{S}{N}\right){\text{单脉冲}} ]
非相参积累达不到这个收益,通常还要额外扣除积累损失。很多初学者用一个很高的相参积累增益,作用距离算出来特别乐观,就是因为忘了积累过程本身有杂波抑制损失、通道失配和量化损失。做系统预算时,把这些损失写进总的 (L) 里,结果才会贴近实测。
5.5 波形参数不是孤立选的
最后想说一个很实际的体会:雷达方程、距离分辨率和多普勒盲速这三组公式,在教材里分开讲,在工程里却是绑定在一起的。脉冲宽度 (\tau) 同时影响盲区、距离分辨率和发射能量;脉冲重复频率 (f_r) 同时影响最大不模糊距离、多普勒不模糊范围和盲速位置;接收机带宽 (B) 又和脉宽、脉压增益直接相关。哪怕只是估一部雷达的“能不能看到目标”,也要先把这几个参数摆在一起权衡,而不是只看某一个公式。
我自己复习《雷达原理》第三版时,有一个笨办法:把每个公式里的每个变量都追问一遍“改了它,哪些指标会好,哪些指标会坏”。雷达方程里增益提上去,距离变远,但孔径和成本变大;重复频率提上去,盲速变大,但距离模糊变严重;带宽提上去,分辨率变好,但噪声功率也跟着变大,灵敏度又会下降。把这些权衡关系想明白,公式就不再是死记硬背的东西了,而是可以用来和同事讨论波形设计、解释仿真结果、甚至判断一部雷达到底行不行的一组趁手工具。