基于空时最优处理器的低空风切变风速估计方法专利检索-机载气象雷达机载设施专利检索查询-专利查询网

积极推动地理标志专门立法

2022-03-10 地理标志，立法，知识产权
保护知识产权是对创新最大的激励

2022-03-10 保护知识产权，创新，激励
谢商华：加快制定知识产权基本法

2022-03-10 知识产权基本法
擦亮“双奥之城”品牌

2022-03-10 双奥，知识产权
让冰雪运动“热”力全开

2022-03-10 冰雪运动，知识产权
携手共奋进　走好强国路

2022-03-10 强国，知识产权
坚持创新引领　方能稳中求进

2022-03-10 创新，稳中求进，知识产权
答好“两张卷” 奋进新征程

2022-03-10 知识产权
专家解读政府工作报告中的创新和知识产权相关部署

2022-03-10 政府工作报告，创新，知识产权
今年政府工作报告指出：加强知识产权保护和运用

2022-03-10 政府工作报告，知识产权保护

基于空时最优处理器的低空风切变风速估计方法

阅读：1136发布：2020-06-27

IPRDB可以提供基于空时最优处理器的低空风切变风速估计方法专利检索，专利查询，专利分析的服务。并且本发明涉及一种基于空时最优处理器的低空风切变风速估计方法。本发明属于机载气象雷达信号处理领域。一种基于空时最优处理器的低空风切变风速估计方法，包括以下步骤：1)从关联的机载设备中读取飞机状态参数、雷达天线扫描参数等作为先验信息；2)选取待检测距离单元，利用空时插值法校正地杂波的距离依赖性；3)利用先验信息补偿飞机速度；4)估计待检测距离单元的地杂波协方差矩阵；5)利用空时最优处理器抑制地杂波并进行风速估计；6)重复步骤2)‑5)，依次对范围内所有距离单元的回波数据进行风速估计。本发明针对强杂波背景下风速估计精度变差的问题，在机载气象雷达中引入空时信息，提出了一种基于空时最优处理器的低空风切变风速估计方法。仿真结果证明了该方法的有效性。，下面是基于空时最优处理器的低空风切变风速估计方法专利的具体信息内容。

权利要求

1. 一种基于空时最优处理器的低空风切变风速估计方法，其特征在于：针对待检测距离单元的回波数据，所述的风速估计方法包括按顺序进行的下列步骤： 1) 从关联的机载设备中读取飞机状态参数、雷达天线扫描参数等作为先验信息； 2) 选取待检测距离单元，利用空时插值法校正地杂波的距离依赖性； 3) 利用先验信息补偿飞机速度； 4) 估计待检测距离单元的地杂波协方差矩阵； 5) 利用空时最优处理器抑制地杂波并进行风速估计，该步骤还包括:根据低空风切变的分布特性，设定多普勒谱宽与中屯、频率，构造风切变信号的协方差矩阵，包括 5a)、风切变目标Rs的参数化模型为其中，P为回波信号功率，fs为信号空间频率，fo为信号多普勒中屯、频率，Of为信号多普勒谱宽，且有上式中瓜为KXK维矩阵，且

，Bs为NXN维全1矩阵；加）、风切变目标Rs的参数化模型可由其回波模型推导得到：某一距离单元内低空风切变回波模型可写为其中Zn,k为接收信号幅度;f为信号归一化多普勒频率，且通常认为/~zV(乂回波信号相关函数为由信号相关函数即可推得风切变目标的参数化模型;其中N为阵元数，K为脉冲数； 6) 重复步骤2)-5)，依次对范围内所有距离单元的回波数据进行风速估计。

说明书全文

基于空时最优处理器的低空风切变风速估计方法

技术领域

[0001] 本发明属于机载气象雷达低空风切变检测技术领域，特别是涉及一种基于空时最优处理器的低空风切变风速估计方法。

背景技术

[0002] 机载气象雷达可对雷雨、风切变、湍流等灾害天气进行探测与预警，是飞机实时感知航路气象的重要设备，其对危险天气的探测能力对保证飞机行驶的安全性至关重要。低空风切变常指高度600米以下风向、风速突然变化的气象现象，通常是强气流由空中冲击到地面后向四周扩散的过程。低空风切变具有发生突然、时间短、尺度小、强度大和不易检测等特点，当飞机起飞和着陆进入强低空风切变区域时，会面临由风向、风速变化带来的升力损失，若处理不得当或缺乏足够的调节空间就可能遭遇飞行事故。因此，低空风切变检测技术的研究已是航空运输领域的一项重要课题。

[0003] 机载气象雷达下视检测低空风切变时，有用信号会被强地杂波所覆盖。在低空风切变风速估计的工作流程中，地杂波抑制一直是最大的难点。传统的在地杂波背景下估计风场风速的方法包括杂波图法、基于参数化模型的谱估计方法、利用模式分析的扩展Prony 方法、零陷滤波器法等，韩雁飞等提出了基于压缩感知的低空风切变风速估计方法，能够在脉冲数较少时实现风速的精确估计，但该方法在应用时需要确定杂波谱宽，且并未涉及风速估计过程中的强杂波抑制问题。这些方法的本质在于寻找合适的凹口在抑制地杂波的同时保留风切变信号，但在强杂波背景下，地杂波难以完全消除，残余杂波仍会覆盖掉风切变信号的多普勒信息，影响风速估计结果的准确性。

[0004] 相对于传统机载气象雷达的单天线体制，相控阵天线体制具有灵活度高、扫描速度快、易于波束赋形的优点，其回波信号中包含目标空间采样信息。目前国外已经对机载相控阵气象雷达开展了相关的研究工作。2013年9月，美国国家大气研究中心（National Center for Atmospheric Research,NCAR)表示，美国正在研制的下一代机载气象雷达采用双极化相控阵体制。空时自适应处理（Space Time Adaptive Processing,STAP)是机载相控阵雷达杂波抑制与目标检测的关键技术，空时最优处理器可在对有用信号输出功率不变的条件下，尽可能地抑制杂波、干扰及噪声。在我国，STAP技术在机载预警雷达上的应用已经颇为成熟，但此种应用多针对点目标，不能直接用于低空风切变等分布式目标的检测与估计，且机载气象雷达工作在前视状态，前视阵的杂波谱在距离上不平稳，因此传统的空时最优处理器不能直接应用于机载气象雷达中。

发明内容

[0005] 为了解决上述问题，本发明的目的在于提供一种在高杂噪比、低信噪比时，仍然可以精确估计风场风速的基于空时最优处理器的低空风切变风速估计方法。

[0006] 为了达到上述目的，本发明提供的基于空时最优处理器的低空风切变风速估计方法，包括按顺序进行的下列步骤：

[0007] 1)从关联的机载设备中读取飞机状态参数、雷达天线扫描参数等作为先验信息；

[0008] 2)选取待检测距离单元，利用空时插值法校正地杂波的距离依赖性；

[0009] 3)利用先验信息补偿飞机速度；

[0010] 4)估计待检测距离单元的地杂波协方差矩阵；

[0011] 5)利用空时最优处理器抑制地杂波并进行风速估计；

[0012] 6)重复步骤2)-5)，依次对范围内所有距离单元的回波数据进行风速估计。

[0013] 在步骤1)中，所述的从关联的机载设备中读取飞机状态参数、雷达天线扫描参数等作为先验信息是读取飞机速度、飞行高度、天线扫描的主瓣方位角和俯仰角等作为先验信息。

[0014] 在步骤2)中，所述的选取待检测距离单元，利用空时插值法校正地杂波的距离依赖性是利用空时插值法，构造某一距离单元与待检测距离单元的插值变换矩阵T，使用该距离单元数据与插值变换矩阵T相乘，使该距离单元内地杂波的空时二维分布与待检测距离单元趋于一致。依次处理所有距离单元数据，可得到与待检测距离单元地杂波分布特性相同的多个独立同分布（Independent and Identically Distributed, IID)样本。

[0015] 在步骤3)中，所述的利用先验彳目息补偿飞机速度是利用上述先验彳目息对接收彳目号进行相位补偿，剔除机载平台运动带来的多普勒偏移，以消除机载平台的运动性对参数估计的影响。

[0016] 在步骤4)中，所述的估计待检测距离单元的地杂波协方差矩阵是利用步骤2)中得到的多个独立同分布样本求取其协方差矩阵的统计平均值。

[0017] 在步骤5)中，所述的利用空时最优处理器抑制地杂波并进行风速估计是指，构造可自适应抑制机载前视阵地杂波并积累待检测单元内风切变信号的最优处理器;根据低空风切变的分布特性，设定多普勒谱宽与中心频率，构造风切变信号的协方差矩阵，并联合步骤4)中计算得到的待检测距离单元内地杂波协方差矩阵求解最优处理器的权矢量;按多普勒中心频率做循环，更新最优处理器权矢量，当输出信号功率最大时，此时设定的多普勒中心频率即为待检测距离单元内风切变信号的多普勒频率估计值，进而能够计算待检测距离单元内的风场速度。

[0018] 在步骤6)中，所述的重复步骤2)_5)，依次对范围内所有距离单元的回波数据进行风速估计是指更新待检测距离单元，再次利用步骤2)_5)估计其风切变信号多普勒中心频率。将范围内所有距离单元处理完毕后，得到低空风切变场风速随距离的变化曲线。

[0019] 本发明提供的基于空时最优处理器的低空风切变风速估计方法是将相控阵体系引入到机载气象雷达中，利用空时插值法校正前视阵地杂波的距离依赖性，利用最优处理器自适应地抑制地杂波并匹配低空风切变信号，估计风切变信号的多普勒频率后得到风场速度估计。在强杂波条件下，本发明方法仍能有效地抑制地杂波并精确地估计风场风速。

附图说明

[0020] 图1为本发明提供的基于空时最优处理器的低空风切变风速估计方法的工作原理图；[0021 ]图2为机载前视阵回波仿真示意图；

[0022]图3为本发明提供的基于空时最优处理器的低空风切变风速估计方法的基本流程图

[0023]图4为补偿飞机速度后的仿真信号空时二维谱；

[0024]图5为空时插值法补偿地杂波距离依赖性效果图；

[0025]图6为空时处理器频响特性；

[0026] 图7为本发明方法多普勒估计结果；

[0027] 图8为本发明方法风场速度估计结果。

具体实施方式

[0028] 下面结合附图和具体实施例对本发明提供的基于空时最优处理器的低空风切变风速估计方法进行详细说明。

[0029] 如图1所示，本发明提供的基于空时最优处理器的低空风切变风速估计方法包括按顺序进行的下列步骤：

[0030] 1)从关联的机载设备中读取飞机状态参数、雷达天线扫描参数等作为先验信息：

[0031] 从机载导航设备中读取飞机速度、飞行高度、天线扫描的主瓣方位角和俯仰角等作为信号处理的先验信息。

[0032] 2)选取待检测距离单元，利用空时插值法校正地杂波的距离依赖性：

[0033] 本发明假设天线体制为均匀线阵。如图2所示为机载前视阵回波仿真示意图，设定载机速度为V;脉冲重复频率为fr;阵元间隔为d;天线主瓣波束水平方位角为θ〇,俯仰角为外;杂波单元水平方位角和俯仰角表示为Θ和Ρ ;杂波单元空间角频率和时间角频率分别表示为巧(0，f)和崎俄的，且有

[0035]雷达回波信号x = s+c+n，其中X为低空风切变信号，c为地杂波，η为高斯白噪声，矩阵形式可写为

[0036] χ=[χι Χ2 ··· xl]T (2)

[0037] 其中 Xfh? # 禮…4,( 将每个距离单元的接收数据排成NK X 1的列向量N为阵元数，K为脉冲数，L为最大距离单元数。

[0038] 空时插值法通过一个插值变换矩阵使得训练距离单元的地杂波空时分布与待检测距离单元一致，从而达到校正地杂波距离相关性的目的，具体实现步骤如下：

[0039] a)选取第i号距离单元为待检测距离单元，构造待检测距离单元的空时导向矢量矩阵Vl。首先将水平方位角均匀分成Μ份，得到对应化丨'^^^然后构造每个水平方位角 0 m和距离单元俯仰角只对应的空时导向矢量ν(6〇 = ν,(1)®ν、(ίϋ4*

[0041] 式中j表示虚数单位，可以得到待检测距离单元的空时导向矢量矩阵为

[0042] Vi=[v(9i) ν(θ2) ··· ν(θΜ)] (4)

