CN102447445A - 一种音频参数均衡的方法和音频参数均衡器 - Google Patents
一种音频参数均衡的方法和音频参数均衡器 Download PDFInfo
- Publication number
- CN102447445A CN102447445A CN2010105020696A CN201010502069A CN102447445A CN 102447445 A CN102447445 A CN 102447445A CN 2010105020696 A CN2010105020696 A CN 2010105020696A CN 201010502069 A CN201010502069 A CN 201010502069A CN 102447445 A CN102447445 A CN 102447445A
- Authority
- CN
- China
- Prior art keywords
- spectral coefficient
- mdct
- frequency band
- frequency
- gain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Tone Control, Compression And Expansion, Limiting Amplitude (AREA)
Abstract
本发明提供了一种音频参数均衡的方法和音频参数均衡器,所述方法包括:通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。通过本发明能够减少参数均衡处理中的运算量,且实现简单。
Description
技术领域
本发明涉及音频处理技术领域,特别是涉及一种音频参数均衡的方法和一种音频参数均衡器。
背景技术
音频均衡器通过调整不同频带信号的增益实现摇滚,流行等各种音效。按照调整参数范围的不同,均衡器分为图形均衡器和参数均衡器两种。对于图形均衡器,是一些并联的带通滤波器,它们接收相同的输入信号,但每个滤波器各频带的中心频率和带宽为设定值,用户只能通过调整各频带的增益值来控制对各频带信号增强或削弱的幅度。另一种是参数均衡器也是一些带通滤波器的组合,除了频带增益外,它还允许用户调整各频带的中心频率和带宽,具有最大的灵活性,图形均衡器的控制细致程度不如参数均衡器,可以说图形均衡器是参数均衡器的一个子集。
基于无限脉冲响应(IIR,Infinite Impulse Response)或有限冲激响应(FIR,Finite Impulse Response)带通滤波器的参数均衡器的设计思想一般是根据用户设定的中心频率和带宽参数,同步计算和更新各个带通滤波器的系数。但是带通滤波器的设计算法运算量很大。当用户尝试性的连续调整均衡器参数以获得最佳效果时,如果运算速度不够,滤波器系数的计算和更新过程不能和用户调整过程同步,最后设定的参数可能并不符合用户的要求。此外,这种方法缺乏控制机制,如用户将带宽参数调整到较小值时,由滤波器设计算法得到的滤波器阶数可能相当高,则计算量比较大,设计过程中改变参数时都要重新计算,无法实时实现。
总之,需要本领域技术人员迫切解决的一个技术问题就是:如何能够一种音频参数均衡技术,能够减少参数均衡处理中的运算量,且实现简单。
发明内容
本发明所要解决的技术问题是提供一种音频参数均衡的方法和音频参数均衡器,能够减少参数均衡处理中的运算量,且实现简单。
为了解决上述问题,本发明公开了一种音频参数均衡的方法,包括:
通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;
根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;
通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。
进一步,通过MDCT变换获得N个MDCT谱系数;其中,N为MDCT变换的块长;
优选的,所述根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,包括:
通过将第b个频带的中心频率Fc[b]除以滤波带宽FB,获得与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],0≤k[b]≤N-1;
通过将第b个频带的带宽BW(b)除以滤波带宽FB,获得与第b个频带相对应的MDCT谱系数的个数num[b];
依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
其中,b=1、2、3....m;m为用户设置的频带个数。
优选的,第b个频带对应的MDCT谱系数的序号k的范围为:[k[b]-l[b],k[b]+r[b]];l[b]、r[b]分别为k[b]号MDCT谱系数左、右两侧谱系数的个数;
当num[b]为奇数时,l[b]=r[b];当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|<|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]大1;当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|>|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]小1;其中,l[b]+r[b]+1=num[b],Fca(k[b]-1)为第k[b]-1号MDCT谱系数对应的窄带信号中心频率,Fca(k[b]+1)为第k[b]+1号MDCT谱系数对应的窄带信号中心频率。
优选的,所述根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整,包括:
依据增益gain[k]分别对各个谱系数进行增益调整。
优选的,采用线性插值法,并根据第k号MDCT谱系数对应的谱线到第k[b]号MDCT谱系数对应的谱线的距离确定所述wb[k],其中:当k<k[b]-l[b]或者k>k[b]+r[b]时,wb[k]取值为零;当k=k[b]时,wb[k]取值为一;当k[b]-l[b]≤k<k[b]时,wb[k]在[0,1)范围内正比例取值;当k[b]<k≤k[b]+r[b]时,wb[k]在[0,1)范围内反比例取值。
此外,本发明还公开了一种音频参数均衡器,包括:
MDCT变换模块,用于通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
频带划分模块,用于根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;
增益调整模块,用于根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;
IMDCT变换模块,用于通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。
进一步,MDCT变换模块通过MDCT变换获得N个MDCT谱系数;其中,N为MDCT变换的块长;
优选的,所述频带划分模块包括:
中心频率序号获取子模块,用于通过将第b个频带的中心频率Fc[b]除以滤波带宽FB,获得与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],0≤k[b]≤N-1;
谱系数个数获取子模块,通过将第b个频带的带宽BW(b)除以滤波带宽FB,获得与第b个频带相对应的MDCT谱系数的个数num[b];
谱系数范围确定子模块,用于依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
其中,b=1、2、3....m;m为用户设置的频带个数。
优选的,第b个频带对应的MDCT谱系数的序号k的范围为:[k[b]-l[b],k[b]+r[b]];l[b]、r[b]分别对应k[b]号谱系数左、右两侧谱系数的个数;
当num[b]为奇数时,l[b]=r[b];当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|<|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]大1;当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|>|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]小1;其中,l[b]+r[b]+1=num[b],Fca(k[b]-1)为第k[b]-1号MDCT谱系数对应的窄带信号中心频率,Fca(k[b]+1)为第k[b]+1号MDCT谱系数对应的窄带信号中心频率。
优选的,所述根增益调整模块包括:
增益调整子模块,用于依据增益gain[k]分别对各个谱系数进行增益调整。
优选的,采用线性插值法,并根据第k号MDCT谱系数对应的谱线到第k[b]号MDCT谱系数对应的谱线的距离确定所述wb[k],其中,当k<k[b]-l[b]或者k>k[b]+r[b]时,wb[k]取值为零;当k=k[b]时,wb[k]取值为一;当k[b]-l[b]≤k<k[b]时,wb[k]在[0,1)范围内正比例取值;当k[b]<k≤k[b]+r[b]时,wb[k]在[0,1)范围内反比例取值。
与现有技术相比,本发明具有以下优点:
本发明首先采用MDCT变换将输入时域音频信号转换成频域MDCT谱系数,然后根据用户设定的中心频率和带宽定位各频带对应的MDCT谱系数的范围,再根据频带增益和各系数距离中心频率的距离确定各MDCT谱系数的幅度增益值,对MDCT谱系数进行增益调整后进行IMDCT变换得到输出时域信号。在本发明中,由于各MDCT谱系数可以独立的调整幅度,MDCT谱系数作为一个集合对应一个带通滤波器。某带通滤波器中心频率和带宽的改变,由对应分组中选择的MDCT谱系数的序号和个数的改变来模拟,因此该方法具有很强的灵活性。相比传统方法,用户设定不同的中心频率和带宽时,不需要进行带通滤波器的设计和更新,只需相应调整各MDCT谱系数的增益值。
此外,不同于逐点运算的传统方案,本发明可以采用MDCT变换一次处理多个样本点,由于MDCT变换具有成熟的快速算法,因此,本方案具有运算量小的优点。进一步,MDCT变换的运算量只和MDCT变换长度N有关,因此无论用户如何调整参数,本方案的运算复杂度基本不变。
附图说明
图1是本发明一种音频参数均衡的方法实施例一的流程图;
图2是一种音频信号数据MDCT变换的示意图;
图3是一种音频信号数据IMDCT变换的示意图;
图4是本发明一种音频参数均衡的方法实施例二的流程图;
图5是本发明一种音频参数均衡器实施例的流程图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明的核心构思之一在于:在实现参数均衡器时,当用户调整频带中心频率和带宽时,不需要计算和更新带通滤波器系数。而是采用一种对多个窄带滤波器的输出(MDCT谱系数)进行分组模拟带通滤波器的方法,即根据用户设定的中心频率和带宽定位各频带对应的MDCT谱系数的范围,再根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整,最后通过IMDCT变换得到输出时域信号。
参照图1,示出了本发明一种音频参数均衡的方法实施例一的流程图,包括:
步骤11,通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
改进的离散余弦变换(Modified Discrete Cosine Transform,MDCT)是一种与傅立叶变换相关的变换,MDCT变换的计算公式为:
上式中,X[k]是通过MDCT变换得到的频域MDCT谱系数,通过MDCT变换获得N个MDCT谱系数,N为MDCT变换的块长,为输入的时域音频信号x[n]的加窗数据,n,k分别为时域音频信号、频域MDCT谱系数的序号。其中,
本发明实施例中使用的窗函数为正弦窗h[n],
h[n]=sin[π(n+1/2)/2N] [3]
一次MDCT变换的输入是前后两个信号数据块拼接而成的2N点加窗数据,相邻两次变换的输入数据会共用一个数据块的N点数据,如图2所示,为一种音频信号数据MDCT变换的示意图,其示出了两次变换输入数据的重叠,图中各数据块均含有N点数据。
MDCT变换本质上是一种基于时域混叠消除的滤波器组。各MDCT谱系数X[k]均可看成一个窄带带通滤波器的输出。设输入信号x[n]的采样率为Fs,则各MDCT谱系数对应的窄带信号带宽(对应MDCT变换的各窄带滤波器带宽)均为:
其中,Fs为输入时域音频信号的采样率。各MDCT谱系数对应的窄带信号中心频率为:
由[4]可知,在采样率Fs一定时,MDCT变换的块长N越大,FB越小。一般要求FB<30Hz,较小的带宽有利于准确的实现用户设定的各种参数设置。此外,为了符合快速算法的需要,N/2应取为2的整数次幂。
步骤12,根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;
参数均衡器可供调节的频带(也称频段)个数一般在5~30之间,而最常见是10段均衡器。显然,对于图形均衡器,用户仅能改变各频段的增益值。传统方案中,一旦用户改变各频段的中心频率和带宽,就需要重新设计带通滤波器,更新滤波器系数。而本方案中,由于N个MDCT谱系数分别对应N个带宽为FB的窄带滤波器的输出,可以采用对X[k]进行频带划分的方法,即对X[k]进行分组,一个组中的若干MDCT谱系数作为一个集合对应一个带通滤波器。某带通滤波器中心频率和带宽的改变,由对应分组中选择的MDCT谱系数的序号和个数的改变来模拟。具体的,所述步骤12包括如下子步骤:
子步骤121,计算与第b个频带的中心频率相对应的MDCT谱系数的序号k[b];
由[5]可知,与中心频率最对应的MDCT谱系数的序号的计算方法为:将第b个频带的中心频率Fc[b]除以滤波带宽FB,考虑到除得的结果可能为小数,在具体实施中,可以修正为:
上式中,b=1、2、3....m,m为用户设置的频带个数,k[b]的取值范围是0~N-1,INT为取整运算,Fc[b]为用户设置的第b个频带的中心频率,FB为由[4]计算得到的MDCT变换的滤波器带宽。
子步骤122,计算与第b个频带相对应的MDCT谱系数的个数num[b];
-般的,通过将第b个频带的带宽BW(b)除以滤波带宽FB得到num[b],考虑到除得的结果可能为小数,在具体实施中,可以通过四舍五入修正为:
上式中,b的取值范围是1~m,INT为取整运算,BW(b)为第b个频带的带宽。
子步骤123,依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
具体的,第b个频带对应的MDCT谱系数的序号范围为:[k[b]-l[b],k[b]+r[b]];l[b]、r[b]分别对应k[b]号频谱系数左、右两侧谱系数的个数,第b个频段包含的总的谱系数个数满足num[b]=r[b]+l[b]+1。左、右两侧谱系数的取值具体包括如下三种情况:
1、当num[b]为奇数时,第b个频带对应的各个MDCT谱系数的范围为:以第k[b]号MDCT谱系数为中心,左、右两边分别取相同数目个MDCT谱系数。则,
但是,当num[b]为偶数时,以序号为k[b]的MDCT谱系数为中心,-侧的谱系数应比另一侧多一个,则具体由与k[b]号最邻近的两个MDCT谱系数对应的中心频率与Fc[b]的接近程度决定。设定Fca(k[b]-1)为第k[b]-1号MDCT谱系数对应的窄带信号中心频率,Fca(k[b]+1)为第k[b]+1号MDCT谱系数对应的窄带信号中心频率,均由公式[5]计算得到。
2、当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|<|Fca(k[b]+1)-Fc(b)|时,由于左边的第一个谱系数更靠近中心频率,因此左边多一个谱系数,左边最远端的谱系数相对来说离中心频率也会近些。第b个频带对应的各个MDCT谱系数的范围为:以第k[b]号MDCT谱系数为中心,左边比右边多取一个MDCT谱系数。则,
3、当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|>|Fca(k[b]+1)-Fc(b)|时,类似的,第b个频带对应的各个MDCT谱系数的范围为:以第k[b]号MDCT谱系数为中心,左边比右边少取一个MDCT谱系数。则,
频带分组划分的功能是由各频段中心频率Fc[b]确定其对应的MDCT谱系数的序号k[b]。再结合带宽参数BW(b)确定该频段对应的MDCT谱系数的范围。显然,本发明实施例可以简便的实现各种形式的带宽设置,如非对称形式的,也允许频带之间有重叠部分。
步骤13,根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;
在确定了各频段对应的MDCT谱系数的范围后,需要根据用户设定的各频段的增益值确定对应的MDCT谱系数分组中各谱系数的幅度增益值。基本的原则是增益值从k[b]号MDCT谱系数开始,向两侧逐渐衰减。具体的,所述步骤13包括如下子步骤:
子步骤131,依据用户设置的各频带增益,依次计算各MDCT谱系数的增益gain[k];
其中,g[b]为用户设置的第b个频带的增益(单位取db);wb[k]为与第b个频带增益对应的第k号MDCT谱系数的加权值,k=0~N-1。wb[k]与各谱线到中心频率谱线的距离成反比例取值,即:距离越小,wb[k]越大,取值越接近1;距离越大,wb[k]越小,取值越接近0。上式的意义是统计各频带增益参数对某一谱系数的影响。
wb[k]为处于0~1之间的数值,如果第k号谱系数在第b个频带对应的谱系数范围(也即[k[b]-l[b],k[b]+r[b]])之外,则wb[k]=0;如果第k号谱系数在此范围内,则
在本发明的-个实施例中,可以采用采用线性插值法,并根据第k号MDCT谱系数对应的谱线到中心频率对应的谱线的距离确定wb[k]。基本原则是各频带中心频率谱线处(对应k[b]号谱系数)的加权值wb[k]=1,各频带两端点谱线处的加权值wb[k]=0,两端点与中心频率间的谱系数采用的加权值根据各谱线到中心频率谱线的距离用线性插值方法得到。wb[k]的取值如下:
也就是说,当MDCT谱系数的序号k在第b个频段包含MDCT谱系数范围之外时,即k<k[b]-l[b]或者k>k[b]+r[b],没有增益,则wb[k]=0;当k为中心频率对应的MDCT谱系数时,即k=k[b],增益应该最大,wb[k]=1;当k[b]-l[b]≤k<k[b]时,根据k取值线性递增,wb[k]在[0,1)范围内正比例取值;当k[b]<k≤k[b]+r[b]时,根据k取值线性递增,wb[k]在[0,1)范围内反比例取值。
子步骤132,依据增益gain[k]分别对各个谱系数进行增益调整。
结合[8]、[9],由用户设定的各频带增益值g[b]可计算得各谱系数的增益db值gain[k]。gain[k]与幅度增益因子的关系为。
gain_fac[k]=10gain[k]/20
各谱系数的增益调整即是乘上幅度增益因子,即
Y[k]=X[k]*gain_fac[k]
其中,Y[k]为增益调整后的MDCT谱系数。
步骤14,根据用户设置的参数,对输入信号的MDCT谱系数进行相应幅度调整后,对谱系数Y[k]进行改进的离散余弦反变换(Inverse ModifiedDiscrete Cosine Transform,IMDCT)变换即得到输出数字信号y[n]。
IMDCT变换的计算公式为:
IMDCT变换的运算过程含有加窗和重叠加的过程。
合成窗g[n]需与步骤11中的分析窗h[n]对应,对于正弦窗,
g[n]=h[n] (n=0~2N-1)。
而后将前一个合成窗的后N点与当前合成窗的前N点相加得到N点输出,如图3所示,为一种音频信号数据IMDCT变换的示意图。通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号,即实现了对音频参数的均衡。
下面,举一个具体的例子对本发明进行详细说明。在本发明实施例中,输入信号x[n]的采样率Fs=44100Hz,N取典型值为1024,依据公式[4],得到FB≈22Hz。具体的,参照图4,示出了本发明一种音频参数均衡的方法实施例二的流程图,包括:
步骤41,通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
步骤42,计算与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],b=1,2,....m,m为用户设置的频带个数;
由于均衡器可供调节的频段个数一般在5~30之间,而最常见是10段均衡器,在本发明实施例中,假设用户设定10个频带,各个频带的中心频率和带宽参数的设置如下表所示:
表1
表1中,第一行为中心频率,第二行为频段带宽。这里带宽理解为各频段调节的作用范围,且以中心频率对称。如第7号带的中心频率为6kHz,带宽为4kHz,则该带的控制范围为[4kHz~8kHz],该带增益的调整不影响这个范围之外的频率成分,而这个范围内的频率成分,距离6kHz越远,该带增益调整带来的影响越小。
具体的,计算与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],此处以7号频带(band)为例,由公式[6]得:
即对第7个band的增益调整对应着以第272根MDCT谱线为中心的若干MDCT谱系数的幅度的调整。
步骤43,计算与第b个频带相对应的MDCT谱系数的个数num[b];
其中,b分别取1,2,....m;
例如,针对第7个band,由公式[7]得:
即与第7个频带相对应的MDCT谱系数的个数为182。
步骤44,依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
其中,b分别取1,2,....m。在本实施例中,num[b]=182,为偶数,则以k[b]为中心,一侧的谱系数应比另一侧多一个,具体由与k[b]最邻近的两个MDCT谱系数对应的中心频率与Fc[b]的接近程度决定。
本实施例中Fc[b]=6000Hz,k[b]=272。由公式[5],得到Fca(271)=5973Hz,Fca(273)=6017Hz。显然右侧的Fca(273)更接近Fc[b],因此第7号band对应的MDCT谱系数的序号范围是[272-90,272+91],即右边比左边多一个。最终确定下来与用户设定值对应的第7号band包含的MDCT谱系数的序号范围是[182,363]。也即对7号频段的增益调整将转化为对序号为[182,363]这182个MDCT谱系数的调整。同理,也可得其它频带MDCT谱系数的范围。
步骤45,依次计算第k号MDCT谱系数的幅度增益gain[k],k=0~N-1;
步骤46,依据增益gain[k]获得各个谱系数幅度增益因子,并分别对各个谱系数进行增益调整;
步骤47,通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。
在本发明实施例中,重点针对用户设置的第7个频带进行说明,其他频带的增益调整与上述类似,此处不再赘述,具体方法可参考第7个频带。
参照图5,示出了本发明一种音频参数均衡器实施例的结构图,包括:
MDCT变换模块51,用于通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
频带划分模块52,用于根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;
增益调整模块53,用于根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;
IMDCT变换模块54,用于通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。
进一步,MDCT变换模块通过MDCT变换获得N个MDCT谱系数;其中,N为MDCT变换的块长;
则各MDCT谱系数对应的窄带信号带宽其中,Fs为输入时域音频信号的采样率。
在本发明的优选实施例中,所述频带划分模块52包括:
中心频率序号获取子模块521,用于通过将第b个频带的中心频率Fc[b]除以滤波带宽FB,获得与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],0≤k[b]≤N-1;
谱系数个数获取子模块522,通过将第b个频带的带宽BW(b)除以滤波带宽FB,获得与第b个频带相对应的MDCT谱系数的个数num[b];
谱系数范围确定子模块523,用于依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
其中,b=1、2、3....m;m为用户设置的频带个数。
进一步,第b个频带对应的MDCT谱系数的序号k的范围为:[k[b]-l[b],k[b]+r[b]];l[b]、r[b]分别对应k[b]号谱系数左、右两侧谱系数的个数;
当num[b]为奇数时,第b个频带对应的各个MDCT谱系数的范围为:以第k[b]号MDCT谱系数为中心,左、右两边分别取相同数目个MDCT谱系数,l[b]=r[b];
当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|<|Fca(k[b]+1)-Fc(b)|时,第b个频带对应的各个MDCT谱系数的范围为:以第k[b]号MDCT谱系数为中心,左边比右边多取一个MDCT谱系数,
当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|>|Fca(k[b]+1)-Fc(b)|时,第b个频带对应的各个MDCT谱系数的范围为:以第k[b]号MDCT谱系数为中心,左边比右边少取一个MDCT谱系数,
其中,l[b]+r[b]+1=num[b],Fca(k[b]-1)为第k[b]-1号MDCT谱系数对应的窄带信号中心频率,Fca(k[b]+1)为第k[b]+1号MDCT谱系数对应的窄带信号中心频率。
在本发明的优选实施例中,所述增益调整模块53包括:
增益调整子模块532,用于依据增益gain[k]分别对各个谱系数进行增益调整。
优选的,采用线性插值法,并根据第k号MDCT谱系数对应的谱线到第k[b]号MDCT谱系数对应的谱线的距离确定所述wb[k],其中,
当k<k[b]-l[b]或者k>k[b]+r[b]时,wb[k]取值为零;
当k=k[b]时,wb[k]取值为一;
当k[b]-l[b]≤k<k[b]时,wb[k]在[0,1)范围内正比例取值;
当k[b]<k≤k[b]+r[b]时,wb[k]在[0,1)范围内反比例取值。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上对本发明所提供的一种音频参数均衡的方法和音频参数均衡器,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (12)
1.一种音频参数均衡的方法,其特征在于,包括:
通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;
根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;
通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。
3.如权利要求2所述的方法,其特征在于,所述根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,包括:
通过将第b个频带的中心频率Fc[b]除以滤波带宽FB,获得与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],0≤k[b]≤N-1;
通过将第b个频带的带宽BW(b)除以滤波带宽FB,获得与第b个频带相对应的MDCT谱系数的个数mum[b];
依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
其中,b=1、2、3....m;m为用户设置的频带个数。
4.如权利要求3所述的方法,其特征在于,
第b个频带对应的MDCT谱系数的序号k的范围为:[k[b]-l[b],k[b]+r[b]];l[b]、r[b]分别为k[b]号MDCT谱系数左、右两侧谱系数的个数;
当num[b]为奇数时,l[b]=r[b];
当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|<|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]大1;
当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|>|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]小1;
其中,l[b]+r[b]+1=num[b],Fca(k[b]-1)为第k[b]-1号MDCT谱系数对应的窄带信号中心频率,Fca(k[b]+1)为第k[b]+1号MDCT谱系数对应的窄带信号中心频率。
6.如权利要求5所述的方法,其特征在于,
采用线性插值法,并根据第k号MDCT谱系数对应的谱线到第k[b]号MDCT谱系数对应的谱线的距离确定所述wb[k],其中:
当k<k[b]-l[b]或者k>k[b]+r[b]时,wb[k]取值为零;
当k=k[b]时,wb[k]取值为一;
当k[b]-l[b]≤k<k[b]时,wb[k]在[0,1)范围内正比例取值;
当k[b]<k≤k[b]+r[b]时,wb[k]在[0,1)范围内反比例取值。
7.一种音频参数均衡器,其特征在于,包括:
MDCT变换模块,用于通过MDCT变换将输入的时域音频信号转换成频域MDCT谱系数;
频带划分模块,用于根据用户设置的各频带的中心频率和带宽确定各频带对应的MDCT谱系数的范围,对MDCT谱系数进行频带分组划分;
增益调整模块,用于根据用户设置的各频带的增益对各频带内的MDCT谱系数进行相应增益调整;
IMDCT变换模块,用于通过IMDCT将增益调整后的MDCT谱系数转换为时域音频信号。
9.如权利要求8所述的音频参数均衡器,其特征在于,所述频带划分模块包括:
中心频率序号获取子模块,用于通过将第b个频带的中心频率Fc[b]除以滤波带宽FB,获得与第b个频带的中心频率相对应的MDCT谱系数的序号k[b],0≤k[b]≤N-1;
谱系数个数获取子模块,通过将第b个频带的带宽BW(b)除以滤波带宽FB,获得与第b个频带相对应的MDCT谱系数的个数num[b];
谱系数范围确定子模块,用于依据k[b]和num[b]确定第b个频带对应的MDCT谱系数的范围;
其中,b=1、2、3....m;m为用户设置的频带个数。
10.如权利要求9所述的音频参数均衡器,其特征在于,
第b个频带对应的MDCT谱系数的序号k的范围为:[k[b]-l[b],k[b]+r[b]];l[b]、r[b]分别对应k[b]号谱系数左、右两侧谱系数的个数;
当num[b]为奇数时,l[b]=r[b];
当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|<|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]大1;
当num[b]为偶数,且|Fca(k[b]-1)-Fc(b)|>|Fca(k[b]+1)-Fc(b)|时,l[b]比r[b]小1;
其中,l[b]+r[b]+1=num[b],Fca(k[b]-1)为第k[b]-1号MDCT谱系数对应的窄带信号中心频率,Fca(k[b]+1)为第k[b]+1号MDCT谱系数对应的窄带信号中心频率。
12.如权利要求11所述的音频参数均衡器,其特征在于,
采用线性插值法,并根据第k号MDCT谱系数对应的谱线到第k[b]号MDCT谱系数对应的谱线的距离确定所述wb[k],其中,
当k<k[b]-l[b]或者k>k[b]+r[b]时,wb[k]取值为零;
当k=k[b]时,wb[k]取值为一;
当k[b]-l[b]≤k<k[b]时,wb[k]在[0,1)范围内正比例取值;
当k[b]<k≤k[b]+r[b]时,wb[k]在[0,1)范围内反比例取值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010502069.6A CN102447445B (zh) | 2010-09-30 | 2010-09-30 | 一种音频参数均衡的方法和音频参数均衡器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010502069.6A CN102447445B (zh) | 2010-09-30 | 2010-09-30 | 一种音频参数均衡的方法和音频参数均衡器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102447445A true CN102447445A (zh) | 2012-05-09 |
CN102447445B CN102447445B (zh) | 2016-04-20 |
Family
ID=46009596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010502069.6A Active CN102447445B (zh) | 2010-09-30 | 2010-09-30 | 一种音频参数均衡的方法和音频参数均衡器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102447445B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105720937A (zh) * | 2014-12-01 | 2016-06-29 | 宏达国际电子股份有限公司 | 电子装置和声音信号的分析与播放方法 |
CN106559056A (zh) * | 2015-09-25 | 2017-04-05 | 电信科学技术研究院 | 参数处理方法、音频信号播放方法、装置及音频均衡器 |
CN106973024A (zh) * | 2017-03-23 | 2017-07-21 | 电子科技大学 | 60GHz毫米波通信系统中均衡器的低复杂度的实现方法 |
CN110971213A (zh) * | 2019-12-11 | 2020-04-07 | 易兆微电子(杭州)有限公司 | 一种10段参数均衡器 |
CN111045634A (zh) * | 2018-10-12 | 2020-04-21 | 北京微播视界科技有限公司 | 音频处理方法和装置 |
US10693430B2 (en) | 2018-08-01 | 2020-06-23 | Realtek Semiconductor Corp. | Audio signal processing method and audio equalizer |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1227024A (zh) * | 1996-07-31 | 1999-08-25 | 夸尔柯姆股份有限公司 | 将用户选择的频率响应模式施加到提供给蜂窝式电话扬声器的音频信号上的方法和设备 |
US6341166B1 (en) * | 1997-03-12 | 2002-01-22 | Lsi Logic Corporation | Automatic correction of power spectral balance in audio source material |
US20110096942A1 (en) * | 2009-10-23 | 2011-04-28 | Broadcom Corporation | Noise suppression system and method |
-
2010
- 2010-09-30 CN CN201010502069.6A patent/CN102447445B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1227024A (zh) * | 1996-07-31 | 1999-08-25 | 夸尔柯姆股份有限公司 | 将用户选择的频率响应模式施加到提供给蜂窝式电话扬声器的音频信号上的方法和设备 |
US6341166B1 (en) * | 1997-03-12 | 2002-01-22 | Lsi Logic Corporation | Automatic correction of power spectral balance in audio source material |
US20110096942A1 (en) * | 2009-10-23 | 2011-04-28 | Broadcom Corporation | Noise suppression system and method |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105720937A (zh) * | 2014-12-01 | 2016-06-29 | 宏达国际电子股份有限公司 | 电子装置和声音信号的分析与播放方法 |
CN106559056A (zh) * | 2015-09-25 | 2017-04-05 | 电信科学技术研究院 | 参数处理方法、音频信号播放方法、装置及音频均衡器 |
CN106559056B (zh) * | 2015-09-25 | 2019-05-07 | 电信科学技术研究院 | 参数处理方法、音频信号播放方法、装置及音频均衡器 |
US10326422B2 (en) | 2015-09-25 | 2019-06-18 | China Academy Of Telecommunications Technology | Parameter processing method, audio signal playing method and device, and audio equalizer |
CN106973024A (zh) * | 2017-03-23 | 2017-07-21 | 电子科技大学 | 60GHz毫米波通信系统中均衡器的低复杂度的实现方法 |
CN106973024B (zh) * | 2017-03-23 | 2020-01-03 | 电子科技大学 | 60GHz毫米波通信系统中均衡器的低复杂度的实现方法 |
US10693430B2 (en) | 2018-08-01 | 2020-06-23 | Realtek Semiconductor Corp. | Audio signal processing method and audio equalizer |
CN111045634A (zh) * | 2018-10-12 | 2020-04-21 | 北京微播视界科技有限公司 | 音频处理方法和装置 |
CN111045634B (zh) * | 2018-10-12 | 2023-07-07 | 北京微播视界科技有限公司 | 音频处理方法和装置 |
CN110971213A (zh) * | 2019-12-11 | 2020-04-07 | 易兆微电子(杭州)有限公司 | 一种10段参数均衡器 |
CN110971213B (zh) * | 2019-12-11 | 2020-08-14 | 易兆微电子(杭州)股份有限公司 | 一种10段参数均衡器 |
Also Published As
Publication number | Publication date |
---|---|
CN102447445B (zh) | 2016-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102447445A (zh) | 一种音频参数均衡的方法和音频参数均衡器 | |
EP2122489B1 (en) | Frequency-warped audio equalizer | |
CN101405791B (zh) | 用于产生音频子带值的装置和方法以及用于产生时域音频采样的装置和方法 | |
CN101241150B (zh) | 信号处理装置、方法以及信号生成方法 | |
CA2777657C (en) | Reverberator and method for reverberating an audio signal | |
Goodwin | Residual modeling in music analysis-synthesis | |
CN104851429B (zh) | 叉积增强的基于子带块的谐波换位 | |
JP3528258B2 (ja) | 符号化音声信号の復号化方法及び装置 | |
CN102132340A (zh) | 参数立体声转换系统和方法 | |
CN102419981A (zh) | 音频信号时间尺度和频率尺度缩放处理方法及设备 | |
CN105745706A (zh) | 用于扩展频带的装置、方法和程序 | |
CN103137136A (zh) | 声音处理装置 | |
Furtado et al. | On the design of high-complexity cosine-modulated transmultiplexers based on the frequency-response masking approach | |
CN103714825A (zh) | 基于听觉感知模型的多通道语音增强方法 | |
CN107483019A (zh) | 一种白噪声产生方法和装置 | |
US20210152936A1 (en) | Information processing device, mixing device using the same, and latency reduction method | |
Bruschi et al. | Linear-Phase Octave Graphic Equalizer | |
CN102543091A (zh) | 一种模拟音效的生成系统及方法 | |
JP2017182099A (ja) | 音声合成装置、音声合成方法および音声合成プログラム | |
CN106601249B (zh) | 一种基于听觉感知特性的数字语音实时分解/合成方法 | |
JP5763487B2 (ja) | 音声合成装置、音声合成方法および音声合成プログラム | |
CN102208901B (zh) | 一种mp3解码过程中实现均衡器的方法 | |
Kamada et al. | Constant Q-value Filter Banks with Complex First Order System | |
Chhetri et al. | Implementation of Audio Effect Generator in FPGA | |
JP2014109667A (ja) | 音声合成装置、音声合成方法および音声合成プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: Wuxi District of Jiangsu city of Wuxi province Qingyuan Road 214000 No. 18 Taihu International Science Park sensor network university science and Technology Park 530 building A1001 Applicant after: WUXI ZHONGGAN MICROELECTRONIC CO., LTD. Address before: 214028 national integrated circuit design (21-1), Changjiang Road, New District, Jiangsu, Wuxi, China, China (610) Applicant before: Wuxi Vimicro Co., Ltd. |
|
COR | Change of bibliographic data | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |