Method and system for speech recognition using continuous density hidden Markov models

Technical Field

0001The present invention relates to computer speech recognition, and more particularly, to a computer speech recognition system that utilizes continuous hidden Markov models.

Background of the Invention

0002The area of speech recognition is challenged by the need to produce a speaker-independent continuous speech recognition system which has a minimal recognition error rate. The focus in realizing this goal is on the recognition algorithm that is utilized by the speech recognition system. The recognition algorithm is essentially a mapping of the speech signal, a continuous-time signal, to a set of reference patterns representing the phonetic and phonological descriptions of speech previously obtained from training data. In order to perform this mapping, signal processing techniques such as fast fourier transforms (FFT), linear predictive coding (LPC), or filter banks are applied to a digital form of the speech signal to extract an appropriate parametric representation of the speech signal. A commonly-used representation is a feature vector containing for each time interval, the FFT or LPC coefficients that represent the frequency and/or energy bands contained in the speech signal. A sequence of these feature vectors is mapped to the set of reference patterns which identify linguistic units, words and/or sentences contained in the speech signal.

0003Often, the speech signal does not exactly match the stored reference patterns. The difficulty in finding an exact match is due to the great degree of variability in speech signal characteristics which are not completely and accurately captured by the stored reference patterns. Probabilistic models and statistical techniques have been used with more success in predicting the intended message than techniques that seek an exact match. One such technique is Hidden Markov Models (HMMs). These techniques are more adept for speech recognition since they determine the reference pattern that will more likely match the speech signal rather than finding an exact match.

0004A HMM consists of a sequence of states connected by transitions. A HMM can represent a particular phonetic unit of speech, such as a phoneme or word. Associated with each state is an output probability indicating the likelihood that the state matches a feature vector. For each transition, there is an associated transition probability indicating the likelihood of following the transition. The transition and output probabilities are estimated statistically from previously spoken speech patterns, referred to as "training data." The recognition problem is one of finding the state sequence having the highest probability of matching the feature vectors representing the input speech signal. Primarily, this search process involves enumerating every possible slate sequence that has been modeled and determining the probability that the state sequence matches the input speech signal. The utterance corresponding to the state sequence with the highest probability is selected as the recognized speech utterance.

0005Most HMM-based speech recognition systems are based on discrete HMMs utilizing vector quantization. A discrete HMM has a finite set of output symbols and the transition and output probabilities are based on discrete probability distribution functions (pdfs). Vector quantization is used to characterize the continuous speech signal by a discrete representation referred to as a codeword. A feature vector is matched to a codeword using a distortion measure The feature vector is replaced by the index of the codeword having the smallest distortion measure. The recognition problem is reduced to computing the discrete output probability of an observed speech signal as a table look-up operation which requires minimal computation.

0006However, speech signals are continuous signals. Although it is possible to quantitize continuous signals through codewords, there may be serious degradation associated with such quantization resulting in poor recognition accuracy. Recognition systems utilizing continuous density HMMs do not suffer from the inaccuracy associated with quantization distortion. Continuous density HMMs are able to directly model the continuous speech signal using estimated continuous density probability distribution functions, thereby achieving a higher recognition accuracy. However, continuous density HMMs require a considerable amount of training data and require a longer recognition computation which has deterred their use in most commercial speech recognition systems. Accordingly, a significant problem in continuous speech recognition systems has been the use of continuous density HMMs for achieving high recognition accuracy.

0007Komori Y Et Al: "An efficient output probability computation for continuous HMM using rough and detail models", European Conference on Speech Communication and Technology (Eurospeech), vol. 2, 18 September, 1995, pages 1087-1090, XP000670819 describes a computation algorithm of the output probability for a continuous HMM speech recognizer using a rough and detail HMM combination. This algorithm first estimates the state output probabilities using the rough HMMs and re-estimates those of the probable states using detail HMMs in order to reduce the amount of the output probability computation. Two basic types for the rough HMM and the detail HMM combination are described, firstly a combination of small and large number of mixture HMMs and secondly a combination of context-independent and dependent HMMs. Both rough HMMs and detail HMMs are trained independently by an ordinary expectation-maximization algorithm using the speech database, the labels and the acoustic parameters. The HMM is a continuous HMM of a diagonal co-variance with multi-mixture Gaussian density functions. Two models are adopted, a context independent phone HMM (25 models) and a right context dependent HMM (243 models).

0008EP-A-0 706 171 describes a speech recognition method and apparatus that uses continuous mixture HMMs for probability processing. For the same phoneme of HMMs, two types of HMMs, a phonemic context-independent HMM and a phonemic context-dependent HMM in which a phoneme is dependent on a plurality of surrounding phonemes are employed. The output probabilities of phonemes and their states and mixtures which are likely to contribute to the recognition results are first determined by use of the phonemic context-independent HMM. The output probabilities of the phonemes with the highest probability values and their states and mixtures determined by use of the phonemic context-independent HMM are re-calculated by use of the phonemic context-dependent HMM.

0009It is, therefore, the object of the present invention to provide an improved speech recognition system which improves the modeling of the speech signal to continuous density HMM corresponding to a linguistic expression.

0010This object is solved by the subject-matter of the independent claims. Preferred embodiments are the subject-matter of the dependent claims.

0011In the preferred embodiment, the recognition system utilizes a context-independent and several context-dependent HMMs to represent the speech unit of a phoneme in different contextual patterns. The output and transition probabilities for each of these HMMs are estimated from the training data. The output probabilities associated with like states corresponding to the same modeled phoneme are clustered forming senones. A weighting factor for each context-dependent senone which indicates the robustness of the output probability in predicting unseen data is also generated. In the preferred embodiment, the weighting factor is estimated through deleted interpolation of all the data points in the training data. Alternatively, the weighting factor can be estimated from a parametric representation of the data points or from randomly-generated data points generated from a parametric representation of the data points.

