DIGITAL SIGNAL PROCESSING SYSTEM AND DESIGN METHOD THEREOF

20190220498 ยท 2019-07-18

    Inventors

    Cpc classification

    International classification

    Abstract

    A digital signal processing system and its design method are disclosed. The digital signal processing system is a digital differentiator with a frequency response coefficient h.sub.m and a frequency response function H.sup.(K)(w), where

    [00001] H ( K ) ( w ) .Math. .Math. is .Math. .Math. one .Math. .Math. of .Math. .Math. { - jw K , 0 w < jw K , - w 0 , { jw K , 0 w < - jw K , - w 0 , .Math. - w K , - w < , and .Math. .Math. { jw K , 0 w < - jw K , - w 0 ,

    or combination thereof, where m of h.sub.m has a range 0mM1 and M is the sampling point quantity.

    Claims

    1. A method for designing a digital signal processing system, comprising (1) selecting an order K of the digital signal processing system, where K is an integer, (2) for odd-numbered K and when K - 1 2 is an odd number, setting the frequency response function H ( K ) ( w ) .Math. .Math. to .Math. { - jw K , 0 w < jw K , - w 0 , for odd-numbered K and when K - 1 2 is an even number, setting H.sup.(K)(w) to { jw K , 0 w < - jw K , - w 0 , for even-numbered K and when K 2 is an odd number, setting H.sup.(K)(w) to w.sup.K, w<; and for even-numbered K and when K 2 is an even number, setting H.sup.K(w) to { jw K , 0 w < - jw K , - w 0 , where w is frequency, and H.sup.(K)(w) is the value of Fourier transform (FT); (3) based on the frequency response function, setting the frequency response coefficient h.sub.m, where m has a range 0mM1 and M is the sampling point quantity; and (4) based on the frequency response function, determining a type of the digital signal processing system.

    2. The method according to claim 1, wherein, when the order K of the digital signal processing system is an odd number, K - 1 2 is an odd number, and M is an even odd, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( M 2 .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } .

    3. The method according to claim 1, wherein, when the order K of the digital signal processing system is an odd number, K - 1 2 is an odd number, and M is an even number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. } .

    4. The method according to claim 1, wherein, when the order K of the digital signal processing system is an odd number, K - 1 2 is an even number, and M is an odd number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( M 2 .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } .

    5. The method according to claim 1, wherein, when the order K of the digital signal processing system is an odd number, K - 1 2 is an even number, and M is an even number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. } .

    6. The method according to claim 1, wherein, when the order K of the digital signal processing system is an even number, K 2 is an odd number, and M is an odd number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } .

    7. The method according to claim 1, wherein, when the order K of the digital signal processing system is an even number, K 2 is an odd number, and M is an even number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) } .

    8. The method according to claim 1, wherein when the order K of the digital signal processing system is an even number, K 2 is an even number, and M is an odd number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 - ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } .

    9. The method according to claim 1, wherein, when the order K of the digital signal processing system is an even number, K 2 is an even number, and M is an even number, the frequency response coefficient h.sub.m is 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) } .

    10. A digital signal processing system having an order K where K is an integer, comprising a frequency response function H.sup.(K)(w), which is a Fourier transform (FT) function of order K at frequency w and a frequency response coefficient h.sub.m corresponding to the frequency response function, where 0mM1, M is the sampling point quantity, and H.sup.(K)(w) is one of { - jw K , 0 w < jw K , - w 0 , { jw K , 0 w < - jw K , - w 0 , - w K , - w < , and .Math. .Math. { jw K , 0 w < - jw K , - w 0 , or combination thereof.

    11. The digital signal processing system according to claim 10, further comprising at least a multiplier between an input terminal and an output terminal of the digital signal processing system, wherein each multiplier's amplification factor matches a corresponding frequency response coefficient h.sub.m.

    12. The digital signal processing system according to claim 10, wherein the digital signal processing system is one of a direct, series-connected, and linear-phased system or a combination thereof.

    13. The digital signal processing system according to claim 11, wherein the digital signal processing system is one of a direct, series-connected, and linear-phased system or a combination thereof.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    [0014] FIG. 1 is an amplitude response diagram for a digital differentiator according to prior art design method where the dots are points to be interpolated.

    [0015] FIG. 1A is an amplitude response diagram for an all-band digital diffeentiator with M=51, K=3 according to the present invention.

    [0016] FIG. 1B is an amplitude response diagram for an all-band digital differentiator with M=101, K=7 according to the present invention.

    [0017] FIG. 2A is an amplitude response diagram for an all-band digital differentiator with M=50, K=3 according to the present invention.

    [0018] FIG. 2B is an amplitude response diagram for an all-band digital differentiator with M=1001, K=7 according to the present invention.

    [0019] FIG. 3A is an amplitude response diagram for an all-band digital differentiator with M=51, K=1 according to the present invention.

    [0020] FIG. 3B is an amplitude response diagram for an all-band digital diffeentiator with M=151, K=5 according to the present invention.

    [0021] FIG. 4A is an amplitude response diagram for an all-band digital differentiator with M=50, K=1 according to the present invention.

    [0022] FIG. 4B is an amplitude response diagram for an all-band digital differentiator with M=80, K=5 according to the present invention.

    [0023] FIG. 5A is an amplitude response diagram for an all-band digital differentiator with M=51, K=2 according to the present invention.

    [0024] FIG. 5B is an amplitude response diagram for an all-band digital differentiator with M=85, K=6 according to the present invention.

    [0025] FIG. 6A is an amplitude response diagram for an all-bend digital differentiator with M=50, K=2 according to the present invention.

    [0026] FIG. 6B is an amplitude response diagram for an all-band digital differntiator with M=150, K=6 according to the present invention.

    [0027] FIG. 7A is an amplitude response diagram for an all-band digital diffeentiator with M=51, K=2 according to the present invention.

    [0028] FIG. 7B is an amplitude response diagram for an all-band digital differentiator with M=201, K=6 according to the present invention.

    [0029] FIG. 8A is an amplitude response diagram for an all-band digital diffeantiator with M=50, K=2 according to the present invention.

    [0030] FIG. 8B is an amplitude response diagram for an all-band digital differentiator with M=200, K=6 according to the present invention.

    [0031] FIG. 9 is an amplitude response diagram for a partial-band digital differentiator with M=201, N.sub.1=81, K=3 according to the present invention.

    [0032] FIG. 10 is an amplitude response diagram for a partial-band digital differentiator with M=451, N.sub.1=91, K=7 according to the present invention.

    [0033] FIG. 11 is an amplitude response diagram for a partial-band digital differentiator with M=100, N.sub.1=35, K=3 according to the present invention.

    [0034] FIG. 12 is an amplitude response diagram for a partial-band digital differentiator with M=151, N.sub.1=57, K=1 according to the present invention.

    [0035] FIG. 13 is an amplitude response diagram for a partial-band digital differentiator with M=100, N.sub.1=30, K=5 according to the present invention.

    [0036] FIG. 14 is an amplitude response diagram for a partial-band digital differentiator with M=101, N.sub.1=41, K=2 according to the present invention.

    [0037] FIG. 15 is an amplitude response diagram for a partial-band digital differentiator with M=200, N.sub.1=60, K=6 according to the present invention.

    [0038] FIG. 16 is an amplitude response diagram for a partial-bend digital differentiator with M=201, N.sub.1=86, K=4 according to the present invention.

    [0039] FIG. 17 is an amplitude response diagram for a partial-band digital differentiator with M=100, N.sub.1=48, K=8 according to the present invention.

    [0040] FIG. 18 is a transfer function block diagram for a linear-phased digital differentiator with M=7.

    [0041] FIG. 19 is an operational block diagram of a direct type FIR differantiator.

    DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

    [0042] The following descriptions are exemplary embodiments only, and are not intended to limit the scope, applicability or configuration of the invention in any way. Rather, the following description provides a convenient illustration for implementing exemplary embodiments of the invention. Various changes to the described embodiments may be made in the function and arrangement of the elements described without departing from the scope of the invention as set forth in the appended claims.

    [0043] When discrete Fourier transform H(m) may be expressed as

    [00017] H ( m ) = H r ( 2 .Math. .Math. .Math. m M ) .Math. e j .Math. .Math. .Math. .Math. H ( m ) ,

    H(m) is the phase function, the amplitude response

    [00018] .Math. H r ( 2 .Math. .Math. .Math. m M ) .Math.

    is a real-coefficient function, and

    [00019] H r ( 2 .Math. .Math. .Math. m M ) = { H r ( 0 ) , m = 0 H ( 2 .Math. ( M - m ) M ) , 1 m M - 1

    As the FIR digital signal processing system to be designed is linear-phased, the phase angle H(m) is as follows if H(m) is symmetric:

    [00020] .Math. .Math. H ( m ) = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ;

    the phase angle H(m) is as follows if H(m) is symmetric:

    [00021] .Math. .Math. H ( m ) = { 2 - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - ( 2 ) + ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ;

    wherein m is the largest integer less than or equal to m. Additionally, if H(m) is symmetric, H.sub.m=H(m) and

    [00022] .Math. .Math. H m = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ;

    whereas if H(m) is anti-symmetric, H.sub.m=jH(m) and

    [00023] .Math. .Math. H m = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .

    Assuming that

    [00024] M = 2 .Math. .Math. M , = M - 1 2 , and .Math. .Math. _ M = ( M - 1 2 ) .Math. 2 .Math. .Math. M = .Math. .Math. M ,

    then

    [00025] W * .Math. H = .Math. [ 1 1 1 .Math. 1 1 .Math. 1 1 1 e j .Math. .Math. M e j .Math. .Math. 2 .Math. .Math. M .Math. e j .Math. M - 1 2 .Math. M e j .Math. M + 1 2 .Math. M .Math. e j ( M - 2 ) .Math. M e j ( M - 1 ) .Math. M 1 e 2 .Math. .Math. j .Math. .Math. M e 2 .Math. .Math. j .Math. .Math. 2 .Math. .Math. M .Math. e 2 .Math. .Math. j .Math. M - 1 2 .Math. M e 2 .Math. .Math. j .Math. M + 1 2 .Math. M .Math. e 2 .Math. .Math. j ( M - 2 ) .Math. M e 2 .Math. .Math. j ( M - 1 ) .Math. M .Math. .Math. 1 e ( M - 2 ) .Math. j .Math. .Math. M e ( M - 2 ) .Math. j .Math. .Math. 2 .Math. .Math. M .Math. e ( M - 2 ) .Math. j .Math. M - 1 2 .Math. M e ( M - 2 ) .Math. j .Math. M + 1 2 .Math. M .Math. e ( M - 2 ) .Math. j ( M - 2 ) .Math. M e ( M - 2 ) .Math. j ( M - 1 ) .Math. M 1 e ( M - 1 ) .Math. j .Math. .Math. M e ( M - 1 ) .Math. j .Math. .Math. 2 .Math. .Math. M .Math. e ( M - 1 ) .Math. j .Math. M - 1 2 .Math. M e ( M - 1 ) .Math. j .Math. M + 1 2 .Math. M .Math. e ( M - 1 ) .Math. j ( M - 2 ) .Math. M e ( M - 1 ) .Math. j ( M - 1 ) .Math. M ] .Math. [ H 0 .Math. e 0 .Math. .Math. j .Math. .Math. _ M H 1 .Math. e - j .Math. .Math. _ M H 2 .Math. e - 2 .Math. .Math. j .Math. .Math. _ M .Math. H M - 1 2 .Math. e - M - 1 2 .Math. j .Math. .Math. _ M H M - 1 2 .Math. e M + 1 2 .Math. j .Math. .Math. _ M .Math. H M - 2 .Math. e 2 .Math. .Math. j .Math. .Math. _ M H M - 1 .Math. e j .Math. .Math. _ M ] = .Math. [ H 0 + H 1 .Math. e - j .Math. .Math. _ M + H 2 .Math. e - 2 .Math. .Math. j .Math. .Math. _ M + .Math. + H M - 2 .Math. e 2 .Math. .Math. j .Math. .Math. _ M + H M - 1 .Math. e j .Math. .Math. _ M H 0 + H 1 .Math. e j ( M - _ M ) + H 2 .Math. e j ( 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + H M - 2 .Math. e j ( ( M - 2 ) .Math. M + 2 .Math. .Math. _ M ) + H M - 1 .Math. e j ( ( M - 1 ) .Math. M + _ M ) H 0 + H 1 .Math. e j ( 2 .Math. .Math. M - _ M ) + H 2 .Math. e j ( 2 .Math. 2 .Math. .Math. M - _ M ) + .Math. + H M - 2 .Math. e j ( 2 .Math. ( M - 2 ) .Math. M + 2 .Math. .Math. _ M ) + H M - 1 .Math. e j ( 2 .Math. ( M - 1 ) .Math. M + _ M ) .Math. H 0 + H 1 .Math. e j ( ( M - 2 ) .Math. M - _ M ) + H 2 .Math. e j ( 2 .Math. ( M - 2 ) .Math. M - 2 .Math. _ M ) + .Math. + H M - 2 .Math. e j ( ( M - 2 ) .Math. ( M - 2 ) .Math. M + 2 .Math. .Math. _ M ) + H M - 1 .Math. e j ( ( M - 1 ) .Math. ( M - 2 ) .Math. M + _ M ) H 0 + H 1 .Math. e j ( ( M - 1 ) .Math. M - _ M ) + H 2 .Math. e j ( 2 .Math. ( M - 1 ) .Math. M - 2 .Math. _ M ) + .Math. + H M - 2 .Math. e j ( ( M - 2 ) .Math. ( M - 1 ) .Math. M + 2 .Math. .Math. _ M ] + H M - 1 .Math. e j ( ( M - 1 ) .Math. ( M - 1 ) .Math. M + _ M ) ] ( 1 )

    If the digital processing system to be designed is a FIR digital differentiator, H.sub.m and W*H is identical to the above, and details are therefore omitted.

    [0044] If H(m) is symmetric and M is an odd number, W*H may be obtained using equation (1) as follows:

    [00026] W * .Math. H = .Math. [ H 0 + ( H 1 .Math. e - j .Math. .Math. _ M + H M - 1 .Math. e j .Math. .Math. _ M ) + ( H 2 .Math. e - 2 .Math. .Math. j .Math. .Math. _ M + H M - 2 .Math. e 2 .Math. .Math. j .Math. .Math. _ M ) + .Math. + ( H M - 1 2 .Math. e - j .Math. .Math. M - 1 2 .Math. _ M + H M + 1 2 .Math. e j .Math. M - 1 2 .Math. .Math. _ M ) H 0 + ( H 1 .Math. e j ( M - _ M ) + H M - 1 .Math. e j [ ( M - 1 ) .Math. M + _ M ] ) + .Math. + ( H M - 1 2 .Math. e j ( M - 1 2 .Math. M - M - 1 2 .Math. _ M ) + H M + 1 2 .Math. e j [ M + 1 2 .Math. M + M - 1 2 .Math. _ M ] ) H 0 + ( H 1 .Math. e j ( 2 .Math. .Math. M - _ M ) + H M - 1 .Math. e j [ 2 .Math. ( M - 1 ) .Math. M - _ M ] ) + .Math. + ( H M - 1 2 .Math. e j ( 2 .Math. M - 1 2 .Math. M - M - 1 2 .Math. _ M ) + H M + 1 2 .Math. e j [ 2 .Math. M - 1 2 .Math. M + M - 1 2 .Math. _ M ] ) .Math. H 0 + ( H 1 .Math. e j ( ( M - 2 ) .Math. M - _ M ) + H M - 1 .Math. e j [ ( M - 2 ) .Math. ( M - 1 ) .Math. M + _ M ] ) + .Math. + ( H M - 1 2 .Math. e j [ M - 1 2 .Math. ( M - 2 ) .Math. M - M - 1 2 .Math. .Math. _ M ] + H M + 1 2 .Math. e j [ M - 1 2 .Math. ( M - 2 ) .Math. M + M - 1 2 .Math. _ M ] ) H 0 + ( H 1 .Math. e j ( ( M - 1 ) .Math. M - _ M ) + H M - 1 .Math. e j [ ( M - 1 ) .Math. ( M - 1 ) .Math. M + _ M ] ) + .Math. + ( H M - 1 2 .Math. e j [ M - 1 2 .Math. ( M - 1 ) .Math. M - M - 1 2 .Math. .Math. _ M ] + H M + 1 2 .Math. e j [ M + 1 2 .Math. ( M - 1 ) .Math. M + M - 1 2 .Math. _ M ] ) ] = .Math. [ 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos .Math. .Math. _ M + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos .Math. .Math. 2 .Math. .Math. _ M + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. cos .Math. M - 1 2 .Math. _ M 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. cos ( M - 1 2 .Math. M - M - 1 2 .Math. _ M ) 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( 2 .Math. .Math. M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. cos ( 2 .Math. M - 1 2 .Math. M - M - 1 2 .Math. _ M ) .Math. 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( ( M - 2 ) .Math. M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. ( M - 2 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. cos ( M - 1 2 .Math. ( M - 2 ) .Math. M - 1 2 .Math. M - M - 1 2 .Math. _ M ) 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( ( M - 1 ) .Math. M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. ( M - 1 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. cos ( M - 1 2 .Math. ( M - 1 ) .Math. M - 1 2 .Math. M - M - 1 2 .Math. _ M ) ] ( 2 )

    When M is an even number, W*H may be obtained using equation (1) as follows:

    [00027] [ H 0 + ( H 1 .Math. e - j .Math. .Math. _ M + H M - 1 .Math. e j .Math. .Math. _ M ) + .Math. + ( H M 2 - 1 .Math. e - j ( .Math. M 2 - 1 ) .Math. _ M + H M 2 - 1 .Math. e j ( M 2 - 1 ) .Math. .Math. _ M ) H 0 + ( H 1 .Math. e j ( M - _ M ) + H M - 1 .Math. e j ( ( M - 1 ) .Math. M + _ M ) ) + .Math. + ( H M 2 - 1 .Math. e j [ ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ] + H M 2 + 1 .Math. e j [ ( M 2 + 1 ) .Math. M + ( M 2 - 1 ) .Math. _ M ] ) H 0 + ( H 1 .Math. e j ( 2 .Math. .Math. M - _ M ) + H M - 1 .Math. e j ( 2 .Math. ( M - 1 ) .Math. M + _ M ) ) + .Math. + ( H M 2 - 1 .Math. e j [ ( M 2 - 1 ) .Math. 2 .Math. .Math. M - ( M 2 - 1 ) .Math. _ M ] + H M 2 + 1 .Math. e j [ ( M 2 + 1 ) .Math. 2 .Math. .Math. M + ( M 2 - 1 ) .Math. _ M ] ) .Math. H 0 + ( H 1 .Math. e j ( ( M - 2 ) .Math. M - _ M ) + H M - 1 .Math. e j ( ( M - 2 ) .Math. ( M - 1 ) .Math. M + _ M ) ) + .Math. + ( H M 2 - 1 .Math. e j [ ( M 2 - 1 ) .Math. ( M - 2 ) .Math. M - ( M 2 - 1 ) .Math. .Math. _ M ] + H M 2 + 1 .Math. e j [ ( M 2 - 1 ) .Math. ( M - 2 ) .Math. M + ( M 2 - 1 ) .Math. _ M ] ) H 0 + ( H 1 .Math. e j ( ( M - 1 ) .Math. M - _ M ) + H M - 1 .Math. e j ( ( M - 1 ) .Math. ( M - 1 ) .Math. M + _ M ) ) + .Math. + ( H M 2 - 1 .Math. e j [ ( M 2 - 1 ) .Math. ( M - 1 ) .Math. M - ( M 2 - 1 ) .Math. .Math. _ M ] + H M 2 + 1 .Math. e j [ ( M 2 + 1 ) .Math. ( M - 1 ) .Math. M + ( M 2 - 1 ) .Math. _ M ] ) ] + [ H M 2 .Math. e j .Math. M 2 .Math. _ M H M 2 .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) H M 2 .Math. e j ( M 2 .Math. 2 .Math. .Math. M + M 2 .Math. _ M ) .Math. H M 2 .Math. e j ( M 2 .Math. ( M - 2 ) .Math. M + M 2 .Math. _ M ) H M 2 .Math. e j ( M 2 .Math. ( M - 1 ) .Math. M + M 2 .Math. _ M ) ] = [ 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos .Math. .Math. _ M + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos .Math. .Math. 2 .Math. .Math. _ M + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M 2 - 1 M ) K .Math. cos ( M 2 - 1 ) .Math. _ M 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M 2 - 1 M ) K .Math. cos ( ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( 2 .Math. .Math. M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M 2 - 1 M ) K .Math. cos ( 2 .Math. ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) .Math. 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( ( M - 2 ) .Math. M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. ( M - 2 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M 2 - 1 M ) K .Math. cos ( ( M - 2 ) .Math. ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) 2 .Math. ( 2 .Math. .Math. M ) K .Math. cos ( ( M - 1 ) .Math. M - _ M ) + 2 .Math. ( 2 .Math. .Math. .Math. 2 M ) K .Math. cos ( 2 .Math. ( M - 1 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. ( 2 .Math. .Math. .Math. M 2 - 1 M ) K .Math. cos ( ( M - 1 ) .Math. ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) ] + [ ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j .Math. M 2 .Math. _ M ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. 2 .Math. .Math. M + M 2 .Math. _ M ) .Math. ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 2 ) .Math. M + M 2 .Math. _ M ) ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 1 ) .Math. M + M 2 .Math. _ M ) ] ( 3 )

    [0045] If H(m) is anti-symnmetric and M is an odd number, W*H may be obtained usinig equation (1) as follows:

    [00028] W * .Math. H == [ - 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin .Math. .Math. _ M .Math. - 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin .Math. .Math. 2 .Math. .Math. _ M - .Math. - 2 .Math. .Math. j ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. sin .Math. M - 1 2 .Math. _ M 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. sin ( M - 1 2 .Math. M - M - 1 2 .Math. _ M ) 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( 2 .Math. .Math. M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. sin ( 2 .Math. M - 1 2 .Math. M - M - 1 2 .Math. _ M ) .Math. 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( ( M - 2 ) .Math. M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. ( M - 2 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. sin ( M - 1 2 .Math. ( M - 2 ) .Math. M - M - 1 2 .Math. _ M ) 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( ( M - 1 ) .Math. M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. ( M - 1 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. M - 1 2 M ) K .Math. sin ( M - 1 2 .Math. ( M - 1 ) .Math. M - M - 1 2 .Math. _ M ) ] ( 4 )

    When M is an even number, W*H may be obtained using equation (1) as follows:

    [00029] = [ - 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin .Math. .Math. _ M .Math. - 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin .Math. .Math. 2 .Math. .Math. _ M - .Math. - 2 .Math. .Math. j ( 2 .Math. .Math. ( M 2 - 1 ) M ) K .Math. sin ( M 2 - 1 ) .Math. _ M 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. ( M 2 - 1 ) M ) K .Math. sin ( ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( 2 .Math. .Math. M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. 2 .Math. .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. ( M 2 - 1 ) M ) K .Math. sin ( 2 .Math. ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) .Math. 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( ( M - 2 ) .Math. M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. ( M - 2 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. ( M 2 - 1 ) M ) K .Math. sin ( ( M - 2 ) .Math. ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) 2 .Math. .Math. j ( 2 .Math. .Math. M ) K .Math. sin ( ( M - 1 ) .Math. M - _ M ) + 2 .Math. .Math. j ( 2 .Math. .Math. .Math. 2 M ) K .Math. sin ( 2 .Math. ( M - 1 ) .Math. M - 2 .Math. .Math. _ M ) + .Math. + 2 .Math. .Math. j ( 2 .Math. .Math. ( M 2 - 1 ) M ) K .Math. sin ( ( M - 1 ) .Math. ( M 2 - 1 ) .Math. M - ( M 2 - 1 ) .Math. _ M ) ] + [ ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j .Math. M 2 .Math. _ M ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. 2 .Math. .Math. M + M 2 .Math. _ M ) .Math. ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 2 ) .Math. M + M 2 .Math. _ M ) ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 1 ) .Math. M + M 2 .Math. _ M ) ] ( 5 )

    [0046] Then, .sub.n=1.sup.M n.sup.K sin n is expanded as follows:

    [00030] .Math. n = 1 M .Math. n K .Math. sin .Math. .Math. n .Math. .Math. = sin .Math. .Math. + 2 K .Math. sin .Math. .Math. 2 .Math. .Math. + 3 K .Math. sin .Math. .Math. 3 .Math. .Math. + 4 .Math. .Math. K .Math. .Math. sin .Math. .Math. 4 .Math. .Math. + .Math. + M K .Math. sin .Math. .Math. M .Math. .Math. = sin .Math. .Math. + sin .Math. .Math. 2 .Math. .Math. + sin .Math. .Math. 3 .Math. .Math. + sin .Math. .Math. 4 .Math. .Math. + .Math. + sin .Math. .Math. M .Math. .Math. + ( 2 K + 1 - 1 ) .Math. sin .Math. .Math. 2 .Math. .Math. + ( 3 K + 1 - 1 ) .Math. sin .Math. .Math. 3 .Math. .Math. + ( 4 K + 1 - 1 ) .Math. sin .Math. .Math. 4 .Math. .Math. + .Math. + ( M K + 1 - 1 ) .Math. sin .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. sin .Math. .Math. n .Math. .Math. ++ .Math. ( 2 K - 1 ) .Math. sin .Math. .Math. 2 .Math. .Math. + ( 2 K - 1 ) .Math. sin .Math. .Math. 3 .Math. .Math. + ( 2 K - 1 ) .Math. sin .Math. .Math. 4 .Math. .Math. + .Math. + ( 2 K - 1 ) .Math. sin .Math. .Math. M .Math. .Math. + [ ( 3 K - 1 ) - ( 2 K - 1 ) ] .Math. sin .Math. .Math. 3 .Math. .Math. + [ ( 4 K - 1 ) - ( 2 K - 1 ) ] .Math. sin .Math. .Math. 4 .Math. .Math. + .Math. + [ ( M K - 1 ) - ( 2 K - 1 ) ] .Math. sin .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. sin .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. sin .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. sin .Math. .Math. 3 .Math. .Math. + ( 3 K - 2 K ) .Math. sin .Math. .Math. 4 .Math. .Math. + ( 3 K - 2 K ) .Math. sin .Math. .Math. 4 .Math. .Math. + ( 3 K - 2 K ) .Math. sin .Math. .Math. 5 .Math. .Math. + .Math. + ( 3 K - 2 K ) .Math. sin .Math. .Math. M .Math. .Math. + [ ( 4 K - 2 K ) - ( 3 K - 2 K ) ] .Math. sin .Math. .Math. 4 .Math. .Math. + [ ( 5 K - 2 K ) - ( 3 K - 2 K ) ] .Math. sin .Math. .Math. 5 .Math. .Math. .Math. + .Math. + [ ( M K - 2 K ) - ( 3 K - 2 K ) ] .Math. sin .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. sin .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. sin .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. .Math. n = 1 M .Math. sin .Math. .Math. n .Math. .Math. + ( 4 K - 3 K ) .Math. sin .Math. .Math. 4 .Math. .Math. + ( 4 K - 3 K ) .Math. sin .Math. .Math. 5 .Math. .Math. + .Math. + ( 4 K - 3 K ) .Math. .Math. sin .Math. .Math. M .Math. .Math. + [ ( 5 K - 3 K ) - ( 4 K - 3 K ) ] .Math. sin .Math. .Math. 5 .Math. .Math. + [ ( 6 K - 3 K ) .Math. ( 4 K - 3 K ) ] .Math. sin .Math. .Math. 6 .Math. .Math. + .Math. .Math. + [ ( M K - 3 K ) - ( 4 K - 3 K ) ] .Math. sin .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. sin .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. sin .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. .Math. n = 3 M .Math. sin .Math. .Math. n .Math. .Math. + ( 4 K - 3 K ) .Math. .Math. n = 4 M .Math. sin .Math. .Math. n .Math. .Math. + [ ( 5 K - 3 K ) - ( 4 K - 3 K ) ] .Math. sin .Math. .Math. 5 .Math. .Math. + [ ( 6 K - 3 K ) .Math. ( 4 K - 3 K ) ] .Math. sin .Math. .Math. 6 .Math. + .Math. + [ .Math. ( M K - 3 K ) - ( 4 K - 3 K ) ] .Math. sin .Math. .Math. M .Math. .Math.

    Therefore,

    [0047] [00031] .Math. n = 1 M .Math. n K .Math. sin .Math. .Math. n .Math. .Math. = .Math. n = 1 M .Math. sin .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. sin .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. .Math. n = 3 M .Math. sin .Math. .Math. n .Math. .Math. + .Math. + ( M K - ( M - 1 ) K ) .Math. .Math. n = M - 1 M .Math. sin .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. ( n K - ( n - 1 ) K ) .Math. .Math. m = n M .Math. sin .Math. .Math. m .Math. .Math. .Math. and .Math. .Math. n = 1 M .Math. n K .Math. sin .Math. .Math. n .Math. .Math. = .Math. n = 1 M .Math. .Math. m = n M .Math. sin .Math. .Math. m .Math. .Math. ( n K - ( n - 1 ) K )

    Also,

    [0048] [00032] .Math. n = 1 N .Math. sin .Math. .Math. n .Math. .Math. = cos .Math. 2 - cos ( N + 1 2 .Math. ) 2 .Math. .Math. sin .Math. 2 .Math. n = N 1 N 2 .Math. sin .Math. .Math. n .Math. .Math. = cos ( N 1 - 1 2 ) .Math. - cos ( N 2 + 1 2 ) .Math. 2 .Math. .Math. sin .Math. 2 = sin .Math. N 2 + N 1 2 .Math. .Math. .Math. sin .Math. N 2 - N 1 + 1 2 .Math. sin .Math. 2 .Math. If .Math. .Math. m M = m .Math. .Math. M - _ M , then .Math. m = N 1 N 2 .Math. sin .Math. .Math. m ( m .Math. .Math. M - _ M ) = .Math. m = N 1 N 2 .Math. sin .Math. .Math. m .Math. .Math. m M = cos ( ( N 1 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( N 2 + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) = sin ( 1 2 .Math. ( N 1 + N 2 ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( N 2 - N 1 + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M )

    Especially,

    [0049] [00033] .Math. l = n M .Math. sin .Math. .Math. l ( m .Math. .Math. M - _ M ) = .Math. l = n M .Math. sin .Math. .Math. l .Math. .Math. m M = cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( M + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) = sin ( 1 2 .Math. ( M + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( M - n + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) .Math. n = 1 M .Math. n K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) = .Math. n = 1 M .Math. n K .Math. sin .Math. .Math. n .Math. .Math. m M = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. .Math. l = n M .Math. sin .Math. .Math. l .Math. .Math. m M ] = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( M + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( 1 2 .Math. ( M + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( M - n + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ]

    Then .sub.n=1.sup.M n.sup.K cos n is expanded as follows:

    [00034] .Math. n = 1 M .Math. n K .Math. cos .Math. .Math. n .Math. .Math. = cos .Math. .Math. + 2 K .Math. cos .Math. .Math. 2 .Math. .Math. + 3 K .Math. cos .Math. .Math. 3 .Math. .Math. + 4 K .Math. cos .Math. .Math. 4 .Math. .Math. + .Math. + M K .Math. cos .Math. .Math. M .Math. .Math. = cos .Math. .Math. + cos .Math. .Math. 2 .Math. .Math. + cos .Math. .Math. 3 .Math. .Math. + cos .Math. .Math. 4 .Math. .Math. + .Math. + cos .Math. .Math. M .Math. .Math. + ( 2 K + 1 - 1 ) .Math. cos .Math. .Math. 2 .Math. .Math. + ( 3 K + 1 - 1 ) .Math. cos .Math. .Math. 3 .Math. .Math. + ( 4 K + 1 - 1 ) .Math. cos .Math. .Math. 4 .Math. .Math. + .Math. + ( M K + 1 - 1 ) .Math. cos .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. cos .Math. .Math. n .Math. .Math. ++ .Math. ( 2 K - 1 ) .Math. cos .Math. .Math. 2 .Math. .Math. + ( 2 K - 1 ) .Math. cos .Math. .Math. 3 .Math. .Math. + ( 2 K - 1 ) .Math. cos .Math. .Math. 4 .Math. .Math. + .Math. + ( 2 K - 1 ) .Math. cos .Math. .Math. M .Math. .Math. + [ ( 3 K - 1 ) - ( 2 K - 1 ) ] .Math. cos .Math. .Math. 3 .Math. .Math. + [ ( 4 K - 1 ) - ( 2 K - 1 ) ] .Math. cos .Math. .Math. 4 .Math. .Math. + .Math. + [ ( M K - 1 ) - ( 2 K - 1 ) ] .Math. cos .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. cos .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. cos .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. cos .Math. .Math. 3 .Math. .Math. + ( 3 K - 2 K ) .Math. cos .Math. .Math. 4 .Math. .Math. + ( 3 K - 2 K ) .Math. cos .Math. .Math. 5 .Math. .Math. + .Math. + ( 3 K - 2 K ) .Math. cos .Math. .Math. M .Math. .Math. + [ ( 4 K - 2 K ) - ( 3 K - 2 K ) ] .Math. cos .Math. .Math. 4 .Math. .Math. + [ ( 5 K - 2 K ) - ( 3 K - 2 K ) ] .Math. cos .Math. .Math. 5 .Math. .Math. .Math. + .Math. + [ ( M K - 2 K ) - ( 3 K - 2 K ) ] .Math. cos .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. cos .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. cos .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. .Math. n = 3 M .Math. cos .Math. .Math. n .Math. .Math. + ( 4 K - 3 K ) .Math. cos .Math. .Math. 4 .Math. .Math. + ( 4 K - 3 K ) .Math. cos .Math. .Math. 5 .Math. .Math. + .Math. .Math. + ( 4 K - 3 K ) .Math. cos .Math. .Math. M .Math. .Math. + [ ( 5 K - 3 K ) - ( 4 K - 3 K ) ] .Math. cos .Math. .Math. 5 .Math. .Math. + [ ( 6 K - 3 K ) .Math. ( 4 K - 3 K ) ] .Math. cos .Math. .Math. 6 .Math. .Math. + .Math. .Math. + [ .Math. ( M K - 3 K ) - ( 4 K - 3 K ) ] .Math. cos .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. cos .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. cos .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. .Math. n = 3 M .Math. cos .Math. .Math. n .Math. .Math. + ( 4 K - 3 K ) .Math. .Math. n = 4 M .Math. cos .Math. .Math. n .Math. .Math. + [ ( 5 K - 3 K ) - ( 4 K - 3 K ) ] .Math. cos .Math. .Math. 5 .Math. .Math. + [ ( 6 K - 3 K ) .Math. ( 4 K - 3 K ) ] .Math. cos .Math. .Math. 6 .Math. .Math. + .Math. + [ ( M K - 3 K ) - ( 4 K - 3 K ) ] .Math. cos .Math. .Math. M .Math. .Math.

    Similarly,

    [0050] [00035] .Math. n = 1 M .Math. n K .Math. cos .Math. .Math. n .Math. .Math. = .Math. n = 1 M .Math. cos .Math. .Math. n .Math. .Math. + ( 2 K - 1 ) .Math. .Math. n = 2 M .Math. cos .Math. .Math. n .Math. .Math. + ( 3 K - 2 K ) .Math. .Math. n = 3 M .Math. cos .Math. .Math. n .Math. .Math. + .Math. + ( M K - ( M - 1 ) K ) .Math. .Math. n = M - 1 M .Math. cos .Math. .Math. M .Math. .Math. = .Math. n = 1 M .Math. ( n K - ( n - 1 ) K ) .Math. .Math. m = n M .Math. cos .Math. .Math. m .Math. .Math. .Math. and .Math. .Math. n = 1 M .Math. n K .Math. cos .Math. .Math. n .Math. .Math. = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. .Math. m = n M .Math. cos .Math. .Math. m .Math. .Math. ] .Math. Also , .Math. .Math. n = 1 N .Math. cos .Math. .Math. n .Math. .Math. = sin ( N + 1 2 ) .Math. 2 .Math. .Math. sin .Math. 2 - 1 2 .Math. n = N 1 N 2 .Math. cos .Math. .Math. n .Math. .Math. = sin ( N 2 + 1 2 ) .Math. - sin ( N 1 - 1 2 ) .Math. 2 .Math. .Math. sin .Math. 2 = cos .Math. ( N 2 + N 1 ) 2 .Math. .Math. .Math. sin .Math. N 2 - N 1 + 1 2 .Math. sin .Math. 2 .Math. If .Math. .Math. m M = m .Math. .Math. M - _ M , then .Math. m = n M .Math. cos .Math. .Math. m ( m .Math. .Math. M - _ M ) = .Math. m = n M .Math. cos .Math. .Math. m .Math. .Math. m M = sin ( M + n ) .Math. m M - sin ( n - 1 2 ) .Math. m M 2 .Math. .Math. sin .Math. 1 2 .Math. m M = cos .Math. ( M + n ) 2 .Math. m M .Math. sin .Math. M - n + 1 2 .Math. m M sin .Math. 1 2 .Math. m M .Math. .Math. .Math. m = n M .Math. cos .Math. .Math. m ( m .Math. .Math. M - _ M ) = sin ( M + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) = cos ( 1 2 .Math. ( M + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( M - n + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) . .Math. .Math. Especially , .Math. .Math. n = 1 M .Math. n k .Math. cos .Math. .Math. n ( m .Math. .Math. M - _ M ) = .Math. n = 1 M .Math. n k .Math. cos .Math. .Math. n .Math. .Math. m M = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. .Math. l = n M .Math. cos .Math. .Math. l .Math. .Math. m M ] = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] = .Math. n = 1 M .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( 1 2 .Math. ( M + n ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M - n + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] .

    [0051] For the design of an all-band digital differentiator of odd-numbered order (K is odd-numbered), if

    [00036] K - 1 2

    is an odd number, the frequency response function of the order-K differentiator is:

    [00037] H ( K ) ( w ) = { - jw K , 0 w < jw K , - w 0

    The frequency response function of M sampling points H(m) may be expressed as

    [00038] H ( m ) = H r ( 2 .Math. .Math. .Math. .Math. m M ) .Math. e j .Math. .Math. .Math. .Math. H ( m ) .Math. ,

    where

    [00039] H r ( 2 .Math. .Math. .Math. .Math. m M ) = { - ( 2 .Math. .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( 2 .Math. .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .

    The phase angle is:

    [00040] .Math. .Math. H ( m ) = { 2 - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - 2 + ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    assuming

    [00041] H m = { - j ( 2 .Math. .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. j ( 2 .Math. .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .Math. .Math. and .Math. .Math. .Math. .Math. H m = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1

    [0052] If M is an odd number, the following may be obtained from equation (1):

    [00042] .Math. W * .Math. H = [ - 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n .Math. .Math. _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( M - _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( 2 .Math. .Math. M - _ M ) .Math. 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 2 ) .Math. .Math. M - _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 1 ) .Math. .Math. M - _ M ) ] .Math. and .Math. h = [ h 0 h 1 .Math. h M - 1 ] T = 1 M .Math. W * .Math. H , .Math. .Math. where h m = .Math. 1 M [ 2 .Math. .Math. j .Math. .Math. .Math. n = 1 M - 1 2 .Math. ( - j ) .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M [ 2 .Math. ( 2 .Math. .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. n K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K - ( n - 1 ) K ) .Math. .Math. m = n M - 1 2 .Math. sin .Math. .Math. m ( m .Math. .Math. M - _ M ) ] } = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. .Math. M ) K .Math. .Math. .Math. n = 1 M - 1 2 .Math. .Math. [ ( n K - .Math. ( n - 1 ) K ) .Math. .Math. sin ( 1 2 .Math. ( M - 1 2 + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( M + 1 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } , .Math. 0 m M - 1. .Math.

    [0053] FIG. 1A is an amplitude response diagram for an all-band digital differentiator with M=51, K=3 according to the present invention. The digital differentiator resolves the problem raised by FIG. 1. As, when w=0 and w= (i.e., the normalized frequency is 1), the frequency response function of this type of order-K digital differentiator is 0. Therefore, all-band design would have significant error at high frequencies, as can be seen from FIG. 1A as well.

    [0054] The following Table 1 lists the digital signal processing system coefficients for M=51, order-3 differentiator with all-band design. If the digital signal processing system is embodied in a digital differentiator, Table 1 shows the digital differentiator coefficients.

    TABLE-US-00001 TABLE 1 Coefficient by the present Filter coefficient invention Filter coefficient h(0) 0.607551798225166 h(50) h(1) 0.609854469008660 h(49) h(2) 0.614503545413603 h(48) h(3) 0.621588534444675 h(47) h(4) 0.631248967997969 h(46) h(5) 0.643681252397096 h(45) h(6) 0.659148755404512 h(44) h(7) 0.677996129394913 h(43) h(8) 0.700669392011971 h(42) h(9) 0.727744079964486 h(41) h(10) 0.759965039064996 h(40) h(11) 0.798303432210820 h(39) h(12) 0.844039913190771 h(38) h(13) 0.898888709931315 h(37) h(14) 0.965187695786594 h(36) h(15) 1.046198692736652 h(35) h(16) 1.146599393008319 h(34) h(17) 1.273323945342434 h(33) h(18) 1.437072223661419 h(32) h(19) 1.655180760114309 h(31) h(20) 1.957452574498282 h(30) h(21) 2.398775578461656 h(29) h(22) 3.086430801945935 h(28) h(23) 4.197297752378359 h(27) h(24) 3.875843891007819 h(26) h(25) 0 h(25)

    [0055] FIG. 1B is an amplitude response diagram for an all-band digital differentiator with M=101, K=7 according to the present invention. As, when w=0 and w= (i.e., the normalized frequency is 1), the frequency response function of this type of order-K digital differentiator is 0. Therefore, all-bend design would have significant error at high frequencies, as can be seen from FIG. 1B as well.

    [0056] For even-numbered M, due to the anti-symmetric design, the following may be obtained from equation (2):

    [00043] W * .Math. H = [ .Math. - 2 .Math. .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n .Math. .Math. _ M + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j .Math. M 2 .Math. _ M 2 .Math. .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( M - _ M ) + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( 2 .Math. .Math. M - _ M ) + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. 2 .Math. .Math. M + M 2 .Math. _ M ) .Math. 2 .Math. .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 2 ) .Math. .Math. M + M 2 .Math. _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 1 ) .Math. .Math. M + M 2 .Math. _ M ) ] .Math. .Math. .Math. and .Math. .Math. .Math. .Math. .Math. h = [ h 0 h 1 .Math. h M - 1 ] T = 1 M .Math. W * .Math. H .Math. .Math. where .Math. .Math. h m = 1 M .Math. { 2 .Math. .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( - j ) .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) } = 1 M .Math. { 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. .Math. m + M - 1 2 .Math. } = 1 M .Math. { 2 .Math. ( 2 .Math. .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( 1 2 .Math. ( M 2 + n - 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. .Math. m + M - 1 2 .Math. } , .Math. .Math. 0 m M - 1. .Math.

    [0057] FIG. 2A is an amplitude response diagram for an all-band digital differentiator with M=50, K=3 according to the present invention. FIG. 2B is an amplitude response diagram for an all-bend digital differentiator with M=100, K=7 according to the present invention.

    [0058] For odd-numbered K and

    [00044] K - 1 2

    is an even number, the frequency response function of an order-K differentiator is:

    [00045] H ( K ) ( w ) = { jw K , 0 w < - jw K , - w 0

    The sampling frequency function H(m) is expressed as

    [00046] H ( m ) = H r ( 2 .Math. .Math. .Math. .Math. m M ) .Math. e j .Math. .Math. .Math. .Math. H ( m ) .Math. where H r ( 2 .Math. .Math. .Math. .Math. m M ) = { ( 2 .Math. .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - ( 2 .Math. .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    and the phase angle is

    [00047] .Math. .Math. H ( m ) = { 2 - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - 2 + ( M - 1 2 ) .Math. ( 2 .Math. .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    assuming that

    [00048] H m = { j ( 2 .Math. .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - j ( 2 .Math. .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .Math. .Math. and .Math. .Math. H m = { j ( 2 .Math. .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - j ( 2 .Math. .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .

    [0059] For odd-numbered M, the all-band design may be obtained from equation (1):

    [00049] W * .Math. H = [ .Math. - 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n .Math. .Math. _ M 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( M - _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( 2 .Math. .Math. M - _ M ) .Math. 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) 2 .Math. .Math. j .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) ] .Math. .Math. .Math. and .Math. .Math. .Math. h = [ h 0 h 1 .Math. h M - 1 ] T = 1 M .Math. W * .Math. H , where .Math. .Math. h m = .Math. 1 M [ 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M [ 2 .Math. ( 2 .Math. .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. n K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( 1 2 .Math. ( M - 1 2 + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( M + 1 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } , .Math. .Math. .Math. 0 m M - 1.

    [0060] FIG. 3A is an amplitude response diagram for an all-band digital differentiator with M=51, K=1, and FIG. 3B is an amplitude response diagram for an all-band digital diffeentiator with M=151, K=5.

    [0061] For even-numbered M, due to the anti-symmetric design, the following may be obtained from equation (2):

    [00050] .Math. W .Math. .Math. H = [ - 2 .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. sin .Math. .Math. n .Math. .Math. _ M + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. M 2 .Math. _ M 2 .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. sin .Math. .Math. n ( M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) 2 .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. sin .Math. .Math. n ( 2 .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. 2 .Math. M + M 2 .Math. _ M ) .Math. 2 .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. sin .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 2 ) .Math. M + M 2 .Math. _ M ) 2 .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. sin .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 1 ) .Math. M + M 2 .Math. _ M ) .Math. and .Math. .Math. h = [ h 0 .Math. .Math. h 1 .Math. .Math. .Math. .Math. .Math. h M - 1 ] T = 1 M .Math. W .Math. H , where h m = 1 M .Math. { 2 .Math. j .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. sin .Math. .Math. n ( m .Math. .Math. m - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) } = 1 M .Math. { 2 .Math. j ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. n K .Math. .Math. sin .Math. .Math. n ( m .Math. .Math. m - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. } = 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. n K .Math. .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. ] = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( 1 2 .Math. ( M 2 + n - 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + .Math. ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. } , .Math. 0 m M - 1.

    [0062] FIG. 4A is an amplitude response diagram for an all-band digital differentiator with M=50, K=1, and FIG. 4B is an amplitude response diagram for a all-band digital differentiator with M=80, K=5.

    [0063] For all-band digital differentiator with even-numbered order where K is even-numbered and

    [00051] K 2

    is an odd number, the frequency response function is:

    [00052] H K ( w ) = { - w K , 0 w < - w K , - w 0

    and the sampling frequency response function H(m) is expressed as

    [00053] H ( m ) = H r ( 2 .Math. .Math. .Math. m M ) .Math. e j .Math. .Math. .Math. .Math. H ( m ) .Math. .Math. where .Math. .Math. H r ( 2 .Math. .Math. .Math. m M ) = { - ( 2 .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - ( 2 .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    and the phase angle is

    [00054] .Math. .Math. H ( m ) = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( M - 1 2 ) .Math. ( 2 .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    assuming that

    [00055] H m = { - ( 2 .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - ( 2 .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 , .Math. and .Math. .Math. .Math. H m = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( ( M - 1 ) M ) .Math. ( 2 .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .

    [0064] For odd-numbered M, due to the symmetric design, the following may be obtained from equation (1):

    [00056] .Math. W .Math. .Math. H = [ 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n .Math. .Math. _ M 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( M - _ M ) 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( 2 .Math. M - _ M ) .Math. 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) 2 .Math. .Math. k = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. k ( ( M - 1 ) .Math. M - _ M ) ] .Math. and .Math. h = [ h 0 .Math. .Math. h 1 .Math. .Math. .Math. .Math. .Math. h M - 1 ] T = 1 M .Math. W .Math. H .Math. where h m = 1 M [ 2 .Math. .Math. k = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. k ( m .Math. .Math. m - _ M ) ] = 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. k = 1 M - 1 2 .Math. n K .Math. .Math. cos .Math. .Math. k ( m .Math. .Math. m - _ M ) ] = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K .Math. - ( n - 1 ) K ) .Math. cos .Math. .Math. ( 1 2 .Math. ( M - 1 2 + n ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M + 1 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] .Math. } , .Math. 0 m M - 1.

    [0065] FIG. 5A is an amplitude response diagram for an all-band digital differntiator with M=51, K=2, and FIG. 5B is an amplitude response diagram for an all-band digital differentiator with M=85, K=6.

    [0066] For even-numbered M, due to the symmetric design, the following may be obtained from equation (2):

    [00057] .Math. W .Math. .Math. H = [ 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n .Math. .Math. _ M + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. M 2 .Math. _ M 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( 2 .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. 2 .Math. M + M 2 .Math. _ M ) .Math. 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 2 ) .Math. M + M 2 .Math. _ M ) 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 1 ) .Math. M + M 2 .Math. _ M ) ] .Math. and .Math. .Math. h = [ h 0 .Math. .Math. h 1 .Math. .Math. .Math. .Math. .Math. h M - 1 ] T = 1 M .Math. W .Math. H , where h m = 1 M [ 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( m .Math. .Math. m - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) ] = 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. n K .Math. .Math. cos .Math. .Math. n ( m .Math. .Math. m - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m + M - 1 2 ) .Math. ] = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K .Math. - ( n - 1 ) K ) .Math. cos .Math. .Math. ( 1 2 .Math. ( M 2 - 1 + n ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] .Math. + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) }

    [0067] FIG. 6A is an amplitude response diagram for an all-band digital differentiator with M=50, K=2, and FIG. 6B is an amplitude response diagram for an all-band digital differentiator with M=150, K=6.

    [0068] For even-numbered K and

    [00058] K 2

    is an even number, the frequency response function of an order-K differentiator is:

    [00059] H ( K ) ( w ) = { jw K , 0 w < - jw K , - w 0 ,

    the sampling frequency response function H(m) is expressed as

    [00060] H ( m ) = H r ( 2 .Math. .Math. .Math. m M ) .Math. e jH ( m ) , .Math. where .Math. H r ( 2 .Math. .Math. .Math. m M ) = { ( 2 .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - ( 2 .Math. ( M - m ) M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    the phase angle is

    [00061] .Math. H ( m ) = { 2 - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - 2 + ( M - 1 M ) .Math. ( 2 .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 ,

    assuming that

    [00062] H m = { j ( 2 .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. - j ( 2 .Math. M - m M ) K , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 , .Math. and .Math. .Math. H m = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. ( M - 1 M ) .Math. ( 2 .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .

    [0069] For odd-numbered M, due to the symmetric design, the following may be obtained from equation (1):

    [00063] .Math. W .Math. .Math. H = [ 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n .Math. .Math. _ M ) 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( M - _ M ) 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( 2 .Math. M - _ M ) .Math. 2 .Math. .Math. n = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) 2 .Math. .Math. k = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. k ( ( M - 1 ) .Math. M - _ M ) ] .Math. and .Math. h = [ h 0 .Math. .Math. h 1 .Math. .Math. .Math. .Math. h M - 1 ] T = 1 M .Math. W .Math. H .Math. where h m = 1 M [ 2 .Math. .Math. k = 1 M - 1 2 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. k ( m .Math. .Math. m - _ M ) ] = 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. k = 1 M - 1 2 .Math. n K .Math. .Math. cos .Math. .Math. k ( m .Math. .Math. m - _ M ) ] = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K .Math. - ( n - 1 ) K ) .Math. cos .Math. .Math. ( 1 2 .Math. ( M - 1 2 + n ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M + 1 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] .Math. } , 0 m M - 1.

    [0070] FIG. 7A is an amplitude response diagram for an all-band digital differentiator with M=51, K=2, and FIG. 7B is an amplitude response diagram for an all-band digital differentiator with M=201, K=6.

    [0071] For even-numbered M, due to the anti-symmetric design, the following may be obtained from equation (2):

    [00064] .Math. W .Math. .Math. H = [ 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n .Math. .Math. _ M + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. M 2 .Math. _ M 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. M + M 2 .Math. _ M ) 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( 2 .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. 2 .Math. M + M 2 .Math. _ M ) .Math. 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 2 ) .Math. M + M 2 .Math. _ M ) 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( M 2 .Math. ( M - 1 ) .Math. M + M 2 .Math. _ M ) ] .Math. and .Math. .Math. h = [ h 0 .Math. .Math. h 1 .Math. .Math. .Math. .Math. h M - 1 } T = 1 M .Math. W .Math. H , where h m = 1 M [ 2 .Math. .Math. n = 1 M 2 - 1 .Math. ( 2 .Math. .Math. n M ) K .Math. .Math. cos .Math. .Math. n ( m .Math. .Math. m - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) ] = 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. n K .Math. .Math. cos .Math. .Math. n ( m .Math. .Math. m - _ M ) + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m + M - 1 2 ) .Math. ] = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K .Math. - ( n - 1 ) K ) .Math. cos .Math. .Math. ( 1 2 .Math. ( M 2 - 1 + n ) .Math. ( m .Math. .Math. M - _ M ) ) sin ( 1 2 .Math. ( M 2 - n ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] .Math. + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) }

    [0072] FIG. 8A is an amplitude response diagram for an all-band digital diffeentiator with M=50, K=2, and FIG. 8B is an amplitude response diagram for an all-band digital differentiator with M=200, K=6.

    [0073] The following is a design method for a partial-bend digital differentiator. For partial-band design, it is assumed that

    [00065] N 1 < .Math. M - 1 2 .Math. .

    The design specification is:

    [00066] H m = { - j ( 2 .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , N 1 0 , j ( 2 .Math. ( M - m ) M ) K , m = N 1 + 1 , .Math. .Math. , M - 1 , and .Math. .Math. H m = { - ( M - 1 2 ) .Math. ( 2 .Math. .Math. .Math. m M ) , m = 0 , 1 , .Math. .Math. , .Math. M - 1 2 .Math. 0 , ( M - 1 2 ) .Math. ( 2 .Math. ( M - m ) M ) , m = .Math. M - 1 2 .Math. + 1 , .Math. .Math. , M - 1 .

    [0074] For the design of a partial-band differentiator with an odd-numbered order, the following may be obtained from equation (1) when K is an odd number:

    [00067] .Math. W * .Math. H = [ - 2 .Math. j .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. sin .Math. .Math. n .Math. .Math. _ M ) 2 .Math. j .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( M - _ M ) 2 .Math. j .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( 2 .Math. M - _ M ) .Math. 2 .Math. j .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) 2 .Math. j .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) ] .Math. and .Math. .Math. h = [ h 0 h 1 .Math. h M - 1 ] T = 1 M .Math. W * .Math. H . .Math. Therefore , .Math. h m = .Math. 1 M [ 2 .Math. j .Math. .Math. n = 1 N 1 .Math. .Math. ( - j ) .Math. ( 2 .Math. .Math. n M ) K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 N 1 .Math. .Math. n K .Math. sin .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 N 1 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. .Math. m = n N 1 .Math. .Math. sin .Math. .Math. m ( m .Math. .Math. M - _ M ) ] } = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. .Math. n = 1 N 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( 1 2 .Math. ( N 1 + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( N 1 - n + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } , .Math. 0 m M - 1.

    [0075] FIG. 9 is an amplitude response diagram for a partial-band digital differentiator with M=201, N.sub.1=81, K=3, and FIG. 10 is an amplitude response diagram for a partial-band digital differentiator with M=451, N.sub.1=91, K=7.

    [0076] FIG. 11 is an amplitude response diagram for a partial-band digital differentiator with M=100, =N.sub.1 35, K=3, FIG. 12 is an amplitude response diagram for a partial-band digital differentiator with M=151, N.sub.1=57, K=1, and FIG. 13 is an amplitude response diagram for a partial-band digital differentiator with M=100, N.sub.1=30, K=5.

    [0077] For the design of a partial-band differentiator with an even-numbered order, the following may be obtained from equation (2) when K is an even number.

    [00068] W * .Math. H = [ 2 .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. cos .Math. .Math. n .Math. .Math. _ M ) 2 .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. cos .Math. .Math. n ( M - _ M ) 2 .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. cos .Math. .Math. n ( 2 .Math. M - _ M ) .Math. 2 .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. cos .Math. .Math. n ( ( M - 2 ) .Math. M - _ M ) 2 .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. cos .Math. .Math. n ( ( M - 1 ) .Math. M - _ M ) ] and .Math. .Math. h = [ h 0 h 1 .Math. h M - 1 ] T = 1 M .Math. W * .Math. H , therefore h m = .Math. 1 M [ 2 .Math. .Math. n = 1 N 1 .Math. .Math. ( 2 .Math. .Math. n M ) K .Math. cos .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M [ 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 N 1 .Math. .Math. n K .Math. cos .Math. .Math. n ( m .Math. .Math. M - _ M ) ] = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 N 1 .Math. .Math. [ ( n K - ( n - 1 ) K ) .Math. .Math. m = n N 1 .Math. .Math. cos .Math. .Math. m ( m .Math. .Math. M - _ M ) ] } = .Math. 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. .Math. n = 1 N 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( 1 2 .Math. ( N 1 + n ) .Math. ( m .Math. .Math. M - _ M ) ) .Math. sin ( 1 2 .Math. ( N 1 - n + 1 ) .Math. ( m .Math. .Math. M - _ M ) ) sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } , .Math. 0 m M - 1.

    [0078] FIG. 14 is an amplitude response diagram for a partial-band digital differentiator with M=101, N.sub.1=41, K=2, FIG. 15 is an amplitude response diagram for a partial-band digital differentiator with M=200, N.sub.1=60, K=6, FIG. 16 is an amplitude response diagram for a partial-band digital differentiator with M=201, N.sub.1=86, K=4, and FIG. 7 is an amplitude response diagram for a partial-band digital differentiator with M=100, N.sub.1=48, K=8.

    [0079] The following table is the frequency response functions and frequency response coefficients for all-band design according to the present invention.

    [0080] The all-band design reference table:

    TABLE-US-00002 Ideal freq. response Design method H.sup.(K)(w) h.sub.m, 0 m M 1 [00069] K .Math. .Math. is .Math. .Math. odd .Math. - .Math. numbered , K - 1 2 .Math. .Math. is .Math. .Math. odd .Math. - .Math. numbered M is odd- num- bered [00070] { - jw K , 0 w < jw K , - w 0 [00071] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( M 2 .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } M is even- num- bered [00072] { - jw K , 0 w < jw K , - w 0 [00073] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. } [00074] K .Math. .Math. is .Math. .Math. odd .Math. - .Math. numbered , K - 1 2 .Math. .Math. is .Math. .Math. even .Math. - .Math. numbered M is odd- num- bered [00075] { jw K , 0 w < - jw K , - w 0 [00076] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( M 2 .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } M is even- num- bered [00077] { jw K , 0 w < - jw K , - w 0 [00078] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j .Math. 2 .Math. m + M - 1 2 .Math. } [00079] K .Math. .Math. is .Math. .Math. even .Math. - .Math. numbered , K 2 .Math. .Math. is .Math. .Math. odd .Math. - .Math. numbered M is odd- num- bered w.sup.K, w < [00080] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } M is even- num- bered w.sup.K, w < [00081] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) } [00082] K .Math. .Math. is .Math. .Math. even .Math. - .Math. numbered , K 2 .Math. .Math. is .Math. .Math. even .Math. - .Math. numbered M is odd- num- bered [00083] { jw K , 0 w < - jw K , - w 0 [00084] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M - 1 2 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } M is even- num- bered [00085] { jw K , 0 w < - jw K , - w 0 [00086] 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 M 2 - 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( M 2 - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) - sin ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] + ( 2 .Math. .Math. M 2 M ) K .Math. e j ( m .Math. M 2 .Math. M + M 2 .Math. _ M ) }

    [0081] The partial-band design reference table, assuming

    [00087] N 1 < .Math. M - 1 2 .Math. :

    TABLE-US-00003 [00088] Frequency .Math. .Math. response H m .Math. { - j ( 2 .Math. .Math. .Math. m M ) K , m = 0 , 1 , .Math. .Math. , N 1 0 , j ( 2 .Math. ( M - m ) M ) K , m = N 1 + 1 , .Math. .Math. , M - 1 Design method h.sub.m, 0 m M 1 [00089] K .Math. .Math. is .Math. .Math. odd .Math. - .Math. numbered , h m = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 N 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. cos ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - cos ( ( N 1 + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] } [00090] K .Math. .Math. is .Math. .Math. even .Math. - .Math. numbered , h m = 1 M .Math. { 2 .Math. ( 2 .Math. M ) K .Math. .Math. n = 1 N 1 .Math. [ ( n K - ( n - 1 ) K ) .Math. sin ( ( N 1 + 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) - sin ( ( n - 1 2 ) .Math. ( m .Math. .Math. M - _ M ) ) 2 .Math. sin .Math. 1 2 .Math. ( m .Math. .Math. M - _ M ) ] }

    [0082] The present invention also provides a digital differentiator where its frequency response coefficient h.sub.m is determined by the above tables. Then the type of the digital differentiator such as direct series-connected, or linear-phased digital differentiator may be determined. The setup of the digital differentiator is as such completed.

    [0083] For example, for an order-N FIR filter having impulse response h[n], its output signal is y[n]=h[n]*x[n] where * is convolution operation. This output signal can also be expressed as


    y[n]=h[n]*x[n]=.sub.k=0.sup.Nh[k]x[nk](eq 1)

    [0084] (eq 1) means the output signal y[n] is the difference equation between the system impulse response h[n] and the input signal x[n]. FIG. 19 is an operational block diagram of a direct type FIR differentiator 1, which has an input terminal 11 and an output terminal 12. z.sup.1 is a delay element 13 that delays the input signal x[n] for an unit of time and outputs signal x[n1]. The arrow lines indicate the propagation of signals. The h[n] above each arrow line is a multiplier 14 providing multiplication operation, and the outputs from the multipliers 14 are summed by adders 15 as the output signal y[n]. Under the given specification, the h[n] of various orders of the digital differentiator 1 is designed according to the previous tables. After h[n] is determined, the amplification factor of the multiplier 14 of the corresponding order is determined, thereby completing the design of the digital differentiator 1.

    [0085] While certain novel features of this invention have been shown and described and are pointed out in the annexed claim, it is not intended to be limited to the details above, since it will be understood that various omissions, modifications, substitutions and changes in the forms and details of the device illustrated and in its operation can be made by those skilled in the art without departing in any way from the claims of the present invention.