本发明专利技术公开了一种固态双偏振天气雷达距离副瓣识别与抑制方法，包括分别逐库计算水平通道和垂直通道数据当前距离库和相邻距离库的相关系数；分别逐库计算水平通道和垂直通道相关系数的标准差；如果当前库的相关系数均大于指定门限，并且当前库的相关系数的标准差均小于指定门限，则判定当前距离库回波受到距离副瓣影响。否则，判定当前距离库不受距离副瓣影响；如果判定当前库受距离副瓣影响，则根据当前库的相关系数对水平通道和垂直通道的回波功率分别进行修正。否则，直接输出原始的I/Q数据。该方法降低了距离副瓣对固态双偏振天气雷达的不利影响，有助于提升固态双偏振雷达对灾害天气的探测效果。灾害天气的探测效果。灾害天气的探测效果。
【技术实现步骤摘要】
一种固态双偏振天气雷达距离副瓣识别与抑制方法
[0001]本专利技术涉及雷达信号处理研究领域，特别是一种固态双偏振天气雷达距离副瓣识别与抑制方法。
技术介绍
[0002]脉冲压缩(脉压)是固态天气雷达广泛采用的技术之一，通常采用线性调频(LFM)或者非线性调频(NLFM)脉压信号形式。LFM信号有其固有的缺点，那就是匹配滤波器输出波形的距离副瓣高。为了降低信号的距离副瓣，需要在压缩滤波器中引进某种加权，这将在一定程度上引入脉压输出主瓣的损耗，降低系统的灵敏度。而采用NLFM则可以避免因为加权所引起的失配损失，又可以得到较低的旁瓣电平。脉压信号波形设计的重点是保证其具有良好的自相关性能，由于信号的自相关函数和功率谱之间是一一对应关系，因此LFM/NLFM信号的设计就归结为按要求的信号功率谱来设计信号，一般通过窗函数来实现。
[0003]在处理距离副瓣方面，鉴福升等(2001)提出了一种对窗函数进行组合设计与优化的思路，能够快速灵活地设计出所需信号，但是最优波形的设计是其中的难点。陶广源等(2003)研究了相位编码的脉冲压缩方法，采用迭代加权最小二乘算法来设计优化副瓣抑制滤波器，从而使脉压输出的主瓣副瓣比(主副比)达到了较高的实用水平，该方法同样关注的是波形设计。Kurdzo等(2014)基于遗传算法设计了一个优化的非线性调频波形，相比于窗函数滤波器，该技术可以在保证距离副瓣抑制效果的同时减少信噪比的损失，缺点是需要更大的带宽。
[0004]在当前的固态天气雷达中，LFM脉冲压缩方法还有相当大范围的应用，再加上天气观测的特殊需求，这也会导致雷达回波中可能出现相当数量的距离副瓣回波。这种虚假回波引起数据质量下降，影响对天气现象的认识。因此，如何在波形设计的基础上，从雷达回波角度去识别和抑制距离副瓣，就显得尤为重要。
技术实现思路
[0005]本专利技术要解决的技术问题是针对上述现有技术的不足，而提供一种固态双偏振天气雷达距离副瓣识别与抑制方法，该固态双偏振天气雷达距离副瓣识别与抑制方法能从雷达回波角度识别和抑制距离副瓣，能解决当前固态天气雷达受脉冲压缩导致的距离副瓣影响较大，造成观测天气时出现虚假回波从而引起数据质量下降的技术问题。
[0006]为解决上述技术问题，本专利技术采用的技术方案是：
[0007]一种固态双偏振天气雷达距离副瓣识别与抑制方法，包括如下步骤。
[0008]步骤1、计算相关系数CC_Gate_H和CC_Gate_V：固态双偏振天气雷达具有水平通道回波和垂直通道回波，水平通道回波和垂直通道回波均具有沿射线方向的N个距离库；其中，N≥500；对每个距离库均分别计算相关系数CC_Gate_H和CC_Gate_V；其中，CC_Gate_H为当前第i个距离库中水平通道回波和相邻第j个距离库中水平通道回波的相关系数；CC_Gate_V为当前第i个距离库中垂直通道回波和相邻第j个距离库中垂直通道回波的相关系
数；其中，i≠j。
[0009]步骤2、计算标准差CC_Gate_Std_H和CC_Gate_Std_V：对每个距离库均分别计算标准差CC_Gate_Std_H和CC_Gate_Std_V；其中，CC_Gate_Std_H为当前距离库和左右相邻a个距离库水平通道回波相关系数的标准差；其中，2≤a≤8；CC_Gate_Std_V为当前距离库和相邻a个距离库垂直通道回波相关系数的标准差。
[0010]步骤3、距离副瓣影响判断，具体包括如下步骤：
[0011]步骤3
‑
1、设定相关系数门限CC_Thr和标准差门限CC_Std_Thr。
[0012]步骤3
‑
2、判断：当前库的CC_Gate_H和CC_Gate_V均大于相关系数门限CC_Thr，并且当前库的CC_Gate_Std_H和CC_Gate_Std_V均小于标准差门限CC_Std_Thr，则判定当前距离库回波受到距离副瓣影响；否则，判定当前距离库不受距离副瓣影响。
[0013]步骤4、距离副瓣抑制：当步骤3判断当前距离库回波受到距离副瓣影响时，则需根据当前库的相关系数CC_Gate_H和CC_Gate_V，对水平通道和垂直通道的回波功率或信噪比分别进行修正，以降低或消除距离副瓣对当前距离库回波的影响。
[0014]步骤1中，相关系数CC_Gate_H和CC_Gate_V计算时，根据当前第i个距离库和相邻第j个距离库的脉冲电压以及信号功率值的比值计算得到。
[0015]相关系数CC_Gate_H和CC_Gate_V的计算公式为：
[0016][0017][0018]式中，M表示每个距离库中水平通道回波或垂直通道回波脉冲电压的采集总次数。
[0019]m表示每个距离库中水平通道回波或垂直通道回波脉冲电压的采集顺序号，1≤m≤M。
[0020]表示当前第i个距离库中水平通道第m次采集的回波脉冲电压的共轭。
[0021]V
j_H
(m)表示相邻第j个距离库中水平通道第m次采集的回波脉冲电压。
[0022]表示当前第i个距离库中水平通道回波信号功率的估算值。
[0023]表示相邻第j个距离库中水平通道回波信号功率的估算值。
[0024]表示当前第i个距离库中垂直通道第m次采集的回波脉冲电压的共轭。
[0025]V
j_V
(m)表示相邻第j个距离库中垂直通道第m次采集的回波脉冲电压。
[0026]表示当前第i个距离库中垂直通道回波信号功率的估算值。
[0027]表示相邻第j个距离库中垂直通道回波信号功率的估算值。
[0028]j的取值，需根据雷达实际使用的信号波形脉压输出后的主瓣与副瓣的相对距离
位置及距离库采样率共同确定。
[0029]j＝i+1或j＝i+2或j＝i+3；第j个距离库与第i个距离库之间的间隔要大于等于雷达距离分辨力。
[0030]步骤2中，a＝6，每组标准差CC_Gate_Std_H和CC_Gate_Std_V计算时，均使用当前距离库、当前距离库相邻的前三个距离库以及当前距离库相邻的后三个距离库，共计7个距离库的相关系数进行标准差计算。
[0031]步骤4中，距离副瓣抑制时，水平通道和垂直通道使用相同的修正相关系数，修正相关系数为当前距离库计算得到的相关系数CC_Gate_H和CC_Gate_V的均值。
[0032]步骤4中，距离副瓣抑制时，根据当前库的相关系数CC_Gate_H和CC_Gate_V，对水平通道和垂直通道的信噪比分别进行修正，以降低或消除距离副瓣对当前距离库回波的影响；其中，信噪比修正方法为：当前库信噪比减去当前库修正相关系数的平方与信噪比的乘积。
[0033]步骤3中，相关系数门限CC_Thr为0.5，标准差门限CC_Std_Thr为0.1。
[0034]本专利技术具有如下有益效果：
[0035]1、能从雷达回波角度识别和抑制距离副瓣，能解决当前本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种固态双偏振天气雷达距离副瓣识别与抑制方法，其特征在于：包括如下步骤：步骤1、计算相关系数CC_Gate_H和CC_Gate_V：固态双偏振天气雷达具有水平通道回波和垂直通道回波，水平通道回波和垂直通道回波均具有沿射线方向的N个距离库；其中，N≥500；对每个距离库均分别计算相关系数CC_Gate_H和CC_Gate_V；其中，CC_Gate_H为当前第i个距离库中水平通道回波和相邻第j个距离库中水平通道回波的相关系数；CC_Gate_V为当前第i个距离库中垂直通道回波和相邻第j个距离库中垂直通道回波的相关系数；其中，i≠j；步骤2、计算标准差CC_Gate_Std_H和CC_Gate_Std_V：对每个距离库均分别计算标准差CC_Gate_Std_H和CC_Gate_Std_V；其中，CC_Gate_Std_H为当前距离库和左右相邻a个距离库水平通道回波相关系数的标准差；其中，2≤a≤8；CC_Gate_Std_V为当前距离库和相邻a个距离库垂直通道回波相关系数的标准差；步骤3、距离副瓣影响判断，具体包括如下步骤：步骤3
‑
1、设定相关系数门限CC_Thr和标准差门限CC_Std_Thr；步骤3
‑
2、判断：当前库的CC_Gate_H和CC_Gate_V均大于相关系数门限CC_Thr，并且当前库的CC_Gate_Std_H和CC_Gate_Std_V均小于标准差门限CC_Std_Thr，则判定当前距离库回波受到距离副瓣影响；否则，判定当前距离库不受距离副瓣影响；步骤4、距离副瓣抑制：当步骤3判断当前距离库回波受到距离副瓣影响时，则需根据当前库的相关系数CC_Gate_H和CC_Gate_V，对水平通道和垂直通道的回波功率或信噪比分别进行修正，以降低或消除距离副瓣对当前距离库回波的影响。2.根据权利要求1所述的固态双偏振天气雷达距离副瓣识别与抑制方法，其特征在于：步骤1中，相关系数CC_Gate_H和CC_Gate_V计算时，根据当前第i个距离库和相邻第j个距离库的脉冲电压以及信号功率值的比值计算得到。3.根据权利要求2所述的固态双偏振天气雷达距离副瓣识别与抑制方法，其特征在于：相关系数CC_Gate_H和CC_Gate_V的计算公式为：相关系数CC_Gate_H和CC_Gate_V的计算公式为：式中，M表示每个距离库中水平通道回波或垂直通...
【专利技术属性】
技术研发人员：邵世卿，施春荣，周静艳，陈超，陆小虎，
申请(专利权)人：中船重工鹏力南京大气海洋信息系统有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。