语音增强

摘要

一种用于增强语音的方法，包括：提取音频信号的中心通道；使中心通道的谱变平；以及将变平的语音通道与音频信号混合，从而增强音频信号中的任何语音。还公开了一种用于从具有多个通道的音频信号提取声音的中心通道的方法，用于使音频信号的谱变平的方法，以及用于检测音频信号中的语音的方法。还公开了一种语音增强器，包括：中心通道提取器；平谱器；语音置信度生成器；以及混合器，用于与具有所检测的语音的置信度成比例地混合变平的语音通道与原始音频信号，从而增强音频信号中的任何语音。

权利要求

1.一种用于从具有多个通道的音频信号提取声音的中心通道的方法，所述方法包括：将(1)所述音频信号的第一通道减去候选中心通道的比例α与(2)所述音频信号的第二通道的共轭减去所述候选中心通道的所述比例α相乘；

近似最小化α；以及

通过将所述候选中心通道乘以近似最小化的α，产生所提取的中心通道。

2.一种用于使音频信号的谱变平的方法，所述方法包括：将假定的语音通道分成感知频带；

确定所述感知频带中的哪个具有最大的能量；以及增加具有较少能量的感知频带的增益，从而使所述音频信号中的任何语音的谱变平。

3.如权利要求2所述的方法，其中所述增加包括：将具有较少能量的感知频带的增益增加，直到最大。

4.一种用于检测音频信号中的语音的方法，所述方法包括：测量所述音频信号的候选中心通道中的谱波动；

测量所述音频信号减去所述候选中心通道的谱波动；以及比较所述谱波动，从而检测所述音频信号中的语音。

5.一种用于增强语音的方法，所述方法包括：提取音频信号的中心通道；

使所述中心通道的谱变平；以及将变平的语音通道与所述音频信号混合，从而增强所述音频信号中的任何语音。

6.如权利要求5所述的方法，还包括：生成检测所述中心通道中的语音的置信度；并且其中所述混合包括与具有所检测的语音的置信度成比例地混合所述变平的语音通道与所述音频信号。

7.如权利要求6所述的方法，其中：所述置信度从最低可能的概率到最高可能的概率变化，并且所述生成包括

进一步将生成的置信度限制到高于所述最低可能的概率并且低于所述最高可能的概率的值。

8.如权利要求5所述的方法，其中所述提取包括：使用权利要求1所述的方法，提取音频信号的中心通道。

9.如权利要求5所述的方法，其中所述变平包括：使用权利要求2所述的方法，使所述中心通道的谱变平。

10.如权利要求5所述的方法，其中所述生成包括：使用权利要求3所述的方法，生成检测所述中心通道中的语音的置信度。

11.如权利要求5所述的方法，其中所述提取包括：使用权利要求1所述的方法，提取音频信号的中心通道；其中所述变平包括：使用权利要求2所述的方法，使所述中心通道的谱变平；以及其中所述生成包括：

使用权利要求3所述的方法，生成检测所述中心通道中的语音的置信度。

12.一种计算机可读存储介质，其中设置了用于执行权利要求1-11中任一权利要求所述的方法的计算机程序。

13.一种计算机系统，包括：CPU；

权利要求12所述的存储介质；以及耦合所述CPU和所述存储介质的总线。

14.一种语音增强器，包括：中心通道提取器，用于提取音频信号的中心通道；

平谱器，用于使所述中心通道的谱变平；

语音置信度生成器，用于生成检测所述中心通道中的语音的置信度；以及混合器，用于与具有所检测的语音的置信度成比例地混合变平的语音通道与原始音频信号，从而增强所述音频信号中的任何语音。

说明书

语音增强

发明内容

[0001] 这里描述了如下方法和装置：该方法和装置用于从具有多个通道的音频信号提取声音的中心通道，用于使音频信号的谱变平，用于检测音频信号中的语音，以及用于增强语音。用于从具有多个通道的音频信号提取声音的中心通道的方法可包括：将(1)该音频信号的第一通道减去候选中心通道的比例α与(2)该音频信号的第二通道的共轭减去该候选中心通道的比例α相乘；近似最小化α；以及通过将该候选中心通道乘以近似最小化的α，产生所提取的中心通道。

