Computer program product for automatic recognition of a consistent message using multiple complimentary sources of information

Abstract

A general approach is provided for the combined use of several sources of information in the automatic recognition of a consistent message. For each message unit (e.g., word) the total likelihood score is assumed to be the weighted sum of the likelihood scores resulting from the separate evaluation of each information source. Emphasis is placed on the estimation of weighing factors used in forming this total likelihood. This method can be applied, for example, to the decoding of a consistent message using both handwriting and speech recognition. The present invention includes three procedures which provide the optimal weighing coefficients.

Description

TECHNICAL FIELD

The present invention relates generally to the field of message recognition, and more particularly, to a system and method for using several sources of information in the automatic recognition of a consistent message.

BACKGROUND ART

The simultaneous use of multiple sources of information is known to improve the performance of identity recognition systems. For example, it was shown in UK Patent GB 2,229,305 A to Fleming that verification procedures used to identify an individual can be strengthened by the use of multiple biometric inputs. The Fleming disclosure describes an identity verification system based on collecting, through appropriate biometric tests, various user specific characteristics (e.g., face, fingerprint, and voice characteristics), and combining the results of these test through a weighted sum of each individual score.

The recognition of a person, however, is quite different from the recognition of a message conveyed by that person. Even though voice and/or handwriting inputs, for example, may be used in both cases, user identification requires only a relatively short voice pattern (such as a password), or a relatively concise handwriting pattern (such as a signature), for the sole purpose of ascertaining the identity of the user. A user identification procedure typically must compare the evidence provided by the user to a unique set of templates stored for this user. A score is assigned to the evidence and a comparison of the score to a threshold is performed. The comparison is either successful or unsuccessful (i.e., this is a binary decision).

In contrast, message recognition is a far more complex task. Message recognition generally requires long utterances or a substantial amount of handwriting, involving complex syntactic and semantic structures constructed from a large vocabulary of words. The output of a message recognizer is a sequence of words representing the message being created. A message recognition procedure typically consists of decoding the current set of input data (such as, for example, the utterance spoken for the current word or the sequence of characters handwritten for the current word) by comparing the input against a plurality of potential candidate words in the vocabulary. Each candidate word must be assigned a score, and the decoded word is determined by which score is the highest.

Of particular interest are procedures for adjusting internal parameters used by the recognition system to the characteristics of each individual user. As described in F. Jelinek, "The Development of an Experimental Discrete Dictation Recognizer", Proc. IEEE, Vol. 73, No. 11, pp. 1616-1624 (November 1985), this adjustment, called training, has the potential to improve the accuracy performance during recognition (also referred to as decoding).

Moreover, the simultaneous use of multiple sources of information is known to improve the performance of automatic message recognition systems. For example, it was demonstrated in U.S. Pat. No. 3,192,321 to Nassimbene and E. D. Petajan, "Automatic Lipreading to Enhance Speech Recognition", IEEE Trans. Acoust., Speech, Signal Processing, Vol. ASSP-33, No 1, pp. 40-47, (January 1985), that using lipreading information can enhance automatic speech recognition. However, the procedures used in these disclosures provide a trivial straightforward strategy and are based on a simple comparison of output words from an automatic speech recognizer and automatic lipreading recognizer.

Finally, it was shown in U.S. patent application Ser. No. 07/676,601 filed Mar. 28, 1992 now abandoned, naming J. Bellegarda and D. Kanevsky as co-inventors, and entitled "Improved Message Recognition through the Integration of Speech and Handwriting Information," which is hereby incorporated by reference in its entirety herein, that the integration of speech and handwriting information can lead to the improved recognition of a consistent message. This application addresses the problem of merging the evidence obtained from two sources during decoding, but not during training.

As sensor technology develops, it is likely that in the near future, more and more measurements from many different sources of information will be available to the automatic message recognition system. The challenge is therefore to integrate the evidence provided by these many sources to perform more accurate and robust decoding.

Thus, given the fundamental difference between identification of an individual and identification of a consistent message, which makes it clear that the same procedure cannot be used in both cases, and the demonstrated ability of multiple sources of information to improve the performance of a recognition system, it is readily apparent that a procedure for using multiple sources of information in the automatic recognition of a consistent message is needed.

DISCLOSURE OF THE INVENTION

The present invention provides a general approach for the combined use of several sources of information in the automatic recognition of a consistent message. For each message unit (e.g., word) the total likelihood score is assumed to be the weighted sum of the likelihood scores resulting from the separate evaluation of each information source. Emphasis is placed on the estimation of weighting factors used in forming this total likelihood. This method can be applied, for example, to the decoding of a consistent message using both handwriting and speech recognition. The present invention comprises three procedures which provide the optimal weighting coefficients.

BRIEF DESCRIPTION OF DRAWINGS

The foregoing and other features and advantages of the invention will be apparent from the following more particular description of preferred embodiments of the invention, as illustrated in the accompanying drawings, in which:

FIG. 1 is a general block diagram of a message recognition system.

FIG. 2 is a more detailed block diagram of the message recognition system shown in FIG. 1;

FIG. 3a illustrates a handwritten character being processed to generate a feature vector for a point (P);

FIG. 3b illustrates a four-dimensional handwriting feature vector generated for the point (P) by collecting a plurality of spatial attributes;

FIG. 3c illustrates, in the frequency domain, a speech sample being processed to obtain a feature vector therefrom, specifically, the determination of amplitudes associated with a plurality of frequency bands, not necessarily all of equal bandwidth;

FIGS. 4a-c illustrates an example of a warping of two distinct messages S.sub.1 and S.sub.2 ;

FIG. 5 is a flowchart of the procedure for generating weighted coefficient .alpha.;

FIG. 6 is a flowchart of the deleted interpolation procedure for generating .alpha.; and

FIG. 7 is a flowchart of the expectation-maximization procedure for generating .alpha..

BEST MODE FOR CARRYING OUT THE INVENTION

FIG. 1 illustrates a message recognition system 100, constructed and operated in accordance with a preferred embodiment of the present invention, for the automatic recognition of a consistent message. The present invention provides a system and method for combining multiple sources of information (e.g., voice print, written word, facial movement from one or more users, gestures, typing) in order to recognize a consistent message. For each sub-message unit, such as a word, the present invention forms the total likelihood score as the weighted sum of several likelihood scores, one for each distinct source of information. Recognition is based on the maximization of this total likelihood (i.e., a maximum a posteriori decoding scheme is implemented). Specifically, procedures are provided for implementing a training strategy which yields the optimal weighting coefficients for a given set of circumstances.

As employed herein, a message is considered to be information generated by a human subject. The information is entered into the system 100 through, for example, an acoustic signal, a handwritten stroke signal, a signal from a lipreading device, a signal from a steno device, or through a combination of such signals. The information may represent individual characters, including numbers, letters of the alphabet, punctuation and the like. The information may also represent words, phrases, and paragraphs.

As employed herein, the expression "consistent message" means the same message treated by several different recognition procedures, either simultaneously or otherwise. In other words, the meaning of the material gathered at each source of information has to be identical. Consequently, since the present invention only requires consistency in the message being created, the present invention is able to consider a situation where, for example, the spoken data would be provided by one user (e.g., an executive dictating a memo) and the handwritten data would be provided by another (e.g., a secretary correcting errors just made). Moreover, the present invention can accept as an input a consistent message spoken by two different speakers.

Referring in detail to FIG. 1, blocks 110a, 110b, 110c indicate different source models of information. They consist of information of each source and the training device that receives and operates the corresponding source of information. As such, blocks 110a, 110b, 110c indicate that multiple sources of information can be used as an input to the recognition system 100. The information can come from a speech source model 110a, a handwriting source model 110b, or any other source model (labeled *) 110c that can provide a consistent message to system 100 (e.g., facial movements, touch patterns, brain impulses).

In accordance with the invention there is provided an integrated use of all sources of information. By integrated use it is intended that (a) all sources of information are used simultaneously, and a decoding decision is based upon an overall likelihood score, either joint or merged, that an assemblage of characters, phonemes, movements, patterns, etc., represents a specific word; or that (b) the sources are used sequentially, but with each being used to correct any errors generated by the others. As noted above, the sources of information need not be employed by the same user, so long as both provide a consistent message. In some applications, it may be desirable to have one source provided by one user while the others are provided by another.

The following is a description of the simultaneous and sequential operation of the present invention with reference to an Automatic Handwriting Recognizer (AHR) and an Automatic Speech Recognizer (ASR) as the source models 110a and 110b, respectively. Note that only two source models are being used for brevity.

In regard to the simultaneous use of acoustic and handwriting information, it is noted that both contribute to an overall decoding decision. Ambiguous characters in the handwritten information are recognized primarily on the basis of the corresponding spoken utterance and, conversely, the acoustically confusable words are recognized primarily on the basis of the handwriting information. For this technique words are uttered at the same time as they are written and, as a result, natural pauses appear in the spoken utterance. Hence, it is most advantageous to integrate isolated word speech recognition with mixed cursive and discrete handwriting recognition. If, however, the environment is too noisy for this technique it may be more appropriate to replace the traditional ASR by an automatic lip reading recognizer.

In regard to the sequential use of AHR and ASR, it is noted that AHR may be employed as an enhancement to ASR. For this approach, ASR is considered to be the primary recognition means and AHR is employed primarily as a natural and convenient interface to correct ASR errors that are displayed either on a transparent electronic tablet or on a display terminal. This technique is also suitable for entering new words, such as proper nouns, which are not in the ASR vocabulary but present in an AHR dictionary.

Alternatively, ASR may be employed as an enhancement to AHR. For this approach, AHR is considered the primary recognition means and ASR is used for the purpose of correcting AHR errors. By pointing on AHR errors in some pre-defined manner the speech transducer is activated so that the user may pronounce a correct version of the wrongly decoded handwritten word.

Referring again to FIG. 1, block 130 denotes parameters (.theta..sub.A, .theta..sub.H, . . . .theta..sub.*) that were produced from each of the source models 110 as a result of training. The training of each source parameter .theta..sub.i is done either individually, via procedures for training of each source model, or globally, using weighted sum formulae of likelihood scores with weighted coefficients .alpha..sub.i that are computed in block 150 (described below). The production of parameters generally is described in Bahl et al., "A Maximum Likelihood Approach to Continuous Speech Recognition", IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. PAMI-5, No. 2, pp. 179-180 (March 1983) which is hereby incorporated by reference in its entirety.

In particular, .theta..sub.i is a set of parameters (parameter vectors) that control a decoder (described below) for the i-th source of information. For example, if i corresponds to a speech source, then .theta..sub.i controls a specific model that is used in an automatic speech recognizer (ASR) in order to decode speech of a user. If the ASR is speaker-dependent then .theta..sub.i reflects voice characteristics of the specific person who trained the ASR.

The computed statistical parameters .theta..sub.i are used as an input into block 140. Block 140 computes scores (e.g., likelihood) F.sub..theta. for each source model. Block 150 in turn uses the computed scores and generates coefficients .alpha..sub.i. .alpha..sub.i is trained on the basis of one of the three procedures described in detail below. The value of .alpha..sub.i can be returned to block 130 via line 155 to further refine the statistical parameters .theta..sub.i. This loop is optional and does not have to be implemented.

Block 160 computes the weighted expression based on the computation of .alpha..sub.i and F.sub.i computed in blocks 150 and 140, respectively. Block 170 is a decoding module in which a search of the best candidate word (i.e., the word that maximizes the weighted sum) is performed.

As stated above, the automatic recognizer 100 of the present invention accepts multiple (N) sources of information. In particular, for 1.ltoreq.i.ltoreq.N, let F.sub.i (W;S.sub.i) be some likelihood score attached to a particular sub-message unit W, and let S.sub.i be the observed front-end sequence from the i-th source of information. Typically, W represents a word or a string of characters from some vocabulary, and S.sub.i represents a sequence of feature vectors obtained through suitable signal processing of the raw sensor data.

For each of the N sources of information taken alone, the standard decoding procedure is to select from a list of candidate units the unit W which maximizes the score F.sub.i (*,S.sub.i), where F.sub.i (*,S.sub.i) represents a function in which a different word W is substituted from the list of candidate words in the place of the * to find the word that maximizes the score F.sub.i (W;S.sub.i). Thus, during decoding the automatic recognizer 100 considers a few possible candidate words (or character strings) from the list of alternative words. This alternative list of words is produced either from some approximate fast method (see e.g., "A Fast Approximate Acoustic Match for Large Vocabulary Speech Recognition", Proc. 1989 Eurospeech Conference, Paris (September 1989)) or from perhaps the whole vocabulary which is stored in the recognizer.

However, in the combined use of several information sources considered in accordance with the present invention, the consistent message recognition system 100 maximizes a weighted sum containing all the scores F.sub.i (*,S.sub.i) simultaneously.

Thus, let .psi.(W;S) be the overall score to maximize, then: ##EQU1## is the expression for the weighted sum of scores considered here. For normalization purposes, it will be assumed that ##EQU2## Clearly, if each of the scores F.sub.i (*,S.sub.i) attains its maximum on the same sub-message unit, this unit is the best estimate for the piece of message jointly conveyed by all the S.sub.i. The problem becomes more difficult when these scores attain their maximum on different candidate units in the vocabulary, because the final decision will depend on the weighing factors .alpha..sub.i. In particular, these weighing factors should reflect, for example, the intrinsic quality of each source of information for each user, the distribution of likelihood scores over the set of candidate units, and the environment under which each S.sub.i evidence was produced.

.alpha..sub.i could be controlled manually by a user, according to his or her own intuitive estimate of the relevance of the associated source of information and the reliability of the corresponding recognition process. However, the present invention provides a number of procedures which lead to the automatic determination of explicit or approximate expressions for the weighing factors .alpha..sub.i, (i=1, . . . ,N).

Before proceeding, however, the main trade-offs involved in the maximization of an expression such as (1) should be examined. It seems that obtaining the weighing coefficients .alpha..sub.i is intimately related to obtaining a set of parameter vectors .theta..sub.i governing each individual score F.sub.i (W.vertline.S.sub.i), assuming a simultaneous maximization over all .alpha..sub.i and .theta..sub.i. In practice, however, such a maximization can lead to degenerative solutions for the .alpha..sub.i since the expression (1) is linear in .alpha..sub.i. In addition, forming a non-linear expression in lieu of (1) is rarely feasible from a computational point of view.

In contrast, however, the present invention assumes that, for i=1, . . . , N, some parameters .theta..sub.i.sup.(0) have been found by maximizing each individual score F.sub..theta.i (T.sub.i .vertline.S.sub.i) on some given training text T.sub.i (possibly common to all sources S.sub.i). It then becomes possible to find some suitable coefficients .alpha..sub.i .sup.(0) from (1) using algorithms to be discussed below. The next step is to compute a new set of parameters .theta..sub.i.sup.(1), i=1, . . . , N, to maximize the expression (1) given these coefficients .alpha..sub.i.sup.(0). Note that usually this can be done effectively if the computation of each .theta..sub.i.sup.(0) could be done effectively. For example, if each individual F.sub..theta. is a polynomial in .theta. with positive coefficients, the Baum-Welch technique can be used for the computation of a local maximum of F.sub..theta. (see L. E. Baum,. "An Inequality and Associated Maximization Technique in Statistical Estimation for Probabilistic functions of Markov Processes," Inequalities, Vol. 3, pp. 1-8, 1972). Since the weighted sum (1) will also be a polynomial with positive coefficients, the same efficient technique is applicable for the computation of the parameters .theta..sub.i.sup.(1).

The above process can then be iterated until a stable set of .alpha..sub.i can be isolated. This strategy has the advantage of decoupling the estimation of the weighting coefficients from the estimation of the vector parameters for each information source. As this latter estimation depends on the type of modelling considered for each source, the present disclosure will only consider the estimation of the weighing coefficients .alpha..sub.i.

Referring now to FIG. 2, a detailed illustration of the message recognizer of a preferred embodiment of the present invention is shown. The system 200 can be broken down into two different stages: training and decoding. Generally, the training stage accepts user training data and generates a weighted coefficient .alpha.. The decoding stage accepts a consistent message from multiple sources and produces a recognized message as its output.

Next, the message recognition system 200 will be described in detail.

System 200 includes a user interface 205 coupled to a display terminal 207 for displaying information to a user of the system 200. A pointing device 209, such as a mouse and/or a keyboard, may be provided for inputting additional information to the processor 205.

In accordance with the present invention the system 200 includes a set of interfaces 215 for providing a plurality of inputs to the automatic recognition system 200. The data (information) entered into the message recognition system 200 via interfaces 215 is also referred to as raw input data. The system 200 processes the raw input data into either digital or analog transformations of the input mediums.

The multiple input user interface is comprised of, for example, a handwriting transducer 215b, typically an electronic tablet and stylus, a speech transducer 215a, typically embodied as a microphone coupled to an analog to digital converter, and a generic transducer 215c, designated with a * to symbolize that the message recognition system 200 contemplates being capable of receiving a variety of different sources of information.

It is noted that the handwriting transducer 215b may be embodied in any means suitable for detecting a motion of a writing implement during use by a human subject. Also, the speech transducer 215a may include any means suitable for generating a set of patterns expressive of speech. In addition, a lipreading transducer (not shown) could also be coupled to the message recognition system 200 as described in Petajan, "Automatic Lipreading to Enhance Speech Recognition", CH2145-1, pps. 40-47, IEEE (1985) Petajan discusses a solid state camera that images a subject's lips and nostrils during speech. Moreover, touch patterns, brain impulses, and/or steno information, to name a few, could also be used as an input to the message recognition system 200.

With the use of a conventional electronic tablet for the handwriting transducer 215b data is typically acquired at 70 points per second. The sampled handwriting is applied to an associated handwriting feature vector (FV.sub.H) processor 220a. As employed herein a feature vector is a set of characteristics of the input accumulated during a specific time interval. Processor 220b generates one or more handwriting feature vectors for each input stroke. This is preferably accomplished, as indicated in FIGS. 3a and 3b, in reference to a center of gravity (CG) of a stroke. The input stroke, by example a lower-case "a", is captured by the electronic tablet. Each feature vector includes elements expressive of an x and a y displacement of a sample point (P) on the stroke, relative to the center of gravity, a displacement (b) of P relative to a baseline (BL), and angular information (theta) relating to a line (L) tangent to P. The resulting feature vector, for this presently preferred embodiment, has four elements and thus represents a four dimensional (4-d) vector.

With the use of a microphone and A/D converter as the speech transducer 215a the speech is typically sampled at a 20 kHz rate and is Fast Fourier transformed. The amplitude (a) within 20 frequency bands is determined as (a.sub.1, . . . a.sub.20) for each speech sample. The following articles are cited as teaching aspects of speech sampling and transformation that are suitable for use by the system 200: J. R. Cohen, "Application of an Auditory Model to Speech Recognition", J. Acoustic Soc. of America, Vol. 85, No. 6, pp. 2623-2629 (June 1989); A. Nadas, D. Nahamoo, and M. A. Picheny, "Speech Recognition Using Noise-Adaptive Prototypes", IEEE Trans. Acoustics, Speech, Signal Processing, Vol. 37 No. 10, pp. 1495-1503 (October 1989); and "Adaptive Labelling: Normalization of Speech by Adaptive Transformations Based on Vector Quantization", in Proc. 1988 Int. Conf. on Acoustics, Speech, Signal Processing, New York, N.Y. pp. 40-43 (April 1988).

FIG. 3c illustrates a speech sample being processed to obtain a speech feature vector (FV.sub.s) therefrom. The determination of amplitudes associated with the plurality of frequency bands yields a twenty dimensional (20-d) frequency domain feature vector. The output of the speech transducer 215a is applied to an associated FV.sub.s processor 220a that generates the series of 20-d speech feature vectors, wherein each element of the vector is expressive of an amplitude within one of the 20 sampling bands.

Once again, it is contemplated that other samples will be processed to obtain other feature vectors therefrom. Those skilled in the art would be in a position to generate the required feature vectors from other samples not discussed above.

During training of the system 200 the feature vectors are each applied to an associated word model training processor. Specifically, the speech vectors are applied to a speech model trainer (MT.sub.s) processor 230a, the handwriting feature vectors are applied to a handwriting model trainer (MT.sub.H) processor 230b, and the * feature vectors are applied to a * model trainer (MT.sub.*). Each of the model trainer processors may be coupled (not shown) to the user interface 205 for interacting with a user during training. This interaction may involve the entry of training sets of sentences by the user and the display of training results.

Each model trainer processor functions to determine a set of optimal parameters for an associated model, specifically a speech word model (WM.sub.s) 240a, a handwriting word model (WM.sub.H) 240b, and a * word model (WM.sub.*) 240c. In addition, a function of each of the training processors 230a and 230b is to determine optimal weights as shown in block 235. The optimal weights are represented as .alpha..sub.i (i=1, 2, . . . , N). .alpha..sub.i governs the relative contribution of speech, handwriting, and * to the final message recognition. For example, if during training on a specific set of training data speech is found to provide a higher likelihood, or score, than handwriting, the weight associated with the speech likelihood is increased and the weight associated with the handwriting likelihood is decreased. If the two likelihood scores are found to be comparable the weights are set to be the same.

In a presently preferred embodiment of the invention the speech trainer 230a operates in accordance with a hidden Markov model technique to identify phonemes and individual words. By example, for English speech there are four general classes of phonemes: nasals, stops, vowels and fricatives. The training and use of Markov models for speech recognition is disclosed in a journal article entitled "A Maximum Likelihood Approach to Continuous Speech Recognition", IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol. PAMI-5, No. 2, March 1983 pp. 179-190 by L. Bahl et al. The technique described in this article is one suitable technique for implementing the speech language model 260 (discussed below) and model trainer 230. In general, a hidden Markov model or source is comprised of a collection of states connected to one another by transitions which produce symbols from a finite alphabet. For training, one state is chosen as an initial state and another state is chosen as a final state. A training algorithm assigns probabilities to all strings of transitions from the initial state to the final state, as well as all outputs generated in the process.

The role of the Markov word model for speech is to represent the pronunciation of a word. The parameters of the Markov model are trained by relating the acoustic signal produced by the speaker to the word uttered. There is one Markov word model per word in the vocabulary.

The handwriting model trainer 230b preferably operates in accordance with an elastic stroke matching technique disclosed by C. C. Tappert in an article entitled "Adaptive On-Line Handwriting Recognition", IEEE International Conference on Pattern Recognition Proceedings (1984). As employed herein elastic matching clusters elastically deformed input strokes, as represented by their associated stroke vectors, into stroke templates to form character prototypes.

Coupled to the output of each of the models 240 is a Likelihood Estimator (LE), specifically, a speech LE 250a handwriting LE 250b, and * LE 250c. The general purpose of each of the LEs, in conjunction with their associated word models 240, is to assign probabilities to strings of words in response to input speech or handwriting signals. Each LE can be aided by a Language Model (LM), specifically a speech LM 260c, a handwriting LM 260c and * LM 260c, whose role is to reduce the search space by a priori discarding unlikely strings of words or characters. The LM is an optional feature of the present invention and does not have to be implemented to make the present invention operational.

After training, and during decoding of, for example, on-line speech and/or handwriting, the LEs 250 receive the feature vectors from blocks 220, and, in conjunction with the associated word models, and in some cases the associated LM scores, generate estimations of probable matching words. For example, the speech estimator 250a outputs a list (L.sub.1) of word estimations and the handwriting estimator 250b also outputs a list (L.sub.2) of word estimations.

The output of the Likelihood Merging (LM) block 280 is a list (L) of possible words corresponding to the acoustic, stroke, and/or * inputs. The list of words is ordered by probability. In general, the members of the list L are a function of SL.sub.1, SL.sub.2, and SL.sub.* and the associated weights .alpha..sub.i. The list of possible words is applied to a decoder 290 that selects the candidate word having the highest probability. This selected candidate word is considered to be a recognized message and becomes available to use by other applications including the user interface 205. These other applications may include, for example, a word processing or spreadsheet program.

It is within the scope of the invention to provide automatic pruning or trimming of the scaled lists SL.sub.1, SL.sub.2, and SL* such that L contains only those candidate words that are common to both SL.sub.1, SL.sub.2, and SL.sub.*. This method of trimming is readily accomplished by inspection. An initial trimming of the two word lists simplifies the task of decoder 290 in that fewer candidate words need be considered for likelihood merging. Similar to the LM scores, this is only an option and does not have to be implemented.

FIG. 5 illustrates a flowchart corresponding to a training procedure 500 for .alpha..sub.i. Procedure 500 represents block 150 and 235 described generally above in connection with FIG. 1 and FIG. 2, respectively.

Generally, the coefficients .alpha..sub.i are used to assign a weight to each source of information. .alpha..sub.i tends be higher for those scores function Fi(*, Si) that frequently provide a consistently correct decoding of the consistent message. For example, in noisy situations handwriting can be a more reliable source of information than speech and therefore .alpha..sub.i corresponding to the score function for handwriting recognition should be given a higher weight than if it the speech sample was produced in quiet conditions. However, in the same noisy condition if the handwriting sample was produced by a writer whose hands shake (from, for example, a disease) then even in the noisy situation, the speech sample might now be more relevant than the handwriting sample and .alpha..sub.i would be modified accordingly.

Similar considerations appear if a lipreading sample is entered into the system 200. Light conditions (dark v. bright) could affect the value of .alpha..sub.i that corresponds to lipreading source of information.

Referring in detail to FIG. 5, Block 530 provides a warping of the multiple messages. Input data from blocks 220a, 220b, 220c is aligned against each other in accordance with their time production, as well as other general characteristics. Specifically, the multiple inputs forming the consistent message are initially converted into feature vectors S.sub.i. The warping is done based on alignment of the feature vector sequence with the text of the message. A procedure called Viterbi alignment can be used, for example. Among other places, this algorithm is discussed in F. Jelinek, "The Development of an Experimental Discrete Dictation Recognizer", Proc. IEEE, Vol. 73, No. 11, pp. 1616-1624, November 1985. Once the alignment between the text and the evidence corresponding to each source of information is produced, the overall alignment between all evidence can be deduced.

An example of the warping procedure is shown in FIGS. 4a-c. For illustrative purposes only two sources of information are used: a spoken message and a written message. S.sub.1 represents the feature vector generated from the spoken source information and S.sub.2 represents the feature vector generated from the written source information. Both feature vectors are broken into units from 1 to .sigma.. Next, both feature vectors are aligned with respect to each other, thus creating a direct correlation between the units in the set of speech source information and the set of written source information.

Next, procedure 500 reorders the vectors S.sub.i based on the results of warping step 530, as shown in block 535. The sources of information (speech, handwriting, etc.) are ordered by some predefined criteria. For example, the information can be ordered based on the quality of a given source or by the level of its eventual contribution to the overall decoding. The information can also be ordered randomly.

As mentioned earlier, each .alpha..sub.i should vary in accordance with the degree of confidence that can be placed in the associated source of information and corresponding recognition process. This implies that they should reflect the effects of ambient noises, user-dependent characteristics such as age and sex, and recording equipment quality. Assuming, as in all statistical approaches, that the experimental conditions observed during training will also be encountered during decoding, such effects can be captured while enrolling each user.

Take for example the situation when a handwriting sample, speech sample, and lipreading sample are entered into the message recognition system 100 in relation to a consistent message T, and the quality of the three inputs from best to worst is speech, handwriting, and lipreading. In operation, reordering block 535 would order the three samples as: S.sup.(1) corresponding to the speech sample; S.sup.(2) corresponding to the handwriting sample; and S.sup.(3) corresponding to the lipreading sample.

Blocks 540, 542, 545 label the multiple sources of reordered information The N sources of reordered information are labeled S.sup.(1), S.sup.(2), . . . S.sup.(N), where the numbers that mark each source are referred to as ordered numbers. Next, each i-th source of labeled information is partitioned into K subunits: S.sub.1.sup.(i), S.sub.2.sup.(i), . . . , S.sub.k.sup.(i) (i=1, 2, . . . , N), as indicated in blocks 550, 552, and 555. These subunits of information correspond to pieces of messages T.sub.1, T.sub.2, . . . T.sub.K ; r-th subunit S.sub.r.sup.(i) corresponds to T.sub.r (r=1, 2, . . . , K).

In block 560, subunits from different sources of information with the same ordering number are aligned against each other in accordance with the results obtained in block 530. The alignment for each pair of subunits is induced by global alignment of all whole pieces of different sources of information in block 530.

Subunits of different sources of information corresponding to the same piece of a message are then grouped together in block 562 (regrouping block 562). This regrouping gives K groups of aligned subunits in each group of N aligned subunits.

As shown in block 565, each group of units S.sub.i.sup.(1), S.sub.i.sup.(2), . . . , S.sub.i.sup.(N) (i=1, 2, . . . , K) from the regrouping block 562 gives rise to a likelihood score L.sub.i (also referred to herein as a unit score). Specifically, each i-th group of pieces of source evidence S.sub.i(1), S.sub.i.sup.(2), . . . , S.sub.i.sup.(N) gives rise to computation of likelihood score:

L.sub.i =L(T.sub.i .vertline.S.sub.1.sup.(1), S.sub.i.sup.(2), . . . , S.sub.i.sup.(N)). This computation can be performed using standard methods (e.g., Viterbi and forward-backward algorithms). For example, one can divide a text T.sub.i in subunits U.sub.1, U.sub.2, . . . , U.sub.r and produce strings of labels for each piece of information evidence: f.sub.1.sup.(1) . . . f.sub.d.sup.(1) ; f.sub.1.sup.(2) . . . f.sub.e.sup.(2) ; . . . f.sub.1.sup.(N). L is computed approximately as: L (U.sub.1, U.sub.2, . . . , U.sub.r .vertline.f.sub.1.sup.(1) . . . f.sub.d.sup.(1) ; f.sub.1.sup.(2) . . . .sub.fe.sup.(2) ; . . . f.sub.1.sup.(N)) .fwdarw.L(U.sub.1, U.sub.2, . . . , U.sub.r .vertline..English Pound..sub.1.sup.(1) . . . .English Pound..sub.d.sup.(1) ; .English Pound..sub.d.sup.(2) . . . .English Pound..sub.e.sup.(2) ; . . . .English Pound..sub.f.sup.(N)).

The likelihood scores from block 565 produce a global score F.sub.(t) (T.vertline.S.sup.(1) . . . S.sup.(N)), as indicated in block 567. The global score can be approximately computed as a product of L.sub.1, L.sub.2, . . . , L.sub.k (as generated in block 565). However, in general more complex approximations can be considered.

The computation described in block 567 is stored, as indicated in block 570. Block 575 checks a list of orderings (that was previously prepared) to determine whether all the orderings have been seen. If they have not, then t is increased by one (block 585) and the computation of .alpha. is returned to block 535 for reordering. Otherwise, the procedure 500 proceeds to block 580.

The weighted coefficients .alpha..sup.(t) are computed using the scores F.sub.1 (*), F.sub.2 (*), . . . , F.sub.N (*) that are stored in block 570. After initial values for weighted coefficients are given and scores from block 567 are generated, gamma is computed as in equation (3) above and the derivatives are computed for each i, where i=1, 2, . . . , N. After all derivatives are computed the weighted coefficients are computed. The number of loops for computation of depends on a maximal number M or some other predetermined criteria.

Finally, if the number of iterations exceed some predetermined number N.sub.-- max or some equivalent criteria is met (e.g., .alpha. has not changed significantly with each new computation), then training of the weighted coefficients .alpha. is completed. Otherwise, the computation is returned to the reordering block 535. Before the procedure returns to block 535, however, the stored global scores are flushed, as shown in block 592, and t is set to zero.

By allowing the procedure to returned to reordering block 535, the source information can be reordered in block 535 and/or regrouped in block 562, thus providing a new set of data to calculate .alpha..

The following procedures estimate .alpha..sub.i as part of the training procedure, using data collected for training. The present invention provides three main classes of training procedures: (1) Deleted Interpolation, (2) Expectation-Maximization, and (3) Joint Maximum Likelihood. Only one of these procedures is used at any given time. Which procedure is used depends on a trade-off based on cost-performance.

The operation of these procedures will be next described in detail below.

(A) Deleted Interpolation

Let the training text T be split into K disjoint pieces T.sub.1, T.sub.2, . . . , T.sub.K, each representing a training subtext. Assuming independence, computing .psi., where .psi. stands for the overall score to maximize in (1), over the training text T is equivalent to computing the following expression: ##EQU3## where S.sub.i,k represent the evidence gathered through the i-th source of information using the subtext T.sub.k. Thus, from (1), the optimal values of .alpha..sub.i are those which maximize .psi.(.alpha..sub.i, . . . ,.alpha..sub.N).

Note that this expression is non-linear in .alpha..sub.i, thus getting around the degenerescence problem mentioned earlier. Using a well-known forward-backward algorithm (see L.E. Baum, "An Inequality and Associated Maximization Technique in Statistical Estimation for Probabilistic Functions of Markov Processes," Inequalities, Vol. 3, pp. 1-8 (1972), which is hereby incorporated by reference in its entirety), a local maximum of this expression can be found by starting from an initial guess .alpha..sub.i.sup.(0) (e.g., 1/N) and following the iterative procedure below, where r.gtoreq.0. ##EQU4##

Note that, providing a good starting point is selected, This procedure effectively converges to the global maximum.

The deleted interpolation procedure is shown in FIG. 6. Initially, as shown in block 610, .alpha. is initialized (either randomly or from heuristics); upper index (0) means initialization, and sub indexes 1, 2, . . . N account for the sources of information (N sources). In addition, score functions F.sub.1, F.sub.2, F.sub.3, . . . F.sub.N are also generated in block 610; each score function measures the likelihood of the corresponding source of information to have some textual data (T.sub.k) coupled with given pieces of source data (S.sub.i,k). Numbers i=1, 2, . . . N index the sources of information and k=1, 2, . . . K index different textual data that were used in training. Scores F.sub.1, F.sub.2, . . F.sub.N remain fixed in this procedure whereas .alpha..sub.1, .alpha..sub.2, .alpha..sub.N are changing iteratively.

Next, the computation of weight factors D.sub.i (i=1,2, . . . N) is done, as shown in block 620. Weight factors D.sub.i represent quantitative information measuring the contribution of previously computed coefficients i and score functions F.sub.i. These weight factors depend on the speed with which the total score function in (3) grows at their projection on the i-th coordinate of .alpha..sub.i (the higher this speed the larger the value of D.sub.i), as well as the current value of .alpha..sub.i (the larger .alpha..sub.i the larger D.sub.i).

Block 630 checks to determine if all sources of information have been used to compute the weight factors D.sub.i. If it hasn't then the procedure increments i and loops back to block 620. If all sources have been considered then the process flows to block 650.

Block 650 increments the iteration index r. Block 660 performs the next iterative step for computing .alpha..sub.i. The iterative values for .alpha..sub.i.sup.(r) corresponding to the iterative step r are computed as normalization of weight factors D.sub.i (i=1,2, . . . N). This procedure continues either until iterative values of .alpha..sub.i stabilize or a predetermined threshold value M, shown in block 670, is exceeded.

(B) Expectation-Maximization

More generally, a version of the general Expectation-Maximization (EM) procedure which is appropriate for mixtures such as defined in (1) can be used. This procedure is described in greater detail in R. A. Redner and H. F. Walker, "Mixture Densities, Maximum Likelihood and the EM Algorithm, " Siam Review, Vol. 26, No. 2, April 1984, which is hereby incorporated by reference in its entirety herein. With the same definitions as above, the recursive procedure is as follows: (a) Choose some initial .alpha.i.sup.(0) ; for example, set them all equal to 1/N and (b) compute for k=1,2, . . . ,K: ##EQU5##

Then update the weighing factors as: ##EQU6## where r.gtoreq.0. Again, the convergence of this procedure, while guaranteed, depends on the initial conditions selected for .alpha..sub.i.

The expectation-maximization procedure is shown in FIG. 7. Initially, as shown in block 710, the values for P (probability) in equation (6) are initialized, and in block 720, .alpha. is initialized.

Block 730 provides an iterative expectation computation of score p.sup.(r) based on the probability scores that were estimated in previous iterations. Block 740 checks to determine whether all probabilities of scores (one for each source of information) have been computed. If they have not, then i is incremented and block 730 is repeated. Otherwise procedure 700 proceeds to block 760.

Block 760 increments the iteration index r. Block 770 performs the iterative computation of .alpha.. .alpha. is computed as the average of the probability scores P that were computed in block 730.

Finally, block 780 checks to determine whether the number of iterations has exceeded a predetermined threshold.

(C) Joint Maximum Likelihood

Alternatively, we can arrive at a different procedure by starting from the expression for the joint likelihood. Note first that if only one source of information were present, say S.sub.i, then the task of a maximum a posteriori decoder would be to find the string of sub-message units W which maximizes P(W.vertline.S.sub.i). Using Bayes' rule and the fact that the maximization is performed on W, not S.sub.i, this operation can be written: ##EQU7## Thus, generally speaking, F.sub.i (W;S.sub.i) =P(S.sub.i .vertline.W) P(W).

When all the sequences S.sub.i corresponding to the N sources of information are available, the task of the maximum a posteriori decoder is to find the string of sub-message units W which maximizes P(W.vertline.S.sub.i, . . . S.sub.N). Using Bayes' rule and the chain rule, together with the fact that the maximization is performed on W, not on any of the S.sub.i, it can be written: ##EQU8## This expression has been obtained through the chain rule starting from S.sub.N and proceeding to S.sub.1. Among many other possibilities, N-1 additional ways can be selected to expand (9), starting from S.sub.N.sub.-1, S.sub.N-2, and so on in a round robin fashion. Putting all of these expressions together, the following expression is generated: ##EQU9## where by convention S.sub.0 is null. Relating this to (1), it becomes clear that the weighing factors .alpha..sub.i can be expressed in the following symmetric fashion: ##EQU10## for any 1.ltoreq.i.ltoreq.N. Furthermore, this is an exact expression for each .alpha..sub.i.

Each multiplicative term in (11), such as P(S.sub.j .vertline.S.sub.i, . . . , S.sub.j-1,W), for example, can be understood as follows. In an ideal situation, assuming that a sufficient amount of training data is available, the quantity P(S.sub.j .vertline.S.sub.i, . . . , S.sub.j-1,W) is obtained by training the model corresponding to the j-th source of information for all possible front-end sequences S.sub.j, when these sequences are observed at the same time as fixed front-end sequences produced by the other sources of information S.sub.i, . . . S.sub.j-1, using some common training text W. In practice, however, it is doubtful that a sufficient amount of training data could ever be collected for every source under every possible condition appearing in (11).

One possibility is to reduce the size of the sub-message unit to be recognized, so as to limit the number of potential front-end sequences that can be observed. Another way of achieving the same result would be to constrain the size of the front-end label alphabet by introducing equivalence classes on this alphabet. Alternatively, one may introduce some approximation in (11) so as to alleviate the heavy conditioning. For example, the model corresponding to the i-th source of information may be trained in the presence of only one (the most relevant) additional source of information. Assuming a judicious indexing of the sources of information, this leads to: ##EQU11## whereby convention S.sub.N+1 =S.sub.1. This expression is much easier to implement than (11), at the expense of some loss of optimality in the derivation of the .alpha.i.

There are numerous applications which could benefit from integrating multiple sources of information using the above three procedures. They can be roughly classified into three distinct categories: (1) Multiple Measurements, (2) Multiple Models, and (3) Divide and Conquer Methods.

The multiple measurement category includes speech and handwriting recognition, as described in U.S. patent application Ser. No. 07/676,601, now abandoned (cited above), as well as various enhancements thereof. For example, various syntactic features (nouns, verbs, etc.) of the language could be easily interpreted as multiple sources of information; lipreading information could be useful for improving the recognition of speech in noise; and hand movement information could perhaps make handwriting recognition more robust.

Ultimately, electrical impulses could be measured directly from the brain to reduce the amount of distortion introduced by the various transducers involved. Until such sensors can be made unintrusive enough, however, multiple solutions are expected to be especially useful to provide disabled people with more reliable means of communication.

The multiple model category assumes that for the same source of information several modelling strategies are available to carry out the message recognition. Examples include using both discrete and continuous parameter hidden Markov models to represent one set of acoustic data, as well as using both elastic matching and hidden Markov modelling to recognize a sequence of handwritten characters.

This means that we can compute several scores L.sub.i (W.vertline.S) for the same source information S. But this problem can be cast in the form of (1) by just letting L.sub.i (W.vertline.S)=F.sub.i (W.vertline.S.sub.i) where S.sub.i =S. This allows one to interpret the multiple use of the same data as multiple sources of information.

The divide and conquer method category is especially relevant to the construction of good statistical language models for use in various recognition and translation tasks. The standard approach is to derive, for any given sequence of words W.sub.1 . . . W.sub.M, the n-gram probabilities P.sub.(Wi-1, . . . Wi-n+1), i=1, . . . ,M. These probabilities are used for predicting the current word knowing the n-1 previous words. The main difficulty is that finding n-gram probabilities for n>3 requires too enormous a set of training texts to be practical.

In such a case, one could conceivably replace the estimation of the n-gram probabilities P.sub.(Wi-1, . . . Wi-n+1) by estimation of the following (N=n-2)-term weighted sum: ##EQU12##

In which each term P(W.sub.j .vertline.W.sub.j-1,W.sub.j-2), a 3-gram probability, is easy to calculate. It is clear that each of these terms can be interpreted as a separate source of information, and we can therefore set F.sub.j (W.sub.j ;S.sub.j)=P (W.sub.j .vertline.W.sub.j-1,W.sub.j-2) where S.sub.j =W.sub.j-1,W.sub.j-2. Thus, the procedures described previously can also be used for estimating the .alpha..sub.j coefficients in (13).

While the invention has been particularly shown and described with reference to preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention.

Claims

1. A computer program product for use with a computer system, comprising:

a computer usable medium having computer readable program code means embodied therein for causing a computer to transform a consistent message into a message recognizable by a computer, said computer program product having

first computer readable program code means for causing the computer to transform the consistent message generated by a human in at least two formats into electrical signal representations of the consistent message;

second computer readable program code means for causing the computer to produce from said electrical signal representations of the consistent message a set of parameters for each said format;

third computer readable program code means for causing the computer to generate a likelihood score of recognition from each said set of parameters;

fourth computer readable program code means for causing the computer to use said sets of parameters to train a weighting coefficient for each of the at least two formats of the consistent message, wherein said training of a weighting coefficient further comprises

partitioning means for partitioning information in each said format into one or more subunits, wherein each subunit corresponds to a piece of the message,

grouping means for grouping said subunits from each said format into a plurality of groups, wherein each said subunit in a group corresponds to the same piece of the message,

first determining means for determining a likelihood score of recognition for each of said groups of subunits,

second determining means for determining a global score for the information based on said likelihood score of recognition, and

third determining means for using said global score to determine said trained weighting coefficients;

fifth computer readable program code means for causing the computer to generate a weighted expression based on said trained weighting coefficients and said likelihood scores of recognition; and

sixth computer readable program code means for causing the computer to select a candidate message unit that maximizes said weighted expression to transform said electrical signal representations of the consistent message into a computer recognizable message.

2. A computer program product according to claim 1, wherein said trained weighting coefficients are used to produce a new set of parameters and said new set of parameters is applied to generate a new set of trained weighting coefficients.

3. A computer program product according to claim 2, wherein said new set of parameters is repeatedly produced until a stable set of trained weighting coefficients is isolated.

4. A computer program product according to claim 1, wherein the consistent message is transformed into said electrical signal representations using a transducer.

5. A computer program product according to claim 1, wherein said third determining means and said fourth determining means comprise means for generating and training said trained weighting coefficients using a deleted interpolation procedure.

6. A computer program product according to claim 5, wherein said third determining means comprises:

(a) means for initializing said trained weighting coefficients; and

(b) means for generating a weight factor for each of the at least two formats of the consistent message, whereby each of said weight factors represent quantitative information measuring the contribution of previously generated weighting coefficients; and

wherein said fourth determining means comprises means for generating said trained weighting coefficients based on a normalization of said weight factors for the at least two formats of the consistent message.

7. A computer program product according to claim 6, further comprising means for iteratively generating said weight factors and said trained weighting coefficients until said trained weighting coefficients stabilize or until a predetermined threshold is met.

8. A computer program product according to claim 1, wherein said third determining means and said fourth determining means comprise means for generating and training said trained weighting coefficients using an expectation-maximization procedure.

9. A computer program product according to claim 8, wherein said third determining means comprises:

(a) means for initializing a set of expectation probability scores and said trained weighting coefficients for each of the at least two formats of the consistent message;

(b) means for generating an expectation probability score based on previously generated expectation probability scores for each of the at least two formats of the consistent message; and

(c) means for generating said trained weighting coefficients as the average of said expectation probability scores; and

wherein said fourth determining means comprises means for iteratively generating said expectation probabilities and said training weight coefficients until a predetermined threshold is met.

10. A computer program product according to claim 1, wherein said third determining means and said fourth determining means comprise means for generating and training said trained weighting coefficients using a joint maximum likelihood procedure.

11. A computer program product according to claim 1, wherein said third determining means comprises:

(a) means for producing at least two sets of feature measurements from said at least two formats and organizing said set of measurements into at least two sets of feature vectors;

(b) means for providing a warping of said consistent message formed from said at least two formats;

(c) means for ordering said set of feature vectors based on some predefined criteria to form an ordered set of feature vectors;

(d) means for partitioning said ordered set of feature vectors to form a partitioned set of feature vectors;

(e) means for aligning said partitioned set of feature vectors in accordance with said warping step to form an aligned set of feature vectors;

(f) means for grouping said aligned set of feature vectors to form a grouped set of feature vectors, wherein each group of units from said grouped set of feature vectors gives rise to likelihood score; and

(g) means for producing a global score from said set of likelihood scores; and

wherein said fourth determining means comprises means for generating said trained weighting coefficients from said global score.

12. A computer program product according to claim 11, wherein said third determining means further comprises:

means for storing said global score;

means for reordering said set of feature vectors; and

wherein said fourth determining means further comprises means for generating said training weighting coefficients using said stored global score.

13. A computer program product according to claim 12, wherein said third determining means further comprises means for determining whether a predetermined set of criteria is met.

14. A program storage device readable by a machine, tangibly embodying a program of instructions executable by the machine to perform steps for transforming a consistent message into a message recognizable by a computer, said method steps comprising:

transforming the consistent message generated by a human in at least two formats into electrical signal representations of the consistent message;

producing from said electrical signal representations of the consistent message a set of parameters for each said format;

generating a likelihood score of recognition from each said set of parameters;

using said sets of parameters to train a weighting coefficient for each of the at least two formats of the consistent message, wherein said training comprises

partitioning information in each said format into one or more subunits, wherein each subunit corresponds to a piece of the message,

grouping said subunits from each said format into a plurality of groups, wherein each said subunit in a group corresponds to the same piece of the message,

determining a likelihood score of recognition for each of said groups of subunits,

determining a global score for the information based on said likelihood score, and

determining said trained weighting coefficients using said global score,

generating a weighted expression based on said trained weighting coefficient and said likelihood scores of recognition; and

selecting a candidate message unit that maximizes said weighted expressions to transform said electrical signal representations of the consistent message into a computer recognizable message.

15. A program storage device according to claim 14, wherein said trained weighting coefficients are used to produce a new set of parameters and said new set of parameters is applied to generate a new set of trained weighting coefficients.

16. A program storage device according to claim 15, wherein said new set of parameters is repeatedly produced until a stable set of trained weighting coefficients is isolated.

17. A program storage device according to claim 14, wherein the consistent message is transformed into said electrical signal representations using a transducer.

18. A program storage device according to claim 14, wherein said step of determining said trained weighting coefficients comprises a deleted interpolation procedure.

19. A program storage device according to claim 18, wherein said deleted interpolation procedure comprises the steps of:

(a) initializing said trained weighting coefficients;

(b) generating a weight factor for each of the at least two formats of the consistent message, whereby each of said weight factors represent quantitative information measuring the contribution of previously generated weighting coefficients; and

(c) generating said trained weighting coefficients based on a normalization of said weight factors for the at least two formats of the consistent message.

20. A program storage device according to claim 19, wherein steps of generating said weight factors and of generating said trained weighting coefficients are repeated until iterative values of said trained weighting coefficients stabilize or until a predetermined threshold is met.

21. A program storage device according to claim 14, wherein said step of determining said trained weighting coefficients comprises an expectation-maximization procedure.

22. A program storage device according to claim 21, wherein said expectation-maximization procedure comprises the steps of:

(a) initializing a set of expectation probability scores and said trained weighting coefficients for each of the at least two formats of the consistent message;

(b) generating an expectation probability score based on previously generated expectation probability scores for each of the at least two formats of the consistent message;

(c) generating said trained weighting coefficients as the average of said expectation probability scores; and

(d) repeating said step of generating said expectation probabilities and said step of generating said training weight coefficients until a predetermined threshold is met.

23. A program storage device according to claim 14, wherein said step of determining said trained weighting coefficients comprises a joint maximum likelihood procedure.

24. A program storage device according to claim 14, wherein said joint maximum likelihood procedure comprises the steps of:

(a) producing at least two sets of feature measurements from said at least two formats and organizing said set of measurements into at least two sets of feature vectors;

(b) providing a warping of said consistent message formed from said at least two formats;

(c) ordering said set of feature vectors based on some predefined criteria to form an ordered set of feature vectors;

(d) partitioning said ordered set of feature vectors to form a partitioned set of feature vectors;

(e) aligning said partitioned set of feature vectors in accordance with said warping step to form an aligned set of feature vectors;

(f) grouping said aligned set of feature vectors to form a grouped set of feature vectors, wherein each group of units from said grouped set of feature vectors gives rise to likelihood score;

(g) producing a global score from said set of likelihood scores; and

(h) generating said trained weighting coefficients from said global score.

25. A program storage device according to claim 24, wherein said joint maximum likelihood procedure further comprises the steps of:

storing said global score;

reordering said set of feature vectors; and

generating said training weighting coefficients using said stored global score.

26. A program storage device according to claim 25, wherein said joint maximum likelihood procedure further comprises the step of determining whether a predetermined set of criteria is met.