基于压缩感知的机载相控阵雷达低空风切变风速估计方法技术

一种基于压缩感知的机载相控阵雷达低空风切变风速估计方法。其包括：利用空时插值法构造参考距离单元与待检测距离单元的变换矩阵，从而获得构造待检测单元的杂波协方差矩阵的独立同分布样本，获得杂波协方差矩阵的估计值，进而实现杂波抑制；利用雷达主瓣宽度作为先验信息，构造风切变场的广义空间导向矢量；利用信号谱宽作为先验信息，构造风切变场的广义时间导向矢量；利用广义空间导向矢量和广义时间导向矢量，构造基于风速的风切变场空时基字典，构造稀疏基矩阵；对步骤1)进行杂波抑制后的回波信号进行观测，利用稀疏基矩阵对回波信号进行恢复，实现风速估计。本方法在脉冲数较少且信噪比较低时仍然能够获得较为精确的风场速度估计结果。

【技术实现步骤摘要】
基于压缩感知的机载相控阵雷达低空风切变风速估计方法
本专利技术属于雷达信号处理
，特别设计一种基于压缩感知的机载相控阵雷达低空风切变风速估计方法。技术背景低空风切变通常是指高度在600m以下风向、风速突然变化的气象现象，是对航空运输安全威胁最大的天气现象之一。由于风切变现象具有时间短、尺度小、强度大等特点，从而带来了探测难、预报难等一系列问题。当飞机在起降阶段进入强风切变区域时，由于缺乏足够的调节空间，若操作不当，极易造成飞行事故。因此低空风切变探测技术的研究已成为现代航空运输领域的一个重要课题。机载气象雷达可以对雷雨、风切变、湍流等气象现象进行探测与预警，是飞机实时检测航路气象信息的重要设备。但当该雷达处于下视工作状态时，风切变回波常常处于强杂波背景中。杂波抑制的好坏程度直接影响了对风切变场风速估计结果的准确性。相较于传统单天线体制气象雷达，相控阵雷达具有灵活度高、扫描速度快、易于波束形成等特点。由于在接收回波中增加了目标信号的空域信息，使得相控阵雷达较传统单天线体制雷达在强杂波背景下杂波抑制效果更佳，能够更好地实现对目标的检测,相控阵雷达因此受到了广泛重视。目前世界上一些先进的气象雷达研究机构已经开始了新一代相控阵气象雷达的研究。空时自适应处理(Space-timeAdaptiveProcessing，STAP)是应用在相控阵雷达杂波抑制上的关键技术，STAP技术通过在空时域自适应地形成凹口，能够在有效抑制地杂波的同时，提高雷达的目标检测能力。目前，STAP技术已相对成熟地应用在点目标的检测上，但是风切变属于分布式目标，传统的STAP技术不能直接应用。目前，针对机载雷达的风切变检测的方法主要是针对单天线体制雷达的，但应用在相控阵体制雷达中的风速估计方法鲜有提及。传统机载单天线体制气象雷达风速估计方法主要包括基于时域分析的脉冲对法(PulsePairProcessing，PPP)和基于频域分析的快速傅立叶变换法(FastFourierTransform，FFT)等。虽然这些方法在脉冲数较多、信噪比较高时风速估计性能较好，但在脉冲数较少、信噪比较低的时候，其风速估计性能下降。
技术实现思路
为了解决上述问题，本专利技术的目的在于提供一种能够提高参数估计精度的基于压缩感知的机载相控阵雷达低空风切变风速估计方法。为了达到上述目的，本专利技术提供的基于压缩感知的机载相控阵雷达低空风切变风速估计方法包括按顺序进行的下列步骤：1)利用空时插值法构造参考距离单元与待检测距离单元的变换矩阵，从而获得构造待检测单元的杂波协方差矩阵的独立同分布样本，获得杂波协方差矩阵的估计值，进而实现杂波抑制；2)利用雷达主瓣宽度作为先验信息，构造风切变场的广义空间导向矢量；3)利用信号谱宽作为先验信息，构造风切变场的广义时间导向矢量；4)利用上述广义空间导向矢量和广义时间导向矢量，构造基于风速的风切变场空时基字典，构造稀疏基矩阵；5)对上述步骤1)进行杂波抑制后的回波信号进行观测，利用上述稀疏基矩阵对回波信号进行恢复，实现风速估计。在步骤2)中，所述的利用雷达主瓣宽度作为先验信息，构造风切变场的广义空间导向矢量的方法是：利用雷达主瓣宽度，构造基于方位和俯仰的角信号密度函数，以此构造区别于传统点目标的分布式目标空间导向矢量。在步骤3)中，所述的利用信号谱宽作为先验信息，构造风切变场的广义时间导向矢量的方法是：利用信号谱宽，构造气象回波的频谱扩展函数，并以此构造区别于传统点目标的气象目标时间导向矢量。在步骤4)中，所述的利用上述广义空间导向矢量和广义时间导向矢量，构造基于风速的风切变场空时基字典，构造稀疏基矩阵的方法是：将待估的风速区间离散化，建立离散的待估风速空间，利用离散的风速空间建立时域基字典，将广义空间导向矢量和时域基字典做Kronecker积，得到基于风速的空切变场的空时基字典，再将空时基字典与随机观测矩阵做乘积，得到稀疏基矩阵。在步骤5)中，所述的对上述步骤1)进行杂波抑制后的回波信号进行观测，利用上述稀疏基矩阵对回波信号进行恢复，实现风速估计的方法是：将观测信号向稀疏基矩阵中各原子对应的一维子空间进行正交投影，计算投影误差，与观测信号投影误差最小，即该原子与观测信号相关性最大的原子即为回波信号的恢复原子；构成该恢复原子的中心风速，即为待估距离单元的中心风速。本专利技术提供的基于压缩感知的机载相控阵雷达低空风切变风速估计方法是利用先验信息建立分布式气象目标的精确空时导向矢量，然后利用其回波的稀疏特性恢复信号，进而实现参数估计。本方法在脉冲数较少且信噪比较低时仍然能够获得较为精确的风场速度估计结果。附图说明图1为本专利技术基于压缩感知的机载相控阵雷达低空风切变风速估计方法流程图。图2为机载前视阵雷达阵模型图。图3为雷达接收信号的空时二维谱。图4为35号距离单元风切变回波的频谱图。图5为采样脉冲64，信噪比5dB，杂噪比40dB时的风速估计结果。图6为采样脉冲64，信噪比0db，杂噪比40dB时的风速估计结果。图7为采样脉冲32，信噪比5db，杂噪比40dB时的风速估计结果。图8为采样脉冲32，信噪比0db，杂噪比40dB时的风速估计结果。具体实施方法下面结合附图和具体实例对本专利技术提供的基于压缩感知的机载相控阵雷达低空风切变风速估计方法进行详细说明。如图1所示，本专利技术提供的基于压缩感知的机载相控阵雷达低空风切变风速估计方法包括按顺序进行的下列步骤：1)利用空时插值法构造参考距离单元与待检测距离单元的变换矩阵，从而获得构造待检测单元的杂波协方差矩阵的独立同分布数据样本，获得杂波协方差矩阵的估计值，进而实现杂波抑制；机载前视阵雷达阵模型如图1所示，设载机速度为VR，载机速度方向与天线的阵面轴夹角为90°，天线阵元数为N，脉冲重复频率为fr，相干处理脉冲数为K，阵元间距d＝0.5λ，其中λ为雷达发射脉冲的波长。在本专利技术中，xl表示第l(l＝1,2,…,L)个待检测单元的NK×1维空时快拍数据，其表达式如下：xl＝sl+cl+nl(1)其中，sl、cl、nl分别表示第l个待检测单元的风切变场回波、杂波及噪声，其中假设杂波无起伏无模糊，噪声为加性高斯白噪声。对于第l个距离单元内的风切变场，雷达对其的采样数据可以写成一个N×K的矩阵Sl。其中，Sl的第n行、第k列元素表示雷达第n(n＝1,2,…N)个阵元、第k(k＝1,2,…K)个脉冲对风场回波的采样，当该距离单元内雷达波束照射范围内共有Q个气象散射点时，其具体表达式如下：其中和分别表示第q(q＝1,2,…,Q)个气象散射点的空间角频率和时间角频率，θq、分别表示该气象散射点相对于雷达的方位角和俯仰角，Rq为第q个散射点与载机的斜距，为天线接收方向图。将上面的Sl展开成为NK×1维列向量，即为风切变场回波快拍sl。则雷达全距离单元内的回波信号可以表示为：X＝[x1x2…xL]T(3)图3表示在杂噪比40dB、信噪比5dB条件下接收信号的空时二维谱。可以看出，前视阵地杂波呈非均匀的椭圆分布，并且地杂波的强度要远大于风切变信号的强度，风切变信号几乎全部被杂波淹没。对于第l号待检测单元，将其水平方位角离散化，得到{θm}|m＝1,2,…M，其中M表示离散化的方位角个数。设该检测单元对应的杂波空时本文档来自技高网...

