CN109102821A - Delay time estimation method, system, storage medium and electronic equipment - Google Patents
Delay time estimation method, system, storage medium and electronic equipment Download PDFInfo
- Publication number
- CN109102821A CN109102821A CN201811049712.7A CN201811049712A CN109102821A CN 109102821 A CN109102821 A CN 109102821A CN 201811049712 A CN201811049712 A CN 201811049712A CN 109102821 A CN109102821 A CN 109102821A
- Authority
- CN
- China
- Prior art keywords
- reference signal
- frequency domain
- signal
- time delay
- energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
- G10L2021/02082—Noise filtering the noise being echo, reverberation of the speech
Abstract
The present invention discloses a kind of delay time estimation method, comprising: obtains reference signal and microphone Mike's signal collected, and carries out Fast Fourier Transform (FFT) to obtain frequency domain reference signal and frequency domain Mike's signal；The frequency domain reference signal is input to sef-adapting filter, to obtain corresponding to the related frequency domain reference signal of the frequency domain reference signal included in the frequency domain Mike signal, the frequency domain Mike signal is for updating the sef-adapting filter；Filter energy is calculated for determining time delay value according to the related frequency domain reference signal of sef-adapting filter output.Solve in the prior art using cross-correlation method, performance sharply declines under environment interference situation, in the case that the more complicated scene of environment or it is double say the highly unstable problem of time delay estimated result.
Description
Technical field
The present invention relates to Audio Signal Processing technical field more particularly to a kind of delay time estimation method, system, storage mediums
And electronic equipment.
Background technique
With the continuous development of information technology, various distributed intelligence hardware are increasingly wide in the application of every field
It is general.Echo cancellor is as the hot spot that link indispensable in smart machine interaction is always that those skilled in the relevant arts study.
For capableing of the TV set-top box of voice control, since the collected signal of microphone can mix phonetic order and TV programme
Own voice needs to eliminate TV speech to distinguish the phonetic order of user's sending.Therefore, accurately estimation
TV speech signal largely affects the quality of echo cancellation to the time delay between microphone.
Echo cancellor is that distal end is prevented by eliminating or removing the audio signal for the distal end picked up in local microphone
A kind of processing method for returning to of sound.Existing typical echo cancel method is the method based on time delay estimation, calculates ginseng
The linear dependence of signal and Mike's signal is examined, and chooses the corresponding time delay of maximum cross-correlation as device time delay, to be used for back
What sound was eliminated.
Although cross-correlation method is a kind of delay time estimation method that principle is relatively simple, but since algorithm itself lacks
Fall into, performance sharply declines under environment interference situation, in the case that the more complicated scene of environment or it is double say time delay estimation knot
Fruit is highly unstable.In addition, the points of FFT are more, and complexity is bigger, CPU usage since the time delay for needing to estimate is longer
It is higher.
Summary of the invention
The embodiment of the present invention provides a kind of delay time estimation method, system electronics and storage medium, at least solving
One of above-mentioned technical problem.
In a first aspect, the embodiment of the present invention provides a kind of delay time estimation method, comprising:
Reference signal and microphone Mike's signal collected are obtained, and carries out Fast Fourier Transform (FFT) to obtain frequency domain ginseng
Examine signal and frequency domain Mike's signal；
The frequency domain reference signal is input to sef-adapting filter, to obtain included in the frequency domain Mike signal
Corresponding to the related frequency domain reference signal of the frequency domain reference signal, the frequency domain Mike signal is for updating the adaptive filter
Wave device；
Filter energy is calculated for true according to the related frequency domain reference signal of sef-adapting filter output
Fixed response time value.
Second aspect, the embodiment of the present invention provide a kind of time delay estimating system, comprising:
Signal acquisition program module for obtaining reference signal and microphone Mike's signal collected, and carries out quick
Fourier transformation is to obtain frequency domain reference signal and frequency domain Mike's signal；
Adaptive-filtering program module, for the frequency domain reference signal to be input to sef-adapting filter, to obtain
State the related frequency domain reference signal for corresponding to the frequency domain reference signal included in frequency domain Mike's signal, the frequency domain Mike
Signal is for updating the sef-adapting filter；
Time delay determines program module, based on the related frequency domain reference signal exported according to the sef-adapting filter
Filter energy is calculated for determining time delay value.
The third aspect, the embodiment of the present invention provide a kind of storage medium, are stored with one or more in the storage medium
Including the program executed instruction, it is described execute instruction can by electronic equipment (including but not limited to computer, server, or
Network equipment etc.) it reads and executes, for executing any of the above-described delay time estimation method of the present invention.
Fourth aspect provides a kind of electronic equipment comprising: at least one processor, and with described at least one
Manage the memory of device communication connection, wherein the memory is stored with the instruction that can be executed by least one described processor, institute
It states instruction to be executed by least one described processor, so that at least one described processor is able to carry out any of the above-described of the present invention
Delay time estimation method.
5th aspect, the embodiment of the present invention also provide a kind of computer program product, and the computer program product includes
The computer program of storage on a storage medium, the computer program includes program instruction, when described program instruction is calculated
When machine executes, the computer is made to execute any of the above-described delay time estimation method.
The embodiment of the present invention is determined by using sef-adapting filter is wrapped in the collected Mike's signal of microphone institute
Contain to the associated related frequency domain reference signal of reference signal, and further by way of calculating sef-adapting filter energy
Determine the time delay value of reference signal.It solves in the prior art using cross-correlation method, performance is sharply under environment interference situation
Decline, in the case that the more complicated scene of environment or it is double say the highly unstable problem of time delay estimated result.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are some embodiments of the invention, for this field
For those of ordinary skill, without creative efforts, it is also possible to obtain other drawings based on these drawings.
Fig. 1 is the flow chart of an embodiment of delay time estimation method of the invention；
Fig. 2 is the flow chart of an embodiment of step S10 of the invention；
Fig. 3 is the flow chart of an embodiment of step S30 of the invention；
Fig. 4 is the flow chart of another embodiment of step S30 of the invention；
Fig. 5 is the flow chart of the another embodiment of step S30 of the invention；
Fig. 6 is the flow chart of another embodiment of delay time estimation method of the invention；
Fig. 7 is the functional block diagram of an embodiment of time delay estimating system of the invention；
Fig. 8 is the functional block diagram of an embodiment of the signal acquisition program module in time delay estimating system of the invention；
Fig. 9 determines the functional block diagram of an embodiment of program module for the time delay in time delay estimating system of the invention；
Figure 10 be invention time delay estimating system in time delay determine program module another embodiment functional block diagram；
Figure 11 be invention time delay estimating system in time delay determine program module another embodiment functional block diagram；
Figure 12 is the structural schematic diagram of an embodiment of electronic equipment of the invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
It should be noted that in the absence of conflict, the features in the embodiments and the embodiments of the present application can phase
Mutually combination.
The present invention can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, member
Part, data structure etc..The present invention can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
In the present invention, the fingers such as " module ", " device ", " system " are applied to the related entities of computer, such as hardware, hardware
Combination, software or software in execution with software etc..In detail, for example, element can with but be not limited to run on processing
Process, processor, object, executable element, execution thread, program and/or the computer of device.In addition, running on server
Application program or shell script, server can be element.One or more elements can be in the process and/or thread of execution
In, and element can be localized and/or be distributed between two or multiple stage computers on one computer, and can be by each
Kind computer-readable medium operation.Element can also according to the signal with one or more data packets, for example, from one with
Another element interacts in local system, distributed system, and/or the network in internet passes through signal and other system interactions
The signals of data communicated by locally and/or remotely process.
Finally, it is to be noted that, herein, relational terms such as first and second and the like be used merely to by
One entity or operation are distinguished with another entity or operation, without necessarily requiring or implying these entities or operation
Between there are any actual relationship or orders.Moreover, the terms "include", "comprise", not only include those elements, and
And further include other elements that are not explicitly listed, or further include for this process, method, article or equipment institute it is intrinsic
Element.In the absence of more restrictions, the element limited by sentence " including ... ", it is not excluded that including described want
There is also other identical elements in the process, method, article or equipment of element.
As shown in Figure 1, the embodiment of the present invention provides a kind of delay time estimation method, this method for electronic equipment (for example,
Intelligent sound box, TV set-top box of voice control etc.) echo cancellor, method includes the following steps:
S10, reference signal and microphone Mike's signal collected are obtained, and carries out Fast Fourier Transform (FFT) to obtain frequency
Domain reference signal and frequency domain Mike's signal；
S20, the frequency domain reference signal is input to sef-adapting filter, to obtain being wrapped in the frequency domain Mike signal
The related frequency domain reference signal corresponding to the frequency domain reference signal contained, the frequency domain Mike signal are described adaptive for updating
Answer filter；
S30, sef-adapting filter energy is calculated according to the related frequency domain reference signal of sef-adapting filter output
Amount is for determining time delay value.
The embodiment of the present invention is determined by using sef-adapting filter is wrapped in the collected Mike's signal of microphone institute
Contain to the associated related frequency domain reference signal of reference signal, and further by way of calculating sef-adapting filter energy
Determine the time delay value of reference signal.It solves in the prior art using cross-correlation method, performance is sharply under environment interference situation
Decline, in the case that the more complicated scene of environment or it is double say the highly unstable problem of time delay estimated result.
Further, since being the block carried out based on frequency domain reference signal and frequency domain Mike signal to sef-adapting filter in frequency domain
A update is updated rather than, the complexity of sef-adapting filter update is reduced, to reduce the complexity of estimation time delay.It is logical
It crosses and adaptive-filtering is carried out to reference signal, obtain and frequency domain reference associated with reference signal included in Mike's signal
Signal most matches the block or corresponding sampled point of (most close), to be delayed, after constantly adaptively adding centainly
Processing can obtain a more stable result.
As shown in Fig. 2, above-mentioned steps S10 in some embodiments of the invention: the acquisition reference signal and microphone
Mike's signal collected, and Fast Fourier Transform (FFT) is carried out to obtain frequency domain reference signal and frequency domain Mike's signal includes:
S11, it obtains pre-stored reference signal and passes through the acquired Mike's signal of microphone；Illustratively, for intelligence
For energy speaker, during playing song, the song of broadcasting is stored as reference signal, passes through the Mike institute of intelligent sound box
Collected user instruction voice is that Mike acquires signal.
S12, the reference signal and Mike's signal are input to low-pass filter and be filtered；Low-pass filtering
Device uses the FIR of 15 ranks, it is contemplated that complexity and stability, it is down-sampled using octuple down-sampled.
S13, respectively to after being filtered reference signal and Mike's signal carry out it is down-sampled；Specifically, in order to prevent
Frequency aliasing, using low-pass filter to after being filtered Mike's signal and reference signal carry out low-pass filtering treatment, it
Carry out again afterwards it is down-sampled, obtain it is down-sampled after signal.The effect for first passing through filter is to filter high-frequency signal, to make
The phenomenon of aliasing will not be generated when frequency spectrum extends to the outside during down-sampled, aliasing can generate high-frequency signal and be converted to low frequency
Band.
S14, fast Fourier change is carried out to down-sampled obtained down-sampled reference signal and down-sampled Mike's signal respectively
It changes to obtain the frequency domain reference signal and frequency domain Mike's signal.To after down-sampled Mike's signal and reference signal by time domain
It is transformed into frequency domain, to reduce the complexity of data processing.Wherein, the frame length of Fast Fourier Transform (FFT) (FFT) is 128 or 256
Either 512 or other sizes, this is not limited by the present invention.
As shown in figure 3, above-mentioned steps S30 in some embodiments of the invention: described according to the sef-adapting filter
The related frequency domain reference signal of output calculates filter energy for determining that time delay value includes:
S31, the energy for calculating every blocking filter block according to the related frequency domain reference signal under frequency domain；For example, in frequency
The value of filter block is many complex points, [a under domain1+j*b1,a2+j*b2,...,an+j*bn], then the energy of the block can be with table
It is shown as (a1 2+b1 2+a2 2+b2 2+...+an 2+bn 2) namely absolute value of a complex number square or it is square of amplitude.
S32, the time delay value is determined according to the maximum value in the energy of every blocking filter.For example, filter block
Energy is once [1,2,4,2,1], then maximum energy was 4 (from left to right the 3rd value), then corresponding time delay is 3 (i.e.
It is 4 index value), indicate that time delay block number is 3.
The embodiment of the present invention is under frequency domain for determining time delay value, required meter by way of calculating filter block energy
The data volume for calculating processing is few, reduces the complexity of estimation time delay value.
As shown in figure 4, above-mentioned steps S30 in some embodiments of the invention: described according to the sef-adapting filter
The related frequency domain reference signal of output calculates filter energy for determining that time delay value includes:
S31 ', inverse Fourier transform is carried out to obtain corresponding relevant time domain with reference to letter to the related frequency domain reference signal
Number；
S32 ', the energy for calculating each sampled point according to the relevant time domain reference signal under time domain；For example, input
A string of sampled points [1,2,3,4], then the energy of each sampled point is 12, 22, 32, 42, that is, correspond to square of sampled point.
S33 ', the time delay value is determined according to the maximum value in the energy of each sampled point.For example, it is assumed that each adopting
Sampling point energy is [1,2,4,2,1], then maximum energy was 4 (from left to right the 3rd value), then corresponding time delay is 3 (i.e.
It is 4 index value), indicate that time delay sampled point is 3, it is with block time delay the difference is that unit is different, as fruit block time delay will turn
Sampled point time delay is changed to, then needing to be multiplied by every piece of size on the basis of block time delay.
In the embodiment of the present invention under time domain, estimation time delay value, Ke Yitong are used for by way of calculating sampled point energy
The estimation of numerous sampled point energy under time domain is crossed to promote the precision estimated time delay value.
As shown in figure 5, above-mentioned steps S30 in some embodiments of the invention: described according to the sef-adapting filter
The related frequency domain reference signal of output calculates filter energy for determining that time delay value includes:
S31 ', the energy for calculating every blocking filter block according to the related frequency domain reference signal under frequency domain；
S32 ', the first time delay value is determined according to the maximum value in the energy of every blocking filter；
S33 ', inverse Fourier transform is carried out to obtain corresponding relevant time domain with reference to letter to the related frequency domain reference signal
Number；
S34 ', the energy for calculating each sampled point according to the relevant time domain reference signal under time domain；
S35 ', the second time delay value is determined according to the maximum value in the energy of each sampled point；
S36 ', the time delay value is determined according to first time delay value and second time delay value.For example, every blocking filter
Block has 512 points, if sampled point time delay estimated result is 1024, block time delay result is 1, and sampled point time delay is transformed into block time delay
1024/512=2, not equal to the estimated value of block time delay, then current estimated result cancels, this moment still exports last moment
Time delay value.If the result of this when of block time delay estimation is also 2, the result pair of block time delay and the estimation of sampled point time delay
It answers, exportable current result.
In the embodiment of the present invention, sampled point delay and block delay are to need to carry out corresponding to peak value searching filter W value
Which point be exactly to be delayed, it should be the result of peak value searching.The embodiment of the present invention comprehensively considers sampled point delay and block prolongs
When, the two mutually as reference, further improves the precision of time delay estimation.
As shown in fig. 6, the flow chart of the FIELD Embodiments for delay time estimation method of the invention, specifically includes following step
It is rapid: down-sampled, Fourier transformation, adaptive-filtering, peak value searching and post-processing.Input two paths of signals (Mike's signal and reference
Signal), export every frame time delay.Each step is introduced respectively below.
1) down-sampled
Low-pass filtering first is carried out to reference signal and Mike's signal, then carries out down-sampled (algorithm complexity can be reduced).
Wherein, to prevent frequency aliasing, a low-pass filter is first passed through, here using the FIR of 15 ranks or 7 ranks, this hair
It is bright that this is not construed as limiting.And sampling consideration complexity and stability, down-sampled or 4 times are down-sampled using 8 times, the present invention is to this
It is not construed as limiting.Mike's signal and reference signal carry out down-sampled together, guarantee that data length per treatment is consistent, including low pass
Filtering and sampling.
The effect for first passing through filter is to filter high-frequency signal, to make when frequency spectrum extends to the outside during down-sampled
The phenomenon of aliasing will not be generated, aliasing can generate high-frequency signal and be converted to low-frequency band.
2) FFT (Fourier transformation)
Down-sampled later Mike's signal and reference signal carry out FFT transform (Fourier transformation) respectively.It is multiple in order to reduce
Miscellaneous degree, in frequency domain processing, it is therefore desirable to which the frame length of FFT, FFT are 128 or 256 or 512 or other sizes, the present invention couple
This is not construed as limiting.
3) adaptive-filtering
Linear convolution mode is replaced used here as cyclic convolution, is realized using overlap-save method, 50% overlapping is used.It is adaptive
The purpose that should be filtered is to estimate part relevant to reference signal in Mike's signal.The input of adaptive-filtering
It is reference signal, output is relevant part in Mike's signal of estimation.
For kth blocking filter and reference signal, reference signal filtering output result are as follows:
The later half element of y (k)=IFFT [X (k) W (k)],
Wherein, X (k) is distal end block signal, and W (k) is filter block coefficient, and why only retains later half element, is
Because later half element is the result of cyclic convolution.Wherein, element is exactly sampled point, and distal end block signal is exactly the reference before referring to
One piece one piece of the form that signal is divided into.
Time domain block error signal are as follows:
E (k)=d (k)-y (k),
Wherein, d (k) indicates Mike's signal.
Area block error signal are as follows:
E (k)=FFT [0e (k)],
Wherein, 00 of the addition half before e (k) is indicated.
E (k) is normalized, is obtained
Wherein, | X (k) | indicate the smoothed energy of reference signal,δIt is the fixed value for preventing filter divergence.
The renewal amount of filter are as follows:
The first half element,
Wherein, μ is the step-size factor factor, because only that the result of the first half is only correct result, later half needs to abandon.
This method means overlap-save method and only retains a part.
Filter update formula are as follows:
W (k+1)=W (k)+FFT [Φ (k) 0],
Wherein 0 indicate to add the 0 of half afterwards in Φ (k).
The update of filter is the crucial step of comparison, but Φ (k) therein and errorIt is related, andWith it is close
The estimated value of end signal, reference signal in proximal end is related, and estimated value of the reference signal in proximal end needs to carry out reference signal filter
Wave.So this is the process of a continuous iteration.
4) peak value searching
Here consider the delay of output sampled point and block delay, need two different branches.Wherein,
The requirements of process of sampled point time delay first carries out IFFT to entire filter and transforms to time domain, and it is every then to calculate filter
The energy of a coefficient chooses point time delay of that the maximum sampled point of energy as estimation.
And the process of block delay is the energy in the every blocking filter block of frequency-domain calculations, then takes maximum piece of energy to be used as block
Time delay.
Here sampled point delay and block delay is to need which point corresponding to peak value searching carried out to filter W value
It is exactly delayed, it should be the result of peak value searching.Point delay and block delay can comprehensively consider or individually consider for estimating
Timing is prolonged.The required precision of block delay is not high, and complexity is low, but error is big；Sampled point delay, error is small, but more complicated.The two
It can also estimate together, can further promote precision mutually as reference.
5) it post-processes
Post-processing mainly accounts in terms of two: (1) the case where time delay is except filter length；(2) in short-term when
The case where prolonging abnormal shake.
Alleviating the current means of first problem is to carry out short-time analysis with continuous 20 frame, obtains filter to every frame
Average energy and energy peak, be then averaged this 20 frame, compares its average energy and peak energy, if peak energy < average
Some multiple proportion of energy, the estimation for being considered as 20 frame are all insincere.20 frames why are carried out once to count, be prevent it is short
When it is interior there is unexpected fragmentary time delay value, this adjusts subsequent AEC filter nonsensical, can also reduce certain complexity.
Alleviating the method that Second Problem uses is more currently to estimate delay and the before delay of the averaged power spectrum of 20 frames,
If it exceeds some threshold value, that is, think that present frame estimated result is insincere, that just exports former frame result.If it is considered to present frame
It is credible, and continuous multiple frames are credible, that is, export the time delay result of present frame and be multiplied by down-sampled multiple adopting as final output
Sampling point time delay result.
And the post-processing about block delay is more relatively easy, because the estimated result floating of block time delay is smaller, excessively
Post-processing will affect accuracy instead, therefore only consider a problem (1), the average energy of energy and entire filter to current block
Amount is compared, if it is less than some threshold value, then it is assumed that current time delay result is invalid.
Frame concept mentioned here is exactly the sampling number of the Mike's signal disposably handled and reference signal or is length
Degree.Since it is desired that output in real time, so a segment of audio cannot be all handled, so needing to export in processing.Here
As the concept of frame with the length of block is before, but filter block as several forms.And the input post-processed is
The time delay result for each frame estimated before.Here be that smooth equal operation is carried out to multiframe result, post-processing be substantially in order to
Make result more stable.Output is also time delay as a result, being only a comprehensive output in conjunction with multiframe.
It should be noted that for the various method embodiments described above, for simple description, therefore, it is stated as a series of
Movement merge, but those skilled in the art should understand that, the present invention is not limited by the sequence of acts described because
According to the present invention, some steps may be performed in other sequences or simultaneously.Secondly, those skilled in the art should also know
It knows, the embodiments described in the specification are all preferred embodiments, and related actions and modules is not necessarily of the invention
It is necessary.In the above-described embodiments, it all emphasizes particularly on different fields to the description of each embodiment, there is no the portion being described in detail in some embodiment
Point, reference can be made to the related descriptions of other embodiments.
As shown in fig. 7, the embodiment of the present invention also provides a kind of time delay estimating system 700, comprising:
Signal acquisition program module 710 for obtaining reference signal and microphone Mike's signal collected, and carries out fast
Fast Fourier transformation is to obtain frequency domain reference signal and frequency domain Mike's signal；
Adaptive-filtering program module 720, for the frequency domain reference signal to be input to sef-adapting filter, to obtain
Correspond to the related frequency domain reference signal of the frequency domain reference signal, the frequency domain wheat included in the frequency domain Mike signal
Gram signal is for updating the sef-adapting filter；
Time delay determines program module 730, the related frequency domain reference letter for being exported according to the sef-adapting filter
Number filter energy is calculated with for determining time delay value.
The embodiment of the present invention is determined by using sef-adapting filter is wrapped in the collected Mike's signal of microphone institute
Contain to the associated related frequency domain reference signal of reference signal, and further by way of calculating sef-adapting filter energy
Determine the time delay value of reference signal.It solves in the prior art using cross-correlation method, performance is sharply under environment interference situation
Decline, in the case that the more complicated scene of environment or it is double say the highly unstable problem of time delay estimated result.
Further, since being the block carried out based on frequency domain reference signal and frequency domain Mike signal to sef-adapting filter in frequency domain
A update is updated rather than, the complexity of sef-adapting filter update is reduced, to reduce the complexity of estimation time delay.It is logical
It crosses and adaptive-filtering is carried out to reference signal, obtain and frequency domain reference associated with reference signal included in Mike's signal
Signal most matches the block or corresponding sampled point of (most close), to be delayed, after constantly adaptively adding centainly
Processing can obtain a more stable result.
As shown in figure 8, in some embodiments of the present invention, signal acquisition program module 710 includes:
Signal acquisition program unit 711, for obtaining pre-stored reference signal and by the acquired Mike of microphone
Signal；
Program unit 712 is filtered, for the reference signal and Mike's signal to be input to low-pass filter
It is filtered；
Down-sampled program unit 713, for respectively to after being filtered reference signal and Mike's signal carry out drop adopt
Sample；
Fourier transformation program unit 714, for respectively to down-sampled obtained down-sampled reference signal and down-sampled wheat
Gram signal carries out Fast Fourier Transform (FFT) to obtain the frequency domain reference signal and frequency domain Mike's signal.
As shown in figure 9, time delay determines that program module 730 includes: in some embodiments of the present invention
Energy balane program unit 731, for calculating every piece of filtering according to the related frequency domain reference signal under frequency domain
The energy of device block；
Time delay determines program unit 732, when determining described for the maximum value in the energy according to every blocking filter
Prolong value.
As shown in Figure 10, in some embodiments of the present invention, time delay determines that program module 730 includes:
Signal conversion procedure unit 731 ', for carrying out inverse Fourier transform to the related frequency domain reference signal to obtain
Corresponding relevant time domain reference signal；
Energy balane program unit 732 ', for calculating each sampling according to the relevant time domain reference signal under time domain
The energy of point；
Time delay determines program unit 733 ', when determining described for the maximum value in the energy according to each sampled point
Prolong value.
As shown in figure 11, in some embodiments of the present invention, time delay determines that program module 730 includes:
First energy balane program unit 731 ', for calculating every piece according to the related frequency domain reference signal under frequency domain
The energy of filter block；
First time delay determines program unit 732 ', determines for the maximum value in the energy according to every blocking filter
One time delay value；
Inverse Fourier transform program unit 733 ', for the related frequency domain reference signal carry out inverse Fourier transform with
Obtain corresponding relevant time domain reference signal；
Second energy balane program unit 734 ', it is each for being calculated under time domain according to the relevant time domain reference signal
The energy of sampled point；
Second time delay determines program unit 735 ', determines for the maximum value in the energy according to each sampled point
Two time delay values；
Time delay determines program unit 736 ', when for determining described according to first time delay value and second time delay value
Prolong value.
In some embodiments, the embodiment of the present invention provides a kind of non-volatile computer readable storage medium storing program for executing, described to deposit
Being stored in storage media one or more includes the programs executed instruction, it is described execute instruction can by electronic equipment (including but
It is not limited to computer, server or the network equipment etc.) it reads and executes, for executing any of the above-described time delay of the present invention
Estimation method.
In some embodiments, the embodiment of the present invention also provides a kind of computer program product, and the computer program produces
Product include the computer program being stored on non-volatile computer readable storage medium storing program for executing, and the computer program includes that program refers to
It enables, when described program instruction is computer-executed, the computer is made to execute any of the above-described delay time estimation method.
In some embodiments, the embodiment of the present invention also provides a kind of electronic equipment comprising: at least one processor,
And the memory being connect at least one described processor communication, wherein the memory is stored with can be by described at least one
The instruction that a processor executes, described instruction is executed by least one described processor, so that at least one described processor energy
Enough execution time delay estimation methods.
In some embodiments, the embodiment of the present invention also provides a kind of storage medium, is stored thereon with computer program,
It is characterized in that, delay time estimation method when which is executed by processor.
The time delay estimating system of the embodiments of the present invention can be used for executing the delay time estimation method of the embodiment of the present invention, and
Reach the realization delay time estimation method technical effect achieved of the embodiments of the present invention accordingly, which is not described herein again.This
Related function module can be realized by hardware processor (hardware processor) in inventive embodiments.
Figure 12 is the hardware configuration signal of the electronic equipment for the execution time delay estimation method that another embodiment of the application provides
Figure, as shown in figure 12, which includes:
One or more processors 1210 and memory 1220, in Figure 12 by taking a processor 1210 as an example.
The equipment of execution time delay estimation method can also include: input unit 1230 and output device 1240.
Processor 1210, memory 1220, input unit 1230 and output device 1240 can by bus or other
Mode connects, in Figure 12 for being connected by bus.
Memory 1220 is used as a kind of non-volatile computer readable storage medium storing program for executing, can be used for storing non-volatile software journey
Sequence, non-volatile computer executable program and module, such as the corresponding program of delay time estimation method in the embodiment of the present application
Instruction/module.Non-volatile software program, instruction and the mould that processor 1210 is stored in memory 1220 by operation
Block, thereby executing the various function application and data processing of server, i.e. realization above method embodiment delay time estimation method.
Memory 1220 may include storing program area and storage data area, wherein storing program area can store operation system
Application program required for system, at least one function；Storage data area can be stored to be created according to using for time delay estimation device
Data etc..In addition, memory 1220 may include high-speed random access memory, it can also include nonvolatile memory,
A for example, at least disk memory, flush memory device or other non-volatile solid state memory parts.In some embodiments,
Optional memory 1220 includes the memory remotely located relative to processor 1210, these remote memories can pass through network
It is connected to time delay estimation device.The example of above-mentioned network includes but is not limited to internet, intranet, local area network, moves and lead to
Letter net and combinations thereof.
Input unit 1230 can receive the number or character information of input, and generates and set with the user of time delay estimation device
It sets and the related signal of function control.Output device 1240 may include that display screen etc. shows equipment.
One or more of modules are stored in the memory 1220, when by one or more of processors
When 1210 execution, the delay time estimation method in above-mentioned any means embodiment is executed.
Method provided by the embodiment of the present application can be performed in the said goods, has the corresponding functional module of execution method and has
Beneficial effect.The not technical detail of detailed description in the present embodiment, reference can be made to method provided by the embodiment of the present application.
The electronic equipment of the embodiment of the present application exists in a variety of forms, including but not limited to:
(1) mobile communication equipment: the characteristics of this kind of equipment is that have mobile communication function, and to provide speech, data
Communication is main target.This Terminal Type includes: smart phone (such as iPhone), multimedia handset, functional mobile phone and low
Hold mobile phone etc..
(2) super mobile personal computer equipment: this kind of equipment belongs to the scope of personal computer, there is calculating and processing function
Can, generally also have mobile Internet access characteristic.This Terminal Type includes: PDA, MID and UMPC equipment etc., such as iPad.
(3) portable entertainment device: this kind of equipment can show and play multimedia content.Such equipment include: audio,
Video player (such as iPod), handheld device, e-book and intelligent toy and portable car-mounted navigation equipment.
(4) server: providing the equipment of the service of calculating, and the composition of server includes that processor, hard disk, memory, system are total
Line etc., server is similar with general computer architecture, but due to needing to provide highly reliable service, in processing energy
Power, stability, reliability, safety, scalability, manageability etc. are more demanding.
(5) other electronic devices with data interaction function.
The apparatus embodiments described above are merely exemplary, wherein described, unit can as illustrated by the separation member
It is physically separated with being or may not be, component shown as a unit may or may not be physics list
Member, it can it is in one place, or may be distributed over multiple network units.It can be selected according to the actual needs
In some or all of the modules achieve the purpose of the solution of this embodiment.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It is realized by the mode of software plus general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, above-mentioned technology
Scheme substantially in other words can be embodied in the form of software products the part that the relevant technologies contribute, the computer
Software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including some instructions to
So that computer equipment (can be personal computer, server or the network equipment etc.) execute each embodiment or
Method described in certain parts of embodiment.
Finally, it should be noted that above embodiments are only to illustrate the technical solution of the application, rather than its limitations；Although
The application is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, each embodiment technical solution of the application that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (12)
1. a kind of delay time estimation method, comprising:
Reference signal and microphone Mike's signal collected are obtained, and carries out Fast Fourier Transform (FFT) to obtain frequency domain reference letter
Number and frequency domain Mike signal；
The frequency domain reference signal is input to sef-adapting filter, it is corresponding included in the frequency domain Mike signal to obtain
In the related frequency domain reference signal of the frequency domain reference signal, the frequency domain Mike signal is for updating the adaptive-filtering
Device；
Filter energy is calculated for when determining according to the related frequency domain reference signal of sef-adapting filter output
Prolong value.
2. according to the method described in claim 1, wherein, the acquisition reference signal and microphone Mike's signal collected,
And Fast Fourier Transform (FFT) is carried out to obtain frequency domain reference signal and frequency domain Mike's signal includes:
Obtain pre-stored reference signal and by the acquired Mike's signal of microphone；
The reference signal and Mike's signal are input to low-pass filter and be filtered；
Respectively to after being filtered reference signal and Mike's signal carry out it is down-sampled；
Fast Fourier Transform (FFT) is carried out to obtain to down-sampled obtained down-sampled reference signal and down-sampled Mike's signal respectively
The frequency domain reference signal and frequency domain Mike's signal.
3. according to the method described in claim 1, wherein, the related frequency domain according to sef-adapting filter output
Reference signal calculates filter energy for determining that time delay value includes:
The energy of every blocking filter block is calculated according to the related frequency domain reference signal under frequency domain；
The time delay value is determined according to the maximum value in the energy of every blocking filter.
4. according to the method described in claim 1, wherein, the related frequency domain according to sef-adapting filter output
Reference signal calculates filter energy for determining that time delay value includes:
Inverse Fourier transform is carried out to obtain corresponding relevant time domain reference signal to the related frequency domain reference signal；
The energy of each sampled point is calculated according to the relevant time domain reference signal under time domain；
The time delay value is determined according to the maximum value in the energy of each sampled point.
5. according to the method described in claim 1, wherein, the related frequency domain according to sef-adapting filter output
Reference signal calculates filter energy for determining that time delay value includes:
The energy of every blocking filter block is calculated according to the related frequency domain reference signal under frequency domain；
The first time delay value is determined according to the maximum value in the energy of every blocking filter；
Inverse Fourier transform is carried out to obtain corresponding relevant time domain reference signal to the related frequency domain reference signal；
The energy of each sampled point is calculated according to the relevant time domain reference signal under time domain；
The second time delay value is determined according to the maximum value in the energy of each sampled point；
The time delay value is determined according to first time delay value and second time delay value.
6. a kind of time delay estimating system, comprising:
Signal acquisition program module for obtaining reference signal and microphone Mike's signal collected, and carries out in quick Fu
Leaf transformation is to obtain frequency domain reference signal and frequency domain Mike's signal；
Adaptive-filtering program module, for the frequency domain reference signal to be input to sef-adapting filter, to obtain the frequency
Correspond to the related frequency domain reference signal of the frequency domain reference signal, the frequency domain Mike signal included in the Mike's signal of domain
For updating the sef-adapting filter；
Time delay determines program module, and the related frequency domain reference signal for being exported according to the sef-adapting filter calculates filter
Wave device energy is for determining time delay value.
7. system according to claim 1, wherein the signal acquisition program module includes:
Signal acquisition program unit, for obtaining pre-stored reference signal and by the acquired Mike's signal of microphone；
Program unit is filtered, for the reference signal and Mike's signal to be input to low-pass filter and be filtered
Processing；
Down-sampled program unit, for respectively to after being filtered reference signal and Mike's signal carry out it is down-sampled；
Fourier transformation program unit, for respectively to down-sampled obtained down-sampled reference signal and down-sampled Mike's signal into
Row Fast Fourier Transform (FFT) is to obtain the frequency domain reference signal and frequency domain Mike's signal.
8. system according to claim 1, wherein the time delay determines that program module includes:
Energy balane program unit, for calculating the energy of every blocking filter block according to the related frequency domain reference signal under frequency domain
Amount；
Time delay determines program unit, determines the time delay value for the maximum value in the energy according to every blocking filter.
9. system according to claim 1, wherein the time delay determines that program module includes:
Signal conversion procedure unit, for carrying out inverse Fourier transform to the related frequency domain reference signal to obtain corresponding phase
Close time domain reference signal；
Energy balane program unit, for calculating the energy of each sampled point according to the relevant time domain reference signal under time domain
Amount；
Time delay determines program unit, determines the time delay value for the maximum value in the energy according to each sampled point.
10. system according to claim 1, wherein the time delay determines that program module includes:
First energy balane program unit, for calculating every blocking filter block according to the related frequency domain reference signal under frequency domain
Energy；
First time delay determines program unit, determines the first time delay for the maximum value in the energy according to every blocking filter
Value；
Inverse Fourier transform program unit, for carrying out inverse Fourier transform to the related frequency domain reference signal to be corresponded to
Relevant time domain reference signal；
Second energy balane program unit, for calculating each sampled point according to the relevant time domain reference signal under time domain
Energy；
Second time delay determines program unit, determines the second time delay for the maximum value in the energy according to each sampled point
Value；
Time delay determines program unit, for determining the time delay value according to first time delay value and second time delay value.
11. a kind of electronic equipment comprising: at least one processor, and connect at least one described processor communication
Memory, wherein the memory be stored with can by least one described processor execute instruction, described instruction by it is described extremely
A few processor executes, so that at least one described processor is able to carry out side described in any one of claim 1-5
The step of method.
12. a kind of storage medium, is stored thereon with computer program, which is characterized in that the realization when program is executed by processor
The step of method described in any one of claim 1-5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811049712.7A CN109102821B (en) | 2018-09-10 | 2018-09-10 | Time delay estimation method, time delay estimation system, storage medium and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811049712.7A CN109102821B (en) | 2018-09-10 | 2018-09-10 | Time delay estimation method, time delay estimation system, storage medium and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109102821A true CN109102821A (en) | 2018-12-28 |
| CN109102821B CN109102821B (en) | 2021-05-25 |
Family
ID=64865653
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811049712.7A Active CN109102821B (en) | 2018-09-10 | 2018-09-10 | Time delay estimation method, time delay estimation system, storage medium and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109102821B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109584898A (en) * | 2018-12-29 | 2019-04-05 | 上海瑾盛通信科技有限公司 | A kind of processing method of voice signal, device, storage medium and electronic equipment |
| CN109862503A (en) * | 2019-01-30 | 2019-06-07 | 北京雷石天地电子技术有限公司 | A kind of method and apparatus of loudspeaker delay adjust automatically |
| CN110349592A (en) * | 2019-07-17 | 2019-10-18 | 百度在线网络技术（北京）有限公司 | Method and apparatus for output information |
| CN110992973A (en) * | 2019-11-29 | 2020-04-10 | 维沃移动通信有限公司 | Signal time delay determination method and electronic equipment |
| CN111294473A (en) * | 2019-01-28 | 2020-06-16 | 展讯通信（上海）有限公司 | Signal processing method and device |
| CN111613238A (en) * | 2020-05-21 | 2020-09-01 | 北京百度网讯科技有限公司 | Method, device and equipment for determining time delay between signals and storage medium |
Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3369077A (en) * | 1964-06-09 | 1968-02-13 | Ibm | Pitch modification of audio waveforms |
| CN1691716A (en) * | 2004-04-23 | 2005-11-02 | 北京三星通信技术研究有限公司 | Echo eliminator |
| CN1953060A (en) * | 2006-11-24 | 2007-04-25 | 北京中星微电子有限公司 | Echo elimination device for microphone and method thereof |
| CN101026659A (en) * | 2006-02-23 | 2007-08-29 | 中兴通讯股份有限公司 | Method for realizing echo time delay positioning |
| CN101119135A (en) * | 2007-07-04 | 2008-02-06 | 深圳市融合视讯科技有限公司 | Step parameter regulation means and equipment for eliminating echo |
| US20100004937A1 (en) * | 2008-07-03 | 2010-01-07 | Thomson Licensing | Method for time scaling of a sequence of input signal values |
| CN103700374A (en) * | 2013-12-25 | 2014-04-02 | 宁波菊风系统软件有限公司 | Method for determining system time delay in acoustic echo cancellation and acoustic echo cancellation method |
| CN104038181A (en) * | 2014-06-05 | 2014-09-10 | 北京航空航天大学 | Self-adapting filter construction method based on NLMS algorithm |
| CN104751853A (en) * | 2013-12-31 | 2015-07-01 | 联芯科技有限公司 | Double-microphone noise inhibiting method and system |
| CN105472191A (en) * | 2015-11-18 | 2016-04-06 | 百度在线网络技术（北京）有限公司 | Method and device for tracking echo time delay |
| CN105810202A (en) * | 2014-12-31 | 2016-07-27 | 展讯通信（上海）有限公司 | Method and device for reducing echoes, and communication equipment |
| CN106936407A (en) * | 2017-01-12 | 2017-07-07 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Area block minimum mean square self-adaption filtering method |
| CN107123430A (en) * | 2017-04-12 | 2017-09-01 | 广州视源电子科技股份有限公司 | Echo cancel method, device, meeting flat board and computer-readable storage medium |
| CN107610713A (en) * | 2017-10-23 | 2018-01-19 | 科大讯飞股份有限公司 | Echo cancel method and device based on time delay estimation |
| CN107645697A (en) * | 2016-07-20 | 2018-01-30 | 森海塞尔通信公司 | Sef-adapting filter unit as Echo Canceller |
| CN107785026A (en) * | 2017-10-18 | 2018-03-09 | 会听声学科技（北京）有限公司 | A kind of delay time estimation method eliminated for set top box indoor echo |
-
2018
- 2018-09-10 CN CN201811049712.7A patent/CN109102821B/en active Active
Patent Citations (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3369077A (en) * | 1964-06-09 | 1968-02-13 | Ibm | Pitch modification of audio waveforms |
| CN1691716A (en) * | 2004-04-23 | 2005-11-02 | 北京三星通信技术研究有限公司 | Echo eliminator |
| CN101026659A (en) * | 2006-02-23 | 2007-08-29 | 中兴通讯股份有限公司 | Method for realizing echo time delay positioning |
| CN1953060A (en) * | 2006-11-24 | 2007-04-25 | 北京中星微电子有限公司 | Echo elimination device for microphone and method thereof |
| CN101119135A (en) * | 2007-07-04 | 2008-02-06 | 深圳市融合视讯科技有限公司 | Step parameter regulation means and equipment for eliminating echo |
| US20100004937A1 (en) * | 2008-07-03 | 2010-01-07 | Thomson Licensing | Method for time scaling of a sequence of input signal values |
| CN103700374A (en) * | 2013-12-25 | 2014-04-02 | 宁波菊风系统软件有限公司 | Method for determining system time delay in acoustic echo cancellation and acoustic echo cancellation method |
| CN104751853A (en) * | 2013-12-31 | 2015-07-01 | 联芯科技有限公司 | Double-microphone noise inhibiting method and system |
| CN104038181A (en) * | 2014-06-05 | 2014-09-10 | 北京航空航天大学 | Self-adapting filter construction method based on NLMS algorithm |
| CN105810202A (en) * | 2014-12-31 | 2016-07-27 | 展讯通信（上海）有限公司 | Method and device for reducing echoes, and communication equipment |
| CN105472191A (en) * | 2015-11-18 | 2016-04-06 | 百度在线网络技术（北京）有限公司 | Method and device for tracking echo time delay |
| CN107645697A (en) * | 2016-07-20 | 2018-01-30 | 森海塞尔通信公司 | Sef-adapting filter unit as Echo Canceller |
| CN106936407A (en) * | 2017-01-12 | 2017-07-07 | 西南电子技术研究所（中国电子科技集团公司第十研究所） | Area block minimum mean square self-adaption filtering method |
| CN107123430A (en) * | 2017-04-12 | 2017-09-01 | 广州视源电子科技股份有限公司 | Echo cancel method, device, meeting flat board and computer-readable storage medium |
| CN107785026A (en) * | 2017-10-18 | 2018-03-09 | 会听声学科技（北京）有限公司 | A kind of delay time estimation method eliminated for set top box indoor echo |
| CN107610713A (en) * | 2017-10-23 | 2018-01-19 | 科大讯飞股份有限公司 | Echo cancel method and device based on time delay estimation |
Non-Patent Citations (3)
| Title |
|---|
| T. KOSAKA ET AL.: "A novel frequency domain filtered-X LMS algorithm for active noise reduction", 《1997 IEEE INTERNATIONAL CONFERENCE ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING》 * |
| 刘庆华: "室内噪声环境下语音信号时延估计", 《声学与电子工程》 * |
| 董蔚: "数字助听器中时延估计与语音增强技术的研究", 《中国优秀硕士学位论文全文库数据库（电子期刊）》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109584898A (en) * | 2018-12-29 | 2019-04-05 | 上海瑾盛通信科技有限公司 | A kind of processing method of voice signal, device, storage medium and electronic equipment |
| CN111294473A (en) * | 2019-01-28 | 2020-06-16 | 展讯通信（上海）有限公司 | Signal processing method and device |
| CN111294473B (en) * | 2019-01-28 | 2022-01-04 | 展讯通信（上海）有限公司 | Signal processing method and device |
| CN109862503A (en) * | 2019-01-30 | 2019-06-07 | 北京雷石天地电子技术有限公司 | A kind of method and apparatus of loudspeaker delay adjust automatically |
| CN110349592A (en) * | 2019-07-17 | 2019-10-18 | 百度在线网络技术（北京）有限公司 | Method and apparatus for output information |
| CN110992973A (en) * | 2019-11-29 | 2020-04-10 | 维沃移动通信有限公司 | Signal time delay determination method and electronic equipment |
| CN111613238A (en) * | 2020-05-21 | 2020-09-01 | 北京百度网讯科技有限公司 | Method, device and equipment for determining time delay between signals and storage medium |
| CN111613238B (en) * | 2020-05-21 | 2023-09-19 | 阿波罗智联(北京)科技有限公司 | Method, device, equipment and storage medium for determining delay between signals |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109102821B (en) | 2021-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109102821A (en) | Delay time estimation method, system, storage medium and electronic equipment | |
| CN107123430B (en) | Echo cancel method, device, meeting plate and computer storage medium | |
| US11349525B2 (en) | Double talk detection method, double talk detection apparatus and echo cancellation system | |
| WO2017099728A1 (en) | System and method for suppression of non-linear acoustic echoes | |
| CN109473118A (en) | Double-channel pronunciation Enhancement Method and device | |
| CN104883462B (en) | A kind of sef-adapting filter and filtering method for eliminating acoustic echo | |
| CN105989850B (en) | A kind of echo cancellation method and device | |
| CN104994249B (en) | Sound method for echo cancellation and device | |
| US20160182770A1 (en) | Display apparatus and echo cancellation method thereof | |
| CN110402540A (en) | Active denoising method, device, chip, active control system and storage medium | |
| CN110211602B (en) | Intelligent voice enhanced communication method and device | |
| Gil-Cacho et al. | Wiener variable step size and gradient spectral variance smoothing for double-talk-robust acoustic echo cancellation and acoustic feedback cancellation | |
| CN109727605B (en) | Method and system for processing sound signal | |
| US20180090153A1 (en) | Signal processing apparatus, signal processing method, and signal processing program | |
| CN110431624A (en) | Residual echo detection method, residual echo detection device, speech processing chip and electronic equipment | |
| CN109961798B (en) | Echo cancellation system, echo cancellation method, readable computer storage medium, and terminal | |
| Chergui et al. | A new pre-whitening transform domain LMS algorithm and its application to speech denoising | |
| CN108597527A (en) | Multichannel audio processing method, device, computer readable storage medium and terminal | |
| BR112014009647B1 (en) | NOISE Attenuation APPLIANCE AND NOISE Attenuation METHOD | |
| CN109584898B (en) | Voice signal processing method and device, storage medium and electronic equipment | |
| Pu et al. | A signal decorrelation PNLMS algorithm for double-talk acoustic echo cancellation | |
| CN108986837A (en) | A kind of filter update method and device | |
| CN114220451A (en) | Audio denoising method, electronic device, and storage medium | |
| EP2712208A1 (en) | Audio processing device, audio processing method, and recording medium on which audio processing program is recorded | |
| CN112687285B (en) | Echo cancellation 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 | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 215123 14 Tengfei Innovation Park, 388 Xinping street, Suzhou Industrial Park, Suzhou, Jiangsu. Applicant after: Sipic Technology Co.,Ltd. Address before: 215123 14 Tengfei Innovation Park, 388 Xinping street, Suzhou Industrial Park, Suzhou, Jiangsu. Applicant before: AI SPEECH Co.,Ltd. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Time delay estimation methods, systems, storage media, and electronic devices Effective date of registration: 20230726 Granted publication date: 20210525 Pledgee: CITIC Bank Limited by Share Ltd. Suzhou branch Pledgor: Sipic Technology Co.,Ltd. Registration number: Y2023980049433 |