CN102353952A - 一种频域相干累加的线谱检测方法 - Google Patents

一种频域相干累加的线谱检测方法 Download PDF

Info

Publication number
CN102353952A
CN102353952A CN2011101487001A CN201110148700A CN102353952A CN 102353952 A CN102353952 A CN 102353952A CN 2011101487001 A CN2011101487001 A CN 2011101487001A CN 201110148700 A CN201110148700 A CN 201110148700A CN 102353952 A CN102353952 A CN 102353952A
Authority
CN
China
Prior art keywords
signal
noise
segment
coherent accumulation
data
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.)
Pending
Application number
CN2011101487001A
Other languages
English (en)
Inventor
朴胜春
刘佳
张海刚
袁冬
石珺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN2011101487001A priority Critical patent/CN102353952A/zh
Publication of CN102353952A publication Critical patent/CN102353952A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明提供的是一种频域相干累加的线谱检测方法。依据信号场和噪声场的时间相关特性,对接收信号进行分段处理,使各段信号间存在时延,时延量满足大于噪声的相关半径、小于信号的相关半径;将每个信号段内部再分为若干子数据段,允许每段的数据有部分的交叠;对信号进行时延补偿;第四步,对各段信号累加,得到相干累加的输出信号;第五步,将各个相干累加的输出,进行功率平均得本发明估计出的功率谱。本发明的主要思想是利用信号的相干性,将存在时延差的信号进行相位补偿,则信号形成同相相加,而噪声因为相位是随机的,不会形成同相相加,因而可以提高对微弱的信号检测能力。

Description

一种频域相干累加的线谱检测方法
技术领域
本发明涉及的是一种水声信号处理方法,具体地说是一种线谱信号检测方法。
背景技术
被动声纳是根据目标航行时发出的辐射噪声来检测目标并判别目标类型的。对线谱进行高质量的谱估计,将为目标识别提供重要依据。线谱的检测和提取一直是国内外研究的重点。目前,线谱检测通常是根据计算的周期图的谱值进行检测的。线谱检测方法分为两大类:利用单个时刻谱值的实时检测法和利用多个时刻的数据推迟决策方法。
线谱检测最经典的方法是直接将在某个方向上接收的信号经滤波、解析变换等处理后计算周期图,然后把输出同设定门限进行比较,判断是否存在窄带信号。该方法可以同时保证高检测概率和较低的虚警概率。但由于单个时刻的信息有限,在低信噪比下不能同时得到高的检测概率和低虚警概率。
为了提高在低信噪比下线谱检测、跟踪和提取能力,可以利用多行的数据,推迟决策,即进行LOFAR处理。往往将多种方法进行组合,大体可以分为四类方法:第一类是采用经典的图像处理方法。优点是在谱线为直线时在低信噪比下有很强的提取能力,但当谱线为曲线时效果差,对于宽的谱线还会出现一根谱线重建成多根谱线的情况。第二类采用了经典的目标检测和跟踪理论。LOFAR图中噪声背景下的谱线得到增强,但只能确定各点是否是谱线点,处理后的图像并不能给出信号的任何相对或时变的幅度信息。第三类将线谱检测和提取问题转化为寻优问题。该算法较好地解决了谱线中间出现断点、随机噪声点干扰等线谱检测中的难题,并提高了低信噪比下的线谱检测和提取能力,但是谱线初始状态对算法性能影响较大。第四类模拟人工从LOFAR图中提取谱线的智能线谱检测和提取方法。由于在实现过程中对线谱识别方法进行了优化,并注重在线谱检测的每个环节上设法降低计算量,使算法的计算量较小,适合于实时处理。
发明内容
本发明的目的在于提供一种可以提高对微弱的信号检测能力的频域相干累加的线谱检测方法。
本发明的目的是这样实现的:
第一步,依据信号场和噪声场的时间相关特性,对接收信号x(t)进行分段处理,使各段信号间存在时延,时延量满足大于噪声的相关半径,小于信号的相关半径;
第二步,将每个信号段xi(t)内部再分为若干子数据段xij(t),允许每段的数据有部分的交叠;
第三步,对信号进行时延补偿,首先将信号变换到频域,
xij(t)ββfρXij(jw)
这里Xij(jw)为xij(t)的频谱,再依据时延量在频域对信号进行相位补偿,得到的变换域信号为
Figure BDA0000066196580000021
第四步,对各段信号累加,得到相干累加的输出信号Bj(jw)
Figure BDA0000066196580000022
第五步,将各个相干累加的输出,进行功率平均
Figure BDA0000066196580000023
P(jw)即为本发明估计出的功率谱。
本发明针对水声信号微弱线谱检测问题,提出了一种频域相干累加的方法以提高检测效果。本发明的主要思想是利用信号的相干性,将存在时延差的信号进行相位补偿,则信号形成同相相加,而噪声因为相位是随机的,不会形成同相相加,因而可以提高对微弱的信号检测能力。
附图说明
图1是本发明的流程图;
图2仿真实验相干累加与功率谱的对比图;
图3实验数据相干累加与功率谱的对比图。
具体实施方式
下面结合附图举例对本发明做更详细的描述:
结合图1,本发明的具体内容如下:
第一步,依据信号场和噪声场的时间相关特性,对接收信号x(t)进行分段处理,使各段信号间存在时延,时延量满足大于噪声的相关半径,小于信号的相关半径。设原始接收信号为x(t),每个数据段的延时量τ要小于信号的时间相关半径τs,大于噪声的时间相关半径τn,即
τn TNτTNτs
则第i段数据为
Figure BDA0000066196580000031
Figure BDA0000066196580000032
其中,
Figure BDA0000066196580000033
第二步,为得到平滑的谱估计,将每个信号段xi(t)内部再分为若干子数据段xij(t),允许每段的数据有部分的交叠。例如,若每一段数据重合一半时,这时的段数
Figure BDA0000066196580000034
式中L是段数,M是每段的长度。
第三步,对信号进行时延补偿,首先将信号变换到频域,
xij(t)ββfρXij(jw)
这里Xij(jw)为xij(t)的频谱,再依据时延量在频域对信号进行相位补偿,得到的变换域信号为
Figure BDA0000066196580000035
第四步,对各段信号累加,得到相干累加的输出信号Bj(jw)
Figure BDA0000066196580000036
第五步,为减小噪声的方差,将各个相干累加的输出,进行功率平均
Figure BDA0000066196580000037
P(jw)即为本发明估计出的功率谱。由于噪声相位是随机的,没有和信号一样形成同相相加,因而可以相对的提高信噪比,实现了对微弱信号检测能力的增强。
采用仿真和实验两种方式来验证本发明的方法。
(一)数值仿真
仿真条件:
采用27.2Hz的正弦型信号,噪声为加性带限高斯白噪声,频带为23Hz-35Hz,定义信噪比SNR为带宽内信号与噪声的功率比:
Figure BDA0000066196580000041
在信噪比为-20dB时,对接收信号x(t)进行分段处理,使各段信号间存在时延,选用窗口长度为5s,将每个信号段xi(t)内部再分为若干子数据段xij(t),使每段的数据有一半的交叠,接下来对信号进行时延补偿,首先将信号变换到频域,再依据时延量在频域对信号进行相位补偿,最后对各段信号累加,得到相干累加的输出信号,输出结果如图2所示,同时图2中给出了拥有同样频率分辨率的功率谱估计的结果。
可以看出,在信噪比为-20dB的条件下,原有的功率谱估计已经不能检测到信号,但是利用本发明中的算法信噪比依然有近10dB。
(二)实验数据处理
实验数据为海上实验中测量得到的船舶辐射噪声,实验中的噪声不再是简单的高斯白噪声,采用5分钟的数据进行相干累加处理,结果如图3所示,同时给出功率谱的对比结果。

Claims (3)

1.一种频域相干累加的线谱检测方法,其特征是:
第一步,依据信号场和噪声场的时间相关特性,对接收信号x(t)进行分段处理,使各段信号间存在时延,时延量满足大于噪声的相关半径、小于信号的相关半径;
第二步,将每个信号段xi(t)内部再分为若干子数据段xij(t),允许每段的数据有部分的交叠;
第三步,对信号进行时延补偿,首先将信号变换到频域,
x ij ( t ) → f X ij ( jw )
这里Xij(jw)为xij(t)的频谱,再依据时延量在频域对信号进行相位补偿,得到的变换域信号为
Y ij ( jw ) = X ij ( jw ) e - jw t i ;
第四步,对各段信号累加,得到相干累加的输出信号Bj(jw)
B j ( jw ) = Σ i = 1 M Y ij ( jw ) ;
第五步,将各个相干累加的输出,进行功率平均
P ( jw ) = 1 N Σ j = 1 N | B j ( jw ) | 2
P(jw)即为估计出的功率谱。
2.根据权利要求1所述的一种频域相干累加的线谱检测方法,其特征是所述时延量满足大于噪声的相关半径、小于信号的相关半径是指:原始接收信号为x(t),每个数据段的延时量τ要小于信号的时间相关半径τs,大于噪声的时间相关半径τn,即
τn<τ<τs
则第i段数据为
xi(t)=x(t-τi)
τi=(i-1)τ
其中,i=1,2,·,M。
3.根据权利要求1或2所述的一种频域相干累加的线谱检测方法,其特征是所述允许每段的数据有部分的交叠是指:若每一段数据重合一半时,这时的段数
Figure FDA0000066196570000021
式中L是段数,M是每段的长度。
CN2011101487001A 2011-06-03 2011-06-03 一种频域相干累加的线谱检测方法 Pending CN102353952A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101487001A CN102353952A (zh) 2011-06-03 2011-06-03 一种频域相干累加的线谱检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101487001A CN102353952A (zh) 2011-06-03 2011-06-03 一种频域相干累加的线谱检测方法

Publications (1)

Publication Number Publication Date
CN102353952A true CN102353952A (zh) 2012-02-15

Family

ID=45577546

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101487001A Pending CN102353952A (zh) 2011-06-03 2011-06-03 一种频域相干累加的线谱检测方法

Country Status (1)

Country Link
CN (1) CN102353952A (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237871A (zh) * 2013-06-08 2014-12-24 中国科学院声学研究所 一种基于相位补偿的时延差估计方法
CN105277934A (zh) * 2015-09-24 2016-01-27 哈尔滨工程大学 一种基于阵列的弱线谱目标被动检测方法
CN105785346A (zh) * 2014-12-26 2016-07-20 中国科学院声学研究所 一种基于相位方差加权的未知目标线谱检测方法及系统
CN108732558A (zh) * 2018-03-26 2018-11-02 东南大学 一种基于拷贝信号分段的匹配滤波计算方法
CN109061591A (zh) * 2018-07-23 2018-12-21 东南大学 一种基于序贯聚类的时频线谱检测方法
CN109655148A (zh) * 2018-12-19 2019-04-19 南京世海声学科技有限公司 一种船舶噪声非平稳低频线谱自主提取方法
CN109743750A (zh) * 2018-12-20 2019-05-10 中国电子科技集团公司电子科学研究院 一种视频泄漏信号去噪算法及装置
CN111220265A (zh) * 2019-12-10 2020-06-02 哈尔滨工程大学 一种船舶辐射噪声谐波线谱增强方法
CN111368679A (zh) * 2020-02-26 2020-07-03 西北工业大学 一种低秩矩阵分解的谱线检测方法
CN112019284A (zh) * 2020-08-27 2020-12-01 中电科仪器仪表有限公司 一种低信噪比下的窄带信号时差计算方法及系统
CN112098984A (zh) * 2020-11-02 2020-12-18 南昌大学 散射信号的目标体检测方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001076086A2 (en) * 2000-03-29 2001-10-11 Time Domain Corporation System and method of using multiple correlator receivers in an impulse radio system
CN101231315A (zh) * 2007-01-24 2008-07-30 涂亚庆 频率估计的多段采样信号融合处理方法
CN101494508A (zh) * 2009-02-26 2009-07-29 上海交通大学 基于特征循环频率的频谱检测方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001076086A2 (en) * 2000-03-29 2001-10-11 Time Domain Corporation System and method of using multiple correlator receivers in an impulse radio system
CN101231315A (zh) * 2007-01-24 2008-07-30 涂亚庆 频率估计的多段采样信号融合处理方法
CN101494508A (zh) * 2009-02-26 2009-07-29 上海交通大学 基于特征循环频率的频谱检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《浙江大学学报(工学版)》 20071231 刘辉涛等 窄带弱信号的线谱检测__相干累加频域批处理自适应线谱增强方法 第41卷, 第12期 *
刘辉涛等: "窄带弱信号的线谱检测――相干累加频域批处理自适应线谱增强方法", 《浙江大学学报(工学版)》, vol. 41, no. 12, 31 December 2007 (2007-12-31) *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104237871B (zh) * 2013-06-08 2017-01-11 中国科学院声学研究所 一种基于相位补偿的时延差估计方法
CN104237871A (zh) * 2013-06-08 2014-12-24 中国科学院声学研究所 一种基于相位补偿的时延差估计方法
CN105785346A (zh) * 2014-12-26 2016-07-20 中国科学院声学研究所 一种基于相位方差加权的未知目标线谱检测方法及系统
CN105785346B (zh) * 2014-12-26 2017-12-29 中国科学院声学研究所 一种基于相位方差加权的未知目标线谱检测方法及系统
CN105277934A (zh) * 2015-09-24 2016-01-27 哈尔滨工程大学 一种基于阵列的弱线谱目标被动检测方法
CN108732558B (zh) * 2018-03-26 2021-08-24 东南大学 一种基于拷贝信号分段的匹配滤波计算方法
CN108732558A (zh) * 2018-03-26 2018-11-02 东南大学 一种基于拷贝信号分段的匹配滤波计算方法
CN109061591A (zh) * 2018-07-23 2018-12-21 东南大学 一种基于序贯聚类的时频线谱检测方法
CN109655148A (zh) * 2018-12-19 2019-04-19 南京世海声学科技有限公司 一种船舶噪声非平稳低频线谱自主提取方法
CN109743750A (zh) * 2018-12-20 2019-05-10 中国电子科技集团公司电子科学研究院 一种视频泄漏信号去噪算法及装置
CN109743750B (zh) * 2018-12-20 2022-12-16 中国电子科技集团公司电子科学研究院 一种视频泄漏信号去噪算法及装置
CN111220265A (zh) * 2019-12-10 2020-06-02 哈尔滨工程大学 一种船舶辐射噪声谐波线谱增强方法
CN111368679A (zh) * 2020-02-26 2020-07-03 西北工业大学 一种低秩矩阵分解的谱线检测方法
CN111368679B (zh) * 2020-02-26 2022-03-22 西北工业大学 一种低秩矩阵分解的谱线检测方法
CN112019284A (zh) * 2020-08-27 2020-12-01 中电科仪器仪表有限公司 一种低信噪比下的窄带信号时差计算方法及系统
CN112098984A (zh) * 2020-11-02 2020-12-18 南昌大学 散射信号的目标体检测方法及装置
CN112098984B (zh) * 2020-11-02 2021-02-23 南昌大学 散射信号的目标体检测方法及装置

Similar Documents

Publication Publication Date Title
CN102353952A (zh) 一种频域相干累加的线谱检测方法
CN104569948B (zh) 海杂波背景下子带自适应glrt‑ltd检测方法
CN110376575B (zh) 一种基于阻尼参数匹配随机共振的低频线谱检测方法
CN106597408B (zh) 基于时频分析和瞬时频率曲线拟合的高阶pps信号参数估计方法
CN103559887B (zh) 用于语音增强系统的背景噪声估计方法
CN102967854A (zh) Frft域海杂波中目标的多重分形检测方法
CN103995950A (zh) 基于空域相关修正阈值的变小波系数局部放电信号消噪方法
CN105223576A (zh) 一种基于单矢量潜标的线谱信号目标自动检测方法
CN101900601B (zh) 复杂多途水声环境下直达声辨识方法
CN106301289A (zh) 利用自适应滤波算法消除泥浆脉冲信号中的泵冲噪声的方法
CN102636775B (zh) 基于模糊逻辑识别的风廓线雷达回波谱重构方法
CN104614769A (zh) 一种压制地震面波的聚束滤波方法
CN103197347A (zh) 一种基于自适应时窗的吸收分析油气预测方法
CN101581787B (zh) 基于最小均方算法的激光测云雷达信号自适应识别方法
CN103905656A (zh) 残留回声的检测方法及装置
CN114785379A (zh) 一种水声janus信号参数估计方法及系统
CN102546052B (zh) 一种自适应门限控制信号参数统计测量方法
CN116165610A (zh) 一种改进奇异值分解的海杂波抑制算法
CN111090089B (zh) 一种基于两类辅助数据的空时自适应检测方法
CN102571671B (zh) 一种基于mspwvd的prbc-lfm复合信号伪码序列盲估计方法
US8127138B1 (en) Method for embedding information in sonar
CN111090088B (zh) 基于主动声纳阵列斜对称特性的泄露目标空时检测方法
CN110542927B (zh) 变窗口加权地震数据尖峰噪声压制方法
Ma et al. Time of arrival estimation for underwater acoustic signal using multi-feature fusion
CN116312623A (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
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120215