面向声品质优化的有源噪声控制算法研究
发布时间：2021-09-24
摘要：传统有源控制算法未考虑控制后的声品质特性，而误差滤波最小均方(Filter-errorleastmeansquare，FeLMS)算法通过对误差滤波器的设置可较好地优化声品质。已有FeLMS算法的优化目标大多仅基于心理声学客观参量，未考虑不同类型噪声的频谱差异。针对该问题，以直升机舱内噪声为研究对象，研究面向声品质优化的有源控制算法。首先将直升机舱内人耳位置处采集的噪声样本作为初级噪声，开展烦恼度主观评价实验，然后对烦恼度实验数据进行多元线性回归得到直升机舱室噪声烦恼度的声品质模型，最后根据模型中烦恼度得分的频率响应曲线进行滤波器设计，以此为基础优化FeLMS算法。计算机仿真和主观评价实验结果表明，进行控制后噪声烦恼度明显降低。相较于传统滤波x最小均方(Filter-xLeastMeanSquare，FxLMS)算法，FeLMS算法在提升声品质方面效果更为明显。
关键词：声学;声品质;有源噪声控制;FeLMS算法;声品质模型;滤波器设计
噪声控制分为无源噪声控制(Passivenoisecontrol，PNC)和有源噪声控制(Activenoisecontrol，ANC)两种手段。无源噪声控制主要包括吸声处理、隔声处理、使用消声器、振动的隔离与降低等。这些控制方法的机理在于使噪声声波与声学材料或结构相互作用而消耗能量，从而达到降噪的目的，且对降低中高频的噪声较为有效。对于低频的噪声，更多是采用有源控制的办法。
在空间某一点，通过初级噪声与次级噪声的相消性干涉达到降噪目的的噪声控制方式称为有源噪声控制[1]。传统的有源噪声控制以相消干涉后噪声的声功率为目标函数，通过使抵消后噪声的声功率最小化从而起到降噪效果。降噪后，由于噪声中低频成分的大幅减少，中高频成分占主导，尽管声压级明显降低，但声品质可能并没有明显改善，甚至听起会来更不舒服。
面向声品质优化的有源噪声控制是一个重要的发展趋势。目前关于声品质控制方面的研究可大致分为两类。第一类是采取传统的有源控制算法，之后对控制后的噪声进行主观声品质评价。刘宗巍等[2]对控制后的响度、尖锐度以及基于这两个评价指标所生成的综合指标的变化进行了评估。Oliveira等[3]对有源控制后噪声的响度与粗糙度的变化进行了评估。Canvet等[4]对有源控制后响度、烦恼度和舒适性的改善进行了评估。Lin等则对降噪后的语义清晰度进行了评价[5]。另一类则是通过对控制算法进行改进，进而改变控制系统结构从而提高声品质。Kuo等提出FeLMS算法[6]和主动噪声均衡(Activenoiseequalizer，ANE)系统[7–8]。ANE系统可以衰减或放大次级噪声，从而为改善声品质提供了方法和思路。FeLMS算法是在FxLMS算法基础上加入滤波器，对抵消后的噪声进行滤波，进而改善声品质。关于声品质控制算法方面的后续研究大多都是在这两种算法基础上进行优化和拓展。Bao等[9]基于A计权曲线对误差信号和参考信号进行滤波。Wang等[10]对ANE系统增益系数进行了研究，根据临界频带将汽车内噪声划分为24个频段，并分别设置增益系数，以降低车内的响度和尖锐度。Sommerfeldt等研究了FeLMS算法对于降低响度的有效性[11]。姜顺明等[12]基于等响曲线对误差信号和参考信号进行滤波。
在基于声品质控制算法的研究中，大部分研究工作中的优化目标只从声学客观参量出发，并没有考虑到不同噪声的频谱差异。针对这一缺陷，本文以直升机舱室噪声为研究对象，设计烦恼度评分实验，建立声品质模型，以烦恼度为控制目标，对烦恼度在Bark域上进行计算并根据结果设计滤波器和基于FeLMS算法的控制系统，通过计算机仿真验证其对直升机舱室声品质改善的效果。
1自适应有源控制算法
自适应滤波是能自动调节其自身单位脉冲响应以达到最优化的维纳滤波器。设计自适应滤波器时可以不必要求预先知道信号与噪声的自相关函数，而且在滤波过程中即使这些自相关函数缓慢变化，也能自适应调节，自动调节到满足最小均方误差的要求。调整自适应滤波器权系数的方法称为自适应算法。
1.1最小均方(Leastmeansquare,LMS)和FxLMS算法
自适应滤波的最终目的就是寻求最佳权矢量WO。LMS算法正是求WO的一种简单有效的递推方法。LMS运用最陡下降法获得权系数递推公式。
1.2FeLMS算法
FeLMS算法是在FxLMS算法的基础上加入了误差滤波器，对误差信号e进行滤波。其算法框图如图2所示。图2中Hw(z)为加入的误差滤波器，r′(n)为参考信号x(n)通过次级通路传递函数C(z)和误差滤波器Hw(z)的信号，e′(n)是误差信号e(n)经过误差滤波器Hw(z)滤波后的信号。
2舱室噪声烦恼度评价与声品质模型构建
获取直升机舱室噪声的声样本和对其进行声品质模型的构建是整个实验的基础，本节对噪声样本的获取和声品质模型构建进行详细描述。
2.1舱室稳态噪声样本采集
选取某型直升机，分别在其开车、巡航、悬停3种工况下测量3分钟。针对直升机舱室降噪的需求，应选择直升机舱内乘员正常坐姿情况下人耳所在位置附近的局部区域作为测点，测量这些测点处的噪声。测试采用多通道B&KPulse3660D数据采集设备记录声压信号，时域数据采样频率为32768Hz，共采集87个声样本。一般情况下，对于稳态噪声，用于主观评价的样本时长为5s较为合适。
2.2舱室噪声声品质评价实验
选取了24位被试进行主观评价实验，其均为在读本科生和研究生且听力正常。地点选择在西北工业大学航海学院五楼会议室，声样本经动圈式高保真监听级头戴耳机(SENNHEISERHD280)播放给被试。
参考评分法结合了评分法和成对比较法的优点，适合于没有经过专业培训的被试，它对被试记忆力要求较低，而且实验结果受其他样本影响较小，实验评价结果精度较高，可用于声样本很多的情况，非常适合本次烦恼度评价实验。因此，最终选择参考评分法进行实验。
通过以上主观评价实验，共获得24份问卷。为保证后续分析结果的有效性，通过误判分析、相关分析和对被试评分范围的考察对数据进行有效性检验及一致性检验。
被试对烦恼度的打分应在1～9之间，若被试的评分较为集中，则可认定为无效的数据，应当予以剔除。误判分析用于判断同一被试两次评价结果的一致性。若同一被试对同一声样本多次重复评价结果存在较大差异，则可认定为无效数据，应当予以剔除。相关分析用于判断同一被试两次评价结果之间的相关性。用同样的量表对同一组受试者进行两次重复测量，计算各项得分之间相关分析可以说明该量表的测量信度，相关系数低说明被试在评价过程中评价尺度处于一个不固定的状态，则可认定数据无效，应当予以剔除。误判分析和相关分析结果如图3和图4所示。
经过上述数据分析，最终共获得17名有效被试烦恼度评价结果。
2.3声品质模型的构建
烦恼感是噪声所具有的典型感知属性，源于噪声引起的各种即时效应，是公众对噪声的理性认识，与噪声刺激直接相关。噪声烦恼度与很多因素相关，声音的响度、尖锐度、粗糙度等都会对声音的烦恼度产生影响。本文以被试的烦恼度得分均值为因变量，选取响度L、尖锐度S、粗糙度R、波动强度FL为自变量，建立直升机舱室噪声声品质模型。选取所有声样本烦恼度评分均值的80%用于训练，20%用于检验。烦恼度得分的多元线性回归模型见式(6)。
相关期刊推荐：《噪声与振动控制》(双月刊)主要内容有噪声与振动控制专题论文;噪声振动治理的新理论、新方法、新经验、新的设计思想及各种工程实例;基础理论讲座;噪声振动测试技术;国内外噪声振动控制元器件、新技术、新材料、新产品以及工厂介绍等。可供相关专业的大专院校师生、各单位噪声振动专业人员参考。
为了验证模型的适用性，将预留的烦恼度主观评分数据与通过所建立模型计算的烦恼度数值对比，进而对声品质模型进行检验。绘制了实验值与预测值对比折线图，如图5所示。
由图5可以看出，实验值与模型预估值大致相同。虽有少量声样本的计算值和实验值相差较大，但整体趋势接近。对实验值和模型预测值进行相关分析，两者的相关系数为0.9186，表明所建立的模型对于不同工况下的直升机舱内噪声具有较好的有效性和适用性。
3基于声品质的有源噪声控制算法及仿真
本节选取主观评价实验中的某一声样本，根据声样本频谱并结合声品质模型设计滤波器。通过计算机仿真比较滤波器加入前后系统的控制效果即烦恼度的变化。
3.1滤波器设计
FeLMS算法中大多采用基于A计权曲线或等响曲线设计的滤波器。但都只考虑到人耳特性，并没有考虑到不同噪声的频谱差异，对差异较大的噪声控制效果相差也较大。选取主观评价实验中的某一声样本，根据声样本频谱并结合烦恼度模型设计滤波器。选取的噪声时频图如图6所示。基于Zwicker[13]模型对选取的声样本进行心里声学客观参量计算，所得计算结果如表1所示。
声样本的响度、粗糙度、尖锐度和波动强度在每个Bark域内的计算值如图7所示。
将声品质客观参量在每个Bark域内的计算数值代入声品质模型中，可得声样本在每个Bark域内的烦恼度得分，如图8(a)所示。为了使烦恼度得分与滤波器的响应幅值匹配，对各Bark域内的烦恼度得分进行数据归一化处理，如图8(b)所示。
考虑到高频降噪在有源控制系统中占比很小，且系统控制效果随频段的增加而衰减。故选取前16个Bark域，即频率在0～3500Hz范围内的烦恼度得分，按照归一化后的分值在频域内进行滤波器拟合，拟合成的滤波器在各频段响应幅值与烦恼度得分大致相等。拟合成的滤波器频响曲线如图9所示。
图9中，虚线代表滤波器频响曲线，黑色柱状线代表烦恼度得分。结合有源噪声控制主要控制中低频段的特点，适当调高滤波器在50Hz和500Hz两个频率附近的响应幅值，滤波器在其他频段的响应幅值与烦恼度得分基本拟合。
3.2仿真及结果对比分析
对选取的噪声样本在MATLAB中进行有源控制仿真，分别采用FxLMS和FeLMS两种算法。其中在FeLMS算法中分别设置以A计权曲线拟合设计的滤波器和以声品质模型拟合设计的滤波器。为了叙述方便，下文将基于A计权曲线的FeLMS算法记为A-FeLMS，将基于声品质模型的FeLMS算法记为Y-FeLMS。分别对基于3种算法的控制系统进行仿真分析与比较。控制后的噪声频谱如图10所示。
可以看出，经基于Y-FeLMS算法的控制系统控制后的噪声在100Hz以下的低频段和1800Hz附近频段处声压级相比前两种算法较大，而在900Hz附近的频段声压级相对较低。对控制后的噪声进行计算，噪声的心理声学客观参量计算结果如表2所示。
从表2中可以看到，经基于FeLMS两种算法的控制系统控制后的降噪量分别为24.1dB和24.5dB，相比于传统FxLMS算法27.2dB的降噪量并没有优势。但在烦恼度方面，基于FxLMS算法的烦恼度下降到0.4707，A-FeLMS下降到0.3901，Y-FeLMS算法下降到0.2457。可以看出Y-FeLMS相比于前两种算法实现了更好的声品质控制。
3.3主观评价验证为验证
控制效果计算的准确性，完全按照之前评价实验的客观条件再次组织声品质评价实验。要求被试对基于3种算法控制后的噪声样本进行评分。主观评价结果如图11所示。
根据图11中的实验值可看出，相比较于FxLMS算法和A-FeLMS算法，经基于Y-FeLMS算法控制后的烦恼度多下降了0.64和1.12个等级。结果表明，根据声品质模型优化的FeLMS算法有助于改善直升机舱室内的声品质。
4结语
通过对噪声样本的声品质客观参量分析，建立声品质模型，并根据所建立的声品质模型设计滤波器，并以此为基础优化FeLMS算法，提升有源噪声控制技术在改善声品质方面的表现。本文研究主要工作和结论如下：
(1)以某型直升机在开车、巡航和悬停3种工况下舱室内的稳态噪声为研究对象，采集并制成87个声样本。通过主观声评价实验和多元线性回归方法建立直升机舱内噪声声品质模型，并对其进行验证。结果表明实验真值和模型预测值相关系数为0.9186，具有较高的可靠性。
(2)根据所建立的声品质模型设计滤波器，对FeLMS算法进行优化，并建立基于烦恼度的YFeLMS算法有源控制系统。仿真结果表明，经过控制后，烦恼度相较于FxLMS算法和基于A计权曲线的A-FeLMS算法多下降0.64和1.12个等级，控制效果有比较明显的提升。——论文作者：来昊，陈克安，王磊
