Skip to content

第1章 雷达方程的发展

雷达距离方程是在第二次世界大战期间发展起来的,用于分析雷达系统性能,并指导雷达研发人员在当时有限的设计选项中进行选择。该主题的最早文献曾受军事安全限制,但在战争结束后发表,并被广泛传播。

预测雷达探测距离的基本方法与早期研究保持一致,这一研究首次发表于1947年美国海军研究实验室的Norton和Omberg的论文中[1,2]。我们将该方程称为原始雷达方程。在本章中,我们回顾雷达方程的发展步骤,并讨论其演化为可应用于现代雷达系统分析与设计的形式。为了避免不同作者使用的符号差异造成混淆,我们将引用文献中的符号替换为一致的符号系统,这些符号在使用时定义,并在本书末尾的附录A中列出。这使得能够直接比较这些方程及其局限性和对当前问题的适用性。

1.1 雷达方程基础

雷达方程的目标是计算在给定的雷达、目标和环境参数条件下,能够实现期望探测性能的最大距离 \(R_{m}\)。本章讨论的雷达方程限定在热噪声为唯一干扰源的环境下,即目标回波信号必须在这种干扰下被成功检测。针对其他环境的方程将在第3章中推导。

雷达距离方程的推导分为三步:

  1. 表达在给定参数下可获得的最大信噪比随距离的变化;
  2. 表达满足探测需求所需的最小信噪比;
  3. 结合这两个表达式,求解满足要求的雷达可实现的最大探测距离。

1.1.1 最大可用信噪比

在二战期间,North在机密研究中[3](后发表于《IEEE会刊》)指出,当接收系统对发射波形使用匹配滤波器时,可以获得最大可能的信噪功率比 \((S / N)_{\max}\)。这个最大比值等于波形的能量比 \(E / N_{0}\) \({ }^{1}\),其中 \(E\) 是回波信号能量,\(N_{0}\) 是竞争热噪声的功率谱密度。第1步所要推导的表达式就是在给定系统参数下的 \(E / N_{0}\)。在[1]中,\(E / N_{0}\) 定义为单脉冲的可用能量比,参照点是接收天线输出端口。本书始终采用这一参考点。

Note

\({ }^{1}\) 注意,本书中的能量比是 \(E / N_{0}\),而不是某些教材中使用的 \(\Re=2E/N_{0}\)

从各向同性发射天线出发,在任意距离 \(R\) 处,发射脉冲的能量密度 \(E_{p}\) 为:

\[ \begin{equation*} E_{p}=\frac{E_{t}}{4 \pi R^{2}} \quad\left(\mathrm{~J} / \mathrm{m}^{2}\right) \tag{1.1} \end{equation*} \]

其中 \(E_{t}\)(焦耳)是发射脉冲能量,\(4 \pi R^{2}\) 是以雷达为中心半径为 \(R\) 的球面面积。假设发射矩形脉冲,脉宽为 \(\tau(\mathrm{s})\),峰值功率为 \(P_{t}(\mathrm{W})\),则 \(E_{t}=P_{t}\tau(\mathrm{J})\)。当发射天线具有增益 \(G_{t}\) 时,波束轴上的能量密度被该增益放大:

\[ \begin{equation*} E_{p}=\frac{E_{t} G_{t}}{4 \pi R^{2}}=\frac{P_{t} \tau G_{t}}{4 \pi R^{2}} \quad\left(\mathrm{~J} / \mathrm{m}^{2}\right) \tag{1.2} \end{equation*} \]

经过雷达散射截面为 \(\sigma(\mathrm{m}^{2})\) 的目标反射后,入射到接收天线的回波能量密度 \(E_{a}\) 为:

\[ \begin{equation*} E_{a}=\frac{E_{p}\sigma}{4 \pi R^{2}}=\frac{P_{t}\tau\sigma}{(4 \pi R^{2})^{2}} \quad\left(\mathrm{~J} / \mathrm{m}^{2}\right) \tag{1.3} \end{equation*} \]

式(1.3)源于雷达散射截面的定义[4]。

具有有效孔径 \(A_{r}\) 的接收天线接收到的能量 \(E\) 为:

\[ \begin{equation*} E=E_{a} A_{r}=\frac{P_{t}\tau G_{t} A_{r}\sigma}{(4 \pi R^{2})^{2}} \quad(\mathrm{~J}) \tag{1.4} \end{equation*} \]

接收天线增益的表达式为:

\[ \begin{equation*} G_{r}=\frac{4 \pi A_{r}}{\lambda^{2}} \tag{1.5} \end{equation*} \]

其中 \(\lambda\) 为波长,由此得到脉冲雷达接收天线输出端口处的最大可用信号能量:

\[ \begin{equation*} E=\frac{P_{t}\tau G_{t} G_{r} \lambda^{2} \sigma}{(4 \pi)^{3} R^{4}} \quad(\mathrm{~J}) \tag{1.6} \end{equation*} \]

到目前为止,我们考虑的是理想情况,即发射机输出直接送入天线,且发射机—目标—接收机路径上不存在任何减少接收信号能量的损耗。为了考虑这些射频损耗,可以在分母中引入因子 \(L_{1}\),得到:

\[ \begin{equation*} E=\frac{P_{t}\tau G_{t} G_{r} \lambda^{2} \sigma}{(4 \pi)^{3} R^{4} L_{1}} \quad(\mathrm{~J}) \tag{1.7} \end{equation*} \]

式(1.7)的推导遵循现代雷达教材中的方法,但不同之处在于它是用脉冲的能量而不是功率来表述发射与接收信号。

为了求出可用能量比,参考直接连接到天线输出端口的接收机,噪声功率谱密度 \(N_{0}(\mathrm{W}/\mathrm{Hz})\) 表示为:

\[ \begin{equation*} N_{0}=kT_{s} \quad (\mathrm{W}/\mathrm{Hz} \text{ 或 J}) \tag{1.8} \end{equation*} \]

其中
\(k=1.38 \times 10^{-23} \mathrm{J}/\mathrm{K}\) 为玻尔兹曼常数;
\(T_{s}=\) 系统噪声温度(开尔文)(见第6章)。

结合式(1.7)和(1.8),得到最大可用信噪功率比:

\[ \begin{equation*} \left(\frac{S}{N}\right)_{\max}=\frac{E}{N_{0}}=\frac{P_{t}\tau G_{t} G_{r} \lambda^{2}\sigma}{(4 \pi)^{3} R^{4} kT_{s} L_{1}} \tag{1.9} \end{equation*} \]

1.1.2 最小所需信噪比

雷达方程推导的第二个表达式给出了获得指定探测性能所需的能量比。在[1]中假设探测由人工操作员完成,其通过阴极射线管显示器观察 \(n\) 个连续回波脉冲:

\[ \begin{equation*} V(n)=\frac{E_{\min}}{N_{0}} \tag{1.10} \end{equation*} \]

其中 \(E_{\min}\) 是每个接收脉冲的最小能量,在最佳观测条件下使 \(n\) 脉冲组在显示器上可见。这与当前定义[4]一致:

可见度因子(脉冲雷达)

单脉冲信号能量与单位带宽噪声功率之比,在最佳带宽和观察环境下,能在显示器上达到规定的探测概率和虚警概率。测量位置为接收机中频部分。

因此最小所需能量为:

\[ \begin{equation*} E_{\min}=N_{0} V(n)=kT_{s}V(n) \quad(\mathrm{J})\tag{1.11} \end{equation*} \]

脉冲数 \(n\) 由脉冲重复频率 \(f_{r}\) 与观测时间 \(t_{o}\) 的乘积决定,而 \(t_{o}\) 取雷达波束在目标上的驻留时间或显示器与人眼观察时间常数中的较小值:

\[ \begin{equation*} n=f_{r} t_{o} \tag{1.12} \end{equation*} \]

因此,通过对宽波束或慢扫描条件下获得的连续脉冲进行积分,可以降低所需的能量比。

在[5]中报道了在最佳A型示波器\({ }^{2}\)\(V(n)\) 的实验值。对于PPI显示\({ }^{3}\)的类似结果在[6-12]中给出。对于由电子检测取代人工目视检测的雷达系统,可见度因子被可探测性因子 \(D(n)\) 所取代,其定义为[4]:

被可探测性因子 \(D(n)\)

