Device for recognizing an input pattern with approximate patterns used for reference patterns on mapping
BACKGROUND OF THE INVENTION 
     This invention relates to a device for recognizing an input pattern by the use of a predetermined number of reference patterns. Although equally well applicable to recognition of various patterns, such as type-printed or hand-printed letters, a device according to this invention will be described in the following mainly in connection with a speech recognition device. 
     A device for recognizing continuous speech sound of one or more spoken words and for encoding the result of recognition is advantageous as a device for supplying data and/or program words to an electronic digital computer and a device for supplying control data to various apparatus. The reasons are as follows. 
     First, the input operation may be carried out by an untrained person because it is only necessary to pronounce the input data rather than manipulating a keyboard or a like facility. Secondly, the input operation is possible with hands and feet used in accomplishing other purposes. In the third place, the input data may be supplied to the computer or the like even from a remote location merely through an ordinary telephone network. 
     Because of these merits, the speech recognition devices are widely demanded and have been developed at various places in the world into practical use. 
     In a speech recognition device, it is preferred in general to carry out pattern matching with the technique of dynamic programming resorted to as described in, for example, U.S. Pat. No. 3,816,722 issued to Hiroaki Sakoe and Seibi Chiba, assignors to the present assignee. In a speech recognition divece of this type, speech sound is subjected to spectrum analysis, sampling, and digitization and is thereby transformed into a time sequence of vectors representative of features of the speech sound at the respective sampling instants (hereafter referred to as a time sequence of feature vectors). The speech sound is representative of one or more continuously spoken words of a preselected vocabulary. The time sequence is representative of a speech sound pattern of the continuously spoken word or words. Prior to recognition of each speech sound pattern supplied to the device (hereafter named an input pattern), which is unknown to the device, at least one standard speech sound pattern for each word of the vocabulary (hereafter called at least one reference pattern) is supplied to the device and memorized therein. Comparison, namely, pattern matching, is carried out between the input pattern and every reference pattern with the dynamic programming technique resorted to. One of the reference pattern that is most similar to the input pattern is selected. The word represented by the selected reference pattern gives the result of recognition. 
     It is to be noted here that the input pattern is subject to a complicated and nonlinear deformation as regards the time axis as a result of variations in the speed of utterance as pointed out in the above-cited patent and also in an article contributed by Hiroaki Sakoe and Seibi Chiba to IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-26, No. 1 (February issue, 1978), pages 43-49, under the title of &#34;Dynamic Programming Algorithm Optimization for Spoken Word Recognition.&#34; Optimum pattern matching is achieved only after nonlinearly compensating for fluctuations or shifts between the time axes of the input pattern and the respective reference patterns. As will be discussed more in detail in the following, a considerable amount of calculation is necessary even with application thereto of the dynamic programming technique. High-priced calculators are indispensable in accomplishing the calculation within a reasonable interval of time. 
     Accurate or reliable speech recognition devices are therefore expensive. Low-priced ones are objectionable as to the performance. Conventional speech or pattern recognition devices are thus still defective in respect of performance-to-price ratio. 
     SUMMARY OF THE INVENTION 
     It is therefore an object of the present invention to provide a pattern recognition device that is excellent in the ratio of performance to price. 
     It is a more specific object of this invention to provide a pattern recognition device of the type described, which is not expensive but retains the capability of carrying out pattern matching and time normalization. 
     According to this invention, there is provided a pattern recognition device for recognizing an input pattern by calculating similarity measures between the input pattern and a predetermined number N of reference patterns. The device comprises input pattern buffer means for memorizing a time sequence of input pattern feature vectors a i  (i=1, 2, . . . , I) representative of the input pattern, reference pattern memory means for memorizing time sequences of reference pattern feature vectors b n   j  (n=1, 2, . . . , N; j=1, 2, . . . , J n ) representative of the respective reference patterns, standard vector memory means for memorizing standard vectors v m  (m=1, 2, . . . , F) preselected in consideration of the reference pattern feature vectors b n   j , and dictionary memory means for memorizing a dictionary given by suffixes for identifying selected ones v m  n ( j ) of the standard vectors v m . The selected standard vectors v m  n ( j ) are those of the standard vectors v m  which are most similar to the reference pattern feature vectors b n   j  of the respective reference patterns. The device further comprises time normalizing means for finding a mapping function j=j(i) for carrying out time normalization between the input pattern and each of the reference patterns by minimizing as a whole the difference between the input pattern feature vectors a i  and the standard vectors v m  ( j ) specified by the dictionary for the said each reference pattern and pattern matching means for carrying out pattern matching between the input pattern and the said each reference pattern based on the mapping function j=j(i) to derive a similarity measure between the input pattern and the said each reference pattern. 
     By substituting approximate patterns composed of the selected standard vectors for the respective reference patterns according to this invention, it is possible to astonishingly reduce the amount of calculation necessary on finding the mapping or warping function. 
    
    
     BRIEF DESCRIPTION OF THE DRAWING 
     FIG. 1 diagrammatically shows an input pattern and a reference pattern for describing the principles of the time normalization technique known in the art and used also in a pattern recognition device according to the instant invention; 
     FIG. 2 is a graph showing an input pattern, a reference pattern, a mapping function, and a normalization window; 
     FIG. 3 is a block diagram of a pattern recognition device according to an embodiment of this invention; 
     FIG. 4 is a block diagram of a pattern recognition device according to a practical embodiment of this invention; 
     FIG. 5 is a time chart for describing operation of the device illustrated in FIG. 4 in a first part of a time normalization step; 
     FIG. 6 is a block diagram of a time normalization unit for use in the device shown in FIG. 4; 
     FIG. 7 is a block diagram of a function calculator for use in the unit depicted in FIG. 6; 
     FIG. 8 is a time chart for describing operation of the device comprising the calculator illustrated in FIG. 7 in a latter part of the time normalization step and in a pattern maching step; and 
     FIG. 9 is a block diagram of a pattern matching unit for use in the device depicted in FIG. 4. 
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     Referring to FIGS. 1 and 2, the technique of pattern matching described in the above-referenced patent and article will briefly be reviewed at first for a better understanding of the present invention. Merely for simplicity of description, an input pattern A to be recognized will be restricted in the description to that representative of a spoken word. For recognition of each of various input patterns A&#39;s, first through N-th reference patterns B n  (n=1, 2, . . . , N) are prepared for a predetermined number of reference words of a vocabulary. Each word of the vocabulary may be pronounced or uttered in more than one manner of pronunciation on preparing the reference patterns B n . In any event, it is possible to differentiate the reference patterns B n  from one another by the affix n. In other words, the respective reference patterns B n  are identified by reference pattern-designating numbers n. When represented by the respective reference patterns B n , the reference words are also identified by the reference pattern-designating numbers n. 
     An input pattern A is represented by a time sequence of first through I-th input pattern feature vectors a i  (i=1, 2, . . . , I) arranged along a time axis i at euqally spaced i-th instants as: 
     
         A=a.sub.1, a.sub.2, . . . , a.sub.I,                       (1) 
    
     where the number I of the feature vectors a i  in the sequence, namely, the length of the sequence, will be called an input pattern length or duration. Like Equation (1), the respective reference patterns B n  are given by time sequences consisting of first through J n  -th reference pattern feature vectors b n   j  (j=1, 2, . . . , J n ), respectively, as: 
     
         B.sup.n =b.sup.n.sub.1, b.sup.n.sub.2, . . . , b.sup.n.sub.Jn, (2) 
    
     where J n  (written as Jn in a simplified fashion in Equations (2)) will be named reference pattern durations. Time axes j n  will be considered for the respective reference patterns B n  because the speeds of pronunciation of the input and the reference words may differ fron one another although the feature vectors a i  and b n   j  are spaced along the time axes i and j n  usually at an equal interval that is decided by a common sampling interval. Each of the reference patterns B n  will be represented by a time sequence of first through J-th reference pattern feature vectors b j  (j=1, 2, . . . , J) as: 
     
         B=b.sub.1, b.sub.2, . . . , b.sub.J. 
    
     The vectors a i , b n   j , and b j  and other vectors that will appear later are represented by the corresponding bold letters in the drawing and will be denoted by the usual letters in the following description unless the vector notations are desirable for some reason or other. 
     The feature vectors a i  and b j  are defined by first through Q-th vector components in a Q-dimensional space. The vector components are obtained by analysing each pattern A or B into a plurality of channels, Q in number, of different frequencies, such as logarithmically equally spaced frequencies, as: 
     
         a.sub.i =(a.sub.i,1, a.sub.i,2, . . . , a.sub.i,Q) 
    
     and 
     
         b.sub.j =(b.sub.j,1, b.sub.j,2, . . . , b.sub.j,Q), 
    
     where a i ,q or a iq  and b j ,q or b jq  (q=1, 2, . . . , Q) represent q-th vector components. It is possible to deal with the vector components either serially or in parallel. 
     In order to quantitatively represent the likelihood or similarity between an input pattern A and a reference pattern B, a time-normalized distance D(A, B) is defined as a similarity measure S(A, B) therebetween by: ##EQU1## where d(i, j) is defined, in turn, by: 
     
         d(i, j)=∥a.sub.i -b.sub.j ∥, 
    
     and is called a vector distance between one each of the input and the reference pattern feature vectors a i  and b j . In connection with the vector distances d(i, j), it may be mentioned here that the suffixes i and j are and will be used, as the case may be, to represent only one each of i and j defined hereinabove. This applies to other affixes used to differentiate similar quantitites from one another. In Equation (3), j=j(i) is a mapping or warping function for mapping from the time axis i of the input pattern A to the time j of the reference pattern B. With the mapping function decided as exemplified in FIG. 2 at 10, namely, selected so as to make Equation (3) hold, it is possible to eliminate the nonlinear shift or fluctuation along the time axis i relative to the time axis j or compensate for the shift between the axes i and j on an i-j plane. In other words, the mapping function j=j(i) is for determining the correspondence between the axes i and j. The time-normalized distance D(A, B) is thereby rendered substantially equal to zero when the reference pattern B corresponds to the input pattern A. Otherwise, the time-normalized distances D(A, B)&#39;s or D(A, B n ) have finite values in whichever manner the mapping functions are selected. 
     The known technique of dynamic programming may be resorted to on solving the minimization problem given by Equation (3). More particularly, the minimum sought for is found by successively calculating recurrence coefficients g(i, j), as called herein, by the use of a recurrence formula: 
     
         g(i, j)=d(i, j)+min g(i-1, j-p),                           (4) 
    
     p=O, 1, 2 
     in an ascending order for each of i and j. A direction coefficient p(i, j) of the mapping function j=j(i) for each value of j is varied from 0 to 2 for each value of i because it is usually sufficient on calculating the time-normalized distance D(A, B) to refer to only three previously calculated recurrence coefficients g(1-1, j), g(i-1, j-1), and g(i-1, j-2), which will be referred to as reference recurrence coefficients. The recurrence formula (4) is calculated, starting at an initial recurrence coefficient g(1, 1) to eventually reach an ultimate recurrence coefficient g(I, J). The initial recurrence coefficient is given by an initial condiciton: 
     
         g(1, 1)=d(1, 1).                                           (5) 
    
     In practice, it is sufficient that the recurrence formula (4) be calculated only within a normalization or adjustment window known in the art to be given by: 
     
         j-r≦i≦j+r, 
    
     where r represents a constant called a width or length of the window. The window is depicted in FIG. 2 between two boundary lines 11 and 12, the boundaries inclusive. 
     Even with the normalization window, it is necessary on calculating the recurrence formula (4) to calculate the vector distances d(i, j) as often as (2r+1)J times. On the other hand, the width of the normalization window r and the reference pattern durations j n  are about six and twenty, respectively, for the examples described in the article referred to hereinabove. The number of dimensions of the feature vectors a i  or b j  is usually about ten. It is thus necessary to calculate sums and differences about 2,600 times in order to minimize Equation (3) for each of N reference patterns B n . 
     Turning to FIG. 3, a pattern recognition device according to an embodiment of this invention is for recognizing an input pattern A by the use of a plurality of reference patterns B n  n=1, 2, . . . , N) and comprises an input pattern buffer 16 for temporarily memorizing a time sequence of input pattern feature vectors a i  given by Equation (1), a reference pattern memory 17 for memorizing time sequences of reference pattern feature vectors b n   j  given by Equations (2), and a standard vector memory 18 for memorizing a plurality of standard vectors v m  (m=1, 2, . . . , F) preselected in consideration of the reference pattern feature vectors b n   j . The number F of the standard vectors v m  may be considerably less than the total number of the reference pattern feature vectors b n   j . The suffixes m to the standard vectors v m  will be denoted also by k in the following as the case may be. 
     During a dictionary forming step designated in FIG. 3 by dashed lines and carried out as an off-line operation prior to supply of the input pattern A to the device, the reference pattern feature vectors b j  of each reference pattern B are compared with the standard vectors v m . Those approximate pattern feature vectors v m (j) are selected from the standard vectors v m  as indicated at 18&#39;, which are most similar to the feature vectors b j , respectively. A time sequence of first through J-th approximate pattern feature vectors v m (j) given by: 
     
         V=v.sub.m(1), v.sub.m(2), . . . , v.sub.m(j), 
    
     defines an approximate pattern V. With approximate patterns V n  thus formed of approximate pattern feature vectors v m  n.sub.(j)  (j=1, 2, . . . , J n ) for the respective reference patterns B n , only the suffixes m n  (j) of the approximate pattern feature vectors v m  n.sub.(j) of the respective approximate patterns V n  are stored as a dictionary M n  or M in a dictionary memory 19. 
     As soon as the input pattern A is supplied in FIG. 3 to the input pattern buffer 16, a time normalization step indicated by solid lines is carried out by a time normalization unit 21. At first, the approximate patterns V n  are successively reproduced in effect with reference to the dictionary M as shown at 21&#39;. The time normalization is carried out between the input pattern A and each reproduced approximate pattern V to find a mapping function j=j(i) defining the correspondence between the time axes i and j of the input pattern A and the reproduced approximate pattern V. Inasmuch as the approximate patterns V n  are nearly identical with the respective reference patterns B n , it is possible to consider that the mapping functions j n  =j n  (i) are found for the respective reference patterns B n . 
     Further referring to FIG. 3, next follows a pattern matching step specified by dash-dot lines. Responsive to the input pattern feature vectors a i  and those reference pattern feature vectors b j (i) of each reference pattern B which are related to the input pattern feature vectors a i  by the mapping function j=j(i) found for the approximate pattern V, a pattern matching unit 22 calculates the time-normalized distance D(A, B) by integrating the vector distances d(i, j)) with respect to i. The matching unit 22 thereby produces a similarity measure signal s representative of the time-normalized distance D(A, B). 
     Referring now to FIGS. 4 and 5, a pattern recognition device according to a practical embodiment of this invention comprises similar parts designated by like reference numerals and is operable by similar signals represented by like reference symbols. The device further comprises an input unit 25 for analysing a speech sound signal SP for an actually spoken word into the input pattern A. The input device 25 may be, for example, a spectrum analyser that is referred to by Louis C. W. Pols in an article he contributed to IEEE Transactions on Computers, Vol. C-20, No. 9 (September issue, 1971), pages 972-978, in FIG. 1, at band filters, log, detection, multiplexer and ADC. As soon as the input pattern A is stored in the input buffer 16 following the dictionary forming step mentioned hereinabove, a control unit 26 of the device produces for use in the time normalization and the pattern matching steps a mode signal e that is delivered to the time normalization unit 21 and indicates at first &#34;1&#34; to specify a distance calculation mode and then alternating &#34;2&#34;  and &#34;3&#34; for specifying a time normalization and a function calculation mode, respectively, during the time that the control unit 26 produces also a reference pattern designating signal n for designating the reference patterns B n  in a predetermined order. While the mode signal e specifies each of the modes, the control unit 26 produces an input pattern feature vector specifying signal i 1  variable in an ascending order from 1 to I to specify the input pattern feature vectors a i . Responsive to the signal i 1 , the input pattern buffer 16 supplies the normalization unit 21 with a first input pattern signal a sequentially representative of the feature vectors a i . 
     When compared with the device illustrated with reference to FIG. 3 and as will become clear as the description proceeds, the device being illustrated carries out the time normalization step in the distance calculation mode and the respective time normalization modes. The approximate patterns V n  are reproduced with reference to the dictionary memory 19 at the beginning of each time normalization mode. The matching step is carried out simultaneously with calculation of the mapping function j=j(i) in each function calculation mode. 
     While the mode signal e specifies in FIGS. 4 and 5 the distance calculation mode and while furthermore the input pattern feature vector specifying signal i 1  specifies each of the input pattern feature vectors a i , the control unit 26 produces a standard vector specifying signal k for specifying the suffixes 1 to F of the standard vectors v k  in a preselected order. In response to the standard vector specifying signals k, the standard vector memory 18 supplies the normalization unit 21 with a standard vector signal v that successively represent the standard vectors v k . Whenever the mode signal e is switched to the time normalization mode, the control unit 26 produces a suffix specifying signal j 1  that is variable in an ascending order preferably within the normalization window while the input pattern feature vector specifying signal i 1  specifies each of the input pattern feature vectors a i . Responsive to the reference pattern designating signal n and the suffix specifying signal j 1 , the dictionary memory 19 supplies the normalization unit 21 with a suffix signal m serially representative of this suffixes m(j) of the approximate pattern feature vectors v m (j) of the approximate pattern corresponding to the reference pattern B being designated by the reference pattern designating signal n. The input pattern feature vector specifying signal i 1  and either the standard vector specifying signal k or the suffix specifying signal j 1  are supplied also to the normalization unit 21. 
     Referring to FIG. 6 in addition to FIGS. 4 and 5, the time normalization unit 21 comprises a distance calculator 31 responsive to the input pattern and the standard vector signals a and v for calculating vector distances d&lt;i, k&gt; between all combinations (i, k) or (i, m) of the input pattern feature and the standard vectors a i  and v k  or v m  to produce an input distance signal d successively representative of the calculated vector distances d&lt;i, k&gt;. When the mode signal e indicates the distance calculation mode, a first multiplexer 32 selects the standard vector specifying signal k to produce the same as a first address signal m 1 . The distance signal d is stored in a distance memory matrix 33 in addresses specified by the input pattern feature vector specifying signal i 1  and the address signal m 1 . 
     Let it now be assumed in FIGS. 4 through 6 that the mode signal e specifies the time normalization mode, that the reference pattern designating signal n designates a reference pattern B, and that the input pattern feature vector specifying signal i 1  specifies the first input pattern feature vector a 1 . At the beginning of the time normalization mode, the suffix specifying signal j 1  specifies the suffix m(1) of the first approximate pattern feature vector v m (1) of the approximate pattern V corresponding to the designated reference pattern B. As described in general hereinabove, the suffix signal m produced by the dictionary memory represents the above-mentioned suffix m(1) under the circumstances. The multiplexer 32 selects the suffix signal m as the address signal m 1 . Accessed by the input pattern feature vector specifying and the address signals i 1  and m 1 , the distance memory matrix 33 produces an output distance signal d&#39; that represents at first a vector distance d&lt;1, m(1)&gt; between the first input pattern feature vector a 1  and the first approximate pattern feature vector v m (1) and is delivered to an adder-minimum finder 36 as the initial recurrence coefficient defined by Equation (5) and thence as an output recurrence coefficient signal g further to a working memory 37. Responsive to the input pattern feature vector and the suffix specifying signals i 1  and j 1 , the recurrence coefficient signal g is stored in the working memory 37. 
     When the input pattern feature vector and the suffix specifying signals i 1  and j 1  specify a general combination (i, j) for one each of the input and the approximate pattern feature vectors a i  and v m (j), the working memory 37 supplies the adder-minimum finder 36 with an input recurrence coefficient signal g&#39; successively representative of the reference recurrence coefficients g(i=1, j-p). The adder-minimum finder 36 finds the minimum of the reference recurrence coefficients g(i-1, j-p) and decides that direction coefficient p(i, j) to produce an input direction coefficient signal p representative thereof, for which the minimum is found. A second multiplexer 38 selects the suffix signal j 1  as a second address signal m 2  while the mode signal e indacates the time normalization mode and produces no second address signal in the distance calculation mode. The direction coefficient signal p is stored in a buffer memory 39 in an address specified by the input pattern feature vector specifying and the second address signals i 1  and m 2 . On the other hand, the suffix signal m supplied to the first multiplexer 32 now represents the suffix m(j). The distance memory matrix 33 therefore makes the output distance signal d&#39; represent a vector distance d&lt;i, m(j)&gt;. The adder-minimum finder 36 thus makes the output recurrence coefficient signal g represent the recurrence coefficient g(i, j) given by the recurrence formula (4). When the input pattern feature vector and the suffix specifying signals i 1  and j 1  eventually specify the combination (I, J), the mode signal e is switched to indicate the function calculation mode. 
     Referring now to FIGS. 7 and 8 in addition to FIGS. 4 and 6, the time normalization unit 21 comprises a function calculator 41 that comprises, in turn, a register 42. At the end of the time normalization mode, the suffix specifying signal j 1  that now specfies by J the suffix m(J) of the J-th approximate pattern feature vector v m (J) of the approximate pattern V corresponding to the reference pattern B being designated by the reference pattern designating signal n is registered in the register 42. As soon as the mode signal e begins to indicate the function calculation mode, the first multiplexer 32 stops production of the first address signal m 1 . Furthermore, the control unit 26 produces a reset pulse R that will be referred to later. During the function calculation mode, the unit 26 delivers a sequence of first clock pulses clk 1  to the register 42 and makes the input pattern feature vector specifying signal i 1  vary in a descending order, from I to 1, to specify the input pattern feature vectors a I , . . . , a 2 , and a 1 . The register 42 is for producing its content as a function signal j 2 , which is selected by the second multiplexer 38 as the second address signal m 2 . 
     When the input pattern feature vector specifying signal i 1  is representative of I, the register 42 makes the function signal j 2  represent J previously stored in the register 42. The value J of the signal j 2  represents the value of the mapping function j=j(i) for i=I, namely, j(I). Responsive to the input pattern feature vector specifying and the second address signals i 1  and m 2  representative of I and J, respectively, the buffer memory 39 produces an output direction coefficient signal p&#39; representative of the direction coefficient p(I, J). A subtractor 43 substracts the output direction coefficient signal p&#39; from the function signal j 2 . The difference is registered in the register 42 in response to build down of a first one of the first clock pulses clk 1 . The function signal j 2  represents j(I-1). 
     When the input pattern feature vector specifying signal i 1  represents a next value I-1, the output direction coefficient signal p&#39; represents p(I-1, j(I-1)). The difference represents j(I-1)-p(I-1, j(I-1)), which difference is registered in the register 42 as the value j(I-2) of the mapping function j=j(i). The processess are repeated until the input pattern feature vector specifying signal i 1  eventually represents 1. The function calculator 41 thus successively calculates: 
     
         j()=j(i-1)-p(i-1, j(i-1)). 
    
     In the meantime, the input pattern buffer 16 supplies the pattern matching unit 22 with a second input pattern signal a&#39; sequentially representative of the input pattern feature vectors a i  from the I-th one to the first one a I  to a 1 . Responsive to the function signal j 2 , the reference pattern memory 17 supplies the matching unit 22 with a reference pattern signal b successively representative of those reference pattern feature vectors b j (i) of the reference pattern B designated by the reference pattern designating signal n, from the J-th one b j (I) to the first one b j (1), which are optimally related to the I-th through the first input pattern feature vectors a I  through a 1 , respectively. 
     Referring to FIG. 9 in addition to FIGS. 4, 6, and 8, the pattern matching unit 22 carries out the pattern matching step while the mode signal e indicates the function calculation mode. Responsive to the reset pulse R mentioned hereinabove, an accumulator 46 of the matching unit 22 is reset to zero. Supplied with the second input pattern and the reference pattern signals a&#39; and b, a subtractor 47 of the matching unit 22 calculates the difference between the specified input pattern feature vector a i  and the optimally related reference pattern feature vector b j (i) to produce a difference signal ab. The most significant bit msb of the difference signal ab is supplied as an add-subtract control signal to the accumulator 46. When the difference is positive, the most significant bit msb representative of &#34;0&#34; makes the accumulator 46 carry out addition of the difference to the content thereof. When the difference is negative, the most significant bit msb of &#34;1&#34; makes the accumulator 46 subtract the difference from the content thereof. In this manner, the accumulator 46 integrates the absolute values of the differences. When the input pattern feature vector specifying signal i 1  varies from I eventually to 1, the accumulator 46 is loaded with a signal representative of the time-normalized distance D(A, B) given by Equation (3). The matching unit 22 produces a similarity measure signal s representative of the time-normalized distance D(A, B). In practice, the subtractor 47 serially calculates the differences between the vector components a iq  of each input pattern feature vector a i  and the corresponding components b j (i),q of the reference pattern feature vector b j (i) related to that input pattern feature vector a i  by the mapping function j=j(i). The differences are integrated in the accumulator 46 in response to a sequence of second clock pulses clk 2  supplied thereto from the control unit 26. 
     Turning back to FIG. 4, the device comprises a decision unit 49 that may be a microprocessor described, for example, on page 7 of &#34;The Am2900 Family Data Book&#34; published 1976 by Advanced Micro Devices, Inc., California, U.S.A. The similarity measure signals s&#39;s or s n  successively delivered to the decision unit 49 from the pattern matching unit 22 during the pattern matching steps for the respective reference patterns B n  are compared with one another. From the values of the reference pattern designating signal n supplied from the control unit 26, the decision unit 49 selects an optimum reference pattern-identifying number n for which the similarity measure signals s n  give the minimum time-normalized distance. 
     As pointed out hereinabove, it is necessary in prior art to calculate the vector distances d(i, j) as often as (2r+1)JN times on finding the mapping functions j n  =j n  (i) of calculating the time-normalized distances D(A, B n ) for N reference patterns. When the number N of the reference patterns B n  is two hundred, (2r+1)JN is as great as 52,000. In marked contrast, the vector distances d&lt;i, k&gt; are calculated FI times by the distance calculator 31 described with reference to FIGS. 5 and 6 and IN times by the matching unit 22. The total is given by (F+N)I, which is equal to only 6,000 when I≈J and F=100. 
     While this invention has thus for been described in conjunction with a few preferred embodiments thereof, it is now possible for those skilled in the art to put this invention into effect in various other ways. For example, minimization of the distances between the input pattern feature vectors a i  and the standard vectors v m (j) specified by the dictionary M for each reference pattern B on finding a mapping function j=j(i) is equivalent to maximization of the similarity between the said input pattern feature and standard vectors a i  and v m (j). The recurrence formula calculator comprising the adder-minimum finder 36 and the working and the buffer memories 37 and 39 may calculate according to the asymmetric dynamic programming algorithm given in the above-cited article, Table I, for P=1, a recurrence formula that is similar to Equation (4).