[0002] 用于使音频信号的谱变平的方法可包括：将假定的语音通道分成感知频带；确定该感知频带中的哪个具有最大的能量；以及增加具有较少能量的感知频带的增益，从而使音频信号中的任何语音的谱变平。增加可包括将具有较少能量的感知频带的增益增加直至最大。

[0003] 用于检测音频信号中的语音的方法可包括：测量音频信号的候选中心通道中的谱波动；测量减去候选中心通道的音频信号的谱波动；以及比较谱波动，从而检测音频信号中的语音。

[0004] 用于增强语音的方法可包括：提取音频信号的中心通道；使中心通道的谱变平；以及将变平的语音通道与该音频信号混合，从而增强该音频信号中的任何语音。该方法还可包括生成检测中心通道中的语音的置信度，并且混合可包括与具有所检测的语音的置信度成比例地混合变平的语音通道与音频信号。置信度可从最低可能的概率到最高可能的概率变化，并且生成可包括进一步将生成的置信度限制到高于最低可能的概率并且低于最高可能的概率的值。提取可包括使用上述方法提取音频信号的中心通道。变平可包括使用上述方法使中心通道的谱变平。生成可包括使用上述方法生成检测中心通道中的语音的置信度。

[0005] 提取可包括使用上述方法提取音频信号的中心通道；变平可包括使用上述方法使中心通道的谱变平；并且生成可包括使用上述方法生成检测中心通道中的语音的置信度。

[0006] 这里教授了一种计算机可读存储介质以及计算机系统，其中该计算机可读存储介质中设置有用于执行任何上述方法的计算机程序，该计算机系统包括CPU、存储介质以及耦合CPU和存储介质的总线。

附图说明

[0007] 图1是根据本发明的一个实施例的语音增强器的功能框图。

[0008] 图2示出了具有1ERB(Equivalent Rectangular Bandwidth，等效矩形带宽)的间距从而导致共40个频带的适当滤波器集合。

[0009] 图3描述了根据本发明的一个实施例的混合过程。

[0010] 图4示出了根据本发明的一个实施例的计算机系统。

具体实施方式

[0011] 图1是根据本发明的一个实施例的语音增强器1的功能框图。语音增强器1包括输入信号17、离散傅里叶变换器10a、10b、中心通道提取器11、平谱器12、声音活动检测器13、可变增益放大器15a、15c，反离散傅里叶变换器18a、18b以及输出信号18。输入信号17分别包括左右通道17a、17b，并且输出信号18类似地分别包括左右通道18a、18b。

[0012] 相应的离散傅里叶变换器18接收输入信号17的左右通道17a、17b作为输入，并且产生变换19a、19b作为输出。中心频道提取器11接收变换19并且产生假想中心通道C20作为输出。平谱器12接收假想中心通道C 20作为输入，并且产生整形的中心通道24作为输出，同时声音活动检测器13接收同一输入C 20，并且一方面产生用于可变增益放大器
14a和14c的控制信号22以及另一方面产生用于可变增益放大器14b的控制信号21作为输出。

[0013] 放大器14a分别接收左通道变换19a和声音活动检测器13的输出控制信号22作为输入和控制信号。同样，放大器14c分别接收右通道变换19b和声音活动检测器输出控制信号22作为输入和控制信号。放大器14b接收平谱器12的经谱整形的中心通道24和输出声音活动检测器控制信号21作为输入和控制信号。

[0014] 混合器15a接收从放大器14输出的、经增益调整的左变换23a以及经增益调整的、经谱整形的中心通道25，并且产生信号26a作为输出。类似地，混合器15b接收来自放大器14c的、经增益调整的右变换23b以及经增益调整的、经谱整形的中心通道25，并且产生信号26b作为输出。

[0015] 反变换器18a、18b接收相应的信号26a、26b，并且产生相应推出的左右通道信号L′18a、R′18b。

[0016] 以下更详细地描述语音增强器1的操作。依次(首先大致概括，然后更详细地)描述根据一个实施例的中心通道提取、平谱、声音活动检测以及混合的过程。

[0017] 中心通道提取

[0018] 假设如下：

[0019] (1)感兴趣的信号17包含语音。

[0020] (2)在多通道信号(即，左和右，或者立体声)的情况下，该语音是中心平坦的(center panned)。

[0021] (3)真正平坦的中心包括源左和右信号的比例α。。

[0022] (4)减去该比例的结果是一对正交信号

[0023] 根据这些假设进行操作，中心通道提取器11从立体声信号17提取中心平坦的内容C 20。对于中心平坦的内容，左右通道二者的相同区域包含该中心平坦的内容。通过从左右通道二者移除相同的部分来提取该中心平坦的内容。

[0024] 可针对剩余左右信号(在多块的帧上或者使用随着新块输入连续地更新的方法)计算LR＊＝0(其中＊表示共轭)，并且调整比例α直至该量足够接近零。

[0025] 平谱

[0026] 听觉滤波器将假定的语音通道中的语音分成感知频带。确定每个数据块具有最大能量的频带。然后改变该块的语音通道的谱形状，以补偿剩余频带中较低的能量。使谱变平：具有较低能量的频带使其增益增加直至某个最大值。在一个实施例中，所有频带可共享最大增益。在替选实施例中，每个频带可具有其自己的最大增益。(在所有频带具有相同能量的退化情况下，则谱已经是平的。可认为谱整形没有发生，或者可认为用相同功能实现了谱整形。)

[0027] 平谱发生而与通道内容无关。可处理非语音，但是稍后在系统中不使用该非语音。与语音相比，非语音具有非常不同的谱，因此针对非语音的变平一般与针对语音的不同。

[0028] 声音活动检测器

[0029] 一旦假设的语音被隔离到单个通道，则针对语音内容对该语音进行分析。它包含语音吗？内容独立于平谱而被分析。语音内容通过测量相邻的数据帧中的谱波动而被确定。(每帧可包括许多数据块，但是在48kHZ采样速率，帧通常是两个、四个或者八个块。)[0030] 在从立体声提取语音通道的情况下，剩余的立体声信号可帮助语音分析。该概念更普遍地应用于任何多通道源中的相邻通道。

[0031] 混合

[0032] 当认为存在语音时，将变平的语音通道以某一比例与原始信号混合，该比例与语音通道实际包含语音的置信度有关。一般而言，当置信度高时，使用较多变平的语音通道。当置信度低时，使用较少变平的语音通道。

[0033] 依次更详细地描述根据一个实施例的中心通道提取、平谱、声音活动检测以及混合的过程。

[0034] 从2通道源提取假想中心和周围通道

[0035] 使用语音增强，期望仅对中心平坦的音频进行提取、处理以及再插入。在立体声混合中，语音通常是中心平坦的。

[0036] 现在描述从2通道混合中提取中心平坦的音频(假想中心通道)。数学证明构成第一部分。第二部分将该证明应用于真实的立体声信号，以推出假想中心。

[0037] 当从原始立体声减去假想中心时，剩下具有正交通道的立体声信号。类似的方法从周围平坦的音频推出假想周围通道。

[0038] 中心通道提取-数学证明

[0039] 给定某一两通道信号，可将通道分成左(L)和右(R)。左右通道均包含独有的信息以及公共的信息。可将公共的信息表示为C(中心平坦的)，并且将独有的信息分别表示为L和R(仅左和仅右)。

[0040] L＝L+C

[0041] R＝R+C (1)

[0042] “独有”意味着L和R彼此正交：

[0043] LR*＝0 (2)

[0044] 如果将L和R分为实部和虚部，则

[0045] LrRr+LiRi＝0 (3)

[0046] 其中Lr是L的实部，Li是L的虚部，并且对于R类似。

[0047] 现在假设通过从L和R减去中心平坦的C而从非正交对(L和R)产生正交对(L和R)。

[0048] L＝L-C (4)

[0049] R＝R-C (5)

[0050] 现在令C＝αC，其中C是假设的中心频道，并且α是比例因子：

[0051] L＝L-αC (6)

[0052] R＝R-αC (7)

[0053] 将方程(6)和(7)代入方程(3)：

[0054] LrRr+LiRi＝(Lr-αCr)(Rr-αCr)+(Li-αCi)(Ri-αCi)

[0055] ＝LrRr-αCr(Lr+Rr)+α2Cr2+LiRi-αCi(Li+Ri)+α2Ci2

[0056] ＝α2[Cr2+Ci2]+α[-Cr(Lr+Rr)-Ci(Li+Ri)]+[LrRr+LiRi]

[0057] ＝0

[0058] (8)

[0059] 方程(8)是二次方程形式：

[0060] α2X+αY+Z＝0 (9)

[0061] 其中通过以下得到根：

[0062]

[0063] 现在令方程(6)和(7)中假设的C如下：

[0064] C＝L+R (11)

[0065] 分为实部和虚部：

[0066] Cr＝Lr+Rr (12)

[0067] Ci＝Li+Ri (13)

[0068] 然后在二次方程(9)中：

[0069] X＝Cr2+Ci2＝(Lr+Rr)2+(Li+Ri)2 (14)

[0070] Y＝-Cr(Lr+Rr)-Ci(Li+Ri)＝-(Lr+Rr)2-(Li+Ri)2＝-X (15)

[0071] Z＝LrRr+LiRi (16)

[0072] 将方程(14)、(15)和(16)代入方程(10)，并且解出α：

[0073]

[0074]

[0075] 选择α的解的负根，并且将α限制到{0，0.5}的范围，避免与周围平坦的信息混淆(尽管该值对本发明不是关键的)。假想中心通道方程然后变为：

[0076]

[0077]

[0078] 其中

[0079]

[0080] (min{}和max{}函数将α限制到{0，0.5}的范围，尽管该值对本发明不是关键的。)

[0081] 可以类似地推出假想周围通道为：

[0082]

[0083]

[0084]

[0085] 其中S是原始立体声对(L，R)中的周围平坦的音频，并且S被假设为(L-R)。再次，选择β的解的负根，并且将β限制到{0，0.5}的范围，避免与中心平坦的信息混淆(尽管该值对本发明不是关键的)。

[0086] 现在推出了C和S，它们可以从原始立体声对(L和R)中被移除，以从原始的两个通道得到四个通道。

[0087] L′＝L-C-S (22)

[0088] R′＝R-C+S (23)

[0089] 其中L′是推出的左，C是推出的中心，R′是推出的右，并且S是推出的周围通道。

[0090] 中心通道提取-应用

[0091] 如上所声明的那样，对于语音增强方法，主要关注是中心频道的提取。在该部分中，将上述技术应用于音频信号的复杂频域表示。

[0092] 假想中心频道提取中的第一步骤是在音频样本块上进行DFT，并获得所得到的变换系数。DFT的块大小取决于采样速率。例如，在48kHZ的采样速率，具有N＝512个样本的块大小将是可接受的。开窗函数w[n](诸如汉明窗)对变换应用之前的样本块进行加权：

[0093] 0≤n＜N (24)

[0094] 其中n是整数，并且N是块中的样本数。

[0095] 方程(25)计算DFT系数为：

[0096]

[0097] 其中x[n，c]是块m的通道c中的样本数n，j是虚数单位(j2＝-1)，并且Xm[k，c]是针对块m中的样本的通道c中的变换系数k。注意通道数是三：左、右、和假想中心(在x[n，c]的情况下，仅左和右)。在以下方程中，左通道被指定为c＝1，假想中心为c＝2(还未推出)，并且右通道为c＝3。另外，快速傅里叶变换(FFT)可以高效地实现DFT。

[0098] 逐频谱单元(frequency bin)地得到左和右的和与差。对实部和虚部进行分组和平方。然后在计算α之前在块之间对每个频谱单元进行平滑。平滑减少了在频谱单元中的功率在数据块之间改变过快时出现的听觉赝像。平滑可通过例如漏积分器(leaky integrator)、非线性平滑器、线性但多极低通平滑器或者更精细的平滑器来完成。

[0099] Bm(k)diff＝(Re{Xm[k，1]}-Re{Xm[k，3]})2+(Im{Xm[k，1]}-Im{Xm[k，3]})2(26a)[0100] Bm(k)sum＝(Re{Xm[k，1]}+Re{Xm[k，3]})2+(Im{Xm[k，1]}+Im{Xm[k，3]})2(26b)[0101] Btemp＝λ1Bm-1(k)diff+(1-λ1)Bm(k)diff

[0102] Bm(k)diff＝Btemp 0＜＜λ1＜1 (26c)

[0103] Btemp＝λ1Bm-1(k)sum+(1-λ1)Bm(k)sum

[0104] Bm(k)diff＝Btemp 0＜＜λ1＜1 (26d)

