技术领域
[0001] 本
发明属于
机载气象雷达低空风切变检测技术领域，特别是涉及一种基于组合空时主通道自适应处理的低空风切变风速估计方法。
背景技术
[0002] 低空风切变属于民航运输领域的灾难性天气之一，它具有持续时间短、作用区域小、瞬间强度大、危害性高等特点。当飞机在
起飞或者进近阶段遭遇风切变时，由于飞行高度较低，飞行员无法及时调整飞机
姿态，将有可能导致空难的发生，因此针对低空风切变的检测和预警就成为当前民航领域的一项重要课题，而风切变风速估计作为整个风切变检测流程的
基础，直接影响了风切变检测的准确程度。
[0003] 机载气象雷达是民航飞机必不可少的重要
电子设备，它是民航飞机的“双眼”。当机载气象雷达工作在下视模式时，风切变的回波
信号会被淹没在地杂波背景中，因此地杂波的抑制效果将直接影响低空风切变的风速估计
精度。传统的地杂波抑制方法有杂波图法、多扫描方法以及零陷
滤波器法等，其本质均在于寻找合适的凹口，在抑制地杂波的同时保留风切变信号。但当阵列存在幅相误差时，杂波谱将在
空域展宽，降低了处理器的
自由度，严重影响杂波抑制的性能，进而影响风速估计结果的准确性。
[0004]
相控阵雷达采用多个天线接收回波信号，相对于传统的单天线雷达，其接收的回波信号中包含目标空域信息，使得空时联合处理成为可能，因此相控阵雷达在强杂波背景下的目标检测和参数估计性能更加优越。空时自适应处理(Space-time Adaptive Processing, STAP)是机载相控阵雷达杂波抑制与目标检测的关键技术，它能够利用雷达回波的空时耦合特性实现杂波抑制和信号匹配。迄今为止，传统的STAP技术多针对点目标进行检测和参数估计，而针对低空风切变等分布式目标检测与参数估计的参考文献相对较少，尤其是幅相误差情况下的机载气象雷达低空风切变检测研究还未有文献报道。
发明内容
[0005] 为了解决上述问题，本发明的目的在于提供一种存在幅相误差时，仍然可以精确估计风速的基于组合空时主通道自适应处理的低空风切变风速估计方法。
[0006] 为了达到上述目的，本发明提供的基于组合空时主通道自适应处理的低空风切变风速估计方法包括按顺序进行的下列步骤：
[0007] 1)从机载设备中读取飞机状态参数、机载气象雷达参数作为先验信息；
[0008] 2)利用步骤1)中读取的先验信息对机载气象雷达回波信号进行
相位补偿，剔除飞机运动带来的速度偏移；
[0009] 3)对机载气象雷达各距离单元相位补偿后的回波信号进行杂波距离依赖性矫正，并估计出待测距离单元的杂波协方差矩阵；
[0010] 4)构造
降维联合空时变换矩阵，对矫正后的机载气象雷达回波信号进行降维自适应处理；
[0011] 5)利用步骤4)获得的降维联合空时变换矩阵和步骤3)估计得到的杂波协方差矩阵求解组合空时主通道自适应降维处理器的最优权矢量，实现杂波抑制和风切变信号匹配；
[0012] 6)利用上述最优权矢量构造代价函数并通过搜索代价函数的最大值来估计出待测距离单元的风速；
[0013] 7)重复步骤3)-6)，依次处理机载气象雷达工作范围内所有距离单元的机载气象雷达回波信号，得到风速随距离的变化曲线。
[0014] 在步骤1)中，所述的飞机状态参数、机载气象雷达参数包括飞机速度V、飞行高度 H、机载气象雷达
波长λ、天线扫描的主瓣方位
角θ0、
俯仰角 和脉冲重复
频率fr。
[0015] 在步骤2)中，所述的利用步骤1)中读取的先验信息对机载气象雷达回波信号进行相位补偿是：利用上述先验信息对机载气象雷达回波信号进行相位补偿，剔除飞机运动带来的速度偏移，以消除飞机的运动对参数估计的影响。
[0016] 在步骤3)中，所述的对机载气象雷达各距离单元相位补偿后的回波信号进行杂波距离依赖性矫正，并估计出待测距离单元的杂波协方差矩阵的方法是：利用杂波配准法矫正杂波距离依赖性，并利用参考距离单元的机载气象雷达回波信号作为样本估计得到待测距离单元的杂波协方差矩阵。
[0017] 在步骤4)中，所述的构造降维联合空时变换矩阵，对矫正后的机载气象雷达回波信号进行降维自适应处理的方法是：利用切比
雪夫权对待测距离单元内的空域波束和时域通道分别进行加权，并通过定义两个
正交的方向向量族构造降维联合空时变换矩阵，对待测距离单元的机载气象雷达回波信号进行自适应降维处理。
[0018] 在步骤5)中，所述的利用步骤4)获得的降维联合空时变换矩阵和步骤3)估计得到的杂波协方差矩阵求解组合空时主通道自适应降维处理器的最优权矢量，实现杂波抑制和风切变信号匹配的方法是：根据LCMV准则，利用步骤4)获得的降维联合空时变换矩阵和步骤3)估计得到的杂波协方差矩阵计算得到组合空时主通道自适应降维处理器的最优权矢量，并利用求解得到的最优权矢量对机载气象雷达回波信号进行滤波。
[0019] 在步骤6)中，所述的利用上述最优权矢量构造代价函数并通过搜索代价函数的最大值来估计出待测距离单元的风速的方法是：循环更新用于重构风场目标空时导向矢量的归一化多普勒中心频率，得到不同的风速匹配结果，最大
输出信号功率对应的速度值即为待测距离单元内的低空风切变信号风速的估计值。
[0020] 本发明提供的基于组合空时主通道自适应处理的低空风切变风速估计方法是立足于新一代相控阵体制的机载气象雷达，将组合空时主通道自适应处理器应用到低空风切变的风速估计中。在存在幅相误差影响的情况下，通过计算降维联合空时变换矩阵构造降维处理器，并利用降维处理器的最优权矢量实现杂波抑制和风场目标信号匹配，进而获得风速的有效估计。
附图说明
[0021] 图1为本发明提供的基于组合空时主通道自适应处理的低空风切变风速估计方法的实现
框图；
[0022] 图2为机载前视阵回波仿真示意图；
[0023] 图3为组合空时主通道自适应降维处理器原理图；
[0024] 图4为机载气象雷达回波信号空时二维谱；
[0025] 图5为第74～78号距离单元滤波器频响图；
[0026] 图6为风速估计结果对比图；
[0027] 图7为幅相误差情况下风速估计结果对比图；
具体实施方式
[0028] 下面结合附图和具体
实施例对本发明提供的基于组合空时主通道自适应处理的低空风切变风速估计方法进行详细说明。
[0029] 如图1所示，本发明提供的基于组合空时主通道自适应处理的低空风切变风速估计方法包括按顺序进行的下列步骤：
[0030] 1)从机载设备中读取飞机状态参数、机载气象雷达参数作为先验信息；
[0031] 从机载设备中读取飞机速度V、飞行高度H、机载气象雷达波长λ、天线扫描的主瓣方位角θ0、俯仰角 和脉冲重复频率fr在内的参数作为
信号处理的先验信息。
[0032] 2)利用步骤1)中读取的先验信息对机载气象雷达回波信号进行相位补偿，剔除飞机运动带来的速度偏移；
[0033] 利用上述先验信息对机载气象雷达回波信号进行相位补偿，剔除飞机运动带来的速度偏移，以消除飞机的运动对参数估计的影响，公式如下：x(l)＝X(l)e-jΩ，其中X(l)表示第l个距离单元未进行相位补偿的机载气象雷达回波信号，
[0034] 3)对机载气象雷达各距离单元相位补偿后的回波信号进行杂波距离依赖性矫正，并估计出待测距离单元的杂波协方差矩阵；
[0035] 利用杂波配准法矫正杂波距离依赖性，并利用参考距离单元的机载气象雷达回波信号估计得到待测距离单元的杂波协方差矩阵 其中，Υ(j)＝Pj,lx(j),(1≤j≤L,j≠l)为第j号参考距离单元矫正后的回波信号，x(j)为该参考距离单元机载气象雷达回波信号，Pj,l为杂波配准补偿矩阵。
[0036] 4)构造降维联合空时变换矩阵，对矫正后的机载气象雷达回波信号进行降维自适应处理；
[0037] 如图2所示，设飞机速度为V，飞行高度为H，在机载平台上沿航向的垂直方向均匀放置N元线阵，阵元间隔d＝0.5λ，其中λ为雷达波长。图中，θ和 为地面散射单元
水平方位角和俯仰角，θ0和 为天线扫描的主瓣方位角和俯仰角，假设杂波无起伏无模糊，噪声为加性高斯白噪声。
[0038] 不失一般性，假设在相干处理时间(Coherent Processing Interval，CPI)内，相干处理脉冲数为K，脉冲重复频率为fr，机载气象雷达在其工作范围内有L个距离单元，第l个距离单元的空时二维快拍数据x(l),(l＝1,2,…,L)可表示为：
[0039] x(l)＝s(l)+c(l)+n (1)
[0040] 其中，c(l)为第l个距离单元的地杂波，n为高斯白噪声，s(l)为该距离单元的低空风切变信号，可表示为：
[0041]
[0042] 其中，Γ为该距离单元的风场回波复幅度，A(ψ0,fl)为该距离单元的风场空时导向矢量，ψ0为风场信号空间锥角， fl为该距离单元内风场目标回波的多普勒中心频率，At(fl)为该距离单元的风场信号的时域导向矢量， 表示Kronecker积，As(ψ0) 为该距离单元的风场信号的空域导向矢量。时域导向矢量At(fl)和空域导向矢量As(ψ0)可以分别表示为：
[0043] At(fl)＝at(fl)⊙g(σf) (3)
[0044] As(ψ0)＝as(ψ0)⊙g(ψ0) (4)
[0045] 其中， 表示该距离单元中多普勒中心频率为fl的点目标的时域导向矢量，⊙为Hadamard积， 表示频率扩展
函数。 表示该距离单元风场中心上点目标的空域导
向矢量， 为风场信号的角度扩展函数，其中，
σθ表示天线扫描的主瓣方位角θ0方向上的角度扩
展， 表示天线扫描的主瓣俯仰角 方向上的角度扩展。
[0046] 如图3所示，假设fll表示待测距离单元待搜索的风场信号归一化多普勒频率(fll取值范围为[-1,1])。利用切比雪夫权对空域波束和时域通道分别进行加权，可以有效削弱空时二维旁瓣杂波，使杂波仅局限在正交的空域主波束和时域主通道内。在二维加权空时响应平面(即ψ-f响应平面)内，分别以ψ＝ψ0和f＝fll为中线的空域主波束和时域主通道可以定义为如下两个正交的方向向量族：
[0047]
[0048] 其中， 和 分别表示所有沿ψ＝ψ0和f＝fll加权导向矢量的集合。f为风场信号占据的归一化多普勒
频率范围，ψ∈[0,π]为风场信号所在的空间锥角范围，B(ψ0,f)和 B(ψ,fll)为不同的加权空时导向矢量，分别可以进一步表示为：
[0049]
[0050]
[0051] 其中，Bt(f)＝{At(f),f∈[-1 1]}表示归一化多普勒频域内所有时间导向矢量的集合，Bs(ψ)＝{As(ψ),ψ∈[0π]}表示空域范围内所有空间导向矢量的集合。Bs(ψ0)和Bt(fll)分别表示为：
[0052]
[0053]
[0054] 式中，Qt＝[qt1 qt1…qtK]T和Qs＝[qs1 qs1…qsN]T分别表示时域和空域的切比雪夫加权矢量，IN和IK分别为N×N和K×K的单位阵， 为NK×N的矩阵，为NK×K的矩阵。 和 组成的矩阵 的秩为N+K-1 (即矩阵
中任意(N+K-1)个列向量相互独立)，由此可以构造降维联合空时变换矩阵
为：
[0055]
[0056] 其中，T'fll由Tfll的前N-1列组成。
[0057] 5)利用步骤4)获得的降维联合空时变换矩阵和步骤3)估计得到的杂波协方差矩阵求解组合空时主通道自适应降维处理器的最优权矢量，实现杂波抑制和风切变信号匹配；
[0058] 根据LCMV准则可以将第l个距离单元内的组合空时主通道自适应降维处理器描述为如下数学优化问题：
[0059]
[0060] 其中，RT(l)＝TlHR(l)Tl为降维杂波协方差矩阵，R(l)为估计得到的该距离单元的杂波协方差矩阵。BT(l)为该距离单元内降维变换后低空风切变信号的空时导向矢量，即：
[0061]
[0062] 其中，B(ψ0,fll)为构造的该距离单元内低空风切变信号的空时导向矢量，fll为上文中假定的归一化多普勒频率。此时，组合空时主通道自适应降维处理器的最优权矢量为：
[0063]
[0064] 6)利用上述最优权矢量构造代价函数并通过搜索代价函数的最大值来估计出待测距离单元的风速；
[0065] 当机载气象雷达回波信号通过组合空时主通道自适应降维处理器后，可以在降维的基础上实现地杂波抑制和风切变信号匹配，此时，可用下式作为代价函数估计待测距离单元风场回波信号的多普勒频率：
[0066]
[0067] 其中，xT(l)＝TlHx(l)，进而得到第l个距离单元的风速估计值为：
[0068]
[0069] 7)重复步骤3)-6)，依次处理机载气象雷达工作范围内所有距离单元的机载气象雷达回波信号，得到风速随距离的变化曲线；
[0070] 利用步骤3)-6)所述方法估计出每个距离单元的风速，并绘出风速随距离的变化曲线，即完成低空风切变场的风速估计过程。
[0071] 本发明提供的基于组合空时主通道自适应处理的低空风切变风速估计方法的效果可以通过以下仿真结果进一步说明。
[0072] 仿真条件描述：低空风切变场分布于飞机前方约8.5-16.5km处，机载气象雷达的天线阵为阵元数N＝8的均匀线阵，阵元间距d＝λ/2，主瓣波束水平方位角为90°，俯仰角为0°，机载气象雷达工作波长为0.05m，脉冲重复频率为7000Hz，最小可分辨距离150m，相干处理脉冲数K＝32，
信噪比5dB，杂噪比40dB，飞机速度为75m/s，飞行高度为600m。
[0073] 图4为仿真的机载气象雷达回波信号空时二维谱，从图4中可见，当机载气象雷达工作在前视模式时，地杂波的空时二维谱为椭圆形，低空风切变的
功率谱总体为一条窄带，并且地杂波的回波功率远大于风场的回波功率，导致风场回波的多普勒信息完全淹没在地杂波的多普勒信息中，这为低空风切变的检测和参数估计带来了困难。
[0074] 图5为利用本发明方法对第74～78号待测距离单元进行空时自适应处理时的滤波器频响图，从图中可以看出，基于组合空时主通道的自适应处理器在零频杂波带上形成了明显的凹口，可以有效地滤除地杂波，并且在归一化多普勒频率分别为0.15、0.22、0.25、 0.3、0.35处形成了有效的信号增益以匹配风场信号。
[0075] 图6为无误差时本发明方法与其他方法的风速估计结果对比图，从图可以看出，利用本发明方法在强杂波下可以获得比较准确的风速估计结果，并且在8.5-16.5km范围内，低空风切变信号的风速随距离变化而表现正负变化的反“S”形分布，这是低空风切变信号最主要的特征。
[0076] 图7为本发明方法在不同幅相误差下的风速估计结果。其中图7(a)为仅存在幅度误差，且幅度误差分别为1％，2％，3％，4％，5％时的风速估计结果图；图7(b)为仅存在相位误差，且相位误差分别为1°，2°，3°，4°，5°时的风速估计结果图；图7(c) 为同时存在不同幅度和相位误差时，本发明方法的风速估计结果对比图。从图中可以看出，随着幅度和相位误差的增大，杂波谱展宽，但本发明方法依然可以较准确地得到整体的风速估计结果。
