CN101556795A - Method and device for computing voice fundamental frequency - Google Patents

Abstract

The invention belongs to the field of signal processing and discloses a method and a device for computing voice fundamental frequency. The method and the device can more accurately estimate the fundamental frequency. In the invention, before the fundamental frequency computation, frequency domain signals in use are reconstructed to generate a reconstruction function with a continuous definitional domain. The function undergoes curve fitting in a definitional domain near each frequency domain peak value according to the corresponding frequency domain signals and effectively inhibits the corresponding frequency domain signals in other definitional domains. When the fundamental frequency is searched, candidate fundamental frequencies and a plurality of frequency multiplications thereof are comprehensively taken into consideration.

Description

The method and apparatus of computing voice fundamental frequency

Technical field

The present invention relates to field of voice signal, particularly the technology of computing voice fundamental frequency.

Background technology

Fast development along with network and multimedia technology, speech processing system has spread in each field such as broadcasting, TV, communication, from broadcasting, TV programme making apparatus all be unable to do without speech processing system to hand-held phone, portable audio/video playback apparatus.

Handle and the encoding and decoding speech field at voice signal, estimate that correctly fundamental frequency is extremely important.

From the angle of voice generation principle, voice derive from the vibration of the operatic tunes, produce sound wave, and the modulation through the sound channel organ obtains voice signal again.The vibration of the operatic tunes can determine the type of the voice signal that produces usually, vowel for example, consonant, fricative or the like.In the voice that reality occurs, vowel is seized of significant proportion.An English word is comprising the vowel of major part usually.On the angle of signal analysis, vowel is mainly by humorous wave component, and promptly its frequency component is made up of a fundamental frequency (also can abbreviate fundamental frequency as) and several its integer multiple frequencies.

Be in 4,161,625 the United States Patent (USP), to disclose a kind of method that from voice signal, obtains fundamental frequency in the patent No..In this patent, by initial voice signal is handled, obtain difference signal, adopt auto-correlation algorithm to obtain fundamental frequency again according to difference signal.

Because in the algorithm of actual speech encoding and decoding and voice signal processing (for example G.729 the encoding and decoding speech standard waits), traditional calculating fundamental frequency algorithm mainly is an auto-correlation algorithm, and promptly the maximum coefficient of autocorrelation by the computing voice signal finds specific value.Because comprise a large amount of noises in the voice signal usually, so use computing method may have certain deviation based on the fundamental frequency of auto-correlation algorithm.

Summary of the invention

The object of the present invention is to provide a kind of method and apparatus of computing voice fundamental frequency, can estimate fundamental frequency more exactly.

The invention discloses a kind of method of computing voice fundamental frequency, may further comprise the steps:

The voice signal of time domain is transformed to discrete frequency domain signal X _i, i=1 wherein, 2 ..., N;

| X _i| in find out each peak value M as local maximum _j, j=1 wherein, 2 ..., L, L are the number of peak value, || expression takes absolute value;

In the related field of definition of described discrete frequency-region signal, L nonoverlapping regional Z of structure _j, each Z _jSize be scheduled to each Z _jCover a D _j, D wherein _jBe M _jIn the pairing value of field of definition;

With each Z _jFor field of definition is constructed continuous function S respectively _j(ω), ω ∈ Z _j, satisfy | S _j(ω _i)-| X _i||＜C1, wherein ω _iBe X _iIn the pairing value of field of definition, C1 is a positive constant;

At each Z _jIn the field of definition that does not have to cover, constructed fuction S ₀(ω), $\mathrm{\ω}\∈\begin{array}{}\end{array}\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]$ And $\mathrm{\ω}\∉Z_j,$ F wherein _sBe sampling rate, satisfy S ₀(ω)＜| X _i|;

Will be by each S _j(ω) and S ₀The S that (ω) is combined into (ω) calculates fundamental frequency as frequency spectrum.

The invention also discloses a kind of equipment of computing voice fundamental frequency, comprising:

Converter unit is used for the voice signal of time domain is transformed to discrete frequency domain signal X _i, i=1 wherein, 2 ..., N;

The peak value computing unit is used for | X _i| in find out each peak value M as local maximum _j, j=1 wherein, 2 ..., L, L are the number of peak value, || expression takes absolute value;

