CN101505443A - 一种虚拟重低音增强方法及*** - Google Patents

一种虚拟重低音增强方法及*** Download PDF

Info

Publication number
CN101505443A
CN101505443A CNA2009100799386A CN200910079938A CN101505443A CN 101505443 A CN101505443 A CN 101505443A CN A2009100799386 A CNA2009100799386 A CN A2009100799386A CN 200910079938 A CN200910079938 A CN 200910079938A CN 101505443 A CN101505443 A CN 101505443A
Authority
CN
China
Prior art keywords
signal
frequency
yield value
audio signal
unit
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
Application number
CNA2009100799386A
Other languages
English (en)
Other versions
CN101505443B (zh
Inventor
张晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuxi Zhonggan Microelectronics Co Ltd
Original Assignee
Vimicro Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Vimicro Corp filed Critical Vimicro Corp
Priority to CN2009100799386A priority Critical patent/CN101505443B/zh
Publication of CN101505443A publication Critical patent/CN101505443A/zh
Priority to US12/605,183 priority patent/US20100232624A1/en
Application granted granted Critical
Publication of CN101505443B publication Critical patent/CN101505443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00Circuits for transducers, loudspeakers or microphones
    • H04R3/04Circuits for transducers, loudspeakers or microphones for correcting frequency response
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R1/00Details of transducers, loudspeakers or microphones
    • H04R1/10Earpieces; Attachments therefor ; Earphones; Monophonic headphones
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R5/00Stereophonic arrangements
    • H04R5/033Headphones for stereophonic communication

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Stereophonic System (AREA)
  • Tone Control, Compression And Expansion, Limiting Amplitude (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明公开了一种虚拟重低音增强方法及***,用以提高重低音增强效果。本发明提供的一种虚拟重低音增强方法包括:通过提取音频信号中的低频信号,将该低频信号转换为频域信号,并确定该频域信号中的基频信号;对所述基频信号进行变调处理,得到多个谐波,并将所述多个谐波进行叠加,将叠加后得到的信号转换为时域信号;将该时域信号与所述音频信号进行合成,得到虚拟重低音增强后的音频信号。

Description

一种虚拟重低音增强方法及***
技术领域
本发明涉及信息处理技术领域,尤其涉及一种虚拟重低音增强方法及***。
背景技术
重低音增强(Bass Enhancement)算法是用来增强声音的低频部分的一种算法,例如增强音乐中的鼓声和贝斯声等低频部分的声音。
一般的音响和耳机都存在一定程度的低频损失,因此重低音的效果已经成为评价音频质量的重要标准之一,因此重低音增强算法很重要。
重低音增强的传统技术是EQ,而目前比较新的一种重低音增强技术是虚拟重低音(Virtual bass)技术,这种虚拟重低音技术主要是依靠产生一种用户心理感知的重低音印象,而并不真正增强低频部分的增益。依据的理论基础是,人感知重低音时,并不主要依靠该低音的基频,而更重要的是依靠该基频产生的各次谐波,即使低音的基频被抑制,只要其各次谐波以及这些谐波的关系依然存在,那么对于人来说,仍然能够感觉到强烈的重低音效果。
可以看出,虚拟重低音增强技术和传统的EQ技术有本质的区别,传统的EQ技术是通过增强音频信号的低频分量的能量来加强重低音效果,而虚拟重低音增强技术是通过增强低音基频的谐波来虚拟增强重低音效果的。
由于对于一般小型的音箱或者耳机来说,低频部分往往有很大的衰减,因此,即使用传统的EQ方法将低频部分增强,仍然不能得到满意的重低音效果,并且存在饱和溢出噪声的问题。而虚拟重低音增强技术增强的是基频的谐波部分,这些谐波的频率相对高些,因此不会受到音箱或者耳机硬件条件的制约,往往能得到更强的重低音效果。
现有的虚拟重低音增强技术一般是通过反馈调制产生低频信号的各次谐波,然而这会产生互调畸变噪声,从而影响了重低音的音质,也就是说,现有的虚拟重低音增强技术的重低音增强效果并不是很好。
发明内容
本发明实施例提供了一种虚拟重低音增强方法及***,用以提高重低音增强效果。
本发明实施例提供的一种虚拟重低音增强方法包括:
提取音频信号中的低频信号,将该低频信号转换为频域信号,并确定该频域信号中的基频信号;
对所述基频信号进行变调处理,得到多个谐波,并将所述多个谐波进行叠加,将叠加后得到的信号转换为时域信号;
将该时域信号与所述音频信号进行合成,得到虚拟重低音增强后的音频信号。
本发明实施例提供的一种虚拟重低音增强***包括:
第一低通滤波单元,用于提取音频信号中的低频信号;
频域转换单元,用于将该低频信号转换为频域信号;
基频确定单元,用于确定该频域信号中的基频信号;
变调单元,用于对所述基频信号进行变调处理,得到多个谐波;
第一合成单元,用于将所述多个谐波进行叠加;
时域转换单元,用于将所述第一合成单元处理后得到的信号转换为时域信号;
第二合成单元,用于将该时域信号与所述音频信号进行合成,得到虚拟重低音增强后的音频信号。
本发明实施例,通过提取音频信号中的低频信号,将该低频信号转换为频域信号,并确定该频域信号中的基频信号;对所述基频信号进行变调处理,得到多个谐波,并将所述多个谐波进行叠加,将叠加后得到的信号转换为时域信号;将该时域信号与所述音频信号进行合成,得到虚拟重低音增强后的音频信号,从而在增强了音频信号的虚拟重低音的同时,避免产生噪声,提高了重低音增强效果。
附图说明
图1为本发明实施例提供的一种虚拟重低音增强***的结构示意图;
图2为本发明实施例提供的斜坡函数示意图;
图3为本发明实施例提供的一种虚拟重低音增强方法的流程示意图。
具体实施方式
本发明实施例提供了一种虚拟重低音增强方法及***,用以可靠地产生音频信号的低频部分的各次谐波,并且不产生噪声,从而提高重低音增强效果。
本发明实施例提供的新的虚拟重低音增强的数字信号处理技术,通过滤波和变采样率技术将音频信号的低频部分单独抽取出来,然后将其变换到频域,再采用频域变调技术将低频部分的基频信号,按照各次谐波对应的频率进行变调,最后将处理后的信号恢复到时域,并与原音频信号合成,从而在增强了重低音效果的同时,不产生噪声。
下面结合附图对本发明实施例提供的技术进行说明。
参见图1,本发明实施例提供的一种虚拟重低音增强***包括:
第一低通滤波单元11,用于按照预设的截止频率对输入的音频信号进行低通滤波处理,得到该音频信号中的低频信号,该低频信号是时域的低频信号。
抽取单元12,用于按照预设的抽取倍数,对第一低通滤波单元11输出的低频信号进行抽取处理,并将处理后的信号发送给频域转换单元13。
频域转换单元13,用于将抽取单元12发送的低频信号转换为频域信号。
基频确定单元14,用于对频域转换单元13发送的频域信号进行分析,确定其中的基频信号。
变调单元15,用于对基频信号进行变调处理,得到多个谐波。
第一合成单元16,用于将变调单元15发送的多个谐波进行叠加,即将多个谐波进行复数相加,也就是说,将用于表示每个谐波的复数的实部和实部相加,将得到的和作为新的复数的实部;将用于表示每个谐波的复数的虚部和虚部相加,将得到的和作为新的复数的虚部;最后将新的复数所表示的信号作为叠加后的信号。
时域转换单元17,用于将第一合成单元16处理后得到的信号转换为时域信号。
插值单元18,用于按照预设的插值倍数,对所述时域转换单元17输出的信号进行插值处理,并将处理后的信号发送给第二低通滤波单元19。
第二低通滤波单元19,用于按照预先设置的截止频率,将插值单元18发送的信号进行低通滤波处理,并将处理后的信号发送给第二合成单元21。
延时处理单元20,用于按照预设时延值对输入的音频信号进行延时处理后发送给第二合成单元21。
第二合成单元21,将延时处理单元20输出的音频信号,与第二低通滤波单元19输出的信号进行合成。
自动增益控制单元22,用于将经过所述第二合成单元21合成的信号进行自动增益控制(AGC)。
较佳地,所述频域转换单元13包括:
分析窗单元,用于将抽取单元12输出的低频信号进行分析窗处理。
FFT单元,用于对经过分析窗单元处理后的信号进行快速傅里叶变换(FFT),得到低频信号的频域信号。
较佳地,所述时域转换单元17包括:
综合窗单元,用于将第一合成单元16输出的频域的低频信号进行综合窗处理。
IFFT单元,用于对经过综合窗单元处理后的信号进行快速傅里叶逆变换(IFFT),得到时域的低频信号。
下面对上述各个单元进行详细说明。
上述第一低通滤波单元11和第二低通滤波单元19的功能实际上是一样的,就是一个简单的低通滤波器(LPF),LPF主要作用是把信号的低频部分滤出来。
LPF的截止频率的设定需要考虑两个方面:一方面不能太小,太小容易使低频成分被衰减;另一方面也不能太大,因为后续还要通过抽取单元12进行进一步处理,如果截止频率太大,会容易造成信号频谱混叠。
通常音频信号中的1千赫兹(KHz)以下的频率部分就包含了几乎所有的低音成分,因此本发明实施例中的截止频率应不小于1KHz。
例如:抽取单元12采用M倍抽取比例,音频信号的采样率为44.1KHz,则抽取单元12抽取得到的信号的采样率降低为约44.1KHz/M,而信号的频率应在44.1KHz/2M以下才不会发生混叠,因此本发明实施例的第一低通滤波单元11和第二低通滤波单元19中预设的截止频率应不大于44.1KHz/2M。
抽取单元12和插值单元18:
抽取单元12和插值单元18主要采用采样率转换技术,抽取单元12的抽取即从输入的信号序列中,每隔M个点抽取一个点,M即是抽取倍数。相应地,插值单元18的插值是在输入的信号序列中,在每个点后***M-1个零,M即是插值倍数,其值与抽取倍数的值相同。
设置抽取单元12和插值单元18进行采样率转换的操作的目的,是通过降低采样率,使得频域转换单元13和时域转换单元17在较低的采样率下工作,因此可以大大降低运算复杂度。
进一步,本发明实施例考虑到虽然降采样可以减小FFT和IFFT的处理点数,但是增加的低通滤波操作会增加额外的运算量,并且降采样的倍数越高,低通滤波器的通带越窄,满足要求的滤波器的阶数越高。因此,需要折衷(tradeoff)。经过试验,本发明实施例可以选择8倍降采样处理,即M=8,例如,输入的音频信号的采样率为44.1KHz,则低通滤波器的截止频率fs必须满足:fs≤44100/2/8,即fs≤2756Hz。
本发明实施例采用的低通滤波器可以采用1.5KHz的截止频率的,64阶的FIR低通滤波器。阻带衰减大于50分贝(dB),低通滤波器由matlab设计。本发明实施例之所以没有选择在同样性能下阶数可以更低的IIR低通滤波器,原因包括:
FIR滤波器虽然有64阶,但是可以和抽取单元12和插值单元18一起构成快速算法***,真正的复杂度相当于64/8=8阶的FIR滤波器,因此算法的复杂度并不高;
而IIR滤波器由于有反馈操作,因此必须逐点运算,不能和变采样单元(即抽取单元12和插值单元18)组成快速算法***,另外,由于IIR滤波对数据精度要求较高,量化误差较大,因此也给设计带来一定困难,用较高精度的数据参与运算,往往增加运算量;
FIR滤波器具有线性相位,即所有频率具有相同的群延时,这点非常关键,不会带来相位失真。更为关键的是,在将增强后的低频信号与原信号叠加的时候,有一个相位对准的问题。如果采用IIR滤波器处理不当,就可能出现因反相而信号对消(phase cancel)的情况,因此会降低音频质量。
对于分析窗单元、综合窗单元、FFT单元和IFFT单元这些基本的数字信号处理单元,由于主要采用的是现有技术,因此本发明实施例不再赘述。
较佳地,基频确定单元14包括:
各频带频率确定单元,用于确定频域转换单元13输出的频域信号中各个频带(bin)对应的频率。
频带选取单元,用于按照预先设置的频率范围选取多个频带,并将这些多个频带中幅度最大的频带所对应的频率作为基频信号的频率,将该频带的幅度作为基频信号的幅度。
首先对各频带频率确定单元的功能进行详细说明。
信号通过FFT得到的频域信号,分为实部和虚部两部分,分别用Real和Imag表示。
假设FFT的点数用fftsize表示,则Real和Imag各为fftsize/2长度的序列。
可以通过实部和虚部得到频域信号用幅度和相位表示的形式,幅度和相位分别可以用magn和phase表示,则有:
Magn = Real * Real + Imag * Imag
Phase = arctg ( Imag Real )
得到了幅度和相位后,可以计算信号在每个bin里的准确频率,bin表示信号经过FFT变换到频域后的每个条带,即频带。例如,如果FFT点数为128,则共有128个bin。
第i个bin里的相位可以用Phase(i)表示,下面以计算第i个bin里的信号的准确频率为例进行说明。
假设上一帧第i个bin对应的相位为Phase_old(i),则可知上一帧和当前帧这两帧在该bin里的相位差Tmp为:
Tmp=Phase(i)-Phase_old(i)
因为第i个bin的标准相位差TmpS为:
TmpS = 2 πi stepsize fftsize
其中,stepsize表示一次处理信号的步长,一般来说,stepsize比fftsize要小,有一定的重叠,这样处理会更加准确。较佳地,本发明实施例中取stepsize为fftsize的四分之一长度,即M=fftsize/stepsize=4
计算Tmp与TmpS的差值TmpD:
TmpD=Tmp-TmpS
将这个差值TmpD规划到正负π之间,得到TmpD′,由此可以计算得到频率偏差FreqD为:
FreqD = TmpD ′ 2 π * M * FreqPerBin
最后,将FreqD与标准频率相加,即可得到第i个bin对应的准确频率值FreqS(i)为:
FreqS(i)=i*FreqPerBin+FreqD
其中,FreqPerBin表示每个频带的带宽。
下面详细介绍一下频带选取单元的功能。
假设原始输入的音频信号的采样率为44.1KHz,经过8倍抽取后,采样率降为约5.5KHz,如果采用256点的FFT,则每个bin里的频率范围约为20Hz,一般来说,基频的频率很低,在80Hz以下,因此只需对频率最低的几个bin进行搜索即可,较佳地,选取频率最低的4个bin,即在频率最低的四个bin里搜索基频信号,查看哪个bin里的幅度Magn最大,Magn最大的bin所对应的频率FreqS即是要寻找的基频信号的频率F,该bin对应的幅度即是基频信号的幅度MF。即:
F_i=arg[Max(Magn(i))],i=0~3
F=FreqS[F_i]
MF=Magn[F_i]
较佳地,所述变调单元15包括:
谐波频率确定单元,用于将基频信号的频率与预先设置的多个整数分别相乘,得到多个谐波的频率。
谐波幅度确定单元,用于将所述基频信号的幅度与预先设置的多个幅度比例因子分别相乘,得到多个谐波的幅度。
变调单元15的主要任务就是通过基频信号的频率F和幅度MF得到基频信号的各个谐波成分。
基频信号的各个谐波的频率都是基频信号的频率F的整数倍。较佳地,本发明实施例中只考虑最低的5次谐波,则各个谐波的频率和幅度分别为:
Fh(k)=kF,k=1,2,3,4,5
MFh(k)=a(k)MF,k=1,2,3,4,5
其中,a(k)表示第k次谐波的幅度比例因子,每个谐波的幅度比例因子是不同的,并且是预先设置好的,一般来说频率越高的谐波,能量越低,因此a(k)是一个大于零的小数,并且随着k的增大而减小。
然后,根据各个谐波的频率计算其准确相位,以第k次谐波为例,假设第k次谐波的频率Fh(k)在第i个bin的范围内,则Fh(k)与该bin内的标准频率的差的归一化值FreqD为:
FreqD=(Fh(k)-i*FreqPerBin)/FreqPerBin
计算相对相位差TmpD为:
TmpD = 2 π FreqD M
TmpD与标准相位差相加得到准确的相位差Tmp:
Tmp = TmpD + 2 πi stepsize fftsize
将准确的相位差Tmp与计算得到Tmp之前累计的相位差Tmp_sum相加得到最终相位Phase(k)为:
Phase(k)=Tmp+Tmp_sum
并且,更新累计相位差Tmp_sum=Phase(k),其中,Tmp_sum的初始值为0。
最后,通过幅度MF(k)和相位Phase(k)就可以计算出该谐波的实部Real(i)和虚部Imag(i):
Real(i)=MF(k)*cos(Phase(k))
Imag(i)=MF(k)*sin(Phase(k))
第一合成单元16将所有的谐波进行叠加后就得到了基频信号的所有谐波成分的频域信号,时域转换单元17对该频域信号进行IFFT变换,就得到了其时域信号。
延时处理单元20:
延时处理单元20,将原始信号延时D个样本,D值即所谓的时延值。延时的目的是为了使得原始的音频信号与变调后的基频信号的相位对齐,以避免相位不对齐时造成信号对消。D值的设置方式如下:
D值的确定需要考虑对音频信号的基频部分进行一系列处理过程中的所有可能的延时,包括:第一低通滤波单元11和第二低通滤波单元19中的滤波器的长度,分析窗和综合窗的长度,以及FFT和IFFT转换所需要占用的时间等等。
假设第一低通滤波单元11和第二低通滤波单元19中的滤波器的长度为L,分析窗和综合窗的长度为W,则D值可以为:
D=L/2*2+W/2*M
其中,L/2为一个LPF的延时,共有两个LPF,所以低通滤波处理造成的延时为L;W/2为分析窗处理和综合窗处理造成的延时,因为这部分延时是在抽取后产生的,因此对于抽取前来说相当于还要增大M倍,因此为W/2*M
将延时了D个样本的原始音频信号与变调后的基频信号相加,可能发生饱和溢出,因此相加后的信号进一步需要进入自动增益控制单元22进行处理。
自动增益控制单元22:
一般的AGC模块用来自动改变信号的增益,将小信号放大,将大信号缩小,音量保持适中。而本发明实施例中的自动增益控制单元22并不是这样,因为对于音乐来说,旋律的起承转合,抑扬顿挫是自身的特点,并不能破坏,本发明实施例中使用AGC的目的是在保证声音不发生饱和失真的前提下,提高重低音的音量。也就是说,自动增益控制单元22为了使得一定时间范围内的幅度最大的信号处于饱和边界上,而该范围内的信号大小关系仍然保留,需要采用快降慢升的方法。
较佳地,自动增益控制单元22包括:
第一增益值单元,用于确定第二合成单元21得到的当前帧音频信号中绝对值最大的信号幅度值,将该信号幅度值与预先设置的目标阈值进行比较,得到第一增益值。
第二增益值单元,用于将第一增益值与上一帧音频信号采用的增益值进行比较,当第一增益值小于上一帧音频信号采用的增益值时,令第二增益值等于第一增益值;当第一增益值大于上一帧音频信号采用的增益值时,令第二增益值等于上一帧音频信号采用的增益值与预先设置的步长的和;其中,第二增益值属于预先设置的阈值范围。
帧内平滑单元,用于利用上一帧音频信号采用的增益值,通过斜坡函数对第二增益值做帧内平滑,得到当前帧音频信号采用的增益值。
输出单元,用于将当前帧音频信号采用的增益值与当前帧合成后的音频信号相乘,得到并输出自动增益控制后的音频信号。
例如,在第二合成单元21输出的当前帧音频信号中找到的绝对值最大的信号幅度值为Vmax,然后将Vmax和目标阈值Ti比较,目标阈值是希望信号幅度能够达到的理想值。将Ti和Vmax相比得到理想的增益值gain_t(即所述的第一增益值)为:
gain_t=Ti/Vmax
由于增益变化过快会带来信号突变噪声,因此,本发明实施例采用快降慢升的增益调整方式,包括:
假设当前帧的上一帧计算得到的最终增益值为gain_old,则:
如果gain_t<gain_old,则gain=gain_t,此操作表示快降,gain即所述的第二增益值,最小可降低到一个低阈值LowLimit。
如果gain_t>gain_old,则gain=gain_old+step,此操作表示慢升,其中step为预先设置的增益增大时的过渡步长,gain最大可增大到一个高阈值HighLimit。
也就是说,gain无论怎么调整,都要满足:gainLowLimit≤gain≤HighLimit。
进一步,用新计算出的增益gain和上一帧的gain_old做帧内平滑,可用如图2所示的斜坡函数来进行加权,斜坡函数定义为b(i)=1-i/N,则:
gainW(i)=b(i)gain_old+(1-b(i))gain,i=0~N-1
其中,gainW(i)为当前帧采用的做了帧内平滑后的样点i的增益,N表示每个帧的长度。
可以看出,由于斜坡函数在开始时对于上一帧的gain_old给与了较大权值,对于当前帧的gain给与了较小权值;而在末尾时正好相反。因此可以有效地平滑增益突变的影响。
最后,用做了帧内平滑的增益gainW(i)去处理第二合成单元21的输出信号,即自动增益控制单元22的输入信号input(i),得到输出信号output(i)为:
output(i)=input(i)*gainW(i),i=0~N-1
参见图3,本发明实施例提供的一种虚拟重低音增强方法包括步骤:
S101、按照预先设置的截止频率,提取音频信号中的低频信号,并按照预设的抽取倍数,对所述低频信号进行抽取处理。
S102、将抽取处理得到的低频信号经过分析窗的处理后,进行FFT,得到低频信号的频域信号,并确定该频域信号中各个频带对应的频率。
S103、按照预先设置的频率范围选取多个频带,并将多个频带中幅度最大的频带所对应的频率作为基频信号的频率,将该频带的幅度作为基频信号的幅度。
S104、将基频信号的频率与预先设置的多个整数分别相乘,得到多个谐波的频率;将基频信号的幅度与预先设置的多个幅度比例因子分别相乘,得到多个谐波的幅度;然后将多个谐波叠加,并将叠加后得到的信号通过IFFT转换为时域信号。
S105、对该时域信号进行综合窗的处理后,按照预设的插值倍数,对经过综合窗的处理后的时域信号进行插值处理;按照预先设置的截止频率,将经过插值处理后得到的信号进行低通滤波处理,得到变调处理后的基频信号。
S106、按照预设时延值对原始输入的音频信号进行延时处理,将经过延时处理后的音频信号与得到变调处理后的基频信号进行合成。
S107、将当前帧音频信号与当前帧变调处理后的基频信号合成后得到的音频信号进行自动增益控制。
其中的抽取步骤和插值步骤都是较佳的步骤,不是必需的步骤。
步骤S107中的自动增益控制包括:
步骤一:确定将当前帧音频信号与当前帧变调处理后的基频信号合成后得到的音频信号中绝对值最大的信号幅度值,将该信号幅度值与预先设置的目标阈值进行比较,得到第一增益值。
步骤二:将第一增益值与上一帧音频信号采用的增益值进行比较,当第一增益值小于上一帧音频信号采用的增益值时,令第二增益值等于第一增益值;当第一增益值大于上一帧音频信号采用的增益值时,令第二增益值等于上一帧音频信号采用的增益值与预先设置的步长的和;其中,第二增益值属于预先设置的阈值范围。
步骤三:利用第二增益值对当前帧合成后的音频信号进行增益控制。
较佳地,步骤三包括:
利用上一帧音频信号采用的增益值,通过斜坡函数对第二增益值做帧内平滑,得到当前帧音频信号采用的增益值,并将该增益值与当前帧合成后的音频信号相乘,得到当前帧自动增益控制后的音频信号。
综上所述,本发明实施例提供了一种新的虚拟重低音增强的数字信号处理方法及***,通过低通滤波和变采样率技术将信号低频部分抽取出来,然后用FFT将其变换到频域,再用频域变调技术,将基频信号按照各次谐波对应得频率进行变调,最后将处理后的基频信号通过IFFT恢复到时域,并与原信号合成,从而在不产生噪声的前提下,增强了音频信号的虚拟重低音。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (15)

1、一种虚拟重低音增强方法,其特征在于,该方法包括:
提取音频信号中的低频信号,将该低频信号转换为频域信号,并确定该频域信号中的基频信号;
对所述基频信号进行变调处理,得到多个谐波,并将所述多个谐波进行叠加,将叠加后得到的信号转换为时域信号;
将该时域信号与所述音频信号进行合成,得到虚拟重低音增强后的音频信号。
2、根据权利要求1所述的方法,其特征在于,将所述低频信号转换为频域信号,并确定该频域信号中的基频信号的步骤包括:
将所述低频信号经过分析窗的处理后,进行快速傅里叶变换,得到所述低频信号的频域信号,并确定该频域信号中各个频带对应的频率;
按照预先设置的频率范围选取多个频带,并将所述多个频带中幅度最大的频带所对应的频率作为基频信号的频率,将该频带的幅度作为基频信号的幅度。
3、根据权利要求1所述的方法,其特征在于,对所述基频信号进行变调处理,得到多个谐波的步骤包括:
将所述基频信号的频率与预先设置的多个整数相乘,得到多个谐波的频率;
将所述基频信号的幅度与预先设置的多个幅度比例因子分别相乘,得到多个谐波的幅度。
4、根据权利要求3所述的方法,其特征在于,在提取了所述低频信号之后,将该低频信号转换为频域信号之前,还包括步骤:
按照预设的抽取倍数,对所述低频信号进行抽取处理;
在将所述叠加后得到的信号转换为时域信号之后,将该时域信号与所述音频信号进行合成之前,还包括步骤:
对所述时域信号进行综合窗的处理后,按照预设的插值倍数,对经过综合窗的处理后的时域信号进行插值处理;
按照预先设置的截止频率,将经过所述插值处理后得到的信号进行低通滤波处理。
5、根据权利要求1至4任一权利要求所述的方法,其特征在于,将所述时域信号与所述音频信号进行合成的步骤包括:
按照预设时延值对所述音频信号进行延时处理;
将经过所述延时处理后的音频信号,与所述时域信号进行合成。
6、根据权利要求1所述的方法,其特征在于,将所述时域信号与所述音频信号进行合成之后,还包括步骤:
将所述时域信号与所述音频信号合成后得到的音频信号进行自动增益控制。
7、根据权利要求6所述的方法,其特征在于,所述自动增益控制包括:
确定所述时域信号与所述音频信号合成后得到的当前帧音频信号中绝对值最大的信号幅度值,将该信号幅度值与预先设置的目标阈值进行比较,得到第一增益值;
将所述第一增益值与上一帧音频信号采用的增益值进行比较,当所述第一增益值小于上一帧音频信号采用的增益值时,令第二增益值等于所述第一增益值;当所述第一增益值大于上一帧音频信号采用的增益值时,令第二增益值等于上一帧音频信号采用的增益值与预先设置的步长的和;其中,所述第二增益值属于预先设置的阈值范围;
利用所述第二增益值对当前帧合成后的音频信号进行增益控制。
8、根据权利要求7所述的方法,其特征在于,利用所述第二增益值对当前帧合成后的音频信号进行增益控制的步骤包括:
利用上一帧音频信号采用的增益值,通过斜坡函数对所述第二增益值做帧内平滑,得到当前帧音频信号采用的增益值,并将该增益值与当前帧合成后的音频信号相乘,得到自动增益控制后的音频信号。
9、一种虚拟重低音增强***,其特征在于,该***包括:
第一低通滤波单元,用于提取音频信号中的低频信号;
频域转换单元,用于将该低频信号转换为频域信号;
基频确定单元,用于确定该频域信号中的基频信号;
变调单元,用于对所述基频信号进行变调处理,得到多个谐波;
第一合成单元,用于将所述多个谐波进行叠加;
时域转换单元,用于将所述第一合成单元处理后得到的信号转换为时域信号;
第二合成单元,用于将该时域信号与所述音频信号进行合成,得到虚拟重低音增强后的音频信号。
10、根据权利要求9所述的***,其特征在于,所述基频确定单元包括:
各频带频率确定单元,用于确定该频域信号中各个频带对应的频率;
频带选取单元,用于按照预先设置的频率范围选取多个频带,并将所述多个频带中幅度最大的频带所对应的频率作为基频信号的频率,将该频带的幅度作为基频信号的幅度。
11、根据权利要求9或10所述的***,其特征在于,所述变调单元包括:
谐波频率确定单元,用于将所述基频信号的频率与预先设置的多个整数分别相乘,得到多个谐波的频率;
谐波幅度确定单元,用于将所述基频信号的幅度与预先设置的多个幅度比例因子分别相乘,得到多个谐波的幅度。
12、根据权利要求9所述的***,其特征在于,该***还包括:
抽取单元,用于按照预设的抽取倍数,对所述第一低通滤波单元输出的低频信号进行抽取处理,并将处理后的信号发送给所述频域转换单元;
插值单元,用于按照预设的插值倍数,对所述时域转换单元输出的信号进行插值处理,并将处理后的信号发送给第二低通滤波单元;
第二低通滤波单元,用于按照预先设置的截止频率,将经过所述插值处理后得到的信号进行低通滤波处理,并将处理后的信号发送给所述第二合成单元。
13、根据权利要求12所述的***,其特征在于,该***还包括:
延时处理单元,用于按照预设时延值对所述音频信号进行延时处理;
所述第二合成单元,将经过所述延时处理后的音频信号,与所述第二低通滤波单元输出的信号进行合成。
14、根据权利要求13所述的***,其特征在于,该***还包括:
自动增益控制单元,用于将经过所述第二合成单元合成的信号进行自动增益控制。
15、根据权利要求14所述的***,其特征在于,所述自动增益控制单元包括:
第一增益值单元,用于确定所述时域信号与所述音频信号合成后得到的当前帧音频信号中绝对值最大的信号幅度值,将该信号幅度值与预先设置的目标阈值进行比较,得到第一增益值;
第二增益值单元,用于将所述第一增益值与上一帧音频信号采用的增益值进行比较,当所述第一增益值小于上一帧音频信号采用的增益值时,令第二增益值等于所述第一增益值;当所述第一增益值大于上一帧音频信号采用的增益值时,令第二增益值等于上一帧音频信号采用的增益值与预先设置的步长的和;其中,所述第二增益值属于预先设置的阈值范围;
帧内平滑单元,用于利用上一帧音频信号采用的增益值,通过斜坡函数对所述第二增益值做帧内平滑,得到当前帧音频信号采用的增益值;
输出单元,用于将当前帧音频信号采用的增益值与当前帧合成后的音频信号相乘,得到并输出自动增益控制后的音频信号。
CN2009100799386A 2009-03-13 2009-03-13 一种虚拟重低音增强方法及*** Active CN101505443B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN2009100799386A CN101505443B (zh) 2009-03-13 2009-03-13 一种虚拟重低音增强方法及***
US12/605,183 US20100232624A1 (en) 2009-03-13 2009-10-23 Method and System for Virtual Bass Enhancement

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100799386A CN101505443B (zh) 2009-03-13 2009-03-13 一种虚拟重低音增强方法及***

Publications (2)

Publication Number Publication Date
CN101505443A true CN101505443A (zh) 2009-08-12
CN101505443B CN101505443B (zh) 2013-12-11

Family

ID=40977465

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100799386A Active CN101505443B (zh) 2009-03-13 2009-03-13 一种虚拟重低音增强方法及***

Country Status (2)

Country Link
US (1) US20100232624A1 (zh)
CN (1) CN101505443B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102354500A (zh) * 2011-08-03 2012-02-15 华南理工大学 一种基于谐波控制的虚拟低音增强处理方法
CN104704855A (zh) * 2012-10-15 2015-06-10 杜比国际公司 用于减小基于换位器的虚拟低音***中的延迟的***及方法
CN105632509A (zh) * 2014-11-07 2016-06-01 Tcl集团股份有限公司 一种音频处理方法和音频处理装置
CN106162443A (zh) * 2016-08-29 2016-11-23 广州声姆音响设备有限公司 一种音响***
CN106210988A (zh) * 2016-08-29 2016-12-07 广州声姆音响设备有限公司 一种音响***的低音补偿方法和电路
CN106653049A (zh) * 2015-10-30 2017-05-10 国光电器股份有限公司 时域中的虚拟低音的相加
CN106658284A (zh) * 2015-10-30 2017-05-10 国光电器股份有限公司 频域中的虚拟低音的相加
CN107959906A (zh) * 2017-11-20 2018-04-24 英业达科技有限公司 音效增强方法及音效增强***
CN108495235A (zh) * 2018-05-02 2018-09-04 北京小鱼在家科技有限公司 一种重低音的分离方法、装置、计算机设备及存储介质
CN108632708A (zh) * 2017-03-23 2018-10-09 展讯通信(上海)有限公司 扬声器输出控制方法及***
CN104952455B (zh) * 2015-06-19 2019-03-15 珠海市杰理科技股份有限公司 实现混响的方法和装置
CN114268886A (zh) * 2021-11-17 2022-04-01 厦门立林科技有限公司 虚拟低音的优化方法、***、智能终端及存储介质
CN115442709A (zh) * 2022-07-29 2022-12-06 荣耀终端有限公司 音频处理方法、虚拟低音增强***、设备和存储介质

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SG144752A1 (en) * 2007-01-12 2008-08-28 Sony Corp Audio enhancement method and system
KR101310231B1 (ko) * 2007-01-18 2013-09-25 삼성전자주식회사 저음 증강 장치 및 방법
JP5074115B2 (ja) * 2007-07-12 2012-11-14 ラピスセミコンダクタ株式会社 音響信号処理装置及び音響信号処理方法
US8971551B2 (en) 2009-09-18 2015-03-03 Dolby International Ab Virtual bass synthesis using harmonic transposition
WO2013181299A1 (en) 2012-05-29 2013-12-05 Creative Technology Ltd Adaptive bass processing system
US9247342B2 (en) 2013-05-14 2016-01-26 James J. Croft, III Loudspeaker enclosure system with signal processor for enhanced perception of low frequency output
US10405094B2 (en) * 2015-10-30 2019-09-03 Guoguang Electric Company Limited Addition of virtual bass
US10893362B2 (en) 2015-10-30 2021-01-12 Guoguang Electric Company Limited Addition of virtual bass
CN111198789B (zh) * 2019-12-20 2023-06-09 北京时代民芯科技有限公司 一种fft硬件实现模块的验证方法
US11617046B2 (en) * 2021-07-02 2023-03-28 Tenor Inc. Audio signal reproduction
CN116486833B (zh) * 2023-06-21 2023-09-22 北京探境科技有限公司 音频增益调整方法、装置、存储介质及电子设备

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1550121A (zh) * 2001-09-21 2004-11-24 在电声变换器中控制音频信号的低音放音的方法和装置
US20070253576A1 (en) * 2006-04-27 2007-11-01 National Chiao Tung University Method for virtual bass synthesis
CN101373961A (zh) * 2007-08-23 2009-02-25 索尼株式会社 信号处理设备、信号处理方法及其程序

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7058188B1 (en) * 1999-10-19 2006-06-06 Texas Instruments Incorporated Configurable digital loudness compensation system and method
EP1473965A2 (en) * 2003-04-17 2004-11-03 Matsushita Electric Industrial Co., Ltd. Acoustic signal-processing apparatus and method
CN1830148B (zh) * 2003-07-29 2010-11-24 松下电器产业株式会社 音频信号频带扩展装置
SG144752A1 (en) * 2007-01-12 2008-08-28 Sony Corp Audio enhancement method and system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1550121A (zh) * 2001-09-21 2004-11-24 在电声变换器中控制音频信号的低音放音的方法和装置
US20070253576A1 (en) * 2006-04-27 2007-11-01 National Chiao Tung University Method for virtual bass synthesis
CN101373961A (zh) * 2007-08-23 2009-02-25 索尼株式会社 信号处理设备、信号处理方法及其程序

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MINGSIAN R. BAI ET. AL.: "Synthesis and Implementation of Virtual Bass System with a Phase-Vocoder Approach", 《JOURNAL AUDIO ENGINEERING SOCIETY》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102354500A (zh) * 2011-08-03 2012-02-15 华南理工大学 一种基于谐波控制的虚拟低音增强处理方法
CN104704855A (zh) * 2012-10-15 2015-06-10 杜比国际公司 用于减小基于换位器的虚拟低音***中的延迟的***及方法
CN104704855B (zh) * 2012-10-15 2016-08-24 杜比国际公司 用于减小基于换位器的虚拟低音***中的延迟的***及方法
CN105632509A (zh) * 2014-11-07 2016-06-01 Tcl集团股份有限公司 一种音频处理方法和音频处理装置
CN105632509B (zh) * 2014-11-07 2019-07-19 Tcl集团股份有限公司 一种音频处理方法和音频处理装置
CN104952455B (zh) * 2015-06-19 2019-03-15 珠海市杰理科技股份有限公司 实现混响的方法和装置
CN106658284B (zh) * 2015-10-30 2020-07-17 国光电器股份有限公司 频域中的虚拟低音的相加
CN106653049A (zh) * 2015-10-30 2017-05-10 国光电器股份有限公司 时域中的虚拟低音的相加
CN106658284A (zh) * 2015-10-30 2017-05-10 国光电器股份有限公司 频域中的虚拟低音的相加
CN106210988A (zh) * 2016-08-29 2016-12-07 广州声姆音响设备有限公司 一种音响***的低音补偿方法和电路
CN106162443B (zh) * 2016-08-29 2019-04-26 广州声姆音响设备有限公司 一种音响***
CN106162443A (zh) * 2016-08-29 2016-11-23 广州声姆音响设备有限公司 一种音响***
CN108632708A (zh) * 2017-03-23 2018-10-09 展讯通信(上海)有限公司 扬声器输出控制方法及***
CN108632708B (zh) * 2017-03-23 2020-04-21 展讯通信(上海)有限公司 扬声器输出控制方法及***
CN107959906A (zh) * 2017-11-20 2018-04-24 英业达科技有限公司 音效增强方法及音效增强***
CN108495235A (zh) * 2018-05-02 2018-09-04 北京小鱼在家科技有限公司 一种重低音的分离方法、装置、计算机设备及存储介质
CN108495235B (zh) * 2018-05-02 2020-10-09 北京小鱼在家科技有限公司 一种重低音的分离方法、装置、计算机设备及存储介质
CN114268886A (zh) * 2021-11-17 2022-04-01 厦门立林科技有限公司 虚拟低音的优化方法、***、智能终端及存储介质
CN114268886B (zh) * 2021-11-17 2023-06-30 厦门立林科技有限公司 虚拟低音的优化方法、***、智能终端及存储介质
CN115442709A (zh) * 2022-07-29 2022-12-06 荣耀终端有限公司 音频处理方法、虚拟低音增强***、设备和存储介质

Also Published As

Publication number Publication date
CN101505443B (zh) 2013-12-11
US20100232624A1 (en) 2010-09-16

Similar Documents

Publication Publication Date Title
CN101505443B (zh) 一种虚拟重低音增强方法及***
CN101577848B (zh) 一种重低音增强方法及***
US10750278B2 (en) Adaptive bass processing system
JP5341128B2 (ja) 補聴器における安定性の改善
US9210506B1 (en) FFT bin based signal limiting
US20030216907A1 (en) Enhancing the aural perception of speech
EP2579252B1 (en) Stability and speech audibility improvements in hearing devices
US9672834B2 (en) Dynamic range compression with low distortion for use in hearing aids and audio systems
US6370255B1 (en) Loudness-controlled processing of acoustic signals
US8306821B2 (en) Sub-band periodic signal enhancement system
KR20040035749A (ko) 사운드 신호의 대역폭 확장 방법
EP3100353B1 (en) An audio compression system for compressing an audio signal
CN106340303A (zh) 一种基于时间频率域的语音降噪方法
CN111739544B (zh) 语音处理方法、装置、电子设备及存储介质
CN111970627A (zh) 音频信号的增强方法、装置、存储介质和处理器
US20010027390A1 (en) Speech decoder and a method for decoding speech
US10199048B2 (en) Bass enhancement and separation of an audio signal into a harmonic and transient signal component
US8144762B2 (en) Band extending apparatus and method
CN1134768C (zh) 由时域频谱相减进行信号噪声减小的电话机、处理器及方法
CN110858487A (zh) 一种音频信号缩放处理方法及装置
CN101221767B (zh) 人声语音加强装置与应用于其上的方法
EP3335218B1 (en) An audio signal processing apparatus and method for processing an input audio signal
US20230051509A1 (en) Method and system for implementing a modal processor
US20230217166A1 (en) Bass enhancement for loudspeakers
KR20010076265A (ko) 디지털 그래퍼메트릭 이퀄라이저

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: WUXI VIMICRO CORPORATION

Free format text: FORMER OWNER: BEIJING ZHONGXING MICROELECTRONICS CO., LTD.

Effective date: 20110328

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 100083 15/F, SHINING BUILDING, NO. 35, XUEYUAN ROAD, HAIDIAN DISTRICT, BEIJING TO: 214028 (CHUANGYUAN BUILDING), NATIONAL INTEGRATED CIRCUIT DESIGN PARK, NO. 21-1, YANGTES RIVER ROAD, WUXI NEW DISTRICT, JIANGSU PROVINCE

TA01 Transfer of patent application right

Effective date of registration: 20110328

Address after: 214028 national integrated circuit design Park, Changjiang Road, New District, Jiangsu,, Wuxi

Applicant after: Wuxi Vimicro Co., Ltd.

Address before: 100083, Haidian District, Xueyuan Road, Beijing No. 35, Nanjing Ning building, 15 Floor

Applicant before: Beijing Vimicro Corporation

C14 Grant of patent or utility model
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 214135 Taihu International Science Park Sensor Network University Science Park 530 Building A1001, 18 Qingyuan Road, Wuxi, Jiangsu Province

Patentee after: WUXI ZHONGGAN MICROELECTRONIC CO., LTD.

Address before: 214028 National Integrated Circuit Design Park 21-1 Changjiang Road, New District, Wuxi City, Jiangsu Province (Chuangyuan Building)

Patentee before: Wuxi Vimicro Co., Ltd.