CN113194385B - 基于步长控制的子带自适应反馈消除方法及系统 - Google Patents
基于步长控制的子带自适应反馈消除方法及系统 Download PDFInfo
- Publication number
- CN113194385B CN113194385B CN202110048056.4A CN202110048056A CN113194385B CN 113194385 B CN113194385 B CN 113194385B CN 202110048056 A CN202110048056 A CN 202110048056A CN 113194385 B CN113194385 B CN 113194385B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency domain
- module
- processing
- sub
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000008030 elimination Effects 0.000 title claims abstract description 12
- 238000003379 elimination reaction Methods 0.000 title claims abstract description 12
- 238000001914 filtration Methods 0.000 claims abstract description 49
- 238000006243 chemical reaction Methods 0.000 claims abstract description 19
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 11
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 11
- 238000002789 length control Methods 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 230000003044 adaptive effect Effects 0.000 description 6
- 230000003321 amplification Effects 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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 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
- G10L21/0224—Processing in the time domain
-
- 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
- G10L21/0232—Processing in the frequency domain
-
- 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/04—Circuits for transducers, loudspeakers or microphones for correcting frequency response
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
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)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
本发明公开了基于步长控制的子带自适应反馈消除方法及系统,本发明的方法包括以下步骤:获取麦克风拾音信号d(n)和扬声器输出延时信号x(n);分别对信号d(n)和x(n)依次进行预加重处理、子带滤波处理和FFT转换;进行频域滤波处理,得到频域差值信号E(f,k);对频域滤波器系数进行更新;频域差值信号E(f,k)经IFFT转换到时域后,依次经子带合成滤波处理和信号去加重处理得到时域信号e(n);时域信号e(n)经移频处理之后得到扬声器输出信号y(n)。本发明采用子带自适应算法,可用最少的滤波器长度达到较高的效果,同时减少大量的计算量,提高处理效率。
Description
技术领域
本发明属于音频处理技术领域,具体涉及基于步长控制的子带自适应反馈消除方法及系统。
背景技术
在专业音响工程中从剧场到会议室,从体育场馆到KTV,在这些音响工程中出现声音反馈产生的啸叫问题,是每种类型的扩声系统都会遇到的。
啸叫产生的原因是音响的声音折回到话筒,话筒输出的信号再送到扩声系统放大并再经音箱送出,而后又折回到话筒,如此循环所致。
传统抑制啸叫抑制的方法一是基于物理的,比如采用吸音材料,选用指向型的麦克风来降低环路增益,减少反馈,二是基于DSP的处理,如移频,移相,自动陷波,自适应反馈消除等本地语音和扩音器之间的信号存在很强的相关性,采用自适应方法收敛效率很慢,甚至会发散。
发明内容
为了解决现有自适应反馈消除收敛效果不佳,在扩声技术中传声增益不足,容易产生啸叫的问题,本发明提供了基于步长控制的子带自适应反馈消除方法。本发明通过改进步长控制,即在语音停顿的瞬间,此时没有本地语音,只存在之前语音的混响激励,此时自适应收敛效果好,增加步长;而在本地语音输入时,减小步长。本发明为了减小本地语音和扩音器之间信号的强相关性,在信号输出前增加移频处理。
本发明通过下述技术方案实现:
基于步长控制的子带自适应反馈消除方法,该方法包括以下步骤:
步骤S1,获取麦克风拾音信号d(n)和扬声器输出延时信号x(n);
步骤S2,对信号d(n)依次进行预加重处理、子带滤波处理和FFT转换之后得到频域D(f,k);对信号x(n)依次进行预加重处理、子带滤波处理和FFT转换之后得到频域信号X(f,k);其中,k=0,1,2,…,N/2-1,N表示FFT的点数;
步骤S3,采用频域滤波器对频域信号D(f,k)和X(f,k)进行频域滤波处理,得到频域差值信号E(f,k);
步骤S4,基于信号E(f,k)、D(f,k)和X(f,k)对频域滤波器系数进行更新;
步骤S5,频域差值信号E(f,k)经IFFT转换到时域后,依次经子带合成滤波处理和信号去加重处理得到时域信号e(n);
步骤S6,时域信号e(n)经移频处理之后得到扬声器输出信号y(n)。
优选的,本发明的步骤S4对频域滤波器系数进行更新的具体过程包括:
步骤S42,通过下式迭代更新参数λ(f,k):
λ(f,k)=beta(f,k)*λ(f,k-1)
式中,beta(f,k)表示迭代系数;λ(f,k)的初始值为1;
步骤S43,基于步骤S42得到的参数λ(f,k)计算得到最新步长α(f,k);
步骤S44,通过得到的最新步长α(f,k)对频域滤波器系数进行更新。
优选的,本发明的步骤S41通过下式计算频域信号的均值:
其中,alpha为平滑系数。
优选的,本发明的步骤S42中beta(f,k)通过下式计算得到:
优选的,本发明的步骤S43通过下式得到最新步长α(f,k):
优选的,本发明的S44通过下式进行滤波器系数更新:
G(f,k)=diag[g0(f,k),g1(f,k),...,gm-1(f,k)]
其中,gm(f,k)为第m个子带系数;L为子带长度;M为子带个数;β和ε均为常数;H(f,k)采用更新后的系数估算值;Hm(f,k)为第m个子带更新后的系数估算值;G(f,k)为全频带系数;为更新前的系数估算值,为更新后的系数估算值。
优选的,本发明的扬声器输出延时信号x(n)是输出信号y(n)经过一帧数据的延时得到。
优选的,本发明的步骤S2中的预加重传递函数为:
H(z)=1-az-1
式中,a为预加重系数。
另一方面,本发明还提出了基于步长控制的子带自适应反馈消除系统,该系统包括信号采集模块、信号预加重处理模块、信号子带滤波模块、FFT模块、频域滤波模块、更新模块、IFFT模块、子带合成滤波模块、信号去加重模块和移频模块;
所述信号采集模块用于获取麦克风拾音信号d(n)和扬声器输出延时信号 x(n);
所述信号d(n)依次通过所述信号预加重处理模块、所述信号子带滤波模块和 FFT模块进行预加重处理、子带滤波处理和FFT转换之后得到频域D(f,k);所述信号x(n)依次通过所述信号预加重处理模块、所述信号子带滤波模块和FFT 模块进行预加重处理、子带滤波处理和FFT转换之后得到频域信号X(f,k);其中, k=0,1,2,…,N/2-1,N表示FFT的点数;
所述频域滤波模块采用频域滤波器对频域信号D(f,k)和X(f,k)进行频域滤波处理,得到频域差值信号E(f,k);
所述更新模块基于信号E(f,k)、D(f,k)和X(f,k)对所述频域滤波器系数进行更新;
所述IFFT转换模块将频域差值信号E(f,k)转换到时域信号,所述IFFT转换模块输出的时域信号依次经所述子带合成滤波模块和所述信号去加重模块处理得到时域信号e(n);
所述移频模块用于对时域信号e(n)进行移频处理之后得到输出信号y(n),所述输出信号y(n)经功放通过扬声器播放输出。
本发明具有如下的优点和有益效果:
1、本发明采用子带自适应算法,可用最少的滤波器长度达到较高的效果,同时减少大量的计算量,提高处理效率。
2、本发明采用了步长控制,本地语音过程中减小步长,而在本地语音停止的瞬间,在短暂的混响周期内会加大步长,让自适应达到一个很好的收敛效果。
3、本发明采用了移频后输出,降低了本地语音和扩音器之间信号的强相关性,提高了处理性能。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:
图1为本发明的方法流程示意图。
图2为本发明的系统原理框图。
图3为本发明的具体应用实施例。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。
实施例1
相较于现有的自适应反馈消除技术,本实施例提出了一种基于步长控制的子带自适应反馈消除方法,通过改进补偿控制,实现在本地语音过程中会减小步长,而在本地语音停止的瞬间,在短暂的混响周期内会加大步长,提高自适应的收敛性能。
具体如图1所示,本实施例的方法包括以下步骤:
步骤一,获取麦克风拾音信号d(n)和扬声器输出延时信号x(n)。
本实施例中通过将麦克风采集到的信号转换为数字音频信号,即为麦克风拾音信号d(n),d(n)=s(n)+r(n)。其中,s(n)为本地语言,r(n)为反馈信号。
本实施例中将扬声器输出信号y(n)经过一帧数据的延时得到扬声器输出延时信号x(n)。
步骤二,对信号d(n)依次进行预加重处理、子带滤波处理和FFT转换之后得到频域D(f,k);对信号x(n)依次进行预加重处理、子带滤波处理和FFT转换之后得到频域信号X(f,k);其中,k=0,1,2,…,N/2-1,N表示FFT的点数。
本实施例的预加重传递函数为:
H(z)=1-az-1
式中,a为预加重系数。
步骤三,采用频域滤波器对频域信号D(f,k)和X(f,k)进行频域滤波处理,得到频域差值信号E(f,k):
步骤四,基于信号E(f,k)、D(f,k)和X(f,k)对频域滤波器系数进行更新。本实施例的具体更新过程如下:
其中,alpha为平滑系数。
步骤4.2,通过下式迭代更新参数λ(f,k):
λ(f,k)=beta(f,k)*λ(f,k-1)
式中,beta(f,k)表示迭代系数;λ(f,k)的初始值为1;
本实施例中,beta(f,k)通过下式计算得到:
步骤4.3,通过下式计算得到最新步长α(f,k):
步骤4.4,通过下式对频域滤波器系数进行更新:
G(f,k)=diag[g0(f,k),g1(f,k),...,gm-1(f,k)]
其中,gm(f,k)为第m个子带系数;L为子带长度,本实施例中L取值为 256;M为子带个数,本实施例中M取值为4;β取值为0.5;ε取值为1e-8;H(f,k)采用更新后的系数估算值Hm(f,k)为第m个子带更新后的系数估算值 (通过将更新后的系数估算值分到M个子带获得);G(f,k)为全频带系数(由m 个子带系数合成);为更新前的系数估算值,为更新后的系数估算值。
步骤五,频域差值信号E(f,k)经IFFT转换到时域后,依次经子带合成滤波处理和信号去加重处理得到时域信号e(n);
步骤六,时域信号e(n)经移频处理之后得到扬声器输出信号y(n)。
实施例2
本实施例提出了一种基于步长控制的子带自适应反馈消除系统,具体如图2 所示,本实施例的系统包括信号采集模块、信号预加重处理模块、信号子带滤波模块、FFT模块、频域滤波模块、更新模块、IFFT模块、子带合成滤波模块、信号去加重模块和移频模块;
信号采集模块用于获取麦克风拾音信号d(n)和扬声器输出延时信号x(n);具体处理过程与上述实施例1相同。
信号d(n)依次通过信号预加重处理模块、信号子带滤波模块和FFT模块进行预加重处理、子带滤波处理和FFT转换之后得到频域D(f,k);信号x(n)依次通过信号预加重处理模块、信号子带滤波模块和FFT模块进行预加重处理、子带滤波处理和FFT转换之后得到频域信号X(f,k);其中,k=0,1,2,…,N/2-1,N 表示FFT的点数;具体处理过程与上述实施例1相同。
频域滤波模块采用频域滤波器对频域信号D(f,k)和X(f,k)进行频域滤波处理,得到频域差值信号E(f,k);具体处理过程与上述实施例1相同。
更新模块基于信号E(f,k)、D(f,k)和X(f,k)对频域滤波器系数进行更新;具体更新过程与上述实施例1相同。
IFFT转换模块将频域差值信号E(f,k)转换到时域信号,IFFT转换模块输出的时域信号依次经子带合成滤波模块和信号去加重模块处理得到时域信号e(n);具体处理过程与上述实施例1相同。
移频模块用于对时域信号e(n)进行移频处理之后得到输出信号y(n),输出信号y(n)经功放通过扬声器播放输出。
实施例3
本实施例为实施例1和实施例2的具体应用,其硬件系统由一个麦克风,一个扬声器,一块DSP音频处理板,根据附图3说明,系统在一个房间内,扬声器通过房间反射到麦克风,形成环路。
麦克风拾取本地语音,经AD电路转换成数字信号,得到d(n),扬声器输出信号y(n)经过延时得到x(n)。
音频系统48K采样,然后软件降采样到16K处理算法,一帧数据缓冲点数为64点,FFT变换点数为256点,拖尾处理以256ms处理,系统分成4个子带计算。
具体实施过程如下:
一,将x(n)和d(n)分别进行预加重处理,预加重传函为:H(z)=1-az-1;本实施例中a取值为0.9。
二,预加重处理后,进行子带滤波处理和FFT转换分别得到X(f,k),D(f,k), k=0,1,2,...N/2-1,其中N为FFT的点数;
三,对X(f,k),D(f,k)进行频域滤波处理得到E(f,k),并计算三个信号的均值;
其中,alpha的值取0.98。
四,对频域滤波器系数进行更新。
五,频域滤波后得到的信号E(f,k),经IFFT转换到时域,经子带合成滤波后得到时域信号e(n),移频处理后得到y(n)发送到扬声器。
重复一到五,实现反馈消除,提升扩声系统增益。
系统先不启用本发明反馈消除,增加增益,直到系统开始出现啸叫,记录此时的增益为G1,然后开避本发明,增加增益,到开如出现啸叫,记录此时的增益为G2,系统提升增益为G=G1-G2;启用本发明系统增益可提升10-15dB, 传统的移频法,自动陷波法提升增益在4-6dB。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、 CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和 /或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/ 或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.基于步长控制的子带自适应反馈消除方法,其特征在于,该方法包括以下步骤:
步骤S1,获取麦克风拾音信号d(n)和扬声器输出延时信号x(n);
步骤S2,对信号d(n)依次进行预加重处理、子带滤波处理和FFT转换之后得到频域D(f,k);对信号x(n)依次进行预加重处理、子带滤波处理和FFT转换之后得到频域信号X(f,k);其中,k=0,1,2,…,N/2-1,N表示FFT的点数;
步骤S3,采用频域滤波器对频域信号D(f,k)和X(f,k)进行频域滤波处理,得到频域差值信号E(f,k);
步骤S4,基于信号E(f,k)、D(f,k)和X(f,k)对频域滤波器系数进行更新;
步骤S5,频域差值信号E(f,k)经IFFT转换到时域后,依次经子带合成滤波处理和信号去加重处理得到时域信号e(n);
步骤S6,时域信号e(n)经移频处理之后得到扬声器输出信号y(n);
所述步骤S4对频域滤波器系数进行更新的具体过程包括:
步骤S42,通过下式迭代更新参数λ(f,k):
λ(f,k)=beta(f,k)*λ(f,k-1)
式中,beta(f,k)表示迭代系数;λ(f,k)的初始值为1;
所述步骤S42中beta(f,k)通过下式计算得到:
步骤S43,基于步骤S42得到的参数λ(f,j)计算得到最新步长α(f,k);
步骤S44,通过得到的最新步长α(f,k)对频域滤波器系数进行更新。
5.根据权利要求1-4任一项所述的基于步长控制的子带自适应反馈消除方法,其特征在于,所述扬声器输出延时信号x(n)是输出信号y(n)经过一帧数据的延时得到。
6.根据权利要求1-4任一项所述的基于步长控制的子带自适应反馈消除方法,其特征在于,所述步骤S2中的预加重传递函数为:
H(z)=1-az-1
式中,a为预加重系数。
7.基于步长控制的子带自适应反馈消除系统,其特征在于,该系统包括信号采集模块、信号预加重处理模块、信号子带滤波模块、FFT模块、频域滤波模块、更新模块、IFFT模块、子带合成滤波模块、信号去加重模块和移频模块;
所述信号采集模块用于获取麦克风拾音信号d(n)和扬声器输出延时信号x(n);
所述信号d(n)依次通过所述信号预加重处理模块、所述信号子带滤波模块和FFT模块进行预加重处理、子带滤波处理和FFT转换之后得到频域D(f,k);所述信号x(n)依次通过所述信号预加重处理模块、所述信号子带滤波模块和FFT模块进行预加重处理、子带滤波处理和FFT转换之后得到频域信号X(f,k);其中,k=0,1,2,…,N/2-1,N表示FFT的点数;
所述频域滤波模块采用频域滤波器对频域信号D(f,k)和X(f,k)进行频域滤波处理,得到频域差值信号E(f,k);
所述更新模块基于信号E(f,k)、D(f,k)和X(f,k)对所述频域滤波器系数进行更新;对频域滤波器系数进行更新的具体过程包括:
通过下式迭代更新参数λ(f,k):
λ(f,k)=beta(f,k)*λ(f,k-1)
式中,beta(f,k)表示迭代系数;λ(f,k)的初始值为1;
所述beta(f,k)通过下式计算得到:
基于得到的参数λ(f,k)计算得到最新步长α(f,k);
通过得到的最新步长α(f,k)对频域滤波器系数进行更新;
所述IFFT模块将频域差值信号E(f,k)转换到时域信号,所述IFFT模块输出的时域信号依次经所述子带合成滤波模块和所述信号去加重模块处理得到时域信号e(n);
所述移频模块用于对时域信号e(n)进行移频处理之后得到输出信号y(n),所述输出信号y(n)经功放通过扬声器播放输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110048056.4A CN113194385B (zh) | 2021-01-14 | 2021-01-14 | 基于步长控制的子带自适应反馈消除方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110048056.4A CN113194385B (zh) | 2021-01-14 | 2021-01-14 | 基于步长控制的子带自适应反馈消除方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113194385A CN113194385A (zh) | 2021-07-30 |
CN113194385B true CN113194385B (zh) | 2023-03-10 |
Family
ID=76972603
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110048056.4A Active CN113194385B (zh) | 2021-01-14 | 2021-01-14 | 基于步长控制的子带自适应反馈消除方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113194385B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113938112B (zh) * | 2021-08-12 | 2023-04-07 | 襄阳达安汽车检测中心有限公司 | 无相移滤波数据的生成方法和无相移滤波器 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0000688D0 (en) * | 1999-01-12 | 2000-03-08 | Motorola Inc | Adapative step size for filter and method therefor |
CN101763858A (zh) * | 2009-10-19 | 2010-06-30 | 瑞声声学科技(深圳)有限公司 | 双麦克风信号处理方法 |
JP2014096027A (ja) * | 2012-11-09 | 2014-05-22 | Nippon Telegr & Teleph Corp <Ntt> | パラメータ推定装置、エコー消去装置、パラメータ推定方法、及びプログラム |
EP2905778A1 (en) * | 2012-10-12 | 2015-08-12 | Huawei Technologies Co., Ltd. | Echo cancellation method and device |
CN106454642A (zh) * | 2016-09-23 | 2017-02-22 | 佛山科学技术学院 | 自适应子带音频反馈抑制方法 |
CN108200522A (zh) * | 2017-11-24 | 2018-06-22 | 华侨大学 | 一种变正则化比例归一化子带自适应滤波方法 |
CN109727604A (zh) * | 2018-12-14 | 2019-05-07 | 上海蔚来汽车有限公司 | 用于语音识别前端的频域回声消除方法及计算机储存介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005096274A1 (fr) * | 2004-04-01 | 2005-10-13 | Beijing Media Works Co., Ltd | Dispositif et procede de codage/decodage audio ameliores |
EP2530840B1 (en) * | 2011-05-30 | 2014-09-03 | Harman Becker Automotive Systems GmbH | Efficient sub-band adaptive FIR-filtering |
US10720138B2 (en) * | 2017-04-24 | 2020-07-21 | Cirrus Logic, Inc. | SDR-based adaptive noise cancellation (ANC) system |
CN109754813B (zh) * | 2019-03-26 | 2020-08-25 | 南京时保联信息科技有限公司 | 基于快速收敛特性的变步长回声消除方法 |
CN111988708A (zh) * | 2020-09-03 | 2020-11-24 | 昆腾微电子股份有限公司 | 一种基于单麦克风的啸叫抑制方法及装置 |
-
2021
- 2021-01-14 CN CN202110048056.4A patent/CN113194385B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0000688D0 (en) * | 1999-01-12 | 2000-03-08 | Motorola Inc | Adapative step size for filter and method therefor |
CN101763858A (zh) * | 2009-10-19 | 2010-06-30 | 瑞声声学科技(深圳)有限公司 | 双麦克风信号处理方法 |
EP2905778A1 (en) * | 2012-10-12 | 2015-08-12 | Huawei Technologies Co., Ltd. | Echo cancellation method and device |
JP2014096027A (ja) * | 2012-11-09 | 2014-05-22 | Nippon Telegr & Teleph Corp <Ntt> | パラメータ推定装置、エコー消去装置、パラメータ推定方法、及びプログラム |
CN106454642A (zh) * | 2016-09-23 | 2017-02-22 | 佛山科学技术学院 | 自适应子带音频反馈抑制方法 |
CN108200522A (zh) * | 2017-11-24 | 2018-06-22 | 华侨大学 | 一种变正则化比例归一化子带自适应滤波方法 |
CN109727604A (zh) * | 2018-12-14 | 2019-05-07 | 上海蔚来汽车有限公司 | 用于语音识别前端的频域回声消除方法及计算机储存介质 |
Non-Patent Citations (2)
Title |
---|
一种新的基于L0的变步长IPNLMS算法;刘遵雄等;《计算机仿真》;20121115(第11期);全文 * |
归一化子带自适应滤波器步长控制;倪锦根等;《电子与信息学报》;20091115(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113194385A (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109727604B (zh) | 用于语音识别前端的频域回声消除方法及计算机储存介质 | |
JP4161628B2 (ja) | エコー抑圧方法及び装置 | |
US7092529B2 (en) | Adaptive control system for noise cancellation | |
CN108172231B (zh) | 一种基于卡尔曼滤波的去混响方法及系统 | |
JP4681163B2 (ja) | ハウリング検出抑圧装置、これを備えた音響装置、及び、ハウリング検出抑圧方法 | |
JP4221222B2 (ja) | 非対称マルチチャンネルフィルタ | |
KR101671389B1 (ko) | 가변 대역 폭 적응 노치 필터, 및 가변 대역 폭 적응 노치 필터를 이용하여 하울링을 제거하는 방법 및 장치 | |
US8538052B2 (en) | Generation of probe noise in a feedback cancellation system | |
JPH06500217A (ja) | 静的干渉キヤンセラ | |
WO2014056328A1 (zh) | 一种回声抵消方法和设备 | |
CN102739886A (zh) | 基于回声频谱估计和语音存在概率的立体声回声抵消方法 | |
JP2012501152A (ja) | 前白色化を伴うlmsアルゴリズムによって適応させられる適応フィルタの更新済みフィルタ係数を決定する方法 | |
CN106454642A (zh) | 自适应子带音频反馈抑制方法 | |
CN101635873B (zh) | 用于自适应白化的自适应长时预测滤波器 | |
JP2004349806A (ja) | 多チャネル音響エコー消去方法、その装置、そのプログラム及びその記録媒体 | |
CN115175063A (zh) | 啸叫抑制方法、装置、音响及扩音系统 | |
CN113194385B (zh) | 基于步长控制的子带自适应反馈消除方法及系统 | |
CN103929704B (zh) | 一种基于变换域的自适应声反馈消除的方法和系统 | |
JP3756828B2 (ja) | 反響消去方法、この方法を実施する装置、プログラムおよびその記録媒体 | |
CN104754465B (zh) | 一种自适应信号增强方法和系统 | |
CN113676157B (zh) | 一种抵消直流偏置的自适应滤波算法 | |
JP2002140100A (ja) | 騒音抑圧装置 | |
JP2002223182A (ja) | 反響消去方法、その装置、そのプログラム及びその記録媒体 | |
JP2004349796A (ja) | 音響エコー消去方法、その装置、プログラム及びその記録媒体 | |
EP1305975A4 (en) | ADAPTIVE MICROPHONE ARRAY SYSTEM WITH CONSERVATION OF BINAURAL NOTES |
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 |