CN114034940A - 一种线谱信号的检测方法 - Google Patents

一种线谱信号的检测方法 Download PDF

Info

Publication number
CN114034940A
CN114034940A CN202111317774.3A CN202111317774A CN114034940A CN 114034940 A CN114034940 A CN 114034940A CN 202111317774 A CN202111317774 A CN 202111317774A CN 114034940 A CN114034940 A CN 114034940A
Authority
CN
China
Prior art keywords
fft
threshold
signal
frequency
line spectrum
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
CN202111317774.3A
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.)
PLA Navy Submarine College
Original Assignee
PLA Navy Submarine College
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 PLA Navy Submarine College filed Critical PLA Navy Submarine College
Priority to CN202111317774.3A priority Critical patent/CN114034940A/zh
Publication of CN114034940A publication Critical patent/CN114034940A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/12Measuring electrostatic fields or voltage-potential
    • G01R29/14Measuring field distribution
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis
    • G01R23/165Spectrum analysis; Fourier analysis using filters

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明涉及一种线谱信号的检测方法。检测方式包括如下步骤:S1对固定时长时间窗t0内的采集得到的低频信号x(t)进行FFT计算,得到其信号的频谱Xfft并去除直流分量,得到XFFT;S2对XFFT进行能量归一化以提高动态观测范围;S3计算XFFT norm的检测门限值,为自适应浮动门限与固定门限相结合;S4比较XFFT norm与Xthreshold的大小,判断疑似线谱频点;S5对T=[T0,T1,L,TN]进行时间积分,若在某时间窗M内疑似线谱频点出现的最大频次Q大于某一设定值R,则判断为疑似目标出现,若连续若干秒时间t1内均出现疑似目标,则判断目标信号出现。它可以实现自适应滤波器的特性变化,准确地确定信号的检测门限。

Description

一种线谱信号的检测方法
技术领域
本发明涉及一种线谱信号的检测方法。
背景技术
舰船电场已成为显著的舰船特征信号,可用于舰船探测、跟踪、定位和打击。俄、美、英、加等国家已在舰船电场研究中取得了相当成就,而国内研究与之差距很大。掌握电场分布的真实特性是电场研究的基础,这对电场测量提出了迫切需求,电场具有明显的频谱特征。
从连续或者离散的输入数据中滤除掉噪声和干扰,提取出其中有用信息的过程称为滤波,实现这个过程的相应装置称为滤波器。一个滤波器可以看作一个系统,这个系统的目的就是从富含噪声的数据中提取有用信号,即期望信号。自适应滤波器的思想就是通过前一时刻的滤波结果自动调整下一时刻的滤波器参数,以适应信号和噪声未知或随时间变化的统记特性,从而实现最优滤波。但由于舰船低频物理场信号的强度随目标和距离的影响较大,在采用能量检测法对低频信号检测时,难以确定信号的检测门限。
发明内容
本发明要解决的技术问题是如何克服难以确定信号的检测门限的问题,因为舰船低频物理场信号的强度随目标和距离的影响较大,在采用能量检测法对低频信号检测时,难以确定信号的检测门限,为解决上述技术问题,本线谱信号的检测方法,其特征是:包括以下步骤,
S1-对固定时长如20s左右的时间窗t0内的采集得到的低频信号x(t)进行FFT计算,得到其信号的频谱Xfft并去除直流分量,得到XFFT
XFFT=[X0,X1,L,XN] 式1)
式1)中,N为频点数,Xi对应的频率为Fi
S2-对XFFT进行能量归一化以提高动态观测范围,舰船由远及近航行时探测节点时,其能量差别可达数十至上百倍,而在20s左右的观测时间内其能量变化不大,得到:
Figure BDA0003344173570000021
S3-计算XFFT norm的检测门限值,为自适应浮动门限与固定门限相结合,其值为:
Figure BDA0003344173570000022
式中,Median(.)为中值滤波运算符(主要根据海洋环境噪声频谱计算自适应浮动门限,中值滤波器的长度L=3、5、...2n-1),Δ=[Δ0,Δ1,L,ΔN]为固定门限向量,相对于传统信号检测算法,由于检测过程中XFFT是实时变化的,从而Xthreshold是变化的,即Xthreshold为自适应门限;
S4-比较XFFT norm与Xthreshold的大小,判断疑似线谱频点
Figure BDA0003344173570000031
式中,T=[T0,T1,L,TN],若Ti大于等于0,则有Yi=0;
S5-对T=[T0,T1,L,TN]进行时间积分,有
Figure BDA0003344173570000032
若在某时间窗M内疑似线谱频点出现的最大频次Q大于某一设定值R,则判断为疑似目标出现,若连续若干秒时间t1内均出现疑似目标,则判断目标信号出现。
优选的,对于连续谱信号,可采用连续累计能量的方法并配合浮动门限进行检测。
本发明一种线谱信号的检测方法的有益效果是:可以实现自适应滤波器的特性变化,准确地确定信号的检测门限从而使得自适应滤波器能够调节其自身传输特性以达到最优。由于信号和噪声的特性是随着时间变化的,所以是无法预知的,所以只使用两种具有固定滤波系数的FIR和IIR滤波器无法跟踪有用信号的变化,没法实现最优滤波。而自适应滤波器不需要信号的先验知识,且计算量小,特别适用于实时处理,它的特性变化是由自适应算法通过调整滤波系数来实现的。
附图说明
下面结合附图对本发明一种线谱信号处理及检测方法作进一步说明:
图1为本发明的流程图;
图2为目标线谱信号图;
具体实施方式
如图1所示,一种线谱信号的检测方法,其特征是:包括以下步骤,
S1-对固定时长如20s左右的时间窗t0内的采集得到的低频信号x(t)进行FFT计算,得到其信号的频谱Xfft并去除直流分量,得到XFFT
XFFT=[X0,X1,L,XN] 式1)
式1)中,N为频点数,Xi对应的频率为Fi
S2-对XFFT进行能量归一化以提高动态观测范围,舰船由远及近航行时探测节点时,其能量差别可达数十至上百倍,而在20s左右的观测时间内其能量变化不大,得到:
Figure BDA0003344173570000041
S3-计算XFFT norm的检测门限值,为自适应浮动门限与固定门限相结合,其值为:
Figure BDA0003344173570000042
式中,Median(.)为中值滤波运算符(主要根据海洋环境噪声频谱计算自适应浮动门限,中值滤波器的长度L=3、5、...2n-1),Δ=[Δ0,Δ1,L,ΔN]为固定门限向量,相对于传统信号检测算法,由于检测过程中XFFT是实时变化的,从而Xthreshold是变化的,即Xthreshold为自适应门限;
S4-比较XFFT norm与Xthreshold的大小,判断疑似线谱频点
Figure BDA0003344173570000051
式中,T=[T0,T1,L,TN],若Ti大于等于0,则有Yi=0;
S5-对T=[T0,T1,L,TN]进行时间积分,有
Figure BDA0003344173570000052
若在某时间窗M内疑似线谱频点出现的最大频次Q大于某一设定值R,则判断为疑似目标出现,若连续若干秒时间t1内均出现疑似目标,则判断目标信号出现。
优选的,对于连续谱信号,可采用连续累计能量的方法并配合浮动门限进行检测。以舰船轴频电场为例、其检测效果如图2所示,图中连续白色线条表示检测到目标线谱信号。
以上描述显示了本发明的主要特征、基本原理,以及本发明的优点。对于本领域技术人员而言,显然本发明不限于上述示范性实施方式或者实施例的细节,且在不背离本发明的精神或者基本特征的情况下,能够以其他的具体形式实现本发明。因此应将上述实施方式或者实施例看作示范性的,且非限制性的。本发明的范围由所附权利要求而非上述说明限定,因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。

Claims (2)