在脉冲雷达中,单脉冲信号能量与单位带宽噪声功率之比,在给定虚警概率条件下,能实现规定的探测概率。测量方法为在中频放大器带宽内,使用与单脉冲匹配的中频滤波器,并通过最佳视频积分实现。

有时在 \(V\)\(D\) 的下标中加注目标模型和相应的探测概率(例如 \(V_{0(50)}\),表示在稳态目标(Case 0)下实现 \(50\%\) 探测概率所需的能量比[6-12])。通常,\(V\) 大于 \(D\),这是因为人工目视检测效率较低,如第4.6节所述。

Note

\({ }^{2}\) A型示波器显示接收机输出电压随距离的变化。
\({ }^{3}\) PPI显示将接收机电压作为亮度强度绘制在以距离和方位角为坐标的极坐标图上。

1.1.3 脉冲雷达的最大探测距离

将式(1.9)中的可用 \(E / N_{0}\) 与所需值 \(V\)\(D\) 相等时,可得到脉冲雷达的最大距离方程:

对于人工目视探测: $$ R_{m}=\left[\frac{P_{t}\tau G_{t} G_{r} \lambda^{2} \sigma}{(4\pi)^{3} k T_{s} V(n) L_{1}}\right]^{1/4} \quad(\mathrm{m}) \tag{1.13} $$

对于电子探测: $$ R_{m}=\left[\frac{P_{t}\tau G_{t} G_{r} \lambda^{2} \sigma}{(4\pi)^{3} k T_{s} D(n) L_{1}}\right]^{1/4} \quad(\mathrm{m}) \tag{1.14} $$

这些方程仅是雷达探测距离预测的起点。它们给出了在自由空间条件下,沿着共用收发天线波束轴线方向,当接收能量比满足要求时的探测距离。其中,人工目视探测使用最佳滤波器,电子探测使用匹配滤波器。


1.2 原始雷达方程

Norton和Omberg [1] 推导的距离方程由式(1.13)表示,其基本假设如下:

  • 噪声功率谱密度:按下式计算
    $$ N_{0}=k T F_{n} \quad(\mathrm{W}/\mathrm{Hz}) \tag{1.15} $$ 其中 \(T=300 \,\mathrm{K}\)(室温),\(F_{n}\) 为接收机噪声系数。

  • 可见度因子:其推导基于在A型示波器上“勉强可见”的脉冲[5],并定义为包含非最佳带宽效应,而不是像[4]中那样使用理想可见度因子,再引入带宽修正因子[6-12]。

  • 传播效应:文献[1]讨论了对流层吸收、折射及地面反射对雷达—目标路径的影响,但这些效应并未作为项写入原始雷达方程。

  • 射频损耗:在[1]的最终方程中基于式(1.13)引入了收发线损 \(L_{t}\)\(L_{r}\)。然而,该形式因包含多个数值换算常数而较为复杂,已不再实用。

在[1]中还引入了数值换算因子,以便直接给出以英里为单位的距离。此外,通过替换常数 \(4\pi, k, T\) 以及采用可见度因子与目标截面的近似值,最终表达式被一个数值常数主导,从而掩盖了推导中的基本关系。这种做法是为了在当时计算工具有限(计算尺或机电计算器)的情况下,尽量将因素合并为单一数值。然而在今天,这既没有必要,也不可取。为了清晰呈现各关系并便于比较不同形式的雷达方程,后续讨论保留推导中的各项,而不引入数值换算因子。

尽管存在这些限制,二战期间提出的原始雷达方程至今仍具意义。它通过信噪能量比揭示了自由空间路径探测距离与雷达及目标参数之间的基本关系。方程提出者基于North的匹配滤波理论并采用能量比作为探测成功的判据,其洞见值得认可,正如Blake所强调的那样,但却在很多当代文献中被忽视。


1.3 Blake的脉冲雷达方程

美国海军研究实验室的Lamont V. Blake [6-13] 在二战研究成果的基础上,引入了更精确的术语定义。他给出的脉冲雷达基本方程为[12, p. 19, Eq. (1.34)]:

\[ \begin{equation*} R_{m}=\left[\frac{P_{t}\tau G_{t} G_{r} \lambda^{2} \sigma F_{t}^{2} F_{r}^{2}}{(4\pi)^{3} k T_{s} D(n) C_{b} L}\right]^{1/4} \quad(\mathrm{m}) \tag{1.16} \end{equation*} \]

式(1.16)的分子在基本单位下维度为 \(\mathrm{J}\cdot\mathrm{m}^{4}\)。相较于式(1.14),新增加的项为:

\[ \begin{aligned} & F_{t}=\text{发射路径的方向—传播因子;} \\ & F_{r}=\text{接收路径的方向—传播因子(见第8章)。} \end{aligned} \]

分母维度为 J,其中新增项为:

\[ \begin{array}{ll} T_{s} & =\text{系统输入温度,单位K(见第6章);} \\ D(n) & =\text{可探测性因子(定义见1.1.2节;参见第4章);} \\ C_{b} & =\text{带宽修正因子(见10.2.3节);} \\ L & =L_{1} L_{2}=\text{系统损耗因子。} \end{array} \]

系统损耗 \(L\) 由两个部分组成:\(L_{1}\) 为式(1.7)中的射频损耗,降低可用能量比;\(L_{2}\) 为信号处理损耗,使得相对于 \(D(n)\) 需要更高的能量比。


1.3.1 Blake方程中各项的意义

Blake雷达方程中的各项定义旨在确保对有无脉内调制的脉冲雷达、采用非相干积累、在自然环境下运行时均能得到准确结果。以下因素值得注意:

  • 发射能量:使用脉冲能量 \(P_{t}\tau\) 使式(1.16)能直接应用于使用脉冲压缩波形的雷达,此时时宽带宽积 \(\tau B \gg 1\),其中 \(B\) 为发射脉冲频谱带宽。这不同于一些方程仅用 \(P_{t}\) 而分母出现 \(N_{0}B\),后者不适用于有脉内调制的情况。

  • 方向—传播因子:Blake将 \(F_{t}\)\(F_{r}\) 直接引入雷达方程,用以描述天线仰角方向图及地面反射的影响,它们是仰角的函数,基本不随距离变化。其定义为[14, p. 35]:

方向—传播因子

在特定条件下某点的电场幅度与在自由空间条件下,当发射波束指向该点时的电场幅度之比。

在式(1.16)中,该因子不包括天线方位方向图和介质衰减的影响,这些被Blake归入损耗因子 \(L\)\(F\) 因子涵盖了目标仰角处相对于轴向定义的 \(G_{t}, G_{r}\) 的天线增益、地面反射对辐射与接收场的贡献,以及低仰角路径衍射对场的修正。对于阵列天线,还包括非正侧向工作带来的增益下降,第10章中将讨论其计算。

  • 噪声温度:Blake引入系统噪声温度 \(T_{s}\) 来表述内部与外部噪声源的综合效应,取代了式(1.15)中的 \(TF_{n}\)。当接收机噪声系数很小且天线指向冷天空时,这一改进尤为重要。需注意,\(T_{s}\) 参考于接收天线输出端,包含来自天线的热噪声、天线与接收机之间的接收线损 \(L_{r}\),以及接收机自身噪声。

  • 能量比的使用:Blake延续了原始雷达方程中使用能量比作为信号判据的做法。

  • 可探测性因子:探测过程既可能是电子探测,也可能是人工目视探测,但任何偏离理想过程的情况,都需要在损耗因子 \(L_{2}\) 中加入相应的补偿项。

  • 带宽因子\(C_{b}\) 由[12, Eq. (8.14a)]给出:

\[ C_{b}=\frac{B_{n}\tau}{4.8}\left(1+\frac{1.2}{B_{n}\tau}\right)^{2} \tag{1.17} \]

其中 \(B_{n}\) 为接收机噪声带宽。当 \(B_{n}=1.2/\tau\) 时,该因子等于1,对应于“最佳”接收机噪声带宽。其来源与应用将在第10.3.3节讨论,而在电子探测雷达中使用 \(C_{b}\) 引起的误差将在1.5.4节讨论。

  • 损耗因子:式(1.16)中的系统损耗因子 \(L\) 由以下组成:
  • \(L_{t}\):发射机输出至天线端口间的传输线损;
  • \(L_{\alpha}\):双程路径上的大气吸收损耗;
  • \(L_{p}\):扫描雷达的波束形状损耗\({ }^{4}\) [13];
  • \(L_{\chi}\):其他接收与信号处理损耗,未单独包含在 \(T_{s}, C_{b}, L_{p}, L_{\alpha}\) 中。

