CN109085595A - 一种使用水听器接收信号估计空中运动声源速度的方法 - Google Patents
一种使用水听器接收信号估计空中运动声源速度的方法 Download PDFInfo
- Publication number
- CN109085595A CN109085595A CN201810822777.4A CN201810822777A CN109085595A CN 109085595 A CN109085595 A CN 109085595A CN 201810822777 A CN201810822777 A CN 201810822777A CN 109085595 A CN109085595 A CN 109085595A
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- sound source
- time
- hydrophone
- 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 34
- 238000001228 spectrum Methods 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 2
- 230000001052 transient effect Effects 0.000 abstract 1
- 230000005855 radiation Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012897 Levenberg–Marquardt algorithm Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 101150061215 outM gene Proteins 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种使用水听器接收信号估计空中运动声源速度的方法,将PCT与Ferguson和Lo提出的使用单水听器接收信号估计空中目标速度的方法结合在一起。首先,利用PCT估计水听器接收信号在整个时间序列上的瞬时频率,再从瞬时频率中计算4个参数即声源频率f0、运动速度v、飞行高度h和CPA时间μc的初始估计值,最后利用Ferguson和Lo建立的瞬时频率参数模型在最小二乘法准则下迭代拟合瞬时频率估计值求解目标速度。有益效果:估计出了较为准确的信号瞬时频率。并以此为基础利用瞬时频率和声源参数的关系实现了使用水听器接收信号对空中运动声源速度的精确估计。上述方法比利用STFT得到瞬时频率的参数估计方法更为准确。
Description
技术领域
本发明属于信号处理、时频分析理论和水声测量等领域涉及一种使用水听器接收信号估计空中运动声源速度的方法,利用时频分析方法以及水听器接收信号瞬时频率 与空中声源参数之间的关系实现使用水听器接收信号对空中运动声源的速度估计。
背景技术
估计空中运动声源的速度是水下对空预警技术的一个重要方面。很大一部分飞机是由螺旋桨驱动的,其辐射噪声主要分为宽带连续谱噪声和线谱噪声两部分,线谱成 分呈谐波关系,主要集中在低频段,强度相比于宽带连续谱噪声较高,跨空气-水界面 远距离传播后更容易被探测。飞机在巡航时基本可以认为是保持匀速直线运动的,当 水面上方的飞机相对于水听器运动时,线谱噪声会发生多普勒频移现象。水听器接收 信号的瞬时频率中包含着声源频率、运动速度、飞行高度以及声源经过与水听器的最 近会遇点(ClosestPoint of Approach,CPA,如图1所示)时间4个空中声源参数的信 息。因此,利用水听器接收到声源辐射噪声的瞬时频率就可以对空中声源的速度进行 估计。
Ferguson和Lo利用线谱辐射噪声的多普勒效应,建立了瞬时频率与速度等参数的模型,提出了一种使用单个水听器的接收信号估计空中目标速度的方法,并通过实验 评估了该方法的性能。这种方法首先利用短时傅里叶变换(Short-Time Fourier Transform,STFT)从水听器接收信号中估计出其瞬时频率,然后利用瞬时频率估计空 中声源的飞行速度(B.G.Ferguson,K.W.Lo.Transiting aircraft parameter estimation usingunderwater acoustic sensor data[J].IEEE Journal of Oceanic Engineering,1999,24(4): 424-435.)。然而,在水下环境噪声的干扰下,STFT的时频分辨率和能量聚集度较低, 估计瞬时频率的精度下降从而影响上述速度估计方法的性能。
为了提高瞬时频率的估计精度,需要选取一种能量聚集度高且对噪声有抑制作用的时频分析方法,以准确地估计水听器接收信号的瞬时频率。Z.K.Peng和G.Meng 等人在传统Chirplet变换的基础上提出了多项式Chirplet变换(Polynomial Chirplet Transform,PCT),PCT具有能量聚集度高的优点,可以有效地抑制噪声的干扰,并且 适用于非线性调频信号(Z.K.Peng,G.Meng,et al.Polynomial chirplet transform with application toinstantaneous frequency estimation[J].IEEE Journal of Oceanic Engineering,2011,60(9):3222-3229.)。而空中运动声源在水下产生的多普勒频移信号正是一种典型的非线性调频信号,所以PCT非常适合用来估计的水听器接收信号的瞬时频率。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种使用水听器接收信号估计空中运动声源速度的方法。
技术方案
一种使用水听器接收信号估计空中运动声源速度的方法,其特征在于步骤如下:
步骤1:放大水听器接收到的信号并记录,记为sr(t),声源为亚音速飞机,并且 声波在跨介质入水的过程中入射角小于临界角12.8°;声源某一特定线谱的原始频率 恒为f0,运动速度为v,高度为h,CPA时间为μc;
步骤2:对记录到的信号sr(t)利用STFT进行预处理,得到接收信号的时频分布 后采用带通滤波,对滤波后的信号进行希尔伯特变换得到相应的复信号z(t);
步骤3:利用PCT对预处理后的信号z(t)进行时频分析,得到预处理后信号z(t)的时频分布:
PCT的定义式
其中:τ表示对信号做PCT的时刻,f表示信号频率,wα(t)表示非负、对称且归 一化实窗;
采用高斯函数,α给定了高斯窗的时间跨度,和分别为频率旋转算子和频移算子;
所述β=(β1,β1,,βN)为多项式系数向量;
通过迭代实现以β为系数的多项式对信号瞬时频率的最优逼近,收敛后得到最优多项式系数,相应的PCT结果将在下一步被用来估计信号z(t)的瞬时频率;
步骤4:在时刻τk对信号能量在频率上的分布进行三点二次插值,能量峰值对应的频率即为信号的瞬时频率估计值以此获得信号在采样时间上的瞬时频率估计序 列其中,τk表示时间序列中的第k个时刻,K表示时间序列的长度;
步骤5:以水下观测空中运动线谱声源的瞬时频率模型在非线性最小二乘准则下拟合逼近上一步得到的瞬时频率估计序列通过迭代算法估计声源参数, 当算法收敛后即得到声源速度的估计值
有益效果
本发明提出的一种使用水听器接收信号估计空中运动声源速度的方法,将PCT与Ferguson和Lo提出的使用单水听器接收信号估计空中目标速度的方法结合在一起。首 先,利用PCT估计水听器接收信号在整个时间序列上的瞬时频率,再从瞬时频率中计 算4个参数即声源频率f0、运动速度v、飞行高度h和CPA时间μc的初始估计值,最 后利用Ferguson和Lo建立的瞬时频率参数模型在最小二乘法准则下迭代拟合瞬时频 率估计值求解目标速度。
本发明的有益效果是:利用适用于非线性调频信号的PCT对多普勒频移信号进行处理,估计出了较为准确的信号瞬时频率。并以此为基础利用瞬时频率和声源参数的 关系实现了使用水听器接收信号对空中运动声源速度的精确估计。上述方法比利用 STFT得到瞬时频率的参数估计方法更为准确。
附图说明
图1是声音从声源跨空气-水介质到达水听器的传播过程示意图;
图2是使用水听器接收信号估计空中运动声源速度的总体流程图;
图3是PCT中多项式系数优化算法的流程图。
图4是利用PCT得到的信号z(t)的时频分布
具体实施方式
现结合实施例、附图对本发明作进一步描述:
本发明解决其技术问题所采用的技术方案总体流程如图2所示,具体包括以下步骤:
1)放大水听器接收到的信号并记录,记为sr(t)。
2)对记录到的信号sr(t)进行预处理,得到z(t)。
3)利用PCT得到预处理后信号z(t)的时频分布。
4)对信号的时频分布做平滑处理后提取瞬时频率。
5)利用非线性最小二乘法得到声源速度的估计值
以下对本发明的每一个步骤做详细说明:
所述步骤1)的具体实现如下:
声音从声源跨空气-水介质到达水听器的传播过程如图1所示,图中给出了空中声源的飞行轨迹和水听器的空间位置。水听器布放深度为d,目标在巡航状态下以恒定 的高度、方向和速度从水听器上方飞过,经过水听器正上方(CPA)。在这个过程中, 水听器持续接收空中声源辐射的声信号,并将其转换成电信号,电信号通过功率放大 器后用数据采集仪器记录下来,记为sr(t)。其中0<t≤T,T表示数据记录的时间长 度。本发明适用的前提是声源为亚音速飞机,并且声波在跨介质入水的过程中入射角 小于临界角(12.8°)。记声源某一特定线谱的原始频率恒为f0,运动速度为v,高度 为h,CPA时间为μc。
所述步骤2)的具体实现如下:
首先利用STFT粗略得到接收信号的时频分布,从中可以大致确定出声源辐射噪声某频率线谱多普勒频移的范围,采用带通滤波器滤除该频率范围以外的干扰,对滤 波后的信号进行希尔伯特变换得到相应的复信号z(t)。
所述步骤3)的具体实现如下:
利用PCT对预处理后的信号z(t)进行时频分析。PCT的定义为
其中τ表示对信号做PCT的时刻,f表示信号频率,wα(t)表示非负、对称且归一化实窗,通常采用高斯函数,有
wα(t)=exp[(logα)t2] (2)
α给定了高斯窗的时间跨度,β=(β1,β1,,βN)为多项式系数向量,N为多项式阶数,和分别为频率旋转算子和频移算子,表达式为
为了得到最优的多项式系数β,通过迭代用多项式拟合逼近信号的瞬时频率,其流程如图3所示:第一步设定多项式系数初值为β0=(0,0,...,0)、并给定多项式阶数N和 判决门限δ;第二步将多项式系数β代入式(1)中得到z(t)关于(τ,f)的时频分布,从中 提取峰值所对应的瞬时频率随τ的变化轨迹IF(τ);第三步用多项式在最小 二乘准则下拟合逼近前一步得到的瞬时频率,获得多项式系数βi;第四步计算
其中τk表示时间序列中的第k个时刻,K表 示时间序列的长度,i表示当前迭代的次数,i-1表示上一次迭代的次数,若ε≤δ则结 束迭代,否则重复第二步至第三步直到ε≤δ;迭代结束后得到的多项式系数β即为多 项式系数。
所述步骤4)的具体实现如下:
将最优多项式系数β带入到式(1)中得到信号z(t)的时频分布,为了提高频率分 辨率、降低噪声的干扰,对其做平滑处理后提取瞬时频率。在给定时刻τk对信号能量 在频率上的分布PCT(τk,f)做关于频率的三点二次插值,得到更加光滑的分布 PCT′(τk,f)。能量峰值对应的频率为信号的瞬时频率估计值,即 重复这一操作,以此获得信号在采样时间上的瞬时频率估计 序列
所述步骤5)的具体实现如下:
用水下观测空中运动线谱声源的瞬时频率模型在非线性最小二乘准则下拟合逼近 上一步得到的瞬时频率估计序列,估计声源参数,即求解以下最小化问题:
其中,为τk时刻水听器瞬时频率的估计值,f′0,v′,h′,μ′c为待定的声源参数,
fk(f′0,v′,h′,μ′c)为τk时刻水听器接收信号瞬时频率的理论值,满足以下模型:
其中
γ=ca/cw,ca和cw分别为空气中和水中的声速,x(μ)=v′(μ-μc),lw(μ)为方程x(μ)2+d2(1+γh′/lw(μ))2=(lw(μ)+γh′)2的正实根。
公式(4)给出的是一个非线性最小二乘优化问题,只能通过迭代求解,为了提高迭代的速度,首先选取几个特殊时刻,粗略估计4个声源参数作为迭代的初始值声源频率的初始估计值为
飞行速度的初始估计值为
其中,通过计算瞬时频率序列前M个值的均值和最后M个值的均值分别可以得到和即
飞行高度的估计值为
其中,是τk时刻瞬时频率估计值的变化率,k为某一时刻在时间序列中排列的位 置。可以利用有限差分由f(τ)近似得到。最小时τk的取值,即
CPA时间的初始估计值为
由于噪声的干扰,瞬时频率估计序列上存在扰动,因此仅利用几个时刻的瞬时频率估计值得到的参数初始估计值的精度不高。为了充分利用完整的瞬时频 率序列信息,把参数初始估计值代入式(4)中,利用通用的非线性最小二乘 方法(如Matlab中的lsqnonlin函数采用的Levenberg-Marquardt算法)求解该问题, 即可得到声源速度的精确估计值
具体实施例子
为了验证利用本发明给出的方法估计空中运动声源参数的有效性,设计仿真实验如下:假设空气中声速ca=340m/s,水中声速cw=1524m/s,水听器到海平面的垂直距 离d=50m,空中声源的线谱辐射噪声频率f0=68Hz,做速度v=150m/s的匀速直线运 动、飞行高度h=300m、CPA时间μc=0s,采样率fs=512Hz,观察时间t=[-1,3]s, 在这段时间内,声波在跨空气-水界面传播时的入射角小于临界角。PCT中多项式阶数N取为20,判决门限δ取为0.001,高斯窗长度取为1024,α取为0.005,计算时 M取为25。由于空中声源是运动状态,它与水听器之间的距离是变化的,所以水听器 接收信号的幅值也随之改变,以致信噪比也是不断改变的。因此假设水听器接收到的 背景噪声大小不变,考虑空中声源辐射噪声的强度按球面扩展规律衰减,在观察时间 两端处信噪比SNR=0dB,噪声为高斯白噪声。根据以上参数得到仿真的水听器接收信 号,并使用本发明给出的方法对仿真信号进行处理,得到水听器接收信号的时频分布 如图4,空中声源参数估计值如表1。
表1空中运动声源参数估计值
从表1中可以看出利用本发明给出的方法可以较为准确地估计出空中运动声源的速度等参数,证明了这种方法的有效性。
Claims (1)
1.一种使用水听器接收信号估计空中运动声源速度的方法,其特征在于步骤如下:
步骤1:放大水听器接收到的信号并记录,记为sr(t),声源为亚音速飞机,并且声波在跨介质入水的过程中入射角小于临界角12.8°;声源某一特定线谱的原始频率恒为f0,运动速度为v,高度为h,CPA时间为μc;
步骤2:对记录到的信号sr(t)利用STFT进行预处理,得到接收信号的时频分布后采用带通滤波,对滤波后的信号进行希尔伯特变换得到相应的复信号z(t);
步骤3:利用PCT对预处理后的信号z(t)进行时频分析,得到预处理后信号z(t)的时频分布:
PCT的定义式
其中:τ表示对信号做PCT的时刻,f表示信号频率,wα(t)表示非负、对称且归一化实窗;
采用高斯函数,α给定了高斯窗的时间跨度,和分别为频率旋转算子和频移算子;
所述β=(β1,β1,...,βN)为多项式系数向量;
通过迭代实现以β为系数的多项式对信号瞬时频率的最优逼近,收敛后得到最优多项式系数,相应的PCT结果将在下一步被用来估计信号z(t)的瞬时频率;
步骤4:在时刻τk对信号能量在频率上的分布进行三点二次插值,能量峰值对应的频率即为信号的瞬时频率估计值以此获得信号在采样时间上的瞬时频率估计序列其中,τk表示时间序列中的第k个时刻,K表示时间序列的长度;
步骤5:以水下观测空中运动线谱声源的瞬时频率模型在非线性最小二乘准则下拟合逼近上一步得到的瞬时频率估计序列通过迭代算法估计声源参数,当算法收敛后即得到声源速度的估计值
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810822777.4A CN109085595B (zh) | 2018-07-25 | 2018-07-25 | 一种使用水听器接收信号估计空中运动声源速度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810822777.4A CN109085595B (zh) | 2018-07-25 | 2018-07-25 | 一种使用水听器接收信号估计空中运动声源速度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109085595A true CN109085595A (zh) | 2018-12-25 |
CN109085595B CN109085595B (zh) | 2022-06-07 |
Family
ID=64838342
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810822777.4A Active CN109085595B (zh) | 2018-07-25 | 2018-07-25 | 一种使用水听器接收信号估计空中运动声源速度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109085595B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110488219A (zh) * | 2019-08-26 | 2019-11-22 | 自然资源部第一海洋研究所 | 一种浅海单水听器运动目标距离估计方法 |
CN112114286A (zh) * | 2020-06-23 | 2020-12-22 | 山东省科学院海洋仪器仪表研究所 | 基于线谱生命周期和单矢量水听器的多目标跟踪方法 |
CN114580479A (zh) * | 2022-03-04 | 2022-06-03 | 安徽工业大学 | 一种基于stft周期性能量峰的船舶螺旋桨转速预报方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3716822A (en) * | 1969-09-25 | 1973-02-13 | Y Tsuruta | Method and device for detecting relative speed by using supersonic wave |
CN1588124A (zh) * | 2004-08-20 | 2005-03-02 | 清华大学 | 利用多普勒变换的被动测速测距法及装置 |
CN101900811A (zh) * | 2010-07-23 | 2010-12-01 | 哈尔滨工程大学 | 一种基于单水听器的被动测距方法 |
CN103020479A (zh) * | 2012-12-28 | 2013-04-03 | 上海交通大学 | 基于非线性调频小波变换的信号瞬时频率估计方法 |
CN103076604A (zh) * | 2012-12-31 | 2013-05-01 | 东南大学 | 一种基于频散特征的低频水声脉冲信号距离的测量方法 |
CN105425212A (zh) * | 2015-11-18 | 2016-03-23 | 西北工业大学 | 一种声源定位方法 |
US20160139245A1 (en) * | 2014-05-22 | 2016-05-19 | The United States Of America As Represented By The Secretary Of The Navy | Sparsity-Driven Passive Tracking of Acoustic Sources |
CN107515390A (zh) * | 2017-09-15 | 2017-12-26 | 哈尔滨工程大学 | 一种基于单矢量传感器的空中目标定位方法 |
WO2018047805A1 (ja) * | 2016-09-09 | 2018-03-15 | 日本電気株式会社 | 移動音源速度推定装置、速度監視システム、移動音源速度推定方法、および移動音源速度推定用プログラムが記憶された記憶媒体 |
-
2018
- 2018-07-25 CN CN201810822777.4A patent/CN109085595B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3716822A (en) * | 1969-09-25 | 1973-02-13 | Y Tsuruta | Method and device for detecting relative speed by using supersonic wave |
CN1588124A (zh) * | 2004-08-20 | 2005-03-02 | 清华大学 | 利用多普勒变换的被动测速测距法及装置 |
CN101900811A (zh) * | 2010-07-23 | 2010-12-01 | 哈尔滨工程大学 | 一种基于单水听器的被动测距方法 |
CN103020479A (zh) * | 2012-12-28 | 2013-04-03 | 上海交通大学 | 基于非线性调频小波变换的信号瞬时频率估计方法 |
CN103076604A (zh) * | 2012-12-31 | 2013-05-01 | 东南大学 | 一种基于频散特征的低频水声脉冲信号距离的测量方法 |
US20160139245A1 (en) * | 2014-05-22 | 2016-05-19 | The United States Of America As Represented By The Secretary Of The Navy | Sparsity-Driven Passive Tracking of Acoustic Sources |
CN105425212A (zh) * | 2015-11-18 | 2016-03-23 | 西北工业大学 | 一种声源定位方法 |
WO2018047805A1 (ja) * | 2016-09-09 | 2018-03-15 | 日本電気株式会社 | 移動音源速度推定装置、速度監視システム、移動音源速度推定方法、および移動音源速度推定用プログラムが記憶された記憶媒体 |
CN107515390A (zh) * | 2017-09-15 | 2017-12-26 | 哈尔滨工程大学 | 一种基于单矢量传感器的空中目标定位方法 |
Non-Patent Citations (6)
Title |
---|
BRIAN G. FERGUSON等: "Transiting Aircraft Parameter Estimation Using Underwater Acoustic Sensor Data", 《IEEE JOURNAL OF OCEANIC ENGINEERING》 * |
LINGJI XU等: "Analysis of moving source characteristics using polynomial chirplet transform", 《THE JOURNAL OF ACOUSTICAL SOCIETY OF AMERICA》 * |
LINGJI XU等: "Parameter estimation of underwater moving sources by using matched Wigner transform", 《APPLIED ACOUSTICS》 * |
Z. K. PENG等: "Polynomial Chirplet Transform With Application to Instantaneous Frequency Estimation", 《IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT》 * |
徐灵基: "水声信号多普勒参数估计及其应用研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 * |
杨龙等: "水下运动声源多普勒频率估计的PCT时频分析方法", 《声学技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110488219A (zh) * | 2019-08-26 | 2019-11-22 | 自然资源部第一海洋研究所 | 一种浅海单水听器运动目标距离估计方法 |
CN112114286A (zh) * | 2020-06-23 | 2020-12-22 | 山东省科学院海洋仪器仪表研究所 | 基于线谱生命周期和单矢量水听器的多目标跟踪方法 |
CN114580479A (zh) * | 2022-03-04 | 2022-06-03 | 安徽工业大学 | 一种基于stft周期性能量峰的船舶螺旋桨转速预报方法 |
CN114580479B (zh) * | 2022-03-04 | 2024-09-06 | 安徽工业大学 | 一种基于stft周期性能量峰的船舶螺旋桨转速预报方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109085595B (zh) | 2022-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109085595B (zh) | 一种使用水听器接收信号估计空中运动声源速度的方法 | |
CN108226933B (zh) | 一种基于条纹干涉结构的深海宽带目标深度估计方法 | |
CN106970371A (zh) | 一种基于Keystone和匹配滤波的目标检测方法 | |
CN105589066B (zh) | 一种利用垂直矢量阵估计水下匀速运动航行器参数的方法 | |
Ferguson et al. | Turboprop and rotary-wing aircraft flight parameter estimation using both narrow-band and broadband passive acoustic signal-processing methods | |
Ferguson | Time-delay estimation techniques applied to the acoustic detection of jet aircraft transits | |
CN107515390B (zh) | 一种基于单矢量传感器的空中目标定位方法 | |
CN113281727B (zh) | 一种基于水平线列阵的输出增强的波束形成方法及其系统 | |
CN113687308B (zh) | 一种基于弯曲波的冰上震源定位方法 | |
CN106772248B (zh) | 高速列车模型在声学风洞内的噪声源定位方法 | |
CN111257943A (zh) | 一种浅海低频声场时频干涉条纹提取方法 | |
CN108646248B (zh) | 一种针对低速运动声源的被动声学测速测距方法 | |
CN108562891A (zh) | 一种深海低信噪比条件下声源深度自主实时跟踪方法 | |
CN110515077A (zh) | 一种高轨舰船目标子孔径isar成像方法 | |
CN110716202B (zh) | 基于重采样的调频连续波声呐测距方法及装置 | |
CN109116359B (zh) | 一种机载雷达低空风切变风场回波风速的估计方法 | |
CN112671489B (zh) | 基于沃特森模型的短波航空移动信道建模方法 | |
Timlelt et al. | Closed-form solution to motion parameter estimation of an acoustic source exploiting Doppler effect | |
Ferguson | Time-frequency signal analysis of hydrophone data | |
CN116794643A (zh) | 一种基于声反射断层扫描成像的逆投影方法 | |
CN115242583B (zh) | 一种基于水平线列阵的信道脉冲响应被动估计方法 | |
CN115902849A (zh) | 一种基于波束输出强度重采样的深海声源深度估计方法 | |
Hutcheson et al. | Experimental Study of Noise Shielding by a NACA 0012 Airfoil | |
Kroeber et al. | Comparison of microphone array measurements of an airfoil with high-lift devices in open and closed wind tunnels | |
Duarte et al. | Ultrasonic time of flight estimation for wind speed measurement based on time-frequency domain using STFT |
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 |