一种高频地波雷达电离层杂波抑制方法技术

一种高频地波雷达电离层杂波抑制方法，本发明专利技术涉及高频地波雷达电离层杂波抑制方法。本发明专利技术为了解决现有方法较多地用来抑制高频雷达中的窄带射频干扰以及瞬态冲击干扰，不能很好的抑制电离层的干扰，且不宜应用于宽波束雷达，以及会造成一阶Bragg峰信息的缺失以及海浪回波能量的减弱的问题。本发明专利技术相对传统杂波抑制算法，利用电离层多聚集区杂波之间的互相关性来完成对消处理，在抑制电离层杂波的同时，保留海浪回波及一阶Bragg峰的信息，减少正常海洋回波能量的下降。经过本发明专利技术方法处理后，雷达回波谱中的强电离层杂波被明显抑制，杂波抑制幅度提高，回波信杂噪比获得显著改善。本发明专利技术用于高频地波雷达探测领域。

【技术实现步骤摘要】
一种高频地波雷达电离层杂波抑制方法
本专利技术涉及高频地波雷达探测领域，具体涉及高频地波雷达电离层杂波抑制方法。
技术介绍
高频地波雷达具有独特的超视距和全天候探测能力，广泛应用于舰船和飞机等运动目标的探测和海洋观测。在实际应用中，由于地波雷达发射天线做不到完全理想地使电磁波沿海面传播，因此虽然大部分能量能够沿水平面发射，但仍不可避免地存在小部分能量向上发射，并经由电离层传播反射回来，被雷达所接收形成了电离层杂波。电离层杂波能量很强，由于电离层电子浓度和高度的时变特性，杂波通常在距离维和多普勒维均有扩展，表现出明显的非平稳信号特征，这些杂波经常使雷达在相应距离上性能急剧下降甚至失效，严重影响了雷达的探测性能。在海态反演应用中，电离层杂波能量有时甚至超过海浪回波的能量，对于海洋信息的提取十分不利，其不确定性也为杂波消除带来了极大的困难。因此，对高频地波雷达电离层杂波抑制一直是研究的难点和亟待解决的问题。目前现有的电离层杂波抑制方法包括自适应波束形成技术、极化滤波技术等。自适应波束形成技术是通过让阵列天线方向图的主瓣指向所需方向，使波束的副瓣与外部干扰的分布相适应来抑制干扰，但它对天线主瓣波束宽度有一定的限制,因而不宜应用于宽波束雷达。极化滤波技术是利用目标和干扰在极化域的特征来进行滤波处理，由于该技术需要较多的水平天线使系统变得复杂，故在高频地波雷达的实际工程应用中很少采用。常用于高频地波雷达杂波抑制的特征值分解法以及基于时频分析、基于(多普勒)频域滤波的干扰检测和抑制方法，则较多地用来抑制高频雷达中的窄带射频干扰以及瞬态冲击干扰(主要是流星余迹干扰和雷电干扰)。采用特征值分解法抑制高频地波雷达电离层杂波的同时，也会造成一阶Bragg峰信息的缺失以及海浪回波能量的减弱。电离层一般包括D层、E层和F层。高频地波雷达所接收的电离层杂波可能存在的传播模式主要包括Es(或E)层直接回波、F层直接回波以及Es(或E)层海面模式回波等，根据其杂波分布情况可分为聚集区和扩散区。经由天波、海洋调制的混合路径产生的电离层杂波，存在于电离层高度之上且存在较大的距离扩展，由于闪烁等原因，其时频分布显得模糊散乱，有时强于海洋回波，这样的杂波区域称为扩展区。对应于天顶方向的Es层直接回波距离分布较集中且能量通常很强，有着聚集的时频分布，形成的集中分布条带状区域称为聚集区，通常出现在雷达距离-多普勒谱(R-D谱)杂波区的最底部。因Es层的强反射作用，在R-D谱的上方存在Es层杂波的一跳和二跳回波等，这些杂波区域能量一般较强且很集中，也属于聚集区。由于电离层Es层的厚度非常薄，通常情况下为百米量级，而高频地波雷达的距离分辨率为数公里，因此Es层杂波及其多跳回波在回波谱图上出现的距离相对集中，能量也主要集中在邻近的几个距离单元上(多为3至5个距离元)。聚集区杂波对同距离的有用信号形成压制性影响，必须对其进行抑制。
技术实现思路
本专利技术的目的是为了解决现有方法较多地用来抑制高频雷达中的窄带射频干扰以及瞬态冲击干扰，不能很好的抑制电离层的干扰，且不宜应用于宽波束雷达，以及会造成一阶Bragg峰信息的缺失以及海浪回波能量的减弱的问题，而提出一种高频地波雷达电离层杂波抑制方法。一种高频地波雷达电离层杂波抑制方法包括以下步骤：步骤一：由雷达距离-多普勒谱(R-D谱)，根据电离层杂波分布特性，按距离由近至远确定对应每个电离层杂波聚集区的杂波中心距离元ki,i＝0,1,…,p-1，对应Es层杂波及其一跳回波、二跳回波等，p为电离层杂波聚集区的个数；步骤二：根据步骤一确定的杂波中心距离元ki，计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关系数步骤三：设置电离层杂波聚集区的相关系数参考阈值Th1，比较由步骤二得到的距离相关系数的值与参考阈值Th1的大小关系，如果则判定该杂波中心距离元对应杂波区域为聚集区，数据有效；否则，判定该杂波中心距离元对应杂波区域不是聚集区，数据无效，舍弃；筛选出所有满足条件的杂波中心距离元ki，按从小到大的顺序依次排列并赋值给dc,c＝0,1,…,m-1，表示聚集区的杂波中心距离元，m为有效数据个数，执行步骤四；步骤四：利用电离层Es层杂波与其多跳回波之间的互相关性，计算电离层Es层聚集区的杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数j＝1,2,…,m-1；步骤五：设置电离层杂波聚集区的互相关系数参考阈值Th2，将步骤四求得的每个互相关系数与参考阈值Th2作比较，如果则认为第j跳回波聚集区数据有效；否则，舍弃第j跳回波聚集区数据；筛选出所有满足条件的杂波中心距离元dj，对应区域为有效聚集区，执行步骤六；步骤六：由步骤五筛选得到的每个满足条件的杂波中心距离元dj，对应的雷达接收回波谱信号记为n＝1,…,N，其中N为采样点数，dj为有效聚集区的杂波中心距离元；使用的上下相邻r个距离元信号构造聚集区信号矩阵使用杂波中心距离元d0对应的信号的上下相邻r个距离元信号构造聚集区信号矩阵r的取值由聚集区内距离元的个数确定(r可以为二分之一的距离元个数)；步骤七：对应每个满足条件的dj，由步骤六得到相应的聚集区信号矩阵取和每个构造多聚集区联合参考矩阵X(和分别取第1至2r+1列)；步骤八：根据步骤七得到的多聚集区联合参考矩阵X，求得X的相关矩阵的最大似然估计R；步骤九：对步骤八求得的最大似然估计R进行特征值分解，得到降序排列的特征值λs和与λs相对应的特征向量vs，s为特征值的个数；步骤十：在R中，电离层杂波的能量所占比重较大，在步骤十得到的特征向量vs中，认为前L个大的特征值所对应的特征向量为电离层杂波的特征向量，取前L(L为经验值)个大的特征值所对应的特征向量为电离层杂波的特征向量，构成矩阵，形成电离层杂波子空间特征向量矩阵V；步骤十一：利用步骤十得到的特征向量矩阵V，构造正交投影矩阵P；步骤十二：通过将步骤七得到的多聚集区联合参考矩阵X，投影到由步骤十一确定的正交投影矩阵P的正交子空间上，来消除待处理距离元上与多聚集区联合参考矩阵相关的电离层杂波，求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y。本专利技术的有益效果为：本专利技术的目的是提供一种基于多聚集区距离相关性的高频地波雷达电离层杂波抑制方法，以解决现有
技术介绍
中存在的问题，更好的实现电离层杂波抑制，提高高频地波雷达的探测性能。本专利技术相对传统杂波抑制算法，利用电离层多聚集区杂波之间的互相关性来完成对消处理，在抑制电离层杂波的同时，保留海浪回波及一阶Bragg峰的信息，减少正常海洋回波能量的下降。经过本专利技术方法处理后，雷达回波谱中的强电离层杂波被明显抑制，杂波抑制幅度提高，回波信杂噪比获得显著改善。实施例一中定量分析本专利技术的抑制效果，以两个一阶Bragg峰为信号，信号以外的其余部分为杂波和噪声基底，分别计算算法处理前后的信杂噪比，计算两者差值可以得到杂波抑制幅度，它反映了算法的真实性能。取121.5km处的回波数据计算，得到处理前回波谱的信噪比为14.74dB，处理后回波谱的信噪比为30.65dB，信噪比改善幅度提高15.91dB，杂波抑制的效果明显。实施例二中定量分析本专利技术的抑制效果，分别计算算法处理前后的信杂噪比。取119km本文档来自技高网...

【技术保护点】
1.一种高频地波雷达电离层杂波抑制方法，其特征在于：所述高频地波雷达电离层杂波抑制方法包括以下步骤：步骤一：由雷达距离‑多普勒谱，确定对应每个电离层杂波聚集区的杂波中心距离元ki,i＝0,1,…,p‑1，p为电离层杂波聚集区的个数；步骤二：根据步骤一确定的杂波中心距离元ki，计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关系数

【技术特征摘要】
1.一种高频地波雷达电离层杂波抑制方法，其特征在于：所述高频地波雷达电离层杂波抑制方法包括以下步骤：步骤一：由雷达距离-多普勒谱，确定对应每个电离层杂波聚集区的杂波中心距离元ki,i＝0,1,…,p-1，p为电离层杂波聚集区的个数；步骤二：根据步骤一确定的杂波中心距离元ki，计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关系数步骤三：设置电离层杂波聚集区的相关系数参考阈值Th1，比较由步骤二得到的距离相关系数的值与参考阈值Th1的大小关系，如果则判定该杂波中心距离元对应杂波区域为聚集区，数据有效；否则，判定该杂波中心距离元对应杂波区域不是聚集区，数据无效，舍弃；筛选出所有满足条件的杂波中心距离元ki，按从小到大的顺序依次排列并赋值给dc,c＝0,1,…,m-1，表示聚集区的杂波中心距离元，m为有效数据个数，执行步骤四；步骤四：计算杂波中心距离元d0与第j跳回波聚集区的杂波中心距离元dj的互相关系数步骤五：设置电离层杂波聚集区的互相关系数参考阈值Th2，将步骤四求得的每个互相关系数与参考阈值Th2作比较，如果则认为第j跳回波聚集区数据有效；否则，舍弃第j跳回波聚集区数据；筛选出所有满足条件的杂波中心距离元dj，对应区域为有效聚集区，执行步骤六；步骤六：由步骤五筛选得到的每个满足条件的杂波中心距离元dj，对应的雷达接收回波谱信号记为其中N为采样点数，dj为有效聚集区的杂波中心距离元；使用的上下相邻r个距离元信号构造聚集区信号矩阵使用杂波中心距离元d0对应的雷达接收回波谱信号的上下相邻r个距离元信号构造聚集区信号矩阵步骤七：取和构造多聚集区联合参考矩阵X；步骤八：根据步骤七得到的多聚集区联合参考矩阵X，求得X的相关矩阵的最大似然估计R；步骤九：对步骤八求得的最大似然估计R进行特征值分解，得到降序排列的特征值λs和与λs相对应的特征向量vs，s为特征值的个数；步骤十：取前L个大的特征值所对应的特征向量为电离层杂波的特征向量，构成矩阵，形成电离层杂波子空间特征向量矩阵V；步骤十一：利用步骤十得到的特征向量矩阵V，构造正交投影矩阵P；步骤十二：通过将步骤七得到的多聚集区联合参考矩阵X，投影到由步骤十一确定的正交投影矩阵P的正交子空间上，来消除待处理距离元上与多聚集区联合参考矩阵相关的电离层杂波，求出待处理聚集区经过正交投影处理后的杂波抑制信号矩阵Y。2.根据权利要求1所述一种高频地波雷达电离层杂波抑制方法，其特征在于：所述步骤二中计算电离层杂波在杂波中心距离元ki和与其相邻的距离元ki+1上的距离相关...

【专利技术属性】
技术研发人员：董英凝，董博远，邓维波，赵彬，
申请(专利权)人：哈尔滨工业大学，
类型：发明
国别省市：黑龙江,23