0012The recognition engine receives an input speech utterance and generates candidate word sequences which will most likely match the feature vectors of the input speech utterance. The word sequences can be composed of various senone alignments corresponding to state sequences of HMMs. The recognition engine determines which senone/state alignment best matches the feature vectors by utilizing an acoustic and language probability score. The acoustic probability score represents the likelihood that the senone alignment corresponds to the feature vectors and the language probability score indicates the likelihood of the utterance corresponding to the senone alignment occurring in the language. The acoustic probability score is based on an output and transition probability analysis. The output probability analysis utilizes the output probabilities of both the context-dependent and context-independent senones by weighing each output probability as a function of the weighting factor. The output probability having the more robust estimate will dominate the analysis thereby improving the output probability analysis. An improvement in the output probability analysis improves the acoustic score and, in turn, the overall recognition accuracy.

Brief Description of the Drawings

0013The foregoing and other features and advantages of the invention will be apparent from the following more particular description of the preferred embodiment of the invention, as illustrated in the accompanying drawings in which like reference characters refer to the same elements throughout the different views The drawings are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the invention.

0014Figure 1 is a block diagram of a speech recognition system employed in the preferred embodiment.

0015Figure 2 is a flow diagram of a training method used in the system of Figure 1.

0016Figure 3 is flow diagram of the method for calculating weighing factors used in the system of Figure 1.

0017Figure 4 is a flow diagram of the preferred embodiment for calculating a new value for lambda as used in the system of Figure 3.

0018Figure 5 is a flow diagram of a first alternate embodiment for calculating a new value for lambda as used in the system of Figure 3.

0019Figure 6 is a flow diagram of a second alternate embodiment for calculating a new value for lambda as used in the system of Figure 3.

0020Figures 7A and 7B depict an example of the hidden Markov models and senone structures associated with a phoneme.

0021Figure 8 is a flow diagram of the speech recognition method used in the system of Figure 1.

Detailed Description of the Invention

0022The preferred embodiment of the present invention recognizes that an improved recognition accuracy can be obtained in speech recognition systems that employ continuous density hidden Markov models by weighing different output probabilities representing the same phonetic unit relative to the degree to which each output probability can predict unseen data. The speech recognition system of the claimed invention receives an input speech utterance, in the form of a continuous signal, and generates the most likely linguistic expression that corresponds to the utterance. The preferred embodiment recognizes a linguistic expression by matching a set of feature vectors that form a parametric representation of the speech signal to a sequence of hidden Markov models (HMMs) which identify possible linguistic expressions. A HMM may represent a phoneme, and a sequence of HMMs may represent words or sentences composed of phonemes.

0023Continuous density probability distribution functions, such as a mixture of Gaussian probability distribution functions, can be utilized to represent the output probability of a state, since they are more accurate at modeling a speech signal. The output probability function is statistically estimated from training data. Often there is an insufficient amount of training data to accurately estimate the output probability function. To account for this problem, context-independent and context-dependent models are constructed for a predetermined set of phonemes. The output probabilities of the context-independent model are then interpolated with the output probabilities of the context-dependent model. This is done through a weighting or interpolation factor which estimates the degree to which the output probability function of the context-dependent HMM can predict data not previously encountered in the training data. Thus, the new modified output probability function of a context-dependent state is a combination of the output probability functions of both models weighed in accordance with the robustness of the estimates. Accordingly, in the preferred embodiment, deleted interpolation is used to smooth the probability space rather than the parameter space.

0024Figure 1 illustrates a speech recognition system 10 which can be used to implement the recognition and training processes in accordance with the preferred embodiment of the invention. The speech recognition system 10 contains an input device 12, such as but not limited to a microphone, which receives an input speech utterance and generates a corresponding analog electric signal. Alternatively, a stored speech utterance that is stored on a storage device can be used as the input speech utterance. The analog electric signal corresponding to the speech utterance is transmitted to analog-to-digital (A/D) converter 14, which converts the analog signal to a sequence of digital samples. The digital samples are then transmitted to a feature extractor 16 which extracts a parametric representation of the digitized input speech signal. This parametric representation captures the acoustic properties of the input speech utterance. Preferably, the feature extractor 16 performs spectral analysis to generate a sequence of feature vectors, each of which contains coefficients representing a spectra of the input speech signal. Methods for performing the spectral analysis are well-known in the art of signal processing and can include fast Fourier transforms (FFT), linear predictive coding (LPC), and cepstral coefficients, all of which can be utilized by feature extractor 16. Feature extractor 16 can be any conventional processor that performs spectral analysis. Spectral analysis may be performed every ten milliseconds to divide the input speech signal into a feature vector which represents twenty-five milliseconds of the utterance. However, this invention is not limited to using feature vectors that represent twenty-five milliseconds of the utterance. Feature vectors representing different lengths of time of a speech utterance can also be used. This process is repeated for the entire input speech signal and results in a sequence of feature vectors which is transmitted to a data processor 38. Data processor 38 can be any conventional computer, such as a desktop personal computer. The data processor contains a switching block 18 which routes the sequence of feature vectors. Switching block 18 can be implemented in hardware or software. However, the speech recognition system is not limited to executing on a data processor. Other types of executable mediums can be used, such as but not limited to, a computer readable storage medium which can be a memory device, compact disc, or floppy disk.

0025During the initial training phase, switching block 18 is switched to direct the feature vectors to the training engine 20. Training engine 20 utilizes the feature vectors to estimate the parameters of the HMMs which will represent the phonemes present in the training data and to compute a set of weighting factors for use by the recognition engine 34. A more detailed description of the method employed by training engine 20 is presented below with reference to Figures 2-6. Briefly, the training engine 20 generates context-independent and context-dependent phoneme-based hidden Markov models by estimating the parameters for these models from the training data. The output distributions for each context-dependent state are clustered forming senones which are stored in senone table storage 30. The senone table storage 30, in general, holds senones for both context-dependent and context-independent HMMs. Senones identifiers for each HMM are stored in the HMM storage 28 The In addition, a weighting factor for each context-dependent senone is calculated and stored in lambda table storage 26 for use by the recognition engine 34. The lambda table storage 26 holds lambda values indexed by context-dependent HMMs. The training engine 20 also utilizes a text transcript that holds a translation 22 of the training data and a dictionary 24 that contains a phonemic description of each word in order to assure that each word is correctly modeled. A more detailed description of the operation of the training engine 20 will be discussed below. The dictionary 24 contains a pronunciation of each word in terms of phonemes. For example, a dictionary entry for "add" might be "/AE DD/"

0026After the initial training phase, switching block 18 is switched to transmit the feature vectors to recognition engine 34. Recognition engine 34 recognizes the sequence of feature vectors as a linguistic expression composed of phonemes that form words, which, in turn, form sentences. A detailed description of the method employed by the recognition engine 34 is presented below with reference to Figure 8. Recognition engine 34 utilizes the context-independent and context-dependent hidden Markov models stored in HMM storage 28, the context-dependent and context-independent senones stored in senone table storage 30, the weighting factors stored in lambda table storage 26, and a language model stored in language model storage 32 and dictionary 24. The language model storage 22 may specify a grammar. In the preferred embodiment, the linguistic expression which is generated from recognition engine 34 is displayed to an output device 36, such as a conventional printer, computer monitor, or the like. However, this invention is not limited to displaying the linguistic expression to an output device. For example, the linguistic expression can be used as input into another program or processor for further processing or may be stored

0027Figures 2-6 are flow charts that illustrate the steps performed in the training phase of the system where the parameters of the HMMs and senones are estimated and the weighting factors are calculated. In short, the training method starts off by receiving input speech utterances, in the form of words, sentences, paragraphs, or the like, and converts them to parametric representations, known as feature vectors. The structure of the hidden Markov models and senones is formed and estimates of the parameters for these data structures are calculated from the training data. The weighting factors are then determined through the technique of deleted interpolation.

0028Referring to Figure 2, the training method commences by receiving a sequence of speech utterances (step 42) which is converted into a sequence of feature vectors (step 44) as previously described above with reference to Figure 1. The complete set of feature vectors is referred to as the "training data." In the preferred embodiment, LPC cepstral analysis is employed to model the speech signal and results in a feature vector that contains the following 39 cepstral and energy coefficients representing the frequency and energy spectra contained in the signal: (1) 12 LPC mel-frequency cepstral coefficients, <i>x</i>_<i>k</i>(<i>t</i>), for 1 &lt;= <i>K</i> &lt;= 12; (2) 12 LPC delta mel-frequency cepstral coefficients Δ<i>x</i>_<i>k</i>(<i>l</i>), for 1&lt;= <i>K</i> &lt;= 12; (3) 12 LPC delta-delta mel-frequency cepstral coefficients ΔΔ<i>x</i>_<i>k</i>(<i>t</i>), for 1 &lt;= <i>K</i> &lt;= 12; and (4) energy, delta energy, and delta-delta energy coefficients. The use of LCP cepstral analysis to model speed signals is well known in the art of speech recognition systems.

0029In step 46, senone and HMM data structures are generated. Senones are well-known data structure in speech recognition systems and a detailed description of senones and the method used to construct them can be found in M. Huang et al., "Predicting Unseen Triphones with Senones," <i>Proc. lCASSP '93</i> Vol. II, pp. 311-314, 1993. In the preferred embodiment, a HMM can be used to model the speech unit of a phoneme. The HMM can also be referred to as an acoustic model. This speech unit is selected in order to accommodate large-vocabulary recognition. Modeling individual words requires a longer training period and additional storage to store the associated parameters. This is feasible for small vocabulary systems but impractical for those which utilize large vocabularies. However, this invention is not limited to phoneme-based HMMs. Other speech units, such as words, diphones, and syllables can be used as the basis of the HMMs.

0030Two types of HMMs can be utilized. A context-dependent HMM can be used to model a phoneme with its left and right phonemic contexts. This type of model captures the contextual dependencies which are usually present in word modeling. A context-independent HMM can be used to model the phoneme in any context that it appears in the training data, therefore making it independent of any particular context. Predetermined patterns consisting of a set of phonemes and their associated left and right phonemic contexts are selected to be modeled by the context-dependent HMM. These selected patterns represent the most frequently occurring phonemes and the most frequently occurring contexts of these phonemes. The training data will provide the estimates for the parameters of these models. The context-independent models will be based on the selected phonemes and modeled within whatever phonemic context appears in the training data. Similarly, the training data will provide the estimates for the parameters of the context-independent models.

0031The use of both a context-independent and context-dependent model is beneficial in achieving an improved recognition accuracy. Each model's robustness is related to the amount of training data used to estimate its parameters which also enables it to predict data not present in the training data. The combination of the two models provides a more robust estimate benefiting from the training of both models. For example, the context-dependent model is beneficial at modeling co-articulatory effects but may be poorly trained due to limited training data. (Although a speaker may try to pronounce words as concatenated sequences of phones, the speaker's articulator cannot move instantaneously to produce unaffected phones. As a result, a phone is strongly inverted by the phone that precedes it and the phone that follows it in a word. The effects are "coarticulatory effects"). By contrast, the context-independent model is highly trainable thereby producing more robust estimates which are less detailed. The combination of these two models, weighed in the appropriate manner, can be used by the recognition engine to produce a more accurate acoustic probability score.

0032Further to account for the between-speaker differences, such as the formant frequencies (<i>i</i>.<i>e</i>., vocal tract resonance frequencies), present in the male and female vocal tracts, the HMM can utilize a mixture of unimodal distributions for the output probability distributions functions (referred to throughout this application as "output pdf"). Preferably, a mixture of Gaussian probability density functions can be used. However, this invention is not constrained to this particular limitation. Mixtures of other well-known continuous density functions can be used, such as the Laplacian and K₀- type density functions.

0033Further, to capture similarity between states of different context-dependent phonemes and to increase the amount of training data available for each senone, the output distributions of like states of different context-dependent phonetic HMM models for the same context-independent phone are clustered together forming senones.

0034Figure 7A illustrates an example of a context-independent HMM structure for the phoneme /aa/ 114. The context-independent HMM includes three states, denoted as state 1(111), state 2 (112), and state 3 (113). The HMM depicted in Figure 7A models the phoneme /aa/ with any left and right phonemes that appear in the training data as designated by the notation (*,*) in Figure 7A. The first position within the parentheses designates the phoneme that precedes the given phoneme and the second position designates the phoneme that follows the given phoneme. Senones are classified within like states (<i>e</i>.<i>g</i>., state 1) for each type of model (<i>e</i>.<i>g</i>., context-dependent vs. context-independent) corresponding to the same phoneme. In this example, the context-independent HMM has senones 10, 55, and 125 corresponding to states 1, 2 and 3 respectively.

0035Figure 7B shows an example of a corresponding context-dependent HMMs for the phoneme /aa/. In Figure 7B there are five context-dependent models which model the phoneme /aa/ in five distinct phonemic contexts (115-119). For example, the context-dependent model /aa/ (/dh/, /b/) 115, models the phoneme /aa/ in a context where the left or preceding phoneme is /dh/ and where the phoneme /b/ succeeds or is to the right of it. The senones are classified within like states in the different HMMs. In state 1, there are two context-dependent senones, denoted as senones 14 and 25. Overall, for the phoneme /aa/, there are 2 context-dependent senones 14 and 35 and 1 context-independcnt senone 10 for state 1; 2 context-dependent senones 25 and 85 and 1 context-independent senone 55 for state 2; and 1 context-dependent senone 99 and context-independent senone 125 for state 3.

0036Accordingly the phoneme-based continuous density HMM used in the preferred embodiment can be characterized by the following mathematical definition: <ul id="ul0001" list-style="none" compact="compact"><li>(1) <i>N</i>, the number of states in the model: preferably, three states are employed. However, the invention is not restricted to three, as many as five may alternatively be used.</li><li>(2) <i>M</i>, the number of mixtures in the output pdf</li><li>(3) <i>A</i> = {<i>a</i>_<i>ij</i>}, the state transition probability distribution, from state <i>i</i> to state <i>j</i>.</li><li>(4) <i>B</i> = {<i>b</i>_<i>i</i> (<i>x</i>)}, output probability distribution, the probability of emitting feature vector <i>x</i> when in state <i>i</i>, where<maths id="math0001" num=""><img file="EP0805434B1_D0001.tif" /></maths>    where <i>N</i> (<i>x</i>, <i>µ</i>_<i>k</i><i>,V</i>_<i>k</i>) denotes the multi-dimensional Gaussian density function defined by mean vector µ_<i>k</i> and covariance matrix <i>V</i>_<i>k</i>, The number <i>M</i> of mixture-components is typically, anywhere from 1 to 50; and <i>c</i>_<i>k</i> is the weight for the <i>k</i>th mixture component in state <i>l</i>. The output probability distribution associated with each state <i>i</i>, is represented by senone, <i>sd</i>_<i>i</i> and can be denoted as <i>p</i>(<i>x</i>_<i>l</i><i> sd</i>_<i>l</i>). </li><li>(5) π = {π_<i>l</i>}, the initial state distribution.</li></ul>

0037For convenicnce, the compact notation α = (<i>A</i>, <i>B</i>, π) is used to denote the complete parameter set of the model which is otherwise known as the parameter space of a HMM.

0038In step 48 of Figure 2, the parameters for the senone, context-dependent HMM, and the context-independent HMM are estimated. The training phase of a HMM consists of estimating these parameters using the training data, a text of the speech 22, and a dictionary of phonemic spellings of words 24. The output and transition probabilities can be estimated by the well-known Baum-Welch or forward-backward algorithm. The Baum-Welch algorithm is preferred since it makes better use of training data. It is described in Huang et al., <i>Hidden Markov Models For Speech Recognition,</i> Edinburgh University Press, 1990. However, this invention is not limited to this particular training algorithm, others may be utilized. Normally, about five iterations through the training data can be made to obtain good estimates of the parameters.

0039In step 50 of Figure 2, weighting or interpolation factors for each content-dependent senone are generated and are denoted by the mathematical symbol, λ. The weighting factors will be used to interpolate the output probabilities of the context-independent HMM with the output probabilities of the context-dependent HMM. The weighting factors indicate the adeptness of the context-dependent output pdf at predicting unseen data. The output pdf is estimated with training data and will closely predict data which resembles the training data. However, it is impossible to estimate the output pdf with training data that represents every possible input speech utterance, or with sufficient training data for it to predict all unseen data correctly. The role of the weighting factor is to indicate the adeptness of the output pdf to predict unseen data which is a function of the training data used to estimate the context-dependent and context-independent models. As the amount of training data for the context-dependent models gets large, λ will approach 1.0 and the output pdf will be heavily weighed. With a small amount of training data for the context-dependent model, λ will approach 0.0 and the output pdf will be weighed less. The optimal value for λ for each context-dependent senone is determined by deleted interpolation.

0040Briefly, the technique of deleted interpolation partitions the training data into two distinct sets. One set is used to estimate the parameters of the model and a second set is used to determine the weighting factor which indicates how well the output pdf can predict unseen training data. This process is iterative where at each iteration the different sets are rotated and a new model and weighting factor is produced. At the end of all the iterations, the average value of the weighting factor is calculated and used in the recognition phase

0041Figures 3-6 illustrate the steps used in computing the weighting factors. Referring to Figure 3, the training data is partitioned into <i>K</i> blocks in step 60. Preferably, there are two blocks of data. However, the invention is not limited to this number of blocks, others may be used dependent on the constraints of training data storage and training time.

0042A weighting factor is calculated for each context-dependent senone (step 62) by first finding <i>sen</i>_<i>Sl</i>, which is the context-independent senone that corresponds to <i>sen</i>_<i>SD</i> (<i>i</i>.<i>e</i>., the context-dependent senone) using the senone table (step 63). The calculation is derived through an iterative process, step 64, which converges when the difference between the new value of λ, denoted as λ_<i>new</i>_, meets a certain threshold. Preferably, the process converges or finishes when | λ - λ_<i>new</i> | &lt; .0001. The process commences by selecting an initial value for λ, step 66. Initially, for the first iteration of a senone, an initial value is preselected by the user. Preferably, the initial value can be an estimated guess as 0.8. For all other iterations, the initial value will be the previously calculated new value, λ = λ_<i>new</i>. In step 68, the process iterates <i>K</i> times. At each iteration, one block of data is selected as the deleted block, and the selected deleted block is one that was not chosen previously, step 70.

0043The process then proceeds to estimate the output probabilities for each context-dependent (denoted as <i>b</i>₁) and context-independent (denoted as <i>b</i>₂) senone using the training data from all the blocks except the deleted block (step 72). These parameters are estimated using the same technique as described above in reference to the estimation of the parameters of the HMMs in the training phase (<i>i</i>.<i>e</i>., Baum-Welch algorithm).

0044Next in step 74, a new value, λ_<i>new</i>_, is computed The computation assumes that "forced alignment" is sued. During training, if the Viterbi algorithm is used, each feature vector in the training data can be identified with a specific senone. This mapping or vectors with senone is known as "forced alignment." λ_<i>new</i>_, is calculated in accord with the following mathematical definition:<maths id="math0002" num=""><img file="EP0805434B1_D0002.tif" /></maths>    where <ul id="ul0002" list-style="none" compact="compact"><li><i>N</i> = number of data points or feature vectors in the deleted block that corresponds to senone <i>sen</i>_<i>SD</i> using forced alignment</li><li><i>x</i>_i = feature vector <i>i</i>, 1 ≤ <i>i</i> ≤ <i>N</i></li><li><i>b</i>₁ (<i>x</i>_<i>i</i>) <i>=</i> context-dependent output pdf as defined by equation (1) above</li><li><i>b</i>₂ (<i>x</i>_<i>i</i>) = context-independent output pdf as defined by equation (1) above</li><li>λ*<i>b</i>₁(<i>x</i>_<i>i</i>) + (1-λ)*<i>b</i>₂(<i>x</i>_<i>i</i>): referred to as the overall probability.</li></ul>

0045A value of λ_<i>new</i> is determined for each of the <i>K</i> iterations. At the completion of all <i>K</i> iterations, in step 76, an average value is computed which can be represented by the following mathematical expression:<maths id="math0003" num=""><img file="EP0805434B1_D0003.tif" /></maths>    where <ul id="ul0003" list-style="none" compact="compact"><li><i>j</i> = index of deleted block;</li><li><i>K</i> = number of blocks;</li><li>λ<maths id="math0004" num=""><math display="inline"><mrow><mfrac linethickness="0" numalign="left" denomalign="left"><mrow><mtext mathvariant="italic">j</mtext></mrow><mrow><mtext mathvariant="italic">new</mtext></mrow></mfrac></mrow></math><img file="EP0805434B1_D0004.tif" /></maths> = estimate of λ using deleted block <i>j</i>; and</li><li><i>N</i>_<i>j</i><i> =</i> number of points in deleted block <i>j</i> that correspond to <i>sen</i>_<i>SD</i> using forced alignment.</li></ul>

0046Steps 66 through 76 are executed again if the value of λ_<i>new</i> does not meet the prescribed threshold. When the process converges for a particular context-dependent senone, the current value of λ_<i>new</i> is stored in lambda table 26 for the particular context-dependent senone.

0047Figure 4 depicts a flowchart of the steps used in computing the new value for the weighting factor, λ_<i>new</i> in accord with equations (2) and (3) above. The new value is computed by summing the contribution of the context-dependent output pdfrelative to the overall probability for each data point in the deleted block. Thus, in step 79, all points in the deleted block that correspond to <i>sen</i>_<i>SD</i> are found using the model generated in step 48 and forced alignment. In step 80, the process iterates for each data point <i>x</i>_<i>i</i> in the deleted block that is aligned with <i>sen</i>_<i>SD</i>. The contribution of the context-dependent output pdf for data point <i>x</i>_<i>i</i> relative to the overall probability is determined in step 82 in accord with the following mathematical definition:<maths id="math0005" num="(4)"><math display="block"><mrow><mfrac><mrow><mtext>(λ</mtext><msub><mrow><mtext mathvariant="italic">*b</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>))</mtext></mrow><mrow><mtext>(λ*</mtext><msub><mrow><mtext mathvariant="italic">b</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>) + (1 - λ)*</mtext><msub><mrow><mtext mathvariant="italic">b</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>))</mtext></mrow></mfrac></mrow></math><img file="EP0805434B1_D0005.tif" /></maths>

0048The sum of the contributions for all data points computed thus far are totaled in step 84. At the completion of the iteration when all the data points in the deleted block that are aligned with <i>sen</i>_<i>SD</i> have been processed, the average of the contributions is computed, λ_<i>new</i>, step 86 in accord with equation (2) above.

0049The above computation of the weighting factor uses the data points in the deleted block. This produces a more accurate computation at the expense of increasing the training time as well as the amount of storage needed by the training engine to perform the computations. In some instances, it may be more advantageous to generate a parametric representation of the data points in the deleted block that correspond and to use the appropriate parameters instead, Another alternative is to use reconstructed data points from the parametric representation of the data points to <i>sen</i>_<i>SD</i>. These alternatives provide a coarse approximation of the data points but has the advantage of computational efficiency.

