CN103728594B - 基于多通道nlms的外辐射源雷达海杂波干扰抑制方法 - Google Patents
基于多通道nlms的外辐射源雷达海杂波干扰抑制方法 Download PDFInfo
- Publication number
- CN103728594B CN103728594B CN201310752143.3A CN201310752143A CN103728594B CN 103728594 B CN103728594 B CN 103728594B CN 201310752143 A CN201310752143 A CN 201310752143A CN 103728594 B CN103728594 B CN 103728594B
- Authority
- CN
- China
- Prior art keywords
- hyperchannel
- nlms
- wave filter
- moment
- passage
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/2813—Means providing a modification of the radiation pattern for cancelling noise, clutter or interfering signals, e.g. side lobe suppression, side lobe blanking, null-steering arrays
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于多通道NLMS的外辐射源雷达海杂波干扰抑制方法,属于雷达信号处理技术领域。本发明将直达波天线接收到的参考信号调频到海杂波的多普勒频率上,构成与海杂波多普勒频率分量相对应的多个参考通道,得到一种多通道归一化最小均方(NLMS)自适应对消滤波器,建立多通道NLMS滤波器的数学模型,构成多个参考信号对消不同频率的海杂波干扰,可以较好的抑制展宽的海杂波干扰,提高对消增益和目标信噪比。
Description
技术领域
本发明涉及一种基于多通道NLMS的外辐射源雷达海杂波干扰抑制方法,属于雷达信号处理技术领域。
背景技术
外辐射源雷达通常有两个接收天线,一个是直达波天线,指向辐射源接收其发射的参考信号;另一个是回波天线,指向搜索区域接收运动目标反射的回波信号。当外辐射源雷达探测海面目标时,其配置如附图1所示,辐射源发射的信号除了被目标反射外,还可经其它路径进入回波天线,如经海面反射进入回波天线主瓣的海杂波干扰、直接从回波天线副瓣进入的直达波干扰等。因此,回波天线接收到的信号中有较强的直达波和海杂波干扰,而目标反射的回波信号与之相比十分微弱,如果不对这些干扰进行有效抑制,就无法检测目标。
由于海面复杂的运动导致回波天线接收到的海杂波干扰是动态变化的,其在多普勒谱上是展宽的,如果使用传统的适用于静态杂波干扰抑制的归一化最小均方(NLMS)自适应对消算法抑制海杂波,将无法有效抑制展宽的海杂波干扰,附图2为实际数据采用NLMS算法对消后互模糊函数图,从图中可以看出,NLMS对消后残余海杂波较多,会造成比较大的虚警概率。采用较大的步长一定程度上可以抑制展宽的海杂波干扰,但是较大的步长会导致目标信噪比有较大的损失,影响目标检测性能。
归一化最小均方(NLMS)滤波器滤波过程可表示为以下步骤:
1)计算估计误差(对消结果),设滤波器阶数为M
e(n)=d(n)-wH(n)u(n)(1)式中d(n)为n时刻回波信号值,wH(n)为n时刻滤波器权向量的共轭转置(1×M维),u(n)为直达波向量(M×1维),wH(n)u(n)为n时刻滤波器对期望值的估计,e(n)为n时刻滤波器的估计误差(对消输出)。
2)更新滤波器系数
式中为更新步长,λ'为一个小的正常数(通常取0.001),作用是防止除法溢出,e*(n)为对消输出e(n)的共轭,w(n)为n时刻NLMS滤波器的权系数向量(M×1维),w(n+1)为n+1时刻NLMS滤波器的权系数向量。
发明内容
本发明的目的是针对外辐射源雷达探测海面目标时,传统NLMS算法不能有效抑制海杂波干扰的问题,提出一种基于多通道NLMS的外辐射源雷达海杂波干扰抑制方法,能较好地抑制海杂波干扰,提高对消增益和目标信噪比。
本发明具体包括如下步骤:
步骤1,将直达波天线接收到的参考信号调频到海杂波的多普勒频率上,构成与海杂波多普勒频率分量相对应的多个参考通道,得到一种多通道归一化最小均方(NLMS)自适应对消滤波器,建立多通道NLMS滤波器的数学模型如下:
式中T表示矩阵转置,H表示共轭转置,*表示复共轭,W(n)为多通道NLMS滤波器系数矩阵,K为通道个数,wi(n)为n时刻多通道NLMS滤波器第i(i=1,2,...,K)个通道的抽头权向量(M×1维)。U(n)为多通道NLMS滤波器输入矩阵,fs为采样频率,fi(f1=0)为调频率, 为n时刻多通道NLMS滤波器第i个通道的抽头输入向量(M×1维,M为滤波器阶数)。u1(n)为原始参考信号,d(n)为n时刻回波值,e(n)为n时刻多通道NLMS滤波器的输出值也即对消结果,wi(n+1)为多通道NLMS滤波器第i个通道系数更新方程,为n时刻K个参考通道输入值的范数平方和,为多通道NLMS滤波器第i个通道的更新步长,λ'是防止除法溢出取的较小值(通常取0.001)。
步骤2,利用步骤1建立的多通道NLMS滤波器模型进行多通道NLMS对消。
步骤2.1,参数初始化
根据杂波的长度,确定对消滤波器阶数M,然后分析确定海杂波多普勒频谱展宽范围包含的多普勒单元个数,设定展宽的单元个数作为参考通道个数K,并根据积累时间确定对应单元的频率fi(海杂波主多普勒频率,也是参考通道调频率),然后将直达波天线接收到的参考信号进行调频,构成K个参考通道,其中第一个通道为原始参考信号,其余参考通道信号为原始参考信号的调频信号,U(n)=[u1(n);u2(n);...;uK(n)]=0为K个参考通道初始输入矩阵,W(n)=[w1(n);w2(n);...;wK(n)]=0为K个通道自适应滤波器初始权系数矩阵。
步骤2.2,更新输入矩阵U(n)
将n时刻K个参考通道新的输入值其中u1(n)为n时刻原始参考信号值,fs为信号采样频率,fi(f1=0)为多通道NLMS滤波器第i(i=1,2,...,K)个通道的调频率,加入到输入向量ui(n)(i=1,2,...,K)中,构成新的输入矩阵
式中
为n时刻多通道NLMS滤波器第i个通道的抽头输入向量(M×1维)。
步骤2.3,计算n时刻的多通道NLMS滤波器估计误差,也即对消结果
式中为n时刻多通道滤波器第i个通道权系数向量的共轭转置,为n时刻K个通道估计值之和,d(n)为滤波器在n时刻的期望值也即回波信号在n时刻的值,e(n)为估计误差也即对消结果;
步骤2.4,更新多通道NLMS滤波器的权系数矩阵W(n+1)
首先对每个通道的权向量进行更新,更新方程如下,
进而得到多通道自适应滤波器权系数矩阵,
W(n+1)=[w1(n+1),w2(n+1),...,wK(n+1)](10)
式中ui(n)为n时刻多通道NLMS滤波器第i个通道的输入向量,是多通道NLMS滤波器第i个通道的更新步长,λ'为防止除法溢出取的较小值(通常取0.001),wi(n)为n时刻多通道NLMS滤波器第i个通道的权系数向量;为n时刻K个通道输入值范数的平方和。wi(n+1)为n+1时刻多通道NLMS滤波器第i个通道更新的权系数向量。得到每个通道的权向量后,便构成更新的权系数矩阵W(n+1)。
步骤2.5,迭代次数加1,循环执行步骤2.2~步骤2.4,当迭代次数等于需要滤波的数据长度时,循环结束。从而得到利用多通道NLMS方法对消后的数据。
有益效果
本发明通过将参考信号调频到对应的海杂波主导分量的多普勒频率上,构成多个参考信号对消不同频率的海杂波干扰,可以较好的抑制展宽的海杂波干扰,提高对消增益和目标信噪比。
附图说明
图1是背景技术中外辐射源雷达海面目标探测***配置示意图;
图2是背景技术中实际数据采用NLMS算法抑制海杂波后的互模糊函数图;
图3是本发明的多通道NLMS方法结构图;
图4是具体实施方式中实际数据采用多通道NLMS算法抑制海杂波后的互模糊函数图。
具体实施方式
下面结合实施例对本发明做进一步说明。
实施例
实际采集数据验证使用的数据为青岛海边采集数据,取其中一组数据验证算法的有效性。取实际采集数据中直达波和回波各5000000点,基带采样频率为fs=10MHz,通过分析确定海杂波主导分量的多普勒频率主要集中在±2Hz和±4Hz上,然后调频原始参考信号,得到与海杂波主导分量多普勒频率对应的4个参考通道,加上原始参考通道,构成5个参考通道。
本发明的多通道NLMS算法用于本实施例中,包括以下步骤:
1)滤波器参数初始化,通道个数K=5,滤波器阶数M=100,则
2)更新抽头输入矩阵U(n),将n时刻5个参考通道输入值加入到抽头输入向量 中,式中f1=0Hz,f2=2Hz,f3=-2Hz,f4=4Hz,f5=-4Hz,构成n时刻新的抽头输入矩阵:
3)计算n时刻的多通道NLMS滤波器输出值e(n)
式中d(n)为n时刻的回波通道的值。
4)更新多通道滤波器的抽头权系数矩阵W(n+1)
式中λ'=0.001,得到新的权系数矩阵W(n+1)
W(n+1)=[w1(n+1),w2(n+1),...,w5(n+1)](15)
5)取n=n+1,循环执行步骤2)~步骤4),经5000000(数据长度)次迭代后得到全部数据的对消结果。
图2和图4分别为实际采集数据采用NLMS算法和采用本发明提出的多通道NLMS算法对消后的互模糊函数图。其中采用NLMS算法的对消增益为22.0963dB,目标信噪比为14.8445dB,图中残余海杂波较多。采用本发明提出的多通道NLMS算法的对消增益为24.2153dB,目标信噪比为15.1492dB,从图中可以看出,展宽的海杂波得到有效的抑制。因此,采用本发明提出的多通道NLMS算法可以有效抑制展宽的海杂波干扰,提高对消增益和目标信噪比。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施实例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.基于多通道NLMS的外辐射源雷达海杂波干扰抑制方法,其特征在于:具体包括如下步骤:
步骤1,将直达波天线接收到的参考信号调频到海杂波的多普勒频率上,构成与海杂波多普勒频率分量相对应的多个参考通道,得到一种多通道NLMS自适应对消滤波器,建立多通道NLMS滤波器的数学模型如下:
式中T表示矩阵转置,H表示共轭转置,*表示复共轭,W(n)为多通道NLMS滤波器系数矩阵,K为通道个数,wi(n)为n时刻多通道NLMS滤波器第i,i=1,2,...,K,个通道的M×1维抽头权向量;M为滤波器阶数;U(n)为多通道NLMS滤波器输入矩阵,fs为采样频率,fi且f1=0为调频率,
步骤2,利用步骤1建立的多通道NLMS滤波器模型进行多通道NLMS对消;
步骤2.1,参数初始化
根据杂波的长度,确定对消滤波器阶数M,然后分析确定海杂波多普勒频谱展宽范围包含的多普勒单元个数,设定展宽的单元个数作为参考通道个数K,并根据积累时间确定对应单元的海杂波主多普勒频率fi,然后将直达波天线接收到的参考信号进行调频,构成K个参考通道,其中第一个通道为原始参考信号,其余参考通道信号为原始参考信号的调频信号,U(n)=[u1(n);u2(n);...;uK(n)]=0为K个参考通道初始输入矩阵,W(n)=[w1(n);w2(n);...;wK(n)]=0为K个通道自适应滤波器初始权系数矩阵;
步骤2.2,更新输入矩阵U(n)
将n时刻K个参考通道新的输入值其中u1(n)为n时刻原始参考信号值,fs为信号采样频率,fi为多通道NLMS滤波器第i个通道的调频率,加入到输入向量ui(n)中,构成新的输入矩阵
步骤2.3,计算n时刻的多通道NLMS滤波器估计误差,也即对消结果
式中为n时刻多通道滤波器第i个通道权系数向量的共轭转置,为n时刻K个通道估计值之和,d(n)为滤波器在n时刻的期望值即回波信号在n时刻的值,e(n)为估计误差即对消结果;
步骤2.4,更新多通道NLMS滤波器的权系数矩阵W(n+1)
首先对每个通道的权向量进行更新,更新方程如下,
进而得到多通道自适应滤波器权系数矩阵,
W(n+1)=[w1(n+1),w2(n+1),...,wK(n+1)](10)
式中ui(n)为n时刻多通道NLMS滤波器第i个通道的输入向量,是多通道NLMS滤波器第i个通道的更新步长,λ'为防止除法溢出的取值,wi(n)为n时刻多通道NLMS滤波器第i个通道的权系数向量;为n时刻K个通道输入值范数的平方和;wi(n+1)为n+1时刻多通道NLMS滤波器第i个通道更新的权系数向量;得到每个通道的权向量后,便构成更新的权系数矩阵W(n+1);
步骤2.5,迭代次数加1,循环执行步骤2.2~步骤2.4,当迭代次数等于需要滤波的数据长度时,循环结束;从而得到利用多通道NLMS方法对消后的数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310752143.3A CN103728594B (zh) | 2013-12-31 | 2013-12-31 | 基于多通道nlms的外辐射源雷达海杂波干扰抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310752143.3A CN103728594B (zh) | 2013-12-31 | 2013-12-31 | 基于多通道nlms的外辐射源雷达海杂波干扰抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103728594A CN103728594A (zh) | 2014-04-16 |
CN103728594B true CN103728594B (zh) | 2016-03-30 |
Family
ID=50452749
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310752143.3A Expired - Fee Related CN103728594B (zh) | 2013-12-31 | 2013-12-31 | 基于多通道nlms的外辐射源雷达海杂波干扰抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103728594B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926570B (zh) * | 2014-04-22 | 2016-05-04 | 北京理工大学 | 基于多通道rls的外辐射源雷达海杂波干扰抑制方法 |
CN105182300B (zh) * | 2015-09-11 | 2018-06-19 | 上海无线电设备研究所 | 天基外辐射源雷达的地杂波干扰实时抑制方法 |
CN107656246A (zh) * | 2017-08-18 | 2018-02-02 | 上海无线电设备研究所 | 一种天基海空动目标广域探测*** |
CN109061597B (zh) * | 2018-08-23 | 2022-04-15 | 哈尔滨工业大学 | 基于盲源分离与时频脊波域滤波的电离层杂波抑制方法 |
CN110146848B (zh) * | 2019-05-22 | 2023-06-23 | 西安电子科技大学 | 基于分数阶最小均方的调频连续波雷达自干扰消除方法 |
CN110333485B (zh) * | 2019-06-21 | 2021-10-01 | 大亚湾核电运营管理有限责任公司 | 一种海杂波抑制方法、装置及终端设备 |
CN110646769B (zh) * | 2019-09-03 | 2021-07-20 | 武汉大学深圳研究院 | 一种适用于lte外辐射源雷达的时域杂波抑制方法 |
CN110967675B (zh) * | 2019-10-07 | 2022-04-22 | 西安电子科技大学 | 一种基于神经网络的卫星外辐射源目标被动探测方法 |
CN111337895B (zh) * | 2020-01-14 | 2021-11-09 | 北京理工大学 | 多通道海杂波空时相关性分析方法 |
CN115061099A (zh) * | 2022-07-28 | 2022-09-16 | 南京华成微波技术有限公司 | 一种雷达非平稳随队干扰抵消方法、装置和终端设备 |
CN116953654A (zh) * | 2023-09-21 | 2023-10-27 | 北京中科睿信科技有限公司 | 一种机载外辐射源雷达杂波对消方法及设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3639500A1 (de) * | 1986-11-20 | 1988-06-01 | Forschungsgesellschaft Fuer An | Radarempfaenger fuer bewegte radargeraete mit gruppenantenne mit zweidimensional wirkender clutterunterdrueckung |
CN101806887A (zh) * | 2010-03-19 | 2010-08-18 | 清华大学 | 基于空间跟踪滤波器的海杂波抑制和目标检测方法 |
CN101907703A (zh) * | 2010-07-20 | 2010-12-08 | 中国人民解放军空军雷达学院 | 天波超视距雷达空时联合自适应抗干扰方法 |
US8456350B2 (en) * | 2010-12-17 | 2013-06-04 | Raytheon Company | Methods and apparatus for sea state measurement via radar sea clutter eccentricity |
-
2013
- 2013-12-31 CN CN201310752143.3A patent/CN103728594B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE3639500A1 (de) * | 1986-11-20 | 1988-06-01 | Forschungsgesellschaft Fuer An | Radarempfaenger fuer bewegte radargeraete mit gruppenantenne mit zweidimensional wirkender clutterunterdrueckung |
CN101806887A (zh) * | 2010-03-19 | 2010-08-18 | 清华大学 | 基于空间跟踪滤波器的海杂波抑制和目标检测方法 |
CN101907703A (zh) * | 2010-07-20 | 2010-12-08 | 中国人民解放军空军雷达学院 | 天波超视距雷达空时联合自适应抗干扰方法 |
US8456350B2 (en) * | 2010-12-17 | 2013-06-04 | Raytheon Company | Methods and apparatus for sea state measurement via radar sea clutter eccentricity |
Also Published As
Publication number | Publication date |
---|---|
CN103728594A (zh) | 2014-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103728594B (zh) | 基于多通道nlms的外辐射源雷达海杂波干扰抑制方法 | |
Kay | Waveform design for multistatic radar detection | |
CN105527610B (zh) | 基于分数阶时延估计的多天线联合优化杂波抑制方法 | |
CN102830387B (zh) | 一种基于数据预处理的协方差矩阵正交化波束形成方法 | |
CN107153178B (zh) | 外辐射源雷达参考信号含有多径干扰时的目标检测方法 | |
CN104865586A (zh) | 一种阵列天线导航接收机抗干扰及多径抑制方法 | |
CN104360355B (zh) | 抗干扰方法和装置 | |
CN107462877B (zh) | 一种基于先验知识的天波雷达海杂波抑制方法 | |
CN106842140B (zh) | 一种基于和差波束降维的主瓣干扰抑制方法 | |
CN103984676A (zh) | 一种基于协方差矩阵重构的正交投影自适应波束形成方法 | |
CN103197300B (zh) | 一种基于gpu的外辐射源雷达直达波杂波对消实时处理方法 | |
CN104007453A (zh) | 概率搜索辅助的频域空域联合抗干扰方法 | |
CN103760529B (zh) | 一种无源探测高效级联空时自适应处理方法 | |
CN105302936A (zh) | 基于相关计算和协方差矩阵重构的自适应波束形成方法 | |
CN102736069B (zh) | 一种直达波干扰抑制方法 | |
CN103885045B (zh) | 基于子阵划分的循环联合自适应波束形成方法 | |
CN103344945B (zh) | 一种直达波及其多径干扰的抑制方法 | |
CN110865345B (zh) | 一种快速自适应脉冲压缩方法 | |
CN101907702A (zh) | 用于mimo雷达的两维多脉冲对消器 | |
CN104777491A (zh) | 一种盲波束宽带干扰抑制方法及装置 | |
Meller et al. | Processing of noise radar waveforms using block least mean squares algorithm | |
CN103926570B (zh) | 基于多通道rls的外辐射源雷达海杂波干扰抑制方法 | |
Tian et al. | Underwater multi-target passive detection based on transient signals using adaptive empirical mode decomposition | |
Kelly et al. | RFI suppression and sparse image formation for UWB SAR | |
CN105572642B (zh) | 一种基于两级架构的空时自适应处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160330 Termination date: 20161231 |