【技术保护点】
一种基于压缩感知的机载相控阵雷达低空风切变风速估计方法，其特征在于，其包括按顺序进行的下列步骤：1)利用空时插值法构造参考距离单元与待检测距离单元的变换矩阵，从而获得构造待检测单元的杂波协方差矩阵的独立同分布数据样本，获得杂波协方差矩阵的估计值，进而实现杂波抑制；2)利用雷达主瓣宽度作为先验信息，构造风切变场的广义空间导向矢量；3)利用信号谱宽作为先验信息，构造风切变场的广义时间导向矢量；4)利用上述广义空间导向矢量和广义时间导向矢量，构造基于风速的风切变场空时基字典，构造稀疏基矩阵；5)对上述步骤1)进行杂波抑制后的回波信号进行观测，利用上述稀疏基矩阵对回波信号进行恢复，实现风速估计。

【技术特征摘要】
1.一种基于压缩感知的机载相控阵雷达低空风切变风速估计方法，其特征在于，其包括按顺序进行的下列步骤：1)利用空时插值法构造参考距离单元与待检测距离单元的变换矩阵，从而获得构造待检测距离单元的杂波协方差矩阵的独立同分布数据样本，获得杂波协方差矩阵的估计值，进而实现杂波抑制；2)利用雷达主瓣宽度作为先验信息，构造风切变场的广义空间导向矢量；广义空间导向矢量的具体形式为：其中，为雷达主瓣方向俯仰角，θi为雷达主瓣方向方位角，表示方位角为θi，俯仰角为的方位上的点目标的空间导向矢量；为确定性角信号密度函数，表示风切变场在雷达主瓣方向方位角θi和俯仰角上的扩展；是以分布源中心DOA为对称中心的单峰对称函数，其形式为：其中，σθ表示θi方向上的角度扩展，表示方向上的角度扩展；σω为σθ和的集合，表示目标在空间上的扩展；N表示天线阵元数；d表示阵元间距，d＝0.5λ，λ表示雷达发射脉冲的波长；3)利用信号谱宽作为先验信息，构造风切变场的广义时间导向矢量；风切变场的广义时间导向矢量的具体形式为：其中，fd表示风切变回波信号的多普勒频率，表示雷达主瓣方向方位角为θi和俯仰角为时，速度为v的点目标的时间导向矢量，g(σf)表示频率扩展函数，其中σf表示风切变信号回波的谱宽；K表示相干处理脉冲数，fr表示脉冲重复频率，λ表示雷达发射脉冲的波长；4)利用上述广义空间导向矢量和广义时间导向矢量，构造基于风速的风切变场空时基字典，构造稀疏基矩阵；5)对上述步骤1)进行杂波抑制后的回波信号进行观测，利用上述稀疏基矩阵对回波信号进行恢复，实现风速估计。2.根据权利要求1所述的...

【专利技术属性】
技术研发人员：李海，周盟，吴仁彪，
申请(专利权)人：中国民航大学，
类型：发明
国别省市：天津;12