CN112260662A - Adaptive filtering method, computer equipment and device - Google Patents
Adaptive filtering method, computer equipment and device Download PDFInfo
- Publication number
- CN112260662A CN112260662A CN202010968471.7A CN202010968471A CN112260662A CN 112260662 A CN112260662 A CN 112260662A CN 202010968471 A CN202010968471 A CN 202010968471A CN 112260662 A CN112260662 A CN 112260662A
- Authority
- CN
- China
- Prior art keywords
- current frame
- reference signal
- length
- filter coefficient
- 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.)
- Granted
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0043—Adaptive algorithms
Landscapes
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
Description
技术领域technical field
本申请涉及滤波领域,特别是涉及一种自适应滤波方法、计算机设备和装置。The present application relates to the field of filtering, and in particular, to an adaptive filtering method, computer equipment and apparatus.
背景技术Background technique
众所周知,音频通信已经成为人们工作和生活中常用的信息交流方法。当人们使用手机开启免提电话或视频会议终端进行视频会议时,由于声音外放,因此会导致扬声器所播放的声音再次被麦克风采集到,从而导致回声。回声信号会严重影响通信质量并且降低语音识别的准确性,因此有必要采取有效措施来抑制回声,提高通信质量。As we all know, audio communication has become a common method of information exchange in people's work and life. When people use a mobile phone to turn on a hands-free phone or a video conferencing terminal for a video conference, the sound played by the speaker will be picked up by the microphone again due to the external sound, which will cause echoes. The echo signal will seriously affect the communication quality and reduce the accuracy of speech recognition, so it is necessary to take effective measures to suppress the echo and improve the communication quality.
通常,采用自适应滤波算法(回声消除算法)对所采集到的信号进行滤波,以消除回声信号。本申请的发明人在长期的研发过程中,发现自适应滤波算法中,尤其是在滤波系数更新的算法中存在计算量不均匀的问题,例如,计算量陡增。当运算系统负载过大时,容易导致自适应滤波算法失效乃至系统崩溃。Usually, an adaptive filtering algorithm (echo canceling algorithm) is used to filter the collected signal to cancel the echo signal. During the long-term research and development process, the inventor of the present application found that the adaptive filtering algorithm, especially the algorithm for updating the filter coefficients, has the problem of uneven calculation amount, for example, the calculation amount increases sharply. When the load of the computing system is too large, it is easy to cause the failure of the adaptive filtering algorithm and even the collapse of the system.
发明内容SUMMARY OF THE INVENTION
本申请主要解决的技术问题是提供一种自适应滤波方法、计算机设备和装置,能够避免滤波过程中计算量不均匀,提高滤波效果。The main technical problem to be solved by the present application is to provide an adaptive filtering method, computer equipment and device, which can avoid the uneven calculation amount in the filtering process and improve the filtering effect.
为解决上述技术问题,本申请采用的一个技术方案是:提供一种自适应滤波方法,包括:获取当前帧输入信号、当前帧参考信号以及当前帧滤波系数;将当前帧参考信号与历史帧参考信号进行组合,以获得当前帧参考信号集;基于当前帧参考信号和当前帧滤波系数对当前帧输入信号进行自适应滤波,基于当前帧参考信号集与当前帧输入信号对当前帧滤波系数按照预设的划分长度进行更新,以获得后续帧的滤波系数;其中,当前帧参考信号集和历史帧参考信号集中包含数据相同的至少两个部分,历史帧和当前帧之间的间隔帧数与划分长度的乘积等于至少两个部分的差值与当前帧输入信号长度的乘积。In order to solve the above-mentioned technical problems, a technical solution adopted in the present application is to provide an adaptive filtering method, including: obtaining the current frame input signal, the current frame reference signal and the current frame filter coefficient; The signal is combined to obtain the current frame reference signal set; the current frame input signal is adaptively filtered based on the current frame reference signal and the current frame filter coefficient, and the current frame filter coefficient is determined according to the preset based on the current frame reference signal set and the current frame input signal. The set division length is updated to obtain the filter coefficient of the subsequent frame; wherein, the current frame reference signal set and the historical frame reference signal set contain at least two parts with the same data, and the interval frame number between the historical frame and the current frame is the same as the division. The product of the lengths is equal to the product of the difference between the at least two parts and the length of the input signal of the current frame.
其中,将当前帧参考信号与历史帧参考信号进行组合,以获得当前帧参考信号集包括:将当前帧参考信号与历史帧参考信号按照预设更新信号块长度进行组合,以获得多个当前帧参考信号分块;当前帧参考信号集包括多个当前帧参考信号分块;其中,更新信号块长度由划分长度与当前帧输入信号长度所确定。Wherein, combining the current frame reference signal and the historical frame reference signal to obtain the current frame reference signal set includes: combining the current frame reference signal and the historical frame reference signal according to the preset update signal block length to obtain a plurality of current frames The reference signal is divided into blocks; the current frame reference signal set includes a plurality of current frame reference signal blocks; wherein, the length of the update signal block is determined by the division length and the length of the current frame input signal.
其中,更新信号块长度不小于划分长度与当前帧输入信号长度之和减1的值。Wherein, the length of the update signal block is not less than the sum of the division length and the length of the input signal of the current frame minus 1.
其中,更新信号块长度等于2的幂数。Among them, the length of the update signal block is equal to the power of 2.
其中,将当前帧参考信号与历史帧参考信号按照预设更新信号块长度进行组合,以获得多个当前帧参考信号分块的步骤之后包括:将当前帧滤波系数按照划分长度进行划分,以获得多个当前帧滤波系数分块,其中当前帧滤波系数分块与当前帧参考信号分块一一对应;对当前帧滤波系数分块以及对应的当前帧参考信号分块进行卷积,以获得卷积信号;基于当前帧卷积信号与当前帧输入信号的差,更新当前帧滤波系数,以获得当前帧之后的后续帧的滤波系数。The step of combining the current frame reference signal and the historical frame reference signal according to the preset update signal block length to obtain a plurality of current frame reference signal blocks includes: dividing the current frame filter coefficients according to the division lengths to obtain A plurality of current frame filter coefficient blocks, wherein the current frame filter coefficient blocks are in one-to-one correspondence with the current frame reference signal blocks; the current frame filter coefficient blocks and the corresponding current frame reference signal blocks are convolved to obtain a volume product signal; based on the difference between the convolution signal of the current frame and the input signal of the current frame, update the filter coefficient of the current frame to obtain the filter coefficient of the subsequent frame after the current frame.
其中,当前帧和历史帧的参考信号集包含数据相同的至少两个参考信号分块,分别为第一参考信号分块和第二参考信号分块;历史帧和当前帧之间的间隔帧数与划分长度的乘积等于第二参考信号分块和第一参考信号分块的差值与当前帧输入信号长度的乘积。Wherein, the reference signal sets of the current frame and the historical frame include at least two reference signal sub-blocks with the same data, which are the first reference signal sub-block and the second reference signal sub-block respectively; The product with the division length is equal to the product of the difference between the second reference signal block and the first reference signal block and the length of the input signal of the current frame.
其中,划分长度为每帧输入信号长度的二分之一;或划分长度为每帧输入信号长度的二倍。Wherein, the division length is half the length of the input signal of each frame; or the division length is twice the length of the input signal of each frame.
其中,将当前帧参考信号与历史帧参考信号按照预设更新信号块长度进行组合,以获得多个当前帧参考信号分块的步骤包括:当当前帧输入信号长度为160,划分长度为80时,将当前帧输入参考信号和历史帧参考信号以256长度进行组合,以获得多个当前帧参考信号分块;或当当前帧输入信号长度为160,划分长度为320时,将当前帧输入参考信号和历史帧参考信号以512长度进行组合,以获得多个当前帧参考信号分块。The step of combining the current frame reference signal and the historical frame reference signal according to the preset update signal block length to obtain multiple current frame reference signal blocks includes: when the current frame input signal length is 160 and the division length is 80 , combine the current frame input reference signal and the historical frame reference signal with a length of 256 to obtain multiple current frame reference signal blocks; or when the current frame input signal length is 160 and the division length is 320, the current frame input reference signal The signal and the historical frame reference signal are combined with a length of 512 to obtain multiple current frame reference signal blocks.
为解决上述技术问题,本申请采用的另一个技术方案是:提供一种计算机设备,包括处理器,该处理器用于执行指令以实现上述的自适应滤波方法。In order to solve the above technical problem, another technical solution adopted in the present application is to provide a computer device including a processor for executing instructions to implement the above adaptive filtering method.
为解决上述技术问题,本申请采用的另一个技术方案是:提供一种具有存储功能的装置,装置存储有程序,程序被执行时能够实现上述的自适应滤波方法。In order to solve the above technical problem, another technical solution adopted in the present application is to provide a device with a storage function, the device stores a program, and the above-mentioned adaptive filtering method can be implemented when the program is executed.
本申请的有益效果是:区别于现有技术的情况,本申请提供了一种自适应滤波方法,通过设定滤波系数的更新中的划分长度,并使得划分长度与当前帧输入信号长度满足约束关系。使得在每次参考系数的更新过程中可以将当前帧输入信号全部进行计算,从而避免了部分信号数据滞留。进而,每次滤波系数更新的计算量相等,避免了计算量不均匀的问题,还可以有效的提高更新算法的更新效果。The beneficial effects of the present application are: different from the situation in the prior art, the present application provides an adaptive filtering method, by setting the division length in the update of the filter coefficient, and making the division length and the current frame input signal length satisfy constraints relation. This makes it possible to calculate all the input signals of the current frame in the updating process of each reference coefficient, thereby avoiding the retention of part of the signal data. Furthermore, the calculation amount of each filter coefficient update is equal, which avoids the problem of uneven calculation amount, and can effectively improve the update effect of the update algorithm.
附图说明Description of drawings
图1是根据本申请一实施例的自适应滤波方法的流程示意图;1 is a schematic flowchart of an adaptive filtering method according to an embodiment of the present application;
图2是根据本申请一实施例的滤波系数更新方法的流程示意图;2 is a schematic flowchart of a filter coefficient update method according to an embodiment of the present application;
图3是根据本申请一实施例的划分长度确定方法的流程示意图;3 is a schematic flowchart of a method for determining a division length according to an embodiment of the present application;
图4是根据本申请一实施例的自适应滤波器的结构示意图;4 is a schematic structural diagram of an adaptive filter according to an embodiment of the present application;
图5是根据本申请一实施例的计算机设备的结构示意图;5 is a schematic structural diagram of a computer device according to an embodiment of the present application;
图6是根据本申请一实施例的具有存储功能的装置的结构示意图。FIG. 6 is a schematic structural diagram of a device with a storage function according to an embodiment of the present application.
具体实施方式Detailed ways
为使本申请的目的、技术方案及效果更加清楚、明确,以下参照附图并举实施例对本申请进一步详细说明。In order to make the objectives, technical solutions and effects of the present application clearer and clearer, the present application will be further described in detail below with reference to the accompanying drawings and examples.
本申请的实施例可以用于音频设备中音频的回声消除处理,其中音频设备可以包括但不限于移动通讯设备、电脑、话筒等。应当理解的是,本申请的系统及方法的应用场景仅仅是本申请的一些示例或实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图将本申请应用于其他类似情景,例如视频的自适应滤波或其他数据的自适应滤波。The embodiments of the present application can be used for echo cancellation processing of audio in an audio device, where the audio device may include, but is not limited to, a mobile communication device, a computer, a microphone, and the like. It should be understood that the application scenarios of the system and method of the present application are only some examples or embodiments of the present application. For those of ordinary skill in the art, without creative work, they can also use these drawings according to the drawings. Apply this application to other similar scenarios, such as adaptive filtering of video or adaptive filtering of other data.
自适应滤波算法可以根据环境的改变使用自适应算法来改变滤波系数或结构。其中,滤波系数可以自动连续的适应与给定信号,以获得期望响应。自适应滤波方法中滤波系数的更新方法可以包括:最小均方差法(Least Mean Square,即LMS)、递推最小二乘法(Recursive Least Square,即RLS)、归一化最小均方差法(Normalized Least MeanSquares,NLMS)或仿射投影算法(Affine Projection Algorithm,APA)。由于在更新过程中,需要对滤波系数进行划分成多个滤波系数分块后再进行计算。Adaptive Filtering Algorithms Adaptive algorithms can be used to change filter coefficients or structures according to changes in the environment. Among them, the filter coefficients can be automatically and continuously adapted to the given signal to obtain the desired response. The updating method of the filter coefficients in the adaptive filtering method may include: Least Mean Square (LMS), Recursive Least Square (RLS), Normalized Least Square (Normalized Least Square) MeanSquares, NLMS) or Affine Projection Algorithm (Affine Projection Algorithm, APA). In the update process, the filter coefficient needs to be divided into a plurality of filter coefficient blocks before calculation.
例如,每帧输入信号的长度可以是160,滤波系数的划分长度为128。那么在进行更新过程,每次分块计算时,输入信号仅有128个数据进行计算,还会有32个数据点滞留。当滞留点数目达到128时,那一时刻会执行两次计算,即当前帧输入信号计算以及滞留点的计算。此时,运算量会翻倍。若自适应滤波方法的算法较为复杂,预算量陡增时会导致无法在规定的时间内完成所有数据的计算,从而导致算法失效乃至系统崩溃。For example, the length of the input signal of each frame may be 160, and the division length of the filter coefficients may be 128. Then in the update process, each time the block is calculated, the input signal has only 128 data for calculation, and 32 data points remain. When the number of stuck points reaches 128, two calculations will be performed at that moment, namely the calculation of the input signal of the current frame and the calculation of the stuck points. At this point, the amount of computation will be doubled. If the algorithm of the adaptive filtering method is relatively complex, when the budget increases sharply, the calculation of all data cannot be completed within the specified time, resulting in the failure of the algorithm and even the collapse of the system.
为解决上述问题,本申请公开了一种自适应滤波方法。参阅图1,图1是根据本申请一实施例的自适应滤波方法的流程示意图。需注意的是,若有实质上相同的结果,本实施例并不以图1所示的流程顺序为限。如图1所示,该方法包括:To solve the above problems, the present application discloses an adaptive filtering method. Referring to FIG. 1 , FIG. 1 is a schematic flowchart of an adaptive filtering method according to an embodiment of the present application. It should be noted that, if there is substantially the same result, the present embodiment is not limited to the sequence of the processes shown in FIG. 1 . As shown in Figure 1, the method includes:
步骤S110:获取当前帧输入信号、当前帧参考信号以及当前帧滤波系数。Step S110: Obtain the current frame input signal, the current frame reference signal and the current frame filter coefficient.
在一实施例中,当前帧输入信号和当前帧参考信号可以是指当前时间段内获取得到的输入信号和参考信号。In an embodiment, the current frame input signal and the current frame reference signal may refer to the input signal and the reference signal acquired in the current time period.
输入信号可以是通过音频采集设备采集并传输。其中音频传输一般是以10ms或是20ms为一帧进行传输。例如,在输入的音频信号的采样率为16k时,一帧输入信号数据块的离散点数为160点或320点,该情况下可以将每帧输入信号长度记为160。The input signal can be collected and transmitted through an audio capture device. Among them, the audio transmission is generally carried out in a frame of 10ms or 20ms. For example, when the sampling rate of the input audio signal is 16k, the discrete points of a frame of input signal data block are 160 points or 320 points, in this case, the length of each frame of input signal can be recorded as 160.
参考信号可以是从网络端获取的。在不同的场景下进行自适应滤波所采用的参考信号会有所不同。例如,在回声消除处理中,参考信号可以是扬声器所播放的前一帧的声音信号。The reference signal may be obtained from the network side. The reference signal used for adaptive filtering in different scenarios will be different. For example, in the echo cancellation process, the reference signal may be the sound signal of the previous frame played by the speaker.
自适应滤波方法中的滤波系数会在滤波过程中不断更新。在一实施例中,处理当前帧输入信号时所采用的当前帧滤波系数可以为根据前一帧输入信号和前一帧参考信号更新后的滤波系数。滤波系数可以为一个包含N个数据的一维矩阵。The filter coefficients in the adaptive filtering method are continuously updated during the filtering process. In one embodiment, the current frame filter coefficient used when processing the current frame input signal may be the filter coefficient updated according to the previous frame input signal and the previous frame reference signal. The filter coefficient can be a one-dimensional matrix containing N pieces of data.
步骤S120:将当前帧参考信号与历史帧参考信号进行组合,以获得当前帧参考信号集。Step S120: Combine the current frame reference signal and the historical frame reference signal to obtain the current frame reference signal set.
在一实施例中,将所获取的当前帧参考信号与已保存的历史帧参考信号进行组合,以获得当前帧参考信号集。其中,历史帧参考信号是指当前帧之前时刻所获取的参考信号。历史帧可以是指当前帧的前一时刻,或者是当前帧之前的任意时刻。In one embodiment, the acquired current frame reference signal and the saved historical frame reference signal are combined to obtain the current frame reference signal set. The historical frame reference signal refers to a reference signal obtained at a moment before the current frame. The historical frame may refer to the previous moment of the current frame, or any moment before the current frame.
每一帧参考信号均具有与其对应的参考信号集。Each frame of reference signal has its corresponding reference signal set.
步骤S130:基于当前帧滤波系数和当前帧参考信号对当前帧输入信号进行自适应滤波。Step S130: Perform adaptive filtering on the current frame input signal based on the current frame filter coefficient and the current frame reference signal.
在一实施例中,滤波过程可以为:利用当前帧滤波系数与当前帧参考信号进行卷积,获得滤波卷积信号。将当前帧输入信号减去滤波卷积信号,即为滤波后的信号。In one embodiment, the filtering process may be as follows: using the current frame filter coefficient and the current frame reference signal to perform convolution to obtain a filtered convolution signal. The filtered convolution signal is subtracted from the current frame input signal, that is, the filtered signal.
在一实施例中,可以将滤波后的信号传送至扬声器或网络端等进行输出。In one embodiment, the filtered signal can be sent to a speaker or a network terminal for output.
步骤S140:基于当前帧参考信号集与当前帧输入信号对当前帧滤波系数按照预设的划分长度进行更新,以获得当前帧之后的后续帧的滤波系数。Step S140: Based on the current frame reference signal set and the current frame input signal, the filter coefficient of the current frame is updated according to the preset division length to obtain the filter coefficient of the subsequent frame after the current frame.
在一实施例中,滤波系数更新方法可以包括获取划分长度;基于划分长度,更新当前帧滤波系数。其中划分长度与每帧自适应滤波方法获取的输入信号的长度满足约束条件。其中,约束条件可以是指划分长度与每帧输入信号长度之间需要满足的关系。划分长度可以是预设定的,也可以是从别处获取得到的。获得当前帧之后的后续帧的滤波系数可以是指获得当前时刻后一时刻的滤波系数。In one embodiment, the filter coefficient update method may include obtaining the division length; and updating the current frame filter coefficient based on the division length. The division length and the length of the input signal obtained by the adaptive filtering method of each frame satisfy the constraints. The constraint condition may refer to a relationship that needs to be satisfied between the division length and the length of the input signal of each frame. The division length can be preset or obtained from elsewhere. Obtaining filter coefficients of subsequent frames after the current frame may refer to obtaining filter coefficients at a moment after the current moment.
在一实施例中,约束条件可以是:当前帧参考信号集和历史帧参考信号集中包含数据相同的至少两个部分,历史帧和当前帧之间的间隔帧数与划分长度的乘积等于至少两个部分的差值与当前帧输入信号长度的乘积。具体的,当前帧参考信号集与当前时刻前的N时刻(例如,前一时刻,前二时刻等)的历史帧参考信号集中包含数据相同的至少两个部分。划分长度与N的乘积等于这两个部分的差值与当前帧输入信号长度的乘积。其中这两个部分的差值是指这两个部分在对应的参考信号集中编号的差值。输入信号长度是指每帧输入信号中包含的数据量。In one embodiment, the constraints may be: the current frame reference signal set and the historical frame reference signal set contain at least two parts with the same data, and the product of the number of interval frames and the division length between the historical frame and the current frame is equal to at least two parts. The product of the difference of the parts and the length of the input signal of the current frame. Specifically, the current frame reference signal set and the historical frame reference signal set at N moments before the current moment (eg, the previous moment, the previous two moments, etc.) include at least two parts with the same data. The product of the division length and N is equal to the product of the difference between these two parts and the length of the input signal of the current frame. The difference between the two parts refers to the difference between the numbers of the two parts in the corresponding reference signal set. The input signal length refers to the amount of data contained in each frame of the input signal.
本申请所公开的自适应滤波方法,通过设定滤波系数的更新中的划分长度,并使得划分长度与当前帧输入信号长度满足约束关系。使得在每次参考系数的更新过程中可以将当前帧输入信号全部进行计算,从而避免了部分信号数据滞留。进而,每次滤波系数更新的计算量相等,避免了计算量不均匀的问题,还可以有效的提高更新算法的更新效果。In the adaptive filtering method disclosed in the present application, the division length in the updating of the filter coefficient is set, and the division length and the length of the input signal of the current frame satisfy the constraint relationship. This makes it possible to calculate all the input signals of the current frame in the updating process of each reference coefficient, thereby avoiding the retention of part of the signal data. Furthermore, the calculation amount of each filter coefficient update is equal, which avoids the problem of uneven calculation amount, and can effectively improve the update effect of the update algorithm.
滤波系数的更新过程实质在于寻找一个有限脉冲响应,使得该响应与参考信号卷积所得到的卷积信号与目标回声信号的差最小。即,滤波后的输出信号中回声信号最小。具体地,可以参见公式(1):The essence of the filter coefficient update process is to find a finite impulse response, so that the difference between the convolution signal obtained by convolving the response with the reference signal and the target echo signal is the smallest. That is, the echo signal is the smallest in the filtered output signal. Specifically, you can refer to formula (1):
其中:d(i)为输入信号,x(i)为参考信号,为脉冲响应,即滤波系数。Where: d(i) is the input signal, x(i) is the reference signal, is the impulse response, that is, the filter coefficient.
参阅图2,图2是根据本申请一实施例的滤波系数更新方法的流程示意图。需注意的是,若有实质上相同的结果,本实施例并不以图2所示的流程顺序为限。如图2所示,该方法包括:Referring to FIG. 2, FIG. 2 is a schematic flowchart of a method for updating a filter coefficient according to an embodiment of the present application. It should be noted that, if there is substantially the same result, the present embodiment is not limited to the sequence of the processes shown in FIG. 2 . As shown in Figure 2, the method includes:
步骤S210:获取划分长度。Step S210: Obtain the division length.
在一实施例中,划分长度P可以是预设定的,也可以是其他模块或网络端获取得到的。划分长度P可以为大于零的正整数。划分长度P与每帧自适应滤波方法获取的输入信号的长度L满足约束条件。In an embodiment, the division length P may be preset, or may be obtained by other modules or the network. The division length P may be a positive integer greater than zero. The division length P and the length L of the input signal obtained by the adaptive filtering method of each frame satisfy the constraints.
划分长度P的确定方法的更多细节可以参见图3及其相关描述。For more details of the method for determining the division length P, please refer to FIG. 3 and its related descriptions.
步骤S220:获取更新信号块长度。Step S220: Obtain the update signal block length.
在每次进行滤波系数更新时,要求进行计算的信号块长度M为2的幂数。在一实施例中,更新信号块长度M可以是指进行更新计算的参考信号的数据块的数据点数,也可以是指进行更新计算的输入信号的数据块的数据点数。Each time the filter coefficient is updated, the signal block length M required to be calculated is a power of 2. In an embodiment, the block length M of the update signal may refer to the number of data points in the data block of the reference signal for which the update calculation is performed, or may refer to the number of data points of the data block of the input signal for which the update calculation is performed.
在一实施例中,更新信号块长度M可以是基于划分长度P以及每帧输入数据长度L确定的。具体的,更新信号块长度M需要满足条件M≥P+L-1。In an embodiment, the update signal block length M may be determined based on the division length P and the input data length L of each frame. Specifically, the update signal block length M needs to satisfy the condition M≥P+L-1.
步骤S230:获得多个参考信号分块。Step S230: Obtain multiple reference signal blocks.
在一实施例中,可以将当前帧参考信号与已保存的历史帧参考信号按照更新信号块长度M进行组合。从而,可以获得多个用于滤波系数更新计算的参考信号分块。其中,每个参考信号分块中的数据量等于更新信号块长度M。In an embodiment, the current frame reference signal and the saved historical frame reference signal may be combined according to the update signal block length M. Thus, a plurality of reference signal blocks for filter coefficient update calculation can be obtained. Wherein, the amount of data in each reference signal block is equal to the length M of the update signal block.
当前帧参考信号集包括多个当前帧参考信号分块。The current frame reference signal set includes a plurality of current frame reference signal blocks.
步骤S240:获得多个滤波系数分块。Step S240: Obtain a plurality of filter coefficient blocks.
在一实施例中,滤波系数可以基于划分长度P进行划分,从而获得若干滤波器系数分块。具体参见公式(2):In an embodiment, the filter coefficients may be divided based on the division length P, thereby obtaining several filter coefficient blocks. See formula (2) for details:
其中,为滤波系数分块,N为滤波系数的长度。in, is the filter coefficient block, and N is the length of the filter coefficient.
滤波系数的数据量N越大,滤波效果会越好。但是在滤波系数更新过程中,每一次更新计算均需要将滤波系数与参考信号进行卷积后,再与输入信号进行比较,因此更新计算中对于滤波系数的数据量会有限定。通过设定划分长度,将滤波系数分成多个滤波系数分块后在进行更新计算,可以使得滤波系数的数据量不受上述限制。The larger the data amount N of the filter coefficients, the better the filtering effect. However, in the filter coefficient update process, each update calculation needs to convolve the filter coefficient with the reference signal, and then compare it with the input signal, so the data amount of the filter coefficient in the update calculation will be limited. By setting the division length, the filter coefficient is divided into a plurality of filter coefficient blocks and then updated and calculated, so that the data amount of the filter coefficient is not limited above.
在一实施例中,滤波系数分块与参考信号分块是一一对应的。即,每一滤波系数分块均有与其对应的参考信号分块。In one embodiment, the filter coefficient blocks and the reference signal blocks are in a one-to-one correspondence. That is, each filter coefficient block has its corresponding reference signal block.
步骤S250:对滤波系数分块以及参考信号分块进行傅里叶变换。Step S250: Perform Fourier transform on the filter coefficient blocks and the reference signal blocks.
在滤波器系数更新的方法中,可以先将数据由时域转换成频域,再进行计算。信号从时域变换到频域可以通过傅里叶变换实现。In the method of updating the filter coefficients, the data can be converted from the time domain to the frequency domain first, and then the calculation is performed. The transformation of a signal from the time domain to the frequency domain can be accomplished by Fourier transform.
在一实施例中,将滤波系数分块进行傅里叶变换,获得频域滤波系数分块。由于待变换的信号块长度为M,傅里叶变换矩阵可以为M*M的矩阵。滤波系数分块的数据量会小于更新信号块长度M,因此在进行傅里叶变换时可以利用补充0向量的方式使得频域相滤波系数分块的数据量满足更新信号块长度M。具体参见公式(3):In one embodiment, the filter coefficient blocks are subjected to Fourier transform to obtain frequency domain filter coefficient blocks. Since the length of the signal block to be transformed is M, the Fourier transform matrix can be an M*M matrix. The data volume of the filter coefficient block will be less than the length M of the update signal block, so when performing Fourier transform, the method of supplementing 0 vectors can be used to make the data volume of the frequency domain phase filter coefficient block meet the update signal block length M. See formula (3) for details:
其中,wp为频域相滤波系数分块,F表示M*M的傅里叶变换矩阵,0为1-(M-P)个零向量。Among them, w p is the frequency domain phase filter coefficient block, F represents the Fourier transform matrix of M*M, and 0 is 1-(MP) zero vectors.
在一实施例中,将参考信号分块进行傅里叶变换,并将参考信号分块向量转化成为矩阵,从而可以与滤波系数分块进行卷积计算。具体参见公式(4):In an embodiment, the reference signal block is subjected to Fourier transform, and the reference signal block vector is converted into a matrix, so that the convolution calculation can be performed with the filter coefficient block. See formula (4) for details:
Xp=diag{F[x[(n+1)L-pP-M+1],…,x[(n+1)L-pP]]} (4)X p =diag{F[x[(n+1)L-pP-M+1],...,x[(n+1)L-pP]]} (4)
其中,XP为频域滤波系数分块矩阵;diag为对角函数;x[(n+1)L-pP]为编号为“(n+1)L-pP”的参考信号数据。Wherein, X P is the frequency domain filter coefficient block matrix; diag is the diagonal function; x[(n+1)L-pP] is the reference signal data numbered "(n+1)L-pP".
步骤S260:获得卷积信号。Step S260: Obtain a convolution signal.
在一实施例中,可以构建信号转换矩阵,用于去掉计算矩阵中的冗余数据或错误数据,仅对部分数据进行计算。具体地,第一转换矩阵为其中IL为L*L个单位向量。通过设置第一转换矩阵,可以将计算的数据中的部分数据去除,从而提高运算效率。In one embodiment, a signal transformation matrix may be constructed to remove redundant data or erroneous data in the calculation matrix, and only perform calculation on part of the data. Specifically, the first transformation matrix is where IL is L*L unit vectors. By setting the first transformation matrix, part of the data in the calculated data can be removed, thereby improving the operation efficiency.
在一实施例中,可以将有效滤波系数分块与有效参考信号分块进行卷积计算,获得卷积信号。并且,需要将所获得的卷积信号从频域转化为时域。具体地,卷积信号矩阵可以参见公式(5):In an embodiment, a convolution calculation may be performed on the effective filter coefficient blocks and the effective reference signal blocks to obtain a convolution signal. And, the obtained convoluted signal needs to be transformed from the frequency domain to the time domain. Specifically, the convolution signal matrix can refer to formula (5):
步骤S270:构建输入信号分块。Step S270: Construct input signal blocks.
在一实施例中,可以通过添加0向量构建具有M个数据的输入信号分块。具体的参见公式(6):In one embodiment, an input signal block with M data may be constructed by adding a 0 vector. For details, see formula (6):
d=[0 dn]T (6)d=[0 d n ] T (6)
其中,d为输入信号分块;dn为第n帧输入信号,具体参见公式(7)Among them, d is the input signal block; d n is the input signal of the nth frame, see formula (7) for details.
dn=[d[nL+1],…,d[(n+1)L]] (7)d n =[d[nL+1],...,d[(n+1)L]] (7)
步骤S280:利用卷积信号和输入信号分块,对滤波系数进行更新。Step S280: Use the convolution signal and the input signal to block, and update the filter coefficient.
在一实施例中,将输入信号分块与卷积信号矩阵相减,获得差值e,参见公式(8):In one embodiment, the input signal block is subtracted from the convolution signal matrix to obtain the difference e, see formula (8):
e=d-y (8)e=d-y (8)
在一实施例中,可以将差值进行傅里叶变换后,与处理过的参考信号分块相乘,再与当前帧滤波系数相加,获得下一滤波系数。具体参见公式(9):In one embodiment, the difference value may be subjected to Fourier transform, multiplied by the processed reference signal in blocks, and then added to the filter coefficient of the current frame to obtain the next filter coefficient. See formula (9) for details:
其中为下一滤波系数,为当前帧滤波系数,F为傅里叶变换矩阵,为Xp的共轭矩阵,G为约束变换矩阵,Δ为归一化系数。具体地,G和Δ参见公式(10)和(11):in is the next filter coefficient, is the filter coefficient of the current frame, F is the Fourier transform matrix, is the conjugate matrix of X p , G is the constraint transformation matrix, and Δ is the normalization coefficient. Specifically, G and Δ refer to equations (10) and (11):
Δ=1/(βΔ+(1-β)|Xp|2) (11)Δ=1/(βΔ+(1-β)|X p | 2 ) (11)
其中,β为平滑参数,可以取0.92。为第二转换矩阵。第二转换矩阵可以用于去掉计算矩阵中的冗余数据或错误数据,仅对部分数据进行计算。Among them, β is a smoothing parameter, which can be taken as 0.92. is the second transformation matrix. The second transformation matrix can be used to remove redundant data or erroneous data in the calculation matrix, and only perform calculation on part of the data.
由上述滤波系数的更新方法可见,滤波系数的划分长度会直接影响更新过程的运算量。参阅图3,图3是根据本申请一实施例的划分长度确定方法的流程示意图。需注意的是,若有实质上相同的结果,本实施例并不以图3所示的流程顺序为限。如图3所示,该方法包括:It can be seen from the above update method of the filter coefficient that the division length of the filter coefficient will directly affect the calculation amount of the update process. Referring to FIG. 3 , FIG. 3 is a schematic flowchart of a method for determining a division length according to an embodiment of the present application. It should be noted that, if there is substantially the same result, the present embodiment is not limited to the sequence of the processes shown in FIG. 3 . As shown in Figure 3, the method includes:
步骤S310:确定约束关系。Step S310: Determine the constraint relationship.
为了使得卷积操作能够有效,需要满足条件M≥P+L-1。同时,为了提高运算效率,可以将更新当前帧参考信号时所计算的数据(例如,Xp)在后续参考系数更新时重复使用,从而能够减少运算量。In order for the convolution operation to be effective, the condition M≥P+L-1 needs to be satisfied. Meanwhile, in order to improve the operation efficiency, the data (eg, X p ) calculated when the reference signal of the current frame is updated can be reused in the subsequent update of the reference coefficient, so that the amount of operation can be reduced.
在一实施例中,当前帧和历史帧的参考信号集中包含数据相同的至少两个参考信号分块,分别为第一参考信号分块和第二参考信号分块。即,前述当前帧参考信号集和历史帧参考信号集中包含数据相同的至少两个部分即可以是指数据相同的至少两个参考信号分块。第一参考信号分块和第二参考信号分块中的所采用的参考信号数据需要满足如公式(12)的条件,才可以在多帧的滤波系数更新中重复使用。In an embodiment, the reference signal sets of the current frame and the historical frame include at least two reference signal blocks with the same data, which are the first reference signal block and the second reference signal block, respectively. That is, the aforementioned reference signal set of the current frame and the reference signal set of the historical frame containing at least two parts with the same data may refer to at least two reference signal blocks with the same data. The reference signal data used in the first reference signal sub-block and the second reference signal sub-block needs to satisfy the conditions as in the formula (12) before they can be reused in the filter coefficient update of multiple frames.
其中,n1为历史帧滤波系数与当前帧滤波系数相差的帧数。Wherein, n1 is the number of frames by which the filter coefficient of the historical frame differs from the filter coefficient of the current frame.
公式(12)可以理解为,更新当前帧滤波系数(即第n+n1帧)时的第p1项参考信号与更新第n1帧前的滤波系数(即,第n帧)时的第p2项参考信号分块相同。Formula (12) can be understood as the p1th reference signal when the filter coefficient of the current frame (ie, the n+ n1th frame) is updated and the p2th reference signal when the filter coefficient before the n1th frame (ie, the nth frame) is updated. The item reference signal blocks are the same.
同时,p1和p2还需要满足p1>p2。Meanwhile, p1 and p2 also need to satisfy p 1 >p 2 .
简化公式(12)即可以获得划分长度与输入信号长度的约束关系。约束条件为:历史帧和当前帧之间的间隔帧数与划分长度的乘积等于第一参考信号分块和第二参考信号分块的差值与当前帧输入信号长度的乘积。参见公式(13):By simplifying the formula (12), the constraint relationship between the division length and the input signal length can be obtained. Constraints are: the product of the number of interval frames between the historical frame and the current frame and the division length is equal to the product of the difference between the first reference signal block and the second reference signal block and the length of the input signal of the current frame. See equation (13):
n1L=(p1-p2)P (13)n 1 L=(p 1 -p 2 )P (13)
步骤S320:确定参考信号分块的重复使用度。Step S320: Determine the degree of reuse of reference signal blocks.
在一实施例中,可以预先根据实际应用场景或计算设备的配置确定重复使用度。具体的,可以通过确定参考信号分块能够重复使用的数量最大的历史帧与当前帧的间隔帧数(例如,n1),以及能够重复使用的参考信号分块的数量(例如p1-p2的值)。In one embodiment, the degree of reuse may be determined in advance according to the actual application scenario or the configuration of the computing device. Specifically, the number of frames between the historical frame and the current frame with the largest number of reference signal blocks that can be reused (for example, n1), and the number of reference signal blocks that can be reused (for example, the value of p1-p2) can be determined. ).
步骤S330:确定划分长度以及更新信号块长度。Step S330: Determine the division length and update the signal block length.
在一实施例中,将前述确定的参考信号的重复使用度代入约束关系中,基于输入信号的长度可以获得滤波系数的划分长度。此外,可以根据条件M≥P+L-1,确定M的取值范围。并且由于M需要为2的幂数,因此M的取值可以为最接近P+L-1的2的幂数。In one embodiment, the previously determined reusability of the reference signal is substituted into the constraint relationship, and the division length of the filter coefficient can be obtained based on the length of the input signal. In addition, the value range of M can be determined according to the condition M≥P+L-1. And since M needs to be a power of 2, the value of M can be a power of 2 that is closest to P+L-1.
例如,设定p1-p2=2,n1=1时,即更新当前帧滤波系数使用的参考信号分块与更新前一历史帧滤波系数使用的参考信号分块仅有两项不相同。每帧输入信号长度L可以为160,将上述值代入约束关系中,获得L=2P,即划分长度P可以为80。M的取值需要满足3P+1=241<M,M可以为256。综上,在当前帧输入信号长度为160,划分长度为80时,将当前帧输入参考信号和历史帧参考信号以256长度进行组合,以获得多个当前帧参考信号分块。此方案中滤波系数单次更新可以将当前帧输入信号全部代入,从而有效解决数据滞留的问题。For example, when setting p 1 -p 2 =2, n 1 =1, that is, the reference signal block used for updating the filter coefficients of the current frame is different from the reference signal block used for updating the filter coefficients of the previous historical frame only by two items . The length L of the input signal of each frame may be 160, and the above value is substituted into the constraint relationship to obtain L=2P, that is, the division length P may be 80. The value of M needs to satisfy 3P+1=241<M, and M can be 256. To sum up, when the current frame input signal length is 160 and the division length is 80, the current frame input reference signal and the historical frame reference signal are combined with a length of 256 to obtain multiple current frame reference signal blocks. In this scheme, a single update of the filter coefficients can substitute all the input signals of the current frame, thereby effectively solving the problem of data retention.
又例如,设定p1-p2=1,n1=2时,即更新当前帧滤波系数使用的参考信号分块与更新当前之前的第二帧滤波系数使用的参考信号分块仅有一项不相同。每帧输入信号长度L可以为160,将上述值代入约束关系中,获得2L=P,即P可以为320。M的取值需要满足3L-1=479<M,M可以为512。综上,在当前帧输入信号长度为160,划分长度为320时,将当前帧输入参考信号和历史帧参考信号以512长度进行组合,以获得多个当前帧参考信号分块。此方案中滤波系数单次更新既可以将当前帧输入信号全部代入,从而有效解决数据滞留的问题。并且,能够降低整个自适应滤波器计算的复杂度。For another example, when setting p 1 -p 2 =1, n 1 =2, there is only one reference signal block used for updating the filter coefficients of the current frame and the reference signal block used for updating the filter coefficients of the second frame before the current frame Are not the same. The length L of the input signal of each frame can be 160, and the above value is substituted into the constraint relationship to obtain 2L=P, that is, P can be 320. The value of M needs to satisfy 3L-1=479<M, and M can be 512. To sum up, when the length of the current frame input signal is 160 and the division length is 320, the current frame input reference signal and the historical frame reference signal are combined with a length of 512 to obtain multiple current frame reference signal blocks. In this scheme, a single update of the filter coefficients can substitute all the input signals of the current frame, thereby effectively solving the problem of data retention. Also, the computational complexity of the entire adaptive filter can be reduced.
请参阅图,图4是根据本申请一实施例的自适应滤波器的结构示意图。该实施例中,自适应滤波器400包括获取模块410,组合模块420,滤波模块430以及更新模块440。Please refer to the drawings. FIG. 4 is a schematic structural diagram of an adaptive filter according to an embodiment of the present application. In this embodiment, the adaptive filter 400 includes an acquiring
在一实施例中,获取模块410用于获取当前帧输入信号、当前帧参考信号和当前帧滤波系数。In one embodiment, the obtaining
在一实施例中,组合模块420用于将当前帧参考信号与历史帧参考信号进行组合,以获得当前帧参考信号集。组合模块420还可以用于将当前帧参考信号与历史帧参考信号进行组合,以获得当前帧参考信号集包括:将当前帧参考信号与历史帧参考信号按照预设更新信号块长度进行组合,以获得多个当前帧参考信号分块;当前帧参考信号集包括多个当前帧参考信号分块;其中,更新信号块长度由划分长度与当前帧输入信号长度所确定。更新信号块长度不小于划分长度与每帧输入信号长度之和减1的值。更新信号块长度等于2的幂数。In one embodiment, the combining
在一实施例中,滤波模块430用于基于当前帧滤波系数和当前帧参考信号对当前帧输入信号进行自适应滤波。In one embodiment, the
在一实施例中,更新模块440用于基于当前帧参考信号和当前帧滤波系数对当前帧输入信号进行自适应滤波,基于当前帧参考信号集与当前帧输入信号对当前帧滤波系数按照预设的划分长度进行更新,以获得后续帧的滤波系数;其中,当前帧参考信号集和历史帧参考信号集中包含数据相同的至少两个部分,历史帧和当前帧之间的间隔帧数与划分长度的乘积等于至少两个部分的差值与当前帧输入信号长度的乘积。In one embodiment, the
更新模块440还可以用于将当前帧滤波系数按照划分长度进行划分,以获得多个当前帧滤波系数分块,其中当前帧滤波系数分块与当前帧参考信号分块一一对应;对当前帧滤波系数分块以及对应的当前帧参考信号分块进行卷积,以获得卷积信号;基于当前帧卷积信号与当前帧输入信号的差,更新当前帧滤波系数。其中,当前帧和历史帧的参考信号集包含数据相同的至少两个参考信号分块,分别为第一参考信号分块和第二参考信号分块;历史帧和当前帧之间的间隔帧数与划分长度的乘积等于第二参考信号分块和第一参考信号分块的差值与当前帧输入信号长度的乘积。其中,划分长度为每帧输入信号长度的二分之一;或划分长度为每帧输入信号长度的二倍。其中,将当前帧参考信号与历史帧参考信号按照预设更新信号块长度进行组合,以获得多个当前帧参考信号分块的步骤包括:当当前帧输入信号长度为160,划分长度为80时,将当前帧输入参考信号和历史帧参考信号以256长度进行组合,以获得多个当前帧参考信号分块;或当当前帧输入信号长度为160,划分长度为320时,将当前帧输入参考信号和历史帧参考信号以512长度进行组合,以获得多个当前帧参考信号分块。The
请参阅图5,图5是根据本申请一实施例的计算机设备的结构示意图。该实施方式中,计算机设备500包括处理器510。Please refer to FIG. 5 , which is a schematic structural diagram of a computer device according to an embodiment of the present application. In this embodiment,
处理器510还可以称为CPU(Central Processing Unit,中央处理单元)。处理器510可能是一种集成电路芯片,具有信号的处理能力。处理器510还可以是通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。The
计算机设备500可以进一步包括存储器(图中未示出),用于存储处理器510运行所需的指令和数据。
处理器510用于执行指令以实现上述自适应滤波方法任一实施例及任意不冲突的组合所提供的方法。The
请参阅图6,图6是根据本申请一实施例的具有存储功能的装置的结构示意图。本申请实施例的具有存储功能的装置600存储有指令,该指令被执行时实现本申请自适应滤波方法任一实施例以及任意不冲突的组合所提供的方法。其中,该指令可以形成程序文件以软件产品的形式存储在上述具有存储功能的装置中,以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)或处理器(processor)执行本申请各个实施方式方法的全部或部分步骤。而前述的具有存储功能的装置600包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质,或者是计算机、服务器、手机、平板等终端设备。Please refer to FIG. 6 , which is a schematic structural diagram of a device with a storage function according to an embodiment of the present application. The
在本申请所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。In the several embodiments provided in this application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. For example, the apparatus embodiments described above are only illustrative. For example, the division of units is only a logical function division. In actual implementation, there may be other division methods, for example, multiple units or components may be combined or integrated. to another system, or some features can be ignored, or not implemented. On the other hand, the shown or discussed mutual coupling or direct coupling or communication connection may be through some interfaces, indirect coupling or communication connection of devices or units, and may be in electrical, mechanical or other forms.
另外,在本申请各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically alone, or two or more units may be integrated into one unit. The above-mentioned integrated units may be implemented in the form of hardware, or may be implemented in the form of software functional units.
以上所述仅为本申请的实施方式,并非因此限制本申请的专利范围,凡是利用本申请说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本申请的专利保护范围内。The above description is only an embodiment of the present application, and is not intended to limit the scope of the patent of the present application. Any equivalent structure or equivalent process transformation made by using the contents of the description and drawings of the present application, or directly or indirectly applied to other related technologies Fields are similarly included within the scope of patent protection of this application.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010968471.7A CN112260662B (en) | 2020-09-15 | 2020-09-15 | A method for adaptive filtering, computer equipment and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010968471.7A CN112260662B (en) | 2020-09-15 | 2020-09-15 | A method for adaptive filtering, computer equipment and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112260662A true CN112260662A (en) | 2021-01-22 |
CN112260662B CN112260662B (en) | 2025-04-01 |
Family
ID=74232899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010968471.7A Active CN112260662B (en) | 2020-09-15 | 2020-09-15 | A method for adaptive filtering, computer equipment and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112260662B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007174632A (en) * | 2005-11-25 | 2007-07-05 | Nippon Telegr & Teleph Corp <Ntt> | Apparatus, method and program for estimating transmission system, and recording medium |
CN104883462A (en) * | 2014-02-28 | 2015-09-02 | 中国科学院声学研究所 | Adaptive filter and filtering method used for eliminating acoustics echo |
CN105472191A (en) * | 2015-11-18 | 2016-04-06 | 百度在线网络技术(北京)有限公司 | Method and device for tracking echo time delay |
CN106059967A (en) * | 2016-07-01 | 2016-10-26 | 晶晨半导体(上海)有限公司 | Single carrier channel estimating method |
US20170178617A1 (en) * | 2015-12-17 | 2017-06-22 | Harman Becker Automotive Systems Gmbh | Active noise control by adaptive noise filtering |
CN107039043A (en) * | 2017-06-08 | 2017-08-11 | 腾讯科技(深圳)有限公司 | The method and device of signal transacting, the method and system of multi-conference |
CN107123430A (en) * | 2017-04-12 | 2017-09-01 | 广州视源电子科技股份有限公司 | Echo cancellation method, device, conference tablet and computer storage medium |
CN109559756A (en) * | 2018-10-26 | 2019-04-02 | 北京佳讯飞鸿电气股份有限公司 | Filter factor determines method, echo cancel method, related device and equipment |
CN111031448A (en) * | 2019-11-12 | 2020-04-17 | 西安讯飞超脑信息科技有限公司 | Echo cancellation method, echo cancellation device, electronic equipment and storage medium |
CN111445917A (en) * | 2020-03-17 | 2020-07-24 | 浙江大华技术股份有限公司 | Echo cancellation method, device and computer storage medium |
-
2020
- 2020-09-15 CN CN202010968471.7A patent/CN112260662B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007174632A (en) * | 2005-11-25 | 2007-07-05 | Nippon Telegr & Teleph Corp <Ntt> | Apparatus, method and program for estimating transmission system, and recording medium |
CN104883462A (en) * | 2014-02-28 | 2015-09-02 | 中国科学院声学研究所 | Adaptive filter and filtering method used for eliminating acoustics echo |
CN105472191A (en) * | 2015-11-18 | 2016-04-06 | 百度在线网络技术(北京)有限公司 | Method and device for tracking echo time delay |
US20170178617A1 (en) * | 2015-12-17 | 2017-06-22 | Harman Becker Automotive Systems Gmbh | Active noise control by adaptive noise filtering |
CN106059967A (en) * | 2016-07-01 | 2016-10-26 | 晶晨半导体(上海)有限公司 | Single carrier channel estimating method |
CN107123430A (en) * | 2017-04-12 | 2017-09-01 | 广州视源电子科技股份有限公司 | Echo cancellation method, device, conference tablet and computer storage medium |
CN107039043A (en) * | 2017-06-08 | 2017-08-11 | 腾讯科技(深圳)有限公司 | The method and device of signal transacting, the method and system of multi-conference |
CN109559756A (en) * | 2018-10-26 | 2019-04-02 | 北京佳讯飞鸿电气股份有限公司 | Filter factor determines method, echo cancel method, related device and equipment |
CN111031448A (en) * | 2019-11-12 | 2020-04-17 | 西安讯飞超脑信息科技有限公司 | Echo cancellation method, echo cancellation device, electronic equipment and storage medium |
CN111445917A (en) * | 2020-03-17 | 2020-07-24 | 浙江大华技术股份有限公司 | Echo cancellation method, device and computer storage medium |
Non-Patent Citations (3)
Title |
---|
MAHBUB U等: "Single-channel acoustic echo cancellation in noise based on gradient-based adaptive filtering", 《EURASIP JOURNAL ON AUDION, SPEECH, AND MUSIC PROCESSING》, 3 May 2014 (2014-05-03) * |
台宏达;: "语音通信系统中的自适应回声消除技术", 计算机仿真, no. 09, 15 September 2015 (2015-09-15) * |
蔡宇;洪缨;原建平;侯朝焕;郝程鹏;: "语音系统中的子带自适应回声消除技术", 仪器仪表学报, no. 07, 15 July 2013 (2013-07-15) * |
Also Published As
Publication number | Publication date |
---|---|
CN112260662B (en) | 2025-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109686381B (en) | Signal processor for signal enhancement and related method | |
CN111445917B (en) | Echo cancellation method, device and computer storage medium | |
CN102132491B (en) | Method for determining updated filter coefficients of an adaptive filter adapted by an lms algorithm with pre-whitening | |
CN110265054B (en) | Speech signal processing method, device, computer readable storage medium and computer equipment | |
CN105391879A (en) | Echo residue-free double-end communication robust acoustic echo elimination method | |
JP2003503871A (en) | Acoustic echo and noise removal | |
CN111489761A (en) | Residual echo estimator, non-transitory computer readable medium and application processor | |
CN112017680B (en) | Dereverberation method and device | |
CN109215672B (en) | Method, device and equipment for processing sound information | |
CN114242100B (en) | Audio signal processing method, training method, device, equipment and storage medium thereof | |
EP1314247B1 (en) | Partitioned block frequency domain adaptive filter | |
Albu et al. | New iterative kernel algorithms for nonlinear acoustic echo cancellation | |
CN114495960A (en) | Audio noise reduction filtering method, noise reduction filtering device, electronic equipment and storage medium | |
CN112260662A (en) | Adaptive filtering method, computer equipment and device | |
CN108141202B (en) | A block frequency domain adaptive filter device including an adaptive module and a correction module | |
KR20220157475A (en) | Echo Residual Suppression | |
CN112803921B (en) | Adaptive filter, method, medium, and electronic device | |
CN113362844B (en) | Low-complexity decorrelation self-adaptive acoustic echo cancellation method and device | |
Enzner et al. | Hybrid-frequency-resolution adaptive Kalman filter for online identification of long acoustic responses with low input-output latency | |
CN115706757A (en) | Acoustic echo cancellation using control parameters | |
CN113783551A (en) | Filter coefficient determination method, echo cancellation method and device | |
WO2016074422A1 (en) | Audio signal processing method and related device | |
JP7495684B2 (en) | Echo cancellation method, device and electronic device | |
JP6343585B2 (en) | Unknown transmission system estimation device, unknown transmission system estimation method, and program | |
CN119785809A (en) | Speech dereverberation method, device, computer equipment, readable storage medium and program product |
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 |