外辐射源雷达直达波恢复方法
【技术领域】
[0001] 本发明涉及雷达技术领域，尤其涉及一种外辐射源雷达直达波恢复方法。
【背景技术】
[0002] 直达波恢复是调频广播外辐射源雷达的关键技术。随着系统在城市附近或山区 环境的应用，严重的环境杂波使得外辐射源雷达直达波恢复变得非常困难。传统的直达波 恢复采用波束形成的主瓣对准直达波信号方向，利用直达波天线副瓣实现对多径信号的抑 制。在严重多径情况下，多径杂波从主瓣进入，造成信号失真。采用常数模盲均衡技术同样 可以实现直达波提纯，包括同时进入主瓣的多径情况。在严重多径的情况下，常数模方法的 收敛性能较差。而空间分集盲均衡方法可以适应严重多径的情况，增强直达波的恢复能力。 但是随着空间分集维数的增加，算法的计算量大大增加，不利于该方法的工程实现。
【发明内容】
[0003] 本发明所要解决的技术问题在于，提供一种外辐射源雷达直达波恢复方法，该方 法针对空间分集的盲均衡方法，设计了计算量较小的间隔迭代方法，实现了强多径杂波环 境下直达波的恢复。
[0004] 为了解决上述技术问题，本发明提供了一种外辐射源雷达直达波恢复方法包括：
[0005] SlOl :外辐射源雷达的参考通道采用空间分集方法接收多个分集通道的分集信 号；
[0006] S102:针对每一分集通道的分集信号采用一路均衡器进行处理，从而得 到多路均衡器输出信号；其中，朽(々）二式中，Pi( k)表示第i路均 衡器输出信号，X1GO表示当前时刻第i个分集通道的分集信号，fjk)为N fXl向 量，表示第i路均衡器的向量，f/ (k)表示第i路均衡器的第1个均衡器系数， …，...........
........... i = 0, 1，·，·，Ν-1，Ν表示均衡器的 个数；
[0007]S103 :将多路均衡器输出信号进行求和，得到直达波信号；其中，
ζ (k)表不直达波信号；
[0008] S104 :根据所述直达波信号计算得到误差项；其中，e (k) = z (k) [ γ -Z2 (k)]， 式中，e(k)表示误差项，γ表示外辐射源雷达发射的原始信号a(k)的峰度，且
[0009] S105:根据所述误差项计算得到均方误差，并判断所述均方误差是否小于预设门 限；其中，所述均方误差为〇2(k+l) = λ σ2(1〇+θ(1〇#(1〇,式中，λ为遗忘因子，初始值 σ 2(O) = O ;
[0010] S106 :若所述均方误差不小于所述预设门限，则返回S104计算下一误差项 e(k+l)；
[0011] S107 :若所述均方误差小于所述预设门限，则每隔N，采样点，对每一路均衡器 的向量进行迭代；其中，向量迭代公式为fi(k+l) =fi(k) + y ·Χ广(k)e(k)，式中，XiQO =[XiGO Xi(Ic-I)…Xi(Ic-I)…Xi (k-Nf+l)]T，XiQO表示第i个分集通道的输入向量， X1 (k-Ι)表示第i个分集通道的第1个分集信号，μ为步长因子。
[0012] 实施本发明，具有如下有益效果：本发明的外辐射源雷达直达波恢复方法针对空 间分集的盲均衡方法，设计了计算量较小的间隔迭代方法，实现了强多径杂波环境下直达 波的恢复，易于工程实现。
【附图说明】
[0013] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以 根据这些附图获得其他的附图。
[0014] 图1是本发明提供的外辐射源雷达直达波恢复方法的一个实施例的流程示意图；
[0015] 图2是图1的原理框图；
[0016] 图3是对本发明提供的外辐射源雷达直达波恢复方法的仿真验证方法的流程示 意图；
[0017] 图4是图3的原理框图；
[0018] 图5是S201的原理框图；
[0019] 图6是SD-CM算法输出信号模值；
[0020] 图7是SD-CMA算法距离多普勒平面的示意图；
[0021] 图8是间隔为2的IISD-CM算法输出信号模值；
[0022] 图9是间隔为2的IISD-CM算法的距离多普勒图
[0023] 图10是间隔为3的IISD-CM算法输出信号模值；
[0024] 图11是间隔为3的IISD-CM算法的距离多普勒图。
【具体实施方式】
[0025] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0026] 图1是本发明提供的外辐射源雷达直达波恢复方法的一个实施例的流程示意图， 图2是图1的原理框图，如图1和图2所示，包括：
[0027] SlOl :外辐射源雷达的参考通道采用空间分集方法接收多个分集通道的分集信 号；
[0028] S102:针对每一分集通道的分集信号采用一路均衡器进行处理，从而得 到多路均衡器输出信号；其中，尺⑷=Ii ，式中，P1GO表示第i路均 衡器输出信号，X1GO表示当前时刻第i个分集通道的分集信号，fjk)为NfXl向 量，表示第i路均衡器的向量，f/ (k)表示第i路均衡器的第1个均衡器系数，
i=0，l，"·，Ν-1，Ν表示均衡器的 个数；
[0029] S103 :将多路均衡器输出信号进行求和，得到直达波信号；其中，
(k)表不直达波信号； _ . I-U
[0030] S104 :根据所述直达波信号计算得到误差项；其中，e (k) = Z (k) [ γ -Z2 (k)]， 式中，e(k)表示误差项，γ表示外辐射源雷达发射的原始信号a(k)的峰度，且
[0031] S105:根据所述误差项计算得到均方误差，并判断所述均方误差是否小于预设门 限；其中，所述均方误差为〇2(k+l) = λ σ2(1〇+θ(1〇#(1〇,式中，λ为遗忘因子，初始值 σ 2(〇) = 〇 ；
[0032] S106 :若所述均方误差不小于所述预设门限，则返回S104计算下一误差项 e(k+l)；
[0033] S107 :若所述均方误差小于所述预设门限，则每隔义个采样点，对每一路均衡器 的向量进行迭代；其中，向量迭代公式为fi(k+l) =fi(k) + y ·Χ广(k)e(k)，式中，XiQO =[XiGO Xi(Ic-I)…Xi(Ic-I)…Xi (k-Nf+l)]T，XiQO表示第i个分集通道的输入向量， X1 (k-Ι)表示第i个分集通道的第1个分集信号，μ为步长因子。
[0034] 下面通过计算机仿真信号验证本发明实施例的正确性：
[0035] 如图3和图4所示，仿真验证方法的步骤包括：
[0036] S201 :外辐射源雷达通过空间分集方法将原始信号发送至多个分集通道。
[0037] 具体的，如图5所示，所述S201具体包括：
[0038] 外辐射源雷达通过空间分集方法将原始信号发送至多个分集通道，其中，第i个 分集通道上的分集信号为= + ，式中，a(k)为外辐射源雷达发 射的原始信号，H1GO为第i个分集通道接收端的噪声，hjk)为kXl信道响应向量，且
?示第i个分集通道第1个 信道系数。
[0039] S202:采用本发明提供的外辐射源雷达直达波恢复方法恢复出直达波信号。
[0040] S203:根据监测通道信号和所述直达波信号采用最小均方算法进行杂波对消，得 到对消后的信号。
[0041] 具体的，所述S203具体包括：
[0042] 产生监测通道信号y (k)，其中，(幻=C ? ￡/(幻十g(々）+ /?、.(々），式中，ns (k)为接 收机噪声，c为多径杂波，且C = [f