0050Figures 5 and 6 depict these alternate embodiments for the calculation of the weighting factors. Figure 5 depicts the first alternate embodiment. Referring to Figure 5, a parameteric representation for the data points in the deleted block is generated as shown in step 90. In this case, the parametric representation is a mixture of Gaussians. This representation can be made using the Baum-Welch algorithm as described above. The parameters generated include the mean, µ_<i>j</i>, and weight, <i>c</i>_<i>j</i>, for each mixture component <i>j</i>. The computation of the new value of lambda, λ_<i>new</i><i>,</i> can be made in accord with the following mathematical definition for a deleted block µ_<i>j</i>:<maths id="math0006" num=""><img file="EP0805434B1_D0006.tif" /></maths>    where <dl id="dl0001" compact="compact"><dt><i>M</i></dt><dd>= number of normal mixture components;</dd><dt><i>C</i>_<i>j</i></dt><dd>= weight of the <i>jth</i> normal mixture component; note that<maths id="math0007" num=""><img file="EP0805434B1_D0007.tif" /></maths></dd><dt>µ_<i>j</i> =</dt><dd>mean of the <i>j</i>th normal mixture component;</dd></dl>

0051Steps 92 - 98 perform this calculation in the following manner. Step 92 iterates for each mixture and determines the contribution of the context-dependent output probability over the overall probability for the mixture having the corresponding mean and weight parameters. For mixture components, this is represented mathematically as:<maths id="math0008" num="(6)"><math display="block"><mrow><mfrac><mrow><msub><mrow><mtext mathvariant="italic">c</mtext></mrow><mrow><mtext mathvariant="italic">j</mtext></mrow></msub><mtext>*(λ*</mtext><msub><mrow><mtext mathvariant="italic">b</mtext></mrow><mrow><mtext mathvariant="italic">1</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">j</mtext></mrow></msub><mtext>)</mtext></mrow><mrow><msub><mrow><mtext mathvariant="italic">λ*b</mtext></mrow><mrow><mtext mathvariant="italic">1</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">µ</mtext></mrow><mrow><mtext mathvariant="italic">j</mtext></mrow></msub><mtext>)</mtext><mtext mathvariant="italic">+</mtext><mtext>(</mtext><mtext mathvariant="italic">l-λ</mtext><mtext>)*</mtext><msub><mrow><mtext mathvariant="italic">b</mtext></mrow><mrow><mtext mathvariant="italic">2</mtext></mrow></msub><msub><mrow><mtext>(µ</mtext></mrow><mrow><mtext mathvariant="italic">j</mtext></mrow></msub><mtext>)</mtext></mrow></mfrac></mrow></math><img file="EP0805434B1_D0008.tif" /></maths>

0052In step 96, the sum of these contributions is formed for all the mixture components. In step 98, the final sum resulting from step 96 is stored as the value of λ_<i>new</i> for the current <i>sen</i>_<i>SD</i> and the deleted block.

0053Referring to Figure 3, at the completion of the <i>K</i> iterations, the process proceeds to calculate the average value for λ_<i>new</i> in step 76 in accord with equation (3) above. The process continues as described above with reference to Figure 3 until the process converges and the current average value λ_<i>new</i> is stored in the lambda table 26 for the particular context-dependent senone.

0054In the second alternative embodiment of the calculation of the weighting factors, a select number of data points are used, which are randomly generated from a parametric representation for the senone. Figure 6 depicts this second alternate embodiment which can be described mathematically for a deleted block according to equation (2) set forth above, except that {<i>x</i>_<i>i</i>} = generated data points and <i>N</i> = the number of gencrated data points.

0055This alternative embodiment differs from the preferred embodiment, shown in Figure 3, in the determination of the new value of λ_<i>new</i><i>,</i> (step 74). The flow sequence remains as shown in Figure 3. Referring to Figure 6, in step 100, a parametric representation is generated for the data points in the deleted block. The parametric representation can consist of a mixture of Gaussians. This parametric representation can be derived using the Baum-Welch algorithm on the training data in the deleted block. From this parametric representation, a prescribed number of data points is reconstructed using a random number generator with the mean and weight parameters, as shown in step 102. The number of data points that are reconstructed is a trade-off between the desired accuracy of λ_<i>new</i> and the computational requirements. A higher number of data points improves the accuracy of the λ_<i>new</i> at the cost cf greater computational requirements. A suitable number of reconstructed data points per mixture is 100.

0056In step 104, steps 106 and 108 are performed for each data point in the set, step 104. In step 106, the contribution of the context-dependent output probability over the overall probability for the data point is determined. This can be represented mathematically as:<maths id="math0009" num="(9)"><math display="block"><mrow><mfrac><mrow><mtext>(λ*</mtext><msub><mrow><mtext mathvariant="italic">b</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>)) </mtext></mrow><mrow><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">λ*b</mtext></mrow><mrow><mtext>1</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>) + (1 - </mtext><mtext mathvariant="italic">λ</mtext><mtext>)*</mtext><msub><mrow><mtext mathvariant="italic">b</mtext></mrow><mrow><mtext>2</mtext></mrow></msub><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">x</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>))</mtext></mrow></mfrac></mrow></math><img file="EP0805434B1_D0009.tif" /></maths> In step 108, the sum of these contributions is formed for all the data points in the set. At the completion of the iteration through all the data points in the set, the average of all of the contributions is returned as the value of λ_<i>new</i>, (step 110). Referring to Figure 3, at the completion of the <i>K</i> iterations, the process proceeds to calculate the average value for λ_<i>new</i> in step 76 in accord with equation (3) above. The process continues as described above with reference to Figure 3 until the process converges and the current average value <b>λ</b>_<i>new</i> is stored in the lambda table 26 for the particular context-dependent senone.

0057Once the training data has been generated and stored in the appropriate storage locations, the recognition system is ready to execute. The primary task of the speech recognition system is to detect the linguistic message which is embodied in the input speech signal. This task is a multi-level decoding problem since it requires matching the sequence of feature vectors to a sequence of phonemes, matching the sequence of phonemes to a sequence of words, and matching the sequence of words to a sentence. This is performed by forming all possible linguistic expressions that have been modeled and calculating the probability that the expression matches the sequence of feature vectors. Since a linguistic expression is composed of a sequence of phonemes, the determination can involve calculating the likelihood that the phonemes forming the expression match the feature vectors and that the expression is likely to occur (i.e., grammatically correct). The probability that the phonemes forming an expression match the feature vectors can be referred to as an acoustic score and the probability that the expression occurs can be referred to as the language score. The language score takes into consideration the syntax and semantics of the language, such as the grammar of the language, and indicates whether the sequence of words corresponding to the sequence of phonemes form a grammatically-correct linguistic expression

0058In the preferred embodiment, phonemes are represented by HMMs where the output pdf of like states are clustered forming senones. The process of matching a feature vector to a phoneme then entails matching a feature vector to the senones associated with the states of a HMM representing the phoneme. Thus, the linguistic expression can be composed of senones corresponding to states of a sequence of HMMs.

0059In the preferred embodiment of the invention, the task of the recognition engine can be one of finding the word sequence <i>W</i> which maximizes the probability <i>P</i>(<i>W</i>/<i>X</i>). The probability <i>P</i>(<i>W</i>/<i>X</i>) represents the probability of the linguistic expression <i>W</i> occurring given the input speech signal <i>X</i>. <i>W</i> can be a word string denoted as <i>W</i> = <i>w</i>₁, <i>w</i>₂, ..., <i>w</i>_<i>n</i>, where w_<i>i</i> denotes individual words, each word is represented by a sequence of phonemes, <i>w</i>₁ = <i>p</i>₁, <i>p</i>₂, .... <i>p</i>_<i>q</i>, and <i>X</i> is the input speech signal represented by a sequence of feature vectors, denoted as <i>X</i> = <i>x</i>₁, <i>x</i>₂, ..., <i>x</i>_<i>n</i>. This maximization problem can be solved using a modified version of the well-known Bayes formula which is described mathematically as:<maths id="math0010" num="(10)"><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">W</mtext><mtext>/</mtext><mtext mathvariant="italic">X</mtext><mtext>) = (</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">X</mtext><mtext>/</mtext><mtext mathvariant="italic">W</mtext><mtext>)</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">W</mtext><mtext>))/</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">X</mtext><mtext>).</mtext></mrow></math><img file="EP0805434B1_D0010.tif" /></maths><i>P</i>(<i>X</i>/<i>W</i>) is the probability that the input speech signal X matches the word string <i>W,</i> and is referred to as the acoustic score. <i>P</i>(<i>W</i>) is the probability that the word string <i>W</i> will occur, and is referred to as the language score. Since <i>P</i>(<i>X</i>) is independent of <i>W</i>, maximizing <i>P</i>(<i>W</i>/<i>X</i>) is equivalent to maximizing the numerator, namely <i>P</i>(<i>X</i>/<i>W</i>)<i>P</i>(<i>W</i>) over all word sequences <i>W</i>.

0060The recognition task considers various word sequences in trying to determine the best match. For each word sequence that is considered by the recognition task, it computes an acoustic score and a language score. A language score indicates how likely the word sequence is, in the language, and is indicated by the term <i>P</i>(<i>W</i>) in the above equation (10). An acoustic score indicates how well the sequence of acoustic vector features matches the acoustic model for the word sequence <i>W</i>. The acoustic score is indicated by the term <i>P</i>(<i>X</i>/<i>W</i>) in the above formula

0061In computing the acoustic score for a given word sequence, the recognition task considers various senone alignments. A senone alignment is a mapping from the sequence of acoustic feature vectors to senones which assigns a unique senone to each acoustic feature vector. Only the senone alignments which would result in the word sequence under consideration are considered by the recognition task. An acoustic score for the word sequence under the constraints of each senone alignment is computed. The acoustic score for the word sequence is the best acoustic score over all possible senone alignments.

0062Mathematically, this can be expressed as:<maths id="math0011" num=""><math display="block"><mrow><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">X</mtext><mtext>/</mtext><mtext mathvariant="italic">W</mtext><mtext>) </mtext><mtext mathvariant="italic">=</mtext><mtext> Max [</mtext><mtext mathvariant="italic">P</mtext><mtext>(</mtext><mtext mathvariant="italic">X</mtext><mtext>/(</mtext><msub><mrow><mtext mathvariant="italic">W, A</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>))]</mtext></mrow></math><img file="EP0805434B1_D0011.tif" /></maths>       over <i>i</i> = 1 to <i>q</i>    where <ul id="ul0004" list-style="none" compact="compact"><li><i>A</i>₁ ... <i>A</i>_<i>q</i>: all possible senone alignments for word sequence <i>W</i></li></ul>

0063The computation of the acoustic score for word sequence <i>W</i> under the constraint of a given senone alignment <i>A</i> can be further expressed as<maths id="math0012" num=""><img file="EP0805434B1_D0012.tif" /></maths> where senone alignment <i>A</i> aligns or maps the <i>i</i>th acoustic feature vector <i>x</i>_<i>i</i> to the context dependent senone <i>sd</i>_<i>i</i>. <i>P</i>(<i>A</i>) represents the state transition probability for the senone sequence <i>sd</i>_<i>i</i> ... <i><u>sd</u></i>_{<i><u>n</u></i>}<i> P</i>(<i>x</i>_<i>i</i>/<i>sd</i>_<i>i</i>) represents the probability that feature version <i>x</i>_<i>i</i> matches context-dependent senone <i>sd</i>_<i>l</i>.

0064The essence of the acoustic score is the computation of the output probability <i>p</i>(<i>x</i>|<i>sd</i>). This represents the likelihood that the feature vector, <i>x</i>, matches the senone, <i>sd</i>, which corresponds to a context-dependent HMM state. However, a poorly estimated output pdf can contribute to inaccuracies in the computation of the acoustic score. This usually occurs due to insufficient training data. The robustness of the distribution increases with the use of more training data to estimate output pdf.

0065One way to reduce this problem is to utilize several HMMs which model the same phenomes at several levels of detail. The output pdf for a particular state can then be constructed by using the output pdfs at various levels of detail and combining them. The combination is done based on the ability to predict data not seen during training. A robust output pdf which is more adept at predicting unseen data will receive a higher weight while a poorly estimated output pdf will receive a lower weight in the combined output pdf. In the preferred embodiment, several context-dependent HMMs and a context-independent HMM are utilized to model a phoneme. A weighting factor, λ, for each senone corresponding to a context-dependent state which was computed previously in the training phase is used to indicate the weight each senone is given. The larger λ is (approaches 1.0) the more the context-dependent senone dominates and the less the context-independent senone is weighed When λ is small (approaches 0.0), the context-independent senone dominates. Thus, the computation of the output probability, <i>p</i>(<i>x</i>|<i>sd</i>), can be represented by the following mathematical definition:<maths id="math0013" num="(12)"><math display="block"><mrow><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>|</mtext><mtext mathvariant="italic">sd</mtext><mtext>) = λ * </mtext><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>|</mtext><msub><mrow><mtext mathvariant="italic">sd</mtext></mrow><mrow><mtext mathvariant="italic">d</mtext></mrow></msub><mtext>) + (1 - λ) * </mtext><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>|</mtext><msub><mrow><mtext mathvariant="italic">sd</mtext></mrow><mrow><mtext mathvariant="italic">i</mtext></mrow></msub><mtext>),</mtext></mrow></math><img file="EP0805434B1_D0013.tif" /></maths>    where <ul id="ul0005" list-style="none" compact="compact"><li>λ is the weighting factor between 0 and I for senone <i>sd</i>;</li><li><i>x</i> is the feature vector;</li><li><i>sd</i>_<i>d</i> is the senone associated with a state of a context-dependent HMM;</li><li><i>sd</i>_<i>i</i> is the senone associated with the corresponding state of a context-independent HMM;</li><li><i>p</i>(<i>x</i>|<i>sd</i>_<i>d</i>) is the probability of the feature vector x matching senone <i>sd</i>_<i>d</i>; and</li><li><i>p</i>(<i>x</i>|<i>sd</i>_<i>i</i>) is the probability of the feature vector x matching senone <i>sd</i>_<i>i</i>.</li></ul>

0066Thus, the output probability, <i>p</i>(<i>x</i>|<i>sd</i>), is linearly interpolated as a function of the output probabilities of context-dependent and context-independent scnoncs. The weighting or interpolation factor λ indicates the degree to which each senone is interpolated.

0067Figure 8 depicts the operation of the speech recognition method. Referring to Figure 8, the method commences by receiving an input speech utterance, (step 122), which is converted to feature vectors (step 124), which was previously detailed above with reference to Figure 1. In step 126, the method performs steps 128-136 for each word sequence that can represent the input speech utterance. The word sequence can consist of a variety of different senone alignments, where each senone alignment corresponds to a sequence of HMM states. In steps 128-134 a combined recognition score for each possible senone alignment which can represent the word sequence is determined. The combined recognition score can be determined in accord with the modified Bayes formula as denoted above in equation (10). The combined recognition score consists of an acoustic score and a language score. The acoustic score is determined in step 130, the language score is determined in step 132, and the combined score is computed in step 134. The senone alignment having the highest combined recognition score is then selected to represent the word sequence, step 136. In step 138, the method recognizes the input speech utterance as the word sequence having the highest combined recognition score.

0068In step 130, the acoustic score can be determined as described above in accord with equation (11) where the output probability is computed as described above in equation (12).

0069In step 132, the method computes a language score based on the language models representing linguistic expressions stored in language model storage 32. Language models use knowledge of the structure and semantics of a language in predicting the likelihood of the occurrence of a word considering the words that have been previously uttered. The language model can be a bigram language model where the language score is based on the probability of one word being followed by a particular second word. Alternatively, the language model may be based on N-grams other than bigrams or each on subword language probabilities. In addition, other lexical knowledge such as syntax and grammatical rules can be employed to create the language model. Methods for creating and using language models are well-known in the art and are described in more detail in the Huang et al. book referred to above.

0070The above detailed invention improves the recognition capability of a speech recognition system by utilizing multiple continuous density output probabilities corresponding to the same speech event in different contexts. This improves the mapping of the feature vectors to the hidden Markov models since it improves the model's performance in predicting speech events that the model was not trained with. An improvement at this level is extremely beneficial since the mapping at this level is the foundation upon which the recognition process further builds on.

0071However, it should be noted that this invention is not restricted to a speech recognition system. Any application which requires the matching of a speech utterance to a linguistic expression can utilize the claimed invention. The speech utterance can be any form of acoustic data, such as but not limited to, sounds, speech waveforms, and the like. An example of such an application is a speech synthesis system which utilizes probabilistic models to generate a speech waveform from a text string representing a linguistic expression.

0072Although the preferred embodiment of the invention has been described hereinabove in detail, it is desired to emphasize that this is for the purpose of illustrating the invention and thereby to enable those skilled in this art to adapt the invention to various different applications requiring modifications to the apparatus described hereinabove; thus, the specific details of the disclosures herein are not intended to be necessary limitations on the scope of the present invention, which is defined by the appended claims.