1.一种线谱信号的检测方法,其特征是:包括以下步骤,
S1-对固定时长如20s左右的时间窗t0内的采集得到的低频信号x(t)进行FFT计算,得到其信号的频谱Xfft并去除直流分量,得到XFFT
XFFT=[X0,X1,L,XN] 式1)
式1)中,N为频点数,Xi对应的频率为Fi
S2-对XFFT进行能量归一化以提高动态观测范围,舰船由远及近航行时探测节点时,其能量差别可达数十至上百倍,而在20s左右的观测时间内其能量变化不大,得到:
Figure FDA0003344173560000011
S3-计算XFFT norm的检测门限值,为自适应浮动门限与固定门限相结合,其值为:
Figure FDA0003344173560000012
式中,Median(.)为中值滤波运算符(主要根据海洋环境噪声频谱计算自适应浮动门限,中值滤波器的长度L=3、5、...2n-1),Δ=[Δ0,Δ1,L,ΔN]为固定门限向量,相对于传统信号检测算法,由于检测过程中XFFT是实时变化的,从而Xthreshold是变化的,即Xthreshold为自适应门限;
S4-比较XFFT norm与Xthreshold的大小,判断疑似线谱频点
Figure FDA0003344173560000021
式中,T=[T0,T1,L,TN],若Ti大于等于0,则有Yi=0;
S5-对T=[T0,T1,L,TN]进行时间积分,有
Figure FDA0003344173560000022
若在某时间窗M内疑似线谱频点出现的最大频次Q大于某一设定值R,则判断为疑似目标出现,若连续若干秒时间t1内均出现疑似目标,则判断目标信号出现。
2.根据权利要求1所述的一种线谱信号的检测方法,其特征是:对于连续谱信号,采用连续累计能量的方法并配合浮动门限进行检测。
CN202111317774.3A 2021-11-09 2021-11-09 一种线谱信号的检测方法 Pending CN114034940A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111317774.3A CN114034940A (zh) 2021-11-09 2021-11-09 一种线谱信号的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111317774.3A CN114034940A (zh) 2021-11-09 2021-11-09 一种线谱信号的检测方法

Publications (1)

Publication Number Publication Date
CN114034940A true CN114034940A (zh) 2022-02-11

Family

ID=80136833

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111317774.3A Pending CN114034940A (zh) 2021-11-09 2021-11-09 一种线谱信号的检测方法

Country Status (1)

Country Link
CN (1) CN114034940A (zh)

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵文春等: "基于轴频电场线谱特征的目标检测及识别", 兵 工 学 报, vol. 41, no. 6, pages 1165 - 1171 *

Similar Documents

Publication Publication Date Title
CN112395992B (zh) 一种基于改进小波阈值的电力谐波信号去噪方法
CN109557429A (zh) 基于改进小波阈值去噪的gis局部放电故障检测方法
CN106441288A (zh) 一种加速度计的自适应小波降噪方法
CN112380934B (zh) 基于小波熵与稀疏度的电缆局放信号自适应小波去噪方法
CN111024209B (zh) 一种适用于矢量水听器的线谱检测方法
CN113052000B (zh) 一种船舶机械设备早期微弱故障信号特征诊断方法
CN110085247A (zh) 一种针对复杂噪声环境的双麦克风降噪方法
CN111769810B (zh) 一种基于能量峭度谱的流体机械调制频率提取方法
CN109991520A (zh) 一种电缆振荡波局放检测系统用波速计算新方法
CN114089279A (zh) 一种基于均匀同心圆麦克风阵列的声目标定位方法
CN117213569A (zh) 一种超声波流量检测方法
JP4016115B2 (ja) 船舶に伴う交流水中電界信号又は交流磁気信号抽出方法及び装置
CN109584256B (zh) 一种基于霍夫直线检测的脉冲星色散值估计方法
CN111239565A (zh) 基于分层式去噪模型的充油套管局部放电脉冲信号处理方法及系统
CN114034940A (zh) 一种线谱信号的检测方法
CN117335772A (zh) 一种天文瞬时干扰抑制方法
CN114722885B (zh) 托辊搬运小车运行异常智能检测方法及系统
CN115293219A (zh) 一种融合小波和峭度的脉冲信号去噪方法
CN112652321B (zh) 一种基于深度学习相位更加友好的语音降噪系统及方法
CN113436645A (zh) 复杂噪声环境下的机电系统故障在线监测声学处理方法
Li et al. Adaptive de-noising for PD online monitoring based on wavelet transform
Shen et al. Speech Noise Reduction by EMD-LMS
CN117542369A (zh) 一种对于电力二次电缆寻径的降噪方法及系统
CN111257859B (zh) 一种波束域自相关水下目标识别方法
CN113358927B (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