一种基于谱熵的啸叫检测方法技术

一种基于谱熵的啸叫检测方法，包括：对待检测信号进行分帧加窗；对分帧加窗后的信号进行频谱分析；划分子带，并计算每个子带的能量；谱熵计算；检测判决，设置谱熵阈值，若当前帧信号的谱熵小于谱熵阈值T0时，则判决为啸叫帧，否则判决为正常信号帧。本发明专利技术的一种基于谱熵的啸叫检测方法，克服了现有啸叫检算法需要动态调整门限并且在不同声场环境下鲁棒性较差的缺点。在不同的声场环境下均有较好的检测效果。实验结果表明，同基于PAPR的啸叫检测算法相比，基于谱熵的检测算法在保证较高检出率的同时，还可明显降低虚警率。

【技术实现步骤摘要】

本专利技术涉及一种啸叫检测方法。特别是涉及一种基于谱熵的啸叫检测方法。
技术介绍
声反馈现象是指在扩声系统中，麦克风拾音将声信号转变为电信号，经功率放大器放大输出后，声音通过声场折回到麦克风后再经功率放大器放大输出，如此反复循环形成正反馈。根据奈奎斯特准则，信号在同时满足相位和增益条件时，会在频点ω0处产生自激震荡：|G(w,t)F(w,t)|≥1∠G(w,t)F(w,t)＝n2π，n为整数G(w,t)为正向路径传递函数，F(w,t)为反馈路径传递函数。当声场的传递函数满足上述相位和增益条件时，将导致扩声系统输出信号的幅值不断地增加,进而产生刺耳的啸叫。通过检测声反馈中出现的啸叫频点，进行陷波处理，降低啸叫频点处增益，破坏啸叫产生的增益条件，从而达到啸叫抑制的目的。啸叫检测是陷波器法的关键，只有及时准确的检测出啸叫成分的频率，才可准确设计相对应中心频率和陷波深度的陷波滤波器。通过级联陷波器进行滤波，抑制啸叫的发生。图1给出啸叫检测过程。由于啸叫本质是单一频率的正弦信号，在频域存在较大的频率分量且频域能量有不断增加的一个过程。正是基于此，相关学者提出一系列相应的啸叫检测算法。主要包括：PAPR(Peak-to-AveragePowerRatio)、PHPR(Peak-to-HarmonicPowerRatio)、PNPR(Peak-to-NeighboringPowerRatio)、IPMP(InterframePeakMagnitudePersistence)、IMSD(InterframeMagnitudeSlopeDeviation)。检测到啸叫成分后，需设计相应的陷波滤波器，降低啸叫频点处增益，破坏啸叫产生的增益条件，达到抑制啸叫的目的。最常用的陷波滤波器是二阶IIR滤波器，因为IIR滤波器可以用较少阶数获得较好的选择特性，所用存储单元少，运算次数少，较为经济而且高效。2阶IIR滤波器的系统函数：H(z)=b0z-1+b1z-1+b2z-11+a1z-1+a2z-1.]]>目前已有的啸叫检测算法需要动态调整门限并且在不同声场环境下鲁棒性较差。
技术实现思路
本专利技术所要解决的技术问题是，提供一种在不同的声场环境下均有较好的检测效果的基于谱熵的啸叫检测方法。本专利技术所采用的技术方案是：一种基于谱熵的啸叫检测方法，包括如下步骤：1)对待检测信号进行分帧加窗；2)对分帧加窗后的信号进行频谱分析；3)划分子带，并计算每个子带的能量；4)谱熵计算根据步骤3)计算的子带能量，相应的概率密度函数和谱熵为Px(i,m)=Sx(i,m)/Σj=0M-1Sx(i,j)]]>m＝0,1,2,…，M-1Hx(i)=Σm=0M-1Px(i,m)log1Px(i,m)]]>m＝0,1,2,…，M-1其中，Px(i,m)表示第i帧信号的第m个子带的概率密度函数，Hx(i)表示第i帧信号的谱熵；5)检测判决设置谱熵阈值T0，若当前帧信号的谱熵小于谱熵阈值T0时，则判决为啸叫帧，否则判决为正常信号帧。步骤1)所述的加窗是加窗函数。步骤2)所述的频谱分析是利用FFT分析计算得到能量谱：X(i,k)=Σn=0N-1x(n)w(n)e-j2πnkN]]>k＝0，1,2,…,N-1Rx(i,k)＝|X(i,k)|2k＝0，1,2,…,N-1x(n)为待检测信号，w(n)为所加窗函数，N为进行FFT的数据长度，e为自然底数，j表示虚数，X(i,k)为第i帧信号的第k个频点的频谱，Rx(i,k)为第i帧信号的第k个频点的能量谱。步骤3)是根据步骤2)得到的能量谱，将整个频带分为若干个子带，再分别计算每个子带的能量：Sx(i,m)=Σk∈BmRx(i,k)]]>m＝0,1,2,…，M-1Sx(i,m)表示在第i帧信号的第m个子带的能量，M表示划分子带的个数，Bm表示第m个子带对应的所有频点。本专利技术的一种基于谱熵的啸叫检测方法，克服了现有啸叫检算法需要动态调整门限并且在不同声场环境下鲁棒性较差的缺点。在不同的声场环境下均有较好的检测效果。实验结果表明，同基于PAPR的啸叫检测算法相比，基于谱熵的检测算法在保证较高检出率的同时，还可明显降低虚警率。附图说明图1是啸叫检测过程示意图；图2a是啸叫信号下基于PAPR的检测效果图；图2b是啸叫信号下每帧的PAPR以及PAPR阈值示意图；图2c是啸叫信号下基于谱熵的检测效果图；图2d是啸叫信号下每帧的谱熵以及谱熵阈值示意图；图3a是正常音乐信号下基于PAPR的检测效果图；图3b是正常音乐信号下每帧的PAPR以及PAPR阈值示意图；图3c是正常音乐信号下基于谱熵的检测效果图；图3d是正常音乐信号下每帧的谱熵以及谱熵阈值示意图；图4是本专利技术的一种基于谱熵的啸叫检测方法的流程图。具体实施方式下面结合实施例和附图对本专利技术的一种基于谱熵的啸叫检测方法做出详细说明。本专利技术的一种基于谱熵的啸叫检测方法，包括如下步骤：1)对待检测信号进行分帧加窗，所述的加窗是加窗函数(例如汉宁窗、海明窗、矩形窗等)。时域离散信号x(n)是无限长的，用FFT做能量谱估计，必须把x(n)限制在一定的时间区域之内，即进行数据截断。数据的截断相当于加窗分帧处理，每帧的长度FrameLen，每帧偏移量ShiftLen。同时，为了尽量减小频谱泄露和谱间干扰的影响，不让数据突然截断，应加缓变的窗(例如汉宁窗、海明窗等)，使得加窗后的谱旁瓣能量更小，卷积后造成的泄露更小。2)对分帧加窗后的信号进行频谱分析，所述的频谱分析是利用FFT分析计算得到能量谱：X(i,k)=Σn=0N-1x(n)w(n)e-j2πnkN]]>k＝0，1,2,…,N-1Rx(i,k)＝|X(i,k)|2k＝0，1,2,…,N-1x(n)为待检测信号，w(n)为所加窗函数，N为进行FFT的数据长度，e为自然底数，j表示虚数，X(i,k)为第i帧信号的第k个频点的频谱，Rx(i,k)为第i帧信号的第k个频点的能量谱。3)划分子带，并计算每个子带的能量，是根据步骤2)得到的能量谱，将整个频带分为若干个子带，再分别计算每个子带的能量：Sx(i,m)=Σk∈BmRx(i,k)]]>m＝0,1,2,…，M-1Sx(i,m)表示在第i帧信号的第m个子带的能量，本文档来自技高网...

