一种基于黎曼距离的广义旁瓣对消算法制造技术

本发明专利技术为一种基于黎曼距离的广义旁瓣对消算法，涉及一种用于混合传播模式高频超视距雷达的海杂波抑制处理方法。为解决混合传播模式下，高频超视距雷达接收到的非平稳海杂波的抑制问题，本发明专利技术提供一种基于黎曼距离的广义旁瓣对消处理方法。本发明专利技术在传统广义旁瓣对消处理方法基础上，基于训练样本协方差矩阵之间的黎曼距离，挑选有效的训练样本，进行非平稳杂波抑制处理。本发明专利技术用于混合传播模式高频超视距雷达对海杂波以及其他非平稳性杂波进行抑制，以提高信杂噪比，增大目标发现概率。

A generalized sidelobe cancellation algorithm based on Riemann distance

【技术实现步骤摘要】
一种基于黎曼距离的广义旁瓣对消算法
本专利技术涉及一种用于高频超视距雷达杂波抑制的空域自适应处理算法，属于混合传播模式高频超视距雷达杂波抑制领域。
技术介绍
混合传播模式高频超视距雷达是一种新体制高频雷达系统。其在传统单站高频地波超视距雷达的基础上，丰富了电磁波传播路径，结合了地波海面传播和天波沿电离层传播这两种传播模式。其发射信号既有沿海面传播的地波模式，也有经电离层反射的天波模式，电磁波照射到目标后，回波信号有通过地波模式，沿海面传播被接收，也有反射入电离层，由天波模式被接收。这种新的高频超视距雷达体制在近几年得到了较多的关注，其结合了天波雷达的长探测距离和地波雷达高探测性能的优点，同时提高了生存能力。但其也不可避免的引入了两种体制超视距雷达的缺点，而其中一个最严重的问题是非平稳环境下的混合杂波，包括电离层杂波和被电离层污染的展宽的一阶海杂波。在多普勒维和方位维扩展的杂波，会将在Bragg峰附近的低速运动目标淹没，造成目标无法被探测。广义旁瓣对消算法是一种成功被应用在高频超视距雷达领域的空域自适应杂波处理算法。方法的核心在于利用方位扩展杂波(海杂波)的全向特性，利用旁瓣内的杂波信息，估计主瓣杂波信息，进行对消处理，到达杂波抑制的目的。同时，其利用单凹口滤波器滤除主瓣内的目标信息，在保护目标不被对消的同时，实现杂波抑制。传统的广义旁瓣对消算法，对于训练样本，要求其和待测杂波样本是独立同分布的，但这在混合传播模式高频超视距雷达中难以满足。因为与传统地波超视距雷达的单一海面传播路径相比，混合传播模式引入了天波模式的电离层信道传播，导致回波信号中的海杂波可能被电离层杂波污染，导致其特性非平稳。同时，电离层杂波也会通过接收机被接收。因而其背景环境的非平稳性会导致对杂波协方差矩阵的估计误差，使得杂波抑制算法性能下降。因此在混合传播模式高频超视距雷达背景下，利用广义旁瓣对消算法进行杂波抑制处理时，如何挑选有效的训练样本是一个亟待解决的难点。
技术实现思路
本专利技术的目的是为了解决在混合传播模式高频超视距雷达背景下的非平稳杂波抑制问题，本专利技术提供一种基于黎曼距离的广义旁瓣对消处理算法。本专利技术的基于黎曼距离的广义旁瓣对消算法，包括如下步骤：步骤一：对混合传播模式高频超视距雷达的所有通道回波信号进行距离处理、多普勒处理，得到距离-速度-通道三维数据块；步骤二：对各个通道数据块进行数字波束形成处理，构成主波束数据块，在指定的角度单元，对各个通道数据块进行单凹口滤波器辅助波束形成处理，构成辅助波束数据块；步骤三：选定训练样本挑选单元大小，选取指定多普勒单元和角度单元的若干距离单元主波束数据，构成训练样本挑选数据块，计算每个距离单元对应的待处理局域单元的协方差矩阵，构成协方差数据块；步骤四：计算指定距离单元的协方差矩阵与若干其他距离单元的协方差矩阵之间的黎曼距离，构成黎曼距离数组，从小至大挑选其中若干个黎曼距离对应的距离单元作为距离维训练样本单元；步骤五：根据选定的训练样本数据，计算指定距离单元的自适应权向量，以及辅助波束的自适应输出，将主波束与辅助波束的输出做差，得到广义旁瓣对消处理后的输出结果；步骤六：遍历所有感兴趣距离单元、多普勒单元和角度单元，得到经过杂波抑制处理后的距离-速度-角度三维数据结果。优选的是，所述步骤一中：对混合传播模式高频超视距雷达的所有通道回波信号进行距离处理、多普勒处理，得到距离-速度-通道三维数据块为{dataCh}，其维度为R×D×Ch，其中，R为数据的距离维度，D为数据的速度维度，Ch为通道数。优选的是，所述步骤二中，对各个通道数据块进行数字波束形成处理，构成主波束数据块，在指定的角度单元，对各个通道数据块进行单凹口滤波器辅助波束形成处理，构成辅助波束数据块的方法为：对各个通道数据块{dataCh}进行数字波束形成处理，得到距离-速度-角度三维主波束数据块为{dataDBF}，其维度为R×D×A，其中A为数据的角度维度；对{dataCh}在指定角度单元θ，进行单凹口滤波器辅助波束形成处理，将Ch元线阵分为Ch′个子阵，每个子阵包含的阵元数为Ch-Ch′+1，则所需单凹口滤波器的阶数为Ch-Ch′，在θ角度的单凹口滤波器系数为H(θ)，辅助波束数据块为{dataCBF}，其中dataCBF＝dataCh×H(θ)，其维度为R×D×(Ch-Ch′+1)。优选的是，所述步骤三中，选定训练样本挑选单元大小，选取指定多普勒单元和角度单元的若干距离单元主波束数据，构成训练样本挑选数据块，计算每个距离单元对应的待处理局域单元的协方差矩阵，构成协方差数据块的方法为：待挑选的训练样本单元所包含的多普勒单元数为localD，角度单元数为localA，则对指定多普勒单元d和指定角度单元θ，选取setR个距离单元构成的训练样本挑选数据块为{dataL}，其维度为setR×localD×localA，其中，setR代表要挑选的距离单元数目；计算某一距离单元r对应的待处理局域单元的协方差的方法为先将dataLr进行列向量化，得到一个一维列向量vr，其长度为localD×localA；然后计算所述一维列向量vr对应的协方差矩阵为其中·H表示共轭转置，dataLr为{dataL}中，第r个距离单元的数据；对每一个距离单元进行计算，可以得到协方差数据块为{Rd×θ}，其维度为R×(localD×localA)×(localD×localA)。优选的是，所述步骤四中，计算指定距离单元的协方差矩阵与若干其他距离单元的协方差矩阵之间的黎曼距离，构成黎曼距离数组，从小至大挑选其中若干个黎曼距离对应的距离单元作为距离维训练样本单元的方法为：选定距离单元r，其对应的协方差矩阵为Rr，依次计算其与其他训练样本距离单元对应的协方差矩阵Rj,j∈{r-setR/2，r+setR/2}，之间的黎曼距离为：其中，||·||F为Frobenius范数。所要挑选的距离维训练样本数为nSelect，则从RieDis(j)中从小至大选择nSelect个数据对应的距离单元数组为Rselect(k)，其中k＝1,2,......,nSelect。优选的是，所述步骤五中，根据选定的训练样本数据，计算指定距离单元的自适应权向量，以及辅助波束的自适应输出，将主波束与辅助波束的输出做差，得到广义旁瓣对消处理后的输出结果的方法为：对于指定距离单元r，指定多普勒单元d和指定角度单元θ，由选定的训练样本数据得到的主波束输出为DBFout(Rselect,1)＝dataDBF(Rselect,d,θ)，辅助波束输出为CBFout(Ch-Ch′+1,Rselect)＝dataCBF(Rselect,d,Ch-Ch′+1)；辅助波束的自相关矩阵为RCBF＝CBFout·CBFoutH，主波束与辅助波束的互相关矩阵为rDBF＝CBFout·DBFoutH，则指定距离单元的自适应权向量可表示为在选定距离单元r，多普勒本文档来自技高网...

【技术保护点】
1.一种基于黎曼距离的广义旁瓣对消算法，其特征在于，所述方法包括如下步骤：/n步骤一：对混合传播模式高频超视距雷达的所有通道回波信号进行距离处理、多普勒处理，得到距离-速度-通道三维数据块；/n步骤二：对各个通道数据块进行数字波束形成处理，构成主波束数据块，在指定的角度单元，对各个通道数据块进行单凹口滤波器辅助波束形成处理，构成辅助波束数据块；/n步骤三：选定训练样本挑选单元大小，选取指定多普勒单元和角度单元的若干距离单元主波束数据，构成训练样本挑选数据块，计算每个距离单元对应的待处理局域单元的协方差矩阵，构成协方差数据块；/n步骤四：计算指定距离单元的协方差矩阵与若干其他距离单元的协方差矩阵之间的黎曼距离，构成黎曼距离数组，从小至大挑选其中若干个黎曼距离对应的距离单元作为距离维训练样本单元；/n步骤五：根据选定的训练样本数据，计算指定距离单元的自适应权向量，以及辅助波束的自适应输出，将主波束与辅助波束的输出做差，得到广义旁瓣对消处理后的输出结果；/n步骤六：遍历所有感兴趣距离单元、多普勒单元和角度单元，得到经过杂波抑制处理后的距离-速度-角度三维数据结果。/n

