LFM-频率编码复合低截获波形信号处理方法

罗美方, 郝志梅, 王强

罗美方, 郝志梅, 王强. LFM-频率编码复合低截获波形信号处理方法[J]. 电子科技大学学报, 2017, 46(2): 335-339. DOI: 10.3969/j.issn.1001-0548.2017.02.003
引用本文: 罗美方, 郝志梅, 王强. LFM-频率编码复合低截获波形信号处理方法[J]. 电子科技大学学报, 2017, 46(2): 335-339. DOI: 10.3969/j.issn.1001-0548.2017.02.003
LUO Mei-fang, HAO Zhi-mei, WANG Qiang. Signal Processing Technique for LFM/FSK Hybrid LPI Waveform[J]. Journal of University of Electronic Science and Technology of China, 2017, 46(2): 335-339. DOI: 10.3969/j.issn.1001-0548.2017.02.003
Citation: LUO Mei-fang, HAO Zhi-mei, WANG Qiang. Signal Processing Technique for LFM/FSK Hybrid LPI Waveform[J]. Journal of University of Electronic Science and Technology of China, 2017, 46(2): 335-339. DOI: 10.3969/j.issn.1001-0548.2017.02.003

LFM-频率编码复合低截获波形信号处理方法

基金项目: 

航空科学基金 2014ZC07002

详细信息
    作者简介:

    罗美方 (1984-), 女, 高级工程师, 主要从事雷达低截获技术方面的研究

  • 中图分类号: TN95

Signal Processing Technique for LFM/FSK Hybrid LPI Waveform

  • 摘要: LFM-频率编码复合低截获波形把LFM技术和FSK技术组合在一起,获得大的时带积,从而增强了雷达的低截获性能,但LFM-频率编码复合低截获波形的相关距离旁瓣较高,影响雷达检测性能。针对这一问题,该文对旁瓣抑制方法展开了算法研究和工程应用设计,通过离线设计旁瓣抑制滤波器,代替传统匹配滤波或相关处理,从旁瓣抑制滤波器库中选取相应的旁瓣抑制滤波器应用于回波信号采样数据,完成脉压的同时抑制距离旁瓣,实现较为简单,适合工程化应用。仿真结果表明了该设计的有效性,旁瓣性能得到改善。
    Abstract: Hybrid LFM-frequency encoding low probability of interception (LPI) waveform combines the techniques of LFM and FSK, which can obtain a large time-bandwith product, thus enhancing the capability of LPI. However, hybrid LFM/FSK LPI signal suffers from high range sidelobes, resulting in the decrease of radar detectability. To this end, this paper deals with the problem of sidelobe suppression. Specifically, the specific sidelobe suppression filter selected from its library which is designed off-line, instead of matched filter or correlation, is exploited on the echo signal. This operation can accomplish the compression pulse while achieving low sidelobes, and its engineering implementation is simple. Simulation results illustrate the effectiveness of the devised filter, showing its ability to obtain lower sidelobes in comparison with the traditional process.
  • 从低截获概率 (LPI) 雷达理论来说,低截获概率雷达对于信号的要求是峰值功率低、频带宽、调制形式复杂、具有大时宽带宽积等[1-2]。如何通过各种调制技术和方法来有效地降低雷达信号时频域的峰值功率和截获因子是低截获概率雷达波形设计的重点。采用大时宽带宽积信号可以在保证雷达威力和距离分辨力的同时使雷达辐射的峰值功率大大降低,因此,大时宽带宽积信号波形设计是解决途径之一。

    目前,常用的大时宽带宽积信号有线性调频 (LFM)、相位编码、频率编码等均为单一调制信号:频率调制、相位调制,同时也各自存在不足。为了提高雷达信号的抗截获性能,各种混合调制波形逐渐得到了广泛的应用。文献[3]提出了FSK/PSK技术,把FSK和PSK调制技术组合在一起,使用频率分量产生了瞬时扩展,能够获得大时宽带宽积[3-5]。在实际工程中,要在极短的码宽内产生较长的相位跳变具有较大的难度,本文的LFM-频率编码复合信号将频率编码技术和LFM技术组合在一起:子脉冲内线性调频、子脉冲间频率跳变,相对易于实现;该信号采用组合调制的方法,增加波形结构的复杂度,提高信号的随机程度,从而增加截获接收机的截获、识别难度。

    LFM-频率编码信号具有较低的截获概率,在现代LPI雷达系统中具有良好的应用前景,但同时也存在着较高的距离旁瓣问题。因此,需要重点研究复合波形的信号处理方法,有效抑制脉压距离旁瓣,提升复合波形的性能。

    抑制旁瓣的有效方法是对信号进行加权 (加窗),但是该方法是以信噪比损失及距离分辨力变坏为代价。

    基于最小二乘的失配滤波器主要应用于相位编码脉压副瓣的优化[6-7]。本文结合实际应用条件,对传统失配滤波器进行改进,提出一种快速脉压距离旁瓣抑制实现方法,通过预先离线获取旁瓣抑制滤波器、建立旁瓣抑制滤波器库,并依据雷达参数从滤波器库选取处理滤波器,在线实时处理目标回波信号,有效抑制距离旁瓣。

    LFM-频率编码复合信号的每个子脉冲具有相同斜率的LFM信号,每个子脉冲间进行频率编码,LFM-频率编码复合调制信号的数学表达式为:

    $$ \begin{array}{l} s(t) = \sum\limits_{m = 0}^{M - 1} {A{u_c}(t - m\tau )} \times \\ \sum\limits_{m = 0}^{M - 1} {\exp [{\rm{j}}2{\rm{\pi }}{c_m}\Delta ft + {\rm{j\pi }}k{{(t - m\tau )}^2}]} \end{array} $$ (1)

    式中,A为信号幅度;${u_c}(t) = \left\{ \begin{array}{l} 1\;\;\;\;\;0 < t \le \tau \\ 0\;\;\;\;\;其他 \end{array} \right.$,$\tau $为码元宽度;Δf为跳频间隔;M为码长;cm(m=0, 1, …M-1) 为频率编码序列;$k{\rm{ = }}{B_{{\rm{LFM}}}}/\tau $为调频斜率;${B_{{\rm{LFM}}}}$为线性调频信号的带宽。脉宽$T = M\tau $,信号的带宽可以表示为:

    $$ B = (M - 1)\Delta f + {B_{{\rm{LFM}}}} $$ (2)

    LFM-频率编码复合波形具有较高的距离旁瓣,可能影响雷达对目标的探测,必须抑制旁瓣降低影响。

    最小二乘 (least square, LS) 失配滤波器是一种常用的直接失配旁瓣抑制方法,以旁瓣电平最小化为目标函数,从滤波器开始运行直到当前时间,滤波器系数的更新总是令总的平方误差趋于最小[6-14]

    对于长度为N的${\boldsymbol{s}} = {[s(0),s(1), \cdots ,s(N - 1)]^{\rm{T}}}$,假设其失配滤波器的阶数为P($P \ge N$),首先为了使信号长度为P,同时使得输出具有对称性,将编码信号两端补零,则信号矢量S和滤波器系数矢量h分别表示为:

    $$ {{\boldsymbol{S}}^{\rm{T}}} = \left[ {0, \cdots 0,s(0),s(1), \cdots ,s(N - 1),0, \cdots ,0} \right] $$ (3)
    $$ {{\boldsymbol{h}}^{\rm{T}}} = \left[ {h(0),h(1), \cdots ,h(P - 1)} \right] $$ (4)

    脉压滤波器输出为:

    $$ y(m) = \sum\limits_{n = 0}^N {s(n){h^*}(n - m)} $$ (5)

    式中,$m = - (P - 1), \cdots ,P - 1$。

    式 (5) 可以进一步表示为:

    $$ {\boldsymbol{y}} = {{\boldsymbol{h}}^{\rm{H}}}{\boldsymbol{X}} $$ (6)

    式中,hH表示h的共轭转置;XP×(2P-1) 维的系数矩阵:

    $$ {\boldsymbol{X}} = \left( {\begin{array}{*{20}{c}} 0& \cdots &0&{x(0)}&{}&{\begin{array}{*{20}{c}} \cdots &{x(P - 1)} \end{array}}\\ 0&{}&{}& \cdots &{\begin{array}{*{20}{c}} {}&{} \end{array}}&0\\ 0&{x(0)}&{x(1)}& \cdots &{\begin{array}{*{20}{c}} {x(P - 1)}&0 \end{array}}& \cdots \\ {x(0)}&{x(1)}& \cdots &{x(P - 1)}&{\begin{array}{*{20}{c}} 0& \cdots \end{array}}&0 \end{array}} \right) $$ (7)

    脉压处理的理想结果是除主瓣外其余输出全为零,实际情况是不可能的,因此,采取使旁瓣尽可能小的逼近设计。通常,失配滤波器的设计时考虑旁瓣电平和信噪比损失 (增益处理损失LPG),其中,旁瓣电平有峰值旁瓣电平PSL和积分旁瓣电平ISL的计算公式分别为:

    $$ {\rm{PSI}} = 20\log 10\left| {\frac{{{P_s}}}{{{P_{\max }}}}} \right| $$ (8)
    $$ {\rm{ISL}} = 20 \times \log 10\frac{{\sum\limits_i {\left| {{P_{s,i}}} \right|} }}{{\left| {{P_{\max }}} \right|}} $$ (9)
    $$ {\rm{LSNR}} = 20\log 10({P_1}) - 20\log 10({P_2}) $$ (10)

    式中,PmaxPs分别为主副瓣峰值;P1P2为失配滤波前后的主瓣峰值。

    从应用的角度看,旁瓣抑制希望最大旁瓣电平极小化,但旁瓣最大值的位置不固定,很难进行处理,所以采用积分旁瓣电平作为LS失配滤波器的目标函数,从滤波器开始运行直到当前时间,滤波器系数的更新是令总平方误差最小化,同时在主瓣方向输出为常数,于是有:

    $$ \left\{ \begin{array}{l} \min 20 \times \log 10\frac{{\sum\limits_i {\left| {{P_{s,i}}} \right|} }}{{\left| {{P_{\max }}} \right|}} \Rightarrow \left\{ \begin{array}{l} \min {{\boldsymbol{y}}^{\rm{H}}}{\boldsymbol{Wy}}\\ {{\boldsymbol{h}}^{\rm{H}}}{\boldsymbol{S}} = {\boldsymbol{N}} \end{array} \right.\\ {{\boldsymbol{h}}^{\rm{H}}}{\boldsymbol{S}} = {\boldsymbol{N}} \end{array} \right. $$ (11)
    $$ \begin{array}{l} {\boldsymbol{W}} = {\rm{diag}}(w( - P + 1), \cdots ,w(0),w(1), \cdots ,w(P - 1)) = \\ \;\;\;\;\;\;\;\;\;\;{\rm{diag}}(1, \cdots ,0,1, \cdots ,1) \end{array} $$ (12)

    其中,旁瓣加权矩阵W用于调整旁瓣电平和主瓣宽度,如果要增大主瓣宽度,可以减小w(1)、w(-1) 的值。该优化问题的解为:

    $$ {\boldsymbol{h}} = {({\boldsymbol{XW}}{{\boldsymbol{X}}^{\rm{H}}})^{ - 1}}{\boldsymbol{S}} $$ (13)

    从最小二乘失配滤波算法式 (13) 可以看到,LS失配滤波器作为一种自适应滤波器,需要进行矩阵的连乘、共轭转置和求逆等运算,运算量大、实现复杂,不适合工程应用[10-12]

    针对自适应滤波算法及其运算量带来的实际工程应用问题,对LS失配滤波方法进行改进,从式 (13) 可以看到,最小失配滤波器h与回波信号$[x(0),$$x(1), \cdots ,x(M - 1){]^{\rm{T}}}$有关,离线无法得到回波信号,则预先设置滤波器是不可能的,即需要对算法进行调整和改进。利用发射波形${[s(0),s(1), \cdots ,s(N - 1)]^{\rm{T}}}$代替式 (7) 中的回波信号序列${[x(0),x(1), \cdots ,x(M - 1)]^{\rm{T}}}$,构造新的系数矩阵A,从而实现最小二乘失配滤波器w的离线获取。

    调整后的P×(2P-1) 维的系数矩阵A为:

    $$ \begin{array}{l} \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\boldsymbol{A}} = \\ \left( {\begin{array}{*{20}{c}} 0& \cdots &0&{s(0)}&{s(1)}&{\begin{array}{*{20}{c}} \cdots &{s(N - 1)} \end{array}}\\ 0&{}&{}& \cdots &{\begin{array}{*{20}{c}} {}&{} \end{array}}&0\\ 0&{s(0)}&{s(1)}& \cdots &{s(N - 1)}& \cdots \\ {s(0)}&{s(1)}& \cdots &{s(N - 1)}&0&0 \end{array}} \right) \end{array} $$ (14)

    从而得到离线最小二乘失配滤波器w为:

    $$ {\boldsymbol{w}} = {({\boldsymbol{A}}\Lambda {{\boldsymbol{A}}^{\rm{H}}})^{ - 1}}{\boldsymbol{S}} $$ (15)

    式中,信号矢量S如式 (3) 所示。

    经过离线最小二乘滤波器处理结果为:

    $$ {\boldsymbol{y}} = {{\boldsymbol{w}}^{\rm{H}}}{\boldsymbol{X}} $$ (16)

    于是,脉压距离旁瓣抑制滤波器的离线实现流程如图 1所示,按照该流程实现将每个频率编码波形的滤波器实现计算好,将这些滤波器建立成库。在接收到雷达目标回波信号时,根据雷达发射波形信息等,从滤波器库中选取对应的旁瓣抑制滤波器,将选取的旁瓣抑制滤波器应用于回波信号采样数据,完成脉压和旁瓣抑制处理。

    图  1  脉压距离旁瓣抑制滤波器的离线实现流程图

    仿真中考虑的主要参数:码元宽度$\tau $,码长M,LFM带宽BLFM以及跳频间隔Δf。涉及到的LFM-频率编码复合波形包括LFM-Costas和LFM-步进频,其余参数如表 1表 2所示。

    表  1  LFM-Costas频率编码复合波形仿真实验参数
    参数 码元宽度
    /us
    脉冲宽度
    /us
    LFM带宽
    /MHz
    跳频间隔
    /MHz
    1 1 8 9 5
    2 20
    下载: 导出CSV 
    | 显示表格
    表  2  LFM-步进频率编码复合波形仿真实验参数
    参数 码元宽度
    /us
    脉冲宽度
    /us
    LFM带宽
    /MHz
    跳频间隔
    /MHz
    1 1 10 9 5
    2 20
    下载: 导出CSV 
    | 显示表格

    图 2图 3分别为LFM-Costas复合波形和LFM-步进频编码复合波形的失配滤波处理结果。

    图  2  LFM-Costas波形失配滤波结果
    图  3  LFM-步进频波形失配滤波结果

    图 2a所示,在表 1参数1条件下,相关旁瓣在-14.57 dB,失配处理的旁瓣仅为-28.74 dB,改善了14.17 dB。从图 2b可以看到,在表 1参数2条件下,相关旁瓣为-18.05 dB,失配滤波后为-28.26 dB,下降了10.21 dB。

    图 3a图 3b可以看到,表 2中2种参数条件下,相关旁瓣分别在-10.81 dB和-12.22 dB,而失配处理的旁瓣则分别为-35.83 dB和-31.04 dB,相比分别改善了25.02 dB和18.82 dB。

    图 2图 3的结果表明,失配滤波可以有效抑制复合波形的相关距离旁瓣。

    图 4图 5分别为50 kHz多普勒频移下,LFM-Costas复合波形和LFM-步进频编码复合波形的失配滤波处理结果。

    图  4  50 kHz多普勒频移下,LFM-Costas复合波形的处理结果
    图  5  50 kHz多普勒频移下,LFM-步进频编码复合波形的处理结果

    图 4所示,50 kHz的多普勒频移导致LFM-Costas失配处理后的主瓣下降了1.81 dB,旁瓣抬高了1.15 dB;50 KHz的多普勒频移引起LFM-步进频失配处理后的主瓣和旁瓣分别恶化了0.64 dB和2.25 dB,如图 5所示。

    图 4图 5的结果显示,多普勒频移会带来失配滤波性能的下降。

    不同多普勒频移下,失配滤波对复合波形处理后主瓣和旁瓣的变化情况如图 6所示,可以看到,随着多普勒频移的增大,主瓣和旁瓣性能的恶化加剧,在多普勒频移为80 kHz时,复合波形失配滤波处理主瓣和旁瓣性能恶化5 dB左右,峰值旁瓣仍优于-22 dB,在可接受范围之内。

    图  6  不同多普勒频率下,失配滤波处理结果

    针对LFM-频率编码复合低截获波形的高距离旁瓣问题,提出一种快速旁瓣抑制实现方法,通过对最小二乘失配滤波方法的改进,实现旁瓣抑制滤波器的离线获取,降低滤波器在线应用的运算量。该方法离线预先获取旁瓣抑制失配滤波器、建立旁瓣抑制滤波器库,依据雷达/波形参数从旁瓣抑制滤波器库选取处理滤波器,在线实时处理目标回波信号。仿真实验结果表明,本文提出的快速旁瓣抑制方法简单、可行,为LFM-频率编码复合低截获波形的距离旁瓣抑制提供了一种有效的解决途径。

  • 图  1   脉压距离旁瓣抑制滤波器的离线实现流程图

    图  2   LFM-Costas波形失配滤波结果

    图  3   LFM-步进频波形失配滤波结果

    图  4   50 kHz多普勒频移下,LFM-Costas复合波形的处理结果

    图  5   50 kHz多普勒频移下,LFM-步进频编码复合波形的处理结果

    图  6   不同多普勒频率下,失配滤波处理结果

    表  1   LFM-Costas频率编码复合波形仿真实验参数

    参数 码元宽度
    /us
    脉冲宽度
    /us
    LFM带宽
    /MHz
    跳频间隔
    /MHz
    1 1 8 9 5
    2 20
    下载: 导出CSV

    表  2   LFM-步进频率编码复合波形仿真实验参数

    参数 码元宽度
    /us
    脉冲宽度
    /us
    LFM带宽
    /MHz
    跳频间隔
    /MHz
    1 1 10 9 5
    2 20
    下载: 导出CSV
  • [1]

    LYNCH D J. Introduction to RF stealth[M]. New Jersey:SciTech Publishing Inc, 2004:275-350.

    [2]

    PHILLIP E P. Detective and classcifying low probability of intercept radar[M]. 2nd ed. Boston:Artech House, 2009:195-198.

    [3]

    DONOHOE J P, INGELS F M. The ambiguity properties of FSK/PSK signals[C]//Record of the 1990 IEEE International Radar Conference. Arlington, USA:IEEE Press, 1990.

    [4]

    LEMIEUX J A. Analysis of an optimum hybrid radar waveform using frequency hopping and locally optimum signals[C]//Proceedings of the IEEE National Radar Conference. Califomia, USA:IEEE Press, 1991.

    [5]

    SKINNER B J, DONOHOE J P, INGELS F M. Simplified performance estimation of FSK/PSK hybrid signaling radar systems[C]//Proceedings of the 1993 IEEE National aerospace and electronics conference. Dayton, USA:IEEE Press, 1993.

    [6]

    SKINNER B J, DONOHOE J P, INGELS F M. Matched FSK/PSK radar[C]//Record of the 1994 IEEE National Radar Conference. Atlanta, USA:IEEE Press, 1994.

    [7] 徐健. 多相编码脉冲压缩雷达信号的优化设计与仿真[D]. 南京: 河海大学, 2007.

    XU Jian. Polyphase pulse compression codes radar signal optimizing and simulation[D]. Nanjing:Hohai University, 2007.

    [8] 俞璐璐. 多相编码信号研究[D]. 南京: 南京理工大学, 2012.

    YU Lu-lu. Research on polyphase signal[D]. Nanjing:Nanjing University of Science and Technology, 2012.

    [9]

    BERNARD W, SAMUEL D S. Adaptive signal processing[M]. New Jersey, USA:Prentice Hall, 1985:289-327.

    [10]

    SIMON H. Adaptive filter theory[M]. 5th ed. New York, USA:Pearson Education Publishing Inc, 2013:423-456.

    [11]

    ACKROYD M H, GHANI F. Optimum mismatched filters for sidelobe suppression[J]. IEEE Transaction on Aerospace and Electronic Systems, 1973, 3:214-218.

    [12]

    BADEN J M, COHEN M N. Optimal Sidelobe Suppression for Biphase Codes[C]//Proceedings of the 1991 National Telesystems Conference. Altanta GA:IEEE Press, 1991.

    [13]

    ZRNIC B, ZEJAK A, PETROVIC A, et al. Range sidelobe suppression for range sidelobe suppression for pulse compression radars utilizing modified RLS algorithm[C]//Proceedings of the 1998 IEEE 5th International Symposium on Spread Spectrum Techniques and Applications. South Africa:IEEE Press, 1998.

    [14]

    SATO R, SHINRHU M. Simple mismatched filter for binary pulse compression code with small PSL and small S/N loss[J]. IEEE Transactions on Aerospace and Electronic systems, 2003, 39(2):711-718. DOI: 10.1109/TAES.2003.1207277

  • 期刊类型引用(7)

    1. 肖跃,高振斌,赵子杭. 一种复合调制雷达信号设计及性能分析. 激光与光电子学进展. 2022(11): 245-250 . 百度学术
    2. 刘子豪,任丽香,孙颖豪,范花玉. 一种捷变脉冲串的动态规划优化设计方法. 空军预警学院学报. 2021(03): 157-161 . 百度学术
    3. 葛萌萌,余显祥,严正欣,方学立,崔国龙,孔令讲. 脉间波形幅相联合设计抗欺骗干扰方法. 电子科技大学学报. 2021(04): 481-487 . 本站查看
    4. 郝志梅,孙进平,罗美方. 基于复合频率编码的低截获概率波形簇设计. 系统工程与电子技术. 2021(11): 3137-3143 . 百度学术
    5. 付银娟,李勇,徐丽琴,张昆辉. NLFM-Costas射频隐身雷达信号设计及分析. 吉林大学学报(工学版). 2019(03): 994-999 . 百度学术
    6. 付银娟,李勇,卢光跃,张昆辉. 脉间复合调频脉内相位编码雷达信号设计与分析. 哈尔滨工程大学学报. 2019(07): 1347-1353 . 百度学术
    7. 王亚涛,王昊飞. 复合调制雷达波形低截获性能分析. 兵器装备工程学报. 2019(12): 120-124 . 百度学术

    其他类型引用(4)

图(6)  /  表(2)
计量
  • 文章访问数:  6596
  • HTML全文浏览量:  2545
  • PDF下载量:  178
  • 被引次数: 11
出版历程
  • 收稿日期:  2015-10-30
  • 修回日期:  2016-09-25
  • 刊出日期:  2017-03-14

目录

/

返回文章
返回