[0043] b)选取第j号距离单元为参考距离单元，构造该参考距离单元的空时导向矢量矩阵％，可得到第j号与第i号距离单元的插值变换矩阵为其中（·）+表示求伪逆运算。处理后的参考距离单元数据为n=L,lxj，其与待检测距离单元内地杂波的空时二维分布一致。

[0044] c)更新参考距离单元，同样利用空时插值法进行处理，可以得到地杂波分布特性与待检测距离单元相同的多个独立同分布样本。

[0045] 3)利用先验信息补偿飞机速度：

[0046] 利用上述先验信息对接收信号进行相位补偿，剔除机载平台运动带来的多普勒频移，以消除机载平台的运动性对参数估计的影响。

[0047] 4)估计待检测距离单元的地杂波协方差矩阵：

[0048]插值变换处理后，估计待检测距离单元地杂波数据的统计协方差矩阵为

[0050] 5)利用空时最优处理器抑制地杂波并进行风速估计：[0051 ] 可在抑制地杂波的同时匹配低空风切变信号的空时最优处理器可以描述为如下的数学优化问题：

[0053] 其中ω为最优处理器的权矢量，Rc+n为杂波噪声协方差矩阵，Rs为低空风切变信号协方差矩阵。

[0054] 构造代价函数

[0055] Η( ω ,λ)= oHRc+no+A(l-oHRso ) (7)

[0056] 其中，λ为拉格朗日算子，，式(6)可等效为求解下式广义特征值问题：

[0057] Rc+no =XRS ω (8)

[0058] 此时，拉格朗日算子λ可以看做是相应的广义特征值。由于矩阵Rc+n和1是半正定的，所以上式得到的所有广义特征值都是非负实数。当式(8)成立时，代价函数取得最小值 λ。将R;1左乘到式⑶两边，得到

[0060] 则有= 其中ρ{ · }表示求解矩阵最大特征值对应的特征向量。由于地杂波功率远高于噪声功率，可近似认为= έ ,得到

[0061] (10)

[0062]由上式可知，求解最优处理器权矢量ω需要杂波协方差矩阵食^和信号协方差矩阵Rs两个关键信息，其中杂波协方差矩阵灸卩1可由式(5)得到，而在实际检测过程中，风切变信号协方差矩阵Rs是未知的，上式仅代表对目标Rs进行积累。

[0063]风切变目标Rs的参数化模型为

[0065]其中，P为回波信号功率，fs为信号空间频率，fo为信号多普勒中心频率，〇f为信号多普勒谱宽，且有

[0069] 上式中，Bt为KXK维矩阵，且，BS为NXN维全1矩阵。

[0070] 风切变目标Rs的参数化模型可由其回波模型推导得到。低空风切变属于分布式目标，其多普勒频率随距离的分布呈反"S"形，且单距离单元内的风切变信号在多普勒域内存在连续性展宽。基于此种特性，设定待检测距离单元内多普勒中心频率和多普勒谱宽，可以构造信号协方差矩阵。某一距离单元内低空风切变回波模型可写为

[0072] 其中zn,k为接收信号幅度;f为信号归一化多普勒频率，且通常认为/、。此时，回波信号相关函数为

[0073]

[0074]由信号相关函数即可推得风切变目标：fes的参数化模型如式(11)所示。

[0075] 将信号协方差矩阵良带入式（10)中，可以求得最优处理器权矢量ω，此处理器可以自适应地抑制地杂波，并积累多普勒频率处于：内的目标信号。地杂波抑制与信号匹配的结果为

[0076] y= ω Hxi (15)

[0077] 更新求解权矢量ω时设定的多普勒中心频率fQ，得到对不同多普勒频段目标信号进行匹配滤波的最优处理器权矢量，并分别求解处理器输出信号功率。当输出功率最大时，最优处理器可以有效地抑制地杂波并积累低空风切变信号，此时设定的多普勒中心频率即为待检测距离单元内风切变信号的多普勒频率估计值,进而得到该距离单元的风速估计值为

[0079] 谱宽〇f表示了低空风切变场的风速变化率，可取为实际观测中风速谱宽的统计平均值，本步骤本质是固定空间锥角和谱宽〇f，利用空时最优处理器对待检测距离单元的多普勒中心频率做一维搜索。

[0080] 6)重复步骤2)_5)，依次对范围内所有距离单元的回波数据进行风速估计：

[0081] 利用步骤2)_5)所述算法估计出每个距离单元的风场速度，即完成低空风切变场的检测过程。

