基于Gammatone滤波...瞬时频率的水声目标特征提取_连梓旭.pdf
《基于Gammatone滤波...瞬时频率的水声目标特征提取_连梓旭.pdf》由会员分享,可在线阅读,更多相关《基于Gammatone滤波...瞬时频率的水声目标特征提取_连梓旭.pdf(6页珍藏版)》请在咨信网上搜索。
1、第 41 卷第 6 期2022 年 12 月Vol.41,No.6Dec.,2022声学技术Technical Acoustics基于Gammatone滤波器组与子带瞬时频率的水声目标特征提取连梓旭,孙向前(91388部队44分队,广东湛江 524000)摘要:水声目标识别是现代海战中的关键环节,具有重要的军事应用价值。针对海洋环境的复杂性及多变性,文章以模仿人耳听觉感知机理为基础,提出了基于伽玛通(Gammatone)频率瞬时幅频系数的特征提取方法。该方法通过Gammatone滤波器组模拟人耳基底膜对信号的分解、处理过程,在各子带输出信号中提取目标的子带瞬时频率特征,并与传统的瞬时幅度特征相
2、结合,从而更加全面地反映目标的固有属性。利用支持向量机对四类实测水声目标辐射噪声数据进行了识别实验,结果表明,文章的特征提取算法能够较好地抑制噪声,具有较强的鲁棒性。关键词:伽玛通(Gammatone)滤波器组;子带瞬时频率;特征提取;水声目标识别中图分类号:TN911.7 文献标志码:A 文章编号:1000-3630(2022)-06-0827-06Feature extraction of underwater acoustic target based on Gammatone filter-bank and subband instantaneous frequencyLIAN Zix
3、u,SUN Xiangqian(Unit 94,No.91388 Troops of PLA,Zhanjiang524000,Guangdong,China)Abstract:Underwater acoustic target recognition is the key part in modern naval battle and has important military application value.In this paper,by considering the complexity and variability of underwater environment,a f
4、eature extraction method imitating the mechanism of human auditory perception is presented based on Gammatone frequency instantaneous amplitude and frequency coefficients.The Gammatone filter-bank is used to simulate the process of signal decomposition and processing of ear basement membrane,and the
5、 subband instantaneous frequency features of the target extracted from the outputs of Gammatone filter-bank are concatenated with the traditional instantaneous amplitude features to fully reflect the inherent nature of the target.The results of recognition experiments on the measured target data sho
6、w that the proposed feature extraction method can suppress noise and has good robustness.Key words:Gammatone filter-bank;subband instantaneous frequency;feature extraction;underwater acoustic target recognition0引 言水声目标识别是现代海战中的关键技术之一,是水声设备与武器系统智能化的重要标志,具有重要的军事应用价值。在水声目标识别领域,特征提取始终是一个重要的研究课题。然而,海洋环境的
7、复杂性及多变性严重影响水声目标特征提取方法的稳定性,直接导致识别系统性能下降。因此,如何选取鲁棒性好的特征提取方法一直是水声目标识别领域的重要研究内容之一。水声目标辐射噪声的发声机理复杂,其成分主要包含宽带连续谱分量、窄带线谱分量,并且存在明显的幅度调制。其中,宽带连续谱主要由水动力噪声及螺旋桨空化噪声引起,窄带线谱主要由机械噪声引起1。传统的特征提取方法通常是基于谱分析或时频分析理论,但是在实际应用中并未取得令人满意的识别效果。近年来,尽管机器识别取得了巨大的进步,但在大多数情况下,声呐兵的听觉识别效果仍然优于系统的识别效果2。因此,从模拟声呐兵的听觉感知机理出发,研究适用于水声目标的听觉特
8、征提取方法日益受到重视。引用格式:连梓旭,孙向前.基于Gammatone滤波器组与子带瞬时频率的水声目标特征提取J.声学技术,2022,41(6):827-832.LIAN Zixu,SUN Xiangqian.Feature extraction of underwater acoustic target based on Gammatone filter-bank and subband instantaneous frequen-cyJ.Technical Acoustics,2022,41(6):827-832.DOI:10.16300/ki.1000-3630.2022.06.006
9、收稿日期:2021-07-27;修回日期:2021-09-01作者简介:连梓旭(1992),男,黑龙江齐齐哈尔人,硕士,研究方向为水声信号处理。通信作者:连梓旭,E-mail:2022 年声学技术随着听觉心理和生理学的不断发展,听觉特征提取方法取得了巨大的进步,并被广泛应用于语音信号处理中,也有不少学者将其应用于水声目标识别领域。李秀坤等采用伽玛通(Gammatone)滤波器组作为人耳听觉模型并提取水声目标的听觉感知特征2。张揽月等采用美尔频率倒谱系数(Mel Frequency Cepstral Coefficients,MFCC)作为水声目标辐射噪声的特征,证实了听觉特征在水声目标识别中应
10、用的可行性3。马元峰等利用模拟人耳功能的串行压缩伽玛频率滑动(Cascade Compressive Gammachirp,CcGC)滤波器组提取水声目标的听觉谱特征,并分析了听觉特征抑制噪声的机理4。牟林等集成Gammatone听觉滤波器与Meddis模型来模拟耳蜗处理机制,从而提取水声目标的调制谱特征5。吴姚振等结合Gammatone滤波器组与离散小波变换,提取了有效吻合人耳特性的水声目标鲁棒特征6。上述的听觉特征提取方法虽然取得了较好的分类识别效果,但均侧重于水声目标辐射噪声信号的幅度特征,并未深入研究其相位特征。然而,信号的瞬时相位也能反映出目标的部分特性,并且是承载相位信息的主要参数
11、。因此,通过提取水声目标辐射噪声信号的瞬时频率特征,并与幅度特征相结合,能够有效改善水声目标识别系统的性能。本文在已有研究的基础上,基于Gammatone滤波器组与子带瞬时频率,提出Gammatone频率瞬时幅频系数(Gammatone Frequency Instantaneous Amplitude and Frequency Coefficients,GFIAFC)的水声目标特征提取方法。该方法采用Gammatone滤波器组模拟人耳基底膜对信号的频率分解特性,在各子带输出信号中提取目标子带瞬时频率特征,并与传统的瞬时幅度特征相结合,从而更加准确地反映水声目标的固有属性。最后,采用支持向量
12、机(Support Vector Machine,SVM)对四类实测水声目标辐射噪声数据进行分类识别。实验结果表明,该方法在一定信噪比条件下能够有效提升水声目标识别的准确率,具有较强的鲁棒性。1Gammatone滤波器组在人耳的听觉感知系统中,基底膜是进行声信号处理的关键部分。研究表明,基底膜对声信号具有频率分解功能,不同频率的声音会激起基底膜不同部位的振动。从信号处理的角度,可将基底膜的这种频率响应特性抽象为一组带通滤波器,从而模拟基底膜对声信号的频率分解过程。本文用Gammatone 滤波器组来模拟基底膜的频率分解特性。Gammatone滤波器具有简单的时域冲激响应,并且只需要很少的参数就
13、能够较好地模拟基底膜的频率分解功能。Gammatone函数的时域表达式g(t)为g(t)=atn-1exp-2bERB(fc)t cos(2fct+)u(t)(1)式中:u(t)为单位阶跃函数;a为滤波器增益;b为与滤波器带宽有关的参数,通常取1.019;n为滤波器的阶数,研究表明,n=4时Gammatone滤波器就能较好地模拟基底膜的特性7;fc为滤波器的中心频率;为相位,通常设为0;ERB(fc)表示等效矩形带宽(Equivalent Rectangular Bandwidth,ERB),它决定了滤波器的带宽,并且与人耳听觉临界频带有关8。滤波器的ERB与中心频率fc之间的关系可表示为ER
14、B(fc)=24.7+0.108fc(2)Gammatone滤波器组中各滤波器的带宽由式(2)决定,并且每个滤波器的带宽均不同。图1给出了中心频率为1 000 Hz的四阶Gammatone滤波器的时域波形及幅频响应。考虑到计算的复杂性以及水声目标辐射噪声的频率分布范围,本文采用24个滤波器组成24通道的Gammatone滤波器组来模拟人耳基底膜的频率分解功能,其中心频率按照ERB域的转换关系,在505 000 Hz之间均匀分布。频率域与ERB域的转换关系为E=21.4lg(4.37f/1000+1)(3)式中:E为f在ERB域上对应的频率值。各滤波器的中心频率确定后,便可由式(2)计算出各滤波
15、器的ERB9。表1给出了24通道Gammatone滤波器组的中心频率分布情况。图2给出了24通道Gammatone滤波器组的幅频响应。在实际应用中,需要对Gammatone滤波器的时域冲激响应进行离散化处理后才能在计算机中实现。本文采用Slaney提出的方法来实现Gammatone滤波器9。在式(1)中,令w0=2fc,B=2bERB(fc),忽略相位和滤波器增益a,可将式(1)简化为图1四阶Gammatone滤波器组中1 000 Hz中心频率的滤波器时域波形和幅频响应曲线Fig.1Time-domain waveform and amplitude-frequency re-sponse c
16、urves of the filter with a center frequency of 1 000 Hz in the fourth Gammatone filter-bank828第 6 期连梓旭等:基于Gammatone滤波器组与子带瞬时频率的水声目标特征提取g(t)=tn-1exp(-Bt)cos(w0t)u(t)(4)对上式进行拉普拉斯变换,可得Gammatone 函数在s域的传递函数:G(s)=(n-1)!2(s-iw0+B)-n+(s+iw0+B)-n(5)根据s域到z域的映射关系z=esT,并定义a1=cos(w0T),a2=sin(w0T),a3=e-BT,采用冲激响应不
17、变法,可得到4阶Gammatone滤波器在z域上的传递函数:G(z)=T-Ta3a1+(2+1)a2z-11-2a1a3z-1+a23z-2T-Ta3a1-(2+1)a2z-11-2a1a3z-1+a23z-2T-Ta3a1+(2-1)a2z-11-2a1a3z-1+a23z-2T-Ta3a1-(2-1)a2z-11-2a1a3z-1+a23z-2=G1(z)G2(z)G3(z)G4(z)(6)由式(6)可知,在实际应用中,一个 4 阶的Gammatone滤波器通常可由4个2阶的带通滤波器级联而成。2子带瞬时频率瞬时频率是信号瞬时相位的一种表达方式,定义为信号瞬时相位的微分,描述了瞬时相位的变
18、化率,并且已经被广泛应用于信号检测、目标识别等领域的研究中10。然而,基于相位微分的瞬时频率估计方法通常只适用于单频信号。对于宽带信号,由于其存在多个频率分量,并且多分量信号的相位并不具有可加性,因此通过相位微分法可能无法得到正确的瞬时频率。为减轻这一问题带来的影响,本文首先采用24通道Gammatone滤波器组分离出宽带水声信号中的不同频率成分,随后针对各子带输出的窄带信号,计算其子带瞬时频率,从而获得更高的估计精度。对于各子带输出的实信号a(t),首先利用希尔伯特变换构造其解析表达式:x(t)=a(t)+ib(t)(7)b(t)=Ha(t)=1-+a()t-d(8)式中:x(t)为构造的解
19、析信号,b(t)为a(t)的希尔伯特变换。信号的瞬时相位计算公式为(t)=argx(t)=arctanb(t)/a(t)(9)信号的子带瞬时频率可以表示为f=12d(t)dtfs(10)式中:fs为采样频率。然而,直接采用以上方法计算子带瞬时频率依然存在如下问题:(1)计算出的瞬时频率存在负值,但解析信号本身并不存在负频率上的谱,并且进行瞬时频率的估计时有时会出现频率突变的现象;(2)由于离散情况下,需要采用瞬时相位的差分在数值上来近似微分,因此瞬时频率的估计会存在一定的误差11。为解决上述问题,本文采用能量权重滑动平均算法来估计信号的子带瞬时频率12。对于各子带的解析信号x(t),首先将其表
20、示为极坐标形式:x(n)=A(n)expi(n)(11)式中:A(n)为瞬时幅度,(n)为瞬时相位。对于每个采样点,能量权重相位差分向量可定位为y(n)=A2(n)x(n)A(n)x*(n-1)A(n-1)=A2(n)expi(n)-i(n-1)(12)把能量作为权重的优点在于,该式能够有效减表124通道Gammatone滤波器组中心频率分布情况Table 1Centre frequencies of the 24-channel Gammatone filter-bank通道号123456789101112中心频率/HZ50.0086.22127.15173.40225.65284.7035
21、1.41426.79511.96608.20716.94839.80通道号131415161718192021222324中心频率/Hz978.631 135.491 312.741 513.001 739.291 994.972 283.862 610.292 979.123 395.873 866.764 398.82图2 24通道Gammatone滤波器组的幅频响应曲线Fig.2 Amplitude-frequency response curves of the 24-channel Gammatone filter-bank8292022 年声学技术轻频率突变现象,抑制负频率干扰,获
22、得更加精确的瞬时频率估计,且估计的稳健性更好13。然而,当y(n)的实部或虚部在0附近时,其较小的数值变化也可能会导致数量级或数值正负的变化,从而导致其瞬时频率仍然存在突变的现象,并且此突变值还可能导致信号的瞬时频率出现在带限信号的带外。为解决上述问题,将y(n)通过数字低通滤波器,滤波器的输出定义为z(n)。滤波器的滤波算法为z(n)=(1-)y(n)+z(n-1)(13)为滤波系数,表达式为=exp(-fcNfs)(14)式中:fc为滤波器各子带的中心频率,fs为采样频率,N决定了滤波时间,N越大,滤波器输出的子带瞬时频率包络越平滑,但滤波时间相应增长,本文中N取10。由式(14)可以看出
23、,当前滤波输出值主要是由上次滤波输出值决定,本次采样值对当前滤波输出值贡献较小,但具有一定的修正作用14。通过此数字低通滤波器可较好地平滑子带瞬时频率包络,去除瞬时频率突变值,并将子带瞬时频率限制在带限信号的带内。随后,子带瞬时频率可通过式(15)进行估计:f(n)=argz(n)2fs(15)采用能量权重滑动平均法估计出的子带瞬时频率稳健性较好,但对于不同子带,其估计值仍有较大的动态范围,因此将其转换到ERB域上以缩减动态范围,便于后续数据处理。严格意义上讲,从Gammatone滤波器组分离出的各窄带信号可能仍存在多个频率分量,采用相位微分法无法得到正确的瞬时频率。但本文求取的子带瞬时频率承
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Gammatone 滤波 瞬时 频率 目标 特征 提取 连梓旭
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。