WO2017012350A1 - Procédé et dispositif de détermination de divergence d'état de filtre - Google Patents

Procédé et dispositif de détermination de divergence d'état de filtre Download PDF

Info

Publication number
WO2017012350A1
WO2017012350A1 PCT/CN2016/075804 CN2016075804W WO2017012350A1 WO 2017012350 A1 WO2017012350 A1 WO 2017012350A1 CN 2016075804 W CN2016075804 W CN 2016075804W WO 2017012350 A1 WO2017012350 A1 WO 2017012350A1
Authority
WO
WIPO (PCT)
Prior art keywords
filter coefficient
coefficient vector
main reflection
end signal
value
Prior art date
Application number
PCT/CN2016/075804
Other languages
English (en)
Chinese (zh)
Inventor
苏谟特艾雅
刘媛媛
李海婷
Original Assignee
华为技术有限公司
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 华为技术有限公司 filed Critical 华为技术有限公司
Publication of WO2017012350A1 publication Critical patent/WO2017012350A1/fr

Links

Images

Classifications

    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00—Adaptive networks
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B3/00—Line transmission systems
    • H04B3/02—Details
    • H04B3/20—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other
    • H04B3/23—Reducing echo effects or singing; Opening or closing transmitting path; Conditioning for transmission in one direction or the other using a replica of transmitted signal in the time domain, e.g. echo cancellers
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04M—TELEPHONIC COMMUNICATION
    • H04M1/00—Substation equipment, e.g. for use by subscribers
    • H04M1/02—Constructional features of telephone sets
    • H04M1/20—Arrangements for preventing acoustic feed-back

Definitions

  • the present invention relates to the field of communications, and in particular, to a method and apparatus for determining a state of a filter divergence.
  • the call quality is affected by the echo. If the echo cancellation is not performed during the call, especially under the hands-free condition, the normal call cannot be achieved.
  • the so-called echo refers to the situation that the voice signal from the near-end speaker is directly or indirectly collected by the microphone and transmitted back to the far end, that is, the party in the call hears its own voice from the phone used by itself.
  • echo is typically cancelled using adaptive filter techniques.
  • the adaptive filter uses the far-end signal to simulate the path generated by the echo and generates an estimated echo signal, which is then subtracted from the near-end signal and sent to the far end.
  • FIG. 1 is a schematic diagram of the working principle of the echo canceller in the background art of the present invention.
  • One straight line and two curves in the right solid line block diagram in FIG. 1 are used to simulate the sound of the far end from the speaker. And the process of propagation collected by the Mike. Among them, the curve indicates that the sound of the far end is received by the microphone through reflection.
  • the dotted line frame in Fig. 1 is a typical acoustic echo canceller. Assuming that user A and user B are talking, the near-end signal d(n) collected by the microphone can be understood as the voice of user A's own voice and user B's.
  • Echo if user A does not speak, only the echo of user B, and the far-end signal x(n) transmitted by the remote end to synchronize with the near-end signal can be understood as the voice of the user B who is talking with user A, Input two signals into the echo canceller, according to the original adaptive filter
  • the filter coefficient vector can be calculated, and the echo propagation path in the right solid line block diagram in FIG. 1 is simulated according to the filter coefficient vector, and the output signal e(n) processed by the echo canceller is calculated. Further, the output signal e(n) is transmitted back to the remote end, that is, the user B, and the entire echo cancellation process is completed to ensure the quality of the call between the users.
  • FIG. 2 is a schematic diagram of the working principle of the adaptive filter in the background art of the present invention, wherein the input signal of the adaptive filter is the near-end signal d(n) and is synchronized with the near-end signal.
  • the far-end signal x(n) is obtained by adaptive filter echo cancellation processing and the output signal is e(n), and the output signal is sent to the far end. Further, whether the filter state diverges or not determines the quality of the output signal e(n). If the filter state is diverged, it indicates that the output signal is distorted, and the echo canceller fails, which will seriously affect the communication quality.
  • the detection of the filter state is an important component of the echo cancellation system and even the communication system, and is an important criterion for judging the quality of the output signal and the quality of the communication.
  • ERLE Echo Return Loss Enhancement
  • E[d 2 (n)] represents the near-end signal energy collected by the microphone, the near-end signal includes echo and other signals;
  • E[e 2 (n)] represents the output signal energy of the echo canceller, and the output signal Includes other signals and residual echo as described previously.
  • the other signals mentioned here are the voice signals spoken by the user in the actual call, and the ambient noise at the near end. Therefore, the other signals are not affected by default during the processing.
  • the filter state When the filter state converges, ERLE is relatively stable. When the adaptive filter is disturbed and the filter state is diverged, ERLE will be significantly reduced. Specifically, the filter can be determined by setting the threshold. Whether the state of the device converges, when the ERLE is higher than the preset threshold, the filter state converges. When the ERLE is less than the preset threshold, the filter state is diverged.
  • the prior art is based on the energy of the signal to evaluate whether the filter state is divergent.
  • Calculating the ERLE needs to calculate the energy difference between the near-end signal and the output signal, but since the speech signal itself fluctuates greatly, the signal energy fluctuation caused by the signal is also relatively large.
  • the calculation of ERLE in order to smooth the volatility of the signal, additional operations are needed to smooth the signal, and in order to control the dynamic range of the ERLE, the calculation needs to be converted into a logarithmic domain. Therefore, the complexity of the whole calculation process is more complicated. high.
  • the environment of the near-end signal changes, such as noise changes, it will also affect the stability of ERLE, which may cause errors in the judgment of whether the filter state is divergent.
  • the embodiment of the invention provides a method and a device for determining the state of the filter to be diverged, which is used to solve the problem that the filter state judging method existing in the prior art has high computational complexity and may cause a judgment error.
  • an embodiment of the present invention provides a method for determining a state of a filter divergence, including:
  • each pair of sync frames including a near-end signal frame and a far-end signal synchronized with the near-end signal frame frame;
  • M is a preset value, and both N and M are positive integers, and N ⁇ M.
  • the performing the synchronous framing processing on the near-end signal and the far-end signal respectively to obtain the N-pair synchronization frame includes:
  • the near-end signals are sequenced according to the signal timing sequence to obtain N near-end signal frames;
  • the far-end signals synchronized with the near-end signals are sequentially processed according to the signal timing sequence to obtain N far-end signal frames synchronized with the near-end signal frame;
  • Any one of the obtained near-end signal frame and the far-end signal frame synchronized with the any one of the near-end signal frames constitutes a pair of synchronization frames, and a total of N pairs of synchronization frames are obtained.
  • the N pairs of synchronization frames are filtered by a filter to obtain N filter coefficient vector mean values corresponding to the N pairs of synchronization frames, include:
  • the calculating N times of the corresponding N pairs of synchronization frames according to the N filter coefficient vector mean values Filter coefficient status reference, including:
  • the elements in the current filter coefficient vector mean are at least divided into a main reflection area element, a left side element of the main reflection area, and a right side element of the main reflection area, wherein the main reflection area
  • the element refers to all the elements contained in the preset sliding window in the filter coefficient distribution map.
  • the calculating N corresponding to N pairs of synchronization frames according to N filter coefficient vector averages Filter coefficient status reference including:
  • the elements in the current filter coefficient vector mean are at least divided into a main reflection area element, a left side element of the main reflection area, and a right side element of the main reflection area, wherein the main reflection area
  • An element is an element included in the preset sliding window when a sum of squares of all elements included in a preset sliding window in the filter coefficient distribution map is a maximum value, and the pre-predetermined in the filter coefficient distribution map
  • the element on the left side of the sliding window is an element on the left side of the main reflection area
  • the element on the right side of the preset sliding window is an element on the right side of the main reflection area
  • the element in the current filter coefficient vector mean is determined according to the filter coefficient profile At least the main reflection zone element, the left element of the main reflection zone, and the right element of the main reflection zone, including:
  • the method further includes:
  • an embodiment of the present invention provides an apparatus for determining a state of a filter, including:
  • An acquisition module configured to collect a near-end signal and a far-end signal synchronized with the near-end signal
  • a framing module configured to perform synchronous framing processing on the near-end signal and the far-end signal collected by the collecting module, respectively, to obtain N pairs of synchronization frames, each pair of synchronization frames including a near-end signal frame and a a far-end signal frame synchronized with the near-end signal frame;
  • a filtering module configured to filter, by using a filter, the N pairs of synchronization frames obtained by the framing module to obtain N filter coefficient vector mean values corresponding to the N pairs of synchronization frames;
  • a determining module configured to calculate, according to the N filter coefficient vector average values obtained by the filtering module, N filter coefficient state reference quantities corresponding to the N pairs of synchronization frames;
  • M is a preset value, and both N and M are positive integers, and N ⁇ M.
  • the framing module is specifically configured to:
  • the near-end signals collected by the collection module are sequentially processed according to the signal timing sequence to obtain N near-end signal frames;
  • the far-end signals that are collected by the acquiring module and synchronized with the near-end signal are subjected to frame processing in sequence according to signal timing, to obtain N synchronization with the near-end signal frame.
  • Remote signal frame
  • Any one of the obtained near-end signal frame and the far-end signal frame synchronized with the any one of the near-end signal frames constitutes a pair of synchronization frames, and a total of N pairs of synchronization frames are obtained.
  • the filtering module is specifically configured to:
  • the determining module is specifically configured to:
  • Entity of the current filter coefficient vector mean according to the filter coefficient profile
  • the determining module is further configured to:
  • the elements in the current filter coefficient vector mean are at least divided into a main reflection area element, a left side element of the main reflection area, and a right side element of the main reflection area, wherein the main reflection area
  • An element is an element included in the preset sliding window when a sum of squares of all elements included in a preset sliding window in the filter coefficient distribution map is a maximum value, and the pre-predetermined in the filter coefficient distribution map
  • the element on the left side of the sliding window is an element on the left side of the main reflection area
  • the element on the right side of the preset sliding window is an element on the right side of the main reflection area
  • the determining module is further configured to:
  • the device further includes:
  • a resetting module configured to stop updating the filter coefficient vector after the determining, according to the determining module, that the filter state is diverged, and clear the current filter coefficient vector.
  • a method for judging the state of the filter is proposed, which is: collecting a near-end signal and a far-end signal synchronized with the near-end signal, and performing synchronous framing on the near-end signal and the far-end signal, respectively. Processing, obtaining N pairs of synchronization frames, filtering the N pairs of synchronization frames through a filter, and obtaining N filter coefficient vector mean values corresponding to the N pairs of synchronization frames.
  • N filter coefficient state reference quantities corresponding to the N pairs of synchronization frames are respectively calculated, since the filter coefficient state reference quantity is insensitive to changes in signal energy, and is not affected by the near-end signal and The influence of the environment where the far-end signal is located, so the obtained result of the method is stable, and the calculation complexity of the filter coefficient state reference quantity is low, and no additional operation is needed to smooth the fluctuation of the signal.
  • at least M filter coefficient state reference quantities are greater than a preset threshold value among the N filter coefficient state reference quantities, it is determined that the filter state diverges.
  • the embodiment of the present invention uses the filter coefficient state reference quantity to determine whether the filter state is divergent, so as to achieve the purpose of controlling the echo canceler filter coefficient vector update and guaranteeing the call quality, and determining whether the filter state is diverged by the method of the present invention. It can effectively judge whether there is double talk during the call, so as to achieve the purpose of controlling the filter coefficient vector update and obtaining high quality output signal.
  • 1 is a schematic view showing the working principle of an echo canceller in the background art of the present invention
  • FIG. 2 is a schematic diagram showing the working principle of an adaptive filter in the background art of the present invention.
  • FIG. 3 is a flowchart showing an overview of determining a state of a filter divergence according to an embodiment of the present invention
  • 4A and 4B are diagrams showing a filter coefficient distribution diagram according to an embodiment of the present invention.
  • FIG. 5 is a schematic structural diagram of an apparatus for determining a state of a filter divergence according to an embodiment of the present invention
  • FIG. 6 is a schematic structural diagram of an apparatus for determining a state of a filter divergence according to an embodiment of the present invention.
  • the embodiment of the invention provides a method and a device for determining the state of the filter divergence, which are used to solve the problem that the filter state judging method existing in the prior art has high computational complexity and may cause a judgment error.
  • the method and the device are based on the same inventive concept. Since the principles of the method and the device for solving the problem are similar, the implementation of the device and the method can be referred to each other, and the repeated description is not repeated.
  • the call quality is affected by the echo. If the echo cancellation is not performed during the call, especially under the hands-free condition, the normal call cannot be achieved. In the prior art, it is usually used
  • the adaptive filter technique eliminates the echo, and the detection of the filter state is an important criterion for judging the echo cancellation effect. If the filter state is divergent, the echo cancellation effect is poor or even invalid, resulting in distortion of the output signal. If the filter state converges, the echo is mostly eliminated, and the call quality is good.
  • FIG. 3 is a flowchart of an overview of determining a state of divergence of a filter according to an embodiment of the present invention.
  • the specific process of determining a state of divergence of a filter in the present invention is as follows:
  • Step 300 Acquire a near-end signal and a far-end signal synchronized with the near-end signal.
  • the echo canceller combines the near-end signal collected by the microphone with the far-end signal synchronized with the near-end signal, which together serve as the acquired input signal.
  • the near-end signal includes echo and other signals, and other signals mentioned here are voice signals spoken by the user during actual conversation, and environmental noise at the near end.
  • the far-end signal is synchronized with the near-end signal, and the far-end signal is a voice signal from the near-end speaker.
  • synchronization refers to the fact that two signals maintain a certain relative relationship during the change process, and the phenomenon of consistency and unification occurs in time.
  • Step 310 Perform synchronous framing processing on the near-end signal and the far-end signal respectively to obtain N-pair synchronization frames, each pair of synchronization frames including a near-end signal frame and a far-end signal frame synchronized with the near-end signal frame.
  • the near-end signals are sequentially processed according to the signal timing sequence to obtain N near-end signal frames.
  • the far-end signals synchronized with the near-end signals are subjected to frame processing in sequence according to the signal timing, and N far-end signal frames synchronized with the near-end signal frames are obtained.
  • any one of the obtained near-end signal frames and the far-end signal frame synchronized with the any one of the near-end signal frames constitute a pair of synchronization frames, and a total of N pairs of synchronization frames are obtained.
  • the collected near-end signal is composed of 1000 sampling points
  • the preset frame length is 100 sampling points
  • the near-end signals are subjected to frame processing in sequence according to the signal timing, thereby obtaining 10 near-end signal frames
  • the collected far-end signal synchronized with the near-end signal is composed of 1000 sampling points.
  • the far-end signal synchronized with the near-end signal is subjected to frame processing according to the signal timing sequence. Get 10 far-end signal frames synchronized with the near-end signal frame.
  • any one of the obtained near-end signal frames and the far-end signal frame synchronized with the any one of the near-end signal frames constitute a pair of synchronization frames, and a total of 10 pairs of synchronization frames are obtained.
  • Step 320 Filtering the N pairs of synchronization frames by a filter to obtain N filter coefficient vector mean values corresponding to the N pairs of synchronization frames.
  • the filter coefficient vector mean value corresponding to the current pair of synchronization frames is obtained according to the filter coefficient vector update order T and the T-th filter coefficient vector corresponding to the current pair of synchronization frames to update the corresponding T update results.
  • an adaptive filter technique based on an affine projection algorithm is taken as an example, but is not limited to the adaptive filtering technique.
  • the following uses an affine projection algorithm to explain how to update the filter coefficient vector and how to obtain the filter coefficient vector update result.
  • the adaptive filtering technique of the affine projection algorithm it is set to collect one sample point at each time, and then construct a signal vector from multiple samples, assuming that at the nth time, according to the number of elements in the filter coefficient vector L, constructing the far-end signal input vector from the far-end signal x(n):
  • x(nl) [x(nl),...,x(nL-l+1)] T .
  • the far-end signal input matrix is constructed by the far-end signal input vector according to the preset affine projection order p:
  • X(n) [x(n),...,x(n-p+1)].
  • the resulting input matrix of the far-end signal can be expressed as:
  • the near-end signal input vector is constructed from the near-end signal d(n):
  • d(n) [d(n),...,d(n-p+1)] T .
  • a sample point is acquired at each moment, and the filter coefficient vector is updated once at each sampling point to obtain the filter coefficient vector updated at the n+1th time:
  • h(n+1) h(n)+ ⁇ X(n)[X T (n)X(n)] -1 e(n).
  • ⁇ is the update step size of the filter, which is the preset value
  • [*] -1 represents the matrix inversion operation
  • [*] T represents the matrix transposition operation
  • filter coefficient vector updates are not performed for each sample point during signal processing. Instead, the acquired near-end signal and the far-end signal synchronized with the near-end signal are synchronously framed.
  • the filter coefficients of the i-th sync frame can be obtained by using, but not limited to, the following method.
  • Vector mean
  • the filter coefficient vector is updated every 4 sample points, the filter coefficient vector will be updated 40 times, and each filter coefficient vector h k (i) is obtained as the update result of the filter coefficient vector.
  • k denotes an index of the number of times the filter coefficient vector h(i) is updated.
  • Step 330 Calculate N filter coefficient state reference quantities corresponding to the N pairs of synchronization frames according to the N filter coefficient vector mean values.
  • N filter coefficient state reference quantities corresponding to N pairs of sync frames can be obtained by the following two methods, but not limited to the following two methods.
  • the filter coefficient distribution map corresponding to the current filter coefficient vector mean value is determined according to a preset layout rule.
  • the elements in the current filter coefficient vector mean are at least divided into the main reflection area element, the left side element of the main reflection area, and the right side element of the main reflection area, as shown in FIG. 4A.
  • 4A is a diagram showing a distribution of filter coefficients in an embodiment of the present invention.
  • the number of elements in the filter coefficient vector is set by a specific adaptive filter algorithm, and the elements in the filter vector mean are sequentially extracted, and are sequentially arranged from left to right, and each two elements are in the horizontal direction.
  • the spacing of the projections is equal, the window length of the sliding window can be determined by the empirical value, and the number of elements contained in the sliding window is a fixed value, and the point A is the element with the largest absolute value in the filter coefficient vector, and the absolute value is the largest element.
  • the absolute value is represented by h, and the window height is at least twice the maximum absolute value of 2h.
  • the main reflection area element refers to an element included in the preset sliding window when the sum of the squares of all the elements included in the preset sliding window in the filter coefficient distribution map is the maximum value, and the preset sliding window in the filter coefficient distribution map
  • the element on the left is the element on the left side of the main reflection area
  • the element on the right side of the preset sliding window is the main reflection. The right side of the shot.
  • the sum of the squares of the elements on the right side of all the main reflection areas is calculated as the first energy value corresponding to the current filter coefficient vector mean.
  • the third preset algorithm is used to determine the filter coefficient state reference amount corresponding to the current filter coefficient vector mean.
  • the filter coefficient distribution map corresponding to the current filter coefficient vector mean value is determined according to a preset layout rule.
  • the elements in the current filter coefficient vector mean are at least divided into the main reflection area element, the left side element of the main reflection area, and the right side element of the main reflection area, as shown in FIG. 4A.
  • the main reflection area element refers to an element included in the preset sliding window when the sum of the squares of all the elements included in the preset sliding window in the filter coefficient distribution map is the maximum value, and the preset sliding window in the filter coefficient distribution map
  • the element on the left is the element on the left side of the main reflection area
  • the element on the right side of the preset sliding window is the element on the right side of the main reflection area.
  • the third preset algorithm is used to determine the filter coefficient state reference amount corresponding to the current filter coefficient vector mean.
  • the above two methods need to divide at least the elements in the current filter coefficient vector mean into the main reflection area element, the left side element of the main reflection area, and the main reflection area right according to the filter coefficient distribution map. Side element.
  • an element having the largest absolute value among the current filter coefficient vector mean values is determined. Construct a sliding window to determine the window height of the sliding window based on the absolute value of the element with the largest absolute value, with the preset value being the window length of the sliding window.
  • the filter coefficient distribution map corresponding to the filter coefficient vector mean value is determined, as shown in FIG. 4A, the elements in the filter coefficient vector.
  • the number is set by a specific adaptive filter algorithm, and the elements in the filter vector are sequentially taken out and arranged from left to right. The projection of each element in the horizontal direction is equal, and all the elements are connected. Get the filter coefficient distribution map.
  • the element is A point
  • the absolute value corresponding to point A is h
  • the window height is 2h.
  • the preset window is the window length of the sliding window, and the window length is smaller than the filter in the filter coefficient distribution map.
  • the length of the coefficient distribution can be arbitrarily set, but it is generally small, and can also be determined according to empirical values, such as one tenth of the length of the filter coefficient distribution, and the number of elements included in the sliding window is a fixed value.
  • FIG. 4A is a simplified description of the filter coefficient distribution diagram of the present invention.
  • the waveform drawn in FIG. 4A is relatively simple, and the number of elements included in the filter coefficient vector mean is small.
  • the filter coefficient distribution map As shown in Figure 4B.
  • Sliding the sliding window from the first element of the i-th filter coefficient vector mean to the last element on the filter coefficient profile, and determining that the sum of the squares of all the elements included in the sliding window is the maximum value, within the sliding window
  • the included element is the main reflection area element
  • the element on the left side of the sliding window is the left element of the main reflection area
  • the element on the right side of the sliding window is the right element of the main reflection area.
  • the elements contained in the sliding window are the main reflection area elements, and the elements on the left side of the sliding window are the left side elements of the main reflection area.
  • the absolute value of the main reflection area element is large, the main reflection area occupies most of the energy of the filter coefficient vector, the element on the left side of the sliding window is the left element of the main reflection area, and the absolute value of the element on the left side of the main reflection area is small.
  • the energy of the left side of the main reflection area is small, the element on the right side of the sliding window is the element on the right side of the main reflection area, the absolute value of the element on the right side of the main reflection area is larger than the absolute value of the element on the left side of the main reflection area, and the right side of the main reflection area Greater energy.
  • the sum of the squares of the elements on the right side of all the main reflection areas is calculated as the first energy value corresponding to the mean value of the i-th filter coefficient vector, as follows:
  • the first method is used to calculate the filter coefficient state reference amount corresponding to the i-th filter coefficient vector mean according to the first energy value and the smoothed value of the second energy value:
  • another method may be used to calculate a filter coefficient state reference quantity corresponding to the i-th filter coefficient vector mean value, specifically:
  • the second method is used to calculate the filter coefficient state reference amount corresponding to the i-th filter coefficient vector mean according to the first energy value and the smoothed value of the second energy value:
  • another method may be used to calculate a filter coefficient state reference quantity corresponding to the i-th filter coefficient vector mean value, specifically:
  • Step 340 Determine that the filter state diverges when at least M filter coefficient state reference quantities are greater than a preset threshold value among the N filter coefficient state reference quantities.
  • N and M are both positive integers, and N ⁇ M.
  • the update of the filter coefficient vector is stopped, and the current filter coefficient vector is cleared.
  • the method also includes comparing the filter coefficient state reference quantity and the preset threshold value in a certain time window to determine whether the filter state is divergent.
  • the ratio of the filter coefficient state reference quantity corresponding to the preset threshold value corresponding to 10 pairs of synchronization frames is taken. More results:
  • the ST(m) includes a comparison result between the filter coefficient state reference quantity corresponding to the 10 pairs of synchronization frames and the preset threshold value, and ST(i) represents the filter coefficient state reference quantity corresponding to the i-th synchronization frame and the preset. Comparison of threshold values.
  • a strategy Strategy(m) for determining the state of the filter divergence is set, that is, when there are 6 or more filter coefficient state reference quantities greater than a preset threshold value , to determine the filter state divergence, otherwise the filter state converges.
  • the filter state is determined to be diverged, and the filter coefficient vector update is determined according to the determination.
  • the embodiment of the invention provides a method for determining the state of the filter divergence, which is: collecting a near-end signal and a far-end signal synchronized with the near-end signal, and performing synchronous framing processing on the near-end signal and the far-end signal, respectively.
  • N filter coefficient state reference quantities corresponding to the N pairs of synchronization frames are respectively calculated, since the filter coefficient state reference quantity is insensitive to changes in signal energy, and is not affected by the near-end signal and The influence of the environment where the far-end signal is located, so the obtained result of the method is stable, and the calculation complexity of the filter coefficient state reference quantity is low, and no additional operation is needed to smooth the fluctuation of the signal. Finally, when at least M filter coefficient state reference quantities are greater than a preset threshold value among the N filter coefficient state reference quantities, it is determined that the filter state diverges.
  • the embodiment of the present invention uses the filter coefficient state reference quantity to determine whether the filter state is divergent, so as to achieve the purpose of controlling the echo canceler filter coefficient vector update and guaranteeing the call quality, and determining whether the filter state is diverged by the method of the present invention. It can effectively judge whether there is double talk during the call, so as to achieve the purpose of controlling the filter coefficient vector update and obtaining high quality output signal.
  • FIG. 5 is a schematic structural diagram of an apparatus for determining a state of a filter to be diverged according to an embodiment of the present invention.
  • the acquiring module 501 is configured to collect a near-end signal and a far-end signal synchronized with the near-end signal;
  • the framing module 502 is configured to perform synchronous framing processing on the near-end signal and the far-end signal collected by the collecting module 501, respectively, to obtain an N-pair synchronization frame, where each pair of synchronization frames includes a near-end signal frame and a near-end signal.
  • the filtering module 503 is configured to filter the N pairs of synchronization frames obtained by the framing module 502 by using a filter to obtain N filter coefficient vector mean values corresponding to the N pairs of synchronization frames;
  • the determining module 504 is configured to calculate N filter coefficient state reference quantities corresponding to the N pairs of synchronization frames according to the N filter coefficient vector average values obtained by the filtering module 503;
  • N and M are both positive integers, and N ⁇ M.
  • the framing module 502 is specifically configured to:
  • the near-end signals collected by the acquisition module 501 are subjected to frame processing in sequence according to the signal sequence to obtain N near-end signal frames;
  • the far-end signals collected by the acquisition module 501 and synchronized with the near-end signals are subjected to frame processing in sequence according to the signal sequence, and N far-end signal frames synchronized with the near-end signal frame are obtained;
  • Any one of the obtained near-end signal frame and the far-end signal frame synchronized with any one of the near-end signal frames constitutes a pair of synchronization frames, and a total of N pairs of synchronization frames are obtained.
  • the filtering module 503 is specifically configured to:
  • the filter coefficient vector mean values corresponding to the current pair of synchronization frames are obtained according to the filter coefficient vector update order T and the T-th filter coefficient vector corresponding to the current pair of synchronization frames.
  • the determining module 504 is specifically configured to:
  • the filter coefficient distribution map corresponding to the current filter coefficient vector mean value is determined according to a preset layout rule
  • the elements in the current filter coefficient vector mean are at least divided into a main reflection area element, a left side element of the main reflection area, and a right side element of the main reflection area, wherein the main reflection area element refers to a filter coefficient
  • the main reflection area element refers to a filter coefficient
  • the element on the right side of the preset sliding window is the element on the right side of the main reflection area;
  • the determining module 504 is further configured to:
  • the filter coefficient distribution map corresponding to the current filter coefficient vector mean value is determined according to a preset layout rule
  • the elements in the current filter coefficient vector mean are at least divided into The main reflection area element, the left side element of the main reflection area, and the right side element of the main reflection area, wherein the main reflection area element refers to the sum of the squares of all the elements included in the preset sliding window in the filter coefficient distribution map.
  • the main reflection area element refers to the sum of the squares of all the elements included in the preset sliding window in the filter coefficient distribution map.
  • the determining module 504 is further configured to:
  • the device further includes:
  • the reset module 505 is configured to stop updating the filter coefficient vector after the determination module determines that the filter state is diverged, and clear the current filter coefficient vector.
  • each functional module in each embodiment of the present application may be used. Can be integrated in one processing module or individual modules There are two or more modules that can be integrated into one module.
  • the above integrated modules can be implemented in the form of hardware or in the form of software functional modules.
  • the integrated modules if implemented in the form of software functional modules and sold or used as separate products, may be stored in a computer readable storage medium.
  • a computer readable storage medium A number of instructions are included to cause a computer device (which may be a personal computer, server, or network device, etc.) or a processor to perform all or part of the steps of the methods described in various embodiments of the present application.
  • the foregoing storage medium includes: a U disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk, and the like. .
  • FIG. 6 is a schematic structural diagram of a device for determining a state of a filter to be diverged according to an embodiment of the present invention, where the device includes a transceiver 601 and a processor. 602. Memory 603.
  • the transceiver 601, the processor 602, and the memory 603 are connected to each other.
  • the specific connecting medium between the above components is not limited in the embodiment of the present invention.
  • the memory 603, the processor 602, and the transceiver 601 are connected by a bus 604 in FIG. 6.
  • the bus is indicated by a thick line in FIG.
  • the bus can be divided into an address bus, a data bus, a control bus, and the like. For ease of representation, only one thick line is shown in Figure 6, but it does not mean that there is only one bus or one type of bus.
  • the transceiver 601 is configured to collect a near-end signal and a far-end signal synchronized with the near-end signal.
  • the processor 602 is configured to perform synchronous framing processing on the near-end signal and the far-end signal collected by the transceiver 601 to obtain N pairs of synchronization frames, where each pair of synchronization frames includes a near-end signal frame and a near-end signal frame.
  • the memory 603 is used to store the program code executed by the processor 602, and may be a volatile memory, such as a random access memory (English: random-access memory, abbreviation: RAM);
  • the memory 603 may also be a non-volatile memory (English: non-volatile memory), such as read-only memory (English: read-only memory, abbreviation: ROM), flash memory (English: flash memory), hard disk (English: hard Disk drive, abbreviated as: HDD) or solid state drive (English: solid-state drive, SSD), or memory 603 can be used to carry or store desired program code in the form of an instruction or data structure and can be accessed by a computer. Any other medium, but not limited to this.
  • the memory 603 may be a combination of the above memories.
  • the processor 602 in the embodiment of the present invention may be a central processing unit (English: central processing unit, CPU for short).
  • the embodiment of the present invention uses the filter coefficient state reference quantity to determine whether the filter state is diverged, so as to achieve the purpose of controlling the echo canceler filter coefficient vector update and guaranteeing the call quality, because the filter coefficient state reference quantity pair
  • the change of signal energy is not sensitive and is not affected by the environment of the near-end signal and the far-end signal. Therefore, the obtained result of the method has good stability, and the calculation complexity of the filter coefficient state reference quantity is low, therefore, Whether the filter state is diverged by the method of the present invention can effectively determine whether there is double talk during the call, thereby achieving the purpose of controlling the filter coefficient vector update and obtaining a high quality output signal.
  • embodiments of the present invention can be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment, or a combination of software and hardware. Moreover, the invention can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) including computer usable program code.
  • computer-usable storage media including but not limited to disk storage, CD-ROM, optical storage, etc.
  • the present invention is directed to a method, apparatus (system), and computer program according to an embodiment of the present invention.
  • the flow chart and/or block diagram of the product is described. It will be understood that each flow and/or block of the flowchart illustrations and/or FIG.
  • These computer program instructions can be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing device to produce a machine for the execution of instructions for execution by a processor of a computer or other programmable data processing device.
  • the computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture comprising the instruction device.
  • the apparatus implements the functions specified in one or more blocks of a flow or a flow and/or block diagram of the flowchart.
  • These computer program instructions can also be loaded onto a computer or other programmable data processing device such that a series of operational steps are performed on a computer or other programmable device to produce computer-implemented processing for execution on a computer or other programmable device.
  • the instructions provide steps for implementing the functions specified in one or more of the flow or in a block or blocks of a flow diagram.

