Loading [MathJax]/jax/output/HTML-CSS/jax.js
高级搜索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种无重叠频带的自适应非对称并联图示均衡器校正方法

李亚 杨俊杰 冯奇 秦先清

李亚, 杨俊杰, 冯奇, 秦先清. 一种无重叠频带的自适应非对称并联图示均衡器校正方法[J]. 电子与信息学报, 2022, 44(5): 1734-1742. doi: 10.11999/JEIT210220
引用本文: 李亚, 杨俊杰, 冯奇, 秦先清. 一种无重叠频带的自适应非对称并联图示均衡器校正方法[J]. 电子与信息学报, 2022, 44(5): 1734-1742. doi: 10.11999/JEIT210220
LI Ya, YANG Junjie, FENG Qi, QIN Xianqing. An Adaptive Asymmetric Parallel Graphic Equalizer Correction Method without Overlapping Frequency Bands[J]. Journal of Electronics & Information Technology, 2022, 44(5): 1734-1742. doi: 10.11999/JEIT210220
Citation: LI Ya, YANG Junjie, FENG Qi, QIN Xianqing. An Adaptive Asymmetric Parallel Graphic Equalizer Correction Method without Overlapping Frequency Bands[J]. Journal of Electronics & Information Technology, 2022, 44(5): 1734-1742. doi: 10.11999/JEIT210220

一种无重叠频带的自适应非对称并联图示均衡器校正方法

doi: 10.11999/JEIT210220
基金项目: 广州市科技计划(201904010302),广东省教育厅科研项目(2021ZDZX1079),广东技术师范大学人才专项基金(2021SDKYA056)
详细信息
    作者简介:

    李亚:男,1979年生,副教授,研究方向为主动健康技术、智能信息处理、忆阻混沌神经网络等

    杨俊杰:男,1995年生,硕士生,研究方向为音频信号处理

    冯奇:男,1989年生,讲师,研究方向为智能信息处理、压缩感知信号处理等

    秦先清:男,1979年生,总经理,研究方向为音频信号处理

    通讯作者:

    冯奇 qfeng@gpnu.edu.cn

  • 中图分类号: TN912.3

An Adaptive Asymmetric Parallel Graphic Equalizer Correction Method without Overlapping Frequency Bands

Funds: The Science and Technology Program of Guangzhou (201904010302), The Scientific Research Project of Department of Education of Guangdong Province (2021ZDZX1079), The Doctoral Scientific Research Startup Fund of Guangdong Polytechnic Normal University (2021SDKYA056)
  • 摘要: 针对车载音响的声场校正精度低、效率低的问题,该文提出一种无重叠频带的自适应非对称并联图示均衡器校正方法。在车载音响声场动态变化情况下,该方法考虑了有效均衡范围和自适应增益,而不是经典方法中的固定均衡范围和人工给定的增益。通过对实测数据进行实验分析,所提方法比经典方法所用均衡滤波器个数平均减少大约20%,拟合目标增益更准确,校正后频谱曲线更平坦。
  • 音频均衡的目的在于提高收听环境中的声音质量[1,2]。随着汽车的普及,车载音频设备有着广阔的市场前景。由于狭小封闭的车内空间、有限的扬声器安装位置以及汽车行驶时的胎噪路噪,导致车载音频音质下降[3,4]。为了提高车载音响的音质,一般通过调节车载音频处理系统中的核心部件、均衡器来实现[1,2]。通过均衡器直接从音源发声端对音频信号的频率响应曲线进行修正、补偿等调节,可使人耳听到的音频信号的频响曲线更加平直、均匀,使驾乘人员能够听到更加动听、真实的声音。

    音频均衡器是对声频中的某些频段进行提升或衰减的声电处理设备,一般分为图示均衡器和参数均衡器[1,2]。图示均衡器为中心频率与带宽固定而增益可调的滤波器组,参数均衡器为中心频率、带宽和增益均可进行调整的滤波器组[1]。图示均衡器结构简单,直观明了,在专业音响中应用非常广泛。参数均衡器调节灵活,可以完成更细致的修补工作,满足消费者对高品质车载音效的需求,是未来车载音频均衡器发展的方向。

    车载音频系统失真通常为线性失真,对应为系统频率响应的幅度和相位失真。尽管相位失真在某些情况下可以被察觉,但其影响远小于幅度失真的影响。因此,音频均衡器应专注于校正系统的幅度响应,而不是其相位响应[5,6]

    在音频均衡算法设计方面,目前主要有3类算法。第1类是基于智能优化算法:Pepe等人[3]提出了一种基于遗传进化算法的多声道音频均衡技术来调整均衡滤波器系数,该方法可以找到全局最优解,但效率低于常规的优化方法;Prince等人[7]提出了一种基于遗传进化算法的高阶无限脉冲响应(Infinite Impulse Response, IIR)滤波器均衡方法,该方法计算成本较高;Pepe等人[8]提出了一种基于引力搜索算法的音频均衡方法,该方法迭代次数较多,计算量较大。第2类是神经网络算法:Pepe等人[9]提出了一种基于深度神经网络算法来设计音频均衡器;Rämö等人[10,11]提出了一种基于神经网络的1/3倍频程图示均衡器校正方法;Välimäki等人[12]通过训练一个神经网络来拟合图示均衡器的目标增益;Ramírez等人[13]提出了一种端到端的卷积神经网络来设计图示均衡器。尽管神经网络算法具备强大的非线性拟合能力,可以任意精度逼近任何非线性连续函数,但其网络模型复杂度较高,且对硬件设备的计算能力有较高的要求。第3类是滤波器组算法:Rämö等人[14]提出了一种用于音频处理的高阶图示均衡器设计算法,该方法拟合误差在±2 dB以内,不满足高保真音频要求。现有的低阶图示均衡器校正方法有串联图示均衡器[15-17]和并联图示均衡器[18-23]。Liski等人[15-17]提出了一种精确的串联图示均衡器校正方法,该方法拟合误差在1 dB以内,满足高保真音频要求;文献[18-20]提出了一种将串联图示均衡器转化为延迟并联图示均衡器的方法,拟合误差也在1 dB以内;Rämö等人[22]提出了一种高精度的并联图示均衡器校正方法,该方法拟合误差虽然在1 dB以内,但是设计过程较复杂;Bank等人[21]对文献[22]方法进行了改进,减少了设计时间;Tao等人[19]对文献[22]的滤波器阶数进行了优化,降低了滤波器的总阶数;Chen等人[23]提出了一种基于预失真的并联图示均衡器校正方法,该方法拟合误差达到1.7 dB,不满足高保真音频要求。

    综上所述,现有的较低计算成本的先进的图示均衡器校正方法有Liski等人[16-18]提出的精确的串并联图示均衡器和Rämö等人[21,22]提出的高精度并联图示均衡器。这3种设计方法拟合误差都在1 dB以内,符合高保真音频要求。然而,在车载声场环境下,这3种设计方法都不能适应车载音响动态变化的声场,因为这3种方法都是针对静态的情况。

    针对以上问题,本文提出一种能够适应车载音响声场动态变化的高精度并联图示均衡器校正方法。实验结果表明,本文所提方法比文献[16,18,21,22]方法所用均衡滤波器个数平均减少大约20%,拟合的目标增益更准确,校正车载扬声器声场效果更好。

    本文组织如下:第2节介绍图示均衡器的串并转换方法;第3节阐述所提方法;第4节为实验结果分析;第5节为本文结论。

    Liski等人[18]提出了一种将高精度的串联图示均衡器转化为延迟并联图示均衡器的方法。方法如下所述。

    一个2阶滤波器传递函数如式(1)所示[18]

    F(z)=k01+b1z1+b2z21+a1z1+a2z2 (1)

    其中,k0为增益参数,b1b2为前馈系数,a1a2为反馈系数。

    对于S个串联的2阶滤波器组,其传递函数为

    F(z)=G0Ss=11+b1,sz1+b2,sz21+a1,sz1+a2,sz2 (2)

    其中增益参数G0

    G0=Ss=1k0,s (3)

    图示均衡器串联组成示意图如图1所示。其中S是均衡滤波器的个数。单个2阶滤波器结构框图如图2所示。

    图 1  图示均衡器串联组成示意图
    图 2  单个2阶滤波器结构框图

    文献[18]给出的串并转换方法流程如下。首先,将式(2)写成非因果形式

    F(z)=G0Ss=1z2+b1,sz+b2,sz2+a1,sz+a2,s (4)

    假设F(z)的部分分式展开为

    F(z)=d0+Nn=1˜rnzpn (5)

    其中,d0为常数,pn为极点,˜rn为极点对应的留数。

    使用Heaviside cover-up方法计算留数为

    ˜rn=(zpn)F(z)|z=pn (6)

    接下来通过将式(5)转化为标准的因果滤波器形式引入延迟

    F(z)=d0+Nn=1˜rnzpn=d0+Nn=1z1˜rnz1(zpn)=d0+z1Nn=1˜rn1pnz1 (7)

    为了求解d0,联立式(2)和式(7)为

    F(z)=G0Ss=11+b1,sz1+b2,sz21+a1,sz1+a2,sz2 =d0+z1Nn=1˜rn1pnz1 (8)

    z得到d0=G0

    为了得到实值滤波器系数,需要组合复共轭极点对,如式(9)所示

    Fp(z)=˜r11p1z1+˜r21p2z1=(˜r1+˜r2)(˜r1p2+˜r2p1)z1(1p1z1)(1p2z1) (9)

    最后得到的延迟并联图示均衡器传递函数为

    F(z)=d0+z1Ss=1c0,s+c1,sz11+a1,sz1+a2,sz2 (10)

    其中,c0,s=˜r(2s1)+˜r(2s), c1,s=˜r(2s1)p(2s)˜r(2s)p(2s1)。通过以上步骤就将串联的2阶滤波器转化为延迟并联形式。图示均衡器并联组成示意图如图3所示。其中S是均衡滤波器的个数。单个2阶滤波器结构框图如图4所示。

    图 3  图示均衡器并联组成示意图
    图 4  单个2阶滤波器结构框图

    本节提出一种用于车载音频的自适应增益精确并联图示均衡器校正方法。所提设计中目标增益是自适应获得的,所设计的并联图示均衡器校正车载扬声器的声场效果较好。

    第2节介绍的方法需要首先设计精确的串联图示均衡器(Accurate Cascade Graphic Equalizer, ACGE)[16],接着通过其阐述的串并转化方法将ACGE转化为精确的并联图示均衡器(Accurate Parallel Graphic Equalizer, APGE)[18]。APGE方法设计出的图示均衡器与目标增益的误差在±1 dB以内,达到高保真音频要求。然而,该方法的目标增益只是人为地预先设定,并没有考虑到动态变化的情况。在车载音频声场环境下,车载扬声器的声场会随着车辆的高速行驶而不断变化,在这种情况下拟合目标增益效果会变差。

    本文所提出的方法在APGE基础上加入了车载扬声器有效频响范围识别和自适应计算目标增益,称为自适应增益-精确的并联图示均衡器(Adaptive Gain-Accurate Parallel Graphic Equalizer, AG-APGE)。AG-APGE方法减少了均衡滤波器个数,自适应计算的目标增益比人工给定的目标增益更准确,因此图示均衡器响应的拟合效果更精准,最终校正后的曲线更平坦。

    图5是本文算法的应用系统框图。车载音响系统主要由扬声器、功放电路和数字信号处理器(Digital Signal Processor, DSP)等部分组成,形成一个独立的系统。其中,DSP片上系统为整个车载音响系统的核心部分,负责音频的解码、信号处理以及均衡等功能。整个系统的流程为:首先车载音响系统播放特定的白噪声音频,然后通过麦克风收集到车内的环境影响后的失真音频信号。失真信号被传到电脑(Personal Computer, PC)端中的AG-APGE算法模块,该算法模块计算出最优的滤波器组参数后,PC端通过通信协议与车载音响系统的DSP片上系统进行通信,修改DSP系统中对应通道上的均衡滤波器参数,完成参数的自适应调整。

    图 5  车载音响图示均衡器自动校正算法应用系统框图

    设车载音响系统频率响应为H0(k),它表示待校正的频率响应。H0(k)是其传递函数在单位圆上的第k部分归一化弧度制角频率ωk/fs处的值,其中kejωk/fs的简写,fs是采样频率。设校正后频率响应为Hs(k),它表示H0(k)经过由s个均衡滤波器组成的图示均衡器滤波后的频率响应。设目标频率响应为T(k)。设校正频率响应为C(k)

    文献[18]的目标增益是人为预先设定的,不能适应车载音响声场的变化。并且文献[18]的均衡频率范围固定为20~20000 Hz。在车载声场环境下,由于扬声器对音频信号有其特定的频率响应范围,因此图示均衡器设计考虑的频率范围是扬声器的有效频率响应范围而不是固定的20~20000 Hz。AG-APGE方法能够适应车载音响声场的不断变化并且其均衡的频率范围是每个扬声器对应的有效频率响应范围,并不是固定的20~20000 Hz。

    因此,首先需要对车载扬声器进行有效频率响应范围识别。步骤如下:取车载音响系统幅度响应得到|H0(k)|,以其最大值做归一化处理并转化为对数幅度响应刻度,取归一化处理后的系统对数幅度响应的–15 dB(经验值)处对应的频率值,最小的频率值和最大的频率值即为扬声器通带范围的边界。根据扬声器的频率响应范围,可以将扬声器分为低音、中音、高音、中低音、中高音和全频段等类型。车载扬声器分类和常见的频率响应范围如表1所示。设扬声器上限截止频率为ef (end frequency),扬声器下限截止频率为if (initial frequency)。

    表 1  车载扬声器分类及其频响范围表
    车载扬声器分类常见频率响应范围
    低音20~150 Hz
    中音150 Hz~5 kHz
    高音5~20 kHz
    中低音20 Hz~5 kHz
    中高音110 Hz~20 kHz
    全频段50 Hz~20 kHz
    下载: 导出CSV 
    | 显示表格

    当扬声器识别为全频段扬声器时,取目标幅度响应为车载音响系统幅度响应在有效频率响应范围内的算术平均值。即如式(11)所示

    |T(k)|=meanωifωkωef(|H0(k)|) (11)

    其中,mean(·)表示取算术平均值。

    对数校正幅度响应20lg|C(k)|20lg|T(k)|20lg|H0(k)|相减,即如式(12)所示

    20lg|C(k)|=20lg|T(k)|20lg|H0(k)|,ωifωkωef (12)

    接下来,在20lg|C(k)|上计算目标增益点对应的增益值,以此自适应获得S个目标增益设置点。分别对应S个均衡滤波器。

    这部分首先利用3.1节得到的S个目标增益设置点和2阶IIR峰谷滤波器原型设计ACGE,再用第2节介绍的串并转化方法将其转化为APGE[18]

    本文使用的2阶对称IIR峰谷滤波器原型定义为[24]

    F(z)=1+Gβ2cos(ωc)z1+(1Gβ)z21+β2cos(ωc)z1+(1β)z2 (13)
    β={tan(B/2),G=1|G2B1||G2G2B|tan(B2), (14)

    其中,G为线性峰值增益,ωc为中心频率,B为带宽,GB为带宽处线性增益值。

    对于1/3倍频程情况,文献[18]方法定义了20lg(GB,s)=0.4×20lg(Gs)。其中s=1,2,,SGB,sGs分别为第s个频段的均衡滤波器的带宽增益和峰值增益。同时定义了每个频带滤波器的带宽为相邻频带滤波器中心频率的差值,对于1/3倍频程情况,有Bs=(321/32)ωc,s0.4662ωc,s

    1/3倍频程的中心频率fc和带宽fB表2所示。

    表 2  1/3倍频程中心频率和带宽
    频带编号fc(Hz)fB(Hz)频带编号fc(Hz)fB(Hz)
    1209.317800373
    22511.7181000466.2
    331.514.7191250582.8
    44018.6201600745.9
    55023.3212000932.4
    66329.42225001166
    78037.32331501469
    810046.62440001865
    912558.32550002331
    1016074.62663002937
    1120093.22780003730
    12250116.628100004662
    13315146.929125005828
    14400186.530160007459
    15500233.131200009324
    16630293.7
    下载: 导出CSV 
    | 显示表格

    其中,ωc=2πfc/fsB=2πfB/fsfs是采样频率。

    由于相邻滤波器之间的重叠频带影响,会导致实际图示均衡器幅度响应和目标增益误差较大,文献[18]通过使用交互矩阵来减小这种误差。文献[18]通过在目标增益点之间加入额外的频率点,频率点频率值为相邻目标增益点频率值的几何平均值,频率点增益值为相邻目标增益点增益值的算术平均值。具有2S1S列的交互矩阵定义为

    Dk,s=20lg|Fs(k)|gp (15)

    其中,k=1,2,,2S1并且s=1,2,,S。1/3倍频程中gp=17 dB。则均衡滤波器最优增益为

    gopt1=D+t=(DTD)1DTt (16)

    其中,D+D的伪逆矩阵,t为具有2S1个元素的列向量,其奇数行元素值为原始目标增益分贝值,偶数行元素值为相邻增益分贝值的算术平均值。

    为了进一步减小误差,需要将gopt1替代gp=17 dB作为新的原型增益去更新交互矩阵D。得到的新的交互矩阵定义为D1。则均衡滤波器最优增益为

    gopt2=D1+t=(D1TD1)1D1Tt (17)

    通过式(17)最终得到每个频带的滤波器最优分贝增益。将最优分贝增益转化为最优线性增益后带入式(13)和式(14)中,就设计出了S个串联的2阶对称峰谷滤波器组。最后,通过第2节介绍的串并转化方法就将该串联的2阶峰谷滤波器组转化为延迟并联形式。延迟并联形式的2阶带通滤波器中心频率处两侧的幅度响应是非对称的,即滤波器中心频率处两侧幅度响应的陡峭程度不一样。

    将所提方法设计的图示均衡器幅度响应去拟合对数校正幅度响应20lg|C(k)|,用图示均衡器幅度响应去校正车载音响系统对数幅度响应20lg|H0(k)|,最终得到校正后对数幅度响应20lg|HS(k)|,如式(18)和式(19)所示

    20lg|HS(k)|=20lg|H0(k)|+20lg|Ss=1˜Fs(k)| (18)
    20lg|HS(k)|20lg|T(k)| (19)

    其中,˜Fs(k)Fs(k)经过变换后的并联结构的滤波器频率响应,20lg|HS(k)|趋近于20lg|T(k)|

    20lg|HS(k)|20lg|T(k)|比较,评估校正效果。校正效果通过谱平坦度度量(Spectral Flatness Measure, SFM)指标衡量[25]。定义频率响应˜Hs(k)=Hs(k)T(k),则SFM定义为˜Hs(k)模平方的几何平均值与算术平均值之比,即如式(20)所示

    SFMs=NNNk=1|˜Hs(k)|2Nk=1|˜Hs(k)|2 (20)

    其中,N为在有效频率响应范围内取的频率值数量,SFMs表示经过s个图示均衡滤波器校正后的SFM

    本文通过MATLAB平台对采集的50个实测数据进行音频的分析以及算法的实现。所提方法AG-APGE和ACGE[16], APGE[18]以及文献[21]的并联图示均衡器设计(Parallel Graphic Equalizer, PGE)对比。使用的采样频率为192 kHz。

    图6显示了不同图示均衡器校正方法拟合滤波器目标增益的情况。对采集的某个实测数据进行扬声器频率响应范围识别可以得到该车载扬声器为全频段扬声器,频率响应范围是50~14000 Hz。在50~14000 Hz这段频率范围进行自适应校正。本文方法所用均衡滤波器个数是动态变化的,为10~25个。在该实验中,本文方法使用了25个2阶带通滤波器。而ACGE, PGE和APGE分别使用了31个2阶峰谷滤波器、62个2阶带通滤波器和31个2阶带通滤波器。图6(a)图6(b)图6(c)中,3种设计方法都是固定的20~20000 Hz频率范围,并且由于增益是人为设定的,滤波器的实际增益与目标增益有一定的误差,从而拟合的均衡器响应与目标增益误差较大。红圆圈为滤波器目标增益设置点,黑圆圈为实际增益设置点,可以看到这3种设计方法都有6个校正无效的点,并且拟合的最大增益误差分别为0.94 dB, 3.33 dB和1.14 dB。这既浪费了计算资源也不满足高保真音频要求。图6(d)中,由于加入了扬声器频响范围识别和自适应计算增益,因而均衡频率范围控制在50~14000 Hz,滤波器的实际增益与目标增益一致,拟合的最大增益误差仅为0.23 dB。这既节省了计算资源也满足了高保真音频要求。

    图 6  不同方法拟合的均衡器响应

    图7显示了不同方法拟合的均衡器响应比较情况。将4种设计方法得到的均衡器响应放在一起,可以明显看出,只有本文方法AG-APGE高精度拟合目标增益设置点。

    图 7  不同方法拟合的均衡滤波器响应比较

    图8显示了不同图示均衡器校正方法对全频段扬声器进行校正的情况。对采集的某个实测数据进行扬声器频率响应范围识别可以得到该扬声器为全频段扬声器,频率响应范围是50~14000 Hz。在50~14000 Hz这段频率范围进行自适应校正。由于车内的干扰以及音响设备本身的干扰,在全频段范围内白噪声会有不同程度的失真。图8(a)图8(b)图8(c)中,蓝色曲线为车载扬声器实际幅度响应,黑粗实线分别为3种设计方法得到的均衡器幅度响应,黑点划线为目标幅度响应,黑虚线为目标幅度响应的±1 dB误差边界线。将均衡器幅度响应去校正车载扬声器实际幅度响应得到校正后幅度响应,即图中的红实线。从图8(a)图8(b)图8(c)可以看到,3种设计方法校正后的幅度响应部分超出了±1 dB的误差容限,不满足高保真音频要求,且谱平坦度分别为0.994, 0.989和0.993。图8(d)中,经过自适应校正算法处理后,校正后的音频信号频谱曲线趋于平直,基本在目标响应的±1 dB误差范围内,达到较好的校正效果,基本满足高保真音频要求。且谱平坦度为0.996。

    图 8  不同方法的车载扬声器均衡

    图9显示了不同方法校正后的幅度响应对比情况。从中可以明显看出,本文方法比其他3种方法的校正效果好。

    图 9  不同方法的车载扬声器均衡比较

    表3显示了不同设计方法的指标对比情况。对采集的50个数据进行了分析统计,表中的数值是对50个数据求取算术平均值后的结果。ACGE, PGE, APGE和AG-APGE拟合的平均最大增益误差分别为1.08 dB, 2.65 dB, 1.07 dB和0.23 dB。只有本文方法AG-APGE满足高保真音频要求。信号平均原始谱平坦度为0.652,4种方法平均校正后信号谱平坦度分别为0.987, 0.980, 0.986和0.991。由于本文方法拟合目标增益精度最高,从而校正后信号谱平坦度最高。前3种设计方法由于是固定的均衡范围,因而其均衡滤波器个数分别为31, 62和31。而本文方法加入了扬声器频响范围识别,因而其平均均衡滤波器个数最少,为25个。

    表 3  不同方法指标对比
    指标ACGEPGEAPGEAG-APGE
    平均拟合最大误差(dB)1.082.651.070.23
    平均原始SFM0.6520.6520.6520.652
    平均校正后SFM0.9870.9800.9860.991
    平均均衡滤波器个数31623125
    平均设计所需时间0.57 ms1.22 s0.70 ms0.48 ms
    平均操作数
    加法124248124100
    乘法155249125101
    总和279497249201
    下载: 导出CSV 
    | 显示表格

    对于ACGE,每个频带滤波器有4次加法和4次乘法,31个这样的频带加上G0部分1次乘法以及31个频带之间的乘法,一共就有124次加法和155次乘法。对于PGE,每个滤波器有3次加法和4次乘法,62个这样的滤波器加上常数项以及它们之间的加法,一共有248次加法和249次乘法。对于APGE,每个频带滤波器有3次加法和4次乘法,31个这样的频带以及整个结构的加法,一共有124次加法和125次乘法。对于AG-APGE,由于平均使用了25个滤波器,因此一共有100次加法和101次乘法。因此本文方法的操作数最少。表3中设计时间为输入目标增益设置点后计算滤波器系数所需的时间。PGE方法设计滤波器步骤较复杂,因此平均设计所需时间最长,为1.22 s。本文方法平均设计所需时间和ACGE方法以及APGE方法为同一个数量级,但由于本文方法比这两种方法所用平均均衡滤波器个数少大约20%,因此本文方法设计图示均衡器速度最快。

    本文介绍了一种无重叠频带的自适应非对称并联图示均衡器校正方法并用于车载扬声器均衡,提出了一种自适应计算滤波器目标增益和识别车载扬声器频响范围的方法。图示均衡器是基于并联结构设计的,这意味着图示均衡器可以在图形处理器中实现,以提高效率。现有的经典图示均衡器校正方法拟合目标增益精度都很高,拟合误差都在1 dB以内。然而,在车载应用背景下,车载音响的声场会动态变化,导致滤波器拟合的增益与目标增益有一定的误差。并且每个车载扬声器有其特定的有效频响范围,而经典方法的固定均衡频率范围会造成计算资源的浪费。经过具体的仿真和对多种实验结果的评估,证明了所提出的方法在车载应用背景下比经典方法所用均衡滤波器个数平均减少大约20%,同时高精度拟合目标增益,且校正效果较好。因此,所提出的方法比经典方法更适用于车载音频均衡应用。

  • 图  1  图示均衡器串联组成示意图

    图  2  单个2阶滤波器结构框图

    图  3  图示均衡器并联组成示意图

    图  4  单个2阶滤波器结构框图

    图  5  车载音响图示均衡器自动校正算法应用系统框图

    图  6  不同方法拟合的均衡器响应

    图  7  不同方法拟合的均衡滤波器响应比较

    图  8  不同方法的车载扬声器均衡

    图  9  不同方法的车载扬声器均衡比较

    表  1  车载扬声器分类及其频响范围表

    车载扬声器分类常见频率响应范围
    低音20~150 Hz
    中音150 Hz~5 kHz
    高音5~20 kHz
    中低音20 Hz~5 kHz
    中高音110 Hz~20 kHz
    全频段50 Hz~20 kHz
    下载: 导出CSV

    表  2  1/3倍频程中心频率和带宽

    频带编号fc(Hz)fB(Hz)频带编号fc(Hz)fB(Hz)
    1209.317800373
    22511.7181000466.2
    331.514.7191250582.8
    44018.6201600745.9
    55023.3212000932.4
    66329.42225001166
    78037.32331501469
    810046.62440001865
    912558.32550002331
    1016074.62663002937
    1120093.22780003730
    12250116.628100004662
    13315146.929125005828
    14400186.530160007459
    15500233.131200009324
    16630293.7
    下载: 导出CSV

    表  3  不同方法指标对比

    指标ACGEPGEAPGEAG-APGE
    平均拟合最大误差(dB)1.082.651.070.23
    平均原始SFM0.6520.6520.6520.652
    平均校正后SFM0.9870.9800.9860.991
    平均均衡滤波器个数31623125
    平均设计所需时间0.57 ms1.22 s0.70 ms0.48 ms
    平均操作数
    加法124248124100
    乘法155249125101
    总和279497249201
    下载: 导出CSV
  • [1] VÄLIMÄKI V and REISS J D. All about audio equalization: Solutions and frontiers[J]. Applied Sciences, 2016, 6(5): 129. doi: 10.3390/app6050129
    [2] CECCHI S, CARINI A, and SPORS S. Room response equalization—a review[J]. Applied Sciences, 2018, 8(1): 16. doi: 10.3390/app8010016
    [3] PEPE G, GABRIELLI L, SQUARTINI S, et al. Evolutionary tuning of filters coefficients for binaural audio equalization[J]. Applied Acoustics, 2020, 163: 107204. doi: 10.1016/j.apacoust.2019.107204
    [4] DAGAR A, NITISH S S, and HEGDE R. Joint adaptive impulse response estimation and inverse filtering for enhancing in-car audio[C]. 2018 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Calgary, Canada, 2018: 416–420. doi: 10.1109/ICASSP.2018.8462329.
    [5] VAIRETTI G, DE SENA E, CATRYSSE M, et al. An automatic design procedure for low-order IIR parametric equalizers[J]. Journal of the Audio Engineering Society, 2018, 66(11): 935–952. doi: 10.17743/jaes.2018.0049
    [6] RÄMÖ J and VÄLIMÄKI V. Graphic delay equalizer[C]. ICASSP 2019-2019 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Brighton, England, 2019: 8018–8022. doi: 10.1109/icassp.2019.8682949.
    [7] PRINCE S and KUMAR K R S. A novel Nth-order IIR filter-based graphic equalizer optimized through genetic algorithm for computing filter order[J]. Soft Computing, 2019, 23(8): 2683–2691. doi: 10.1007/s00500-018-3640-9
    [8] PEPE G, GABRIELLI L, SQUARTINI S, et al. Gravitational search algorithm for IIR filter-based audio equalization[C]. 2020 28th European Signal Processing Conference (EUSIPCO), Amsterdam, Holland, 2021: 496–500. doi: 10.23919/Eusipco47968.2020.9287421.
    [9] PEPE G, GABRIELLI L, SQUARTINI S, et al. Designing audio equalization filters by deep neural networks[J]. Applied Sciences, 2020, 10(7): 2483. doi: 10.3390/app10072483
    [10] RÄMÖ J, LISKI J, and VÄLIMÄKI V. Third-octave and bark graphic-equalizer design with symmetric band filters[J]. Applied Sciences, 2020, 10(4): 1222. doi: 10.3390/app10041222
    [11] RÄMÖ J and VÄLIMÄKI V. Neural third-octave graphic equalizer[C]. Proceedings of the 22nd International Conference on Digital Audio Effects (DAFx-19), Birmingham, UK, 2019: 2–6.
    [12] VÄLIMÄKI V and RÄMÖ J. Neurally controlled graphic equalizer[J]. IEEE/ACM Transactions on Audio, Speech, and Language Processing, 2019, 27(12): 2140–2149. doi: 10.1109/taslp.2019.2935809
    [13] RAMÍREZ M A M and REISS J D. End-to-end equalization with convolutional neural networks[C]. 21st International Conference on Digital Audio Effects (DAFx-18), Aveiro, Portugal, 2018: 296–303.
    [14] RÄMÖ J and VÄLIMÄKI V. Optimizing a high-order graphic equalizer for audio processing[J]. IEEE Signal Processing Letters, 2014, 21(3): 301–305. doi: 10.1109/lsp.2014.2301557
    [15] LISKI J, RÄMÖ J, and VÄLIMÄKI V. Graphic equalizer design with symmetric biquad filters[C]. 2019 IEEE Workshop on Applications of Signal Processing to Audio and Acoustics (WASPAA), New Paltz, USA, 2019: 55–59. doi: 10.1109/WASPAA.2019.8937168.
    [16] VÄLIMÄKI V and LISKI J. Accurate cascade graphic equalizer[J]. IEEE Signal Processing Letters, 2017, 24(2): 176–180. doi: 10.1109/lsp.2016.2645280
    [17] LISKI J and VÄLIMÄKI V. The quest for the best graphic equalizer[C]. Proceedings of the 20th International Conference on Digital Audio Effects (DAFx-17), Edinburgh, UK, 2017: 95–102.
    [18] LISKI J, BANK B, SMITH J O, et al. Converting series biquad filters into delayed parallel form: Application to graphic equalizers[J]. IEEE Transactions on Signal Processing, 2019, 67(14): 3785–3795. doi: 10.1109/TSP.2019.2919419
    [19] TAO Hao, FANG Yong, LIU Huaping, et al. Optimization method for a high-precision graphic equalizer with lower order filter[C]. 2018 14th IEEE International Conference on Signal Processing (ICSP), Beijing, China, 2018: 22–26. doi: 10.1109/ICSP.2018.8652302.
    [20] BANK B. Converting infinite impulse response filters to parallel form [tips & tricks][J]. IEEE Signal Processing Magazine, 2018, 35(3): 124–130. doi: 10.1109/msp.2018.2805358
    [21] BANK B, BELLOCH J A, and VÄLIMÄKI V. Efficient design of a parallel graphic equalizer[J]. Journal of the Audio Engineering Society, 2017, 65(10): 817–825. doi: 10.17743/jaes.2017.0029
    [22] RÄMÖ J, VÄLIMÄKI V, and BANK B. High-precision parallel graphic equalizer[J]. IEEE/ACM Transactions on Audio, Speech, and Language Processing, 2014, 22(12): 1894–1904. doi: 10.1109/taslp.2014.2354241
    [23] CHEN Zhe, GENG Guosheng, YIN Fuliang, et al. A pre-distortion based design method for digital audio graphic equalizer[J]. Digital Signal Processing, 2014, 25: 296–302. doi: 10.1016/j.dsp.2013.11.007
    [24] ORFANIDIS S J. Introduction to Signal Processing[M]. Upper Saddle River: Prentice Hall, 1995.
    [25] KARJALAINEN M, PIIRILÄ E, JÄRVINEN A, et al. Comparison of loudspeaker equalization methods based on DSP techniques[J]. Journal of the Audio Engineering Society, 1999, 47(1/2): 14–31.
  • 期刊类型引用(1)

    1. 王丹. 多通道广播音频信号处理中的相位校正技术. 电声技术. 2024(10): 117-119 . 百度学术

    其他类型引用(0)

  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  850
  • HTML全文浏览量:  548
  • PDF下载量:  62
  • 被引次数: 1
出版历程
  • 收稿日期:  2021-03-15
  • 修回日期:  2021-10-20
  • 网络出版日期:  2021-11-04
  • 刊出日期:  2022-05-25

目录

/

返回文章
返回