1、2023年第47卷第4期132器 件 与 应 用器 件 与 应 用arts and ApplicationsP文献引用格式：吕雨贞，刘紫赟，沈勇.基于听觉特性优化的级联型均衡滤波器设计 J.电声技术，2023，47（4）：132-136.LV Y Z，LIU Z Y，SHEN Y.Design of cascaded equalization filter based on auditory characteristics optimizationJ.Audio Engineering，2023，47（4）：132-136.中图分类号：TN431.2 文献标识码：A DOI：10.16311/
2、j.audioe.2023.04.038基于听觉特性优化的级联型均衡滤波器设计吕雨贞，刘紫赟，沈 勇*（南京大学声学研究所 近代声学教育部重点实验室，江苏 南京 210093）摘要：受扬声器品质、房间反射等因素的影响，声音信号会出现失真。为了改善声重放效果，使用均衡滤波器对扬声器的频率响应进行均衡。考虑到人耳听觉特性和滤波器参数估计，提出了对数尺度上系统频率响应与目标频率响应之差的代价函数，设计了由若干个 Regalia-Mitra 滤波器级联而成的均衡滤波器，并将其应用于小房间的音箱频率响应均衡。结果表明，设计的均衡滤波器能较好地减少扬声器频率响应与目标频率响应的差异。关键词：扬声器频率响应
3、；人耳听觉特性；均衡滤波器；响应均衡Design of Cascaded Equalization Filter Based on Auditory Characteristics OptimizationLVYuzhen,LIUZiyun,SHENYong*(KeyLaboratoryofModernAcoustics,MOE,andInstituteofAcoustics,NanjingUniversity,Nanjing210093,China)Abstract:Affectedbyfactorssuchasspeakerqualityandroomreflection,thesound
4、signalwillbedistorted.Inordertoimprovethesoundreproductioneffect,anequalizationfilterisusedtoequalizethefrequencyresponseofthespeaker.Consideringtheauditorycharacteristicsofhumanearsandtheestimationoffilterparameters,thecostfunctionofthedifferencebetweenthesystemfrequencyresponseandthetargetfrequenc
5、yresponseonthelogarithmicscaleisproposed,andanequalizationfiltercomposedofseveralRegalia-Mitrafiltersisdesignedandappliedtothefrequencyresponseequalizationofspeakersinsmallrooms.Theresultsshowthatthedesignedequalizationfiltercanreducethedifferencebetweenthespeakerfrequencyresponseandthetargetfrequen
6、cyresponse.Keywords:speakerfrequencyresponse;hearingcharacteristicsofhumanear;equalizingfilter;responseequalization0 引 言由于实际的扬声器可能存在缺陷，加上房间反射等因素的影响，听者听到的声音存在失真。为了改善声重放效果并使听者获得期望的听音感受，利用均衡滤波器对扬声器频率响应进行补偿和优化。均衡滤波器通常由若干个低阶无限冲激响应（Infinite Impulse Response，IIR）滤波器级联而成，低阶 IIR 滤波器有一阶搁架滤波器和二阶峰值滤波器。一阶搁架滤波器的重
7、点参数为拐点频率和增益，二阶峰值滤波器的重点参数为中心频率、增益以及品质因子。通过合理调整级联的低价 IIR 滤波器参数，可以有效降低均衡滤波器的阶数，实现更好的均衡。手动调整滤波器参数值较为低效，误差难以保证，因此设计自动均衡滤波器对参数进行调整。自动设计均衡滤波器一般基于代价函数最小化进行设计，已有的研究通常使用系统频率响应和目标频率响应的幅度差作为代价函数1-2。VAIRETTI G提出一种新的代价函数，即系统频率响应与目标频率响应差3。该代价函数的优点是既考虑了幅度误差，又考虑了相位误差，同时比只使用频率响应幅度误差的代价函数有更高的数学可处理性。实际上，人耳对声音强度的感知并不是线性
8、，基金项目：国家重点研发计划项目(2018YFB1403800)作者简介：吕雨贞（1998），女，博士在读，研究方向为声学。E-mail：。通信作者：沈 勇（1965），男，博士，教授，研究方向为声学。E-mail：。2023年第47卷第4期133Parts and ApplicationS器 件 与 应 用器 件 与 应 用而是近似满足对数关系4-5。考虑到无法对复数直接取对数，同时为保证代价函数的可导性，可以适当压缩频率响应误差，以匹配人耳的感知差异。基于人耳听觉特性，提出在对数尺度上系统与目标频率响应差的代价函数，使用此代价函数计算滤波器参数，从而实现具有更好听音效果的均衡滤波器设计。1
9、 Regalia-Mitra 结构的滤波器文章设计的均衡滤波器由一阶搁架滤波器和二阶峰值滤波器级联而成，结构如图 1 所示，包括一个全通滤波器和一个前馈路径6。图 1 Regalia-Mitra 滤波器Regalia-Mitra 滤波器的传递函数 Fm(z)为()()()()()()()()()1LFLF11LF1HFHF11HFLFHF1112mmmFzAzVAz=+（1）式中：V 为滤波器增益，Am(z)为全通滤波器传递函数。采用线性增益的 Regalia-Mitra 滤波器结构，此时全通滤波器传递函数 Am(z)与增益 V 无关，IIR滤波器传递函数 Fm(z)与增益 V 呈线性关系。分
10、别设置一阶全通滤波器传递函数 Am(z)和增益 V 可得一阶搁架滤波器，分别设置二阶全通滤波器传递函数 Am(z)和增益 V 可得二阶峰值滤波器。低频搁架滤波器和高频搁架滤波器的一阶全通滤波器传递函数分别定义为()1LFLF11LF1HFHF11HFLFHF1112111tan 1tan tan 1tan 1mmmcscscscsFzAzVAzazAza z=+=（2）()()()()()1LFLF11LF1HFHF11HFLFHF1112111tan 1tan tan 1tan 1mmmcscscscsFzAzVAzazAza zazAzazffaffffaff+=+（3）aLF和 aHF分
11、别为7()()()()1LFLF11LF1HFHF11HFLFHF1112111tan 1tan tan 1tan 1mmmcscscscsFzAzVAzazAza zazAzazffaffffaff=+=+=+=+（4）()()HFtan 1tan 1cscsffaff=+（5）式中：fs为采样频率，fc为一阶搁架滤波器的拐点频率。二阶全通滤波器的传递函数定义为()()()12212111ada zzAzda zaz+=+（6）式中：d=-cos(2f0/fs)。a 的表达式为()()()()()()()()()()()()12212dB1201dB,121111tan 1tan 1bsbs
12、NMkMSMmsNs mksmada zzAzda zazffaffENW kHkFkT kFkCFkENW kHk FkT k=+=+=+（7）式中：f0为二阶峰值滤波器的中心频率，fb为二阶峰值滤波器的-3 dB 带宽。2 基于人耳听觉特性的代价函数与增益估计使用对数函数模拟人耳的听觉特性，提出代价函数为()()()()()()()()()()12212dB1201dB,121111tan 1tan 110log110log1sssbsbsNMkMSMmsNs mksmada zzAzda zazffaffENW kHkFkT kFkCFkENW kHk FkT k=+=+=+=+（8）式
13、中：k 为离散频率索引；N 为需要计算的频率个数；ms为系统频率响应；T(k)为目标频率响应；FM(k)为 M 阶均衡滤波器的频率响应；W(k)为频率加权函数，表明不同频率的重要性。M 阶均衡滤波器由 S 个低阶 IIR 滤波器级联，其频率响应 FM(k)为()()()()()()12212dB1201dB,121111tan 1tan 110log110log1sssbsbsNMkMSMmsNs mksmada zzAzda zazffaffENW kHkFkT kFkCFkENW kHk FkT k=+=+=+=+=+（9）式中：C 为均衡滤波器全局增益；s 为级联的 IIR滤波器索引；m
14、s为第 s 个 IIR 滤波器的阶数，一般为 1 或 2；Fms(k)为第 s 个 IIR 滤波器的频率响应。对第 s 个滤波器的参数进行优化，使代价函数最小。此时代价函数可表示为()()()()dB,12110log1ssNs mksmENW kHk FkT k=+（10）式中：Hs-1(k)为前 s-1 个 IIR 滤波器级联的均衡滤波器均衡后的频率响应。2023年第47卷第4期134器 件 与 应 用器 件 与 应 用arts and ApplicationsP由式（1）可得()()()()()()()()()()()()()()()()()()()2*022202*0220 01112
15、12smmmmmFkAVAFkV Fkkk=+=+（11）式中：Fm(k)=1+Am(k)，Fm(k)=1-Am(k)。当全通滤波器传递函数 Am(z)已知时，令代价函数对增益 V 的一阶偏导数为 0 来估计增益 V。为了方便计算，利用 Fm(k)和 Fm(k)正交8，增益 V的估计为()()()()()()()2*022202*0220 01112122smmmmmmkmkkksbsFkAVAFkV FkW kFk Hk T kVW kHkFkW kHk T kCW kHkffQaffak=（12）类似于滤波器增益 V 计算，均衡滤波器全局增益 C 的估计为()()()()()()()2*0
16、220 01112122 sin 2sin 12tan 21kksbsFkAVAFkV FkW kFk Hk T kVW kHkFkkW kHk T kCW kHkffQaffak=|+（13）3 均衡滤波器设计均衡滤波器由 S 个低阶 IIR 滤波器级联，根据式（10）依次设计每个低阶 IIR 滤波器。3.1 低阶 IIR 滤波器参数的初值选取低阶 IIR 滤波器的参数可分为全通滤波器参数和滤波器增益 V。由式（12）可知，利用全通滤波器参数和待均衡频率响应可以估计滤波器增益 V，因此只需设置全通滤波器参数。设置极点网格，一个极点对应一个全通滤波器参数，根据待均衡频率响应构造极点网格对应的滤波器组。计算滤波器组的代价函数，选择其中代价函数最小的低阶 IIR 滤波器参数作为第 s 个低阶 IIR 滤波器的初值。根据目标优化的频率范围，设定搁架滤波器的拐点频率 fc和峰值滤波器的中心频率 f0、Q 值分布，生成极点网格。复极点 p=ej对应峰值滤波器参数d=-cos()和 a=2。正实极点 p=对应低频搁架滤波器参数 aLF=2，负实极点=-对应高频搁架滤波器参数 aHF=2。通过设定