CN106878575B - The estimation method and device of residual echo - Google Patents
The estimation method and device of residual echo Download PDFInfo
- Publication number
- CN106878575B CN106878575B CN201710105912.9A CN201710105912A CN106878575B CN 106878575 B CN106878575 B CN 106878575B CN 201710105912 A CN201710105912 A CN 201710105912A CN 106878575 B CN106878575 B CN 106878575B
- Authority
- CN
- China
- Prior art keywords
- remote signaling
- power spectrum
- residual signals
- residual
- alignment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M9/00—Arrangements for interconnection not involving centralised switching
- H04M9/08—Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04M—TELEPHONIC COMMUNICATION
- H04M3/00—Automatic or semi-automatic exchanges
- H04M3/22—Arrangements for supervision, monitoring or testing
Abstract
The estimation method and device of residual echo provided in an embodiment of the present invention, belong to speech signal processing technology.The described method includes: obtaining remote signaling after residual signals and alignment；Obtain the first power spectrum of remote signaling and the residual signals and the coherent function for being aligned rear remote signaling after the power spectrum of the residual signals, the alignment；Adaptive smooth parameter is calculated；Based on the adaptive smooth parameter and preset rules, the second power spectrum of remote signaling after the alignment is obtained；According to the second power spectrum of remote signaling after the coherent function and the alignment, the estimated value of residual echo is obtained.Residual echo can be accurately estimated with this.
Description
Technical field
The present invention relates to speech signal processing technologies, a kind of estimation method in particular to residual echo and
Device.
Background technique
In reality in loudspeaker and microphone circuit system, due to the non-linear frequency response of power amplifier, loudspeaker vibration
The factors such as non-linear and microphone non-linear frequency response, cause linear echo elimination to be difficult echo cancellor is clean, residual
Echo has pace of change fast, the big feature of energy dynamics range.The residual echo suppression module of existing echo cancellation technology is right
The estimation of residual echo causes residual echo to inhibit insufficient, seriously affects sense of hearing generally there are underestimating.Meeting in engineering
Empirical inhibition is carried out to the residual echo underestimated, although this method can eliminate echo, the language Chong Die with echo
Line point also will receive inhibition, even be eliminated so as to cause voice distortion.
Summary of the invention
In view of this, the estimation method and device for being designed to provide a kind of residual echo of the embodiment of the present invention, to change
The kind above problem.To achieve the goals above, the technical solution adopted by the present invention is as follows:
In a first aspect, the embodiment of the invention provides a kind of estimation methods of residual echo, which comprises obtain residual
Remote signaling after difference signal and alignment；Obtain the first power of remote signaling after the power spectrum of the residual signals, the alignment
Spectrum and the residual signals and the coherent function for being aligned rear remote signaling；Adaptive smooth parameter is calculated；Based on institute
Adaptive smooth parameter and preset rules are stated, the second power spectrum of remote signaling after the alignment is obtained；According to described relevant
The second power spectrum of remote signaling, obtains the estimated value of residual echo after function and the alignment.
Second aspect, the embodiment of the invention provides a kind of estimation device of residual echo, first acquisition units, for obtaining
Take remote signaling after obtaining residual signals and alignment；Second acquisition unit, for obtaining the power spectrum, described of the residual signals
First power spectrum of remote signaling and the residual signals and the coherent function for being aligned rear remote signaling after alignment；It calculates single
Member, for adaptive smooth parameter to be calculated；Third acquiring unit, for being based on the adaptive smooth parameter and presetting
Rule obtains the second power spectrum of remote signaling after the alignment；Estimated value obtaining unit, for according to the coherent function with
And after the alignment remote signaling the second power spectrum, obtain the estimated value of residual echo.
The estimation method and device of a kind of residual echo provided in an embodiment of the present invention, which comprises obtain residual error
Remote signaling after signal and alignment；Obtain the first power spectrum of remote signaling after the power spectrum of the residual signals, the alignment
With the residual signals and the coherent function for being aligned rear remote signaling；Adaptive smooth parameter is calculated；Based on described
Adaptive smooth parameter and preset rules obtain the second power spectrum of remote signaling after the alignment；According to the relevant letter
The second power spectrum of remote signaling, obtains the estimated value of residual echo after the several and described alignment.Can accurately it be estimated with this
Residual echo.
Other features and advantages of the present invention will be illustrated in subsequent specification, also, partly be become from specification
It is clear that by implementing understanding of the embodiment of the present invention.The objectives and other advantages of the invention can be by written theory
Specifically noted structure is achieved and obtained in bright book, claims and attached drawing.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 is a kind of structural block diagram that can be applied to the electronic equipment in the embodiment of the present application；
Fig. 2 is the flow chart of the estimation method for the residual echo that first embodiment of the invention provides；
Fig. 3 is the application environment schematic diagram of the estimation method for the residual echo that first embodiment of the invention provides；
Fig. 4 is the structural block diagram of the estimation device for the residual echo that second embodiment of the invention provides.
In figure: 100- electronic equipment；110- memory；111- storage control；112- processor；300- loudspeaker；
302- microphone；304- delay estimation module；306- alignment module；308- acoustic echo cancellation module；The suppression of 310- residual echo
Molding block；400- device；410- first acquisition unit；420- second acquisition unit；430- computing unit；440- third obtains single
Member；450- estimated value obtaining unit；460- a priori SNR estimation unit；470- gain calculating unit；480- output obtains single
Member.
Specific embodiment
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause
This, is not intended to limit claimed invention to the detailed description of the embodiment of the present invention provided in the accompanying drawings below
Range, but it is merely representative of selected embodiment of the invention.Based on the embodiment of the present invention, those skilled in the art are not doing
Every other embodiment obtained under the premise of creative work out, shall fall within the protection scope of the present invention.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing.Meanwhile of the invention
In description, term " first ", " second " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Referring to Fig. 1, Fig. 1 shows a kind of structural block diagram of electronic equipment 100 that can be applied in the embodiment of the present application.
The electronic equipment 100 can be used as user terminal, be also possible to computer or server, the user terminal can for mobile phone or
Tablet computer.As shown in Figure 1, electronic equipment 100 may include memory 110, storage control 111, processor 112 and residual
The estimation device of echo.
Memory 110, storage control 111 are directly or indirectly electrically connected between each element of processor 112, to realize
The transmission or interaction of data.For example, can realize electricity by one or more communication bus or signal bus between these elements
Connection.The estimation method of residual echo respectively includes at least one can be stored in the form of software or firmware (firmware)
Software function module in memory 110, such as the estimation device software function module that includes or the calculating of the residual echo
Machine program.
Memory 110 can store various software programs and module, such as residual echo provided by the embodiments of the present application
Estimation method and the corresponding program instruction/module of device.Software journey of the processor 112 by operation storage in the memory 110
Sequence and module realize the residual echo in the embodiment of the present application thereby executing various function application and data processing
Estimation method.Memory 110 can include but is not limited to random access memory (Random Access Memory, RAM), only
It reads memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only
Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM),
Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Processor 112 can be a kind of IC chip, have signal handling capacity.Above-mentioned processor can be general
Processor, including central processing unit (Central Processing Unit, abbreviation CPU), network processing unit (Network
Processor, abbreviation NP) etc.；It can also be digital signal processor (DSP), specific integrated circuit (ASIC), ready-made programmable
Gate array (FPGA) either other programmable logic device, discrete gate or transistor logic, discrete hardware components.It can
To realize or execute disclosed each method, step and the logic diagram in the embodiment of the present application.General processor can be micro-
Processor or the processor are also possible to any conventional processor etc..
First embodiment
Referring to figure 2., the embodiment of the invention provides a kind of estimation methods of residual echo, which comprises
Step S200: remote signaling after residual signals and alignment is obtained；
Step S210: the first power spectrum of remote signaling and institute after the power spectrum of the residual signals, the alignment are obtained
State residual signals and the coherent function for being aligned rear remote signaling；
In reality in loudspeaker and microphone circuit system, due to the non-linear frequency response of power amplifier, loudspeaker vibration
The factors such as non-linear and microphone non-linear frequency response, cause linear echo elimination to be difficult echo cancellor is clean, residual
Echo inhibits for eliminating residual echo.
Residual echo has pace of change fast, the big feature of energy dynamics range.The residual of existing echo cancellation technology is returned
Sound suppression module causes residual echo to inhibit insufficient, seriously affects to the estimation of residual echo generally there are underestimating
Sense of hearing.
Empirical inhibition can be carried out to the residual echo underestimated in engineering, although this method can eliminate echo,
It is that the phonological component Chong Die with echo also will receive inhibition, is even eliminated so as to cause voice distortion.
Referring to figure 3., as shown in figure 3, loudspeaker 300 plays reference signal x (n), the echo generated in spatial is believed
Number d (n)；Microphone 302 acquires main signal p (n), p (n)=s (n)+d (n)+b (n), and wherein x (n) is reference signal, and s (n) is
Voice signal, b (n) are noise signal, and p (n) is main signal.N be discrete time scale, n=0,1,2 ....Delay estimation module
304 with reference signal x (n) and main signal p (n) for input signal, export as Delay Estima-tion value K；xhistoryIt (n) is reference signal
The historical data of x (n).Alignment module 306 is according to Delay Estima-tion value K, in xhistory(n) distally believe after finding time unifying in
Number i.e. x (n-K).Acoustic echo cancellation module 308 is based on remote signaling x (n-K) after alignment, obtains estimated echo signal
Further according to the main signal of input, residual signals are obtainedResidual signals x (n-K) and alignment after remote signaling
Input signal of the x (n-K) as the residual echo suppression module 310.Sig (n) is the residual echo suppression module 310
Output signal.
Specifically, the residual signals x (n-K) and alignment after remote signaling x (n-K) using as unit of frame as described residual
Stay the input signal of echo suppression module 310.It needs to input speech signal sub-frame processing.
Based on following formula (1):
In formula (1), e is residual signals, and m is frame number, and n is discrete time, and ω is discrete frequency, and win is adding window, and FFT is
Fast discrete Fourier transformation, LfftFor the length of fast discrete Fourier transformation, LshiftIndicate that frame moves, K is Delay Estima-tion value, x
For remote signaling after alignment, SeeFor the power spectrum of the residual signals, SxxFor the first power of remote signaling after the alignment
Spectrum, SexTo indicate the residual signals and the crosspower spectrum for being aligned rear remote signaling, ρ is constant, CohexFor the residual error
Signal and the coherent function for being aligned rear remote signaling.
Step S220: adaptive smooth parameter is calculated；
Step S230: being based on the adaptive smooth parameter and preset rules, obtains remote signaling after the alignment
Second power spectrum；
Based on following formula (2):
Sxx,new(m, ω)=(1- σ (m, ω)) Sxx,new(m-1,ω)+σ(m,ω)XH(m-K/Lshift,ω)·X(m-K/
Lshift,ω) (2)
Second power spectrum of remote signaling after being aligned, wherein formula (2) is the preset rules, Sxx,nweFor alignment
Second power spectrum of remote signaling afterwards, σ (m, ω) are the adaptive smooth parameter；
σ (m, ω)=min (XH(m-K,ω)·X(m-K,ω)/Sxx(m,ω),1)。
Step S240: it according to the second power spectrum of remote signaling after the coherent function and the alignment, is remained
The estimated value of echo.
Based on Szz(m, ω)=abs (Cohex(m,ω))2·Sxx,new(m, ω) obtains the estimated value of residual echo,
In, SzzFor the estimated value of residual echo, CohexFor the coherent function, Sxx,newFor the second function of remote signaling after the alignment
Rate spectrum.
In order to obtain the output signal for eliminating residual echo, the method also includes:
Step S250: the power spectrum of estimated value and the residual signals based on the residual echo obtains priori letter
It makes an uproar the estimated value of ratio；
Based on following formula (3):
Wherein SzzFor the estimated value of residual echo, β (m, ω)=(ξDD(m,ω)/(1+ξDD(m,ω)))2, η (m, ω) is
Posteriori SNR, ξDD(m, ω) is the prior weight estimated based on Dicision-Directed method, sigprev(m-1,ω)
It is exported for the frequency domain of former frame, θ is constant, and ξ (m, ω) is the estimated value of prior weight.So obtain estimating for prior weight
Evaluation ξ (m, ω).
Step S260: wiener gain is calculated in the estimated value based on the prior weight；
Based on formula (4):
Wherein, G (m, ω) is wiener gain.
Step S270: the Fourier transformation based on the wiener gain, the residual signals, by inverse fourier transform,
Obtain output signal.
Based on formula (5):
In formula (5), E (m, ω) is the Fourier transformation of residual signals in formula (1), and IFFT () is fast Fourier
Inverse transformation, sig (m, n) are the output signal of m frame.
The estimation method of residual echo provided in an embodiment of the present invention is calculated based on adaptive smooth parameter, Ke Yiyou
Effect prevents from underestimating residual echo, and can preferably track the energy spectrum of residual echo.
The estimation method of residual echo provided in an embodiment of the present invention, which comprises obtain residual signals and alignment
Remote signaling afterwards；Obtain the first power spectrum of remote signaling and the residual error after the power spectrum of the residual signals, the alignment
Signal and the coherent function for being aligned rear remote signaling；Adaptive smooth parameter is calculated；Based on the adaptive smooth
Parameter and preset rules obtain the second power spectrum of remote signaling after the alignment；According to the coherent function and described
The second power spectrum of remote signaling, obtains the estimated value of residual echo after alignment.Residual echo can be accurately estimated with this.
Second embodiment
Referring to figure 4., the embodiment of the invention provides a kind of estimation device 400 of residual echo, described device 400 is wrapped
It includes:
First acquisition unit 410, for remote signaling after obtaining residual signals and being aligned.
Second acquisition unit 420, for obtaining first of remote signaling after the power spectrum of the residual signals, the alignment
Power spectrum and the residual signals and the coherent function for being aligned rear remote signaling.
Specifically, second acquisition unit 420 are used for:
Based on following formula (1):
Wherein, e is residual signals, and m is frame number, and n is discrete time, and ω is discrete frequency, and win is adding window, and FFT is quick
Discrete Fourier transform, LfftFor the length of fast discrete Fourier transformation, LshiftIndicate that frame moves, K is Delay Estima-tion value, and x is pair
Remote signaling after neat, SeeFor the power spectrum of the residual signals, SxxFor the first power spectrum of remote signaling after the alignment, Sex
To indicate the residual signals and the crosspower spectrum for being aligned rear remote signaling, ρ is constant, CohexFor the residual signals with
The coherent function of remote signaling after the alignment.
Computing unit 430, for adaptive smooth parameter to be calculated.
Third acquiring unit 440, for being based on the adaptive smooth parameter and preset rules, after obtaining the alignment
Second power spectrum of remote signaling；
Third acquiring unit 440, is used for:
Based on following formula (2):
Sxx,new(m, ω)=(1- σ (m, ω)) Sxx,new(m-1,ω)+σ(m,ω)XH(m-K/Lshift,ω)·X(m-K/
Lshift,ω) (2)
Second power spectrum of remote signaling after being aligned, wherein Sxx,nweFor the second power of remote signaling after alignment
Spectrum, σ (m, ω) are the adaptive smooth parameter；
σ (m, ω)=min (XH(m-K,ω)·X(m-K,ω)/Sxx(m,ω),1)。
Estimated value obtaining unit 450, for the second function according to remote signaling after the coherent function and the alignment
Rate spectrum, obtains the estimated value of residual echo.
Estimated value obtaining unit 450, for being based on Szz(m, ω)=abs (Cohex(m,ω))2·Sxx,new(m, ω), is obtained
Obtain the estimated value of residual echo, wherein SzzFor the estimated value of residual echo, Coh is that the residual signals are remote with after described be aligned
The coherent function of end signal, Sxx,newFor the second power spectrum of remote signaling after the alignment.
Described device 400 further include:
A priori SNR estimation unit 460, for based on the residual echo estimated value and the residual signals
Power spectrum obtains the estimated value of prior weight；
Wiener gain is calculated for the estimated value based on the prior weight in gain calculating unit 470；
It exports obtaining unit 480 and passes through Fu for the Fourier transformation based on the wiener gain, the residual signals
Vertical leaf inverse transformation, obtains output signal.
The above each unit can be by software code realization, at this point, above-mentioned each unit can be stored in memory 110.
The above each unit can equally be realized by hardware such as IC chip.
The technical effect of the estimation device 400 of residual echo provided in an embodiment of the present invention, realization principle and generation and
Preceding method embodiment is identical, and to briefly describe, Installation practice part does not refer to place, can refer in preceding method embodiment
Corresponding contents.
In several embodiments provided herein, it should be understood that disclosed device and method can also pass through
Other modes are realized.The apparatus embodiments described above are merely exemplary, for example, flow chart and block diagram in attached drawing
Show the device of multiple embodiments according to the present invention, the architectural framework in the cards of method and computer program product,
Function and operation.In this regard, each box in flowchart or block diagram can represent the one of a module, section or code
Part, a part of the module, section or code, which includes that one or more is for implementing the specified logical function, to be held
Row instruction.It should also be noted that function marked in the box can also be to be different from some implementations as replacement
The sequence marked in attached drawing occurs.For example, two continuous boxes can actually be basically executed in parallel, they are sometimes
It can execute in the opposite order, this depends on the function involved.It is also noted that every in block diagram and or flow chart
The combination of box in a box and block diagram and or flow chart can use the dedicated base for executing defined function or movement
It realizes, or can realize using a combination of dedicated hardware and computer instructions in the system of hardware.
In addition, each functional module in each embodiment of the present invention can integrate one independent portion of formation together
Point, it is also possible to modules individualism, an independent part can also be integrated to form with two or more modules.
It, can be with if the function is realized and when sold or used as an independent product in the form of software function module
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.It needs
Illustrate, herein, relational terms such as first and second and the like be used merely to by an entity or operation with
Another entity or operation distinguish, and without necessarily requiring or implying between these entities or operation, there are any this realities
The relationship or sequence on border.Moreover, the terms "include", "comprise" or its any other variant are intended to the packet of nonexcludability
Contain, so that the process, method, article or equipment for including a series of elements not only includes those elements, but also including
Other elements that are not explicitly listed, or further include for elements inherent to such a process, method, article, or device.
In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including the element
Process, method, article or equipment in there is also other identical elements.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.It should also be noted that similar label and letter exist
Similar terms are indicated in following attached drawing, therefore, once being defined in a certain Xiang Yi attached drawing, are then not required in subsequent attached drawing
It is further defined and explained.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
It should be noted that, in this document, relational terms such as first and second and the like are used merely to a reality
Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation
In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to
Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, article or equipment
Intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that
There is also other identical elements in process, method, article or equipment including the element.
Claims (4)
1. a kind of estimation method of residual echo, which is characterized in that the described method includes:
Remote signaling after obtaining residual signals and being aligned；
Based on formula
Obtain after the power spectrum of the residual signals, the alignment the first power spectrum of remote signaling and the residual signals with it is described right
The coherent function of remote signaling after neat；Wherein, e is residual signals, and m is frame number, and n is discrete time, and ω is discrete frequency, win
For adding window, FFT is fast discrete Fourier transformation, LfftFor the length of fast discrete Fourier transformation, LshiftIndicate that frame moves, K is
Delay Estima-tion value, x are remote signaling after alignment, SeeFor the power spectrum of the residual signals, SxxFor remote signaling after the alignment
The first power spectrum, SexTo indicate the residual signals and the crosspower spectrum for being aligned rear remote signaling, ρ is constant, Cohex
For the residual signals and the coherent function for being aligned rear remote signaling, X is the fast discrete Fourier for being aligned remote signaling
Transformed value；
Based on formula σ (m, ω)=min (XH(m-K,ω)·X(m-K,ω)/Sxx(m, ω), 1) adaptive smooth is calculated
Parameter, σ (m, ω) are the adaptive smooth parameter；
Based on formula Sxx,new(m, ω)=(1- σ (m, ω)) Sxx,new(m-1,ω)+σ(m,ω)XH(m-K/Lshift,ω)·X
(m-K/Lshift, ω) be aligned after remote signaling the second power spectrum, wherein SXx, newFor remote signaling after the alignment
Second power spectrum；
Based on Szz(m, ω)=abs (Cohex(m,ω))2·Sxx,new(m, ω) obtains the estimated value of residual echo, wherein Szz
For the estimated value of residual echo, CohexFor the coherent function, SXx, newFor the second power spectrum of remote signaling after the alignment.
2. the method according to claim 1, wherein the method also includes:
The power spectrum of estimated value and the residual signals based on the residual echo, obtains the estimated value of prior weight；
Based on the estimated value of the prior weight, wiener gain is calculated；
Fourier transformation based on the wiener gain, the residual signals obtains output signal by inverse fourier transform.
3. a kind of estimation device of residual echo, which is characterized in that described device includes:
First acquisition unit, for remote signaling after obtaining residual signals and being aligned；
Second acquisition unit, for being based on formula
Obtain after the power spectrum of the residual signals, the alignment the first power spectrum of remote signaling and the residual signals with it is described right
The coherent function of remote signaling after neat；Wherein, e is residual signals, and m is frame number, and n is discrete time, and ω is discrete frequency, win
For adding window, FFT is fast discrete Fourier transformation, LfftFor the length of fast discrete Fourier transformation, LshiftIndicate that frame moves, K is
Delay Estima-tion value, x are remote signaling after alignment, SeeFor the power spectrum of the residual signals, SxxFor remote signaling after the alignment
The first power spectrum, SexTo indicate the residual signals and the crosspower spectrum for being aligned rear remote signaling, ρ is constant, Cohex
For the residual signals and the coherent function for being aligned rear remote signaling, X is the fast discrete Fourier for being aligned remote signaling
Transformed value；Computing unit is used for formula σ (m, ω)=min (XH(m-K,ω)·X(m-K,ω)/Sxx(m, ω), 1) it calculates
To adaptive smooth parameter, σ (m, ω) is the adaptive smooth parameter；
Third acquiring unit is based on formula Sxx,new(m, ω)=(1- σ (m, ω)) Sxx,new(m-1,ω)+σ(m,ω)XH(m-K/
Lshift,ω)·X(m-K/Lshift, ω) be aligned after remote signaling the second power spectrum, wherein SXx, newFor the alignment
Second power spectrum of remote signaling afterwards；
Estimated value obtaining unit is based on Szz(m, ω)=abs (Cohex(m,ω))2·Sxx,new(m, ω), obtains residual echo
Estimated value, wherein SzzFor the estimated value of residual echo, CohexFor the coherent function, SXx, newDistally to believe after the alignment
Number the second power spectrum.
4. device according to claim 3, which is characterized in that described device further include:
A priori SNR estimation unit, for the power spectrum of estimated value and the residual signals based on the residual echo,
Obtain the estimated value of prior weight；
Wiener gain is calculated for the estimated value based on the prior weight in gain calculating unit；
It exports obtaining unit and passes through Fourier inversion for the Fourier transformation based on the wiener gain, the residual signals
It changes, obtains output signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710105912.9A CN106878575B (en) | 2017-02-24 | 2017-02-24 | The estimation method and device of residual echo |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710105912.9A CN106878575B (en) | 2017-02-24 | 2017-02-24 | The estimation method and device of residual echo |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106878575A CN106878575A (en) | 2017-06-20 |
| CN106878575B true CN106878575B (en) | 2019-11-05 |
Family
ID=59167889
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710105912.9A Active CN106878575B (en) | 2017-02-24 | 2017-02-24 | The estimation method and device of residual echo |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106878575B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107749303B (en) * | 2017-07-31 | 2020-03-31 | 电信科学技术研究院 | Post-processing method and device for voice signal output by acoustic echo canceller |
| CN108696648B (en) * | 2018-05-16 | 2021-08-24 | 上海小度技术有限公司 | Method, device, equipment and storage medium for processing short-time voice signal |
| CN109120379A (en) * | 2018-08-30 | 2019-01-01 | 武汉虹信通信技术有限责任公司 | A kind of adaptive modulation coding method suitable for the more scenes of wireless communication system |
| CN109547655A (en) * | 2018-12-30 | 2019-03-29 | 广东大仓机器人科技有限公司 | A kind of method of the echo cancellation process of voice-over-net call |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1176034A (en) * | 1995-02-24 | 1998-03-11 | 艾利森公司 | Apparatus and method for canceling acoustic echoes including non-linear distortions in loudspeaker telephones |
| CN101964670A (en) * | 2009-07-21 | 2011-02-02 | 雅马哈株式会社 | Echo suppression method and apparatus thereof |
| CN102196109A (en) * | 2010-03-01 | 2011-09-21 | 联芯科技有限公司 | Residual echo detection method and system |
| CN103327201A (en) * | 2012-03-20 | 2013-09-25 | 联芯科技有限公司 | Elimination method and system of remaining echoes |
| CN103561184A (en) * | 2013-11-05 | 2014-02-05 | 武汉烽火众智数字技术有限责任公司 | Frequency-convertible echo cancellation method based on near-end audio signal calibration and correction |
| CN104427144A (en) * | 2013-09-11 | 2015-03-18 | 联芯科技有限公司 | Linear echo cancellation method and device thereof |
| CN104427143A (en) * | 2013-09-06 | 2015-03-18 | 联芯科技有限公司 | Residual echo detection method and system |
| CN105280195A (en) * | 2015-11-04 | 2016-01-27 | 腾讯科技（深圳）有限公司 | Method and device for processing speech signal |
| WO2016098164A1 (en) * | 2014-12-15 | 2016-06-23 | 三菱電機株式会社 | Echo canceler device and echo cancel method |
| CN105825864A (en) * | 2016-05-19 | 2016-08-03 | 南京奇音石信息技术有限公司 | Double-talk detection and echo cancellation method based on zero-crossing rate |
-
2017
- 2017-02-24 CN CN201710105912.9A patent/CN106878575B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1176034A (en) * | 1995-02-24 | 1998-03-11 | 艾利森公司 | Apparatus and method for canceling acoustic echoes including non-linear distortions in loudspeaker telephones |
| CN101964670A (en) * | 2009-07-21 | 2011-02-02 | 雅马哈株式会社 | Echo suppression method and apparatus thereof |
| CN102196109A (en) * | 2010-03-01 | 2011-09-21 | 联芯科技有限公司 | Residual echo detection method and system |
| CN103327201A (en) * | 2012-03-20 | 2013-09-25 | 联芯科技有限公司 | Elimination method and system of remaining echoes |
| CN104427143A (en) * | 2013-09-06 | 2015-03-18 | 联芯科技有限公司 | Residual echo detection method and system |
| CN104427144A (en) * | 2013-09-11 | 2015-03-18 | 联芯科技有限公司 | Linear echo cancellation method and device thereof |
| CN103561184A (en) * | 2013-11-05 | 2014-02-05 | 武汉烽火众智数字技术有限责任公司 | Frequency-convertible echo cancellation method based on near-end audio signal calibration and correction |
| WO2016098164A1 (en) * | 2014-12-15 | 2016-06-23 | 三菱電機株式会社 | Echo canceler device and echo cancel method |
| CN105280195A (en) * | 2015-11-04 | 2016-01-27 | 腾讯科技（深圳）有限公司 | Method and device for processing speech signal |
| CN105825864A (en) * | 2016-05-19 | 2016-08-03 | 南京奇音石信息技术有限公司 | Double-talk detection and echo cancellation method based on zero-crossing rate |
Non-Patent Citations (1)
| Title |
|---|
| 用于回声消除系统的自适应延时估计算法研究;陈国志，乐彦杰，张磊，胡凯文，昊翔;《科学技术与工程》;20150131;244-249 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106878575A (en) | 2017-06-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106878575B (en) | The estimation method and device of residual echo | |
| Brandt et al. | Integrating time signals in frequency domain–Comparison with time domain integration | |
| CN104994249B (en) | Sound method for echo cancellation and device | |
| CN108293170B (en) | Method and apparatus for adaptive phase distortion free amplitude response equalization in beamforming applications | |
| CN109102821B (en) | Time delay estimation method, time delay estimation system, storage medium and electronic equipment | |
| Jukić et al. | Speech dereverberation using weighted prediction error with Laplacian model of the desired signal | |
| CN102576543A (en) | Multi-input noise suppresion device, multi-input noise suppression method, program, and integrated circuit | |
| CN112017679B (en) | Method, device and equipment for updating adaptive filter coefficients | |
| KR101568937B1 (en) | Apparatus and method for supressing non-linear echo talker using volterra filter | |
| CN112863535A (en) | Residual echo and noise elimination method and device | |
| US9978383B2 (en) | Method for processing speech/audio signal and apparatus | |
| EP2752848A1 (en) | Method and apparatus for generating a noise reduced audio signal using a microphone array | |
| CN105989850B (en) | A kind of echo cancellation method and device | |
| CN106031196A (en) | Signal-processing device, method, and program | |
| CN104200811A (en) | Self-adaption spectral subtraction and noise elimination processing method and device for voice signals | |
| US20130208905A1 (en) | Temporal Interpolation Of Adjacent Spectra | |
| CN111863008A (en) | Audio noise reduction method and device and storage medium | |
| Gil-Cacho et al. | A frequency-domain adaptive filter (FDAF) prediction error method (PEM) framework for double-talk-robust acoustic echo cancellation | |
| CN105144290A (en) | Signal processing device, signal processing method, and signal processing program | |
| EP2689419A1 (en) | Method and arrangement for damping dominant frequencies in an audio signal | |
| Hofmann et al. | Significance-aware filtering for nonlinear acoustic echo cancellation | |
| Lee et al. | Stereophonic acoustic echo suppression incorporating spectro-temporal correlations | |
| Sandoval et al. | Using linear prediction to mitigate end effects in empirical mode decomposition | |
| Zhao et al. | A collaborative spline adaptive filter for nonlinear echo cancellation | |
| CN106847299B (en) | Time delay estimation method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |