CN113541851B - 一种稳态宽带电磁频谱抑制方法 - Google Patents
一种稳态宽带电磁频谱抑制方法 Download PDFInfo
- Publication number
- CN113541851B CN113541851B CN202110821164.0A CN202110821164A CN113541851B CN 113541851 B CN113541851 B CN 113541851B CN 202110821164 A CN202110821164 A CN 202110821164A CN 113541851 B CN113541851 B CN 113541851B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- background signal
- value
- frequency spectrum
- background
- 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04K—SECRET COMMUNICATION; JAMMING OF COMMUNICATION
- H04K3/00—Jamming of communication; Counter-measures
- H04K3/20—Countermeasures against jamming
- H04K3/22—Countermeasures against jamming including jamming detection and monitoring
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04K—SECRET COMMUNICATION; JAMMING OF COMMUNICATION
- H04K3/00—Jamming of communication; Counter-measures
- H04K3/20—Countermeasures against jamming
- H04K3/22—Countermeasures against jamming including jamming detection and monitoring
- H04K3/224—Countermeasures against jamming including jamming detection and monitoring with countermeasures at transmission and/or reception of the jammed signal, e.g. stopping operation of transmitter or receiver, nulling or enhancing transmitted power in direction of or at frequency of jammer
- H04K3/228—Elimination in the received signal of jamming or of data corrupted by jamming
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种稳态宽带电磁频谱抑制方法,包括通过宽带接收机采集宽带背景信号,对采集得到的宽带背景信号进行N阶离散傅立叶变换;将当前每点的值与前次对应位置的值进行比较后,保留两者中的最大值;计算本次保留的最大值与上一次比较保留的最大值的差值,当差值小于设定的阈值时,则停止背景信号采集,判定背景信号频谱稳定并得到背景信号频谱;根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制。通过本发明,可以实现消除背景信号频谱,将新出现的信号凸显出来的作用。
Description
技术领域
本发明涉及通信领域,具体是一种稳态宽带电磁频谱抑制方法。
背景技术
在现有技术中,信号检测技术是基于频谱进行信号提取的,考虑到检测的信号背景环境异常复杂,在宽带频谱内信号个数和信号样式非常繁杂,且容易受到有意或无意信号干扰。目前大多采用的设计实现思路:一种是基于空域自适应滤波的方法,抑制某些方向的来波干扰信号,这样的技术途径往往需要相控阵天线,从而导致系统部署不方便,且当来波方向一致时该方法失效;第二种是用波形串行干扰抵消的方法,该方法能够解决来波方向一致时的干扰抑制问题,通过耦合器获取本端发射信号,使耦合信号通过由衰减器、放大器、移相器和时延线等组成的模拟干扰信道,将相位调整为与干扰信号相差180度的奇数倍后,再通过耦合器输出至接收端,从而实现干扰信号的抑制。但当频谱复杂,干扰背景信号较多时,就无法实现多信号干扰的抑制。考虑到复杂电磁背景下的信号监测时,主要利用的频谱信息,如果能将宽带频谱进行有效抑制,将会极大地提升侦察设备的检测正确性,因此,一种稳态宽带电磁频谱抑制方法亟待研究,以确保对抗条件下的监测设备检测能力。
发明内容
本发明的目的在于克服现有技术的不足,提供一种稳态宽带电磁频谱抑制方法,包括如下步骤:
步骤一,通过宽带接收机采集宽带背景信号,对采集得到的宽带背景信号进行N阶离散傅立叶变换;
步骤二,将当前N阶离散傅立叶变换得到的每点的模值与前次对应的点的模值进行比较后,保留两者中的最大值;
步骤三,计算本次保留的最大值与上一次比较保留的最大值的差值,当差值小于设定的阈值时,则停止背景信号采集,判定背景信号频谱稳定并得到背景信号频谱,执行步骤四,若大于阈值,则执行步骤一;
步骤四,根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制,当接收到的频谱数据的频谱幅度超过背景信号频谱幅度时,则采用当前的频谱值;当接收到的频谱数据的频谱幅度低于背景信号频谱幅度时,采用底噪进行填充,得到抑制背景噪声后的频谱数据。
进一步的,所述的对采集得到的宽带背景信号进行N阶离散傅立叶变换,包括如下过程:
采集宽带背景信号,通过下式计算离散频谱,
其中,x(n),n=0...N-1为采集的宽带背景信号,X(m)为其离散傅立叶变换值。
进一步的,所述的将当前每点的值与前次对应位置的值进行比较后,保留两者中的最大值,采用如下公式:
定义变量βi(m)=Max[βi-1(m),|Xi(m)|],m=0...N-1,i=2,3,......,
其中|.|表示对频谱取模运算,且β1(m)=|X1(m)|。
进一步的,所述的计算本次比较后保留的最大值与上一次比较后保留的最大值的差值,采用如下公式:
ρi(m)=|βi(m)-βi-1(m)|,i=2,3,.......,m=0...N-1,
当ρi(m)<γ时,停止停止背景信号采集,γ为设定的阈值。
进一步的,所述的根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制,当接收到的频谱数据的频谱幅度超过背景信号频谱幅度时,则采用当前的频谱值;当接收到的频谱数据的频谱幅度低于背景信号频谱幅度时,采用底噪进行填充,得到抑制背景噪声后的频谱数据,包括如下过程:
当频谱幅度超过背景信号频谱时,采用当前的频谱值;当频谱幅度低于背景信号频谱时,采用底噪进行填充;即:
其中,X(m)为当前采集到的背景信号的频谱数据,Xn(m)为拟合出来的频谱底噪,βk(m)为前面得到的背景信号频谱;
本发明的有益效果是:本发明可以较好地消除背景信号频谱,将新出现的信号凸显出来,为后续载波检测很好地消除了背景信号的干扰。
附图说明
图1为一种稳态宽带电磁频谱抑制方法的流程示意图;
图2为无目标信号时的背景频谱;
图3为无目标信号时的消除频谱;
图4为包含有目标信号的背景频谱;
图5为包含有目标信号的消除频谱。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
如图1所示,一种稳态宽带电磁频谱抑制方法,包括如下步骤:
步骤一,通过宽带接收机采集宽带背景信号,对采集得到的宽带背景信号进行N阶离散傅立叶变换;
步骤二,将当前N阶离散傅立叶变换得到的每点的模值与前次对应的点的模值进行比较后,保留两者中的最大值;
步骤三,计算本次保留的最大值与上一次比较保留的最大值的差值,当差值小于设定的阈值时,则停止背景信号采集,判定背景信号频谱稳定并得到背景信号频谱,执行步骤四,若大于阈值,则执行步骤一;
步骤四,根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制,当接收到的频谱数据的频谱幅度超过背景信号频谱幅度时,则采用当前的频谱值;当接收到的频谱数据的频谱幅度低于背景信号频谱幅度时,采用底噪进行填充,得到抑制背景噪声后的频谱数据。
所述的对采集得到的宽带背景信号进行N阶离散傅立叶变换,包括如下过程:
采集宽带背景信号,通过下式计算离散频谱,
其中,x(n),n=0...N-1为采集的宽带背景信号,X(m)为其离散傅立叶变换值。
所述的将当前每点(1~N)的值与前次对应位置的值进行比较后,保留两者中的最大值,采用如下公式:
定义变量βi(m)=Max[βi-1(m),|Xi(m)|],m=0...N-1,i=2,3,......,
其中}.|表示对频谱取模运算,且β1(m)=|X1(m)|。
所述的计算本次比较后保留的最大值与上一次比较后保留的最大值的差值,采用如下公式:
ρi(m)=|βi(m)-βi-1(m)|,i=2,3,.......,m=0...N-1,
当ρi(m)<γ时,停止停止背景信号采集,γ为设定的阈值。
所述的根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制,当接收到的频谱数据的频谱幅度超过背景信号频谱幅度时,则采用当前的频谱值;当接收到的频谱数据的频谱幅度低于背景信号频谱幅度时,采用底噪进行填充,得到抑制背景噪声后的频谱数据,包括如下过程:
当频谱幅度超过背景信号频谱时,采用当前的频谱值;当频谱幅度低于背景信号频谱时,采用底噪进行填充;即:
其中,X(m)为当前采集到的背景信号的频谱数据,Xn(m)为拟合出来的频谱底噪,βk(m)为前面得到的背景信号频谱;
具体的,通过宽带接收机采集宽带数据,利用采集得到宽带数据进行N阶离散傅立叶变换,将当前每点(1~N)与前次对应位置值相比较,比较后保留两者的最大值;计算本次最大值与上次最大值的差值,当该值小于一定特定阀值时,停止背景频谱采集,即认为背景频谱收集充分,充分收集完背景频谱数据后,就需要进行背景频谱抑制,抑制后的频谱就能够突现目标信号。
其中,x(n),n=0...N-1为采集的宽带背景信号,X(m)为其离散傅立叶变换值,也中文也称之为背景频谱。由于收集背景频谱需要多次采集,故定义Xi(m),i=1,2,......为第i次采集数据计算得到的背景频谱,。
定义变量βi(m)=Max[βi-1(m),|Xi(m)|],m=0...N-1,i=2,3,......,其中|.|表示对频谱取模运算,且β1(m)=|X1(m)|。
计算ρi(m)=|βi(m)-βi-1(m)|,i=2,3,.......,m=0...N-1,该变量用于衡量背景信号变化程度的指征,当对于i=k,ρk(m)<γ,γ为设定的门限为一常量,即可认为背景信号的频谱特征已经进行了充分地采集,此时背景频谱采集流程完成,进入背景频谱抑制环节。
背景频谱抑制是利用上述处理得到的背景信号频谱βk(m),对当前接收到的频谱数据进行频谱抑制处理。常规的消除方法一般采用直接相减的处理方式,这种处理方式存在一个显而易见的弊端:由于背景信号是时变的,当进行后续检测是,如果背景信号消失,抑制后的频谱会出现较大的凹陷,对后续抑制后频谱显示和载波检测都会带来干扰。为了解决这个问题,采用了底噪补偿的算法,当频谱幅度超过背景信号频谱时,采用当前的频谱值;当频谱幅度低于背景信号频谱时,采用底噪进行填充;即:
其中,X(m)为当前采集到信号的频谱数据,Xn(m)为拟合出来的频谱底噪,βk(m)为前面得到的背景信号频谱。频谱底噪Xn(m)的估计处理方式:取背景频谱的最小值作为噪底基准,再估算背景频谱的抖动方差,然后根据抖动方差生成随机噪声信号,最后叠加到噪底基准上形成一个随机底噪。具体表现形式为Xn(m)=α(m),m=0...N-1,每个α(m)符合N(μ,σ2)的正态分布,其中μ=Min(βk(m)),m=0...N-1,
通过上述一系列处理得到的抑制背景噪声后的频谱数据SPd(m),基于该频谱数据就非常有利于后续的信号检测和目标发现。
以上所述仅是本发明的优选实施方式,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (5)
1.一种稳态宽带电磁频谱抑制方法,其特征在于,包括如下步骤:
步骤一,通过宽带接收机采集宽带背景信号,对采集得到的宽带背景信号进行N阶离散傅立叶变换;
步骤二,将当前N阶离散傅立叶变换得到的每点的模值与前次对应的点的模值进行比较后,保留两者中的最大值;
步骤三,计算本次保留的最大值与上一次比较保留的最大值的差值,当差值小于设定的阈值时,则停止背景信号采集,判定背景信号频谱稳定并得到背景信号频谱,执行步骤四,若大于阈值,则执行步骤一;
步骤四,根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制,当接收到的频谱数据的频谱幅度超过背景信号频谱幅度时,则采用当前的频谱幅度值;当接收到的频谱数据的频谱幅度低于背景信号频谱幅度时,采用底噪进行填充,得到抑制背景噪声后的频谱数据。
3.根据权利要求2所述的一种稳态宽带电磁频谱抑制方法,其特征在于,将当前每点的值与上一次比较后保留的值进行比较后的最大值,采用如下公式:
定义变量βi(m)=Max[βi-1(m),|Xi(m)|],m=0…N-1,i=2,3,……,
其中|.|表示对频谱取模运算,且β1(m)=|X1(m)|。
4.根据权利要求3所述的一种稳态宽带电磁频谱抑制方法,其特征在于,所述的计算本次比较后保留的最大值与上一次比较后保留的最大值的差值,采用如下公式:
ρi(m)=|βi(m)-βi-1(m)|,i=2,3,……,m=0…N-1,
当ρi(m)<γ时,停止背景信号采集,γ为设定的阈值。
5.根据权利要求1所述的一种稳态宽带电磁频谱抑制方法,其特征在于,所述的根据得到的背景信号频谱,对接收到的频谱数据进行频谱抑制,当接收到的频谱数据的频谱幅度超过背景信号频谱幅度时,则采用当前的频谱值;当接收到的频谱数据的频谱幅度低于背景信号频谱幅度时,采用底噪进行填充,得到抑制背景噪声后的频谱数据,包括如下过程:
当频谱幅度超过背景信号频谱时,采用当前的频谱值;当频谱幅度低于背景信号频谱时,采用底噪进行填充;即:
其中,X(m)为接收到的频谱数据,Xn(m)为拟合出来的频谱底噪,βk(m)为前面得到的背景信号频谱;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110821164.0A CN113541851B (zh) | 2021-07-20 | 2021-07-20 | 一种稳态宽带电磁频谱抑制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110821164.0A CN113541851B (zh) | 2021-07-20 | 2021-07-20 | 一种稳态宽带电磁频谱抑制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113541851A CN113541851A (zh) | 2021-10-22 |
CN113541851B true CN113541851B (zh) | 2022-04-15 |
Family
ID=78129034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110821164.0A Active CN113541851B (zh) | 2021-07-20 | 2021-07-20 | 一种稳态宽带电磁频谱抑制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113541851B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102392924A (zh) * | 2011-09-21 | 2012-03-28 | 中国海洋石油总公司 | 一种抑制管道堵塞声信号中宽带背景噪声的方法 |
CN103886865A (zh) * | 2012-12-20 | 2014-06-25 | 索尼公司 | 声音处理装置、声音处理方法和程序 |
CN106328151A (zh) * | 2015-06-30 | 2017-01-11 | 芋头科技(杭州)有限公司 | 一种环噪消除系统及其应用方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3252252B2 (ja) * | 1996-08-21 | 2002-02-04 | 株式会社山武 | パターン照合装置 |
US20030223597A1 (en) * | 2002-05-29 | 2003-12-04 | Sunil Puria | Adapative noise compensation for dynamic signal enhancement |
US8767974B1 (en) * | 2005-06-15 | 2014-07-01 | Hewlett-Packard Development Company, L.P. | System and method for generating comfort noise |
JP5347794B2 (ja) * | 2009-07-21 | 2013-11-20 | ヤマハ株式会社 | エコー抑圧方法およびその装置 |
CN105842561B (zh) * | 2016-03-28 | 2018-06-26 | 北京航空航天大学 | 一种适用于现场电磁干扰检测的背景信号消除方法 |
CN106157310B (zh) * | 2016-07-06 | 2018-09-14 | 南京汇川图像视觉技术有限公司 | 基于混合自适应水平集模型与多通道结合的TFT LCD mura缺陷检测方法 |
-
2021
- 2021-07-20 CN CN202110821164.0A patent/CN113541851B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102392924A (zh) * | 2011-09-21 | 2012-03-28 | 中国海洋石油总公司 | 一种抑制管道堵塞声信号中宽带背景噪声的方法 |
CN103886865A (zh) * | 2012-12-20 | 2014-06-25 | 索尼公司 | 声音处理装置、声音处理方法和程序 |
CN106328151A (zh) * | 2015-06-30 | 2017-01-11 | 芋头科技(杭州)有限公司 | 一种环噪消除系统及其应用方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113541851A (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101729461B (zh) | 消除单频干扰及多频干扰的系统及方法 | |
CN112737711B (zh) | 一种基于自适应噪声基底估计的宽带载波检测方法 | |
CN108627807B (zh) | 一种机载雷达抗干扰方法 | |
TW200919986A (en) | Jamming detector and jamming detecting method | |
CN102645589B (zh) | 脉冲检测方法及系统 | |
CN113972929B (zh) | 一种高动态多普勒下扩频信号的捕获方法 | |
CN109787656A (zh) | Ofdm电力线通信的自动增益控制装置 | |
CN106772254A (zh) | 改进的基于数字自适应干扰对消的收发隔离方法 | |
CN110855382A (zh) | 基于csi进行人体活动数据检测的方法、装置及电子设备 | |
CN113541851B (zh) | 一种稳态宽带电磁频谱抑制方法 | |
CN108768557B (zh) | 一种从宽带接收信号的频域中检测延时差的方法 | |
CN109709581B (zh) | 一种卫星导航信号强干扰循环特征参数快速检测获取方法 | |
CN108718223B (zh) | 一种非合作信号的盲频谱感知方法 | |
CN113972951B (zh) | 一种低信噪比高动态信号载波多普勒频偏的梳状搜索方法 | |
CN113965215B (zh) | 一种使用高速跳频信号进行单脉冲跟踪解调的方法 | |
CN112630737B (zh) | 对雷达中频回波信号的预处理方法 | |
CN112883787B (zh) | 一种基于频谱匹配的短样本低频正弦信号参数估计方法 | |
CN102508060B (zh) | 一种抗外界脉冲干扰的阵列校准预处理方法 | |
CN108132460B (zh) | 基于频域通道均衡的脉冲压缩补偿算法 | |
US8001167B2 (en) | Automatic BNE seed calculator | |
CN116599606B (zh) | 基于信道化加权互相关处理的扩谱信号接收方法及系统 | |
EP2684073A1 (en) | Radio frequency digital receiver system and method | |
CN112800862B (zh) | 一种非平稳信号时频矩阵重构方法及系统 | |
CN116996137B (zh) | 一种基于加权叠加的低信噪比宽带线性调频信号检测方法 | |
CN115361108B (zh) | 一种超宽带高精度接收时间戳获取方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |