RAPIDLY TRAINING A SPEECH RECOGNIZER TO A SUBSEQUENT SPEAKER GIVEN TRAINING DATA OF A REFERENCE SPEAKER

Abstract: Apparatus and method for training the statistics of a Markov Model speech recognizer to a subsequent speaker who utters part of a training text after the recognizer has been trained for the statistics of a reference speaker who utters a full training text. Where labels generated by an acoustic processor in response to uttered speech serve as outputs for Markov models, the present apparatus and method determine label output probabilities at transitions in the Markov models corresponding to the subsequent speaker where there is sparse training data. Specifically, label output probabilities for the subsequent speaker are re-parameterized based on confusion matrix entries having values indicative of the similarity between an ℓth label output of the subsequent speaker and a kth label output for the reference speaker. The label output probabilities based on re-parameterized data are combined with initialized label output probabilities to form "smoothed" label output probabilities which feature smoothed probability distributions. Based on label outputs generated when the subsequent speaker utters the shortened training text, "basic" label output probabilities computed by conventional methodology are linearly averaged against the smoothed label output probabilities to produce improved label output probabilities.

Description: The present invention relates to apparatus and method for training the statistics of a Markov Model speech recognizer to a subsequent speaker after the recognizer has been trained for a reference speaker. One approach to speech recognition involves the use of Hidden Markov Models (HMM). Hidden Markov Models have been discussed in various articles such as: "Continuous Speech Recognition by Statistical Methods" by F. Jelinek, Proceedings of the IEEE, volume 64, number 4, 1976 and "A Maximum Likelihood Approach to Continuous Speech Recognition", by L.R. Bahl, F. Jelinek, and R.L. Mercer IEEE Transactions on Pattern Analysis and Machine Intelligence, volume PAMI-5, Number 2, March 1983. These articles are incorporated herein by reference. In performing speech recognition based on Hidden Markov Models, successive intervals of speech are examined by an acoustic processor with respect to various predefined characteristics of speech. For example, respective amplitudes for each of various energy frequency bands are determined for each time interval. Each respective amplitude represents a component, or feature. Together, the components combine to form a feature vector. The acoustic processor defines a finite set of prototype, or reference, vectors. Each prototype vector has a unique label which identifies it. The feature vector at each successive time interval is compared with each prototype vector. Based on a prescribed distance measure, the closest prototype is selected. Hence, for each time interval a prototype vector (which most closely represents the feature vector of the interval) is selected. As speech is uttered, the acoustic processor provides as output a string of labels. In accordance with Markov models speech recognition, a set of Markov models is defined. Typically, such Markov models have corresponded one-to-one with phonetic elements. For eighty phonetic elements, then, there are eighty respective Markov models. The Markov models corresponding to successive phonetic elements of a word can be concatenated to form a Markov model baseform for the word. Each Markov model is characterized as having a plurality of states and a plurality of transitions. Each transition extends from a state to a state. At least some transitions represent a time interval during which a prototype vector label is selected by the acoustic processor. For each transition, there is a transition probability and, in some cases, output probabilities. Typically associated, the transition probabilities indicate the likelihood of following a given transition in a Markov model. The output probabilities indicate the likelihood of a certain output label (e.g., prototype vector label) being produced during a given transition. For a certain transition A ij extending from state i to state j, there is an associated transition probability P(A ij ) and, where there are 200 different prototype vectors, there are 200 associated output probabilities: Normally, but not necessarily, the skeletal structure of states with connecting transitions (without probability values assigned) is the same for each Markov model. For a given speaker, the various respective Markov models for the different phonetic elements differ typically in the values of the probabilities associated therewith. In order to be operative, the various transition probabilities and output probabilities for each Markov model must be determined. The physical implementation of the Markov model is referred to as a "phone machine", or Markov model phone machine. The phone machine for a corresponding phonetic element includes memory locations for storing the transition probabilities, output probabilities, shape of the phone machine, identifiers indicating which phonetic element is represented thereby, and other such information which characterizes the respective Markov model. The process of determining the transition probabilities and output probabilities so that they may be stored for phone machines is referred to as "training." Typically, a distinct set of transition probabilities and output probabilities must be determined for each speaker. That is, for each speaker, the speech recognizer stores data (e.g., transition probability values and output probability values) for a respective set of phone machines. The conventional approach to training is for a speaker to utter a known sample text into an acoustic processor. The sample text represents a known sequence of phonetic elements and, hence, a known corresponding sequence of phone machines. The acoustic processor generates a string of prototype labels in response to the uttered speech input. From the string of prototype labels generated for the known sample text and from initially set values of the transition probabilities and output probabilities (which may not reflect actual speech characteristics), improved probability values can be determined by applying a forward-backward algorithm, or Baum-Welch algorithm, to produce transition counts and output counts, deriving transition probabilities and output probabilities therefrom, applying the forward-backward algorithm with the derived probabilities to produce updated counts, and so on over a number of iterations. The probability values after the last iteration are referred to herein as "basic" transition probabilities and "basic" output probabilities In order to generate reasonably accurate "basic" probabilities, it is necessary for a speaker to utter a relatively long sample text, extending, for example 20 minutes. In accordance with prior technology, each speaker would be required to utter the 20-minute sample text in order to train the speech recognizer to his/her speech. A required training period of 20 minutes per speaker may be undesirably long and inconvenient. Also however, the amount of computing the speech recognizer must perform for 20 minutes of training text in accordance with the forward-backward algorithm in order to determine "basic" probabilities is excessive. Accordingly, a significant problem in speaker dependent Markov model speech recognition has involved the lengthy period during which each speaker must utter text and the computationally costly process of applying the forward-backward algorithm to the full text for each speaker. In the US Patent 4,718,094 entitled "Speech Recognition System", some of the transitions are grouped together to have common output probabilities applied thereto. Although reducing the required amount of training data, the sample text has remained nonetheless lengthy when multiple speakers are to be recognized. In the article "Speaker Adaptation for a Hidden Markov Model", published in the ICASSP 86 Proceedings - IEEE - IECEJ - ASJ International Conference on Acoustics, Speech, and Signal Processing, Tokyo, 7th - 11th April 1986, vol. 4, pages 2667 - 2670, IEEE, New York, US; by the authors K. Sugawara et al., a method for adapting Markov models to a subsequent speaker is disclosed. During the training process, parameters of an hidden Markov model are calculated iteratively using "forward-backward algorithm". The adaptation method uses the intermediate results of the last iteration. These values are used to reconstruct the state (probabilistic counts at each node) of the last iteration of the training process. Then the newly added training data for adaptation can be processed as the original data. The parameters are calculated as the weighted sum of the restored probabilistic counts and the ones obtained using new training data. Adaptation can be done on the word basis. This means that the user can select an error-prone word and retrain or adapt its corresponding hidden Markov model. SUMMARY OF THE INVENTION It is an object of the present invention to provide apparatus and method for rapidly training a Hidden Markov Model speech recognizer to subsequent speakers after the speech recognizer has been trained to a reference speaker. That is, after transition probabilities and output probabilities for a reference speaker have been determined, training statistics for a subsequent speaker are determined. According to the present invention, the time each subsequent speaker must utter text during training is greatly reduced. In addition, the level of computation required to determine the transition probabilities and label output probabilities for a subsequent speaker is also reduced. The present invention achieves the above objects by using (a) data determined during the training of the speech recognizer to a reference speaker and (b) similarity relations between data generated for the reference speaker and data generated for the subsequent speaker during utterance of the shortened sample text. The computational cost in employing the reference speaker data and similarity relation data is considerably less than applying the forward-backward algorithm to prototype outputs generated in response to utterance of the full training text. In accordance with the invention, it is presumed that, for a subsequent speaker, each transition A ij from state i to state j has a corresponding transition probability P(A ij ) that is to be computed. There is also allocated, for each transition at which a label is produced, a number of label output probabilities to be computed. Where there are 200 distinct labels, there are normally 200 different label output probabilities at each non-null transition. For each transition from a state i to a state j, there is then (1) a transition probability P(A ij ) therefor and (2) for each non-null transition, a probability P(F₂(ℓ)|A ij ) of generating each ℓth label given the transition from state i to state j (of the given model), where F₂(ℓ) represents the ℓth label for a subsequent speaker. The present invention reduces the data required to find the numerous probabilities P(F₂(ℓ)|A ij ) by re-structuring the transition from state i to state j as a plurality of parallel transitions wherein the probability of each of the parallel transitions is of the form: P(F₁(k)|A ij ) P(F₂(ℓ)|F₁(k)) where k identifies one of the N labels for the reference speaker. Taken collectively over all k values, the label output probability for the ℓth label output at transition A ij is: In examining the above equation, it is observed that the P(F₁(k)|A ij ) component corresponds to previously determined data for the reference speaker and that the P(F₂(ℓ)|F₁(k)) component corresponds to a "confusion matrix" which defines a relationship between the labels of the reference speaker and a subsequent speaker. The confusion matrix is an NxN matrix (e.g., a 200x200 matrix for an embodiment in which there are 200 different labels). Each entry represents the probability of a label ℓ spoken by the subsequent speaker corresponding to a label k spoken by the reference speaker. In this expanded model with parallel transitions, the probability of each transition corresponds to P(F₁(k)|A ij ) and the output probabilities correspond to confusion matrix values P(F₂(ℓ)|F₁(k)). From confusion matrix entries and other information, updated re-parameterized label output probabilities are computed. The updated re-parameterized label output probabilities are linearly averaged with initialized label output probabilities in a deleted estimation processor to form smoothed relational probabilities. Basic parameter label output probabilities, determined in a conventional forward-backward algorithm manner, are then linearly averaged with the smoothed relational label output probabilities to form final label output probabilities for the subsequent speaker. The final label output probabilities for the subsequent speaker are generated with the subsequent speaker uttering only approximately 1/3 of the full training text; however, the accuracy achieved with a Markov model speech recognizer employing the final label output probabilities is comparable to the accuracy achieved when full text training is used. Moreover, the overall computational cost in computing the final label output probabilities is considerably less than that required for full training text training. BRIEF DESCRIPTION OF THE DRAWINGS Fig. 1 is a general block diagram of a speech recognizer. Fig. 2 is a graph illustrating an embodiment of a phonetic Markov model. Fig. 3 is a diagram showing two Markov models linked together to form a baseform for the word "THE". Fig.4 is an embodiment of a "phone machine" which stores information relating to a respective Markov model from a predefined set of Markov models. Fig. 5 is an illustration of a stored matrix which contains probabilities of labels for 240 different arc types. Fig. 6 is a graph illustrating an embodiment of a "fenemic" Markov model. Fig. 7 is an illustration showing a trellis representation of concatenated fenemic Markov models. Fig. 8 is a general block diagram of a training apparatus according to the present invention. Fig. 9 is a block diagram of a probability initializer of Fig. 8. Fig. 10 is an illustration of an original Markov model graph between two states i and j. Fig. 11 is an illustration of an expanded Markov model graph between two states i and j, wherein parallel branches correspond to re-parameterized data. Fig. 12 is a block diagram of a re-parameterized forward-backward processor of Fig. 8. Fig. 13 is an illustration of a stored confusion matrix. Fig. 14 is a block diagram of a preferred basic parameter forward-backward processor of Fig. 8. Fig. 15 is a block diagram of a preferred deleted estimation processor, such those of Fig. 8. Fig. 16 is a flowchart diagram of a training methodology according to the present invention. A. Speech Recognition System Environment with Phonetic Phone Machines Referring to Fig. 1, a basic speech recognition system 100 is shown to include an acoustic processor 102, a stack decoder 104, an acoustic match processor 106, and a language model processor 108. Briefly, the acoustic processor 102 defines an alphabet of speech types based on some prescribed characteristics. Each characteristic preferably relates to a value corresponding to an energy frequency band. In the European Patent Application 85 111 905.7 entitled "Nonlinear Signal Processing in a Speech Recognition", each characteristic corresponds to a "neural firing rate" ascribed to a corresponding frequency band of energy. Defining all speech as a spectral space, speech types represent regions in the space. Each speech type covers a distinct cluster of sounds in the space. Each component of the space corresponds to one of the prescribed characteristics. Hence, given values for each characteristic, a point or vector in the space is specified. Each point or vector, in turn, is contained in one of the clusters. A convenient method of determining the appropriate cluster for a given point or vector involves the use of prototypes. For each cluster, a representative prototype vector is defined. In determining which cluster (i.e., speech type) a point or vector belongs in, the distance between the subject point and each prototype vector is measured in a known manner. The subject point belongs to the speech type having the "closest" prototype vector. For one interval of speech after another, respective values for the prescribed characteristics are combined to represent a "feature vector". For each time interval, a determination is made as to which speech type best represents the feature vector generated during the interval. Preferably, the acoustic processor 102 defines 200 different speech types, each having a respective prototype vector. Each speech type is identified by a corresponding numerical identifier, or label. Table 1 includes a preferable list of 200 speech types having labels AA11 through ZX3-. By conforming the labels to sounds represented thereby (AA11 to an "A" sound, ZX3- to a "Z" sound, and so on), human interpretation is facilitated. A more general scheme of labels identifies the 200 labels simply as L1 through L200. In processing successive intervals of a speech waveform, the acoustic processor 102 generates a string of labels y₁y₂y₃...y T . Each label in the string is selected from the alphabet of 200 labels, based on which corresponding speech type is determined the closest for a given interval t. The prior art discloses a variety of techniques for clustering speech into speech types and representing each speech type by a prototype; and then selecting the closest prototype for each successive time interval. The label string generated by acoustic processor 102 enters the stack decoder 104 and the acoustic match processor 106. The purpose of the stack decoder 104 is to determine which "path" of words is a likely path for the input string of labels, based on acoustic match and language model inputs. The language model processor 108 is based on word contexts. In particular, a collection of words has been previously examined to determine the frequency of three-word (i.e., tri-gram) and two-word (bi-gram) sequences, as well as single word (uni-gram) occurrences. Based on the tri-gram, bi-gram, and uni-gram data, the language model processor 108 provides input to the stack decoder 104 regarding the likelihood of a specific word following two preceding words. The other input to the stack decoder 104 is from the acoustic match processor 106. The acoustic match processor 106 produces an acoustic match score for each word in the vocabulary. The acoustic match score is determined based on Markov models implemented in the form of "phone" machines. The structure of a Markov model 120 is illustrated in Fig. 2. The Markov model 120 of Fig. 2 is a phonetic model in that it corresponds to a phonetic element. Phonetic elements (e.g., elements in the International Phonetic Alphabet) are well-known and in common use. The phonetic Markov model 120 in Fig. 2 is represented as a plurality of states and transitions, each transition extending from a state to a state. Each transition may be viewed as corresponding to a time interval during which a label can be generated by the acoustic processor 102 (of Fig. 1). At successive time intervals, speech may be viewed as being produced while passing from a state to a state along a transition. The phonetic Markov model 120 is illustrated with seven states S1 through S7 and thirteen transitions. From the initial state S1, there are three possible paths. The first extends along a transition tr₁ to state S4. The second extends along a transition tr₈ to state S2. And the third path from state S1 extends along transition tr₁₁ to a final state S7. From state S2, transition tr₉ and tr₁₂ may alternatively be followed to states S3 and S7 respectively. From state S3, transition tr₁₀ or tr₁₃ may alternatively be followed to state S7. Transitions tr₁₁, tr₁₂ and tr₁₃ are each represented by a dashed line and are referred to as "null transitions." A null transition corresponds to the departure from one state without outputting a label. States S4, S5, and S6 each have a first transition which extends to a next state and a second transition which "self-loops" back to the original state. With self-looping, one label output after another can be produced at successive time intervals without advancing to a "new" state. Null transitions and self-loop transitions provide flexibility to the Markov model of Fig. 2 and account for variations which may occur in pronunciation at different times or under different conditions. Preferably, in the phonetic environment, each phonetic element is represented by respective phonetic Markov models. Words are represented as word baseforms by linking together Markov models of successive phonetic elements. In Fig.3, one form of the word THE is spelled phonetically as DH-UH1. A phonetic Markov model for DH is shown followed by a phonetic Markov model for UH1 The word "THE" is thus represented as a sequence of two phonetic Markov models. The respective phonetic Markov models are shown to have similar general structures (i.e., seven states and thirteen transitions). However, referring again to the Markov model 120, it is observed that each transition tr₁ through tr₁₃ has a probability associated therewith and that, at each non- null transition, there is allocated an output probability for each label. For each transition, tr₁ through tr₁₀ there are two hundred probabilities; one probability indicating the likelihood of a label L1 (first label of the alphabet) being produced at a given transition, one probability of a label L2 (the second label of the alphabet being produced at the given transition), ..., and one probability indicating the likelihood of a label L200 (the 200th label of the alphabet) being produced at the given transition. This is shown in Fig. 2 as If there are 80 phonetic elements, there is a corresponding set of 80 different Markov models (for a given speaker). Each Markov model in the set has distinct transition probabilities and label output probabilities. The probabilities are determined during a training period (discussed in detail hereinbelow) and are stored in memory in the acoustic match processor 106 (of Fig. 1). The implementation of a Markov model, with probabilities and other related data stored therefor, is referred to as a "phone machine". One simple method of storing the probabilities for a phone machine is shown in Fig. 4. Letting the phone machine for the phonetic element DH correspond to the 3rd phone machine of the set, the probability each of the 13 transitions of the DH phone machine are stored in a memory location identified as p₃[tr i ] where 1≦i≦13. The label output probabilities for each non-null transition are stored in locations identified as p₃[1] through p 1 / 3[200]. Similar storage is provided for the phone machine corresponding to the phonetic element UH1. If UH1 corresponds to the 75th phone machine in the set. all UH1 probabilities are identified by a 75 subscript. For 80 phonetic phone machines for Markov models like that of Fig. 2, the following probabilities are ordinarily determined and stored: 200 [distinct labels generatable] x10 [non-null transitions per phone machine] x80 [distinct phone machines] 160,000 [label probabilities] 13 [transitions per phone] x80 [distinct phones] 1040 [transition probabilities] There are thus a total of 160,000+1040=161,040 probabilities that are known (or estimated). To reduce the computational time required to generate all the probabilities and to reduce the required storage, several methods have been implemented. According to one method, it is assumed that the non-null transitions of each phone machine may be grouped into three sections. The first section includes transitions 1, 2, and 8; the second includes transitions 3, 4, 5, and 9; and the third includes transitions 6, 7, and 10. It is also assumed that, although each transition maintains its respective probability, the label output probabilities at the various transitions in a section are the same. That is, for transitions 1, 2, and 8, there is a single common label output probability distribution. Similarly, the transitions in the second section have a common label output probability distribution and the same for the transitions in the third section. Table 2 shows storage for the DH phone according to this method. In Table 2, "node" corresponds to the state number; "arc" corresponds to a transition between two nodes (i.e., from a state to a state); "count" is a value determined during training; and "prob" is a probability value for a corresponding transition. Three counts are provided for each label output AA11 to ZX3-. Table 2 lists label outputs AE13 through XX23 which have count values above a given threshold (e.g. 0.01). Although the table is condensed for purposes of facilitated viewing, it should be recognized that the table actually stored in the system has three count entries for each of the 200 labels. The first column of counts corresponds to counts generated in transitions corresponding to the first section of the phone machine; the second column corresponds to counts for the second section; and the third column corresponds to counts for the third section. It is noted that each non-null transition is assigned to a corresponding column number 4, 5, or 6 --for example, transition 1 → 2 is assigned to column 4 and transition 6 → 6 is assigned to column 6. Another memory structure for a phone machine is illustrated by Fig. 4 and Fig. 5. In Fig. 4, one phone machine is identified by its numerical label 3 (e.g., for the DH phonetic element phone machine). In memory, the number of states of the phone machine is stored. For each state, the number of arcs (or transitions) is also stored. Also, for each arc, the destination state is indicated and an arc type is specified. As in Table 2, it is assumed that each phone machine has three label output probability distributions associated therewith; that is a label output distribution for the first section, a label output distribution for the second section and a label output distribution for the third section of the phone machine. For 80 phone machines, there are 240 label output probability distributions. In Fig. 5, there are 240 label output distributions stored in a matrix form. Each column of the matrix represents a label output probability distribution (across the 200 labels) and each corresponds to a distinct arc type. In Fig. 4, the arc type entry points to one of the columns in the Fig. 5 matrix. With probabilities stored for each phone machine and with prescribed sequences of phone machines representing words in a vocabulary, the acoustic match processor 106 is able to determine the likelihood of a word given a string of labels (or speech types). That is, given a string of labels, the acoustic processor 106 examines the phone machines for a first word and determines the likelihood of those phone machines generating successive labels in the string. Briefly, the acoustic match score for a word is computed by following the various paths the label string could follow through the phone machines which form the word and applying the probabilities appropriately. By way of short example, suppose the acoustic processor 102 (of Fig. 1) generates labels L200-L199-L15-L14-... The computing of an acoustic match score by the acoustic match processor 106 (of Fig. 1) commences as follows for the word "THE". At the first time interval, label L200 is generated by the acoustic processor 102 (of Fig. 1). Two initial branches are considered: along the first branch, label L200 may be produced along transition tr₁; and along the second branch label L200 may be produced along transition tr₈. The acoustic matcher 106 looks to the phone machine for phonetic element DH to see what the stored probabilities are for (1) taking transition tr₁ and (2) producing label L200 at transition tr₁ (in the first section of the phone machine). If the tr₁ probability is 0.5 and the label output probability for L200 is 0.001, then the likelihood of following the first path and producing the L200 output is (0.5)(0.001)=0.0005. After reaching state S4, there are two possible paths for producing label L199. One path involves following the self-loop transition tr₂ and producing label L199. The other path involves following transition tr₃ to state S5 and producing label L199. Suppose that the probability of following transition tr₂ in the DH phone machine is 0.3 and that the probability of following transition tr₃ is 0.7. Also suppose that the probability of producing the label L199 on transition tr₂ in the DH phone machine is 0.025 and that the probability of producing the label L199 on transition tr₃ in the DH phone machine is 0.004. The probability of producing label L199 along transition tr₂ is computed as (0.3)(0.025)=0.0075 . The probability of producing label L199 along transition tr₃ is computed as (0.7)(0.004)=0.0028 . The likelihood of producing labels L200-L199 along transitions tr₁ and then tr₂ is computed as (0.0005)(0.0075) and the likelihood of producing labels L200-L199 along transitions tr₁ and then tr₃ is computed as (0.0005)(0.0028). The joint probability of following producing label L200 at transition tr₁ and producing label L199 at transition tr₂ or tr₃ is then PROB I = (0.0005)(0.0075)+(0.0005)(0.0028) = 0.00000515 Along the second branch extending from state S1 along transition tr₈, the acoustic matcher 106 (of Fig. 1) retrieves the probability of following transition tr₈ and the probability of producing label L200 at transition tr₈. The product of p(tr₈) x p⁸[200] yields the probability of following transition tr₈ and producing label L200 along that transition. If label L199 is to be produced by the DH phone machine from state S2, it must be along transition tr₉. The related probability of that event is p(tr₉) x p⁹[199]. Hence, the probability of generating L200-L199 by following transitions tr₈ and then tr₉ is PROB II = {p(tr₈) x p⁸[200] }{p(tr₉) x p⁹[199]}. The probability of the DH phone machine producing the labels L200-L199 is then PROB = PROB I + PROB II . For the labels which follow L200 and L199, computations are further performed, and expanded along the various possible branches. At this point, it is observed that the DH phone machine of Fig. 3 may be traversed with any number of label outputs being produced. Along null transition tr₁₁, the DH phone machine may be, in effect, passed over without producing any label output. In this event, the probability of following the null transition tr₁₁ enters the computation and the probability of generating the label string continues with the UH1 phone machine. There are various known techniques for partitioning the label string either based on words or based on phone machines. In the former case, words may be uttered with distinct pauses therebetween, thereby enabling the system to detect which substring of labels corresponds to a word based on intervals of detected silence. For phonetic phone machines and labels which are generated by an acoustic processor 102 at centisecond intervals, there are typically eighty to one hundred labels per word. With such isolated word speech, the labels corresponding to one word after another are identified. In the aforementioned US Patent 4,718,094 probability distributions for phonetic phone machine end-times are determined. Accordingly, labels for one phone machine after another, based on probabilities, are determined. One vocabulary word after another is selected in order to compute an acoustic match score for each word. For each word, the sequence of phone machines --including the probabilities stored therefor-- is retrieved. The possible paths which the identified labels can follow through the sequence of phone machines is determined and, as illustrated hereinabove, probabilities are computed. This approach is called a detailed acoustic match and is described in the aforementioned US Patent. B. Speech Recognition with Fenemic Phone Machines As an alternative to phonetic phone machines, "fenemic" phone machines may be employed in constructing the word baseforms. A "feneme" is simply a label. A "fenemic" phone machine is a phone machine which corresponds to a feneme, or label, rather than to a phonetic element. Where there are 200 labels in the alphabet of labels, there are 200 fenemic phone machines --one for each feneme. Each fenemic phone machine has a structure 160 such as that shown in Fig. 6. The fenemic phone machine has two states S1 and S2. Between states S1 and S2 are a null transition and a non-null transition. At state S1 is a self-loop non-null transition. Each transition has a respective probability determined during training. Each non-null transition has a probability for each label output. For example, the self-loop transition for the fenemic phone machine AA11 has a transition probability p(tr s ) and a probability of producing each label at the self-loop transition: The probability of producing the speech type AA11 (e.g.label L1) for the AA11 fenemic phone machine should be high. However, due to differing speech pronunciations from one time to another, other labels normally also have a probability of being produced during a transition of the AA11 phone machine. As with phonetic phones, word baseforms can be constructed as respective sequences of fenemic phone machines. There are normally on the order of eighty to one hundred fenemic phone machines per word where labels are generated every centisecond. A convenient manner of depicting all the paths through successive fenemic or phonetic, phone machines is in trellis form. A fenemic phone machine trellis 180 is shown in Fig. 7. For a given string of labels generated by an acoustic processor 102 (of Fig. 1), a variety of paths through the trellis may be followed. Starting at the upper left --state S1 at time t1-- the first fenemic phone may be skipped by following the null (dashed) transition vertically down to state S2. State S2 represents the initial state of the second fenemic phone machine. Also from state S1 at time t1, a label output may be produced following the self-loop, which is depicted as a horizontal move from state S1 back to state S1 at time t2. Lastly, a label output may be produced with a move from state S1 to state S2 at time t2. For fenemic phones, each point on the trellis provides similar alternative paths. Still referring to Fig. 7, it is noted that the trellis provides numerous possible routes for a string of labels. For example, for two labels produced between states S1 at time t1 and S3 at time t3, thirteen routes are possible. For example, a first route includes two self-loops at state S1 followed by two null transitions at time t3; a second route includes a self-loop transition at state S1, a transition to state S3 at time t3, and a null transition at time t3; a third route includes a self-loop at state S1, a null transition to state S2, a self-loop at state S2, and a null transition to state S3; a fourth route includes a self-loop at state S1, a null transition to state S2, a transition to state S3 at time t3; a fifth route includes a null transition at time t1, a self-loop at state S2, and a transition to state S3 at time t3; a sixth route includes a null transition at time t1, a transition from state S2 to S3 at time t2, and a null transition at state S3 and so on. The computations are preferably performed by iteratively computing the probabilities of each node (i.e. state-time point) in the trellis. This is described in detail in the above-mentioned F. Jelinek article "Continuous Speech Recognition Statistical Methods." Apparatus and methodology for constructing fenemic word baseforms formed of fenemic phones are set forth in the European Patent Application 86 104 220.8 which is incorporated herein by reference to the extent (if any) required for adequate disclosure, and which is entitled "Feneme-based Markov Models for Words". As with the present invention and other patent applications cited herein, this application pertains to an invention made by members of the IBM Corporation Speech Recognition group. Whether phonetic phone machines or fenemic phone machines are used, acoustic match scores for word baseforms are computed for at least some words in the vocabulary. Words having the higher scores are more likely candidates for the spoken word to be recognized. C. Training the Phone Machines Having described the operation of a Hidden Markov Model speech recognizer which is based on phone machines (such as phonetic phone machines and fenemic phone machines) which have been trained, the present invention which relates to training is now disclosed. In particular, apparatus and method for determining the label output probabilities corresponding to a subsequent speaker following the training for a reference speaker is now disclosed. Fig. 8 is a general block diagram of apparatus 200 for training (Markov model) phone machines according to the invention. Input speech of a training text is shown entering a label generator 202. Label generator 202 looks at one time interval of speech after another and determines the values for each of various prescribed characteristics, or features, relating to the speech. Each characteristic corresponds to one component of a feature vector. At each interval, a feature vector is generated. For each speaker, there are specified a distinct alphabet of prototype vectors. The alphabet of prototype vectors for each speaker is stored in prototype storage 204. For each successive feature vector at successive time intervals, the closest prototype vector thereto is determined. The label for the closest prototype vector is then assigned to the given interval. For successive intervals, a string of labels is generated. Referring to Fig. 2, it is observed that the label generator 202 and prototypes storage 204 correspond to the acoustic processor 102. Referring back to Fig. 8, a probability initializer 206 is shown. The purpose of the probability initializer 206 is to provide starting point probability values from which accurate, informative probability values may be determined. The initialized values normally do not closely reflect real speech and are inadequate for use in speech recognition. A preferred initializer 300 is depicted in Fig. 9. In accordance with the invention, it is assumed that a reference speaker has previously uttered a full training text, that a series of input vectors were generated in response to uttered speech of the reference speaker, and that transition probabilities P(A ij ) and label output probabilities P(F₁(k)|A ij ) have been previously determined. The input vectors for the reference speaker enter a labeller element 302. The labeller element 302 assigns a label to each input vector. Significantly, the labeller element 302 assigns labels based on the prototype vectors defined for a subsequent speaker. The prototype vectors for the subsequent speaker are determined from data only generated in response to the subsequent speaker uttering part of the training text. The defining of prototype vectors and clustering feature vectors into speech types is well-known in the art. The string of labels exitting labeller element 302 is referred to as A₁ = x₁x₂x₃.... The label string A₂ enters an initial probability calculator 306. Also entering the initial probability calculator 306 is alignment information from time alignment storage 308. Storage 308 contains data correlating phone machine transitions with labels. In this regard, it is observed that during the training of the reference speaker, a text corresponding to a given sequence of phone machines is uttered. From Viterbi alignment, a correlation (time alignment) between each phone machine (for the reference speaker) and a label segment in a string A₁, (i.e. a label string based on reference speaker prototypes) is known. The string A₂ is aligned label-by-label the same as is string A₁. If b1 labels of string A₁ are aligned with a first phone machine, b1 labels of string A₂ are also aligned with the first phone machine; and so on. The correlation between labels of A₂ and phone machines are therefore provided in storage 308. By way of example, reference is made to a phone machine PM1. Viterbi alignment during the training period for the reference speaker specifies that the first p labels of string A₁ correlate to phone machine PM1. The number of labels assigned to each successive phone machine transition are stored in the time alignment storage 308. The initial probability calculator 306 divides the string of labels A₂ in the same manner as is string A₁ so that the first p labels of string A₂ correlate to the transition PM1. Based on the alignment between phone machines and label string A₂, initial counts for the subsequent speaker -- which correspond to occurrences of a given label for a given phone machine -- are computed. For example, in the training text, the phone machine PM1 may occur fifteen times. The labels which correlate to each occurrence of PM1 are known due to the alignment process. The number of times each label (i.e., speech type), is aligned so that it occurs during an occurrence of phone machine PM1 is counted. That is, for label L1 aligned with a phone machine PM1, there is a corresponding count; for label L2 aligned with phone machine PM1, there is a corresponding count; and so on. A similar procedure is followed for each phone machine, resulting in each phone machine having a respective count for each label. The initial probability calculator 306 normalizes the counts to produce probabilities of the form P₀(F₂(ℓ)|PMi) where ℓ represents a label and PMi represents the ith phone machine. Assuming an even distribution of the labels over the various transitions of a phone machine, the initial probability calculator 306 provides as output label output probabilities of the form P₀(F₂(ℓ)|A ij ) where P₀ indicates an initial probability. The initial transition probabilities P₀(A ij ) are simply obtained by assigning fixed values which are speaker independent. Referring again to Fig. 8, the initial probabilities are shown entering a re-parameterized forward-backward processor 208 and a basic parameter forward-backward processor 210. The reparameterized forward-backward processor 208 receives as input (i) the initialized probabilities, (ii) a string of labels y₁y₂y₃..., generated -- based on subsequent speaker prototypes -- when the subsequent speaker utters part of the training text, and (iii) the label output probabilities for the reference speaker P(F₁(k)|A ij ) as previously determined during the reference speaker training period. In brief, the re-parameterized forward-backward processor 208 re-specifies an original Markov model to an expanded graph form which includes similarity data between the reference speaker and the subsequent speaker. Stated otherwise, each non-null transition in an original Markov model has a transition probability and probabilities of producing respective labels at the transition; for the expanded graph, output probabilities of the Markov model are redefined. Each output probability for the Markov model P(F₂(ℓ)|A ij ) is re-specified in terms of the sum of similarity measures P(F₂(ℓ)|F₁(k)) over all k values normalized by P(F₁(k)|A ij ) where k identifies a label in the alphabet of N labels. The re-parameterized forward-backward processor 208 produces re-parameterized label output probabilities based on expanded Markov models according to the expression: Fig. 10 illustrates an original Markov graph for a label output probability. Fig. 11 shows an expanded graph based on the expression set forth immediately hereinabove. A specific embodiment of a re-parameterized forward-backward processor 400 is illustrated in Fig. 12. The initialized probabilities P₀(A ij ) and P₀(F₂(ℓ)|A ij ), and the string of labels for the subsequent speaker y₁y₂y₃... are shown entering a forward-backward algorithm processor 402. In response to these inputs, the processor 402 generates a transition count count n+1 (A ij ) and a label output count count n+1 (F₂(ℓ), A ij ). "n" is an iteration index value which is set to zero for the initialized probabilities. The transition count is normalized in a transition probability calculator 404 to provide an updated transition probability value P n+1 (A ij ) which enters a counter multiplexer 406. For values of n < nmax (where nmax is preferably 3), the updated transition probability value enters the processor 402 as input thereto. When n reaches nmax, the updated transition probability is outputted from the re-parameterized forward-backward processor 400 (as shown in Fig. 12) as a re-parameterized transition probability P(A ij ). Referring to Fig. 12, the label output probability count count n+1 (F₂(ℓ), A ij ) is shown entering an expanded graph count calculator 408. Also entering calculator 408 are the label output probabilities previously determined for the reference speaker P(F₁(k)|A ij ); a previously determined value for the label output probability at time n, P n (F₂(ℓ)|A ij ); and a previously determined similarity measure value P n (F₂(ℓ)|F₁(k)). The label output probabilities P(F₁(k)|A ij ) are received as input from an external element which stores data relating to the reference speaker. For the first iteration, n=0. Each P₀ probability thus represents an initialized value. The calculator 408, in effect, performs the following calculation to determine relational counts: The right side of the expression includes elements which are provided as inputs to the calculator 408. By summing the counts count n+1 (F₂(ℓ), F₁(k), A ij ) over all transitions A ij in an adder 410, the following computation is performed: By normalizing the summed counts in an updated confusion matrix entry calculator 412, updated similarity measure values P n+1 (F₂(ℓ)|F₁(k) at time (n+1) are computed. The most recent updated value for each similarity measure value is entered into a confusion matrix storage 414. The storage 414 contains one entry for each (ℓ,k) pair of labels. Hence, for an alphabet of N=200 labels, the matrix has 200x200 entries. A sample matrix is shown in Fig. 13. In the sample matrix, there is a 0.002 probability of label L1 for the subsequent speaker given label K1 of the reference speaker. Similarly, there is a 0.010 probability of label L2 for the subsequent speaker given label K1 of the reference speaker. The probabilities for L3, L4, and L200 (given label K1 of the reference speaker) are 0.001, 0.004, and 0.001, respectively. Labels L5 through L199 also have respective probabilities for label K1. The total of the probabilities along the K1 column (and each other column) equals one. That is, for a given reference speaker label k, one of the labels ℓ for the subsequent speaker is produced. Referring again to Fig. 12, it is observed that the confusion matrix entries are iteratively updated. The most recent entry values on each iteration are provided as input to the calculator 408 for a next iteration (if any). In addition, however, the most recent updated entry values enter an updated label output probability calculator 416. The other input to calculator 416 are the label output probabilities associated with the reference speaker; namely, P(F₁(k)|A ij ) for each k. The calculator 416 produces an updated label output probability After n=nmax iterations, an update counter multiplexer 418 provides the updated label output probability P r (F₂(ℓ)|A ij ) as an output from the re-parameterized forward-backward processor 400. For previous iterations (n < nmax), the label output probability from the calculator 416 is sent as P n (F₂(ℓ)|A ij ) to processor 402 and calculator 408 for a next iteration. By way of further comment, it is noted that the forward-backward algorithm processor 402 operates in accordance with the teachings set forth in the aforementioned IEEE articles. In brief, the forward-backward algorithm involves computations relating to movement from a state i in a Markov model to a state j in the Markov model. In arriving at state i at a time t, starting at an initial state, there is a predetermined probability α t (state(i)). Similarly, there is a probability of traversing from state j at time (t+1) through to the final state of the Markov model β t+1 (state(j)). The α probability involves moving forward to state i whereas β involves computing backward from the final state. The count indicating how often a path from the initial state to the final state via the i → j transition is followed wherein label ℓ is produced is represented by the expression: where t represents time intervals from t=1 to t=T and δ = 1 if y t = F₂(ℓ) and δ = 0 otherwise. Of particular significance, it is noted that all counts count n+1 (F₂(ℓ), A ij ) for the expanded Markov model graph are readily calculated by using counts obtained for the original graph --i.e., counts related to label output probabilities. The processor 402 performs the above computation on successive iterations (as n increments in value). Referring to Fig. 14, an embodiment of a basic parameter forward-backward processor 500 is shown. The basic parameter forward-backward algorithm processor 502 --like processor 402-- receives the initialized probabilities and label string y₁y₂y₃ ... as inputs. Probability values P′ n (A ij ) and P′ n (F₂(ℓ)|A ij ) --computed as described below -- enter as iteratively updated inputs to the processor 502 when they are generated. In the first iteration n=0. The processor 502 determines basic counts count′ n+1 (A ij ) and count′ n+1 (F₂, A ij ) for transitions and label outputs respectively. Calculators 504 and 506 transform the counts into basic transition probabilities and basic label output probabilities, respectively, by normalizing the counts. The basic probabilities P′ n+1 (A ij ) and P′ n+1 (F₂(ℓ)|A ij ) enter an updated counter multiplexer 508. If the count for n is less than a value nmax′, the basic probabilities enter the processor 502 as inputs P′ n (A ij ) and P′ n (F₂(ℓ)|A ij ) for the next iteration of forward-backward processing. When n reaches nmax′, the last update basic probabilities P′(A ij ) and P′(F₂(ℓ)|A ij ) are outputted. Referring again to Fig. 8, it is noted that the probabilities P r (A ij ) and P r (F₂(ℓ)|A ij ) resulting after re-parameterization by the re-parameterized forward-backward processor 208 enter a deleted estimation processor 212. Similarly, the basic parameter probabilities P′(A ij ) and P′(F₂(ℓ)|A ij ) which are output from the processor 210 enter another deleted estimation processor 214. Each deletion estimation processor 212 and 214 operate in a similar manner to linearly average probabilities and to smooth the distribution of probabilities. By smoothing, is meant, rounding peaks and raising valleys in the label output probability distribution. The smoothing is required because, due to sparse training data, sharp contrasts in label output probabilities may occur which, with further training data, would be less extreme. The concept of deleted estimation processing is described in a 1979 work by F. Jelinek and R.L. Mercer entitled "Interpolated Estimation of Markov Source Parameters from Sparse Data " In addition, a work entitled Pattern Recognition in Practice editted by Gelsema et al. (1980) at pages 381-402 discusses deleted estimation processing. Both references are incorporated herein by reference to provide background. In Fig. 15, a preferred deleted estimation processor 600 is shown. The deleted estimation processor 600 has three inputs. The first input is the string of labels y₁y₂y₃.... The second input corresponds to forward-backward (F-B) processed probabilities. For deleted estimation processor 212, the second input corresponds to the output from the re-parameterized forward-backward processor 208. For deleted estimation processor 214, the second input corresponds to the output from the basic forward-backward processor 210. The third input corresponds to contrast probabilities. For deleted estimation processor 212, the third input corresponds to the initialized probabilities produced by the initializer 206. For deleted estimation processor 214, the third input corresponds to the smoothed probabilities computed by the deleted estimation processor 212. The output from a deleted estimation processor includes (i) linear averaged transition probabilities and (ii) linear averaged label output probabilities, where the probability distributions are smoothed. The preferable deleted estimation processor 600 (Fig. 15) includes a training data divider 602 which separates the training text labels of the subsequent speaker into two portions. The first portion of labels and the forward-backward processed probabilities (from processor 208 or 210) enter a first forward-backward processor 604. Processor 604 computes counts and probabilities therefrom. The second portion of labels and the forward-backward processed probabilities (from processor 208 or 210) enter a second identical forward-backward processor 606. Processor 606 computes counts and probabilities therefrom. The respective counts and probabilities from the processors 604 and 606 enter a maximum likelihood smoother processor 608. The other input to processor 608 are contrast probabilities. The processor 608 performs a linear averaging from the inputs thereto. In the case of deleted estimation processor 212, the label output counts and the transition counts from F-B processor 604 are identified as C₁(F₂(ℓ), A ij ) and C₁(A ij ) and the label output probabilities are identified as P₁(F₂(ℓ)|A ij ). The label output counts from F-B processor 606 are identified as C₂ and the label output probabilities are identified as P₂. Preferably, C₁ and P₁ correspond to the odd sentences label string and C₂ and P₂ correspond to the even sentences string, based on the separation performed by the training data divider 602. The computation performed by the maximum likelihood and smoothing processor 608 of the deleted estimation processor 212 is: wherein λ m is a selectable weighting factor and wherein A ij represents a transition from a state i to a state j in a phone machine. Preferably the factor λ m is selected to maximize the following expression with respect to λ m wherein π represents the product over all m values or all A ij εS m values: where P₀ represents an initial probability and wherein S m represents an mth class of M classes of transitions. All (A ij ) events are preferably classified so that S m = {(A ij )|T m-1 < C₁(A ij ) + C₂(A ij ) < T m } where T m is a predefined count range limit separating one arc class from another. In explaining the value M, it is noted that each of the 240 arcs has a count obtained during subsequent speaker training. Prior to smoothing, arcs (i.e., transitions) are classified into M different classes based on count value. For example, arcs having count range limits from 0 to 10 form a first class (m=1), arcs having limits from 10+ to 100 form a second class, and arcs having a limit of 100+ form a third class. (In such case M=3.) The amount of smoothing depends on which class an arc belongs to. M is typically less than 10. In the case of deleted estimation processor 214, the label output counts and transition counts from F-B processor 604 are identified as C*₁ and the label output probabilities are identified as P*₁. The label output counts and transition counts from F-B processor 606 are identified as C*₂s and the label output probabilities are identified as P*₂. Preferably, C*₁ and P*₁ correspond to the odd sentences label string and C*₂ and P*₂ correspond to the even sentences label string, based on the separation performed by the training data divider 602. The computation performed by the maximum likelihood smoother processor 608 of the deleted estimation processor 214 is: Preferably the factor λ* m is selected to maximize the following expression with respect to λ* m . where P s represents a smoothed probability and where S* m is an mth class of M classes of transitions based on transition count values and is defined as S* m ={A ij |T m-1 < C*₁(A ij ) + C*₂(A ij ) < T m }. The label output probabilities computed by the deleted estimation processor 214, P*(F₂(ℓ)|A ij ), are greatly improved values over the basic label output probabilities P′(F₂(ℓ)|A ij ); the re-parameterized label output probabilities P r (F₂(ℓ)|A ij ); and the smoothed label output probabilities output from the first deleted estimation processor 212, namely P s (F₂(ℓ)|A ij ). A deleted estimation process similar to the one described above is also applied to the transition probabilities to obtain P*(A ij ). However, it is noted that the degree of improvement over P′(A ij ), P r (A ij ), or P s (A ij ) is not as significant as for the case of output probabilities. In this regard, the reparameterization methodology of determining label output probabilities from re-parameterized values has been implemented on the IBM MVS system. A series of experiments has been performed for many speakers using a sample training script. When a total of 100 sentences of sample text is used to train subsequent speakers, an average recognition rate of 95% is achieved. Using a shortened training text of 30 sentences --rather than 100 sentences-- results in 89% accuracy. Employing the re-parameterization with the shortened training text of 30 sentences has resulted in an average recognition rate of 92.5%. A 30% improvement in the error rate as compared to that of straight 30 sentences training is achieved by the re-parameterization alone. When the two deleted estimation steps are added to the process, accuracy is raised further to approximately 95%. With 1/3 the training text for a subsequent speaker, the present invention permits average recognition accuracy similar to the accuracy achieved based on full text training. Hence, instead of a 20 minute period for uttering training text, a subsequent speaker need utter only 7 minutes of training text. Moreover, notwithstanding the various computations discussed hereinabove, the present invention achieves computational savings as compared to processing full training text by conventional forward-backward algorithm techniques. Fig. 16 shows a flowchart of an algorithm 700 which implements the present invention. At steps 701 and 702, index value n is set to zero and initialized values for transition, label output, and confusion matrix entry probabilities are set. In step 704, a forward-backward algorithm is applied to generate counts for the original Markov model graph in terms of transition probabilities and label output probabilities. In step 706, the counts count n+1 (F₂(ℓ), A ij ) for the label output probabilities are respecified to provide counts for an expanded Markov model graph (in which the output probabilities correspond to confusion matrix entries) as count n+1 (F₂(ℓ), F₁(k), A ij ). The counts are summed over all i → j transitions to provide counts (step 708) that reflect a similarity measure between an ℓth label for the subsequent speaker and a kth label for the reference speaker. The entry counts are used to update the probability values for the confusion matrix entries (step 710). Updated re-parameterized label output probabilities are computed in step 712. While the updated re-parameterized label output probabilities are being determined, conventional techniques are employed in updating the transition probabilities in step 714. After steps 712 and 714, the value of n is incremented (step 716). If the new value for n is greater than or equal to a prescribed value nmax (where nmax is preferably 3) at step 718, the most recent output at steps 712 and 714 enter a first deleted estimation step 720. If n is less than nmax, steps 704 through 718 repeat. The deleted estimation step 720 provides as outputs smoothed probabilities P s (A ij ) and P s (F₂(ℓ)|A ij ) wherein the label output probabilities are based in part on similarity information between the reference speaker and the subsequent speaker. While the probabilities P s (A ij ) and P s (F₂(ℓ)|A ij ) are being computed, a basic forward-backward algorithm is iteratively applied, starting with the initial probabilities for the first iteration at n=0 at step 730. Over successive iterations, updated transition and label output probabilities are computed in standard fashion in step 732. Steps 734 and 736 effect the iterations until n ≧ nmax is reached, whereupon the last update basic transition probabilities P′(A ij ) and last update basic label output probabilities P′(F₂(ℓ)|A ij ) are provided as output directed to a second deleted estimation step 740. In step 740, the smoothed probabilities from the first deleted estimation step 720 and the last update basic label output and transition probabilities following step 736 enter a deleted estimation processor (such as processor 600 described hereinabove) together with the label string y₁y₂y₃.... At step 740, final label output probabilities P*(F₂(ℓ)|A ij ) and transition probabilities P*(A ij ) are determined. In this regard it is observed that the present invention is applicable in general to Markov model speech recognition. Accordingly, the invention applies whether phonetic phone machines, fenemic phone machines, or other Markov model phone machines are employed.

Claims: A speech recognizer comprising acoustic processor means for outputting in response to a speech input a string of labels one after another at successive time intervals out of an alphabet of labels, a differently designed alphabet being provided for each speaker, and in which words correspond to respective sequences of Markov model phone machines, wherein each phone machine has a plurality of transitions each extending from a state i to a state j and wherein the probability of each transition and the probability of label outputs being produced at this transition differ among speakers, comprising means for training the recognizer to a subsequent speaker after the recognizer is trained with a reference speaker in response to the utterance by said reference speaker of a sampled text, characterized in that the training means further comprises means for computing from labels generated by said acoustic processor means in response to the subsequent speaker uttering part of the sampled text a first set of label output probability values of a set of Markov model phone machines corresponding to the said part of the sampled text; means for computing and storing each entry of a matrix representing the probability of a label l spoken by the subsequent speaker corresponding to a label k spoken by the reference speaker; means for generating a second set of label output probabilities for the said set of Markov model phone machines which are based on the values of the said matrix and the label output probabilities of the reference speaker; and means for linear averaging the first set of label output probabilities against the second set of label output probabilities, given the labels generated in response to the subsequent speaker uttering said parts of the sample text, to produce final label output probabilities for the said set of Markov model phone machines. The speech recognizer of claim 1 characterized in that said linear averaging means includes a deleted estimation processor. The speech recognizer of claim 2 characterized in that said deleted estimation processor includes: means for dividing label inputs generated by the acoustic processor means during utterance of sample text into two portions; means for applying a forward-backward algorithm to the first set of label output probabilities for one of the two portions of labels to produce counts C * ₁ and probabilities P * ₁; means for applying a forward-backward algorithm to the first set of label output probabilities for the other portion of labels to produce counts C * ₂ and probabilities P*₂; smoothing processor means, which receives as input the counts C * ₁ and C * ₂ the probabilities P * ₁ and P * ₂ together with the second set of label output probabilities, for producing final label output probabilities wherein λ* m is selected to maximize the likelihood of the final label output probabilities corresponding to spoken text and wherein A ij represents a transition from a state i to a state j in a phone machine. The speech recognizer of claim 3 characterized in that the smoothing processor means includes means for evaluating λ* m according to the expression: where S * m represents an mth class of transitions such that S* m = {A ij |T m-1 <C*₁(A ij )+C*₂(A ij )<T m }; wherein T m represents a predefined count range limit and wherein each class of transitions is defined between a pair of adjacent range limits. The speech recognizer of claim 1 characterized in that said means for generating the second set of label output probabilities includes: re-parameterizing means for computing smoothed label output probabilities P s (F₂(ℓ)|A ij ) in response to receiving as inputs (i) the known label output probabilities of the reference speaker P(F₁(k)|A ij ), (ii) initialized label output probabilities and transition probabilities P₀(F₂(ℓ)|A ij ) and P₀(A ij ), respectively, and (iii) labels generated by said acoustic processor means in response to the subsequent speaker uttering said part of the sample text. The speech recognizer of claim 5 characterized in that said re-parameterizing means includes: means for computing a respective probability value P [F₂(ℓ)|F₁(k)] for each entry in a confusion matrix indicative of the probability of the ℓth label as spoken by the subsequent speaker corresponding to the kth label as spoken by the reference speaker; means for storing the values for the confusion matrix entries; means for expanding each label output probability of the subsequent speaker P [F₂(ℓ)|A ij ] in terms of confusion matrix entry values and label output probabilities for the reference speaker as: where N is the number of labels in the alphabet of labels. The speech recognizer of claim 6 characterized in that said expanding means includes: means for calculating a transition probability P(A ij ) in response to a transition count input; forward-backward processor means for iteratively producing transition counts count n+1 (A ij ) and label output counts count n+1 (F₂(ℓ), A ij ) in response to a transition probability input P n (A ij ) a previously calculated label output probability P n (F₂(ℓ)|A ij ), and labels generated in response to the subsequent speaker uttering said part of the sample text; entry update calculator means for computing a value P n+1 (F₂(ℓ)|F₁(k)) in response to receiving as input (i) a previous entry value P n (F₂(ℓ)|F₁(k), (ii) label output probabilities for the reference speaker P(F₁(k)|A ij ), (iii) the most recent label output count count n+1 (F₂(ℓ), A ij ) from the forward-backward processor means, and (iv) a previous label output probability value P n (F₂(ℓ)|A ij ) and calculator means for computing updated re-parameterized label output probabilities The speech recognizer of claim 7 further characterized by: initializing means for producing initial values for P₀(A ij ) the transition probabilities and for P₀(F₂(ℓ)|A ij ) the label output probabilities; and deleted estimation processor means for linear averaging the re-parameterized label output probabilities received as input from said re-parameterizing means against the initial label output probabilities received as input from the initializing means, given the label outputs generated by the acoustic processor means in response to the subsequent speaker uttering said part of the sample text, to produce smoothed label output probabilities P s (F₂(ℓ)|A ij ) which have a smoothed probability distribution. A machine-implemented method for training a Markov model speech recognition system to a subsequent speaker after the system is trained to a first speaker, wherein each word corresponds to a sequence of Markov model phone machines and wherein speech is labelled by an acoustic processor based on a speaker-dependent label alphabet, comprising the steps of generating a string of labels in response to the utterance of a sample text by a first speaker and computing, based on the label string generated for the first speaker the probabilities P(F₁(k)|A ij ) of each kth label F₁(k) being produced at a transition A ij from a state i to a state j of a Markov model phone machine for the first speaker, characterized by the steps of (a) computing and storing from labels generated by said acoustic processor in response to the subsequent speaker uttering sampled text each entry of a matrix representing the probability of the lth label F₂(l) spoken by the subsequent speaker corresponding to the kth label F₁(k) spoken by the reference speaker, each entry being of the form P(F₂(ℓ)|F₁(k)) (b) computing probabilities of each lth label F₂(l) being outputted at a transition of A ij of a Markov model phone machine for the subsequent speaker as: where N is the number of labels in the alphabet. The method of claim 9, characterized by further steps of (c) generating a short string of labels in the acoustic processor in response to the utterance by the subsequent speaker of a shortened sample text; (d) initializing transition probability values P n (A ij ) and label probability values P n (F₂(ℓ)|A ij ) for the phone machines corresponding to the subsequent speaker, based on the short string of labels, where index value n = 0; (e) computing, with a forward-backward algorithm, up-dated transition counts count n+1 (A ij ); (f) computing, with a forward-backward algorithm, up-dated label output probability counts as: where α t is a cumulative probability in moving forward to state i from an initial state and where ß t+1 is a cumulative probability of proceeding backward from a final state of a Markov model sequence of states to state j; (g) from counts count n+1 (A ij ) and counts count n+1 (F₂(ℓ), A ij ), computing up-dated transition probability values P n+1 (A ij ) and up-dated values for P n+1 (F₂(ℓ)|A ij ). The method of claim 10, characterized in that computing up-dated P n+1 (F₂(ℓ)|A ij ) values comprises the steps of (h) computing a relational count as: (i) summing the relational counts over all transitions A ij as: (j) from the summed relational counts, estimating up-dated confusion matrix values as: The method of claim 11, characterized by the further steps of (k) incrementing the value of n and repeating steps (e) through (j) over a prescribed number of iterations; the label output probabilities resulting after the prescribed number of iterations representing updated reparameterized label output probabilities. The method of claim 12, characterized by the further step of (l) producing, in response to the entry of (i) updated transition probabilities and reparameterized label output probabilities, (ii) initialized probabilities, and (iii) the short string of labels into a deleted estimation processor linear averaged smoothed enhanced probabilities between the updated reparameterized label output probabilities and the initialized label output probabilities. The method of claim 13, characterized in that the step of producing smoothed probabilities includes the steps of (m) dividing the short string of labels into two portions; (n) applying a forward-backward algorithm, which is based on the updated probabilities, to the first portion of labels to generate counts C₁(F₂(ℓ), A ij ) and C₁(A ij ) and a probability P₁(F₂(ℓ)|A ij ) for each label output; (o) applying a forward-backward algorithm, which is based on the last update probabilities, to the second portion to generate counts C₂(F₂(ℓ), A ij ) and C₂(A ij ) and a probability P₂(F₂(ℓ)|, A ij ) for each label output; and (p) evaluating a factor λ m to maximize the likelihood the expression: where P₀ represents an initial probability and wherein S m represents an mth class of M transitions where transitions are classified based on transition count value. The method of claim 14, wherein the step of producing the smoothed probabilities further comprises the step of: (q) computing: wherein λ m is a selectable weighting factor and wherein A ij represents a transition from a state i to a state j in a phone machine. The method of claim 15, characterized by the further steps of (r) computing basic transition probabilities and label output probabilities for the phone machines of the subsequent speaker which includes the steps of applying a forward-backward algorithm to the short string of labels to compute basic counts and normalizing the basic counts; (s) producing, in response to the entry of (i) smoothed transition probabilities and smoothed label output probabilities (ii) basic transition probabilities and basic label output probabilities, and (iii) the short string of labels into a deleted estimation processor, final label output probabilities which are linearly averaged between the smoothed label output probabilities and the basic label output probabilities. The method of claim 16, characterized in that the step of producing final label output probabilities includes the steps of (t) dividing the short string of labels into two portions; (u) applying a forward-backward algorithm, which is based on the smoothed label output probabilities, to the first portion of labels to generate counts C * ₁(F(ℓ)|A ij ) and C * ₁(A ij ) and a probability P * ₁(F(ℓ)|A ij ) for each label output; (v) applying a forward-backward algorithm, which is based on the smoothed label output probabilities, to the second portion to generate counts C*₂(F(ℓ)|A ij ) and C*₂(A ij ) and a probability P * ₂(F(ℓ)|A ij ) for each label output, and evaluating a factor λ* m to maximize the likelihood of the expression: where P s represents a smoothed probability and where S* m is an mth class of M classes of transitions based on transition count values and is defined as S* m = {A ij |T m-1 <C*₁(A ij )+C*₂(A ij )<T m } where T m-1 and T m are adjacent predefined count range limits which specify a respective class of transitions having counts between the two limits. The method of claim 17, characterized in that the step of producing the final label output (w) computing: probabilities further comprises the step of (w) computing: The method of claim 18 further characterized by the step of producing final transition probabilities which are linearly averaged between the smoothed transition probabilities and the basic transition probabilities.

Inventors: BAHL, LALIT RAI, MERCER, ROBERT LEROY, NAHAMOO, DAVID
Applicants: INTERNATIONAL BUSINESS MACHINES CORPORATION
Application year: 1988
Main classification: G10L5
IPC codes: G10L15, G10L11