CN115015942A - 一种自适应水下目标声学测速装置及方法 - Google Patents
一种自适应水下目标声学测速装置及方法 Download PDFInfo
- Publication number
- CN115015942A CN115015942A CN202210947918.1A CN202210947918A CN115015942A CN 115015942 A CN115015942 A CN 115015942A CN 202210947918 A CN202210947918 A CN 202210947918A CN 115015942 A CN115015942 A CN 115015942A
- Authority
- CN
- China
- Prior art keywords
- frequency
- signal
- spectrum
- module
- frequency band
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005259 measurement Methods 0.000 title claims abstract description 12
- 238000001228 spectrum Methods 0.000 claims abstract description 138
- 238000000691 measurement method Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 55
- 238000012545 processing Methods 0.000 claims description 33
- 230000009467 reduction Effects 0.000 claims description 33
- 238000012937 correction Methods 0.000 claims description 28
- 230000001133 acceleration Effects 0.000 claims description 22
- 238000010586 diagram Methods 0.000 claims description 19
- 238000007670 refining Methods 0.000 claims description 17
- 230000003044 adaptive effect Effects 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 12
- 238000009499 grossing Methods 0.000 claims description 9
- 239000000654 additive Substances 0.000 claims description 7
- 230000000996 additive effect Effects 0.000 claims description 7
- 230000005855 radiation Effects 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 10
- 238000001514 detection method Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000035484 reaction time Effects 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
Images
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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- 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
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
-
- 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种自适应水下目标声学测速装置及方法,包括线谱多频段估计模块,自适应降噪模块,频谱细化模块,频谱校正模块,多普勒频移计算模块,水下目标参数计算模块,分批解算模块。本发明的装置及方法能够提高水下目标声学测速的自适应性能、频谱分辨精度和实时性。对于以任意频率采样所得的输入信号均能输出频谱分辨率为0.01Hz的频谱图,极大提高了通用性与适应性;可自动对目标进行分批处理,将输出结果中速度相同而频率不同的信号来源划分为同一集合,提高了对水下航行器的分辨能力及抗干扰能力。
Description
技术领域
本发明涉及海洋测量技术领域,特别是涉及一种自适应水下目标声学测速装置及方法。
背景技术
声波是目前唯一能在海洋中远距离传播的信息载体,也是水下探测所依赖的媒介。因此,水声探测技术逐步成为了世界各国研究和开发的重点。该技术通过对被动接受的水声辐射信号进行运算处理来获得水下目标的运动参数,如速度、加速度等。而后可以对水下目标的运动参数进行估计,得到关于目标的运动分析结果,并对目标的运动轨迹进行拟合跟踪。低频段水声被动监测技术相较于其他水声探测技术具有:
依托于水听器的被动式监测方法相较于传统的主动式监测方法,具有隐蔽性好,适应性强等特点。
针对水下航行器难以消除的低频段上的辐射噪声进行过滤提取,有效性与针对性更强。
虽然现有的低频段水声被动监测技术对比以往水声探测技术有一定优势,但是相较于计算机技术以及水下航行器技术的充分发展,现有的被动水声探测技术仍存在反应时间长,自动化程度低,通用性弱等缺点。不同于实验室的理想环境,在实际环境中,过长的反应时间会给水下目标更大的逃逸可能。同时现代技术发展的方向也在朝着自动化、标准化前进。故提出一种可解决自适应性能、频谱分辨精度、实时性等方面问题的方法迫在眉睫。
发明内容
针对现有技术存在的问题,本发明的目的在于针对自适应性能、频谱分辨精度、实时性等方面问题,通过对单通道采集信号进行自适应解算,提出一种自适应水下目标声学测速方法及装置。该方法以水听器的输出信号为输入,对信号进行降噪处理,从中提取出包含目标信号噪源的多个频段,并进行频谱细化,得到统一的高精度频谱;最后利用多普勒频移计算公式、速度计算公式、加速度计算公式、分批解算方法等解算出水下目标的相对径向速度、加速度、分批结果等目标参数。
为实现上述目的,本发明提供一种自适应水下目标声学测速装置,所述装置包括线谱多频段估计模块、自适应降噪模块、频谱细化模块、频谱校正模块、多普勒频移计算模块、水下目标参数计算模块和分批解算模块;其中,
所述线谱多频段估计模块,用于接收原始输入信号并将包含目标辐射噪声的频段选取出来;
所述自适应降噪模块,用于接收原始输入信号并自适应地进行降噪处理;
所述频谱细化模块,用于接收目标频段选取结果与降噪后的信号;
所述频谱校正模块,用于接收频谱细化后输出的信号,并将校正后的频谱信号并输出;
所述多普勒频移计算模块,用于接收频谱校正模块输出的信号,对该信号进行自适应的平滑处理,最后将频偏公式带入平滑后的多普勒频移曲线,计算出多普勒频偏并输出;
所述水下目标参数计算模块,用于接收目标信号的多普勒频移,根据速度和加速度计算公式计算并输出目标的速度和加速度;
所述分批解算模块,用于接收目标信号的速度,根据不同信号源的速度,将同一速度不同频率的信号源划分为同一目标,输出分批结果。
进一步,所述线谱多频段估计模块输入为WAV格式原始信号,输出与频谱细化模块相连接;所述自适应降噪模块输入为WAV格式原始信号,输出与频谱细化模块输入相连接;所述频谱细化模块输出与频谱校正模块输入相连接;所述频谱校正模块输出与多普勒频移计算模块输入相连接;所述多普勒频移计算模块输出与水下目标参数计算模块输入相连接;所述水下目标参数计算模块输出为水下目标的速度、加速度并与分批解算模块输入相连接。
进一步,所述线谱多频段估计模块将原始信号通过STFT后得到时频图,并采用自适应的能量门限,得到一系列线谱频率值,利用多普勒频移计算公式拓宽其上下限,得到一系列的第一频段集合;再取时频图作时域积分,并采用自适应的能量门限,得到一系列的线谱频率值,利用多普勒频移计算公式拓宽其上下限,得到一系列的第二频段集合;取第一频段集合和第二频段集合的频段交集作为目标频段选取结果输出。
进一步,自适应降噪模块先将原始信号通过截止频率为100Hz的高通滤波器,再利用变分模态分解的方法、小波阈值的方法或时域积分降噪的方法对信号进行降噪处理,过滤掉加性噪声后输出。
进一步,频谱细化模块用于接收目标频段选取结果与降噪后的信号,以目标频段选取结果作为上下限,将降噪后的信号通过Chirp-Z变换进行频谱细分或对降噪后的信号选用ZOOM-FFT进行频谱细化,输出频谱精度为0.01Hz的频谱,且对于任意频率的采集信号均能输出频谱精度为0.01Hz的信号。
进一步,频谱校正模块用于接收频谱细化后输出的信号,并依次找出该信号中各频段的功率峰值点作为中心频率点,将中心频率点左右0.5Hz范围内的频段选取出来,作为校正后的频谱信号并输出。
进一步,多普勒频移计算模块用于接收频谱校正模块输出的信号,对该信号进行自适应的平滑处理,最后将频偏公式带入平滑后的多普勒频移曲线,计算出多普勒频偏并输出。
另一方面,本发明提供一种自适应水下目标声学测速方法,所述方法包括以下步骤:
S1. 接收WAV格式的原始信号,输入到线谱多频段估计模块中,输出目标频段选取结果;
S2. 将原始信号输入到自适应降噪模块中,输出降噪后的信号;
S3. 将目标频段选取结果与降噪后的信号输入到频谱细化模块中,输出细化后的时频谱;
S4. 将细化后的频谱图输入到频谱校正模块,输出频谱校正后的时频图;
S5. 将频谱校正后的时频图输入到多普勒频移计算模块,输出多普勒频移曲线图;
S6. 将多普勒频移曲线图输入到水下目标参数计算模块,输出水下目标的速度和加速度;
S7. 根据水下目标的速度解算出分批情况,输出分批结果。
进一步,步骤S1中,还包括以下子步骤:
S11. 接收由水听器输出的WAV格式的信号,是一段包含目标辐射噪声在内的时域信号;
S12. 线谱多频段估计模块对接收的信号进行STFT处理,处理完成后保存处理结果,其中
S14. 按照频段选取公式对第一频率值集合进行排序分组,将同一组的频率极值作为该频段的上下限,得到一系列的第一频段集合;
式中为频段内的频率值,为中心频率,为当前时刻目标的运动速度,为目标与声纳夹角,为声速;得到的每一个频段的频率上限为,频率下限为;根据多普勒频移计算公式,取中心频率为,拓宽频段下界,取中心频率为,拓宽频段上界;
S16. 按照频段选取公式对第二频率值集合进行排序分组,将同一组的频率极值作为该频段的上下限,得到一系列的第二频段集合;
S17. 取第一频段集合与第二频段集合的交集作为最终的目标频段。
进一步,步骤S3中,对降噪后的频谱图进行Chirp-Z变换,以目标频段作为Chirp-Z变换的频率选择范围,得到变换后的细化频谱;
Chirp-Z变换的公式为:
式中,为离散序列的Chirp-Z变换,为数据长度,,,表示起始采样点的矢量半径长度,表示起始采样点的相角,表示两相邻采样点之间的角度差,的大小表示螺旋线的伸展率,为虚数单位,e为自然数,k为Z变换中的序列,n为离散序列x(n)的序列。
与现有技术相比,本装置及方法具有如下优点:
对于以任意频率采样所得的输入信号均能输出频谱分辨率为0.01Hz的频谱图,极大提高了通用性与适应性,为后续标准化处理奠定了基础。
可自动对目标进行分批处理,将输出结果中速度相同而频率不同的信号来源划分为同一集合,提高了对水下航行器的分辨能力及抗干扰能力。
附图说明
图1示出了根据本发明实施例中水下目标声学测速装置和方法的示意性流程图;
图2示出了根据本发明实施例中对信号进行STFT处理后的时频图;
图3示出了根据本发明实施例中经过降噪处理后得到的频谱图;
图4示出了根据本发明实施例中经过Chirp-Z变换后得到的时频图;
图5示出了根据本发明实施例中平滑处理后的频谱图;
图6示出了根据本发明实施例中多普勒频移曲线图;
图7示出了根据本发明实施例中该频段所求得全时段移动速度曲线图;
图8示出了根据本发明实施例中该频段所求得全时段加速度曲线图。
具体实施方式
下面将结合附图,对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
以下结合图1-图8对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
如图1所示,根据本发明的第一方面,提供了一种自适应水下目标声学测速装置,该装置包括线谱多频段估计模块1、自适应降噪模块2、频谱细化模块3、频谱校正模块4、多普勒频移计算模块5、水下目标参数计算模块6、分批解算模块7。所述线谱多频段估计模块1输入为WAV格式原始信号,输出与频谱细化模块3相连接;所述自适应降噪模块2输入为WAV格式原始信号,输出与频谱细化模块3输入相连接;所述频谱细化模块输出3与频谱校正模块4输入相连接;所述频谱校正模块4输出与多普勒频移计算模块5输入相连接;所述多普勒频移计算模块5输出与水下目标参数计算模块6输入相连接;所述水下目标参数计算模块6输出为水下目标的速度、加速度并与分批解算模块7输入相连接。图1为上述水下目标声学测速装置和方法的示意性流程图。
其中,线谱多频段估计模块1,用于接收原始输入信号并将包含目标辐射噪声的频段截选出来。将原始信号通过STFT后得到时频图,并采用自适应的能量门限,得到一系列线谱频率值,利用多普勒频移计算公式拓宽其上下限,得到一系列的第一频段集合。再取时频图作时域积分,并采用自适应的能量门限,得到一系列的线谱频率值,利用多普勒频移计算公式拓宽其上下限,得到一系列的第二频段集合。取第一频段集合和第二频段集合的频段交集作为目标频段选取结果输出。
自适应降噪模块2,用于接收原始输入信号并自适应地进行降噪处理。由于水听器中电子元器件所引发的噪声大多集中在100Hz以下的频段内,故先将原始信号通过截止频率为100Hz的高通滤波器,再利用VMD(变分模态分解)的方法对信号进行降噪处理,过滤掉大多数加性噪声后输出。
作为可替换方式,在原始信号通过高通滤波器之后,利用小波阈值的方法对信号进行降噪处理,过滤掉大多数加性噪声后输出。
作为另一种可替换方式,在原始信号通过高通滤波器之后,利用时域积分降噪的方法对信号进行降噪处理,过滤掉大多数加性噪声后输出。
频谱细化模块3,用于接收目标频段选取结果与降噪后的信号,为了得到更高精度的频谱并提高通用性,以目标频段选取结果作为上下限,将降噪后的信号通过Chirp-Z变换进行频谱细分,输出频谱精度为0.01Hz的频谱。且对于任意频率的采集信号均能输出频谱精度为0.01Hz的信号。
作为可替换方式,对降噪后的信号选用ZOOM-FFT进行频谱细化。
频谱校正模块4,用于接收频谱细化后输出的信号,并依次找出该信号中各频段的功率峰值点作为中心频率点,将中心频率点左右0.5Hz范围内的频段截选出来,作为校正后的频谱信号并输出。
作为可替换方式,利用正弦波信号频率估计方法选择中心频率点。
多普勒频移计算模块5,用于接收频谱校正模块输出的信号,对该信号进行自适应的平滑处理,最后将频偏公式带入平滑后的多普勒频移曲线,计算出多普勒频偏并输出。
水下目标参数计算模块6,用于接收目标信号的多普勒频移。根据速度、加速度计算公式计算并输出目标的速度、加速度。
分批解算模块7,用于接收目标信号的速度。根据不同信号源的速度,将同一速度不同频率的信号源划分为同一目标,输出分批结果。
根据本发明的第二方面,提供了一种基于多普勒频移自适应水下目标测速方法,包括如下步骤:
S1.接收WAV格式的原始信号,输入到线谱多频段估计模块1中,输出目标频段选取结果。
S2.将原始信号输入到自适应降噪模块2中,输出降噪后的信号。
S3.将目标频段选取结果与降噪后的信号输入到频谱细化模块3中,输出细化后的时频谱。
S4.将细化后的频谱图输入到频谱校正模块4,输出频谱校正后的时频图。
S5.将频谱校正后的时频图输入到多普勒频移计算模块5,输出多普勒频移曲线图。
S6.将多普勒频移曲线图输入到水下目标参数计算模块6,输出水下目标的速度、加速度。
S7.通过水下目标参数模块和分批解算模块根据水下目标的速度解算出分批情况,输出分批结果。
在图2-图8示出的一个具体实施例中,对SWellEx-96实验中的Event S5事件中垂直阵列(vertical line arrays, VLA)采集到的数据进行处理与分析,实施步骤如下:
1) 接收到由水听器输出的WAV格式的信号,是一段包含目标辐射噪声在内的时域信号。
2) 对接收的信号进行STFT处理。进一步的,处理完成后保存处理结果。
图2为对信号进行STFT处理后的时频图。
4) 按照频段选取公式对第一频率值集合进行排序分组,将同一组的频率极值作为该频段的上下限,得到一系列的第一频段集合。
式中为频段内的频率值,为中心频率,为当前时刻目标的运动速度,为目标与声纳夹角,为声速。假设得到的每一个频段的频率上限为,频率下限为。在此基础上,考虑到受多普勒频移影响,线谱会在频域上会有一定展宽,为获得尽可能准确的频段,同样根据多普勒频移计算公式,取中心频率为,拓宽频段下界,取中心频率为,拓宽频段上界。
6) 按照频段选取公式对第二频率值集合进行排序分组,将同一组的频率极值作为该频段的上下限,得到一系列的第二频段集合。
7) 取第一频段集合与第二频段集合的交集作为最终的目标频段。
最终的目标频段如表1所示。
8)对接收的目标信号使用VMD(变分模态分解)法进行降噪,得到降噪后的频谱结果。
在其他的一些示例中,对接收的目标信号使用小波阈值法进行降噪,得到降噪后的频谱结果。
小波阈值法为:
图3为经过降噪处理后得到的频谱图。
9) 对降噪后的频谱图进行Chirp-Z变换,以目标频段作为Chirp-Z变换的频率选择范围,得到变换后的细化频谱。
Chirp-Z变换的公式为:
式中,为离散序列的Chirp-Z变换,为数据长度, , , 表示起始采样点 的矢量半径长度, 表示起始采样点 的相角, 表示两相邻采样点之间的角度差, 的大小表示螺旋线的伸展率, 为虚数单位,e为自然数,k为Z变换中的序列,n为离散序列x(n)的序列。
在其他的一些示例中,对降噪后的频谱图进行ZOOM-FFT,得到变换后的细化频谱。
图4为经过Chirp-Z变换后得到的时频图。
10) 对经过Chirp-Z变换后得到的频谱图上各频段内。找出各时刻每一个目标频段上的最大值对应频率,作为该时刻线谱的瞬时频率。再对该线谱进行加权平均峰值处理,并进行自适应的平滑处理,得到平滑处理后的不同频段内的线谱的瞬时频率。
由于频段较多,在此仅取表1中频段6作为示例,其余频段计算过程与方法可参考该频段相关计算示例。
图5为平滑处理后的频谱图。
11) 在平滑处理后的线谱上利用以功率峰值点作为各频段的中心频率点。再由各频段的中心频率点及频移公式计算得到各频段内信号源的频偏。
在其他的一些示例中,利用一种正弦波信号频率估计方法选择中心频率。再由各频段的中心频率点及频移公式计算得到各频段内信号源的频偏。
正弦波信号频率估计方法:
图6为多普勒频移曲线图。
12) 由各信号源的频偏及速度计算公式计算依次得到各信号源的移动速度集合。
图7为该频段所求得全时段移动速度曲线图。
13) 根据各信号源的移动速度集合计算得到各信号源的加速度集合。
图8为该频段所求得全时段加速度曲线图。
14) 根据各信号源的移动速度集合及加速度集合对信号源进行分批处理,得到分批处理结果。
根据各信号源移动速度及加速度分析可知,十一个频段内的信号源来自同一水下目标,分批后抗干扰效果显著。
在本说明书的描述中,参考术语“实施例”、“实施例”等的描述意指结合该实施例或实施例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或实施例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或实施例。此外,本领域的技术人员可以在不产生矛盾的情况下,将本说明书中描述的不同实施例或实施例以及其中的特征进行结合或组合。
上述内容虽然已经示出和描述了本发明的实施例,可以理解的是,上述实施例是实施例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型等更新操作。
Claims (10)
1.一种自适应水下目标声学测速装置,其特征在于,所述装置包括线谱多频段估计模块、自适应降噪模块、频谱细化模块、频谱校正模块、多普勒频移计算模块、水下目标参数计算模块和分批解算模块;其中,
所述线谱多频段估计模块,用于接收原始输入信号并选取包含目标辐射噪声线谱成分的频段;
所述自适应降噪模块,用于接收原始输入信号并自适应地进行降噪处理;
所述频谱细化模块,用于接收目标频段选取结果与降噪后的信号;
所述频谱校正模块,用于接收频谱细化后输出的信号,并输出校正后的频谱信号;
所述多普勒频移计算模块,用于接收频谱校正模块输出的信号,对该信号进行自适应的平滑处理,同时将频偏公式带入平滑后的多普勒频移曲线,计算出多普勒频偏并输出;
所述水下目标参数计算模块,用于接收目标信号的多普勒频移,根据速度和加速度计算公式计算并输出目标的速度和加速度;
所述分批解算模块,用于接收目标信号的速度,根据不同信号源的速度,将同一速度不同频率的信号源划分为同一目标,输出分批结果。
2.根据权利要求1所述的自适应水下目标声学测速装置,其特征在于,所述线谱多频段估计模块输入为WAV格式原始信号,输出与频谱细化模块相连接;所述自适应降噪模块输入为WAV格式原始信号,输出与频谱细化模块输入相连接;所述频谱细化模块输出与频谱校正模块输入相连接;所述频谱校正模块输出与多普勒频移计算模块输入相连接;所述多普勒频移计算模块输出与水下目标参数计算模块输入相连接;所述水下目标参数计算模块输出为水下目标的速度、加速度并与分批解算模块输入相连接。
3.根据权利要求1或2所述的自适应水下目标声学测速装置,其特征在于,所述线谱多频段估计模块将原始信号通过STFT后得到时频图,并采用自适应能量门限,得到一系列线谱频率值,利用多普勒频移计算公式拓宽其上下限,得到一系列的第一频段集合;再取时频图作时域积分,并采用自适应的能量门限,得到一系列的线谱频率值,利用多普勒频移计算公式拓宽其上下限,得到一系列的第二频段集合;取第一频段集合和第二频段集合的频段交集作为目标频段选取结果输出。
4.根据权利要求3所述的自适应水下目标声学测速装置,其特征在于,自适应降噪模块先将原始信号通过截止频率为100Hz的高通滤波器,再利用变分模态分解的方法、小波阈值的方法或时域积分降噪的方法对信号进行降噪处理,过滤掉加性噪声后输出。
5.根据权利要求4所述的自适应水下目标声学测速装置,其特征在于,频谱细化模块用于接收目标频段选取结果与降噪后的信号,以目标频段选取结果作为上下限,将降噪后的信号通过Chirp-Z变换进行频谱细分或对降噪后的信号选用ZOOM-FFT进行频谱细化,输出频谱精度为0.01Hz的频谱,且对于任意频率的采集信号均能输出频谱精度为0.01Hz的信号。
6.根据权利要求5所述的自适应水下目标声学测速装置,其特征在于,频谱校正模块用于接收频谱细化后输出的信号,并依次找出该信号中各频段的功率峰值点作为中心频率点,选取出中心频率点左右0.5Hz范围内的频段,作为校正后的频谱信号并输出。
7.根据权利要求6所述的自适应水下目标声学测速装置,其特征在于,多普勒频移计算模块用于接收频谱校正模块输出的信号,对该信号进行自适应的平滑处理,将频偏公式带入平滑后的多普勒频移曲线,即可计算出多普勒频偏并输出。
8.一种自适应水下目标声学测速方法,其特征在于,所述方法包括以下步骤:
S1. 接收WAV格式的原始信号,输入到线谱多频段估计模块中,输出目标频段选取结果;
S2. 将原始信号输入到自适应降噪模块中,输出降噪后的信号;
S3. 将目标频段选取结果与降噪后的信号输入到频谱细化模块中,输出细化后的时频谱;
S4. 将细化后的频谱图输入到频谱校正模块,输出频谱校正后的时频图;
S5. 将频谱校正后的时频图输入到多普勒频移计算模块,输出多普勒频移曲线图;
S6. 将多普勒频移曲线图输入到水下目标参数计算模块,输出水下目标的速度和加速度;
S7. 根据水下目标的速度解算出分批情况,输出分批结果。
9.根据权利要求8所述的自适应水下目标声学测速方法,其特征在于,所述方法中,步骤S1中,还包括以下子步骤:
S11. 接收到由水听器输出的WAV格式的信号,是一段包含目标辐射噪声在内的时域信号;
S12. 线谱多频段估计模块对接收的信号进行STFT处理,处理完成后保存处理结果;
式中STFT的变换公式为:
S14. 按照频段选取公式对第一频率值集合进行排序分组,将同一组的频率极值作为该频段的上下限,得到一系列的第一频段集合;
式中为频段内的频率值,为中心频率,为当前时刻目标的运动速度,为目标与声纳夹角,为声速;得到的每一个频段的频率上限为,频率下限为;根据多普勒频移计算公式,取中心频率为,拓宽频段下界,取中心频率为,拓宽频段上界;
S16. 按照频段选取公式对第二频率值集合进行排序分组,将同一组的频率极值作为该频段的上下限,得到一系列的第二频段集合;
S17. 取第一频段集合与第二频段集合的交集作为最终的目标频段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210947918.1A CN115015942B (zh) | 2022-08-09 | 2022-08-09 | 一种自适应水下目标声学测速装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210947918.1A CN115015942B (zh) | 2022-08-09 | 2022-08-09 | 一种自适应水下目标声学测速装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115015942A true CN115015942A (zh) | 2022-09-06 |
CN115015942B CN115015942B (zh) | 2022-11-01 |
Family
ID=83066220
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210947918.1A Active CN115015942B (zh) | 2022-08-09 | 2022-08-09 | 一种自适应水下目标声学测速装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115015942B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116973901A (zh) * | 2023-09-14 | 2023-10-31 | 海底鹰深海科技股份有限公司 | 时频分析在声呐信号处理中的算法应用 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120169542A1 (en) * | 2006-04-28 | 2012-07-05 | Mathews Michael B | System and method for positioning using hybrid spectral compression and cross correlation signal processing |
US20200066132A1 (en) * | 2018-08-24 | 2020-02-27 | Digital Global Systems, Inc. | Systems, Methods, and Devices for Automatic Signal Detection based on Power Distribution by Frequency over Time |
CN111123286A (zh) * | 2020-01-12 | 2020-05-08 | 哈尔滨理工大学 | 一种基于自校准全光纤多普勒电缆长度测试方法及装置 |
CN111505618A (zh) * | 2020-04-14 | 2020-08-07 | 东南大学 | 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法 |
CN111750972A (zh) * | 2020-05-26 | 2020-10-09 | 北京理工大学 | 一种适用于双频连续激光雷达的目标振动测量方法 |
CN112526474A (zh) * | 2020-11-23 | 2021-03-19 | 哈尔滨工程大学 | 基于全相位傅里叶变换的fmcw雷达距离速度联合估计方法 |
CN114862809A (zh) * | 2022-05-19 | 2022-08-05 | 杭州意能电力技术有限公司 | 一种基于移动终端与图像处理的振动监测方法和装置 |
-
2022
- 2022-08-09 CN CN202210947918.1A patent/CN115015942B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120169542A1 (en) * | 2006-04-28 | 2012-07-05 | Mathews Michael B | System and method for positioning using hybrid spectral compression and cross correlation signal processing |
US20200066132A1 (en) * | 2018-08-24 | 2020-02-27 | Digital Global Systems, Inc. | Systems, Methods, and Devices for Automatic Signal Detection based on Power Distribution by Frequency over Time |
CN111123286A (zh) * | 2020-01-12 | 2020-05-08 | 哈尔滨理工大学 | 一种基于自校准全光纤多普勒电缆长度测试方法及装置 |
CN111505618A (zh) * | 2020-04-14 | 2020-08-07 | 东南大学 | 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法 |
CN111750972A (zh) * | 2020-05-26 | 2020-10-09 | 北京理工大学 | 一种适用于双频连续激光雷达的目标振动测量方法 |
CN112526474A (zh) * | 2020-11-23 | 2021-03-19 | 哈尔滨工程大学 | 基于全相位傅里叶变换的fmcw雷达距离速度联合估计方法 |
CN114862809A (zh) * | 2022-05-19 | 2022-08-05 | 杭州意能电力技术有限公司 | 一种基于移动终端与图像处理的振动监测方法和装置 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116973901A (zh) * | 2023-09-14 | 2023-10-31 | 海底鹰深海科技股份有限公司 | 时频分析在声呐信号处理中的算法应用 |
Also Published As
Publication number | Publication date |
---|---|
CN115015942B (zh) | 2022-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110307994B (zh) | 异声检测装置以及异声检测方法 | |
CN107505135B (zh) | 一种滚动轴承复合故障提取方法及系统 | |
CN108921014B (zh) | 一种基于改进噪声包络信号识别的螺旋桨轴频搜索方法 | |
CN101236250B (zh) | 声音判定方法和声音判定装置 | |
CN106199532B (zh) | 基于混合傅立叶-小波分析的探地雷达信号降噪方法 | |
CN110108467B (zh) | 基于便携式移动设备的主动发声测速方法 | |
CN115015942B (zh) | 一种自适应水下目标声学测速装置及方法 | |
CN106226078A (zh) | 一种基于麦克风阵列的多普勒畸变声学信号的校正方法 | |
Barbini et al. | Phase editing as a signal pre-processing step for automated bearing fault detection | |
CN113074935B (zh) | 一种用于齿轮箱冲击故障特征的声学分离和诊断方法 | |
CN111856401A (zh) | 一种基于互谱相位拟合的时延估计方法 | |
CN117111048B (zh) | 一种基于改进itd的多波束海底地形测量回波降噪方法 | |
CN104729677B (zh) | 一种非平稳噪声信号的时域数字计权方法 | |
CN113361331A (zh) | 基于加窗插值fft的工频干扰消除方法、系统和介质 | |
CN113899444A (zh) | 一种基于汉宁双窗的振弦传感器共振频率测量方法 | |
CN114578290A (zh) | 一种基于频率修正的麦克风阵列声源定位方法 | |
CN113188797A (zh) | 一种基于传声器阵列的轴承故障诊断方法 | |
CN113093106A (zh) | 一种声源定位方法及系统 | |
Melgoza et al. | Comparing radar receiver pulse deinterleaving performance of differing window functions for bandpass FIR filter design | |
CN107124678B (zh) | 一种音频谐波失真的测量系统 | |
CN116929539A (zh) | 一种基于可视化声源技术的电抗器故障诊断方法与系统 | |
CN111611686A (zh) | 一种通信信号时频域的检测方法 | |
CN116754907A (zh) | 一种变压器故障分析方法和系统 | |
CN111623968A (zh) | 一种基于自适应形态滤波的故障特征提取方法 | |
WO2016006079A1 (ja) | ピーク周波数検出装置、方法およびプログラム |
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 |