CN106601266A - 回声消除方法、装置及*** - Google Patents
回声消除方法、装置及*** Download PDFInfo
- Publication number
- CN106601266A CN106601266A CN201710106188.1A CN201710106188A CN106601266A CN 106601266 A CN106601266 A CN 106601266A CN 201710106188 A CN201710106188 A CN 201710106188A CN 106601266 A CN106601266 A CN 106601266A
- Authority
- CN
- China
- Prior art keywords
- subband
- signal
- reference signal
- echo
- primary input
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 57
- 230000008569 process Effects 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims description 30
- 230000000644 propagated effect Effects 0.000 claims description 19
- 230000003044 adaptive effect Effects 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 238000003786 synthesis reaction Methods 0.000 claims description 7
- 238000009499 grossing Methods 0.000 claims description 2
- 230000001902 propagating effect Effects 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000001914 filtration Methods 0.000 description 6
- 230000002093 peripheral effect Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000002592 echocardiography Methods 0.000 description 3
- 235000008434 ginseng Nutrition 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 241000209140 Triticum Species 0.000 description 2
- 235000021307 Triticum Nutrition 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000004549 pulsed laser deposition Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Speech or voice signal processing techniques 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
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明实施例提供的回声消除方法、装置及***，属于语音信号处理技术领域。所述方法包括：根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。提高回声消除过程中的计算效率。
Description
技术领域
本发明涉及语音信号处理技术领域，具体而言，涉及回声消除方法、装置及***。
背景技术
传统回声消除技术主要包含线性回声消除和残留回声消除两个模块。主说话人的语音、扬声器产生的回声和环境噪声被麦克风采集，通过回声消除去除线性回声，然后通过残留回声抑制模块消除残留回声。线性回声消除技术主要使用频域分块自适应滤波处理或子带自适应滤波处理，在回声消除过程中，一般需要双讲检测和噪声估计来解决外界干扰问题，增加***负担。并且存在回声路径较长引起的计算量较大；在外界语音和噪声干扰下，自适应滤波器不稳定性和收敛速度慢的问题。
发明内容
有鉴于此，本发明实施例的目的在于提供一种回声消除方法、装置及***。为了实现上述目的，本发明实施例采取的技术方案如下：
第一方面，本发明实施例提供了一种回声消除方法，所述方法包括：根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。
第二方面，本发明实施例提供了一种回声消除装置，所述装置包括：第一分析滤波器组单元，用于根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；第二分析滤波器组单元，用于根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；声学回声消除模块，用于基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。
第三方面，本发明实施例提供了一种回声消除***，所述***包括：电子设备、麦克风以及扬声器，所述电子设备包括上述的装置，所述电子设备分别与所述麦克风、所述扬声器耦合。
本发明实施例提供的回声消除方法、装置及***，所述方法包括：根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。提高回声消除过程中的计算效率。
本发明的其他特征和优点将在随后的说明书阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明实施例了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为一种可应用于本申请实施例中的回声消除***的结构框图；
图2为本发明第一实施例提供的回声消除方法的流程图；
图3为本发明第二实施例提供的回声消除装置的结构框图；
图4为本发明第二实施例提供的回声消除装置的一种应用环境示意图；
图5为本发明第二实施例提供的回声消除装置的另一种应用环境示意图；
图6为本发明第二实施例提供的回声消除装置的再一种应用环境示意图。
图中：100-电子设备；110-存储器；111-存储控制器；112-处理器；113-外设接口；200-回声消除***；201-麦克风；202-扬声器；300-装置；310-第一分析滤波器组单元；320-第二分析滤波器组单元；330-声学回声消除单元；340-综合滤波器组单元；350-残留回声抑制单元；351-滤波器组延时补偿子单元；352-残留回声抑制子单元。
具体实施方式
下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
请参阅图1，图1示出了一种可应用于本申请实施例中的回声消除***200的结构框图。所述回声消除***200包括电子设备100、麦克风201、扬声器202。所述电子设备100分别与麦克风201、扬声器202耦合。该电子设备100可以作为用户终端，也可以是计算机或服务器，所述用户终端可以为手机或平板电脑。如图1所示，电子设备100可以包括存储器110、存储控制器111、处理器112、外设接口113和回声消除装置。所述外设接口113分别与所述麦克风201、扬声器202电连接。
存储器110、存储控制器111、处理器112、外设接口113各元件之间直接或间接地电连接，以实现数据的传输或交互。例如，这些元件之间可以通过一条或多条通讯总线或信号总线实现电连接。回声消除方法分别包括至少一个可以以软件或固件(firmware)的形式存储于存储器110中的软件功能模块，例如所述回声消除装置包括的软件功能模块或计算机程序。
存储器110可以存储各种软件程序以及模块，如本申请实施例提供的回声消除方法及装置对应的程序指令/模块。处理器112通过运行存储在存储器110中的软件程序以及模块，从而执行各种功能应用以及数据处理，即实现本申请实施例中的回声消除方法。存储器110可以包括但不限于随机存取存储器(Random Access Memory，RAM)，只读存储器(ReadOnly Memory，ROM)，可编程只读存储器(Programmable Read-Only Memory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。
处理器112可以是一种集成电路芯片，具有信号处理能力。上述处理器可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(NetworkProcessor，简称NP)等；还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。其可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
外设接口113将各种输入/输出装置耦合至处理器112以及存储器110。在本实施例中，麦克风201和扬声器202通过外设接口113耦合至处理器112以及存储器110。
麦克风201用于采集用户的语音、扬声器202产生的回声和环境噪声；扬声器202用于输出用户的语音。
第一实施例
请参阅图2，本发明实施例提供了一种回声消除方法，所述方法包括：
步骤S200：根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；
步骤S210：根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；
扬声器202产生参考信号，参考信号设为x(n)；麦克风201采集主输入信号，主输入信号设为p(n)，p(n)＝s(n)+d(n)+b(n)，其中，d(n)为x(n)经过空间传播的回声信号，即d(n)定义为所述参考信号经过第一延时传播的回声，第一延时为K，s(n)为语音信号，b(n)为噪声信号。n为离散时间标度，n＝0,1,2……。
作为一种实施方式，对获取到的参考信号通过分析滤波器组，获得N个子带参考信号。同理可得，对获取到的主输入信号进行通过分析滤波器组，获得N个主输入子带信号。所述分析滤波器组包括复指数调制分析滤波器组或余弦变换分析滤波器组。
步骤S220：基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声；
所述N个子带参考信号对应有N个子带自适应滤波器，每个子带自适应滤波器的长度为L为全带自适应滤波器长度，所述基于所述N个子带参考信号，获得N个子带估计回声信号，包括：
基于获得N个子带估计回声信号，其中，为第i个子带估计回声信号，为第i个子带自适应滤波器的向量，i＝0,1,…,N-1。
所述将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，以至于消除所述参考信号经过第一延时传播的回声，包括：
基于获得N个子带残差信号，其中，为第i个子带估计回声信号，pi(n)为第i个主输入子带信号，ei(n)为第i个子带残差信号，i＝0,1,…,N-1。
每个子带自适应滤波器的系数更新为：基于 更新第i个子带自适应滤波器的系数，其中，μi(n)为第i个子带自适应滤波器的步长因子，Ai(n)为第i个子带自适应滤波器的趋势矩阵，Δ为常数，ei(n)为第i个子带残差信号，i＝0,1,…,N-1。
所述μi(n)为0和两者中的最大值，即：
其中，max(·)表示取较大值，abs(·)表示绝对值。Rpp,i(n)为第i个主输入子带信号平滑后的能量，为第i个子带估计回声信号平滑后的能量，Ree,i(n)为第i个子带残差信号平滑后的能量，ε为常数；
所述Ai(n)由以下表达式(1)得到：
其中，α为常数，取值范围[-1.0,1.0]；j＝0,1,…,Lw-1；|·|为取模运算，Ai(n)为对角方阵。
Rpp,i(n)、Ree,i(n)由以下表达式(2)获得：
其中，λ为常数，λ取值范围为[0.7,1.0]，||·||2表示欧几里得模2范数。
步骤S230：基于所述N个子带残差信号，获得全带残差信号；
步骤S240：基于所述全带残差信号以及获取到的所述参考信号经过延时修正后的延时参考信号，通过处理后获得输出信号。
本发明实施例提供的回声消除方法，所述方法包括：根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。提高回声消除过程中的计算效率。
第二实施例
请参阅图3，本发明实施例提供了一种回声消除装置300，所述装置300包括：
第一分析滤波器组单元310，用于根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号。
第二分析滤波器组单元320，用于根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声。
声学回声消除单元330，用于基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。
所述装置300还包括：
综合滤波器组单元340，用于基于所述N个子带残差信号，获得全带残差信号。
残留回声抑制单元350，用于基于所述全带残差信号以及获取到的所述参考信号经过延时修正后的延时参考信号，通过处理后获得输出信号。
作为一种实施方式，所述残留回声抑制单元350包括滤波器组延时补偿子单元351和残留回声抑制子单元352。
所述滤波器组延时补偿子单元351，用于获取到所述参考信号经过延时修正后的延时参考信号。
请参阅图4，参考信号x(n)经过第一分析滤波器组单元310，输出为N个子带参考信号即x0(n),x1(n),...,xN-1(n)，子带数为N，子带带宽为fs/(2N)，fs为采样频率。同理，主输入信号设为p(n)，p(n)＝s(n)+d(n)+b(n)，其中，d(n)为x(n)经过空间传播的回声信号，即d(n)定义为所述参考信号经过第一延时传播的回声，第一延时为K，s(n)为语音信号，b(n)为噪声信号。n为离散时间标度，n＝0,1,2……。主输入信号p(n)经过第二分析滤波器组单元320，输出为N个主输入子带信号p0(n),p1(n),...,pN-1(n)。
N个子带参考信号即x0(n),x1(n),...,xN-1(n)经过声学回声消除单元330，获得N个子带估计回声信号即p0(n),p1(n),...,pN-1(n)和相减得到N个子带残差信号即e0(n),e1(n),...,eN-1(n)。
N个子带残差信号e0(n),e1(n),...,eN-1(n)经过综合滤波器组单元340，得到全带残差信号即e(n)。x(n)经过滤波器组延时补偿子单元351，得到经过延时修正后的延时参考信号即xcomp(n)。xcomp(n)和e(n)输入到残留回声抑制子单元352中，通过处理后获得输出信号O(n)。
请参阅图5，将第一分析滤波器组单元310、第二分析滤波器组单元320、声学回声消除单元330、综合滤波器组单元340细节展开，H0,H1,...,HN-1为各个子带的分析带通滤波器，F0,F1,...,FN-1为各个子带的合成带通滤波器，↓D表示抽取D，↑D表示插值D；声学回声消除单元330中W0(n),W1(n),...,WN-1(n)为各个子带对应的自适应滤波器。子带自适应滤波器W0(n),W1(n),...,WN-1(n)为有限长度冲击响应(Finite Impulse Response，FIR)滤波器，每个子带自适应滤波器的长度为Lw＝L/N,L为全带自适应滤波器长度。相对于全带时域自适应滤波器计算效率提高N倍，N为子带数。在外界语音和噪声干扰下，利用自适应滤波器的变步长控制滤波器更新，大大降低自适应滤波器的不稳定性和提高自适应滤波器收敛速度。
在本实施例，分析滤波器组如第一分析滤波器组单元310、第二分析滤波器组单元320采用多项结构的调制滤波器组的形式。该多项结构的调制滤波器组使用一个长度为LP的FIR(Finite Impulse Response，有限长冲激响应)的仿酉原型滤波器h(n)，经过调制后，得到N个子带分析滤波器，该N个子带分析滤波器可以表示为Hk(z),k取值为0到N-1的整数。
请参阅图6，T为调制矩阵，对于DFT滤波器组，T表示离散傅里叶变换；对于余弦调制滤波器组，T表示离散余弦变换。L1为抽取矩阵，L2为一个排序矩阵。对于具有N个子带的分析滤波器组，其第K个子带的滤波器多项I型可以表示为公式(3)：
其中，Hk|m(z)为Hk(z)的多项I型第m个分量的Z域表示，D为抽取或插值数(可以使用抽取和插值一致的滤波器组)。
对于某一输入信号，其多项II型可以分解为以下公式(4)：
其中，Xm(z)为X(z)的多项II型第个分量。
将输入信号和分析滤波器组写成如下的向量形式：
根据(5)和(6)，分析滤波器组的输出为：
第k个子带的调制序列Lp同时为D和N的倍数，第k个子带滤波器调制后的滤波器系数为 和关系如下公式(7)：
其中，ID为D×D的单位阵，bi为原型滤波器系数，i＝0,1,...,Lp-1。
中的元素以N为周期，其表达式如下：
其中，为一个周期N内的离散调制向量，对于复指数调制滤波器组，为离散复指数调制向量；对于余弦调制滤波器组，为离散余弦调制向量。分析滤波器组单元的多项表达式(8)为：
H(z)＝T·L2·B·L1(z) (8)
表达式(8)中，为调制矩阵。
根据滤波器组单元完全重构条件同理可得到合成滤波器组单元的结构表达式(9)为：
表达式(9)中，为Lp×Lp反单位阵。
为完全重构后信号。为了便于工程实现，引入公式(10)：
公式(10)中，0D和ID分别为D×D零矩阵和D×D单位阵。
由(10)式，得到重构后的信号为
即V(z)最后个元素。
综上，调制分析滤波器组如本实施例中的第一分析滤波器组单元310、第二分析滤波器组单元320的快速实现由(8)式给出，调制合成滤波器组如本实施例中的综合滤波器组单元340的快速实现由(10)式给出。基于上述的调制滤波器组的快速实现方法，并结合子带自适应滤波器算法，大大提高计算效率。
以上各单元可以是由软件代码实现，此时，上述的各单元可存储于存储器110内。以上各单元同样可以由硬件例如集成电路芯片实现。
本发明实施例提供的回声消除装置300，其实现原理及产生的技术效果和前述方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的***来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
Claims (10)
1.一种回声消除方法，其特征在于，所述方法包括：
根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；
根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；
基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。
2.根据权利要求1所述的方法，其特征在于，所述N个子带参考信号对应有N个子带自适应滤波器，每个子带自适应滤波器的长度为L为全带自适应滤波器长度，所述基于所述N个子带参考信号，获得N个子带估计回声信号，包括：
基于获得N个子带估计回声信号，其中，为第i个子带估计回声信号，为第i个子带自适应滤波器的向量，i＝0,1,…,N-1。
3.根据权利要求2所述的方法，其特征在于，所述将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，以至于消除所述参考信号经过第一延时传播的回声，包括：
基于获得N个子带残差信号，其中，为第i个子带估计回声信号，pi(n)为第i个主输入子带信号，ei(n)为第i个子带残差信号，i＝0,1,…,N-1。
4.根据权利要求2所述的方法，其特征在于，每个子带自适应滤波器的系数更新为：基于更新第i个子带自适应滤波器的系数，其中，μi(n)为第i个子带自适应滤波器的步长因子，Ai(n)为第i个子带自适应滤波器的趋势矩阵，Δ为常数，ei(n)为第i个子带残差信号，i＝0,1,…,N-1。
5.根据权利要求4所述的方法，其特征在于，所述μi(n)为0和两者中的最大值，其中，Rpp,i(n)为第i个主输入子带信号平滑后的能量，为第i个子带估计回声信号平滑后的能量，Ree,i(n)为第i个子带残差信号平滑后的能量，ε为常数；
所述Ai(n)由以下表达式得到：
其中，α为常数，j＝0,1,…,Lw-1。
6.根据权利要求1所述的方法，其特征在于，所述根据获取到的参考信号，获得N个子带参考信号，包括：
对获取到的参考信号通过分析滤波器组，获得N个子带参考信号。
7.根据权利要求1所述的方法，其特征在于，所述方法还包括：
基于所述N个子带残差信号，获得全带残差信号；
基于所述全带残差信号以及获取到的所述参考信号经过延时修正后的延时参考信号，通过处理后获得输出信号。
8.一种回声消除装置，其特征在于，所述装置包括：
第一分析滤波器组单元，用于根据获取到的参考信号，获得所述参考信号对应的N个子带参考信号；
第二分析滤波器组单元，用于根据获取到的主输入信号，获得所述主输入信号对应的N个主输入子带信号，其中，所述主输入信号包括所述参考信号经过第一延时传播的回声；
声学回声消除单元，用于基于所述N个子带参考信号，获得N个子带估计回声信号，再将所述N个主输入子带信号减去所述N个子带估计回声信号，获得N个子带残差信号，从而消除所述参考信号经过第一延时传播的回声。
9.根据权利要求8所述的装置，其特征在于，所述装置还包括：
综合滤波器组单元，用于基于所述N个子带残差信号，获得全带残差信号；
残留回声抑制单元，用于基于所述全带残差信号以及获取到的所述参考信号经过延时修正后的延时参考信号，通过处理后获得输出信号。
10.一种回声消除***，其特征在于，所述***包括：电子设备、麦克风以及扬声器，所述电子设备包括如权利要求8-9任一项所述的装置，所述电子设备分别与所述麦克风、所述扬声器耦合。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710106188.1A CN106601266A (zh) | 2017-02-24 | 2017-02-24 | 回声消除方法、装置及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710106188.1A CN106601266A (zh) | 2017-02-24 | 2017-02-24 | 回声消除方法、装置及*** |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106601266A true CN106601266A (zh) | 2017-04-26 |
Family
ID=58588059
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710106188.1A Pending CN106601266A (zh) | 2017-02-24 | 2017-02-24 | 回声消除方法、装置及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106601266A (zh) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107799123A (zh) * | 2017-12-14 | 2018-03-13 | 南京地平线机器人技术有限公司 | 控制回声消除器的方法和具有回声消除功能的装置 |
| CN108831496A (zh) * | 2018-05-22 | 2018-11-16 | 出门问问信息科技有限公司 | 一种回声抑制方法及装置、存储介质、电子设备 |
| CN108986837A (zh) * | 2018-09-05 | 2018-12-11 | 科大讯飞股份有限公司 | 一种滤波器更新方法及装置 |
| CN109451398A (zh) * | 2018-11-16 | 2019-03-08 | 珠海市杰理科技股份有限公司 | 声反馈消除设备、声反馈消除方法、音频处理*** |
| CN110021289A (zh) * | 2019-03-28 | 2019-07-16 | 腾讯科技（深圳）有限公司 | 一种声音信号处理方法、装置及存储介质 |
| CN111213359A (zh) * | 2017-10-04 | 2020-05-29 | 主动音频有限公司 | 回声消除器和用于回声消除器的方法 |
| CN111988708A (zh) * | 2020-09-03 | 2020-11-24 | 昆腾微电子股份有限公司 | 一种基于单麦克风的啸叫抑制方法及装置 |
| CN112542176A (zh) * | 2020-11-04 | 2021-03-23 | 北京百度网讯科技有限公司 | 信号增强方法、装置及存储介质 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101923862A (zh) * | 2002-08-16 | 2010-12-22 | 艾玛复合信号公司 | 使用自适应滤波器处理子带信号的方法和*** |
| CN103730125A (zh) * | 2012-10-12 | 2014-04-16 | 华为技术有限公司 | 一种回声抵消方法和设备 |
| CN104395957A (zh) * | 2012-04-30 | 2015-03-04 | 创新科技有限公司 | 一种通用可重构回声消除*** |
| CN104427144A (zh) * | 2013-09-11 | 2015-03-18 | 联芯科技有限公司 | 一种线性回声消除方法及其装置 |
| CN106210368A (zh) * | 2016-06-20 | 2016-12-07 | 百度在线网络技术（北京）有限公司 | 消除多通道声回波的方法和装置 |
-
2017
- 2017-02-24 CN CN201710106188.1A patent/CN106601266A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101923862A (zh) * | 2002-08-16 | 2010-12-22 | 艾玛复合信号公司 | 使用自适应滤波器处理子带信号的方法和*** |
| CN104395957A (zh) * | 2012-04-30 | 2015-03-04 | 创新科技有限公司 | 一种通用可重构回声消除*** |
| CN103730125A (zh) * | 2012-10-12 | 2014-04-16 | 华为技术有限公司 | 一种回声抵消方法和设备 |
| CN104427144A (zh) * | 2013-09-11 | 2015-03-18 | 联芯科技有限公司 | 一种线性回声消除方法及其装置 |
| CN106210368A (zh) * | 2016-06-20 | 2016-12-07 | 百度在线网络技术（北京）有限公司 | 消除多通道声回波的方法和装置 |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111213359A (zh) * | 2017-10-04 | 2020-05-29 | 主动音频有限公司 | 回声消除器和用于回声消除器的方法 |
| CN111213359B (zh) * | 2017-10-04 | 2021-09-10 | 主动音频有限公司 | 回声消除器和用于回声消除器的方法 |
| CN107799123A (zh) * | 2017-12-14 | 2018-03-13 | 南京地平线机器人技术有限公司 | 控制回声消除器的方法和具有回声消除功能的装置 |
| CN108831496A (zh) * | 2018-05-22 | 2018-11-16 | 出门问问信息科技有限公司 | 一种回声抑制方法及装置、存储介质、电子设备 |
| CN108986837A (zh) * | 2018-09-05 | 2018-12-11 | 科大讯飞股份有限公司 | 一种滤波器更新方法及装置 |
| CN108986837B (zh) * | 2018-09-05 | 2021-08-17 | 科大讯飞股份有限公司 | 一种滤波器更新方法及装置 |
| CN109451398A (zh) * | 2018-11-16 | 2019-03-08 | 珠海市杰理科技股份有限公司 | 声反馈消除设备、声反馈消除方法、音频处理*** |
| CN110021289A (zh) * | 2019-03-28 | 2019-07-16 | 腾讯科技（深圳）有限公司 | 一种声音信号处理方法、装置及存储介质 |
| CN110021289B (zh) * | 2019-03-28 | 2021-08-31 | 腾讯科技（深圳）有限公司 | 一种声音信号处理方法、装置及存储介质 |
| CN111988708A (zh) * | 2020-09-03 | 2020-11-24 | 昆腾微电子股份有限公司 | 一种基于单麦克风的啸叫抑制方法及装置 |
| CN112542176A (zh) * | 2020-11-04 | 2021-03-23 | 北京百度网讯科技有限公司 | 信号增强方法、装置及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106601266A (zh) | 回声消除方法、装置及*** | |
| Chang et al. | Analysis of conjugate gradient algorithms for adaptive filtering | |
| CN107845389A (zh) | 一种基于多分辨率听觉倒谱系数和深度卷积神经网络的语音增强方法 | |
| CN101227537B (zh) | 一种宽带声学回声消除方法 | |
| CN106353649B (zh) | 一种基于提升小波变换的局部放电信号去噪方法 | |
| CN104427144B (zh) | 一种线性回声消除方法及其装置 | |
| Schertler | Selective block update of NLMS type algorithms | |
| CN102377454B (zh) | 一种回声消除的方法和装置 | |
| CN105825865B (zh) | 噪声环境下的回声消除方法及*** | |
| Xu et al. | An improved method for the design of FIR quadrature mirror-image filter banks | |
| Buchner et al. | A generalization of a class of blind source separation algorithms for convolutive mixtures | |
| CN112201273A (zh) | 一种噪声功率谱密度计算方法、***、设备及介质 | |
| Islam et al. | Performance comparison of modified LMS and RLS algorithms in de-noising of ECG signals | |
| Babatas et al. | Time and frequency based sparse bounded component analysis algorithms for convolutive mixtures | |
| Takahashi et al. | Steady-state mean-square performance analysis of a relaxed set-membership NLMS algorithm by the energy conservation argument | |
| CN107924685A (zh) | 信号处理装置和方法 | |
| Derkx et al. | New constraining method for partitioned block frequency-domain adaptive filters | |
| Comminiello et al. | Functional link expansions for nonlinear modeling of audio and speech signals | |
| CN106847299A (zh) | 延时的估计方法及装置 | |
| Cheng et al. | Wavelet packet transform applied to active noise control system for mixed noise in nonlinear environment | |
| Abutalebi et al. | Affine projection algorithm for oversampled subband adaptive filters | |
| CN1688104B (zh) | 数字信号处理方法及装置 | |
| CN114242103A (zh) | 语音降噪方法、装置、计算机设备和存储介质 | |
| CN106297813A (zh) | 分离的音频分析和处理 | |
| CN105957534A (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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170426 |
|
| WD01 | Invention patent application deemed withdrawn after publication |