CN113542980B - 一种抑制扬声器串扰的方法 - Google Patents

一种抑制扬声器串扰的方法 Download PDF

Info

Publication number
CN113542980B
CN113542980B CN202110825596.9A CN202110825596A CN113542980B CN 113542980 B CN113542980 B CN 113542980B CN 202110825596 A CN202110825596 A CN 202110825596A CN 113542980 B CN113542980 B CN 113542980B
Authority
CN
China
Prior art keywords
frequency
loudspeaker
crosstalk
amplitude
signal
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.)
Active
Application number
CN202110825596.9A
Other languages
English (en)
Other versions
CN113542980A (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.)
Shenzhen Yueersheng Acoustics Co ltd
Original Assignee
Shenzhen Yueersheng Acoustics Co ltd
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 Shenzhen Yueersheng Acoustics Co ltd filed Critical Shenzhen Yueersheng Acoustics Co ltd
Priority to CN202110825596.9A priority Critical patent/CN113542980B/zh
Publication of CN113542980A publication Critical patent/CN113542980A/zh
Application granted granted Critical
Publication of CN113542980B publication Critical patent/CN113542980B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R3/00—Circuits for transducers, loudspeakers or microphones
    • H04R3/02—Circuits for transducers, loudspeakers or microphones for preventing acoustic reaction, i.e. acoustic oscillatory feedback
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208—Noise filtering
    • G10L2021/02082—Noise filtering the noise being echo, reverberation of the speech

Landscapes

  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Human Computer Interaction (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Quality & Reliability (AREA)
  • Computational Linguistics (AREA)
  • Multimedia (AREA)
  • General Health & Medical Sciences (AREA)
  • Otolaryngology (AREA)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

本发明公开了一种抑制扬声器串扰的方法，属于音频信号处理技术领域。本发明采用基于双谱重构的方法，即综合利用声音传播通道的幅频响应和相频响应特性，在麦克风拾取通道对扬声器原信号进行重构，再对消；能很好在对拾取到的串扰信号进行抵消，既包括直达信号，也包括一定程度的多径传播的混响信号。本发明所采用的基于重叠相加的时频域变换方法，减少滤波处理的延时，延时时间大幅缩短。同时，本发明采用基于多项式拟合获取频谱响应系数，解决了耳机腔体产生的回声过程中，由于多次折射和反射等复杂的过程导致的幅频响应特性和相频响应特性随声压大小呈非线性特性的问题。

Description

一种抑制扬声器串扰的方法
技术领域
本发明属于音频信号处理技术领域，尤其涉及一种抑制扬声器串扰的方法。
背景技术
耳机或其它小功率通话设备，由于设计需要或是结构限制，存在扬声器与拾音麦克风串扰的问题，即扬声器播放的声音又被本机的麦克风拾取后回传给对方，形成回声情况。
目前解决串扰的方法，一是基于IIR型滤波器，另一种是基于自适应FIR型滤波器方法。其中，IIR型滤波器，方法实现简单，可在时域上直接滤波输出。但由于IIR型滤波器只反映了声音传输过程中幅频响应的特性，不能反映相频特性的变化，还需加入单独的延时处理环节，增加了复杂度，对于混响也无能为力。而自适应FIR型滤波器方法，采用远端扬声器信号与近端麦克风采集信号相逼近的自适应方式实现滤波，FIR滤波器具有线性相位和稳定性，可以在一定程度上反映幅频响应和相频响应的变化，对混响也有抑制作用。但这一算法对噪声较为敏感，在信噪比较低或噪声始终存在的场合，收敛过程会发散，不能起到很好地滤波作用。
另外，现有技术通用的DSP(Digital Signal Processing，数字信号处理)来说，为保证处理能力和算力瓶颈问题，一般采用2*block(数据块)的延时策略，即用一个数据块的时长去接收数据，另一块的时长内作处理，在第三个输入数据块到来时，将之前处理过的数据块同步输出；这样就会有2个数据块的延时，在fft(fast Fourier transform，快速傅立叶变换)快速傅里叶时频处理中，时间分辨率与频率分辨率是相互矛盾的；所以语音处理一般选择10ms-30ms的数据块作为处理单元，以保证短时稳定特性和长时的时变特性。即使选择10ms的数据块，最终输出也会有20ms的延时；这对延时要求苛刻的场合很难满足要求。
发明内容
基于此，本发明实施例的主要目的在于提供一种抑制扬声器串扰的方法，旨在解决通过IIR型滤波器或自适应FIR型滤波器处理串扰存在的技术问题；采用基于双谱重构的方法，即综合利用声音传播通道的幅频响应和相频响应特性，在麦克风拾取通道对扬声器原信号进行重构，再对消；能很好在对拾取到的串扰信号进行抵消，既包括直达信号，也包括一定程度的多径传播的混响信号，而且延时较小。
本发明实施例提供一种抑制扬声器串扰的方法，在麦克风拾取通道对扬声器的原始信号进行重构，得到重构信号；通过所述重构信号对在所述麦克风拾取通道拾取到的串扰信号进行抵消；所述串扰信号包括直达信号或/和混响信号。
本发明提供的一种抑制扬声器串扰的方法，包括以下步骤：
(1)基于重叠相加(OLA，overlap add)的时频域变换方法，将扬声器的音频输入信号sin，麦克风的输入信号min，分别以10-30ms的数据块为单位，将10ms-30ms的输入数据块分成若干子块，作m(fast Fourier transform，快速傅立叶变换)变换，将信号变换到频域；由功率谱获得数据块时间段内的声压值，得到声压级别序数i，基于声压级别序数i得到扬声器端输入信号的幅频与相频为sam[k][i]和
Figure GDA0004087340310000021
麦克风的输入信号min，经fft后，其幅频与相频为mam(k)和/>
Figure GDA0004087340310000022
其中k为第k个频点；
(2)基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A，求出音频输入信号sin对应的由串扰引起的幅频响应sre_am[k][i]和相频响应
Figure GDA0004087340310000023
(3)音频信号恢复；采用谱减法对麦克风的输入信号min进行恢复；采用如式I所示的复指数形式获得输出频域：
Figure GDA0004087340310000031
其中：mout为麦克风输入信号min经串扰处理后的时域信号；
Figure GDA0004087340310000032
为第k个频点的频域；
j为虚数单位；
将获得的频域
Figure GDA0004087340310000033
变换到时域/>
Figure GDA0004087340310000034
按照步骤(1)中的基于重叠相加的时频域变换方法，将其最后的子块输出，得到最终的麦克风端的恢复信号。
步骤(1)中所述的基于重叠相加的时频域变换方法，具体步骤如下：
将输入信号数据分为10-30ms的数据块；每次作fft处理的是10ms-30ms的数据块，将10ms-30ms的输入数据块分成若干子块(sub_block)，每次将最新输入的子块数据添加到fft的数据序列中，将最后的子块抛弃；这样每次作fft时，会有3个子块的重叠；ifft处理后输出时，每次只输出处理后的最后一个子块的数据，即最新的1组数据；这样，最终的延时是2倍子块数据的时长，即2*sub_block；减少滤波处理的延时。
步骤(2)中所述的基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A，通过下述方式实现：
1)通过实验测试获取扬声器原输入音频信号在不同声压级别序数i下，其在频点k的原频率幅值am[i]，及其对应的由串扰引起的的幅频响应re_am[i]；原频率相位Φ[i]，及其对应的串扰引起的相频响应re_Φ[i]；并组成如下幅频响应对、相频响应对：
幅频响应对：{am[k][i]，re_am[k][i]}；其中，k为信号频点序号；i为声压级别序数；
相频响应对：{Φ[k][i]，re_Φ[k][i]}；其中，k为信号频点序号；i为声压级别序数；
2)建立多项式方程，求解多项式系数：
y1＝a0+a1*x1+a2*x1^2+a3*x1^3+…+an*x1^n；
y2＝a0+a1*x2+a2*x2^2+a3*x2^3+…+an*x2^n；
……
yi＝a0+a1*xi+a2*xi^2+a3*xi^3+…+an*xi^n；
上式[xi，yi]，表示幅频响应对，或相频响应对；
i为声压级别序数，i∈[1，I]，I为整数；
n为多项式次数，n∈[1，N]，N为整数；
将上式多项式组写成如下矩阵形式：
Y＝AHX；
3)采用基于最小二乘法和最陡梯度下降法求解系数。
作为优选的实施方式，步骤2)中
考虑小功率扬声器混响或回声相对较小，其非线性的复杂性相对低些，减少计算的复杂度，作为优选的实施方式，可设所述的N＝3。
步骤(3)中所述的采用基于最小二乘法和最陡梯度下降法求解系数，具体为：
以最小二乘法构建目标函数如式II所示：
Figure GDA0004087340310000041
采用最陡梯度下降及迭代方法求出A；
Ap＝Ap-1+λ*g；
其中：
p为第p次迭代；
λ为迭代步长；
g为第k次迭代的梯度矩阵；
其中第n个系数的梯度值由以下公式得到：
Figure GDA0004087340310000051
当式II满足阈值时，Ap即为所求系数矩阵。
步骤3)中，将所述的式I由欧拉公式展开：
Figure GDA0004087340310000052
步骤3)中，所述的阈值为迭代结束的条件，为迭代次数阈值或均方差阈值中先达到的一种，迭代即终止。具体的，所述迭代次数阈值，为大于设定的迭代次数阈值，迭代随即终止；所述均方差阈值，为小于设定的均方误差阈值(如0.1)，迭代随即终止。
本发明与现有技术相比，具有以下的优点和有益效果：
1.本发明采用基于双谱重构的方法，即综合利用声音传播通道的幅频响应和相频响应特性，在麦克风拾取通道对扬声器原信号进行重构，再对消；能很好在对拾取到的串扰信号进行抵消，既包括直达信号，也包括一定程度的多径传播的混响信号,而且延时较小。
2.本发明所采用的基于重叠相加的时频域变换方法，将现有技术中2个数据块的延迟减小到2个子块的延迟，减少滤波处理的延时，延时时间大幅缩短。
3.本发明采用基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A的方法，对耳机腔体回声这一复杂过程作拟合逼近；解决了耳机腔体产生的回声过程中，由于多次折射和反射等复杂的过程导致的幅频响应特性和相频响应特性随声压大小呈非线性特性的问题。
附图说明
图1为本发明一实施例具有录音功能的蓝牙耳机的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
需要说明，若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后、顶、底……)，则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
目前，常见的解决串扰的方法，一种是基于IIR(Infinite Impulse Response，无限冲激响应)型滤波器进行解决，另一种是基于自适应FIR(Finite Impulse Response，有限冲激响应)型滤波器进行解决。但是，由于IIR型滤波器只反映声音传输过程中幅频响应的特性，不能反映相频响应特性的变化，还需加入单独的延时处理环节，增加了处理的复杂度，也不能处理混响的串扰问题。但是，通过自适应FIR型滤波器解决串扰的算法对噪声较为敏感，在信噪比较低或噪声始终存在的场合，收敛过程会发散，不能起到很好的滤波作用。
基于现有技术的缺陷，本发明实施例的主要目的在于提供一种抑制扬声器串扰的方法，旨在解决通过IIR型滤波器或自适应FIR型滤波器处理串扰存在的技术问题；采用基于双谱重构的方法，即综合利用声音传播通道的幅频响应和相频响应特性，在麦克风拾取通道对扬声器原信号进行重构，再对消；能很好在对拾取到的串扰信号进行抵消，既包括直达信号，也包括一定程度的多径传播的混响信号,而且延时较小。
本发明实施例提供一种抑制扬声器串扰的方法，在麦克风拾取通道对扬声器的原始信号进行重构，得到重构信号；通过所述重构信号对在所述麦克风拾取通道拾取到的串扰信号进行抵消；所述串扰信号包括直达信号或/和混响信号。
本发明提供的一种抑制扬声器串扰的方法，包括以下步骤：
(1)基于重叠相加的时频域变换方法，将扬声器的音频输入信号sin，麦克风的输入信号min，分别以10-30ms的数据块为单位，作fft(fast Fourier transform，快速傅立叶变换)变换，将信号变换到频域；由功率谱获得数据块时间段内的声压值，得到声压级别序数i，基于声压级别序数i得到扬声器端输入信号的幅频与相频为sam[k][i]和
Figure GDA0004087340310000073
麦克风的输入信号min，经fft后，其幅频与相频为mam(k)和/>
Figure GDA0004087340310000074
其中k为第k个频点；
(2)基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A，求出音频输入信号sin对应的由串扰引起的幅频响应sre_am[k][i]和相频响应
Figure GDA0004087340310000075
(3)音频信号恢复；采用谱减法对麦克风的输入信号min进行恢复；采用如式I所示的复指数形式获得输出频域：
Figure GDA0004087340310000071
其中：mout为麦克风输入信号min经串扰处理后的时域信号；
Figure GDA0004087340310000072
为第k个频点的频域；
j为虚数单位；
将获得的频域
Figure GDA0004087340310000077
变换到时域/>
Figure GDA0004087340310000076
按照步骤(1)中的基于重叠相加的时频域变换方法，将其最后的子块输出，得到最终的麦克风端的恢复信号。
现有技术通用的DSP(Digital Signal Processing，数字信号处理)来说，为保证处理能力和算力瓶颈问题，一般采用2*block(数据块)的延时策略，即用一个数据块的时长去接收数据，另一块的时长内作处理，在第三个输入数据块到来时，将之前处理过的数据块同步输出；这样就会有2个数据块的延时。
在fft(fast Fourier transform，快速傅立叶变换)快速傅里叶时频处理中，时间分辨率与频率分辨率是相互矛盾的；所以语音处理一般选择10ms-30ms的数据块作为处理单元，以保证短时稳定特性和长时的时变特性。即使选择10ms的数据块，最终输出也会有20ms的延时；这对延时要求苛刻的场合很难满足要求。
步骤(1)中所述的基于重叠相加的时频域变换方法，具体步骤如下：
将输入信号数据分为10-30ms的数据块；每次作fft处理的是10ms-30ms的数据块，将10ms-30ms的输入数据块分成若干子块(sub_block)，每次将最新输入的子块数据添加到fft的数据序列中，将最后的子块抛弃；这样每次作fft时，会有3个子块的重叠；ifft处理后输出时，每次只输出处理后的最后一个子块的数据，即最新的1组数据；这样，最终的延时是2倍子块数据的时长，即2*sub_block；减少滤波处理的延时；
如图1所示，将10ms的数据块分成4个子块，第一次参与fft的是数据块block1(输入子块0，1，2，3)；ifft后只选择子块3’输出；3’比输入子块3延时2个子块的时长。
以采样率16khZ为例，每个子块2.5ms(40个采样点)；最终时延为5ms(2*2.5ms)，如按照现有技术的2*block延时，则是20ms；本发明方法减少滤波处理的延时，延时时间大幅缩短。
对于步骤(2)中；耳机腔体产生的回声过程，包括了多次折射和反射等复杂的过程，其幅频响应特性和相频响应特性随声压大小呈非线性特性。为解决这一问题，本发明提出了基于多项式拟合的方法，对这一复杂过程作拟合逼近。
步骤(2)中所述的基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A，通过下述方式实现：
1)通过实验，获取扬声器原输入音频信号在不同声压级别序数i下，其在频点k的原频率幅值am[i]，及其对应的由串扰引起的的幅频响应re_am[i]；原频率相位Ф[i]，及其对应的串扰引起的相频响应re_Ф[i]；并组成如下幅频响应对、相频响应对：
幅频响应对：{am[k][i],re_am[k][i]}；其中，k为信号频点序号；i为声压级别序数；
相频响应对：{Ф[k][i],re_Ф[k][i]}；其中，k为信号频点序号；i为声压级别序数；
2)建立多项式方程，求解多项式系数：
y1＝a0+a1*x1+a2*x1^2+a3*x1^3+…+an*x1^n；
y2＝a0+a1*x2+a2*x2^2+a3*x2^3+…+an*x2^n；
……
yi＝a0+a1*xi+a2*xi^2+a3*xi^3+…+an*xi^n；
上式[xi,yi]，表示幅频响应对，或相频响应对；
i为声压级别序数，i∈[1，I]，I为整数；
n为多项式次数，n∈[1，N]，N为整数；考虑小功率扬声器混响或回声相对较小，其非线性的复杂性相对低些，可设N＝3；以减少计算的复杂度；
将上式多项式组写成如下矩阵形式：
Y＝AHX；
3)采用基于最小二乘法和最陡梯度下降法求解系数。
步骤3)中所述的采用基于最小二乘法和最陡梯度下降法求解系数，具体为：
以最小二乘法构建目标函数如式II所示：
Figure GDA0004087340310000091
采用最陡梯度下降及迭代方法求出A；
Ap＝Ap-1+λ*g；
其中：
p为第p次迭代；
λ为迭代步长；
g为第k次迭代的梯度矩阵；
其中第n个系数的梯度值由以下公式得到：
Figure GDA0004087340310000092
当式II满足阈值时，Ap即为所求系数矩阵。
步骤(3)中，将式I由欧拉公式展开：
Figure GDA0004087340310000101
步骤3)中，所述的阈值为迭代结束的条件，为迭代次数阈值或均方差阈值中先达到的一种，迭代即终止。具体的，所述迭代次数阈值，为大于设定的迭代次数阈值，迭代随即终止；所述均方差阈值，为小于设定的均方误差阈值(如0.1)，迭代随即终止。
本发明与现有技术相比，具有以下的有益效果：
1.本发明采用基于双谱重构的方法，即综合利用声音传播通道的幅频响应和相频响应特性，在麦克风拾取通道对扬声器原信号进行重构，再对消；能很好在对拾取到的串扰信号进行抵消，既包括直达信号，也包括一定程度的多径传播的混响信号,而且延时较小。
2.本发明所采用的基于重叠相加的时频域变换方法，将现有技术中2个数据块的延迟减小到2个子块的延迟，减少滤波处理的延时，延时时间大幅缩短。
3.本发明采用基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A的方法，对耳机腔体回声这一复杂过程作拟合逼近；解决了耳机腔体产生的回声过程中，由于多次折射和反射等复杂的过程导致的幅频响应特性和相频响应特性随声压大小呈非线性特性的问题。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种抑制扬声器串扰的方法，其特征在于：包括以下步骤：
(1)基于重叠相加的时频域变换方法，将扬声器的音频输入信号sin，麦克风的输入信号min，分别以10-30ms的数据块为单位，将10ms-30ms的输入数据块分成若干子块，作fft变换，将信号变换到频域；由功率谱获得数据块时间段内的声压值，得到声压级别序数i，基于声压级别序数i得到扬声器端输入信号的幅频与相频为sam[k][i]和
Figure FDA0004087340300000011
麦克风的输入信号min，经fft后，其幅频与相频为mam(k)和/>
Figure FDA0004087340300000012
其中k为第k个频点；
(2)基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A，求出音频输入信号sin对应的由串扰引起的幅频响应sre_am[k][i]和相频响应
Figure FDA0004087340300000013
(3)音频信号恢复；采用谱减法对麦克风的输入信号min进行恢复；采用如式I所示的复指数形式获得输出频域：
Figure FDA0004087340300000014
其中：mout为麦克风输入信号min经串扰处理后的时域信号；
Figure FDA0004087340300000015
为第k个频点的频域；
j为虚数单位；
将获得的频域
Figure FDA0004087340300000016
变换到时域/>
Figure FDA0004087340300000017
按照步骤(1)中的基于重叠相加的时频域变换方法，将其最后的子块输出，得到最终的麦克风端的恢复信号。
2.根据权利要求1所述的抑制扬声器串扰的方法，其特征在于：
步骤(1)中所述的基于重叠相加的时频域变换方法，具体步骤如下：
将输入信号数据分为10-30ms的数据块；每次作fft处理的是10ms-30ms的数据块，将10ms-30ms的输入数据块分成若干子块，每次将最新输入的子块数据添加到fft的数据序列中，将最后的子块抛弃；这样每次作fft时，会有3个子块的重叠；ifft处理后输出时，每次只输出处理后的最后一个子块的数据。
3.根据权利要求1所述的抑制扬声器串扰的方法，其特征在于：
步骤(2)中所述的基于多项式拟合获取频谱响应系数，获得频谱响应系数矩阵A，通过下述方式实现：
1)通过实验测试获取扬声器原输入音频信号在不同声压级别序数i下，其在频点k的原频率幅值am[i]，及其对应的由串扰引起的的幅频响应re_am[i]；原频率相位Ф[i]，及其对应的串扰引起的相频响应re_Ф[i]；并组成如下幅频响应对、相频响应对：
幅频响应对：{am[k][i],re_am[k][i]}；其中，k为信号频点序号；i为声压级别序数；
相频响应对：{Ф[k][i],re_Ф[k][i]}；其中，k为信号频点序号；i为声压级别序数；
2)建立多项式方程，求解多项式系数：
y1＝a0+a1*x1+a2*x1^2+a3*x1^3+…+an*x1^n；
y2＝a0+a1*x2+a2*x2^2+a3*x2^3+…+an*x2^n；
……
yi＝a0+a1*xi+a2*xi^2+a3*xi^3+…+an*xi^n；
上式[xi,yi]，表示幅频响应对，或相频响应对；
i为声压级别序数，i∈[1，I]，I为整数；
n为多项式次数，n∈[1，N]，N为整数；
将上式多项式组写成如下矩阵形式：
Y＝AHX；
3)采用基于最小二乘法和最陡梯度下降法求解系数。
4.根据权利要求3所述的抑制扬声器串扰的方法，其特征在于：所述的N＝3。
5.根据权利要求3所述的抑制扬声器串扰的方法，其特征在于：步骤3)中所述的采用基于最小二乘法和最陡梯度下降法求解系数，具体为：
以最小二乘法构建目标函数如式II所示：
Figure FDA0004087340300000031
采用最陡梯度下降及迭代方法求出A；
Ap＝Ap-1+λ*g；
其中：
p为第p次迭代；
λ为迭代步长；
g为第k次迭代的梯度矩阵；
其中第n个系数的梯度值由以下公式得到：
Figure FDA0004087340300000032
当式II满足阈值时，Ap即为所求系数矩阵。
6.根据权利要求5所述的抑制扬声器串扰的方法，其特征在于：步骤(3)中，将所述的式I由欧拉公式展开：
Figure FDA0004087340300000033
7.根据权利要求5所述的抑制扬声器串扰的方法，其特征在于：步骤3)中所述的阈值为迭代结束的条件，为迭代次数阈值或均方差阈值中先达到的一种，迭代即终止。
CN202110825596.9A 2021-07-21 2021-07-21 一种抑制扬声器串扰的方法 Active CN113542980B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110825596.9A CN113542980B (zh) 2021-07-21 2021-07-21 一种抑制扬声器串扰的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110825596.9A CN113542980B (zh) 2021-07-21 2021-07-21 一种抑制扬声器串扰的方法

Publications (2)

Publication Number Publication Date
CN113542980A CN113542980A (zh) 2021-10-22
CN113542980B true CN113542980B (zh) 2023-03-31

Family

ID=78100813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110825596.9A Active CN113542980B (zh) 2021-07-21 2021-07-21 一种抑制扬声器串扰的方法

Country Status (1)

Country Link
CN (1) CN113542980B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101346895A (zh) * 2005-10-26 2009-01-14 日本电气株式会社 回声抑制方法及设备
JP2015037239A (ja) * 2013-08-13 2015-02-23 日本電信電話株式会社 残響抑圧装置とその方法と、プログラムとその記録媒体
CN108376548A (zh) * 2018-01-16 2018-08-07 厦门亿联网络技术股份有限公司 一种基于麦克风阵列的回声消除方法与系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101346896B (zh) * 2005-10-26 2012-09-05 日本电气株式会社 回声抑制方法及设备
JPWO2009051197A1 (ja) * 2007-10-19 2011-03-03 日本電気株式会社 エコー抑圧方法及び装置
EP3296992B1 (en) * 2008-03-20 2021-09-22 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus and method for modifying a parameterized representation
EP2562751B1 (en) * 2011-08-22 2014-06-11 Svox AG Temporal interpolation of adjacent spectra
CN103067821B (zh) * 2012-12-12 2015-03-11 歌尔声学股份有限公司 一种基于双麦克的语音混响消减方法和装置
CN111933166B (zh) * 2020-07-31 2024-01-26 广州视源电子科技股份有限公司 回声抑制方法、装置、存储介质和电子设备
CN112147590B (zh) * 2020-08-25 2023-04-07 西安电子科技大学 一种基于响应估计频域拟合的通道均衡方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101346895A (zh) * 2005-10-26 2009-01-14 日本电气株式会社 回声抑制方法及设备
JP2015037239A (ja) * 2013-08-13 2015-02-23 日本電信電話株式会社 残響抑圧装置とその方法と、プログラムとその記録媒体
CN108376548A (zh) * 2018-01-16 2018-08-07 厦门亿联网络技术股份有限公司 一种基于麦克风阵列的回声消除方法与系统

Also Published As

Publication number Publication date
CN113542980A (zh) 2021-10-22

Similar Documents

Publication Publication Date Title
CN109727604B (zh) 用于语音识别前端的频域回声消除方法及计算机储存介质
US7620172B2 (en) Method and system for eliminating noises and echo in voice signals
US8000482B2 (en) Microphone array processing system for noisy multipath environments
CN105825864B (zh) 基于过零率指标的双端说话检测与回声消除方法
CN105575397B (zh) 语音降噪方法及语音采集设备
CN104883462B (zh) 一种用于消除声学回声的自适应滤波器及滤波方法
CN107636758A (zh) 声学回声消除系统和方法
CN104835503A (zh) 一种改进gsc自适应语音增强方法
WO2012142270A1 (en) Systems, methods, apparatus, and computer readable media for equalization
CN110956975B (zh) 回声消除方法及装置
CN105810202B (zh) 一种降低回声的方法、装置及通讯设备
CN111213359B (zh) 回声消除器和用于回声消除器的方法
CN110211602B (zh) 智能语音增强通信方法及装置
CN104410761B (zh) 一种仿射投影符号子带凸组合自适应回声消除方法
CN115132215A (zh) 一种单通道语音增强方法
CN106024001A (zh) 一种提高麦克风阵列语音增强性能的方法
CN113542980B (zh) 一种抑制扬声器串扰的方法
JP2003188776A (ja) 音響エコー消去方法、装置及び音響エコー消去プログラム
JP3756828B2 (ja) 反響消去方法、この方法を実施する装置、プログラムおよびその記録媒体
US20050122893A1 (en) Techniques to reduce echo
JP2003309493A (ja) 反響低減方法、反響低減装置、反響低減プログラム
CN115620737A (zh) 语音信号处理装置、方法、电子设备和扩音系统
JP4396449B2 (ja) 残響除去方法及びその装置
CN112151051A (zh) 音频数据的处理方法和装置及存储介质
CN112489669A (zh) 一种音频信号处理方法、装置、设备和介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant