《实验二啸叫检测与抑制系统》由会员分享，可在线阅读，更多相关《实验二啸叫检测与抑制系统（19页珍藏版）》请在金锄头文库上搜索。
1、实验二啸叫检测与抑制系统一、实验目的1、认识扩音系统中的啸叫现象；2、了解啸叫产生的条件；3、分析啸叫信号的频谱；4、掌握啸叫信号检测的方法；5、掌握啸叫抑制的原理；6、实现啸叫抑制。二、实验设备1、音频功率放大系统；2、录音机；3、计算机；4、matlab 软件三、实验内容1、认识啸叫及其危害声反馈现象在日常生活中非常常见。在多动能报告厅，KTV等同时出现扬声 器和麦克风的场合，由于扬声器和麦克风之间存在电声耦合，必然会导致声反馈 现象的产生。声反馈会在反馈回路中产生再生混响，使讲话、唱歌的声音严重失 真，音质受到破坏，清晰度大大降低。严重时甚至会产生自激啸叫，限制扩声 系统传声增益的提升，
2、使整个系统的正常工作受到影响。还可能会烧毁系统放大 器、扬声器中的高音单元，甚至会对人的听力造成损伤。所以，在应用到扩声系 统的场合，啸叫的检测和抑制非常重要。2、啸叫信号产生的原理声音信号首先从麦克风拾入，经过扩声系统的功率放大器放大后由扬声器送 出，经过各种障碍物的多次反射后，又被麦克风拾入，从而形成一个闭合环路。 如果传声器对某些频点的拾音灵敏度过高，导致声音在这些频点的增益是正值， 就形成了一个正反馈过程，声音信号经过多次反复循环放大后，在某些频点的声音强度超过一定的增益上限，就会发生自激振荡，从而产生啸叫。啸叫的产生必须同时具备以下三个条件:(1) 扬声器和拾音设备(麦克风)要处于同
3、一声场中，从而保证扬声器输出的信号能被拾音设备再度采集，形成正反馈;(2) 拾音设备的拾音灵敏度高，系统的传声增益大；(3) 声场存在缺陷共振，即扩声系统的频谱特性不平坦，在某些频点上容 易出现共振。3、啸叫信号的特征为了检测和抑制啸叫信号，需要对啸叫信号的基本特性做一定的分析，对啸 叫信号在时域和频域上进行分析，最终得出以下结论：啸叫信号的时域波形是一个频率恒定的正弦波，其幅值随着时间的增加迅速 增大，直到超出了功放放大区，进入饱和区和截止区时，产生削波现象，如图2 所示。啸叫信号的频谱图中存在单一的，并且固定不变的啸叫频率点，并且啸叫频 率点对应的幅值远远大于语音信号中其他频率点的幅值，如
4、图2所示。图2啸叫时时域波形图3啸叫时频域波形该实验要做的主要工作就是在声反馈建立的初期，在尽可能少的时间内检测 出啸叫对应的的频率点，并设计相应的陷波器将其抑制掉，保证系统工作在正常 的线性放大区，预防谐波信号的产生。4、啸叫抑制的方法及原理对啸叫的抑制方法大体上分为被动型预防和主动型控制两大类。被动型预防 是根据声场特性，从扩声系统的系统设计、设备选型、声场布局以及声场调整等 环节着手做好预防啸叫的工作；主动型控制力求在控制啸叫的同时提高扩声系统 的传声增益。主要有以下几种方法：基于建筑声学特性在音响界，最初解决回声啸叫的办法是通过降低扩声系统的增益来确保系统 的稳定工作，但是，这样就无法
5、完全满足扩声音量的需要。所以，工程技术人员 开始在室内建筑声学、结构声学以及室内装饰装修方面下功夫，来解决声反馈问 题和啸叫现象。但是，这种方法只是被动地抑制啸叫，造价及施工难度太大，不 具有普遍推广意义。(2)均衡法图4均衡器均衡法的原理就是通过衰减或切除声音信号中某些过强的频率来抑制住声 反馈。均衡器首先将声音频率进行分段，然后通过调节各频率段增益来弥补声场 合扬声器的缺陷，达到抑制某一频点啸叫的目的。但是，这个方法需专业音响师 手动调节，要求音响师准确判断啸叫频点并迅速做出反应；其次，在实际操作中， 由于均衡器对啸叫频点的定位不够精确，会严重损伤音质。移频法移频就是移动频率，移频器的工作
6、原理是对声音信号的频点向上或向下 移几个或几十个频点，破坏构成声反馈的条件，最终达到防止啸叫的目的。但是, 移频法导致了整个声音频率范围内的频率失真，对音质损伤较大；而且，移频器 的工作状态与现场有无啸叫点没有关系，这样就会造成新的失真。因此，应该尽 量少用这种方法。压限器它的工作原理是根据输入声音信号的大小自动改变输出信号的增益。当 声音信号强度超过设定的电平值时，压缩器就将压缩比调到-:1,语音信号强 度就不会再继续增大，啸叫也就自然而然地被消除掉了。但是，这种抑制方法会 对声音带来动态损失，不适宜要求较高的场合。(5)窄带陷波器利用频谱分析和判定条件找出声音信号的啸叫频率点，设计参数可调
7、的数字 陷波器对声反馈频点进行抑制。由于窄带陷波器在某个很窄的频带内有很大的衰 减，因此对啸叫频点定位精度高，陷波器带宽窄，可以较为准确地抑制啸叫频点, 而且对音质的损伤比较低，适合于对音质要求较高的场合。5啸叫抑制系统基本性能指标设计的啸叫抑制系统要具备以下条件：第一，陷波器的阻带带宽要越窄越好，这样可以减少声音的动态失真。第二，对啸叫频点进行抑制时，衰减量的陷波深度越深越好，保证对啸叫的 完全抑制；第三，啸叫频率出现在任何频点，系统都可以自动检测到该频点，并对其进 行抑制；该实验的性能指标如下：处理字长：16bit二阶IIR数字陷波器陷波滤波器频率范围：100Hz-5KHz陷波滤波器带宽：
8、1/10倍频程频率定位精度：1Hz响应时间：1s6、系统算法及其实现61啸叫的检测和判定算法6.1.1简单的检判算法传统的啸叫检测是利用FFT算法实现的，用FFT算法对每帧信号xn，n二。丄，N -1进行频谱分析，求得对应的功率谱Xk,k二N/2-1，然 后计算出判决值R,计算公式如下：R = 10lg1（Ppeak / Pa）（1）在1式中，Ppea是信号功率谱中的最大幅值，所对应得频率点是卩卩碱，pav 为每帧信号中所有功率幅值的平均值。如果在连续的一段时间，判决值都比预设 门限值T高，且最大幅值对应的频点Fpeak都为同一个值，则判定Fpeak为啸叫频点。6.1.2细化频谱分析算法用上面
9、的算法对声音信号的频谱进行分析时，分辨率为：D = F /N = 48KHz/1024 = 46.875Hz（2）f s显而易见，该算法的频率分辨率过低，我们可以采用CZT算法来细化频谱分 析。下面简单的介绍一下CZT算法：CZT（Chirp Z-transform）算法的实质就是对z变换采用螺线抽样。z变换的 定义式为：其中，xnn二丄,N -1是有限长序列。若令抽样点zk= AW-k，k = 丄,M -1，沿z平面上的一段螺线作等分角抽样，就可以得到 抽样点处的z变换：X =七xA-nWnk,k = ,1,M -1（4）zknn=其中，A = Aej2叫,W = We-j2叫,A为起始抽样
10、点Z的矢量半径长度，为起始抽样点的相角，W为螺线的伸展率（当W 1，为螺线内缩；当W，为逆时针；当e，为顺时 针），M为要分析的复频谱点数。令0 = Fl/Fs e =（化-F）/Fs，就可以得到巧口化范围的频谱，CZT算法的频率分辨率为（F厂F）/Fs。可见，CZT算法可以 任意选取F） Fh的频率范围进行分析，可以根据所需的频率分辨率来确定频率取 样点数M。CZT算法的原理示意图如图5所示:图5 CZT算法原理图CZT算法细化频谱效果如图6所示，第一幅图是直接利用FFT变换后的频谱图，第二幅图是利用CZT变换后的细化频谱图。在对啸叫信号进行频域分析时，如果使用CZT对整个频域进行频谱分析，
11、那利用C可变换后的细化频谱么需要N*M(Fh - Fl)/M，其中 F 能通过增加频率采样点数M来实现率的公式Dzt就会很大，使运算量大大增加。谱分需要的mH61鲁；为整个频带范围，要想提高分辨率，只啸叫检测系统要求的频率分辨率又很高,综上所述，FFT算法频率分辨率太低法运算量又太大。综合考虑fFt算法和czT算法的优缺点，我们采用一个折中的系统设计的要求，同时CZT算次复乘，(N-1)M 次复加运算量非常大。按照 CZT算法频率分辨办法：先利用FFT算法对整个频域进行粗分析，一旦检测到疑似声反馈频率点Fc 存在，就利用Fc设定CZT算法的F1和F，其中F1 = Fc-512，Fh = Fc+
12、512，然后对F1Fh频带范围进行频谱细化分析，如图7所示。这种算法在保证较高频 率分辨率的同时减小了 CZT算法的运算量，综合性能较好。1024 CZTHzFc-512FcF.+51215k图7 FFT+CZT算法示意图7、声反馈抑制算法实现7.1方案简述本实验采用参数可调的IIR陷波器来消除声反馈频点的自激啸叫。陷波器实 质是一种带阻滤波器，它的阻带很窄，近似于一个特定的频点，因此又称为点阻 滤波器。当声反馈检判算法检测到一个声反馈频点，就设计一个中心频率为该频 点的陷波器对其进行抑制。72陷波器设计数字陷波滤波器的设计要从滤波器的类型、模拟滤波器的设计准则以及模拟 滤波器到数字滤波器的映
13、射方法三个方面来考虑。陷波器的设计可以选用有限冲击响应FIR (Finite Impulse Response)滤波器 或者无限冲击响应IIR (Infinite Impulse Response)滤波器。先比较一下两者的优 缺点：FIR滤波器有严格线性相位，稳定，适合定点DSP实现等优点，但同时它 运算量太大且频率选择性不高；相比之下，IIR滤波器有良好的频率选择性，在 相同的技术指标下，可用比FTR滤波器少得多的阶数实现要求，计算量小，有利 于陷波器设计的实时性。综合考虑两类滤波器的优缺点，我们最终选用采用IIR 滤波器设计陷波器。IIR滤波器设计的基本步骤：第一，将数字滤波器技术指标，比
14、如通带截止 频率，阻带截止频率等，转换成对应的模拟滤波器技术指标；第二，设计模拟滤 波器Ha(s)使其满足上述技术指标;第三，将设计出的模拟滤波器系统函数Ha(s) 按某种方法转换成数字滤波器的系统函数H(z)。在整个设计过程中存在着两种 不同的变换：一是频带变换，即将低通滤波器变换成高通、带通、带阻滤波器：二是变换域变换，即将模拟滤波器变换成数字滤波器。根据这两种变换的先后顺 序，引出以模拟滤波器为基础的两种设计方法，如图9和10所示。模拟域频带变换变换域变换模拟低通滤波器设计低通变带阻模拟变为数字口尺数字 滤波器模拟低通滤波器设计图9 IIR数字滤波器设计方法1变换域变换模拟域频带变换模拟
15、变为数字低通变带阻IIR数字 滤波器图10 IIR数字滤波器设计方法2模拟低通滤波器几种常见的类型有贝塞尔滤波器、巴特沃思滤波器、切比雪 夫I型、II型和椭圆滤波器。它们的性能比较如表1所示。模拟低通滤波器类型滤波器类型贝塞尔滤波器近似线性相位，微小的群延时失真巴特沃什滤波器平滑的通带/阻带增益响应，适度的群延时变化切比雪夫I型滤波器通带具有等波纹，阻带平滑切比雪夫II型滤波器阻带具有等披纹，通带平滑椭圆函数型滤披器通带与阻带都有波纹，对给定任务提供最低阶数设计表1五类模拟低通滤波器性能对比图从模拟滤波器到数字滤波器，有两种常用的映射方法：脉冲响应不变法和双 线性变换法。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位 脉冲响应h(n)模仿模拟滤波器的冲击响应ha(t)。只有当模拟滤波器的频率响应是 带限的，且带限于折叠频率以内，才能使数字滤波器的频率响应在折叠频率以内