其中,\(L_{t}\)\(L_{\alpha}\) 属于减少可用能量比的射频损耗 \(L_{1}\),而 \(L_{p}\)\(L_{\chi}\) 属于增加所需能量比的损耗 \(L_{2}\)。损耗因子的处理,以及部分因子与所需探测概率相关的情况,将在第4章与第10章中详细讨论。

Note

\({ }^{4}\) Blake将其称为“天线方向图损耗”,但“波束形状损耗”是原始且当前公认的术语。其将在第5章中详细讨论。

1.3.2 距离的求解方法

由于 \(F_{t}\)\(F_{r}\)\(L_{\alpha}\) 等项的复杂距离依赖性,式(1.16)没有封闭解。Blake 讨论了迭代与图解方法[12, pp. 379-388]。一种简单方法被嵌入在常称为Blake图的“脉冲雷达距离计算工作表”中(见图1.1)。该图通过对分贝值进行加减运算求解方程,用袖珍计算器即可方便实现。

天线高于周围地面的高度和目标仰角需要输入到Blake图中作记录,因为结果仅在这些特定条件下有效。用户随后在部分(1)输入五个参数,用它们计算系统噪声温度 \(T_{s}\)。标注为“来自图表”的参数将在第6章中讨论,相关图表会附带等效方程以取代换算表。剩余的距离因子输入在表的(2)部分,并在(3)部分转化为分贝形式,最后求和得到 \(40 \log R_{0(nmi)}\) 的分贝值。自由空间距离 \(R_{0}\)(海里)通过取反对数求得。方向—传播因子 \(F\) 单独计算(见第8章),再乘以 \(R_{0}\) 得到初始距离估计 \(R^{\prime}\)。该估计值随后进入两步迭代以修正大气衰减。当目标处的大气衰减系数较大时(例如毫米波雷达或有降水的微波雷达),可能需要额外的迭代步骤。

Blake图的一个复杂之处在于输入的距离因子采用了混合单位。例如,雷达频率输入为 \(f_{\mathrm{MHz}}=(c/\lambda)/10^{6}\),而不是雷达方程中的 \(\lambda\)。由此得到的距离方程常数将多个换算因子、常数及其量纲组合为:

\[ \begin{aligned} & 10 \log \left[\frac{(10^{3} \,\mathrm{W}/\mathrm{kW})(10^{-6}\,\mathrm{s}/\mu \mathrm{s})(3\times 10^{8}\,\mathrm{m}/\mathrm{s})^{2}/(10^{6}\,\mathrm{Hz}/\mathrm{MHz})^{2}}{(4\pi)^{3}(1.38\times 10^{-23}\,\mathrm{W}\cdot\mathrm{s}/\mathrm{K})(185{,}200\,\mathrm{m}/100\,\mathrm{nmi})^{4}}\right] \\ & = 10 \log \left[1.292^{4}\,(\mathrm{K}\cdot \mathrm{nmi}^{4}/\mathrm{W}\cdot\mathrm{s}^{3}\cdot \mathrm{m}^{2})\right] \\ & = 4.45\,\mathrm{dB}\,[\mathrm{ref}(\mathrm{K}\cdot \mathrm{nmi}^{4}/\mathrm{W}\cdot \mathrm{s}^{3}\cdot \mathrm{m}^{2})] \end{aligned} \]

需要注意的是,在Blake图中输入的是中值雷达散射截面 \(\sigma_{50}\) 和相应的可见度因子 \(V_{0(50)}\),以得到 \(P_{d}=50\%\) 时的探测距离 \(R_{50}\)。通常指定的雷达散射截面是其平均值,对于Swerling 1型(见第4.3.3节)波动目标,平均值比中值大1.5 dB。

脉冲雷达距离计算工作表

图1.1 Blake的脉冲雷达距离计算工作表[8]。引用的图形在第4、6和7章中重现,换算分贝与比值的表格已由数字计算取代。

在[6-8]中给出了 \(V_{0(50)}\) 及电子探测对应的 \(D_{0(50)}\) 曲线,但在很多情况下需要计算 \(P_{d}\neq 50\%\) 的情况。Blake后续的研究[9-12]提供了不同 \(P_{d}\) 和不同目标模型下的 \(D\) 曲线。这些曲线要求在式(1.16)和Blake图中使用目标截面的平均值而不是中值。


1.3.3 Blake图的优点

Blake图的主要优点在于它能记录计算探测距离所用的输入值及若干中间结果。缺乏这种规范可能导致关键数据未被记录,后续无法追溯。当使用工程计算器或数字计算方法替代Blake的换算表和人工中间步骤时,计算既快速又准确。Blake图的具体形式从[6]中的最初版本发展到[12]中的最终形式,经历了多次修改,但基本方法保持不变。第1.4.2节介绍了一种保持基本单位输入的图表改进形式。

Blake图形式化了求解距离相关衰减的迭代过程,它被广泛接受为标准的距离计算方法。虽然理想情况下应以基本单位输入参数,但在现代条件下分贝与比值之间的换算不再依赖表格,口袋计算器或计算机即可轻松完成。然而,保持图表的基本格式仍很重要,以保证结果的可追溯性。如果仅将参数输入计算机程序直接输出探测距离,而未记录和检查输入与中间结果,则容易出现错误。


1.3.4 Blake的相干雷达方程

Blake将式(1.16)推广至相干雷达[12, p. 20, Eq. (1.35)]:

相干雷达:

\[ R_{m}=\left[\frac{P_{av} t_{f}\tau G_{t} G_{r} \lambda^{2} \sigma F_{t}^{2} F_{r}^{2}}{(4\pi)^{3} k T_{s} D(1) C_{b} L}\right]^{1/4} \quad(\mathrm{m}) \tag{1.18} \]

这里,用观测波形的能量 \(P_{av} t_{f}\) 替代了单脉冲发射能量 \(P_{t}\tau=P_{av} t_{r}\),其中 \(P_{av}\) 为平均功率,\(t_{f}\) 为相干处理时间(CPI),\(t_{r}\) 为脉冲重复周期(PRI)。检测所需能量比 \(D(1)\) 对应于相干积分器的单个输出。该结果是一个通用雷达方程(类似于Hall的方程,见1.4.1节),适用于连续波(CW)雷达、脉冲多普勒雷达以及非相干脉冲雷达。

将Blake图应用于相干雷达时,只需在输入中将 \(P_{t}\) 替换为以kW为单位的平均功率 \(P_{av}\),并将 \(\tau\) 替换为以\(\mu\)s为单位的 \(t_{f}\)。如果这些单位对CW雷达不便,则可替换为瓦特和毫秒,而不会改变方程常数。

式(1.18)假设探测由多普勒滤波器输出完成,且没有进一步的非相干积累。这一限制可通过将 \(D(1)\) 替换为 \(D(n^{\prime})\) 来克服,从而使方程适用于对 \(n^{\prime}\) 个滤波输出进行非相干积累的雷达。此时,可用能量的增加由分子中的比值 \(P_{av} t_{f}/P_{t}\tau >1\) 表示,而单样本所需能量的降低由分母中的比值 \(D(n^{\prime})/D(1)<1\) 表示。这样,来自所有 \(n\) 脉冲的相干与非相干增益均被包含在探测距离的计算中。

Blake的雷达方程也通过Skolnik雷达手册的前两版[15,16]得到传播。它能够在实际环境条件下、针对任何类型的波形与处理方式的雷达,准确预测任意仰角目标的最大探测距离。对其方程进行现代化修改所需的步骤相对较少,见1.6节。


1.3.5 Blake的双基地雷达方程

Blake将其基本雷达方程扩展至双基地系统[12, Eq. (1.38)]:

双基地雷达:

\[ \left(R_{t} R_{r}\right)_{m}^{1/2}=\left[\frac{P_{t}\tau G_{t} G_{r} \lambda^{2} \sigma_{b} F_{t}^{2} F_{r}^{2}}{(4\pi)^{3} k T_{s} D(n) C_{b} L}\right]^{1/4} \quad(\mathrm{m}) \tag{1.19} \]

其中 \(R_{t}\)\(R_{r}\) 分别是发射机到目标和目标到接收机的路径。\(\left(R_{t}R_{r}\right)^{1/2}\) 是双基地系统中两条路径的几何平均值。第1.3.1节的评论同样适用于该方程。

Blake图也可以用于双基地雷达,以得到几何平均距离。此时需分别计算 \(F_{t}\)\(F_{r}\),再得到双程因子 \(F=(F_{t}F_{r})^{1/2}\) 并填入图表第8行。第9–14行的大气衰减迭代过程采用两条路径的单程衰减和,并考虑其仰角(应与两天线高度一起单独记录在表头)。

Blake的方程和图表主要用于目标在雷达波形非模糊距离内的探测。对于超出非模糊距离的探测,结果仍然成立,但可能因发射信号的遮蔽(见1.6.1节与10.1.2节)或使用变化PRI导致的积累损失而产生误差。

1.4 雷达方程的其他形式

1.4.1 Hall的雷达方程

Hall 在文献[17]中提出了脉冲雷达的雷达方程,早于Blake的工作,并澄清了前人文献中的许多问题:

脉冲雷达: $$ R_{m}=\left[\frac{P_{t}\tau G^{2}\lambda^{2}\sigma F^{4}}{(4\pi)^{3} k T_{0} F_{n}\times 1.2 D_{x}(n) L_{t} L_{\alpha}}\right]^{1/4} \quad(\mathrm{m}) \tag{1.20} $$

该式等价于(1.14),但有以下替换:

  • 公用收发天线的增益 \(G^{2}\) 替代 \(G_{t}G_{r}\)
  • 分子中加入 \(F^{4}\) 作为双程方向—传播因子;
  • 系统温度采用 \(T_{0}F_{n}\),其中 \(T_{0}=290\ \mathrm{K}\),表示天线输出端处、标准温度环境下的接收机;
  • \(D L\)\(1.2 D_{x}(n) L_{t} L_{\alpha}\) 替代,其中 \(D_{x}(n)\) 为有效可探测性因子,由以下五项相乘得到:
\[ \begin{aligned} D(n) & = \text{非相干积累 $n$ 个脉冲的基本(理论)可探测性因子;} \\ L_{m} & = \text{非最佳带宽损耗;} \\ L_{p} & = \text{波束形状损耗;} \\ L_{c} & = \text{合并损耗;} \\ L_{o} & = \text{操作员和退化因子。} \end{aligned} \]

\(D_{x}\) 被用于现代雷达的方程(见1.6节)。式(1.20)中的系数1.2表示最佳接收机带宽 \(B_{n}=1.2/\tau\) 下的匹配损耗(假设未调制脉冲)。与非最佳中频带宽相关的附加损耗包含在 \(L_{m}\) 中。

Hall 随后提出的通用雷达方程适用于单基地与双基地、任意波形的雷达系统[18, Eq. (8)]:

通用情况:

\[ \left(R_{t}R_{r}\right)_{m}^{1/2}=\left[\frac{P_{av}t_{o} G_{t} G_{r}\lambda^{2}\sigma F_{t}^{2}F_{r}^{2}}{(4\pi)^{3} k T_{s} D_{x}(1) L_{\alpha t} L_{\alpha r}}\right]^{1/4} \quad(\mathrm{m}) \tag{1.21} \]

该表达式结合了Blake的式(1.18)与(1.19),替换为:

  • 相干积累时间 \(t_{f}\) 假设等于观测时间 \(t_{o}\)
  • 可探测性因子与损耗的乘积 \(D(1)L\)\(D_{x}(1)L_{\alpha t}L_{\alpha r}\) 替代;
  • 分母采用发射与接收路径的一程大气衰减乘积 \(L_{\alpha t}L_{\alpha r}\)

由于Blake的方法已广泛应用并能给出准确结果,Hall的表达式受到的关注较少。本文收录它们是为了展示能量比在雷达方程中的合理应用。


1.4.2 Barton的雷达方程

1964年出版的《雷达系统分析》[19]最初采用的雷达方程使用峰值发射功率与接收机带宽内的噪声功率,而非能量。这一不足在1988年出版的《现代雷达系统分析》[20]中得到修正。新版方程基于Blake和Hall的研究,后者引入了有效可探测性因子 \(D_{x}(n)\),替代了基本因子 \(D(n)\) 及其相关损耗。Barton给出的脉冲雷达方程为[20, Eq. (1.2.25)]:

脉冲雷达: $$ R_{m}=\left[\frac{P_{t}\tau G_{t} G_{r}\lambda^{2}\sigma F^{4}}{(4\pi)^{3} k T_{s} D_{x}(n) L_{t} L_{\alpha}}\right]^{1/4} \quad(\mathrm{m}) \tag{1.22} $$

其中,降低可用能量比的射频损耗 \(L_{1}\) 分解为 \(L_{t}\)\(L_{\alpha}\)。增加所需能量比的损耗 \(L_{2}\) 被包含在 \(D_{x}\) 中,其定义避免了Blake因子 \(C_{b}\) 与电子探测结合使用时引发的误差。

该方程另有两种相干雷达形式:

  • 整个观测时间的相干积累[20, Eq. (1.2.26)]:用观测时间内的能量 \(P_{av}t_{o}\) 替代单脉冲能量 \(P_{t}\tau\),并用单样本的 \(D_{x}(1)\) 替代 \(D_{x}(n)\)
  • CPI内相干积累 + \(t_{o}\) 内非相干积累[20, Eq. (1.2.27)]:用CPI内的能量 \(P_{av}t_{f}\) 替代单脉冲能量 \(P_{t}\tau\),并用 \(D_{x}(n^{\prime})\) 替代 \(D_{x}(n)\),其中 \(n^{\prime}=t_{o}/t_{f}\)

第二种形式被用于修正后的Blake图,其示例如图1.2所示。

脉冲雷达距离计算工作表

图1.2 修正后的Blake图[20, p. 21]。

用户在这一版本的Blake图中的输入均采用基本单位,避免了原始图表中的多重换算因子。此时距离方程常数仅包含 \((4\pi)^{3}\)\(k\) 以及将米换算为千米的因子:

\[ \begin{aligned} \text{常数} &= \left[(4\pi)^{3} \cdot 1.38\times 10^{-23}(\mathrm{J}/\mathrm{K})(1000\,\mathrm{m}/\mathrm{km})^{4}\right]^{-1} \\ &= 3.65 \times 10^{7} \, (\mathrm{K}\cdot \mathrm{km}^{4}/\mathrm{J}\cdot \mathrm{m}^{4}) \\ &= 75.6\,\mathrm{dB}\,[\mathrm{ref}(\mathrm{K}\cdot \mathrm{km}^{4}/\mathrm{J}\cdot \mathrm{m}^{4})] \end{aligned} \]

输入步骤、中间计算、分贝换算以及大气衰减迭代的过程与原始Blake图相同。Blake的 \(C_{b}\) 在此被替换为匹配因子 \(M\),定义见第10.2.3节。

2005年出版的《雷达系统分析与建模》给出了最通用的雷达方程形式[21, Eq. (1.20)]:

通用形式: $$ R_{m}^{4}=\frac{P_{av} t_{f} G_{t} G_{r}\lambda^{2}\sigma F_{p}^{2}F_{t}^{2}F_{r}^{2}}{(4\pi)^{3} k T_{s} D_{x}(n^{\prime}) L_{t} L_{\alpha}} \quad(\mathrm{m})\tag{1.23} $$

这一形式与图1.2中的修正Blake图一致,但额外加入了因子 \(F_{p}^{2}\),用于考虑发射与接收天线间可能的极化失配。该因子通常取为1,但其存在强调了接收天线的极化可能与目标回波的实际极化不一致,即使目标截面 \(\sigma\) 是按发射极化计算并给定的。例如,当发射与接收均采用右旋圆极化时,通常指定的 \(\sigma\) 实际对应左旋圆极化分量的回波,而接收天线对此不敏感(见10.1.1节)。


1.5 距离计算中的常见陷阱

有必要在此对上述及文献中的各种雷达方程形式进行评论,并指出使用中可能引发的误差来源。


1.5.1 系统噪声温度 \(\boldsymbol{T}_{s}\)

原始方程(1.9)中的 \(kTF_{n}\) 本意是表示接收机输入端的噪声功率谱密度 \(N_{0}\)。Blake则严谨地定义了其对应的 \(kT_{s}\) 及其多个组成部分,以提高噪声计算的准确性[12, p. 152]。这一方法在现代低噪声系数系统中尤为重要,尤其当天线指向冷天空时。第6章将对其公式进行讨论,并给出实例,表明若简单假设环境温度 \(T_{0}=290\,\mathrm{K}\),会导致方程结果出现显著误差。

Blake的公式应作为所有雷达方程中热噪声表达的标准替代形式。这一做法已在Skolnik雷达手册的前两版[15,16]以及Barton[20-22]中采用。其他一些雷达方程表达式错误地依赖简化关系 \(N_{0}=kT_{0}F_{n}\),该关系仅在 \(F_{n}>10\ \mathrm{dB}\) 的雷达中尚可接受。而在现代雷达中,它无法准确建模噪声功率谱密度,导致最小所需信号能量计算出现 \(1-2\ \mathrm{dB}\) 的误差。


1.5.2 信噪能量比的使用

原始方程以及Blake和Hall的方程均正确地基于输入信号能量与噪声功率谱密度的比值,而不是某个不明确(且通常无法测量)的带宽内的功率比。这使得能够使用North的匹配滤波理论来表示最大可用信噪比[3]:

\[ \left(\frac{S}{N}\right)_{\max}=\frac{E}{N_{0}} \tag{1.24} \]

其中 \(E\) 为信号能量,\(N_{0}\) 为滤波器输入处的噪声谱密度。该关系给出了雷达性能的上限,并可通过引入匹配因子来适应实际接收机与处理器性能低于匹配滤波器的情况。所需输入能量比由此因子提高。

若雷达方程分子采用发射功率 \(P_{t}\) 而分母采用接收机带宽 \(B_{n}\) 或某种“有效带宽”,几乎必然导致错误。在发射脉冲压缩或调频连续波(FM-CW)信号时,这种形式更是完全不适用。即使在未调制脉冲情况下,当 \(B_{n}<1/\tau\) 时,分母中的小 \(B_{n}\) 会使计算的 \(R_{m}\) 过大,趋近无穷大(即接收机失效)。当 \(B_{n}\gg 1/\tau\) 时,计算出的 \(R_{m}\) 又可能偏小,因为较小的视频或显示带宽相对于接收机通过的噪声减少了噪声功率。

纵观雷达文献,许多作者在方程分母中直接引入接收机带宽 \(B_{n}\) 的噪声功率。在少数分子中出现脉冲能量的情况下,推导通常基于未调制脉冲且假设 \(B_{n}\tau=1\)


1.5.3 平均功率的使用

采用平均发射功率而非峰值功率的雷达方程(如式(1.18)、(1.21)或(1.23))更值得推荐,因为它们强调了在热噪声环境中探测距离对平均功率的依赖,而不是对峰值功率或其他波形参数。这些参数虽然在确定探测距离之外的应用中具有实际意义,但并非决定探测性能的根本因素。只要对匹配滤波和积分(相干与非相干)的贡献和损耗表述准确,使用峰值功率并非错误。但在现代雷达中,系统带宽已无法直接测量,使得这种表达愈发困难。


1.5.4 带宽修正与匹配因子

当将式(1.17)中的 \(C_{b}\) 应用于电子探测雷达时,会引入误差。在 \(B_{n}\tau=1.2\) 时取 \(C_{b}=1.0\),会错误地暗示此“最佳带宽滤波器”就是匹配滤波器。而原始雷达方程推导自人工在示波管显示器上目视检测的雷达,实验测得的 \(V\) 值中已包含显示器/观察者过程固有的 \(2-3\ \mathrm{dB}\) 损耗,即便在“最佳带宽”下也是如此[19, p. 171]。在式(1.14)或(1.16)中使用的可探测性因子 \(D\) 定义在匹配滤波条件下,因此需要引入匹配因子 \(M\)(不同于 \(C_{b}\)),以表达任何滤波器不匹配效应。该问题在4.6节和10.3节将进一步讨论。


1.5.5 任意目标的可探测性因子

文献及第4章中给出了常用目标模型的可探测性因子精确方程。Rice [23] 推导的稳态目标(Case 0)单脉冲可探测性因子 \(D_{0}(1)\) 可通过数学程序方便求解。Marcum [24] 将理论扩展至非相干积累 \(n\) 个脉冲的 \(D_{0}(n)\)。Swerling [25] 给出了适用于Case 1–4目标的表达式,其幅度统计对应自由度为2或4的卡方分布,并区分幅度慢变或快变情况。本书将其记为 \(D_{1}(n)\ldots D_{4}(n)\)。文献[21]又补充了一个广义模型,给出 \(D_{e}(n,n_{e})\),适用于幅度统计符合自由度为 \(2n_{e}\) 的卡方分布的目标。该模型涵盖Case 0–4,还包括 \(1\leq n_{e}\leq n\) 的目标(包括常见的非整数值)。

第4章将介绍这些模型的精确与近似计算方法。若采用不合适的目标模型,可能造成数dB的误差,特别是当雷达系统利用时间、频率、空间或极化分集时,Swerling模型描述的目标更易出现问题。

1.5.6 方向—传播因子

在雷达方程的分子中引入方向—传播因子 \(F_{t}\)\(F_{r}\),或它们的几何平均值 \(F\),对于任何目标不在天线仰角波束轴上、且波束主瓣未完全离开雷达—目标路径下方地面的情况,都是准确计算探测距离所必需的。大多数教材虽提及这些因子的必要性,但往往将其与雷达方程分开讨论,而非直接嵌入方程中。对于相控阵雷达在非正侧向方向观测目标时,该因子的处理需要特别关注(见第10.2.1节)。


1.5.7 损耗因子

波束形状损耗 \(L_{p}\) 是式(1.16)和(1.18)中损耗因子 \(L\) 的组成部分之一。这假设目标在回波脉冲的整个积累时间内并未始终处于雷达波束轴线上(例如雷达波束扫描过目标时)。第5章表明,除非积累了分布在天线主瓣范围内的多个脉冲,否则 \(L_{p}\) 是所需探测概率 \(P_{d}\) 的函数。

另一个需要谨慎处理的问题是波束形状损耗与嵌入方向—传播因子 \(F\) 的天线方向图之间的关系。Blake的方法将 \(L_{p}\) 应用于常规二维搜索雷达的方位扫描效应,同时用 \(F^{4}\) 描述目标在仰角方向图中的位置效应。对于二维扫描,可取 \(F=1\),同时在没有地面反射效应的区域引入损耗 \(L_{p}^{2}\),以评估整个仰角扇区的平均性能。该问题将在第8章进一步讨论。

在相控阵雷达中,波束轴向增益随非正侧向扫描角度变化。增益 \(G_{t}\)\(G_{r}\) 通常定义在正侧向波束轴线上,非正侧向扫描因子被包含在 \(F\) 中,因此 \(F\) 同时是仰角与方位的函数。另一种方法是定义整个扫描扇区的平均扫描损耗(见10.2.1节),用于描述扫描范围内的增益变化。与波束形状损耗类似,扫描扇区损耗会增加正侧向所需的能量比,其大小依赖于所要求的 \(P_{d}\),除非通过调整各扫描角度下的发射能量来补偿非正侧向波束的增益降低(见10.2.1节)。

方程中的其他损耗 \(L_{x}\) 用于涵盖一系列损耗因子,这些将在第10章中详细讨论。其中一些损耗因子同样是所需 \(P_{d}\) 的函数,因此若要获得准确的探测距离结果,必须与可探测性因子 \(D\) 一并调整。


1.5.8 探测距离计算常见错误总结

在使用雷达方程时常见的错误包括:

  • \(T_{0}=290\,\mathrm{K}\) 替代实际系统噪声温度 \(T_{s}\)
  • 使用信噪功率比并引入不明确的带宽 \(B\),而不是能量比。任何在方程中包含接收机带宽 \(B\) 的雷达方程都会带来混淆,甚至严重错误;
  • 使用峰值功率但未包含相应的处理因子和损耗;
  • 将带宽修正因子 \(C_{b}\) 用于非包含目视检测损耗的情况;
  • 使用的可见度因子或可探测性因子未考虑目标的实际统计特性和雷达分集效应;
  • 未包含合适的方向—传播因子;
  • 未计入所有实际损耗。

1.6 现代雷达系统的雷达方程

1.6.1 需要修改距离方程的因素