【技术保护点】
一种基于谱熵的啸叫检测方法，其特征在于，包括如下步骤：1)对待检测信号进行分帧加窗；2)对分帧加窗后的信号进行频谱分析；3)划分子带，并计算每个子带的能量；4)谱熵计算根据步骤3)计算的子带能量，相应的概率密度函数和谱熵为Px(i,m)=Sx(i,m)/Σj=0M-1Sx(i,j)]]> m＝0,1,2,…，M‑1Hx(i)=Σm=0M-1Px(i,m)log1Px(i,m)]]> m＝0,1,2,…，M‑1其中，Px(i,m)表示第i帧信号的第m个子带的概率密度函数，Hx(i)表示第i帧信号的谱熵；5)检测判决设置谱熵阈值T0，若当前帧信号的谱熵小于谱熵阈值T0时，则判决为啸叫帧，否则判决为正常信号帧。

【技术特征摘要】
1.一种基于谱熵的啸叫检测方法，其特征在于，包括如下步骤：
1)对待检测信号进行分帧加窗；
2)对分帧加窗后的信号进行频谱分析；
3)划分子带，并计算每个子带的能量；
4)谱熵计算
根据步骤3)计算的子带能量，相应的概率密度函数和谱熵为
Px(i,m)=Sx(i,m)/Σj=0M-1Sx(i,j)]]>m＝0,1,2,…，M-1
Hx(i)=Σm=0M-1Px(i,m)log1Px(i,m)]]>m＝0,1,2,…，M-1
其中，Px(i,m)表示第i帧信号的第m个子带的概率密度函数，Hx(i)表示第i帧信号的谱
熵；
5)检测判决
设置谱熵阈值T0，若当前帧信号的谱熵小于谱熵阈值T0时，则判决为啸叫帧，否则判决
为正常信号帧。
2.根据权利要求1所述的一种基于谱熵的啸叫检测方法，其特征在于，步骤1)所述的加
窗是加窗函数。
3.根据权利要求1所述的一种基于谱熵的...

【专利技术属性】
技术研发人员：张涛，李康康，张彩霞，陈有志，
申请(专利权)人：天津大学，
类型：发明
国别省市：天津;12