CN101858939A - Method and device for detecting harmonic signal - Google Patents
Method and device for detecting harmonic signal Download PDFInfo
- Publication number
- CN101858939A CN101858939A CN200910134924A CN200910134924A CN101858939A CN 101858939 A CN101858939 A CN 101858939A CN 200910134924 A CN200910134924 A CN 200910134924A CN 200910134924 A CN200910134924 A CN 200910134924A CN 101858939 A CN101858939 A CN 101858939A
- Authority
- CN
- China
- Prior art keywords
- peak
- signal
- pseudo
- signal frame
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
Abstract
The embodiment of the invention discloses a method and a device for detecting a harmonic signal, which are used for accurately detecting the harmonic signal. The method comprises the following steps of: acquiring a frequency domain signal in a current signal frame; searching for the frequency domain signal to acquire first peak information; matching the first peak information and second peak information to acquire a matched numerical value, wherein the second peak information is the peak information of the frequency domain signal in the signal frame acquired before the current signal frame; and judging whether the matched numerical value exceeds a preset matched threshold, and determining the current signal frame is a harmonic signal frame if does. The embodiment of the invention also provides the device for detecting the harmonic signal. The method and the device can accurately detect the harmonic signal.
Description
Technical field
The present invention relates to the signal Processing field, relate in particular to a kind of harmonic signal detection method and device.
Background technology
In the parametric audio coding technology, usually at first can carry out Characteristics Detection to sound signal, adopt different parameter extraction technology at different characteristics then, the parameter that extracts will be encoded and pass to demoder, and demoder recovers sound signal according to the parameter that receives.The most common signal characteristic has two kinds of harmonic wave (harmonic) and anharmonic waves.
The essence of sound signal is sound wave, vibration by object produces, wherein, the minimum frequency of object vibration is known as " fundamental frequency " (fundamental frequency), and the ripple that is produced by the vibration of each multiple frequency of fundamental frequency is called as " harmonic wave " (harmonic).In the sound signal, if having abundant harmonic components, this signal is called as harmonic signal so, otherwise then is called as the anharmonic wave signal.
In the prior art, a kind of harmonic signal detection method is:
Convert time-domain signal to frequency-region signal, at first calculate the fundamental frequency of frequency-region signal, detect each peak of frequency-region signal then, calculate the distance between adjacent peak value, calculate the ratio of this distance and fundamental frequency then, when ratio during, judge that then signal is a harmonic signal, otherwise be the anharmonic wave signal greater than certain threshold value.
But, in the above-mentioned prior art, owing to need to determine fundamental frequency, so need a large amount of complicated calculations, simultaneously, in the computation process of fundamental frequency, be subjected to influence of various factors easily, for example frequency multiplication or half erroneous judgement frequently be fundamental frequency, thus the accuracy that causes harmonic wave to detect.
A kind of harmonic signal detection method also is provided in the prior art, has been specially:
Time-domain signal is converted to the correction type discrete cosine transform (MDCT of frequency domain, Modified DiscreteCosine Transform) coefficient, calculate the MDCT coefficient of present frame and the cross-correlation coefficient of former frame MDCT coefficient then, the relation that compares cross-correlation coefficient and threshold value, when greater than certain threshold value, judge that then signal is a harmonic signal, otherwise be the anharmonic wave signal.
But, in the above-mentioned prior art, comprise amplitude and phase information in the MDCT coefficient, and phase information can influence the calculating of cross-correlation coefficient, thereby make that the cross-correlation coefficient that finally obtains is not accurate enough, and cause signal to be judged by accident easily.
Summary of the invention
The embodiment of the invention provides a kind of harmonic signal detection method and device, can improve the accuracy that harmonic signal detects.
The harmonic signal detection method that the embodiment of the invention provides comprises: obtain the frequency-region signal in the current demand signal frame; Search for described frequency-region signal and obtain the first peak value information; Described first peak value information and second peak information mated obtain mating numerical value, described second peak information is the peak information of the frequency-region signal in the signal frame that got access to before described current demand signal frame; Judge whether described coupling numerical value surpasses the coupling threshold value that presets,, determine that then described current demand signal frame is the harmonic signal frame if surpass.
The harmonic signal detection method that the embodiment of the invention provides comprises: obtain the frequency-region signal in the current demand signal frame; Described frequency-region signal is carried out phase place elimination processing obtain the first pseudo-power spectrum signal; Calculate the cross-correlation coefficient between the described first pseudo-power spectrum signal and the second pseudo-power spectrum signal, the described second pseudo-power spectrum signal is eliminated the pseudo-power spectrum signal that obtains after the processing for the frequency-region signal in the signal frame that got access to carries out phase place before described current demand signal frame; Judge whether described cross-correlation coefficient surpasses the coefficient threshold value that presets,, determine that then described current demand signal frame is the harmonic signal frame if surpass.
The harmonic signal pick-up unit that the embodiment of the invention provides comprises: acquiring unit is used for obtaining the frequency-region signal of current demand signal frame; Search unit is used to search for described frequency-region signal and obtains the first peak value information; Matching unit is used for described first peak value information and second peak information mated and obtains mating numerical value, and described second peak information is the peak information of the frequency-region signal in the signal frame that got access to before described current demand signal frame; Judging unit is used to judge whether described coupling numerical value surpasses the coupling threshold value that presets; Generation unit is used for determining that described current demand signal frame is the harmonic signal frame when the described coupling numerical value of described judgment unit judges surpasses the coupling threshold value that presets as a result.
The harmonic signal pick-up unit that the embodiment of the invention provides comprises: acquiring unit is used for obtaining the frequency-region signal of current demand signal frame; Phase place is eliminated the unit, is used for that described frequency-region signal is carried out phase place elimination processing and obtains the first pseudo-power spectrum signal; Computing unit, be used to calculate the cross-correlation coefficient between the described first pseudo-power spectrum signal and the second pseudo-power spectrum signal, the described second pseudo-power spectrum signal is eliminated the pseudo-power spectrum signal that obtains after the processing for the frequency-region signal in the signal frame that got access to carries out phase place before described current demand signal frame; Judging unit is used to judge whether described cross-correlation coefficient surpasses the coefficient threshold value that presets; Generation unit is used for determining that described current demand signal frame is the harmonic signal frame when the described cross-correlation coefficient of described judgment unit judges surpasses the coefficient threshold value that presets as a result.
As can be seen from the above technical solutions, the embodiment of the invention has the following advantages:
In the embodiment of the invention, mate to judge whether the current demand signal frame is the harmonic signal frame by peak information the peak information of current demand signal frame and signal frame before, because the harmonic signal frame has successional characteristics usually, correlativity between the harmonic signal frame of front and back is bigger, therefore utilizes the peak information of different signal frames to mate the accuracy that can improve judgement harmonic signal frame effectively;
In addition, when judging the harmonic signal frame, need not to calculate fundamental frequency in the embodiment of the invention, so can reduce the calculated amount in the harmonic signal frame deterministic process effectively;
In the embodiment of the invention, get access to after the frequency-region signal in the signal frame, this frequency-region signal is carried out phase place eliminate processing, thereby do not comprise phase information in the pseudo-power spectrum signal that obtains after feasible the processing, then calculate the current demand signal frame and can not be subjected to the influence of phase information during the cross-correlation coefficient between the front signal frame with it, thereby can judge the harmonic signal frame accurately.
Description of drawings
Fig. 1 is embodiment synoptic diagram of harmonic signal detection method in the embodiment of the invention;
Fig. 2 is another embodiment synoptic diagram of harmonic signal detection method in the embodiment of the invention;
Fig. 3 is a harmonic signal detection method embodiment synoptic diagram again in the embodiment of the invention;
Fig. 4 is the another embodiment synoptic diagram of harmonic signal detection method in the embodiment of the invention;
Fig. 5 is embodiment synoptic diagram of harmonic signal pick-up unit in the embodiment of the invention;
Fig. 6 is another embodiment synoptic diagram of harmonic signal pick-up unit in the embodiment of the invention;
Fig. 7 is a harmonic signal pick-up unit embodiment synoptic diagram again in the embodiment of the invention;
Fig. 8 is the another embodiment synoptic diagram of harmonic signal pick-up unit in the embodiment of the invention.
Embodiment
The embodiment of the invention provides a kind of harmonic signal detection method and device, is used to improve the accuracy that harmonic signal detects.
See also Fig. 1, embodiment of harmonic signal detection method comprises in the embodiment of the invention:
101, obtain frequency-region signal in the current demand signal frame;
In the present embodiment, be that example describes as the current demand signal frame, be understandable that in actual applications, the signal frame that needs can also be carried out the harmonic signal detection does not limit as the current demand signal frame herein with the current signal frame that is detecting.
102, the described frequency-region signal of search obtains the first peak value information;
In the present embodiment, after from the current demand signal frame, getting access to frequency-region signal, thereby then can search for this frequency-region signal and obtain the first peak value information, concrete can search for the peak information that this frequency-region signal obtains this frequency-region signal, as the first peak value information, this peak information specifically can comprise peak and/or peak amplitude, in actual applications with the peak information of this frequency-region signal, concrete peak information can also not limit for other attributes of peak value herein.
103, first peak value information and second peak information are mated obtain mating numerical value;
In the present embodiment, get access to after the first peak value information, this first peak value information and second peak information can be mated, thereby obtain mating numerical value, concrete matching process will be described in detail in subsequent embodiment.
Second peak information in the present embodiment is the peak information of the frequency-region signal in the signal frame that got access to before the current demand signal frame, be specifically as follows: the peak information of the frequency-region signal in the last signal frame of current demand signal frame, or the peak information of the frequency-region signal in preceding several signal frames of current demand signal frame, do not limit herein.
104, judge whether coupling numerical value surpasses the coupling threshold value that presets, if then execution in step 105, if not, then execution in step 106;
105, determine that the current demand signal frame is the harmonic signal frame;
If coupling numerical value surpasses the coupling threshold value that presets, illustrate that then the matching degree between current demand signal frame and the signal frame before is higher, can determine that the current demand signal frame is the harmonic signal frame.
106, determine that the current demand signal frame is not the harmonic signal frame.
If coupling numerical value is less than or equal to the coupling threshold value that presets, illustrate that then the matching degree between current demand signal frame and the signal frame before is lower, can determine that the current demand signal frame is not the harmonic signal frame.
In the present embodiment, mate to judge whether the current demand signal frame is the harmonic signal frame by peak information the peak information of current demand signal frame and signal frame before, because the harmonic signal frame has successional characteristics usually, correlativity between the harmonic signal frame of front and back is bigger, therefore utilizes the peak information of different signal frames to mate the accuracy that can improve judgement harmonic signal frame effectively;
In addition, when judging the harmonic signal frame, need not to calculate fundamental frequency in the present embodiment, so can reduce the calculated amount in the harmonic signal frame deterministic process effectively.
For ease of understanding, with an instantiation harmonic signal detection method in the present embodiment is described below, see also Fig. 2, specifically comprise:
201, obtain the voice time-domain signal;
In the present embodiment, the signal that pending harmonic signal detects can be the voice time-domain signal.
202, time-domain signal is converted to frequency-region signal;
Get access to after the voice time-domain signal, can carry out time-frequency to this time-domain signal and be converted to frequency-region signal, the mode of concrete time-frequency conversion can be correction type discrete cosine transform (MDCT, ModifiedDiscrete Cosine Transform), or discrete cosine transform (DCT, Discrete CosineTransform), or Fourier transform etc., do not limit herein.
203, frequency-region signal is carried out the peak value detection and obtain the first peak value information;
In the present embodiment, after getting access to the frequency-region signal frame, obtain the first peak value information thereby can carry out the peak value detection to frequency-region signal wherein, concrete peak information can comprise peak set and/or peak amplitude set.
The peak set is made up of several peaks, and the peak amplitude set is made up of several peak amplitude equally.
In the present embodiment, the mode that concrete peak value detects can for:
The peak and/or the peak amplitude of at least one point of amplitude maximum in the search frequency-region signal;
Deposit peak in first peak value location sets, and/or deposit peak amplitude the set of in first peak amplitude.
Need to prove that in the present embodiment, the number of the point of concrete required amplitude maximum of getting can be determined by the scope of required coupling in the practical application.
204, first peak value information and second peak information are mated obtain mating numerical value;
Get access to after the first peak value information, this first peak value information and second peak information can be mated, thereby obtain mating numerical value.
Wherein, second peak information is the peak information of the frequency-region signal in the signal frame that got access to before the current demand signal frame, be specifically as follows: the peak information of the frequency-region signal in the last signal frame of current demand signal frame, or the peak information of the frequency-region signal in preceding several signal frames of current demand signal frame, do not limit herein.
In the present embodiment, concrete matching process can for:
To satisfy the number of times of prerequisite between the peak in peak in the first peak value location sets and the set of second peak as coupling numerical value;
And/or,
To satisfy the number of times of prerequisite between the peak amplitude in peak amplitude in the set of first peak amplitude and the set of second peak amplitude as coupling numerical value.
In the present embodiment promptly, the coupling of peak and the coupling of peak amplitude can only be used wherein a kind of in actual applications, or dual mode can be used in combination.
The above-mentioned prerequisite of having mentioned, in the present embodiment, for the coupling of peak, concrete prerequisite can for:
Arbitrary peak in the set of arbitrary peak in the first peak value location sets and second peak equates, then determines to satisfy prerequisite, occurs equal situation each time, and then this coupling numerical value adds 1, until having mated all peaks;
Or,
The difference of the arbitrary peak in arbitrary peak in the first peak value location sets and the set of second peak is less than the difference threshold value that presets, then determine to satisfy prerequisite, the situation of difference less than the difference threshold value that presets appears each time, then this coupling numerical value adds 1, until having mated all peaks.
Matching process for the matching process of peak amplitude and peak is similar, repeats no more herein.
205, judge whether coupling numerical value surpasses the coupling threshold value, if then execution in step 207, if not, then execution in step 206;
Coupling threshold value in the present embodiment can determine according to actual conditions, if examination criteria is strict more, then this coupling threshold value can be provided with big more.
Need to prove, if several signal frames of choosing before the current demand signal frame compare, then criterion can be decided according to actual conditions, coupling numerical value between the peak information of each signal frame before that for example can require the peak information of current demand signal frame and select is all above coupling threshold value ability execution in step 207, only need perhaps also can to require coupling numerical value between the peak information of the peak information of current demand signal frame and a last signal frame to surpass the coupling threshold value can execution in step 207, perhaps also can require the current demand signal frame peak information and selected before execution in step 207 when surpassing the coupling threshold value above half coupling numerical value is arranged in the plurality of matching numerical value between the peak information of each signal frame, concrete criterion does not limit herein.
In the present embodiment,, can choose current demand signal frame harmonic signal frame before object as a comparison for further improving the accuracy that harmonic signal detects.
206, determine that the current demand signal frame is not the harmonic signal frame, and process ends;
If coupling numerical value does not surpass the coupling threshold value, illustrate that promptly the matching degree between current demand signal frame and the signal frame before is lower, determine that then the current demand signal frame is not the harmonic signal frame.
207, determine that the current demand signal frame is the harmonic signal frame, and execution in step 208;
If coupling numerical value surpasses the coupling threshold value, illustrate that promptly the matching degree between current demand signal frame and the signal frame before is higher, determine that then the current demand signal frame is the harmonic signal frame.
208, peak information is upgraded and carries out the testing process of next signal frame.
In the present embodiment, after definite current demand signal frame is the harmonic signal frame, then can replaces the second original peak information, and next signal frame is detected with the peak information of current demand signal frame, concrete testing process and aforementioned testing process are similar, repeat no more herein.
Need to prove, in the present embodiment, between step 202 and the step 203, promptly get access to after the frequency-region signal, can carry out phase place elimination processing to this frequency-region signal and obtain the first pseudo-power spectrum signal, simultaneously the frequency-region signal of the signal frame before the current demand signal frame being carried out phase place equally eliminates to handle and obtains the second pseudo-power spectrum signal, then step 203 can be carried out peak value to the first pseudo-power spectrum signal and detected and to obtain the first peak value information, step 204 can be mated second peak information that this first peak value information and the second pseudo-power spectrum signal carry out obtaining after peak value detects, the process of concrete coupling is similar, repeats no more herein.
In the present embodiment, mate to judge whether the current demand signal frame is the harmonic signal frame by peak information the peak information of current demand signal frame and signal frame before, because the harmonic signal frame has successional characteristics usually, correlativity between the harmonic signal frame of front and back is bigger, therefore utilizes the peak information of different signal frames to mate the accuracy that can improve judgement harmonic signal frame effectively;
Secondly, when judging the harmonic signal frame, need not to calculate fundamental frequency in the present embodiment, so can reduce the calculated amount in the harmonic signal frame deterministic process effectively;
Once more, the harmonic signal frame that can choose in the present embodiment before the current demand signal frame compares object, thereby can improve the accuracy that the harmonic signal to the current demand signal frame detects;
Further, in the present embodiment, after getting access to frequency-region signal, can carry out phase place to this frequency-region signal and eliminate processing, and the signal that phase place is eliminated after handling carries out peak value detection and coupling, therefore can eliminate the influence that phase place detects harmonic signal, thereby further improve accuracy the harmonic signal detection of current demand signal frame.
Foregoing description the harmonic signal detection method embodiment in the embodiment of the invention, for ease of understanding, with several concrete examples above-mentioned flow process is described in detail below:
One, the absolute matching way of peak:
Process concrete in the present embodiment comprises:
(1) voice time-domain signal frame X comprises T sampled point: X={x (t), t=0 ..., T-1} is converted into the frequency-region signal frame Y={y (i) of I sampled point by the MDCT conversion, i=0 ..., I-1};
(2) signal Y is carried out peak value searching: compare I y (i), find maximum amplitude y (j) wherein, write down the position j at its place, and this peak j put into variable Position (0), then y (j) is changed to a minimum constant C, searches out once more when preventing from for the second time to search for maximum amplitude.
Repeat above-mentioned search procedure afterwards, until obtaining N peak (N>0), and note they position Position (1), Position (2) ..., Position (N-1).
N can get 10 in the present embodiment, and C can get-1E10, and concrete N and the value of C can determine according to practical application, do not limit herein.
Concrete search procedure can comprise following program:
For n=0 to N-1 // N search of execution
Position (n)=j; // deposit j in Position (n)
Y (j)=C; // and y (j) is set to minimal value
end
end
Can choose the peak of N peak value by above-mentioned search.
(3) N the peak PrePosition (m) with N peak Position (n) signal frame last with it that search compares, and notes the number matchPos that both equate.
Last signal frame peak: PrePosition (m) m=0,1 ..., N-1.
Concrete matching process is as follows:
For n=0 to N-1 // N coupling of execution
for?m=0?to?N-1
If Position (n) equals PrePosition (m)
MatchPos=matchPos+1; // as Position (n) and PrePosition (m) when equating, matchPos is added 1
end
end
end
(4) comparison match is counted matchPos and the position threshold ThrePos that pre-sets, if matchPos is greater than ThrePos, determine that then the current demand signal frame is the harmonic signal frame, ThrePos in the present embodiment can determine according to actual conditions, if it is strict more to the detection of harmonic signal, what then the numerical value of this ThrePos can be provided with is big more, and the ThrePos in the present embodiment can be set to 8.
Need to prove, in the present embodiment, if matchPos, determines then that the current demand signal frame is the harmonic signal frame greater than ThrePos, if matchPos is less than ThrePos, whether determine that then the current demand signal frame is not the harmonic signal frame, if matchPos equals ThrePos, can be the harmonic signal frame according to concrete mode decision scheme decision current demand signal frame then, in the present embodiment, if matchPos equals ThrePos, determine that then the current demand signal frame is the harmonic signal frame, other mode decision scheme does not limit herein.
(5) upgrade peak information.
In the present embodiment, to the current demand signal frame detect finish after, can upgrade peak information according to testing result so that carry out the detection of next signal frame, concrete renewal process can for:
Upgrade the peak information PrePosition (n) of front signal frame with Position (n):
for?n=0?to?N-1
PrePosition(n)=Position(n);
end
Two, the relative matching way of peak:
The processing procedure of the detailed process in the present embodiment and the absolute matching way of aforementioned peak is similar, difference only is in the previous embodiment it is when Position (n) equals PrePosition (m), matchPos is added 1, and can be Position (n) in the present embodiment with the absolute value of the difference of PrePosition (m) less than certain difference threshold value, then matchPos is added 1, concrete difference threshold value delta can determine that the delta in the present embodiment can be 0.5 according to actual conditions.
The detailed process and the previous embodiment of mating in the present embodiment are similar, repeat no more herein.
Three, peak information combinations matches mode:
Process concrete in the present embodiment comprises:
(1) voice time-domain signal frame X comprises T sampled point: X={x (t), t=0 ..., T-1} is converted into the frequency-region signal frame Y={y (i) of I sampled point by the MDCT conversion, i=0 ..., I-1};
(2) signal Y is carried out peak value searching: compare I y (i), find maximum amplitude y (j) wherein, write down its amplitude, put into variable V alue (0), write down the position j at its place, and this peak j put into variable Position (0), and then y (j) is changed to a minimum constant C, search out once more when preventing to search for maximum amplitude for the second time.
Repeat above-mentioned search procedure afterwards, until obtaining N peak amplitude and peak (N>0), and note their amplitude and position [Position (1), Value (1)], [Position (2), Value (2)] ..., [Position (N-1), Value (N-1)].
N can get 10 in the present embodiment, and C can get-1E10, and concrete N and the value of C can determine according to practical application, do not limit herein.
Concrete search procedure can comprise following program:
For n=0 to N-1 // N search of execution
Position (n)=j; // deposit j in Position (n)
Value (n)=y (j); // deposit y (j) in Value (n)
Y (j)=C; // and y (j) is set to minimal value
end
end
Can choose the peak and the peak amplitude of N peak value by above-mentioned search.
(3) N peak information with Position (n), Value (n) two signal frames preceding with it compares respectively, and record coupling numerical value.
In the present embodiment, being that example describes with preceding two signal frames of current demand signal frame, being understandable that, in actual applications, can be previous signal frame only equally, perhaps can not limit for more a plurality of signal frames of front herein.
Last signal frame peak information: PrePosition1 (m), PreValue1 (m) m=0,1 ..., N-1
Preceding secondary signal frame peak information: PrePosition2 (m), PreValue2 (m) m=0,1 ..., N-1
Concrete matching process is as follows:
for?n=0?to?N-1
for?m=0?to?N-1
If Position (n) equals PrePosition1 (m)
MatchPos1=matchPos1+1; // as Position (n) and PrePosition1 (m) when equating, matchPos1 is added 1
end
If Position (n) equals PrePosition2 (m)
MatchPos2=matchPos2+1; // as Position (n) and PrePosition2 (m) when equating, matchPos2 is added 1
end
if|Value(n)-PreValue1(m)|≤beta
MatchVal1=matchVal1+1; // when the absolute value of the difference of Value (n) and PreValue1 (m) during, matchVal1 is added 1 less than beta
end
if|Value(n)-PreValue2(m)|≤beta
MatchVal2=matchVal2+1; // when the absolute value of the difference of Value (n) and PreValue2 (m) during, matchVal2 is added 1 less than beta
end
end
end
Wherein, the numerical value of beta can determine according to practical application, beta 〉=0, and the beta in the present embodiment can equal 0.5.
Need to prove, in the present embodiment, what peak adopted is absolute coupling, what peak amplitude adopted is relative coupling, this mainly is that the numerical value of considering amplitude tends to more accurate, the mode of absolute coupling may be difficult to match identical peak amplitude, but in actual applications, peak amplitude equally also can adopt absolute coupling, peak equally also can adopt relative coupling, or all adopt absolute coupling, or all adopt coupling relatively, concrete matching way does not limit herein.
(4) comparison match is counted matchPos1, matchPos2, matchVal1, matchVal2 and position threshold ThrePos, amplitude thresholds ThreVal.
In the present embodiment, when matchPos1 more than or equal to ThrePos, and matchPos2 is more than or equal to ThrePos, and matchVal1 is more than or equal to ThreVal, and matchVal2 determines that just the current demand signal frame is the harmonic signal frame during more than or equal to ThreVal.
Need to prove, in actual applications, when not necessarily needing all to set up, 4 above-mentioned comparisons could determine that the current demand signal frame is the harmonic signal frame, equally can be only determine that the current demand signal frame is the harmonic signal frame when one of them or several are set up, concrete decision rule does not limit herein.
(5) upgrade peak information.
If the current demand signal frame is the harmonic signal frame, then use Position (n), Value (n) to upgrade the peak information of front signal frame:
Concrete renewal process can for:
If signal Y is a harmonic signal
for?n=0?to?N-1
PrePosition2(n)=PrePosition1(n);
PrePosition1(n)=PrePosition(n);
PreValue2(n)=PreValue1(n);
PreValue1(n)=PreValue(n);
end
end
In each above-mentioned embodiment, mate to judge whether the current demand signal frame is the harmonic signal frame by peak information the peak information of current demand signal frame and signal frame before, because the harmonic signal frame has successional characteristics usually, correlativity between the harmonic signal frame of front and back is bigger, therefore utilizes the peak information of different signal frames to mate the accuracy that can improve judgement harmonic signal frame effectively;
Secondly, when judging the harmonic signal frame, need not to calculate fundamental frequency in the present embodiment, so can reduce the calculated amount in the harmonic signal frame deterministic process effectively.
See also Fig. 3, introduce another embodiment of harmonic signal detection method in the embodiment of the invention below, specifically comprise:
301, obtain frequency-region signal in the current demand signal frame;
In the present embodiment, be that example describes as the current demand signal frame, be understandable that in actual applications, the signal frame that needs can also be carried out the harmonic signal detection does not limit as the current demand signal frame herein with the current signal frame that is detecting.
302, frequency-region signal is carried out phase place elimination processing and obtain the first pseudo-power spectrum signal;
Frequency-region signal in the present embodiment includes amplitude and phase information, and this frequency-region signal is carried out can obtaining the first pseudo-power spectrum signal after the phase place elimination processing, does not comprise phase information in this first pseudo-power spectrum signal.
303, the cross-correlation coefficient between the calculating first pseudo-power spectrum signal and the second pseudo-power spectrum signal;
In the present embodiment, get access to after the first pseudo-power spectrum signal, can calculate the cross-correlation coefficient between this first pseudo-power spectrum signal and the second pseudo-power spectrum signal.
The second pseudo-power spectrum signal in the present embodiment is eliminated the pseudo-power spectrum signal that obtains after the processing for the frequency-region signal in the signal frame that got access to carries out phase place before described current demand signal frame, be specifically as follows: the frequency-region signal in the last signal frame of current demand signal frame is eliminated through phase place and is handled the pseudo-power spectrum signal that obtains, or the frequency-region signal in preceding several signal frames of current demand signal frame does not limit through the pseudo-power spectrum signal that phase place elimination processing obtains herein.
304, judge whether cross-correlation coefficient surpasses the coefficient threshold value that presets, if surpass, then execution in step 305, if do not surpass, then execution in step 306;
305, determine that the current demand signal frame is the harmonic signal frame;
If cross-correlation coefficient surpasses the coefficient threshold value that presets, illustrate that then the degree of correlation between current demand signal frame and the signal frame before is higher, can determine that the current demand signal frame is the harmonic signal frame.
306, determine that the current demand signal frame is not the harmonic signal frame.
If cross-correlation coefficient does not surpass the coefficient threshold value that presets, illustrate that then the degree of correlation between current demand signal frame and the signal frame before is lower, can determine that the current demand signal frame is not the harmonic signal frame.
In the present embodiment, get access to after the frequency-region signal in the signal frame, this frequency-region signal is carried out phase place eliminate processing, thereby do not comprise phase information in the pseudo-power spectrum signal that obtains after feasible the processing, then calculate the current demand signal frame and can not be subjected to the influence of phase information during the cross-correlation coefficient between the front signal frame with it, thereby can judge the harmonic signal frame accurately.
For ease of understanding, with an instantiation harmonic signal detection method in the present embodiment is described below, see also Fig. 4, the another embodiment of harmonic signal detection method in the embodiment of the invention comprises:
401, obtain the voice time-domain signal;
In the present embodiment, the signal that pending harmonic signal detects can be the voice time-domain signal.
402, time-domain signal is converted to frequency-region signal;
Get access to after the voice time-domain signal, can carry out time-frequency to this time-domain signal and be converted to frequency-region signal, the mode of concrete time-frequency conversion can be MDCT, or DCT, or Fourier transform etc., do not limit herein.
403, frequency-region signal is carried out phase place elimination processing and obtain the first pseudo-power spectrum signal;
In the present embodiment, after getting access to the frequency-region signal frame, obtain the first pseudo-power spectrum signal thereby can carry out phase place elimination processing to frequency-region signal wherein.
404, the cross-correlation coefficient between the calculating first pseudo-power spectrum signal and the second pseudo-power spectrum signal;
Get access to after the first pseudo-power spectrum signal, can calculate the cross-correlation coefficient between this first pseudo-power spectrum signal and the second pseudo-power spectrum signal.
Wherein, the second pseudo-power spectrum signal is that the frequency-region signal process phase place in the signal frame that got access to before the current demand signal frame is eliminated the pseudo-power spectrum signal after handling, be specifically as follows: the frequency-region signal process phase place in the last signal frame of current demand signal frame is eliminated the pseudo-power spectrum signal after handling, or the pseudo-power spectrum signal after the process of the frequency-region signal in preceding several signal frames of current demand signal frame the phase place elimination processing, do not limit herein.
The process of calculating cross-correlation coefficient in the present embodiment is those skilled in the art's a common practise, does not limit herein.
405, judge whether cross-correlation coefficient surpasses the coefficient threshold value, if then execution in step 407, if not, then execution in step 406;
Coefficient threshold value in the present embodiment can determine according to actual conditions, if examination criteria is strict more, then this coefficient threshold value can be provided with big more.
Need to prove, if several signal frames of choosing before the current demand signal frame compare, then criterion can be decided according to actual conditions, for example can require after the current demand signal frame eliminate to be handled through phase place pseudo-power spectrum signal with selected before each signal frame eliminate through phase place that the cross correlation number average surpasses coefficient threshold value ability execution in step 407 between the pseudo-power spectrum signal after handling, cross-correlation coefficient between the pseudo-power spectrum signal after pseudo-power spectrum signal after perhaps also can requiring only to need the current demand signal frame through phase place elimination processing and the elimination of last signal frame process phase place are handled gets final product execution in step 407 above the coefficient threshold value, perhaps also can require after the current demand signal frame eliminate to be handled through phase place pseudo-power spectrum signal with selected before each signal frame eliminate execution in step 407 when having above half cross-correlation coefficient in several cross-correlation coefficients between the pseudo-power spectrum signal after handling through phase place above the coefficient threshold value, concrete criterion does not limit herein.
In the present embodiment,, can choose current demand signal frame harmonic signal frame before object as a comparison for further improving the accuracy that harmonic signal detects.
406, determine that the current demand signal frame is not the harmonic signal frame, and process ends.
If cross-correlation coefficient does not surpass the coefficient threshold value that presets, illustrate that then the degree of correlation between current demand signal frame and the signal frame before is lower, can determine that the current demand signal frame is not the harmonic signal frame.
407, determine that the current demand signal frame is the harmonic signal frame.
If cross-correlation coefficient surpasses the coefficient threshold value that presets, illustrate that then the degree of correlation between current demand signal frame and the signal frame before is higher, can determine that the current demand signal frame is the harmonic signal frame.
In the present embodiment, get access to after the frequency-region signal in the signal frame, this frequency-region signal is carried out phase place eliminate processing, thereby do not comprise phase information in the pseudo-power spectrum signal that obtains after feasible the processing, then calculate the current demand signal frame and can not be subjected to the influence of phase information during the cross-correlation coefficient between the front signal frame with it, thereby can judge the harmonic signal frame accurately.
Foregoing description the flow process of a kind of harmonic signal detection method of providing of the embodiment of the invention, for ease of understanding, with a concrete example above-mentioned flow process is described in detail below:
Process concrete in the present embodiment comprises:
(1) voice time-domain signal frame X comprises T sampled point: X={x (t), t=0 ..., T-1} is converted into the frequency-region signal frame Y={y (i) of I sampled point by the MDCT conversion, i=0 ..., I-1};
(2) MDCT coefficient Y (being frequency-region signal) is calculated pseudo-power spectrum signal Pseudo_power (i):
(3) calculate the pseudo-power spectrum signal PrePseudo_power (i) of last signal frame and the cross-correlation coefficient CorrXY of the pseudo-power spectrum signal Pseudo_power of current demand signal frame (i):
In the present embodiment be last signal frame with the current demand signal frame as calculating object, be understandable that in actual applications, can also adopt several signal frames before the current demand signal frame as calculating object equally, concrete account form is similar.
(4) relation of comparison cross-correlation coefficient CorrXY and threshold value ThreCorr, if CorrXY is greater than ThreCorr, determine that then the current demand signal frame is the harmonic signal frame, ThreCorr in the present embodiment can determine according to actual conditions, if it is strict more to the detection of harmonic signal, what then the numerical value of this ThreCorr can be provided with is big more, and the ThreCorr in the present embodiment can be set to 0.7.
Need to prove, in the present embodiment, if CorrXY, determines then that the current demand signal frame is the harmonic signal frame greater than ThreCorr, if CorrXY is less than ThreCorr, whether determine that then the current demand signal frame is not the harmonic signal frame, if CorrXY equals ThreCorr, can be the harmonic signal frame according to concrete mode decision scheme decision current demand signal frame then, in the present embodiment, if CorrXY equals ThreCorr, determine that then the current demand signal frame is the harmonic signal frame, other mode decision scheme does not limit herein.
In the present embodiment, get access to after the frequency-region signal in the signal frame, this frequency-region signal is carried out phase place eliminate processing, thereby do not comprise phase information in the pseudo-power spectrum signal that obtains after feasible the processing, then calculate the current demand signal frame and can not be subjected to the influence of phase information during the cross-correlation coefficient between the front signal frame with it, thereby can judge the harmonic signal frame accurately.
Introduce the harmonic signal pick-up unit in the embodiment of the invention below, see also Fig. 5, embodiment of harmonic signal pick-up unit comprises in the embodiment of the invention:
Acquiring unit 501 is used for obtaining the frequency-region signal of current demand signal frame;
Concrete, the search unit 502 in the present embodiment can be searched for the peak information that frequency-region signal obtains this frequency-region signal, with the peak information of this frequency-region signal as the first peak value information.
Judging unit 504 is used to judge whether the coupling numerical value that matching unit 503 obtains surpasses the coupling threshold value that presets;
Below the harmonic signal pick-up unit in the embodiment of the invention is described in detail, see also Fig. 6, another embodiment of harmonic signal pick-up unit comprises in the embodiment of the invention:
Acquiring unit 601 is used for obtaining the frequency-region signal of current demand signal frame;
Phase place is eliminated unit 602, and the frequency-region signal that is used for current demand signal frame that acquiring unit 601 is got access to carries out phase place to be eliminated to handle and obtain the first pseudo-power spectrum signal;
Concrete, in the present embodiment, search unit 603 can be searched for phase place and eliminate unit 602 and handle the peak information that the first pseudo-power spectrum signal later obtains this first pseudo-power spectrum signal, with the peak information of this first pseudo-power spectrum signal as the first peak value information.
The peak information of frequency-region signal that can choose the harmonic signal frame before the current demand signal frame in the present embodiment is as match objects.
The process of mating among process that matching unit 604 specifically mates in the present embodiment and the preceding method embodiment is similar, repeats no more herein.
Judging unit 605 is used to judge whether the coupling numerical value that matching unit 604 obtains surpasses the coupling threshold value that presets;
Foundation and the process judged among the foundation that judging unit 605 is specifically judged in the present embodiment and process and the preceding method embodiment are the same, repeat no more herein.
In the present embodiment, the first peak value information comprises the first peak value location sets and/or the set of first peak amplitude, and second peak information comprises set of second peak and/or the set of second peak amplitude;
The first peak value location sets comprises at least one peak, and the set of first peak amplitude comprises at least one peak amplitude, and the set of second peak comprises at least one peak, and the set of second peak amplitude comprises at least one peak amplitude;
And/or,
In the present embodiment, matching unit 604 can mate the peak information of the peak information of current demand signal frame and signal frame before, then judging unit 605 and as a result generation unit 606 can judge whether the current demand signal frame is the harmonic signal frame, because the harmonic signal frame has successional characteristics usually, correlativity between the harmonic signal frame of front and back is bigger, therefore utilizes the peak information of different signal frames to mate the accuracy that can improve judgement harmonic signal frame effectively;
Secondly, in the present embodiment, when judging the harmonic signal frame, need not to calculate fundamental frequency, so can reduce the calculated amount in the harmonic signal frame deterministic process effectively;
Once more, in the present embodiment, the harmonic signal frame that matching unit 604 can be chosen before the current demand signal frame compares object, thereby can improve the accuracy to the harmonic signal detection of current demand signal frame;
Further, in the present embodiment, after acquiring unit 601 gets access to frequency-region signal, phase place is eliminated unit 602 can carry out phase place elimination processing to this frequency-region signal, the peak information of the signal after then 604 pairs of phase places eliminations of matching unit are handled mates, therefore can eliminate the influence that phase place detects harmonic signal, thereby further improve accuracy the harmonic signal detection of current demand signal frame.
Introduce the another kind of harmonic signal pick-up unit in the embodiment of the invention below, see also Fig. 7, in the embodiment of the invention harmonic signal pick-up unit again an embodiment comprise:
Acquiring unit 701 is used for obtaining the frequency-region signal of current demand signal frame;
Phase place is eliminated unit 702, and the frequency-region signal that is used for that acquiring unit 701 is got access to carries out phase place elimination processing and obtains the first pseudo-power spectrum signal;
Computing unit 703, be used to calculate phase place and eliminate the first pseudo-power spectrum signal that unit 702 obtains and the cross-correlation coefficient between the second pseudo-power spectrum signal, this second pseudo-power spectrum signal is eliminated the pseudo-power spectrum signal that obtains after the processing for the frequency-region signal in the signal frame that got access to carries out phase place before described current demand signal frame, be specifically as follows: the frequency-region signal in the last signal frame of current demand signal frame is eliminated through phase place and is handled the pseudo-power spectrum signal that obtains, or the frequency-region signal in preceding several signal frames of current demand signal frame does not limit through the pseudo-power spectrum signal that phase place elimination processing obtains herein.
Judging unit 704 is used to judge whether the cross-correlation coefficient that computing unit 703 calculates surpasses the coefficient threshold value that presets;
Generation unit 705 as a result, are used for determining that the current demand signal frame is the harmonic signal frame when judging unit 704 judges that cross-correlation coefficient surpasses the coefficient threshold value that presets.
Below the harmonic signal pick-up unit in the embodiment of the invention is described in detail, see also Fig. 8, the another embodiment of harmonic signal pick-up unit comprises in the embodiment of the invention:
Acquiring unit 802 is used for obtaining the frequency-region signal of current demand signal frame;
Phase place is eliminated unit 803, and the frequency-region signal that is used for that acquiring unit 802 is got access to carries out phase place elimination processing and obtains the first pseudo-power spectrum signal;
Phase place elimination unit 803 can carry out phase place elimination processing to frequency-region signal by following mode in the present embodiment:
Frequency-region signal is calculated pseudo-power spectrum signal Pseudo_power (i):
Frequency-region signal is Y={y (i), i=0 ..., I-1}, wherein I is the frequency domain sample point; The first pseudo-power spectrum signal is Pseudo_power (i).
Calculate the pseudo-power spectrum signal PrePseudo_power (i) of last signal frame and the cross-correlation coefficient CorrXY of the pseudo-power spectrum signal Pseudo_power of current demand signal frame (i):
Wherein, PrePseudo_power (i) is the pseudo-power spectrum signal (i.e. the second pseudo-power spectrum signal) of last signal frame, and Pseudo_power (i) is the pseudo-power spectrum signal (i.e. the first pseudo-power spectrum signal) of current demand signal frame, and I is the frequency domain sample point.
Judging unit 805 is used to judge whether the cross-correlation coefficient that computing unit 804 calculates surpasses the coefficient threshold value that presets;
In the present embodiment, acquiring unit 802 gets access to after the frequency-region signal in the signal frame, phase place is eliminated 803 pairs of these frequency-region signals in unit and is carried out phase place elimination processing, thereby do not comprise phase information in the pseudo-power spectrum signal that obtains after feasible the processing, then computing unit 804 calculates the current demand signal frames and can not be subjected to the influence of phase information during the cross-correlation coefficient between the front signal frame with it, thereby can judge the harmonic signal frame accurately.
One of ordinary skill in the art will appreciate that all or part of step that realizes in the foregoing description method is to instruct relevant hardware to finish by program, described program can be stored in a kind of computer-readable recording medium, this program comprises the steps: when carrying out
Obtain the frequency-region signal in the current demand signal frame;
Search for described frequency-region signal and obtain the first peak value information;
Described first peak value information and second peak information mated obtain mating numerical value, described second peak information is the peak information of the frequency-region signal in the signal frame that got access to before described current demand signal frame;
Judge whether described coupling numerical value surpasses the coupling threshold value that presets,, determine that then described current demand signal frame is the harmonic signal frame if surpass.
The above-mentioned storage medium of mentioning can be a ROM (read-only memory), disk or CD etc.
More than a kind of harmonic signal detection method provided by the present invention and device are described in detail, for one of ordinary skill in the art, thought according to the embodiment of the invention, part in specific embodiments and applications all can change, in sum, this description should not be construed as limitation of the present invention.
Claims (17)
1. a harmonic signal detection method is characterized in that, comprising:
Obtain the frequency-region signal in the current demand signal frame;
Search for described frequency-region signal and obtain the first peak value information;
Described first peak value information and second peak information mated obtain mating numerical value, described second peak information is the peak information of the frequency-region signal in the signal frame that got access to before described current demand signal frame;
Judge whether described coupling numerical value surpasses the coupling threshold value that presets,, determine that then described current demand signal frame is the harmonic signal frame if surpass.
2. method according to claim 1 is characterized in that, described search frequency-region signal obtains the first peak value information and comprises:
Search for the peak information of described frequency-region signal, with the peak information of described frequency-region signal as described first peak value information.
3. method according to claim 2 is characterized in that,
Described first peak value information comprises the first peak value location sets and/or the set of first peak amplitude, and described second peak information comprises set of second peak and/or the set of second peak amplitude;
The described first peak value location sets comprises at least one peak, described first peak amplitude set comprises at least one peak amplitude, described second peak set comprises at least one peak, and described second peak amplitude set comprises at least one peak amplitude;
Described first peak value information and second peak information are mated, obtain mating numerical value and comprise:
To satisfy the number of times of prerequisite as described coupling numerical value between the peak in peak in the described first peak value location sets and the set of described second peak;
And/or,
To satisfy the number of times of prerequisite as described coupling numerical value between peak amplitude in described first peak amplitude set and the peak amplitude in the set of described second peak amplitude.
4. method according to claim 3 is characterized in that, judges whether satisfy prerequisite between the peak in the set of the peak in the first peak value location sets and second peak in the following way:
When the arbitrary peak in the described first peak value location sets when arbitrary peak in the set of described second peak equates, then determine to satisfy prerequisite;
Or,
When the difference of the arbitrary peak in the set of the arbitrary peak in the described first peak value location sets and described second peak during, then determine to satisfy prerequisite less than the difference threshold value that presets.
5. according to claim 3 or 4 described methods, it is characterized in that judge whether satisfy prerequisite between peak amplitude in first peak amplitude set peak amplitude in gathering with second peak amplitude in the following way:
When the arbitrary peak amplitude in the set of described first peak amplitude when arbitrary peak amplitude during described second peak amplitude is gathered equates, then determine to satisfy prerequisite;
Or,
The difference of the arbitrary peak amplitude in arbitrary peak amplitude in the set of described first peak amplitude and described second peak amplitude are gathered then determines to satisfy prerequisite during less than the difference threshold value that presets.
6. according to claim 3 or 4 described methods, it is characterized in that the peak information of described search frequency-region signal comprises:
Search for the peak and/or the peak amplitude of at least one point of amplitude maximum in the described frequency-region signal;
Deposit described peak in the described first peak value location sets, and/or deposit described peak amplitude in described first peak amplitude set.
7. according to each described method in the claim 1 to 4, it is characterized in that described second peak information is the peak information of the frequency-region signal at least one the harmonic signal frame that gets access to before the described current demand signal frame.
8. method according to claim 1 is characterized in that, described search frequency-region signal obtains the first peak value information and comprises:
Described frequency-region signal is carried out phase place elimination processing obtain the first pseudo-power spectrum signal;
Search for the peak information of the described first pseudo-power spectrum signal;
With the peak information of the described first pseudo-power spectrum signal as described first peak value information;
Described second peak information is that the frequency-region signal at least one the harmonic signal frame that gets access to before the described current demand signal frame carries out phase place and eliminates the peak information of handling the second pseudo-power spectrum signal that obtains.
9. method according to claim 8 is characterized in that, describedly frequency-region signal is carried out phase place eliminate to handle and to obtain the first pseudo-power spectrum signal and specifically realize in the following way:
Described frequency-region signal is Y={y (i), i=0 ..., I-1}, wherein I is the frequency domain sample point; The first pseudo-power spectrum signal is Pseudo_power (i).
10. a harmonic signal detection method is characterized in that, comprising:
Obtain the frequency-region signal in the current demand signal frame;
Described frequency-region signal is carried out phase place elimination processing obtain the first pseudo-power spectrum signal;
Calculate the cross-correlation coefficient between the described first pseudo-power spectrum signal and the second pseudo-power spectrum signal, the described second pseudo-power spectrum signal is eliminated the pseudo-power spectrum signal that obtains after the processing for the frequency-region signal in the signal frame that got access to carries out phase place before described current demand signal frame;
Judge whether described cross-correlation coefficient surpasses the coefficient threshold value that presets,, determine that then described current demand signal frame is the harmonic signal frame if surpass.
11. method according to claim 10 is characterized in that, describedly frequency-region signal is carried out phase place eliminate to handle and to obtain the first pseudo-power spectrum signal and specifically realize in the following way:
Described frequency-region signal is Y={y (i), i=0 ..., I-1}, wherein I is the frequency domain sample point; The first pseudo-power spectrum signal is Pseudo_power (i).
12. method according to claim 11 is characterized in that, the cross-correlation coefficient between the described first pseudo-power spectrum signal and the second pseudo-power spectrum signal calculates in the following way:
The described second pseudo-power spectrum signal is PrePseudo_power (i), and the described first pseudo-power spectrum signal is Pseudo_power (i), and described cross-correlation coefficient is CorrXY.
13., it is characterized in that the described second pseudo-power spectrum signal is that at least one harmonic signal frame that described current demand signal frame gets access to before carries out the pseudo-power spectrum signal that phase place elimination processing obtains according to each described method in the claim 10 to 12.
14. a harmonic signal pick-up unit is characterized in that, comprising:
Acquiring unit is used for obtaining the frequency-region signal of current demand signal frame;
Search unit is used to search for described frequency-region signal and obtains the first peak value information;
Matching unit is used for described first peak value information and second peak information mated and obtains mating numerical value, and described second peak information is the peak information of the frequency-region signal in the signal frame that got access to before described current demand signal frame;
Judging unit is used to judge whether described coupling numerical value surpasses the coupling threshold value that presets;
Generation unit is used for determining that described current demand signal frame is the harmonic signal frame when the described coupling numerical value of described judgment unit judges surpasses the coupling threshold value that presets as a result.
15. harmonic signal pick-up unit according to claim 14, it is characterized in that, described first peak value information comprises the first peak value location sets and/or the set of first peak amplitude, and described second peak information comprises set of second peak and/or the set of second peak amplitude;
The described first peak value location sets comprises at least one peak, described first peak amplitude set comprises at least one peak amplitude, described second peak set comprises at least one peak, and described second peak amplitude set comprises at least one peak amplitude;
Described matching unit comprises:
First matching unit is used for and will satisfies the number of times of prerequisite as described coupling numerical value between the peak in the peak of the described first peak value location sets and the set of described second peak;
And/or,
Second matching unit satisfies the number of times of prerequisite as described coupling numerical value between the peak amplitude in being used for the peak amplitude of described first peak amplitude set gathered with described second peak amplitude.
16., it is characterized in that described harmonic wave pick-up unit also comprises according to claim 14 or 15 described harmonic signal pick-up units:
Phase place is eliminated the unit, and the frequency-region signal that is used for current demand signal frame that described acquiring unit is got access to carries out phase place to be eliminated to handle and obtain the first pseudo-power spectrum signal;
Described search unit is searched for the peak information of the described first pseudo-power spectrum signal, with described peak information as described first peak value information.
17. a harmonic signal pick-up unit is characterized in that, comprising:
Acquiring unit is used for obtaining the frequency-region signal of current demand signal frame;
Phase place is eliminated the unit, is used for that described frequency-region signal is carried out phase place elimination processing and obtains the first pseudo-power spectrum signal;
Computing unit, be used to calculate the cross-correlation coefficient between the described first pseudo-power spectrum signal and the second pseudo-power spectrum signal, the described second pseudo-power spectrum signal is eliminated the pseudo-power spectrum signal that obtains after the processing for the frequency-region signal in the signal frame that got access to carries out phase place before described current demand signal frame;
Judging unit is used to judge whether described cross-correlation coefficient surpasses the coefficient threshold value that presets;
Generation unit is used for determining that described current demand signal frame is the harmonic signal frame when the described cross-correlation coefficient of described judgment unit judges surpasses the coefficient threshold value that presets as a result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910134924XA CN101858939B (en) | 2009-04-10 | 2009-04-10 | Method and device for detecting harmonic signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910134924XA CN101858939B (en) | 2009-04-10 | 2009-04-10 | Method and device for detecting harmonic signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101858939A true CN101858939A (en) | 2010-10-13 |
CN101858939B CN101858939B (en) | 2012-11-07 |
Family
ID=42944957
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910134924XA Expired - Fee Related CN101858939B (en) | 2009-04-10 | 2009-04-10 | Method and device for detecting harmonic signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101858939B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015074210A1 (en) * | 2013-11-21 | 2015-05-28 | 展讯通信(上海)有限公司 | Initial search method and apparatus for user equipment and cell thereof |
CN107843767A (en) * | 2017-10-19 | 2018-03-27 | 广州致远电子有限公司 | Signal processing method and system for spectrum analysis |
CN109163672A (en) * | 2018-08-03 | 2019-01-08 | 广东工业大学 | A kind of microscopic appearance measurement method based on white light interference zero optical path difference position picking algorithm |
CN110659583A (en) * | 2019-04-26 | 2020-01-07 | 深圳市豪视智能科技有限公司 | Signal processing method and device and related product |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08248070A (en) * | 1995-03-08 | 1996-09-27 | Anritsu Corp | Frequency spectrum analyzer |
JP2006113002A (en) * | 2004-10-18 | 2006-04-27 | Nsk Ltd | Anomaly diagnosis system for mechanical equipment |
US7304597B1 (en) * | 2006-05-26 | 2007-12-04 | Lecroy Corporation | Adaptive interpolation for use in reducing signal spurs |
EP2082550B1 (en) * | 2006-11-13 | 2013-01-23 | Telefonaktiebolaget LM Ericsson (publ) | A method for limiting local bandwidth impairment using tone reservation |
CN100545651C (en) * | 2007-09-27 | 2009-09-30 | 上海交通大学 | Ceramic carrier embrittlement detection method in the ternary catalytic converter assembly |
CN100575894C (en) * | 2007-09-27 | 2009-12-30 | 上海交通大学 | Time-frequency domain identification of short-time characteristic acoustical signal and detection method under the complicated noise floor |
-
2009
- 2009-04-10 CN CN200910134924XA patent/CN101858939B/en not_active Expired - Fee Related
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015074210A1 (en) * | 2013-11-21 | 2015-05-28 | 展讯通信(上海)有限公司 | Initial search method and apparatus for user equipment and cell thereof |
US9853808B2 (en) | 2013-11-21 | 2017-12-26 | Spreadtrum Communications (Shanghai) Co., Ltd. | User equipment, initial cell searching method and apparatus thereof |
CN107843767A (en) * | 2017-10-19 | 2018-03-27 | 广州致远电子有限公司 | Signal processing method and system for spectrum analysis |
CN107843767B (en) * | 2017-10-19 | 2020-12-04 | 广州致远电子有限公司 | Signal processing method and system for spectrum analysis |
CN109163672A (en) * | 2018-08-03 | 2019-01-08 | 广东工业大学 | A kind of microscopic appearance measurement method based on white light interference zero optical path difference position picking algorithm |
CN109163672B (en) * | 2018-08-03 | 2020-06-19 | 广东工业大学 | Micro-topography measuring method based on white light interference zero-optical-path-difference position pickup algorithm |
CN110659583A (en) * | 2019-04-26 | 2020-01-07 | 深圳市豪视智能科技有限公司 | Signal processing method and device and related product |
CN110659583B (en) * | 2019-04-26 | 2023-11-28 | 深圳市豪视智能科技有限公司 | Signal processing method and device and related products |
Also Published As
Publication number | Publication date |
---|---|
CN101858939B (en) | 2012-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100476949C (en) | Multichannel voice detection in adverse environments | |
US7912709B2 (en) | Method and apparatus for estimating harmonic information, spectral envelope information, and degree of voicing of speech signal | |
CN105529028A (en) | Voice analytical method and apparatus | |
CN102915728B (en) | Sound segmentation device and method and speaker recognition system | |
CN102789779A (en) | Speech recognition system and recognition method thereof | |
CN108665903A (en) | A kind of automatic testing method and its system of audio signal similarity degree | |
Wang et al. | Digital audio tampering detection based on ENF consistency | |
CN101858939B (en) | Method and device for detecting harmonic signal | |
KR101762723B1 (en) | Method and apparatus for detecting correctness of pitch period | |
KR100827153B1 (en) | Method and apparatus for extracting degree of voicing in audio signal | |
CN111508457A (en) | Music beat detection method and system | |
CN112908344A (en) | Intelligent recognition method, device, equipment and medium for bird song | |
CN113674763A (en) | Whistling sound identification method, system, equipment and storage medium by utilizing line spectrum characteristics | |
CN115758082A (en) | Fault diagnosis method for rail transit transformer | |
US7966179B2 (en) | Method and apparatus for detecting voice region | |
CN112509601B (en) | Note starting point detection method and system | |
Deng et al. | Selection of optimal wavelet basis for signal denoising | |
CN116895288A (en) | Digital audio self-adaptive copy and paste detection method and device based on pseudo Wigner-Ville distribution | |
CN102214219A (en) | Audio/video content retrieval system and method | |
US7012186B2 (en) | 2-phase pitch detection method and apparatus | |
JP3815323B2 (en) | Frequency conversion block length adaptive conversion apparatus and program | |
WO2023093029A1 (en) | Wake-up word energy calculation method and system, and voice wake-up system and storage medium | |
CN116667920A (en) | Transient signal detection method combining principal component analysis and Hilbert-Huang transform | |
CN100424692C (en) | Audio fast search method | |
Zhao et al. | Pipeline leak fault feature extraction based on wavelet packet analysis and application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121107 Termination date: 20190410 |
|
CF01 | Termination of patent right due to non-payment of annual fee |