Abstract:
A digital filter of linear phase non-cyclic type comprises a shift register composed of serially connected delay elements with a specified delay time. The outputs from the individual delay elements are multiplied and added so that the output sampling frequency is twice that of the input sampling frequency which enters the shift register. The number of multiplications is significantly reduced by making the number of output lines from the delay elements to be odd and by placing adders between the shift register and the device for carrying out multiplications. An alternative method is to use read only memories instead of the means for carrying out multiplications.

Description:
This is a continuation of application Ser. No. 698,008 filed Feb. 4, 1985, now abandoned. 
    
    
     BACKGROUND OF THE INVENTION 
     This invention relates to a digital filter and, in particular, to a digital filter of linear phase non-recursive type. 
     A digital audio system generally makes use of an analog low pass filter in order to remove the high frequency components contained in the output of the digital-analog converter. FIG. 1 shows a circuit block diagram in such a case, a frequency characteristic diagram being shown in FIG. 2. In general, a low pass filter is required with a characteristic shown by the dotted line E so that the source signal A is passed but high frequency components B, C and D are attenuated. It is not a desirable method, however, to make use of a filter with such steep characteristics because distortions will be caused in the high frequency part of the audio signals. 
     For this reason, a digital filter is usually disposed before the digital-analog converter as shown in FIG. 3. If the higher harmonics component B which is closest to the normal playback range can be eliminated or if a digital filter with the characteristics shown by the dotted line F of FIG. 4 can be used, the result will be as shown in FIG. 5, there remaining only higher harmonics C which can be eliminated with an analog low pass filter with a relatively gentle frequency characteristic G. Thus, the distortion of waveforms can be minimized and high fidelity reproduction becomes possible. 
     FIG. 6 shows the circuit structure of a conventional digital filter of non-recursive type with its output sampling frequency twice as large as the input sampling frequency. In FIG. 6, delay elements D 1 . . . D 13  have a same delay time T, A 1 . . . A 14  are coefficients, m 1 . . . m 14  are multipliers and a is an adder. If the digital filter is of linear phase type, the following relationships hold among the coefficients: A 1  =A 14 , A 2  =A 13 , A 3  =A 12 , A 4  =A 11 , A 5  =A 10 , A 6  =A 9  and A 7  =A 8 . 
     The usual method of making the output sampling frequency of such a digital filter twice that of the input sampling frequency so that its characteristics will be as shown by the dotted line F of FIG. 4 (and that, for example, a band with its center at 44.1 KHz will be attenuated but another band with its center twice that value, or 88.2 KHz, as well as the source signal A will pass) is to insert zero data between the input data of the digital filter once every time period of T. If this method is used, the input data of even-numbered multipliers m 2 , m 4  . . . m 14  become 0 at certain time t and those of odd-numbered multipliers m 1 , m 3  . . . m 13  become 0 at time (t+T) as shown in FIG. 6. 
     In order to reduce the number of multipliers by a factor of 2, a circuit composition shown in FIG. 7 may now be considered. At time t, the seven multipliers m 1 . . . m 7  respectively execute multiplications of coefficients A 1 , A 3 , A 5 , A 7 , A 9 , A 11  and A 13  with the output data from the delay elements while at time (t+T) they execute multiplications of coefficients A 2 , A 4 , A 6 , A 8 , A 10 , A 12  and A 14  with the output data from the delay elements. Thus, by executing multiplications of odd-numbered and even-numbered coefficients alternately, the same output can be obtained as that by using the digital filter of FIG. 6. However, the number of multiplications cannot be reduced any further. 
     In general, multiplications take more time than additions and subtractions. If the time required for each multiplication is to be reduced, the hardware will become complicated and expensive. For this reason, it is extremely important to reduce the frequency of multiplications in the data processing of a digital audio system. 
     SUMMARY OF THE INVENTION 
     It is therefore an object of this invention to provide a method of reducing the frequency of multiplications as much as possible in the operation of a digital filter of a linear phase non-cyclic type. 
     A digital filter of linear phase non-cyclic type according to one embodiment of the present invention comprises a shift register with a specific number of serially connected delay elements and the outputs from these delay elements are added and multiplied so that the output sampling frequency will be twice the input sampling frequency which characterizes the input to the shift register, but this is achieved according to this embodiment by making the number of output lines of the delay elements from the shift register to be odd and by providing adders before (or at an earlier stage than) the multiplier means for performing multiplications. 
     According to another embodiment of the present invention, the multiplying means mentioned above is replaced by a read only memory (ROM). Described more specifically, another digital filter according to the present invention comprises a serial-type shift register with a longitudinally connected delay elements, a ROM for which address is selected by the output lines of these delay elements and accumulators for adding outputs from this ROM. The ROM is divided into segments and each segment has associated with it an accumulator which adds the output from the associated ROM segment. The filter output is obtained by adding the outputs from the individual accumulators. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is a block diagram describing the demodulating system of a digital audio unit; 
     FIG. 2 is a frequency characteristic diagram for showing modulation by this system; 
     FIG. 3 is a block diagram of the system of FIG. 1 with a digital filter; 
     FIGS. 4 and 5 are frequency characteristic diagrams for explaining the modulation of the system of FIG. 3; and 
     FIGS. 6 and 7 are circuit diagrams of conventional digital filters. 
     Embodiments of the present invention are shown in FIGS. 8-19. FIG. 8 is a circuit diagram according to a first embodiment; FIGS. 9 and 10 are diagrams presented for explaining the principle of the present invention; FIG. 11 is a time chart for explaining the operation according to the first embodiment; FIG. 12 is a circuit diagram according to a second embodiment; FIG. 13 is a circuit diagram according to a third embodiment; FIGS. 14 and 15 are circuit diagrams, respectively, for the delay elements DA and DB of FIG. 13; FIG. 16 is a time chart for explaining the operation according to the embodiment of FIG. 13; FIG. 17 is a circuit diagram according to a fourth embodiment of the present invention; FIG. 18 is a circuit diagram according to a fifth embodiment of the present invention; and FIG. 19 is a circuit diagram according to a sixth embodiment of the present invention. 
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
     A first embodiment of the present invention is shown in FIG. 8 wherein D 1  . . . D 7  are delay elements, φ 1  and φ 2  are two-phase clock signals, AD 1  . . . AD 3  are first-stage adders, m 1  . . . m 4  are multipliers, A 1  . . . A 7  and O are multiplicative coefficients. Symbol &#34;+&#34; in this and other drawings indicates an adder. In order to explain the operation of this embodiment, FIG. 8 is modified to FIGS. 9 and 10. FIG. 9 makes the comparison easier with the prior art example of FIG. 6, being different therefrom in that the final stage delay element D 13 , the multiplier m 14  related to the output thereof and the multiplicative coefficient A 14  are missing so that the number of output lines from the column of delay elements is odd. FIG. 10 results if input data and zero data are alternately inputted from FIG. 9 and multiplication regarding zero data is omitted. In FIG. 10, calculations are carried out in synchronization with the sampling frequency and multiplications are alternately done regarding odd-numbered coefficients (A 1 , A 3 , A 5 , A 7 , A 9 , A 11  and A 13 ) and even-numbered coefficients (A 2 , A 4 , A 6 , A 8 , A 10  and A 12 ). Since the digital filter is linear-phased, the following relationships hold: A 1  =A 13 , A 2  =A 12 , A 3  =A 11 , A 4  =A 10 , A 5  =A 9  and A 6  =A 8 . It is to be noted at this point that equalities hold always between two odd-numbered coefficients and between two even-numbered coefficients. This makes it possible to first add the pairs, each with matching timing for the execution of multiplication and having equal coefficients and then to carry out multiplications. 
     To explain the above more clearly, if the contents stored in the individual delay elements of FIG. 10 are M 1  . . . M 7 , the sum P of the products of odd-numbered coefficients is 
     
         P=M.sub.1 A.sub.l +M.sub.2 A.sub.3 +M.sub.3 A.sub.5 +M.sub.4 A.sub.7 +M.sub.5 A.sub.9 +M.sub.6 A.sub.11 +M.sub.7 A.sub.13. 
    
     Since A 1  =A 13 , A 3  =A 11  and A 5  =A 9 , we obtain 
     
         P=(M.sub.1 +M.sub.7)A.sub.1 +(M.sub.2 +M.sub.6)A.sub.3 +(M.sub.3 +M.sub.5)A5+M.sub.4 A.sub.7.                              (1) 
    
     Similarly, the sum Q of the products of even-numbered coefficients is 
     
         Q=M.sub.1 A.sub.2 +M.sub.2 A.sub.4 +M.sub.3 A.sub.6 +M.sub.4 A.sub.8 +M.sub.5 A.sub.10 +M.sub.6 A.sub.12. 
    
     But since A 2  =A 12 , A 4  =A 10  and A 6  =A 8 , we obtain 
     
         Q=(M.sub.1 +M.sub.6)A.sub.2 +(M.sub.2 +M.sub.5)A.sub.4 +(M.sub.3 +M.sub.4)A.sub.6                                          (2) 
    
     As can be understood by comparing (1) and (2), different terms are added in P and Q. According to the embodiment shown in FIG. 8, therefore, M 1  . . . M 4  and M 5  . . . M 7  are taken out in different timing by means of two-phase clock pulses φ 1  and φ 2 . 
     FIG. 11 is a time chart for the part of the calculation relating to coefficients A 1  and A 2  with regard to FIG. 8. Input signal A to the adder AD 1  is sequentially taken out in the order of . . . M 2 , M 1 , M 0  . . . in synchronization with clock signal φ 1  while the other input signal B is taken out in the order of . . . M 7 , M 6 , M 5  . . . in synchronization with clock signal φ 2 . Since the total output C from the adder AD 1  is (A+B), output is made sequentially in the order . . . (M 2  +M 7 ), (M 1  +M 7 ), (M 1  +M 6 ), (M 0  +M 6 ) , (M 0  +M 5 ) . . . and the multiplicative coefficients corresponding to them are inputted to the multiplier m 1 , alternating between A 1  and A 2  in synchronization with the clock signals φ 1  and φ 2 . As the result, the multiplier m 1  outputs alternately the first term of (1) and the first term of (2). Similar calculations are performed by the other multipliers m 2 , m 3  and m 4  so that an equivalent of P+Q is obtained as the output from the last-stage adder AD. 
     There is shown in FIG. 12 a second embodiment of the present invention which is different from the aforementioned first embodiment in that the column of delay elements are shifted by a common clock signal φ 1  and the input signals to the first-stage adders AD 1  . . . AD 3  are switched by the switches S 1 , S 2  and S 3  which are in synchronization with the timing for multiplication. 
     A third embodiment of the present invention is shown in FIG. 13, with the internal structures of the first and second delay elements DA and DB of FIG. 13 being respectively shown in FIGS. 14 and 15. Both the first and second delay elements DA and DB comprise a memory 11 having a predetermined delay time, a control input line 12, a gate circuit 13 for switching and a cycling circuit 14 so that the memory content can be both kept circulating and outputted serially if the control signal SEL is &#34;1&#34; and that the aforementioned cycling loop is cut off if SEL is &#34;0&#34;, the shift input data from a shift input line 15 being introduced into the memory 11 and the data originally stored in the memory being shifted to the next stage through a shift output line 16. A data output line 17 to a first-stage adder leads from the output terminal of the memory 11 in the case of a first delay element DA and from the input terminal of the memory 11 in the case of a second delay element DB. 
     FIG. 16 is a time chart for explaining the aforementioned third embodiment. In this example, data &#34;1234&#34; is first stored in the memory of a first delay element DA and is being outputted while SEL is &#34;1&#34;. Next, SEL becomes &#34;0&#34; and a new data &#34;5678&#34; is shift-inputted, sending the original data &#34;1234&#34; out. when SEL becomes &#34;1&#34; again, data &#34;5678&#34; is both maintained in the memory and outputted at the same time. In the meantime, data &#34;abcd&#34; is similarly stored in the memory of a second delay element DB and is being outputted while SEL is &#34;1&#34;. When SEL becomes &#34;0&#34; next, the new data &#34;efgh&#34; is shift-inputted, outputting the old data &#34;abcd&#34;. When SEL becomes &#34;1&#34; again, data &#34;efgh&#34; is both maintained in the memory and being outputted at the same time. In this situation, the signals through the two input lines E and F into the first-stage adder AD 1 , those through the output line G out of the adder AD 1  and the multiplicative coefficients A 1  and A 2  of the multiplier m 1  change as shown in FIG. 16. The same result of calculation can thus be obtained. This embodiment is convenient when data are serially processed. 
     A fourth embodiment of the present invention is shown in FIG. 17. This results by omitting the plurality of multipliers m 1  . . . m 14  and the adder a from the prior art design of FIG. 6 and by using a ROM 21 and an accumulator 22 instead according to a well-known method in this technology. All results of multiplications are stored in this ROM 21 so that when a ROM address is specified by a multiplier and a multiplicand, their product can be fetched immediately. 
     A fifth embodiment of the present invention is shown in FIG. 18 wherein the ROM of FIG. 17 is divided into three parts. Of the outputs X 1  . . . X 18  from the 18 delay elements SR 1  . . . SR 18 , X 1  . . . X 6  become the address of a first ROM 4A, X 7  . .. X 12  that of a second ROM 4B and X 13  . . . X 18  that of a third ROM 4C. These segmented ROMs 4A, 4B and 4C are respectively provided with accumulators 8A, 8B and 8C, each of which comprises a parallel adder ADD, a register R and a shifter S. The outputs Y A , Y B  and Y C  from the accumulators 8A, 8B and 8C are added by an adder 9 to yield an output Y. 
     The output Y A  from accumulator 8A is 
     
         Y.sub.A =α.sub.13 X.sub.6 +α.sub.14 X.sub.5 +. . . +α.sub.18 X.sub.1                                   (3) 
    
     
         but X.sub.n =Σ2.sup.j X.sub.n,j 
    
     where n=1, 2, . . . 6 and the summation is from j=0 to N-1 (j and N, respectively, representing the bit position of the data and the length of the data). Therefore, 
     
         Y.sub.A =Σ2.sup.j (α.sub.13 X.sub.6,j +α.sub.14 X.sub.5,j +. . . +α.sub.17 X.sub.2,j +α.sub.18 X.sub.1,j)(4) 
    
     Since X n ,j (n=1, . . . 6 and j=0, . . . N-1) can be either &#34;0&#34; or &#34;1&#34; and the α are constants, the sum inside the parentheses of (4) can take only 2 6  =64 different values. This means that a capacity of 64 words is all that is required of ROM 4A. The same holds also for the other two ROMs. Thus, the total capacity required of the three ROMs 4A, 4B and 4C of FIG. 18 is 3×64=192 words. 
     FIG. 19 discloses a sixth embodiment of the present invention according to which pairs of shift register outputs having same multiplicative coefficients such as X 1  and X 18 , X 2  and X 17 , X 3  and X 16 , X 4  and X 15 , X 5  and X 14 , X 6  and X 13 , X 7  and X 12 , X 8  and X 11 , and X 9  and X 10  are added first by first-stage adders AD 1  . . . AD 9  and these sums become the addresses of the divided ROMs 4D, 4E and 4F. The circuit composition after the ROM output is the same as shown in FIG. 18. 
     The output Y D  from the accumulator 8D which relates to ROM 4D is 
     
         Y.sub.D =α.sub.7 (X.sub.12 +X.sub.7)+α.sub.8 (X.sub.11 +X.sub.8)+α.sub.9 (X.sub.10 +X.sub.9)               (5) 
    
     
         but X.sub.n =Σ2.sup.j Xn,j 
    
     where n=7, . . . 12 and the summation is from j=0 to N-1 (j and N having the same meaning as explained above). Therefore, 
     
         Y.sub.D =Σ2.sup.j [α.sub.7 (X.sub.12,j +X.sub.7,j +C.sub.j-1)+α.sub.8 (X.sub.11,j +X.sub.8,j +C.sub.j -1)+α.sub.9 (X.sub.10,j +X.sub.9,j +C.sub.j-1)]     (6) 
    
     where C is the carry of the serial adder. Since the expression inside each pair of parentheses in (6) is either &#34;0&#34; or &#34;1&#34; and the α are constants, the expression inside the bracket in (6) can take only 2 3  =8 different values. This means that a capacity of 8 words is all that is required of ROM 4D. Since the same holds also for the other two ROMs, the total capacity required of the three ROMs is 3×8=24 words. 
     This is to be compared with an address of 2 18  =262144 words that would be required of a ROM according to the prior art design. This number is reduced to 192 by the present invention. If first-stage adders are provided, this can further be reduced to 24. 
     The present invention has been described above by a limited number of embodiments. The above description, however, should be considered as illustrative and not as limiting. For example, the number of serially connected delay elements may be varied and the ROM of FIG. 17 may be divided into a number of segments other than three. The scope of the present invention is limited only by the following claims.