[发明专利]外辐射源雷达直达波恢复方法

说明书

技术领域

本发明涉及雷达技术领域，尤其涉及一种外辐射源雷达直达波恢复方法。

背景技术

直达波恢复是调频广播外辐射源雷达的关键技术。随着系统在城市附近或山区环境的应用，严重的环境杂波使得外辐射源雷达直达波恢复变得非常困难。传统的直达波恢复采用波束形成的主瓣对准直达波信号方向，利用直达波天线副瓣实现对多径信号的抑制。在严重多径情况下，多径杂波从主瓣进入，造成信号失真。采用常数模盲均衡技术同样可以实现直达波提纯，包括同时进入主瓣的多径情况。在严重多径的情况下，常数模方法的收敛性能较差。而空间分集盲均衡方法可以适应严重多径的情况，增强直达波的恢复能力。但是随着空间分集维数的增加，算法的计算量大大增加，不利于该方法的工程实现。

发明内容

本发明所要解决的技术问题在于，提供一种外辐射源雷达直达波恢复方法，该方法针对空间分集的盲均衡方法，设计了计算量较小的间隔迭代方法，实现了强多径杂波环境下直达波的恢复。

为了解决上述技术问题，本发明提供了一种外辐射源雷达直达波恢复方法包括：

S101：外辐射源雷达的参考通道采用空间分集方法接收多个分集通道的分集信号；

S102：针对每一分集通道的分集信号采用一路均衡器进行处理，从而得到多路均衡器输出信号；其中，式中，p_i(k)表示第i路均衡器输出信号，x_i(k)表示当前时刻第i个分集通道的分集信号，f_i(k)为N_f×1向量，表示第i路均衡器的向量，f_i^l(k)表示第i路均衡器的第l个均衡器系数，i＝0,1，…,N-1，N表示均衡器的个数；

S103：将多路均衡器输出信号进行求和，得到直达波信号；其中，z(k)表示直达波信号；

S104：根据所述直达波信号计算得到误差项；其中，e(k)＝z(k)[γ-z²(k)]，式中，e(k)表示误差项，γ表示外辐射源雷达发射的原始信号a(k)的峰度，且

S105：根据所述误差项计算得到均方误差，并判断所述均方误差是否小于预设门限；其中，所述均方误差为σ²(k+1)＝λσ²(k)+e(k)e^*(k)，式中，λ为遗忘因子，初始值σ²(0)＝0；

S106：若所述均方误差不小于所述预设门限，则返回S104计算下一误差项e(k+1)；

S107：若所述均方误差小于所述预设门限，则每隔N_g个采样点，对每一路均衡器的向量进行迭代；其中，向量迭代公式为f_i(k+1)＝f_i(k)+μ·X_i^*(k)e(k)，式中，X_i(k)＝[x_i(k) x_i(k-1) … x_i(k-l) … x_i(k-N_f+1)]^T，X_i(k)表示第i个分集通道的输入向量，x_i(k-l)表示第i个分集通道的第l个分集信号，μ为步长因子。

实施本发明，具有如下有益效果：本发明的外辐射源雷达直达波恢复方法针对空间分集的盲均衡方法，设计了计算量较小的间隔迭代方法，实现了强多径杂波环境下直达波的恢复，易于工程实现。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明提供的外辐射源雷达直达波恢复方法的一个实施例的流程示意图；

图2是图1的原理框图；

图3是对本发明提供的外辐射源雷达直达波恢复方法的仿真验证方法的流程示意图；

图4是图3的原理框图；

图5是S201的原理框图；

图6是SD-CMA算法输出信号模值；

图7是SD-CMA算法距离多普勒平面的示意图；

图8是间隔为2的IISD-CMA算法输出信号模值；

图9是间隔为2的IISD-CMA算法的距离多普勒图