雷达回波模拟仿真（三）：杂波（以LFM为例）matlab

  就杂波来源而言，它包括来自地面、大海、云雨等天气、以及一些飞行的动物的回波。特殊的例如电离的流星尾迹回波或是被称为箔条所引起的回波，通常也被认为是杂波。从类型上来讲，杂波可以大体分为面杂波和体杂波两部分。其中，面杂波主要指地杂波和海杂波，通常用单位面积的杂波横截面积（ σ 0 \sigma^0 σ0）,即为单位面积的杂波横截面积，计算式如下： σ 0 = σ c / A c \sigma^0=\sigma_c/A_c σ0=σc​/Ac​ 其中， σ c \sigma_c σc​是面积 A c A_c Ac​上的杂波截面积。 σ 0 \sigma^0 σ0是一个无量纲的参数。 地杂波，除个别特殊的孤立杂波外，是一种分布散射的现象。由于地杂波是散射回波，它与实际地面场景有密切联系，不同的地形地貌杂波散射强度不同。   同样的，海杂波的特性主要决定于海面形状。因为海的形状粗糙度受风的影响较大，所以风力也是决定海杂波大小的一个因素。此外，水表面的张力也会影响粗糙度，从而影响到海面的雷达回波。研究海杂波有两个难点，一是建立海回波测量值与决定海况的环境因素之间的定量关系非常困难，二是海面随时间连续地发生变化，动态性的模拟也是困难之处。

  机载PD雷达，由于其向下俯视性和运动性的特点，其回波的频谱要比地基等固定或低速雷达的频谱宽得多，再加上雷达天线方向图和波束照射区域的运动变化等原因，其杂波表现出再时域上随机起伏。因此，研究机载雷达杂波，必须要从该类杂波的频域入手，分析其频谱特效。   根据杂波散射单元所处的天线波束主旁瓣的位置不同，将机载雷达杂波频谱分为主瓣、旁瓣和高度线三个杂波部分。其杂波频谱较为详细的分析图如图所示   图中， V R V_R VR​为载机速度， V T V_T VT​为目标速度， ϕ 0 \phi_0 ϕ0​为主瓣波束方向与载机速度矢量间的夹角， ϕ T \phi_T ϕT​为目标速度矢量与雷达和目标连线的夹角， ϕ \phi ϕ为地面杂波散射单元A与载机速度矢量间夹角， θ B \theta_B θB​为雷达天线主瓣波束宽度。从图中可以看出，就杂波幅度而言，机载PD雷达主瓣杂波幅度通常最高，高度线杂波次之，旁瓣杂波一般较低；而就频谱宽度而言，主瓣杂波和高度线杂波的宽度较窄，而旁瓣杂波会很宽。

  雷达通过电磁散射回波来发现目标，这种电磁散射的物理特性是雷达工作的基本原理之一。其中，除双基地雷达外，大部分雷达关注的物理量是电磁波的后向散射。对于研究的地、海杂波而言，后向散射系数 σ 0 \sigma^0 σ0是定量研究雷达杂波的重要概念。后向散射系数的定义是地海表面单位面积所对应的雷达横截面积，其值表示在散射波方向上的等效散射总功率与入射方向上入射功率密度和照射面积乘积的比值。   雷达的后向散射系数与场景地形、工作频率、掠射角、海况、极化的方式、风速等因素有关。滨海地区的复杂环境，既有地杂波，又有海杂波，而这两种杂波的形成机理和影响因素有很大的不同，所以必须充分考虑每种情况的后向散射特性才能准确把握和研究。

  美国MIT林肯实验室收集了分布于北美多个地区的总共42个观测站的地杂波测量数据，并建立了地杂波数据库，得出了准确度较高的雷达杂波强度数据。这些数据针对各种不同的地形，具有较好的“地面真实性”。对这些地杂波数据进行分析的结果表明： （1）在低掠射角(0°-10°)， σ 0 \sigma^0 σ0均值几乎线性增加；在某一小的地形区域内， σ 0 \sigma^0 σ0总是在某一均值附近变化；而在其他条件相同的前提下，不同频段的测量值也不同。 （2）在中掠射角(10°-70°)，Nathanson提出了一个近似的计算方法： σ 0 = γ s i n ψ \sigma^0=\gamma sin\psi σ0=γsinψ，其中， ψ \psi ψ为雷达波束主瓣掠射角(擦地角)。 （3）大掠射角(70°-90°)， σ 0 \sigma^0 σ0的值可以很大，这个范围内的散射是由那些小平面将入射能量反射回雷达造成的。另外，接近垂直入射的 σ 0 \sigma^0 σ0的值也会受到天线方向图和天线增益的影响。当散射平面为光滑平面时， σ 0 = G \sigma^0=G σ0=G   由以上的分析，我们可以推导求解一些地杂波后向散射系数的模型。这些模型是后向散射系数与雷达参数和环境参数在内的相关影响因素之间的关系模型。与地形相关的散射系数模型主要有：修正的常量 γ \gamma γ模型、修正的Morchin模型、Barton模型等。通过比较，并结合机载雷达工作掠射角范围的实际情况，考虑采用适用掠射角范围较广，角度值接近0°-90°的修正Morchin模型。其数学表达式如下： σ 0 = A σ c 0 s i n ψ λ + u c o t 2 β 0 e x p [ − t a n 2 ( B − ψ ) t a n 2 β 0 ] \sigma^0= \frac {A\sigma^0_c sin\psi}{\lambda}+ucot^2\beta_0exp[-\frac{tan^2(B-\psi)}{tan^2\beta_0}] σ0=λAσc0​sinψ​+ucot2β0​exp[−tan2β0​tan2(B−ψ)​] 其中， ψ \psi ψ为地形单元擦地角，地貌类型为沙漠，且当 ψ < ψ c \psi