Landscapes

  • Filters That Use Time-Delay Elements (AREA)

Abstract

L'invention concerne un procédé et un dispositif de détermination de divergence d'état de filtre, ceux-ci étant utilisés pour résoudre les problèmes selon lesquels, dans l'état de la technique, un procédé de détermination d'état de filtre présente une relativement grande complexité de calcul et une possibilité d'erreur de détermination. Le procédé consiste : à acquérir un signal d'extrémité proche et un signal d'extrémité distante, synchronisé avec le signal d'extrémité proche (300), et à effectuer respectivement un traitement de verrouillage de trame synchrone sur le signal d'extrémité proche et le signal d'extrémité distante, de manière à obtenir N paires de trames synchrones (310) ; à filtrer les N paires de trames synchrones par l'intermédiaire d'un filtre, de manière à obtenir une valeur moyenne de N vecteurs de coefficients de filtre correspondant aux N paires de trames synchrones (320) ; à calculer respectivement N quantités de référence d'état de coefficient de filtre correspondant aux N paires de trames synchrones en fonction de la valeur moyenne des N vecteurs de coefficients de filtre (330) ; lorsqu'il existe au moins M quantités de référence d'état de coefficient de filtre supérieures à une valeur seuil prédéfinie parmi les N quantités de référence d'état de coefficient de filtre, à déterminer que le filtre est dans un état divergent (340). Par conséquent, il est possible de déterminer efficacement si un filtre est dans un état divergent et de commander la mise à jour d'un vecteur de coefficients de filtre.
PCT/CN2016/075804 2015-07-21 2016-03-07 Procédé et dispositif de détermination de divergence d'état de filtre WO2017012350A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201510433689.1 2015-07-21
CN201510433689.1A CN105141287A (zh) 2015-07-21 2015-07-21 一种判断滤波器状态发散的方法及装置

Publications (1)

Publication Number Publication Date
WO2017012350A1 true WO2017012350A1 (fr) 2017-01-26

Family

ID=54726524

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2016/075804 WO2017012350A1 (fr) 2015-07-21 2016-03-07 Procédé et dispositif de détermination de divergence d'état de filtre

Country Status (2)

Country Link
CN (1) CN105141287A (fr)
WO (1) WO2017012350A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10596720B2 (en) 2015-09-24 2020-03-24 SWISS KRONO Tec AG Method for reducing the emission of volatile organic compounds made of wood materials, and wood materials

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105141287A (zh) * 2015-07-21 2015-12-09 华为技术有限公司 一种判断滤波器状态发散的方法及装置
CN110364174A (zh) * 2019-07-18 2019-10-22 申瓯通信设备有限公司 一种回声消除系统及回音消除方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1242651A (zh) * 1998-03-30 2000-01-26 摩托罗拉公司 用于具有解码器的装置的回声消除器及其方法
JP2003324370A (ja) * 2002-04-26 2003-11-14 Matsushita Electric Works Ltd エコーキャンセラ
CN101026392A (zh) * 2006-02-22 2007-08-29 阿尔卡特朗讯 对滤波器的调整进行控制的方法
CN101119135A (zh) * 2007-07-04 2008-02-06 深圳市融合视讯科技有限公司 一种用于回音消除的步长参数调整方法及设备
CN102377453A (zh) * 2010-08-06 2012-03-14 联芯科技有限公司 一种控制自适应滤波器更新的方法、装置及回声消除器
US8199927B1 (en) * 2007-10-31 2012-06-12 ClearOnce Communications, Inc. Conferencing system implementing echo cancellation and push-to-talk microphone detection using two-stage frequency filter
CN103259563A (zh) * 2012-02-16 2013-08-21 联芯科技有限公司 自适应滤波器发散检测方法及回声消除系统
CN105141287A (zh) * 2015-07-21 2015-12-09 华为技术有限公司 一种判断滤波器状态发散的方法及装置

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1242651A (zh) * 1998-03-30 2000-01-26 摩托罗拉公司 用于具有解码器的装置的回声消除器及其方法
JP2003324370A (ja) * 2002-04-26 2003-11-14 Matsushita Electric Works Ltd エコーキャンセラ
CN101026392A (zh) * 2006-02-22 2007-08-29 阿尔卡特朗讯 对滤波器的调整进行控制的方法
CN101119135A (zh) * 2007-07-04 2008-02-06 深圳市融合视讯科技有限公司 一种用于回音消除的步长参数调整方法及设备
US8199927B1 (en) * 2007-10-31 2012-06-12 ClearOnce Communications, Inc. Conferencing system implementing echo cancellation and push-to-talk microphone detection using two-stage frequency filter
CN102377453A (zh) * 2010-08-06 2012-03-14 联芯科技有限公司 一种控制自适应滤波器更新的方法、装置及回声消除器
CN103259563A (zh) * 2012-02-16 2013-08-21 联芯科技有限公司 自适应滤波器发散检测方法及回声消除系统
CN105141287A (zh) * 2015-07-21 2015-12-09 华为技术有限公司 一种判断滤波器状态发散的方法及装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10596720B2 (en) 2015-09-24 2020-03-24 SWISS KRONO Tec AG Method for reducing the emission of volatile organic compounds made of wood materials, and wood materials