[0082]本发明提供的基于空时最优处理器的低空风切变风速估计方法的效果可以通过以下仿真结果进一步说明。

[0083]仿真条件描述:低空风切变场分布于飞机前方约8.5-16.5km处，天线阵为阵元数为N=8的均匀线阵，阵元间距(1 = λ/2,主瓣波束水平方位角为60°，俯仰角为0°，机载气象雷达工作波长为〇.〇5m，脉冲重复频率为7000Hz，最小可分辨距离150m，相干处理脉冲数Κ = 64,信噪比5dB，杂噪比40dB，飞机速度为75m/s，飞行高度为600m，归一化多普勒谱宽〇f = 0.05〇

[0084] 图3为本发明提供的基于空时最优处理器的低空风切变风速估计方法的基本流程图。其关键步骤在于:利用空时插值法校正地杂波的距离依赖性，构造符合待检测距离单元地杂波二维分布特征的多个IID样本;设计自适应空时处理器，在抑制地杂波的同时对风切变信号进行积累，运用目标检测的思想，搜索风切变信号的多普勒频率。

[0085] 图4为补偿飞机速度后的仿真信号空时二维谱示意图，为便于观察，此时设置风切变信号与地杂波功率相同。如图所示，地杂波的空时二维谱呈现椭圆形分布，低空风切变信号的空时二维谱表现为主瓣波束空间锥角方向上的一条"窄带"。本发明提供方法利用空时最优处理器自适应地在地杂波方向形成凹口，在风切变信号方向形成增益，进而得到风场多普勒频率的有效估计结果。

[0086]图5(a)给出了利用空时插值法校正地杂波距离依赖性前的杂波功率谱，由于机载前视阵雷达的地杂波分布具有距离依赖性，其空时二维谱有明显发散现象，直接使用回波信号估计地杂波协方差矩阵，会致使最优处理器的地杂波抑制性能严重下降。图5(b)给出了利用空时插值法向第20号距离单元做杂波谱补偿的结果，可以看出，补偿后的地杂波空时二维谱明显变窄。

[0087] 图6以20号距离单元为例，显示了 f 0 = 0.5, 〇f = 0.05时的最优处理器频响特性。可以看出，最优处理器沿地杂波分布方向形成自适应凹口，同时对中心频率为fo的宽带信号进行功率积累。按fQ循环，更新最优处理器的权矢量，得到最大输出功率对应的又，即可计算此距离单元的风场风速。

[0088] 图7为本发明方法多普勒估计结果。如图所示，在8.5-16.5km范围内，信号的多普勒频率随距离呈现反"S"形分布，进一步计算风速梯度和告警因子后可完成低空风切变的检测过程。

[0089] 图8为本发明方法与传统方法的风速估计结果对比。在同等杂噪比与信噪比条件下，本发明方法处理效果明显优于自适应频域置零法与固定凹口法等传统方法。结果证明，本发明提供的基于空时最优处理器的低空风切变风速估计方法可以在强杂波条件下，有效地抑制地杂波，并精确估计风场速度。

标题	发布/更新时间	阅读量
一种机载气象雷达视频处理方法-专利编号CN106918803A	2020-05-14	457
基于FPGA的机载气象雷达-专利编号CN106597449A	2020-05-12	593
一种机载气象雷达系统教学试验台-专利编号CN110085090A	2020-05-19	1120
基于FPGA的机载气象雷达-专利编号CN106597449B	2020-05-11	322
一种机载气象雷达教学模拟系统-专利编号CN102129795B	2020-05-16	726
一种飞机机载气象雷达测试系统-专利编号CN104345302A	2020-05-17	576
一种机载气象雷达教学模拟系统-专利编号CN102129795A	2020-05-18	972
一种机载气象雷达雷暴识别方法-专利编号CN108896995A	2020-05-13	1149
具有旋转天线的机载气象雷达-专利编号CN104111459B	2020-05-20	402
具有旋转天线的机载气象雷达-专利编号CN104111459A	2020-05-15	769

基于空时最优处理器的低空风切变风速估计方法

基于空时最优处理器的低空风切变风速估计方法

技术领域

背景技术

发明内容

附图说明

具体实施方式

IPRDB

热门服务

关于我们

友情链接

联系方式