Reconfiguration unit is used in the related field of definition of described discrete frequency-region signal, L nonoverlapping regional Z of structure _j, each Z _jSize be scheduled to each Z _jCover a D _j, D wherein _jBe M _jIn the pairing value of field of definition; With each Z _jFor field of definition is constructed continuous function S respectively _j(ω), ω ∈ Z _j, satisfy | S _j(ω _i)-| X _i||＜C1, wherein ω _iBe X _iIn the pairing value of field of definition, C1 is a positive constant; At each Z _jIn the field of definition that does not have to cover, constructed fuction S ₀(ω), $\mathrm{\ω}\∈\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]$ And $\text{\ω\∉}{Z}_j,$ F wherein _sBe sampling rate, satisfy S ₀(ω)＜| X _i|;

The fundamental tone computing unit is used for by each S _j(ω) and S ₀The S that (ω) is combined into (ω) calculates fundamental frequency as frequency spectrum.

Embodiment of the present invention compared with prior art, the key distinction and effect thereof are:

Before calculating fundamental frequency, earlier used frequency-region signal is reconstructed, generate the reconstruction of function that field of definition is continuous, carry out curve fitting by corresponding frequency-region signal near the field of definition of this function each frequency domain peak value, in other field of definition, corresponding frequency-region signal is effectively suppressed.Because candidate's fundamental frequency and multiple frequency thereof are usually expressed as peak value, so, can improve accuracy and antijamming capability that fundamental frequency calculates by keeping near the frequency-region signal in the field of definition each peak value, significantly weakening frequency-region signal in other field of definition.The frequency-region signal that obtains by conversion disperses, by can more accurately advancing in the frequency spectrum of reconstruction of function representative to search for fundamental frequency to the serialization of field of definition.

Further, when pitch search, take all factors into consideration candidate's fundamental frequency and a plurality of frequency multiplication thereof, can make Search Results more accurate.

Further, can be by the reconstruction of function value in other field of definition outside near the field of definition peak value be made as 0, thus weaken irrelevant frequency component to greatest extent, further improve accuracy and antijamming capability that fundamental frequency calculates.

Description of drawings

Fig. 1 is the method flow diagram according to a kind of computing voice fundamental frequency of first embodiment of the invention;

Fig. 2 is the equipment structure chart according to a kind of computing voice fundamental frequency of third embodiment of the invention.

Embodiment

For making the purpose, technical solutions and advantages of the present invention clearer, embodiments of the present invention are described in further detail below in conjunction with accompanying drawing.

First embodiment of the present invention relates to a kind of method of computing voice fundamental frequency, as shown in Figure 1.

In step 110, the voice signal of importing is transformed into frequency domain from time domain.Specifically, the time domain voice signal of supposing input is x _i, i=1,2 ..., N then can be converted into discrete frequency domain signal X by Fast Fourier Transform (FFT) (FastFourier Transform is called for short " FFT ") _i, i=1 wherein, 2 ..., N.

Need to prove that in this step, the conversion from the time-domain signal to the frequency-region signal is embodied as example with FFT and describes, but in actual applications, also can realize by other modes.Such as, can be by modes such as discrete cosine transform (Discrete Cosine Transform, be called for short " DCT ") or modified discrete cosine transforms, the voice signal of time domain is transformed to discrete frequency-region signal.

Then, in step 120, from each frequency domain signal X _iAbsolute value in, find out each peak value M as local maximum _j, j=1 wherein, 2 ..., L, L are the number of peak value.Such as, earlier according to X ₁, X ₂...., X _N, obtain the absolute value Y of each frequency-region signal ₁, Y ₂...., Y _N, Y wherein _i=| X _i|, i=1,2 ..., N.Then, search for local maximum Y again _i, as search for all and satisfy Y _i＞max (Y _I+1, Y _I-1) Y _iAs peak value M _j, the choosing method of this local maximum is actually among 3 o'clock and chooses maximal value, certainly, in actual applications, also can (among as 5 or 6) choose maximal value in more point.

Then, in step 130, reconstruct the continuous frequency spectrum of voice signal according to the frequency domain peak value of electing.Because each frequency-region signal that obtains after the FFT conversion disperses, and the field of definition continuous functions of frequency spectrum promptly can't be provided, this calculating to fundamental tone has caused difficulty, therefore needs reconstruct and serialization frequency spectrum.Specific as follows:

At first, entire spectrum is divided into two types.One type frequency spectrum is corresponding with the frequency component of fundamental frequency or its certain multiple, and the frequency spectrum of another kind of type then is and the irrelevant pairing frequency spectrum of frequency component of fundamental frequency.Because the frequency component of fundamental frequency and its certain multiple all shows as a local maximum usually on frequency spectrum, can think that therefore in step 120 selected peak value represented the frequency component of fundamental frequency and its certain multiple.Then be considered to the frequency component that has nothing to do in other parts of whole frequency axis.

Secondly, respectively this frequency spectrum of two types is carried out function reconstruct.Specifically, in the related field of definition of discrete frequency-region signal, L of structure respectively with L the corresponding nonoverlapping regional Z of peak value _j, each Z _jSize be scheduled to each Z _jCover a D _j, D wherein _jBe M _jIn the pairing value of field of definition.With each Z _jFor field of definition is constructed continuous function S respectively _j(ω), ω ∈ Z _j, satisfy | S _j(ω _i)-| X _i||＜C1, wherein ω _iBe X _iIn the pairing value of field of definition, C1 is a positive constant.At each Z _jIn the field of definition that does not have to cover, constructed fuction S ₀(ω), $\mathrm{\ω}\∈\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]$ And $\mathrm{\ω}\∉Z_j,$ F wherein _sBe sampling rate, satisfy S ₀(ω _i)＜| X _i|.Below to continuous function S _jMake (ω) is described further.

In the present embodiment, by with the frequency-region signal absolute value of peak value correspondence with and former and later two frequency-region signal absolute values carry out the binomial interpolation, realize continuous function S _jStructure (ω).Such as, first peak value M ₁Corresponding frequency-region signal absolute value Yi is ω in the pairing value of field of definition _i, its former and later two frequency-region signal absolute value (Y then _I-1, Y _I+1) be (ω in the pairing value of field of definition _I-1, ω _I+1).Suppose that this interpolation polynomial represented by following second degree trinomial expression:

f(x)＝ax ²+bx+c

Then by Substitution method can solve corresponding coefficient a, b, c}:

$[a,b,c]=[Y_{i-1},Y_i,Y_{i+1}]{\left(\begin{array}{ccc}{\mathrm{\ω}}_{i-1}^2& {\mathrm{\ω}}_i^2& {\mathrm{\ω}}_{i+1}^2\\ {\mathrm{\ω}}_{i-1}& {\mathrm{\ω}}_i& {\mathrm{\ω}}_{i+1}\\ 1& 1& 1\end{array}\right)}^{-1}$

Therefore, can obtain: S ₁(ω)=a ₁ω ²+ b ₁ω+c ₁

In like manner, can construct each S by the binomial interpolation method _j(ω) function, i.e. S _j(ω)=a _jω ²+ b _jω+c _j

Need to prove, because in the present embodiment, peak value M _jBe actually the maximal value of among 3 o'clock, choosing, if so first peak value M ₁Corresponding frequency-region signal absolute value is Y _i, then regional Z ₁Reference position be Y _I-1In the pairing value of field of definition, end position is Y _I+1In the pairing value of field of definition, promptly $Z_1=[\frac{F_s}{N}{\mathrm{\ω}}_{i-1},\frac{F_s}{N}{\mathrm{\ω}}_{i+1}],$ In like manner, can obtain each Z _jThe zone.The field of definition of the matched curve of each peak value correspondence also can be taked other arbitrarily rational length

For each Z _jThe field of definition that does not have covering, because Pitch Information is not contained in these zones, therefore can be simply with these parts arbitrary function S ₀(ω) replace, $\mathrm{\ω}\∈\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]$ And $\mathrm{\ω}\∉Z_j,$ Function S ₀(ω) only need satisfy S ₀(ω _i)＜| X _i| this condition gets final product.Such as, adopt null function, that is: S ₀(ω)=0.

Because in this step, respectively this frequency spectrum of two types has been carried out function reconstruct, so entire spectrum all to be reconstructed into be a field of definition continuous functions, that is:

$S\left(\mathrm{\ω}\right)=\begin{array}{}\end{array}\left\{\begin{array}{c}{S}_1\left(\mathrm{\ω}\right)=a_1{\mathrm{\ω}}^2+b_1\mathrm{\ω}+c_1,\mathrm{\ω}\∈Z_1\\ S_2\left(\mathrm{\ω}\right)=a_2{\mathrm{\ω}}^2+b_2\mathrm{\ω}+c_2,\mathrm{\ω}\∈Z_2\\ \·\\ \·\\ \·\\ S_L\left(\mathrm{\ω}\right)=a_L{\mathrm{\ω}}^2+b_L\mathrm{\ω}+c_L,\mathrm{\ω}\∈Z_L\\ S_0\left(\mathrm{\ω}\right)=0,\mathrm{\ω}\∈[0,\frac{F_s}{2}]\mathrm{and\ω}{\∉Z}_j\end{array}\right.$

Then, enter step 140, calculate fundamental frequency.Specifically, because in step 130, drawn a field of definition continuous functions S (ω), can directly derive fundamental frequency according to the function characteristics of this function S (ω).Such as, to search in the scope (as from 50 hertz to 500 hertz) that may exist of fundamental tone, the criterion of search is to find the frequency that satisfies following formula:

${\mathrm{\ω}}_p=\underset{\mathrm{\ω}}{\arg \max }\underset{k=1}{\overset{N\left(\mathrm{\ω}\right)}{\mathrm{\Σ}}}{\left|S\left(\mathrm{k\ω}\right)\right|}^2$

Wherein, N (ω) is to be the harmonic wave number of fundamental frequency with ω, ω _pBe fundamental frequency.Need to prove that above-mentioned formula is an object lesson as search criteria, in actual applications, also can adopt other formula, as will square changing to 4 powers or 1 power etc. in the above-mentioned formula.Above-mentioned ω _pThe essence of correlation formula is to take all factors into consideration candidate's fundamental frequency and a plurality of frequency multiplication thereof when pitch search, and concrete formula form can have other variation, can make Search Results more accurate like this

Because in the present embodiment, before calculating fundamental frequency, earlier used frequency-region signal is reconstructed, generate the reconstruction of function that field of definition is continuous, carry out curve fitting by corresponding frequency-region signal near the field of definition of this function each frequency domain peak value, in other field of definition, corresponding frequency-region signal is effectively suppressed.Because candidate's fundamental frequency and multiple frequency thereof are usually expressed as peak value, so, can improve accuracy and antijamming capability that fundamental frequency calculates by keeping near the frequency-region signal in the field of definition each peak value, significantly weakening frequency-region signal in other field of definition.The frequency-region signal that obtains by conversion disperses, by can more accurately advancing in the frequency spectrum of reconstruction of function representative to search for fundamental frequency to the serialization of field of definition.

What deserves to be mentioned is, in the present embodiment, to each Z _jThe function S that the field of definition that does not have to cover is constructed ₀(ω) be: S ₀(ω)=0, thereby weaken irrelevant frequency component to greatest extent, further improve accuracy and antijamming capability that fundamental frequency calculates.And in actual applications, also can be with function S ₀(ω) be changed to a very little value, can search fundamental frequency comparatively exactly equally.

Second embodiment of the present invention relates to a kind of method of computing voice fundamental frequency, and the present embodiment and first embodiment are roughly the same, and its difference is, in the first embodiment, is constructing continuous function S _jIn the time of (ω), be by with the frequency-region signal absolute value of peak value correspondence with and former and later two frequency-region signal absolute values carry out the binomial interpolation and realize; And in the present embodiment, can be by fitting to the segmentation straight line, or come match with cubic polynomial, realize continuous function S _jStructure (ω).

Method embodiment of the present invention can be realized in software, hardware, firmware or the like mode.No matter the present invention be with software, hardware, or the firmware mode realize, instruction code can be stored in the storer of computer-accessible of any kind (for example permanent or revisable, volatibility or non-volatile, solid-state or non-solid-state, medium fixing or that change or the like).Equally, storer can for example be programmable logic array (Programmable Array Logic, be called for short " PAL "), random access memory (Random Access Memory, be called for short " RAM "), programmable read only memory (Programmable Read Only Memory, be called for short " PROM "), ROM (read-only memory) (Read-Only Memory, be called for short " ROM "), Electrically Erasable Read Only Memory (Electrically Erasable Programmable ROM, be called for short " EEPROM "), disk, CD, digital versatile disc (Digital Versatile Disc is called for short " DVD ") or the like.

The 3rd embodiment of the present invention relates to a kind of equipment of computing voice fundamental frequency, as shown in Figure 2, comprising: converter unit is used for the voice signal of time domain is transformed to discrete frequency domain signal X _i, i=1 wherein, 2 ..., N; The peak value computing unit is used for | X _i| in find out each peak value M as local maximum _j, j=1 wherein, 2 ..., L, L are the number of peak value, || expression takes absolute value; Reconfiguration unit is used in the related field of definition of discrete frequency-region signal, L nonoverlapping regional Z of structure _j, each Z _jSize be scheduled to each Z _jCover a D _j, D wherein _jBe M _jIn the pairing value of field of definition; With each Z _jFor field of definition is constructed continuous function S respectively _j(ω), ω ∈ Z _j, satisfy | S _j(ω _i)-| X _i||＜C1, wherein ω _iBe X _iIn the pairing value of field of definition, C1 is a positive constant; At each Z _jIn the field of definition that does not have to cover, constructed fuction S ₀(ω), $\mathrm{\ω}\∈\begin{array}{}\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]\end{array}$ And $\mathrm{\ω}\∉Z_j,$ F wherein _sBe sampling rate, satisfy S ₀(ω i)＜| X _i|; The fundamental tone computing unit is used for by each S _j(ω) and S ₀The S that (ω) is combined into (ω) calculates fundamental frequency as frequency spectrum.

The fundamental tone computing unit calculates fundamental frequency in the following manner: fundamental tone may have scope search, the criterion of search is to find the frequency that satisfies following formula:

${\mathrm{\ω}}_p=\underset{\mathrm{\ω}}{\arg \max }\underset{k=1}{\overset{N\left(\mathrm{\ω}\right)}{\mathrm{\Σ}}}{\left|S\left(\mathrm{k\ω}\right)\right|}^2$

Wherein, N (ω) is to be the harmonic wave number of fundamental frequency with ω, ω _pBe fundamental frequency.

Converter unit can adopt modes such as FFT, discrete cosine transform, modified discrete cosine transform, and the voice signal of time domain is transformed to discrete frequency-region signal.

Reconfiguration unit one of can be in the following ways realized S _jStructure (ω): with the frequency-region signal absolute value of peak value correspondence with and former and later two frequency-region signal absolute values carry out the binomial interpolation or fit to the segmentation straight line, or come match with cubic polynomial.

Need to prove, each unit of mentioning in the present embodiment all is a logical block, physically, a logical block can be a physical location, it also can be the part of a physical location, can also realize that the physics realization mode of these logical blocks itself is not most important with the combination of a plurality of physical locations, the combination of the function that these logical blocks realized is the key that just solves technical matters proposed by the invention.

In addition, for outstanding innovation part of the present invention, present embodiment will not introduced not too close unit with solving technical matters relation proposed by the invention, and this does not show that there is not other unit in this equipment embodiment.

Though by with reference to some preferred embodiment of the present invention, the present invention is illustrated and describes, those of ordinary skill in the art should be understood that and can do various changes to it in the form and details, and without departing from the spirit and scope of the present invention.

Claims (10)

1. the method for a computing voice fundamental frequency is characterized in that, may further comprise the steps:

The voice signal of time domain is transformed to discrete frequency domain signal X _i, i=1 wherein, 2 ..., N;

| X _i| in find out each peak value M as local maximum _j, j=1 wherein, 2 ..., L, L are the number of peak value, || expression takes absolute value;

In the related field of definition of described discrete frequency-region signal, L nonoverlapping regional Z of structure _j, each Z _jSize be scheduled to each Z _jCover a D _j, D wherein _jBe M _jIn the pairing value of field of definition;

With each Z _jFor field of definition is constructed continuous function S respectively _j(ω), ω ∈ Z _j, satisfy | S _j(ω _i)-| X _i||＜C1, wherein ω _iBe X _iIn the pairing value of field of definition, C1 is a positive constant;

At each Z _jIn the field of definition that does not have to cover, constructed fuction S ₀(ω), $\mathrm{\ω}\∈\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]$ And $\mathrm{\ω}\∉Z_j,$ F wherein _sBe sampling rate, satisfy S ₀(ω _i)＜| X _i|;

Will be by each S _j(ω) and S ₀The S that (ω) is combined into (ω) calculates fundamental frequency as frequency spectrum.

2. the method for computing voice fundamental frequency according to claim 1 is characterized in that, and is described with the step of S (ω) as frequency spectrum calculating fundamental frequency, realizes by following substep:

Fundamental tone may have scope search, the criterion of search is to find the frequency that satisfies following formula:

${\mathrm{\ω}}_p=\underset{\mathrm{\ω}}{\arg \max }\underset{k=1}{\overset{N\left(\mathrm{\ω}\right)}{\mathrm{\Σ}}}{\left|S\left(\mathrm{k\ω}\right)\right|}^2$

Wherein, N (ω) is to be the harmonic wave number of fundamental frequency with ω, ω _pBe the result of calculation of fundamental frequency.

3. the method for computing voice fundamental frequency according to claim 2 is characterized in that, is transformed in the step of discrete frequency-region signal at the voice signal with time domain, adopts one of following mapping mode:

Fast fourier transform, discrete cosine transform, modified discrete cosine transform.

4. the method for computing voice fundamental frequency according to claim 3 is characterized in that,

Described at | X _i| in find out each peak value M as local maximum _jStep comprise following substep:

Calculate Y _i=| X _i|;

Search for all and satisfy Y _i＞max (Y _I+1, Y _I-1) Y _iAs peak value M _j

5. the method for computing voice fundamental frequency according to claim 4 is characterized in that,

At described structure continuous function S _jIn the step (ω), realize one of in the following ways S _jStructure (ω):

With the frequency-region signal absolute value of peak value correspondence with and former and later two frequency-region signal absolute values carry out the binomial interpolation or fit to the segmentation straight line, or come match with cubic polynomial.

6. the method for computing voice fundamental frequency according to claim 5 is characterized in that,

Described S ₀(ω)=0.

7. the equipment of a computing voice fundamental frequency is characterized in that, comprising:

Converter unit is used for the voice signal of time domain is transformed to discrete frequency domain signal X _i, i=1 wherein, 2 ..., N;

The peak value computing unit is used for | X _i| in find out each peak value M as local maximum _j, j=1 wherein, 2 ..., L, L are the number of peak value, || expression takes absolute value;

Reconfiguration unit is used in the related field of definition of described discrete frequency-region signal, L nonoverlapping regional Z of structure _j, each Z _jSize be scheduled to each Z _jCover a D _j, D wherein _jBe M _jIn the pairing value of field of definition; With each Z _jFor field of definition is constructed continuous function S respectively _j(ω), ω ∈ Z _j, satisfy | S _j(ω _i)-| X _i||＜C1, wherein ω _iBe X _iIn the pairing value of field of definition, C1 is a positive constant; At each Z _jIn the field of definition that does not have to cover, constructed fuction S ₀(ω), $\mathrm{\ω}\∈\left[\begin{array}{cc}0& \frac{F_s}{2}\end{array}\right]$ And $\mathrm{\ω}\∉Z_j,$ F wherein _sBe sampling rate, satisfy S ₀(ω _i)＜| X _i|;

The fundamental tone computing unit is used for by each S _j(ω) and S ₀The S that (ω) is combined into (ω) calculates fundamental frequency as frequency spectrum.

8. the equipment of computing voice fundamental frequency according to claim 7 is characterized in that,

Described fundamental tone computing unit calculates fundamental frequency in the following manner:

Fundamental tone may have scope search, the criterion of search is to find the frequency that satisfies following formula:

${\mathrm{\ω}}_p=\underset{\mathrm{\ω}}{\arg \max }\underset{k=1}{\overset{N\left(\mathrm{\ω}\right)}{\mathrm{\Σ}}}{\left|S\left(\mathrm{k\ω}\right)\right|}^2$

Wherein, N (ω) is to be the harmonic wave number of fundamental frequency with ω, ω _pBe the result of calculation of fundamental frequency.

9. the equipment of computing voice fundamental frequency according to claim 8 is characterized in that,

One of the following mapping mode that adopts described converter unit realizes the voice signal of time domain is transformed to discrete frequency-region signal:

Fast fourier transform, discrete cosine transform, modified discrete cosine transform.

10. the equipment of computing voice fundamental frequency according to claim 9 is characterized in that,

Described reconfiguration unit is realized one of in the following ways S _jStructure (ω):

With the frequency-region signal absolute value of peak value correspondence with and former and later two frequency-region signal absolute values carry out the binomial interpolation or fit to the segmentation straight line, or come match with cubic polynomial.

Images