雷达技术与反雷达对抗的发展在应用雷达方程时引入了新的问题。这些问题大多可通过在雷达方程分子中引入合适的距离相关响应因子 \(F_{\text{rdr}}\) 来解决。\(F_{\text{rdr}}\) 包含以下几个组成部分。在低可探测目标的情况下,这些因子尤为重要,因为其减小的散射截面可能使得目标在出现遮蔽效应和STC作用的距离之外无法被探测。


1.6.1.1 遮蔽效应 (Eclipsing)

接收信号的遮蔽效应出现在两类波形中:(1)采用固态发射机、占空比超过约 \(1\%\) 的低PRF雷达;(2)中高PRF雷达。\({ }^{5}\)

Note

\({ }^{5}\) 对于低PRF波形,目标在非模糊距离内;对于高PRF波形,目标在非模糊速度内;对于中PRF波形,目标同时超出非模糊距离和速度。

  • 固态发射机:固态射频功率放大器无法像电子管那样以低占空比换取高峰值功率\({ }^{6}\)。因此,固态发射机的设计倾向于更高的占空比,通常 \(D_{u}=5\%\sim 20\%\),以便在器件峰值功率较低的情况下获得所需的平均功率。对于低PRF雷达,脉冲重复周期 \(t_{r}\) 必须足够长以避免距离模糊,而更高的 \(D_{u}\) 需要比电子管发射机常用的更长脉冲。例如,非模糊距离 \(R_{u}=450\,\mathrm{km}\) 需要 \(t_{r}\geq 3\,\mathrm{ms}\)。若 \(D_{u}=20\%\),则发射脉宽 \(\tau \geq 600\,\mu\mathrm{s}\)。此时,来自 \(R<R_{\min}=\tau c/2 \approx 90\,\mathrm{km}\) 的回波会被遮蔽:回波前段与发射脉冲重叠,导致信噪比损失;来自 \(R>R_{u}-R_{\min}\) 的回波同样会被尾端遮蔽。此类情况的雷达方程需在 \(F_{\text{rdr}}\) 中加入距离相关的遮蔽因子 \(F_{\text{ecl}}\)

Note

\({ }^{6}\) 占空比(或占空因子)是脉冲宽度与脉冲重复周期之比。

部分固态雷达采用双脉冲组形式,即每个长脉冲伴随一个频率偏移的短脉冲,用于覆盖长脉冲遮蔽的近程区。在这种情况下,只有长脉冲的平均功率贡献于远程目标的探测,雷达方程应体现这一点。近程脉冲无法避免非模糊距离远端的遮蔽。

对于无法在 \(R_{\min}\) 内探测低可探测目标的雷达,遮蔽因子可能完全阻止目标在该范围内被探测。

  • 中高PRF波形:MPRF与HPRF波形长期用于机载雷达,也出现在部分地基雷达中。这类波形隐含的要求是探测 \(R>R_{u}\) 的目标。最大探测距离的计算需在方程中引入遮蔽因子,该因子在特定PRF下是距离的确定函数。当在每个波束位置使用多种PRF(无论是在波束驻留期间还是扫描间)时,该因子通常以统计方式表示(见10.2.1节)。

1.6.1.2 灵敏度时间控制 (STC)

许多低PRF雷达使用STC,以避免在中、短程对大目标或杂波导致接收机或处理器饱和,并抑制对小型移动目标(如鸟类和车辆)的不必要探测。STC通过在某个STC范围 \(R_{\text{stc}}\) 内降低接收灵敏度来实现,该范围的选择使目标雷达散射截面阈值 \(\sigma_{\text{min}}\) 得以设定,只有超过此阈值的目标才会被探测。衰减作用施加在接收机之前的RF部分,或接收机的早期级,以限制预期回波输入的幅度。图1.5展示了与距离相关的STC因子 \(F_{\text{stc}}\) 的影响。采用STC的雷达即使在近程,也会失去对低于 \(\sigma_{\text{min}}\) 的低可探测目标的探测能力。


1.6.1.3 波束驻留因子

在某些情况下,雷达波束移动过快,接收波束轴线已不再靠近发射脉冲照射长程目标的角度,而回波在延迟 \(t_{d}<t_{r}\) 时到达。

  • 如果雷达连续以角速率 \(\omega\) 扫描,使得两脉冲之间的角度变化 \(\Delta\theta=\omega t_{r}\) 显著但小于波束宽度,则需要引入与距离相关的波束驻留因子。如果 \(\Delta\theta\) 大于波束宽度,回波信号将丢失。
  • 电子扫描雷达通常分步扫描,波束在发射位置驻留,直至最远目标的回波到达。在某些情况下,接收波束可能会返回到先前脉冲照射目标的角度以接收回波,但这种操作方式不常见,且返回驻留的时间有限。

在这两种情况下,扫描雷达都会对超出波束驻留时间延迟的信号产生歧视,因此雷达方程中引入波束驻留因子 \(F_{\text{bd}}\) 作为 \(F_{\text{rdr}}\) 的组成部分来描述这一效应。


1.6.1.4 频率捷变或分集

脉冲到脉冲的频率变化(捷变)或脉冲组到组的频率变化(分集)用于规避干扰,并平均化目标RCS(例如降低闪烁损耗)。在这些技术中,接收机频率会调整到紧前一发射脉冲的频率。其对RCS平均化的有利作用体现在目标模型中,用于计算可探测性因子 \(D_{e}(n,n_{e})\)。而对于在频率改变后到达的回波,则需要在雷达方程中引入频率分集因子 \(F_{\text{fd}}\) 作为 \(F_{\text{rdr}}\) 的组成部分,以反映该类信号的抑制。


1.6.1.5 透镜因子

Weil [26] 引入的透镜损耗来源于对流层折射随仰角的变化。尽管Blake在[12, pp. 188-192]中讨论过,但该因素有时被忽略,因为它未在雷达方程中作为单独项出现。为避免遗漏,这里将双程透镜因子 \(F_{\text{lens2}}\) 纳入 \(F_{\text{rdr}}\)。需注意它不能与分母中的大气衰减 \(L_{\alpha}\) 合并,因为它并非增加系统噪声温度的耗散性损耗。


1.6.2 适用于现代雷达的方程

Modern Radar System Analysis Software, Version 3.0 (MRSAS3) [22] 是一套运行在 Mathcad \({}^{\text{©}}\) 下的完整工作表(通常被称为程序),可计算雷达探测距离与测量精度。该版本开发于雷达技术变革与目标RCS缩小日益显著的时期,前文1.6.1已指出传统雷达方程的局限性。该软件通过在式(1.23)中增加随距离变化的响应因子 \(F_{\text{rdr}}\) 来改进,系统的电压响应平方修正了随距离变化的可用能量:

\[ E=\frac{P_{av} t_{f} G_{t} G_{r} \lambda^{2} \sigma F_{p}^{2} F_{t}^{2} F_{r}^{2} F_{\text{rdr}}^{2}}{(4\pi)^{3} R^{4} L_{t} L_{\alpha}} \quad(\mathrm{J})\tag{1.25} \]

其中:

\[ \begin{aligned} F_{\text{rdr}} &= F_{\text{ecl}} F_{\text{stc}} F_{\text{bd}} F_{\text{fd}} F_{\text{lens2}} \\ F_{\text{ecl}} &= \text{遮蔽因子;} \\ F_{\text{stc}} &= \text{灵敏度时间控制因子;} \\ F_{\text{bd}} &= \text{波束驻留因子;} \\ F_{\text{fd}} &= \text{频率分集因子;} \\ F_{\text{lens2}} &= \text{双程透镜因子。} \end{aligned} \]

这五个因子在第7章和第10章中详细讨论,这里简要说明:

  • 遮蔽因子:长期以来在高/中PRF机载雷达中被认为重要。随着高占空比固态发射机应用于地基雷达,遮蔽成为普遍问题。\(F_{\text{ecl}}\) 的详细见10.1.2。
  • STC因子:在接收机前或接收机早期级施加可变电压增益 \(F_{\text{stc}}\leq 1\)(见10.1.2)。假设竞争噪声主要出现在STC之后的级,因而输出信噪比与 \(F_{\text{stc}}^{2}\) 成正比。
  • 波束驻留因子:通常为1,但当波束运动导致目标回波接收时天线增益下降时,该因子减小(见10.1.2)。
  • 频率分集因子:只要接收机仍锁定在回波频率上,该因子取1。它确保模糊距离回波在频率分集或捷变下被抑制。
  • 透镜因子:Weil指出该因子会降低低仰角长程目标的可用能量,因为对流层折射导致波束仰角扩展(见7.4节)。它不同于增加系统噪声温度的大气吸收损耗 \(L_{\alpha}\),因此应作为独立因子表达。