[0105] 其中Re{}是实部，Im{}是虚部，并且λ1是漏积分器系数。漏积分器具有低通滤波效果，并且λ1的典型值是0.9。然后使用方程(19)推出块m的提取系数α：

[0106]

[0107] 然后使用方程(18)推出块m的假象中心通道：

[0108] Xm[k，2]＝αm(k)(Xm[k，1]+Xm[k，3]) (28)

[0109] 平谱

[0110] 随后是本发明的平谱的实施例的描述。假设主要是语音的单个通道，通过离散傅里叶变换(DFT)或者相关变换将语音信号变换到频域。然后通过对变换频谱单元进行平方将幅值谱变换为功率谱。

[0111] 然后可能在临界的或者听觉滤波器的标度上将频谱单元分组成频带。将语音信号划分为临界频带模仿人类听觉系统(特别是耳蜗)。这些滤波器呈现了近似圆形指数形状，并且在等效矩形带宽(ERB)标度上均匀地被隔开。ERB标度仅是在心理声学中使用的测量，其近似听觉滤波器的带宽和间距。图2示出了具有1ERB的间距从而导致共40个频带的适当滤波器集合。将音频数据分段还有助于消除在每个频谱单元基础上工作时可能发生的听觉赝像。然后，相对于时间对临界频带的功率进行平滑，也就是说，跨越相邻的块进行平滑。

[0112] 得到经平滑的临界频带中的最大功率，并且针对剩余(非最大)频带计算相应的增益，以使得它们的功率更接近最大功率。增益补偿类似于基膜(basilar membrane)的压缩(非线性)特性。这些增益被限制到最大值以避免饱和。为了将这些增益施加于原始信号，它们必须被变换回到DFT格式。因此，每个频带的功率增益首先被变换回到频谱单元功率增益，然后每个频谱单元的功率增益通过取每个频谱单元的平方根而被转换为幅值增益。原始信号变换频谱单元然后可以乘以算出的每个频谱单元的幅值增益。经谱变平的信号然后被从频域变换回到时域。在假想中心的情况下，其在被返回到时域之前首先与原始信号混合。图3描述了该过程。

[0113] 上述平谱系统没有考虑输入信号的特性。如果非语音信号被变平，则音质中感知的变化会很严重。为了避免非语音信号的处理，上述方法可以与声音活动检测器13耦合。当声音活动检测器13指示存在语音时，使用变平的语音。

[0114] 假设已经如上所述地将要变平的信号转换到频域。为了简单起见，省略了以上使用的通道标记。DFT系数被转换为功率，并且然后从DFT域转换到临界频带。

[0115] (29)

[0116] 0≤p＜P

[0117] 其中H[k，p]是P个临界带通滤波器。

[0118] 然后在块之间对每个频带中的功率进行平滑，类似于发生在大脑皮层的时间积分。平滑可通过例如漏积分器、非线性平滑器、线性但多极低通平滑器或者更精细的平滑器来完成。这种平滑还有助于消除可能引起增益在块之间波动过快从而引起可听见的泵送(pumping)的瞬时行为。然后得到峰值功率。

[0119] Em[p]＝λ2Em-1[p]+(1-λ2)Cm[p] 0＜＜λ2＜1 (30a)

[0120]

[0121] 其中Em[p]是平滑的、临界频带的功率，λ2是漏积分器系数，并且Emax是峰值功率。漏积分器有低通滤波的效果，另外，λ2的典型值是0.9。

[0122] 接下来得到每个频带的功率增益，其中限制了最大增益以避免过补偿：

[0123]

[0124] 0＜γ＜1 (31b)

[0125] 其中Gm[p]是要施加于每个频带的功率增益，Gmax是允许的最大功率增益，并且γ确定谱的平整(leveling)度。实际上，γ接近一。如果指定了执行处理的系统以及对增益量的任何其它全局限制，则Gmax取决于动态范围(或者净空(headroom))。Gmax的典型值是20dB。

[0126] 每个频带的功率增益接下来被转换为每个频谱单元的功率，并且取平方根来获得每个频谱单元的幅值增益：

[0127] (32)

[0128] 0≤k＜K

[0129] 其中Ym[k]是每个频谱单元的幅值增益。

[0130] 接下来基于声音活动检测器输出21、22修改幅值增益。接下来描述根据本发明的一个实施例的声音活动检测的方法。