Also Published As

Publication number Publication date
CN105141287A (zh) 2015-12-09

Similar Documents

Publication Publication Date Title
US10205830B2 (en) Echo cancellation data synchronization control method, terminal, and storage medium
Enzner et al. Acoustic echo control
US10250740B2 (en) Echo path change detector
WO2018188282A1 (fr) Procédé et dispositif d'annulation d'écho, ordinateur tablette de conférence et support d'informations informatique
EP2785032B1 (fr) Annulateur d'écho et procédé d'annulation d'écho
CN106486135B (zh) 近端语音检测器、语音系统、对语音进行分类的方法
WO2020181766A1 (fr) Procédé et dispositif de traitement de signal vocal, appareil et support de stockage lisible
US8644522B2 (en) Method and system for modeling external volume changes within an acoustic echo canceller
CN103391381A (zh) 回声消除方法及装置
WO2010083641A1 (fr) Procédé et appareil de détection de double parole
US8199907B2 (en) Echo canceller and echo cancelling method and program
van Waterschoot et al. Double-talk-robust prediction error identification algorithms for acoustic echo cancellation
WO2020252629A1 (fr) Procédé de détection d'écho acoustique résiduel, dispositif de détection d'écho acoustique résiduel, puce de traitement vocal et dispositif électronique
WO2017012350A1 (fr) Procédé et dispositif de détermination de divergence d'état de filtre
CN109559756B (zh) 滤波系数确定方法、回声消除方法、相应装置及设备
CN112489670A (zh) 时延估计方法、装置、终端设备和计算机可读存储介质
CN107786755B (zh) 一种双端通话检测方法和装置
Pfeifenberger et al. Acoustic Echo Cancellation with Cross-Domain Learning.
TR201815047T4 (tr) Bir uzak uç konuşmacı sinyali ve bir birleşik sinyal arasındaki akustik bir bağlamanın belirlenmesi.
CN111355855B (zh) 回声处理方法、装置、设备及存储介质
WO2022194011A1 (fr) Procédé et appareil d'annulation d'écho, dispositif et support de stockage
CN110148421A (zh) 一种残余回声检测方法、终端和装置
CN111989934B (zh) 回声消除装置、回声消除方法、信号处理芯片及电子设备
CN111883155B (zh) 回声消除方法、装置及存储介质
Szwoch et al. A low complexity double-talk detector based on the signal envelope

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16827041

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 16827041

Country of ref document: EP

Kind code of ref document: A1