由此得到的现代雷达方程[22, Eq. (7.67)]为:

现代雷达: $$ R_{m}^{4}=\frac{P_{av} t_{f} G_{t} G_{r}\lambda^{2}\sigma F_{p}^{2} F_{t}^{2} F_{r}^{2} F_{\text{rdr}}^{2}}{(4\pi)^{3} k T_{s} D_{x}(n^{\prime}) L_{t} L_{\alpha}} \quad(\mathrm{m})\tag{1.26} $$


1.6.3 探测距离的计算方法

在Blake图中,通过简单迭代过程求解存在距离相关大气衰减时的探测距离。引入额外的距离相关响应因子以及随距离可能变化的信号处理模式后,需要更稳健的方法。其基本思路与1.1节相同:分别求解可用与所需能量比作为距离的函数,当两者相等时即得到最大探测距离。该方法既可采用图解方式[21, Section 1.6],也可通过计算机实现[22]。文献[22]中的中间结果图表展示了此过程及距离相关项对可用与所需能量比的影响。


1.6.3.1 探测距离计算示例

表1.1列出了一个非相干二维低空监视雷达的参数。图1.3展示了对于 \(1.0\,\mathrm{m}^{2}\) 目标、恒定仰角 \(1^{\circ}\) 入射、在 \(P_{d}=50\%\) 下的可用与所需能量比 \(E/N_{0}\)\(D_{x}\)。计算使用 Modern Radar System Analysis Software [22],在用户设定的最大距离(此例为150 km)内按100个等距步长向内进行。探测距离 \(R_{m}\) 的确定方法为:检查 \(E/N_{0}\)\(D_{x}\) (以分贝表示)的差值,直到差值为正,再通过插值确定确切的 \(R_{m}\)。此例中 \(R_{m}=132\,\mathrm{km}\)。示例雷达建模为在 \(R_{p}=50\,\mathrm{km}\) 内采用MTI处理,但在该范围内MTI损耗引起的 \(D_{x}\) 增加并未影响 \(R_{m}\),只要 \(R_{p}<R_{m}<R_{u}\)(此例 \(R_{u}=135\,\mathrm{km}\))。

表1.1 示例二维雷达参数

参数 数值 参数 数值
雷达频率 \(f_{0}\) 3.0 GHz 波长 \(\lambda\) 0.10 m
峰值功率 \(P_{t}\) 100 kW 平均功率 \(P_{av}\) 110.8 W
脉冲宽度 \(\tau=\tau_{n}\) \(1.0\,\mu s\) 脉冲重复频率 \(f_{r}\) 1108 Hz
发射线损 \(L_{t}\) 1.0 dB 天线增益 \(G\) 40.0 dB
方位波束宽度 \(\theta_{a}\) \(1.3^{\circ}\) 仰角波束宽度 \(\theta_{e}\) \(2.0^{\circ}\)
方位扫描扇区 \(A_{m}\) \(360^{\circ}\) 帧时间 \(t_{s}\) 6.0 s
每次驻留脉冲数 \(n\) 24 系统温度 \(T_{s}\) 987 K
探测概率 \(P_{d}\) 0.50 虚警概率 \(P_{fa}\) \(10^{-6}\)
基本可探测性因子 \(D\) 2.7 dB 匹配因子 \(M\) 0.8 dB
波束形状损耗 \(L_{p}\) 1.2 dB 其他损耗 \(L_{x}\) 3.3 dB
可探测性因子 \(D_{x}\) 8.0 dB 衰减 \(L_{\alpha}\)(在 \(R_{m}\) 处) 1.8 dB
方向—传播因子 \(F\) 1.0 距离相关因子 \(F_{\text{rdr}}\) 0 dB
目标RCS \(\sigma\) \(1.0\,\mathrm{m}^{2}\) 热噪声下探测距离 \(R_{m}\) 132 km

图1.3 示例雷达能量比与距离关系:实线为信号能量,比虚线为所需能量。MTI处理在 \(R<R_{p}=50\,\mathrm{km}\) 内增加 \(L_{x}\)

1.6.3.2 固态雷达示例

一种可能的固态版本示例雷达被建模,具有相同的平均功率和非模糊距离,但占空比为 \(20\%\),得到的脉宽 \(\tau=180\,\mu s\)。在 \(R<27\,\mathrm{km}\)\(108\,\mathrm{km}<R<162\,\mathrm{km}\) 时出现遮蔽效应,其影响如图1.4所示的可用能量比曲线。此时探测发生在 \(R_{m}=114\,\mathrm{km}\),目标从遮蔽区中出现时几乎完全获得脉冲压缩增益。短程的遮蔽在此情况下对 \(R_{m}\) 无影响。

图1.4 固态雷达(\(D_{u}=20\%\))的能量比与距离关系。


1.6.3.3 含STC的雷达示例

图1.5展示了在示例雷达的非模糊距离前半段内应用STC的结果。上曲线表示对目标散射截面 \(\sigma=1\,\mathrm{m}^{2}\) 的期望效果,此时在 \(R<R_{m}\) 的所有距离上,可用能量都保持在所需值之上。实线曲线表示较小目标(\(\sigma=0.032\,\mathrm{m}^{2}=-15\,\mathrm{dBsm}\))的情况,显示该目标在任何距离上都无法达到 \(P_{d}=50\%\) 的要求。这种性能在空管雷达中是理想的,例如希望探测 \(\sigma \geq 1\,\mathrm{m}^{2}\) 的目标,同时排除几乎所有鸟类。但在现代防空雷达中则不可接受。

图1.5 含STC的示例雷达能量比与距离关系。


1.6.4 垂直覆盖图

Blake图或其他雷达方程解给出的是监视雷达在垂直平面覆盖中的单个点,其方向—传播因子 \(F(\theta_{t})\) 在特定目标仰角 \(\theta_{t}\) 下有效。因为Blake图中的距离 \(R^{\prime}\)\(F\) 成正比,所以一种简单绘制覆盖的方法是将 \(F(\theta)\) 作为仰角 \(\theta\) 的函数计算,并将探测距离包络表达为 \(R^{\prime}(\theta)\)。若大气衰减在波束宽度内变化不大,这种方法即可满足;否则,就需要在扇区内多个仰角上重复迭代衰减计算。[22]的处理方法是:对每个绘制的仰角,利用 \(F, L_{\alpha}, L_{\mathrm{rdr}}\) 随距离变化的函数计算可用能量比,并采用足够小的步长生成平滑的覆盖曲线。所需能量比 \(D_{x}\) 随距离变化(例如由波束形状、MTI损耗和STC引起)的计算,也能得到Blake图难以直接获得的结果。示例雷达的覆盖如图1.6所示。

图1.6的计算条件为:表1.1中的雷达,在地表rms高度偏差 \(\sigma_{h}=3\,\mathrm{m}\) 的环境下运行,仅产生一个低于 \(0.1^{\circ}\) 仰角的显著反射瓣。若无此效应,该反射瓣将延伸到非模糊距离 \(R_{u}=135\,\mathrm{km}\) 以外,但在该区间 \(D_{x}\) 增大超过了可用能量比,因此覆盖被截断在 \(R_{u}\)。同样,由于 \(R<R_{p}=50\,\mathrm{km}\) 时使用了MTI处理,\(D_{x}\) 增大导致波束顶部曲线在穿过该范围时出现轻微凸起。在波束中心 \(\theta=1^{\circ}\) 时,\(P_{d}=50\%\) 的曲线与图1.3中显示的探测距离一致。

图1.6 示例雷达垂直覆盖图:\(P_{d}=90\%\)(内曲线)与 \(P_{d}=50\%\)(外曲线)。

在单一仰角进行的距离计算虽有许多用途,但无法传达垂直覆盖图所展示的信息。垂直覆盖图以直观可解释的形式提供了评估监视雷达性能所需的数据。另一方面,图1.3–1.5所示的曲线则用于分析在特定仰角下控制探测距离的因素。两种表现形式及其他中间结果都是必要的,以避免估算雷达距离时常见的错误。


1.6.5 所需的探测概率

雷达方程中使用的可探测性因子 \(D\) 取决于所需的探测概率。该概率通常任意设为 \(P_{d}=0.5\sim 0.9\)。然而,实际需求取决于接收雷达输出的系统运行模式。单次告警很少足以表明目标存在,因为这种告警可能由噪声或干扰的随机虚警引起。通常需要第二次探测来确认,并建立航迹文件,这意味着单次扫描 \(P_{d}\geq 0.5\)(通常 \(\geq 0.8\))才能可靠地获取并维持航迹。图1.7显示了在不同单次扫描 \(P_{d}\) 下,为获得给定累积航迹获取概率所需的扫描次数。

图1.7 累积航迹获取概率随扫描次数的变化,单次扫描探测概率范围 \(50\%\sim 90\%\)。典型的跟踪扫描系统要求在5次连续扫描中至少有3次探测(引自[21])。

图1.8给出了在不同单次扫描 \(P_{d}\) 下的航迹保持概率。由这两幅图可见,若雷达要支持可靠的跟踪扫描,需满足 \(P_{d}\geq 0.8\)

许多现代雷达(包括多功能阵列雷达和与相控阵跟踪雷达配合的常规搜索雷达)可以调配资源,使得即使搜索模式下的 \(P_{d}\) 远低于0.8,也能建立并维持航迹。在初次探测后安排立即验证驻留,并可能相对于常规搜索功能增加发射能量,可使目标存在得到确认,并高概率建立跟踪,且仅需少量增加时间与能量需求。结果是最大跟踪距离由累积探测概率 \(P_{c}\)(而非单次探测概率 \(P_{d}\))达到可接受值所决定。此时雷达方程中可使用 \(P_{d}\approx 0.3\) 或更低的数值,从而显著扩展雷达作用距离[27]。

图1.8 航迹保持概率随扫描次数的变化,不同单次扫描探测概率下的典型跟踪扫描系统。若连续两次探测失败则丢失航迹(引自[21])。


1.7 雷达方程发展的总结

本章回顾了雷达方程从二战起源到现代雷达应用的发展历程,强调了Lamont Blake的贡献,并建议以他的方程表达方式作为当前与未来雷达性能分析的准确起点。文中介绍了对Blake方程的修正,使其适应现代雷达系统更复杂的特性与工作模式,但其基本方法仍然是可靠估算雷达探测距离的基石。正如Marcum在讨论原始雷达方程时指出的[24, p. 1]:

使用上述方程时可能遇到的陷阱几乎是无穷的,其中许多困难早已被人们认识。

自Marcum发表研究以来,又有新的问题被发现。现代计算技术克服了许多局限,并提供了高精度结果,但同时也可能掩盖其他问题,并鼓励使用一些导致误差远超表面精度的处理方法。避免这些错误正是Blake仔细分析的核心目的。这同样也是后续发展([20-22])以及本章和全书其余部分所追求的目标。

References

[1] Norton, K. A. and A. C. Omberg, "The Maximum Range of a Radar Set," Proc. IRE, Vol. 35, No. 1, January 1947, pp. 4-24. Reprinted in [2].

[2] Barton, D. K., (ed.), The Radar Equation, Dedham, MA: Artech House, 1974.

[3] North, D. O., "An Analysis of the Factors Which Determine Signal/Noise Discrimination in Pulsed Carrier Systems," RCA Laboratories Technical Report PTR-6C, June 25, 1943. Reprinted: Proc. IEEE, Vol. 51, No. 67, July 1963, pp. 1015-1027.

[4] IEEE Standard 100, The Authoritative Dictionary of IEEE Standards Terms, 7th ed., New York: IEEE Press, 2000.

[5] Haeff, A. W., "Minimum Detectable Radar Signal and Its Dependence Upon Parameters of Radar Systems," Proc. IRE, Vol. 34, No. 11, November 1946, pp. 857-861.

[6] Blake, L. V., "Interim Report on Basic Pulse-Radar Maximum Range Calculations," Naval Research Laboratory Report 1106, November 1960.

[7] Blake, L. V., "Recent Advancements in Basic Radar Range Calculation," IRE Trans. on Military Electronics, Vol. MIL-5, No. 2, April 1961, pp. 154-164. Reprinted in [2].

[8] Blake, L. V., "A Guide to Basic Pulse-Radar Maximum-Range Calculation, Part 1-Equations, Definitions, and Aids to Calculation," Naval Research Laboratory Report 5868, December 28, 1962.

[9] Blake, L. V., "A Guide to Basic Pulse-Radar Maximum-Range Calculation, Part 1-Equations, Definitions, and Aids to Calculation," Naval Research Laboratory Report 6930, December 23, 1969. Summarized in [2].

[10] Blake, L. V., "A Guide to Basic Pulse-Radar Maximum-Range Calculation, Part 2-Derivation of Equations, Bases of Graphs, and Additional Explanations," Naval Research Laboratory Report 7010, December 31, 1969. [11] Blake, L. V., "Prediction of Radar Range," Chapter 2 in [15, 16].

[12] Blake, L. V., Radar Range-Performance Analysis, Lexington, MA: D. C. Heath, 1980; Dedham, MA: Artech House, 1986.

[13] Blake, L. V., "The Effective Number of Pulses per Beamwidth for a Scanning Radar," Proc IRE, Vol. 41, No. 6, June 1953, pp. 770-774, with addendum: December 1953, p. 1785. Reprinted in [2].

[14] Kerr, D. E., (ed.), Propagation of Short Radio Waves, Vol. 13 in MIT Radiation Laboratory Series, New York: McGraw-Hill, 1951. Reprinted (CD ROM edition): Norwood, MA: Artech House, 1999.

[15] Skolnik, M. I., (ed.), Radar Handbook, New York: McGraw-Hill, 1970.

[16] Skolnik, M. I., (ed.), Radar Handbook, 2nd ed., New York: McGraw-Hill, 1990.

[17] Hall, W. M., "Prediction of Pulse Radar Performance," Proc. IRE, Vol. 44, No. 2, February 1956, pp. 224-231. Reprinted in [2].

[18] Hall, W. M., "General Radar Equation," Space/Aeronautics R and D Handbook, 1962-63. Reprinted in [2].

[19] Barton, D. K., Radar System Analysis, Englewood Cliffs, NJ: Prentice-Hall, 1964; Dedham, MA: Artech House 1976.

[20] Barton, D.. K., Modern Radar System Analysis, Norwood, MA: Artech House, 1988.

[21] Barton, D. K., Radar System Analysis and Modeling, Norwood, MA: Artech House, 2005.

[22] Barton, D. K., Modern Radar System Analysis Software and User's Manual, Version 3.0, Norwood, MA: Artech House, 2007.

[23] Rice, S. O., "Mathematical Analysis of Random Noise," Bell Sys. Tech. J., Vol. 23, No. 3, July 1944, pp. 282-332 and Vol. 24, No. 1, January 1945, pp. 461-556. Reprinted: Selected Papers on Noise and Stochastic Processes, (N. Wax, ed.), New York: Dover Publ., 1954.

[24] Marcum, J. I., "A Statistical Theory of Target Detection by Pulsed Radar, RAND Corp. Research Memo RM-754, December 1, 1947. Reprinted: IRE Trans. on Information Theory, Vol.IT-6, No. 2, April 1960, pp. 59-144; excerpted: Detection and Estimation, (S. S. Haykin, ed.), Stroudsburg, PA: Halstad Press, 1976, pp. 57-121.

[25] Swerling, P., "Probability of Detection for Fluctuating Targets," RAND Corp. Research Memo RM-1217, March 17, 1954. Reprinted: IRE Trans. on Information Theory, Vol. IT-6, No. 2, April 1960, pp. 269-308; Detection and Estimation, (S. S. Haykin, ed.), Stroudsburg, PA: Halstad Press, 1976, 122-158.

[26] Weil, T. A., "Atmospheric Lens Effect, Another Loss for the Radar Range Equation," IEEE Trans. on Aerospace and Electronic Systems, Vol. AES-9, No. 1, January 1973, pp. 51-54. Reprinted in [2].

[27] Barton, D. K., "Maximizing Firm-Track Range on Low-Observable Targets," IEEE International Conf. Radar-2000, Washington, DC, May 8-12, 2000, pp. 24-29.