[0131] 声音活动检测

[0132] 谱通量测量信号的功率谱改变的速度，比较相邻的音频帧之间的功率谱。(帧是多个音频数据块。)谱通量指示音频分类中的声音活动检测或者语音/非语音确定。通常，使用额外的指示器，并且集中(pool)结果来判决音频是否确实是语音。

[0133] 一般而言，语音的谱通量稍微高于音乐的谱通量，也就是说，与语音谱相比，音乐谱在帧之间倾向于更稳定。

[0134] 在立体声的情况下，当提取了假想中心通道时，DFT系数首先被分为中心和侧边音频(原始立体声减去假想中心)。这与传统的中间/侧边(mid/side)立体声处理的不同之处在于中间/侧边处理通常是(L+R)/2、(L-R)/2；而中心/侧边处理是C、L+R-2C。

[0135] 信号如上所述地被转换到频域，DFT系数被转换为功率并且然后从DFT域被转换到临界频带域。然后使用临界频带功率计算中心和侧边二者的谱通量：

[0136] (33a)

[0137] 0≤p＜P

[0138]

[0139] 0≤p＜P

[0140] 其中 是假想中心的临界频带型式， 是剩余信号(左右之和减去中心)的临界频带型式，并且H[k，p]是如之前所述的P个临界带通滤波器。

[0141] 根据之前的2J个数据块产生两个帧缓冲器(用于中心和侧边的幅值)：

[0142]

[0143]

[0144] 下一步骤根据当前和之前的帧的平均功率计算中心通道的权重W。这是在限定的频带范围上完成的：

[0145]

[0146] 频带的范围限于语音的主要带宽(近似100-8000Hz)。然后计算中心和侧边二者的未加权的谱通量：

[0147]

[0148]

[0149] 其中Fx(m)是中心的未加权的谱通量，并且Fs(m)是侧边的未加权的谱通量。

[0150] 然后如下计算谱通量的偏置估计：

[0151] 如果 并且W(m)＞Wmin (37a)

[0152]

[0153] 否则，

[0154] FTot(m)＝0 (37c)

[0155] 其中FTot(m)是总通量估计，并且Wmin是允许的最小权重。Wmin取决于动态范围，但是典型值将是Wmin＝-60dB。

[0156] 通过用简单1阶IIR(Infinite Impulse Response，无限脉冲响应)低通滤波器对FTot(m)的值进行低通滤波，计算谱通量的最终的、经平滑的值。该滤波器取决于信号的采样速率和块大小，但是在一个实施例中，对于fs＝48kHz，可以通过具有0.025＊fs的标准化截止的一阶低通滤波器来定义，其中fs是数字系统的采样速率。

[0157] FTot(m)然后被夹在0≤FTot(m)≤1的范围：

[0158] FTot(m)＝min{max{0.0，FTot(m)}，1.0} (38)

[0159] (根据该实施例，min{}和max{}函数将FTot(m)限制到{0，1}的范围。)[0160] 混合

[0161] 基于声音活动检测器的输出，将变平的中心通道与原始音频信号混合。

[0162] 将用于平谱的每个频谱单元的幅值增益Ym[k](如上所示)施加于假想中心通道Xm[k，2](如上所推出的那样)：

[0163] Xtemp＝Ym[k]Xm[k，2]

[0164] Xm[k，2]＝Xtemp (39)

[0165] 当声音活动检测器13检测到语音时，令FTot(t)＝1；当其检测到非语音时，令FTot(m)＝0。0和1之间的值是可能的，在该情况下声音活动检测器13对语音的存在做出软判决。

[0166] 对于左通道，

[0167] Xtemp＝(1-FTot(m))Xm[k，1]+FTot(m)Xm[k，2]

[0168] Xm[k，1]＝Xtemp

[0169] 0≤FTot(m)≤1 (40a)

[0170] 类似地，对于右通道，

[0171] Xtemp＝(1-FTot(m))Xm[k，3]+FTot(m)Xm[k，2]

[0172] Xm[k，3]＝Xtemp

[0173] 0≤FTot(m)≤1 (40b)

[0174] 实际上，FTot可被限于更窄的值范围。例如，0.1≤FTot(m)≤0.9保留了最终混合中变平的信号和原始信号两者中的小量。

[0175] 然后将每个频谱单元的幅值增益施加于原始输入信号，该原始输入信号然后经由反DFT而被转换回到时域：

[0176]

[0177] 0≤n＜N

[0178] c＝1，3 (41)

[0179] 其中 是原始立体声输入信号x的增强型式。

[0180] 图4示出了根据本发明的一个实施例的计算机4。计算机4包括存储器41、CPU42以及总线43。总线43通信地耦合存储器41和CPU 42。存储器41存储用于执行任何上述方法的计算机程序。

[0181] 已经描述了本发明的多个实施例。然而，本领域的普通技术人员应理解，如何不同地修改所述实施例而不背离本发明的精神和范围。例如，尽管描述中包括离散傅里叶变换，但是本领域的普通技术人员应理解，从时域变换到频域的各种可替选方法，反之亦然。现有技术

[0182] Schaub，A. 和 P.Straub，P.，“Spectral sharpening for speech enhancementnoise reduction”，Proc.ICASSP 1991，多伦多，加拿大，1991年5月，第993-996页。

[0183] Sondhi，M.，“New methods of pitch extraction，Audio andElectroacoustics”，IEEE Transactions，1968年6月，第16卷，第2期，第262-266页。

[0184] Villchur，E.，“Signal Processing to Improve Speech Intelligibility for theHearing Imparied”，第99次Audio Engineering Society Convention，1995年9月。

[0185] Thomas，I. 和 Niederjohn，R.，“Preprocessing of Speech for AddedIntelligibility in High Ambient Noise”， 第 34 次 Audio Engineering SocietyConvention，1968年3月。

[0186] Moore，B. 等 人，“A Model for the Prediction of Thresholds，Loudness，andPartial Loudness”，J.Audio Eng.Soc.，第 45 卷，第 4 章，1997 年 4 月。Moore，B.和Oxenham，A.，“Psychoacoustic consequences of compression inthe peripheral auditory system”，The Journal of the Acoustical Society ofAmerica，2002年12月，第112卷，第6期，第2962-2966页。

[0187] 现有技术-平谱

[0188] 美国专利

[0189] US 6732073 B1 Spectral enhancement of acoustic signals to provideimproved recognition of speech

[0190] US 06993480 B1 Voice intelligibility enhancement system

[0191] US 2006/0206320 A1 Apparatus and method for noise reduction and[0192] speech enhancement with microphones and loudspeakers

[0193] US 07191122 Speech compression system and method

[0194] US 2007/0094017 Frequency domain format enhancement

[0195] 国际专利

[0196] WO 2004/013840 A1 Digital Signal Processing Techniques For ImprovingAudio Clarity And Intelligibility

[0197] WO 2003/015082 Sound Intelligibility Enhancement Using APsychoacoustic Model And An Oversampled Filterbank

[0198] 论文

[0199] Sallberg，B.等人；“Analog Circuit Implementation for SpeechEnhancement Purposes Signals”；Systems and Computers，2004。

[0200] Conference Record of the Thirty-Eighth Asilomar Conference。

[0201] Magotra，N. 和 Sirivara，S.；“Real-time digital speech processing strategiesfor the hearing impaired”；Acoustics，Speech，and Signal Processing，1997.ICASSP-97.，1997，第1211-1214页，第2卷。

[0202] Walker，G.，Byrne，D. 和 Dillon，H. ；“The effects ofmultichannelcompression/expansion amplification on the intelligibility of nonsensesyllables in noise”；The Journal of the Acoustical Society of America，
1984年9月，第76卷，第3期，第746-757页。

[0203] 现有技术-中心提取

[0204] AdobeAudition具有声音/乐器提取功能

[0205] http://www.adobeforums.com/cig-bin/webx/.3bc3a3e5

[0206] Winamp的“center cut”

[0207] http://www.hydrogenaudio.org/forums/lofiversion/index.php/t17450.html[0208] 现有技术-谱通量

[0209] Vinton，M. 和 Robinson C；“Automated Speech/Other Discrimination forLoudness Monitoring”，第118次AES Convention，2005。

[0210] Scheirer E. 和 Slaney M.，“Construction and evaluation of a robustmultifeature speech/music discriminator”，IEEE Transactions onAcoustics，Speech，and Signal Processing(ICASSP′97)，1997，第1331-1334页。