【技术特征摘要】
1.一种基于黎曼距离的广义旁瓣对消算法，其特征在于，所述方法包括如下步骤：
步骤一：对混合传播模式高频超视距雷达的所有通道回波信号进行距离处理、多普勒处理，得到距离-速度-通道三维数据块；
步骤二：对各个通道数据块进行数字波束形成处理，构成主波束数据块，在指定的角度单元，对各个通道数据块进行单凹口滤波器辅助波束形成处理，构成辅助波束数据块；
步骤三：选定训练样本挑选单元大小，选取指定多普勒单元和角度单元的若干距离单元主波束数据，构成训练样本挑选数据块，计算每个距离单元对应的待处理局域单元的协方差矩阵，构成协方差数据块；
步骤四：计算指定距离单元的协方差矩阵与若干其他距离单元的协方差矩阵之间的黎曼距离，构成黎曼距离数组，从小至大挑选其中若干个黎曼距离对应的距离单元作为距离维训练样本单元；
步骤五：根据选定的训练样本数据，计算指定距离单元的自适应权向量，以及辅助波束的自适应输出，将主波束与辅助波束的输出做差，得到广义旁瓣对消处理后的输出结果；
步骤六：遍历所有感兴趣距离单元、多普勒单元和角度单元，得到经过杂波抑制处理后的距离-速度-角度三维数据结果。


2.根据权利要求1所述的基于黎曼距离的广义旁瓣对消算法，其特征在于，所述步骤一中，对混合传播模式高频超视距雷达的所有通道回波信号进行距离处理、多普勒处理，得到距离-速度-通道三维数据块为{dataCh}，其维度为R×D×Ch，其中，R为数据的距离维度，D为数据的速度维度，Ch为通道数。


3.根据权利要求2所述的基于黎曼距离的广义旁瓣对消算法，其特征在于，所述步骤二中，对各个通道数据块进行数字波束形成处理，构成主波束数据块，在指定的角度单元，对各个通道数据块进行单凹口滤波器辅助波束形成处理，构成辅助波束数据块的方法为：
对各个通道数据块{dataCh}进行数字波束形成处理，得到距离-速度-角度三维主波束数据块为{dataDBF}，其维度为R×D×A，其中A为数据的角度维度；
对{dataCh}在指定角度单元θ，进行单凹口滤波器辅助波束形成处理，将Ch元线阵分为Ch′个子阵，每个子阵包含的阵元数为Ch-Ch′+1，则所需单凹口滤波器的阶数为Ch-Ch′，在θ角度的单凹口滤波器系数为H(θ)，辅助波束数据块为{dataCBF}，其中dataCBF＝dataCh×H(θ)，其维度为R×D×(Ch-Ch′+1)。


4.根据权利要求3所述的基于黎曼距离的广义旁瓣对消算法，其特征在于，所述步骤三中，选定训练样本挑选单元大小，选取指定多普勒单元和角度单元的若干距离单元主波束数据，构成训练样本挑选数据块，计算每个距离单元对应的待处理局域单元的协方差矩阵，构成协方差数据块的方法为：
待挑选的训练样本单元所包含的多普勒单元数为localD，角度单元数为localA，则对指定多普勒单元d和指定角度单元θ，选取setR个距离单元构成的训练样本挑选数据块为{dataL}，其维度为setR×localD×localA，其中，setR代表要挑选的距离单元数目；
计算某一距离...

【专利技术属性】
技术研发人员：杨强，张佳智，杨勇，张鑫，张薇，赵梦晓，李佳明，
申请(专利权)人：杨强，张佳智，
类型：发明
国别省市：黑龙;23