Maximally Flat Lowpass Digital Differentiators
|
|
- Linette Cannon
- 5 years ago
- Views:
Transcription
1 Maximally Flat Lowpass Digital Differentiators Ivan W. Selesnick August 3, 00 Electrical Engineering, Polytechnic University 6 Metrotech Center, Brooklyn, NY 0 selesi@taco.poly.edu tel: fax: Abstract This paper describes the design of Type III and Type IV linear-phase FIR lowpass digital differentiators according to the maximally flat criterion. We introduce a two-term recursive formula that enables the simple stable computation of the impulse response coefficients. The same recursive formula is valid for both Type III and Type IV solutions. Revised manuscript submitted to IEEE Transactions on Circuits and Systems II: Analog and Digital Signal Processing. Original manuscript submitted February 6, 00. Manuscript received This work was supported in part by the NSF under CAREER grant CCR-9875.
2 List of Figures Type IV maximally flat lowpass differentiators Type III maximally flat lowpass differentiators List of Tables Numerical values of the weights c(n) for 0 K 3, 0 n
3 Introduction This paper describes the design of linear-phase FIR lowpass digital differentiators according to the maximally flat criterion. We introduce a two-term recursive formula that enables the simple stable computation of the impulse response coefficients. The same recursive formula is valid for both Type III and Type IV solutions. To avoid the undesirable amplification of noise in digital differentiation, lowpass differentiators can be used in place of full-band ones. They find application, for example, in bar code readers [6]. Maximally flat digital full-band differentiators have been described by several authors [3, 7, 8, 9, ] where the approximation is accurate at a single point, usually ω = 0, π/, or π. When the input signal lies in a known frequency band centered around the frequency ω o, it is desirable to use a differentiator which is maximally flat at ω o rather than at ω = 0. Solutions of that type have been introduced in [, 0,, ]. The design of lowpass differentiators can be performed with least squares, the Remez algorithm [8] and other methods [3, 3, ]; however, the design of maximally flat lowpass differentiators has not been previously addressed. The frequency response of the ideal (full-band) digital differentiator is H F B (e jω ) = j ω, ω < π. () The maximally flat FIR approximation to the ideal differentiator satisfies the derivative constraints, H(e jω ) = 0, ω = 0 () d dω H(ejω ) =, ω = 0 (3) d k dω k H(ejω ) = 0, ω = 0, k B () where B is the number of free parameters. Notice that the approximation is performed at a single frequency point ω = 0. Kumar and Dutta Roy described how the solution can be obtained from the maximally flat lowpass FIR filter [9]. The frequency response of the ideal lowpass digital differentiator is H LP (e jω j ω ω < ω c ) = 0 ω c < ω < π where ω c is the cutoff frequency. This paper describes the maximally flat linear-phase FIR approx- (5) 3
4 imation, which satisfies the constraints H(e jω ) = 0, ω = 0 (6) d dω H(ejω ) =, ω = 0 (7) d k dω k H(ejω ) = 0, ω = 0, k L (8) d k dω k H(ejω ) = 0, ω = π, 0 k M. (9) In this case, it is not a single point approximation problem and the solution can not be obtained directly from the maximally flat lowpass filter. Type III and Type IV linear-phase FIR solutions to this problem are derived in Section and Section 3 respectively. As in [], the method used in this paper is based on power series expansion. To compute the impulse response coefficients in a simple way, Section gives a two-term recursive formula which is valid for both Type III and IV solutions. This recursive formula was obtained using the algorithms of [5] for automatic hypergeometric-type sum simplification. A Matlab program to compute the maximally flat differentiator based on the two-term recurrence is available from the author. The derivation of the maximally flat lowpass differentiator below uses the transformation of variables used, for example, by Herrmann in the derivation of maximally flat lowpass Type I FIR filters [6]. Similar to [], the expression obtained in (58) below for the transfer function can be used as the basis of a structure for an efficient implementation of the filter. Type IV Solution The derivation of the solution will depend on a transformation that maps polynomials on the real interval [0, ] to polynomials on the upper half of the unit circle. Let P (x) be a polynomial, P (x) = p 0 + p x + p x + + p n x n. (0) If we define ( ) z + z H(z) = P, () then H(z) will be a Type I transfer function and its frequency response is given by ( H(e jω ) = P ) cos(ω). ()
5 Now, note that a type IV transfer function H (z) can always be written as H (z) = ( z ) H(z). (3) (This is because a Type IV transfer function always has a zero at z =. See [3, 7], for example, for a description of the four types of FIR linear-phase digital filter.) The frequency response of H (z) is given by where Note that even though the term e j ω H (e jω ) = ( e jω ) H(e jω ) () ( = j e j ω ω ) sin H(e jω ) (5) = j e j ω ( cos(ω) P ) cos(ω) (6) ( = j e j ω P ) cos(ω) (7) P (x) := x P (x). (8) appears in the expression for H (e jω ) no fractional delay element is needed for the implementation of this digital filter. (H(z) is a Type I transfer function which requires no fractional delays, and (3) does not introduce any fractional delay.) Our goal is to design P (x) so that H (e jω ) approximates the ideal lowpass differentiator. That means it is desired that ( P ) cos(ω) ω for ω 0 0 for ω π. (9) Equivalently, or in terms of P (x), arccos( x) for x 0 P (x) 0 for x arccos( x) for x 0 P (x) x 0 for x. (0) () These approximation condition can be written as P (x) = A(x) ( x) M () 5
6 and arccos( x) x P (x) = B(x) x L+ (3) where M, and L + are the degrees of tangency at x = and x = 0. The parameters M and L are the parameters that define the solution. A(x) and B(x) are power series in x. Adding () and (3), we get Solving for A(x) we get, arccos( x) x = B(x) x L+ + A(x) ( x) M. () A(x) = Let the Taylor series of the first term be given by: ( ) arccos( x) ( x) M B(x) x L+. (5) x arccos( x) = ( x) M x c(n) x n. (6) Note that the polynomial A(x) of minimal degree which satisfies (5) will be of degree L. (Suppose the degree of A(x) were less than L, then it would be impossible to choose a power series B(x) so that the right-hand side of (5) has the same degree as the left-hand side.) With A(x) being minimal degree, we have A(x) = L c(n) x n, (7) and P (x) = ( x) M L c(n) x n, (8) H (z) = ( z ) H(z) (9) = ( ) z + z ( z ) P (30) = ( ) z + + z M ( L ( ) z + z n z ) c(n). (3) Note that () the degree of P (x) is M + L, () the degree of H (z) is M + L +. That makes h (n) an impulse response of length M + L +. We need only find the Taylor series coefficients c(n). To this end, we define A (x) := ( x) M = c (n) x n (3) 6
7 where and where A (x) := ( ) M + n c (n) = n c (n) = arccos( x) x = n + n! (33) c (n) x n (3) ( ) n which is derived in the Appendix. Then A(x) = A (x) A (x) and c(n) = (35) n c (k) c (n k). (36) In Section, we give a recursive formula for computing the coefficients c(n). k=0 3 Type III Solution Note that a type III transfer function H 3 (z) can always be written as H 3 (z) = ( z ) ( + z ) H(z) (37) where H(z) is a Type I transfer function. (This is because a Type III transfer function always has a zero at z = and z =.) The frequency response of H 3 (z) is given by ( H 3 (e jω ) = j e j ω ω ) ( sin e j ω ω ) cos H(e jω ) (38) = j e jω cos(ω) + ( cos(ω) P ) cos(ω) (39) ( = j e j ω P3 ) cos(ω) (0) where P 3 (x) := x x P (x). () To design P (x) so that H 3 (e jω ) approximates the ideal lowpass differentiator, we ask that ( P 3 ) cos(ω) ω for ω 0 0 for ω π. () 7
8 Equivalently, or in terms of P (x), arccos( x) for x 0 P 3 (x) 0 for x arccos( x) for x 0 P (x) x x 0 for x. (3) () These approximation condition can be written as and P (x) = A(x) ( x) M (5) arccos( x) x x P (x) = B(x) x L+ (6) where M, and L + are the degrees of tangency at x = and x = 0. Adding (5) and (6), we get arccos( x) x x = B(x) x L+ + A(x) ( x) M. (7) As above, the polynomial A(x) of minimal degree which satisfies (7) will be of degree L. Solving for A(x) of minimal degree L we get A(x) = L c(n) x n (8) where c(n) are the coefficients in the Taylor series of Then arccos( x) x ( x) M+ = c(n) x n. (9) H 3 (z) = ( z ) ( + z ) H(z) (50) = ( ) z + z ( z ) ( + z ) P (5) = ( ) z + + z ( z ) ( + z M L ( ) z + z n ) c(n). (5) 8
9 Note that () P (x) is of degree M + L again, and () H 3 (z) is of degree M + L +. That makes h 3 (n) an impulse response of length M + L + 3. To determine c(n), define A (x) as A (x) := ( x) M+ = c (n) x n (53) where ( M + c (n) = + n ) n (5) and define A (x) as in (3). Then A(x) = A (x) A (x) and c(n) = n k=0 c (k) c (n k). The coefficients c (n) are the same for Type III and Type IV solutions, while c (n) has a different form. To evaluate the binomial coefficient for fractional values of the upper entry, we can use the Gamma function Γ, ( M + c (n) = + n ) n = Γ(M + + n) Γ(M + (55) ) Γ(n + ). The Gamma function also satisfies Γ(α + ) = α Γ(α), so we can write a recursive equation for evaluating c (n), In Section we unify the Type III and Type IV solutions. c (n) = n (M + n /) c (n ). (56) Combined Formula Let K denote the number of zeros a transfer function has at z =. A Type IV transfer function always has an even number of zeros at z = ; for the solution given in Section 3 we have K = M. On the other hand, a Type III transfer function always has an odd number of zeros at z = ; for the solution given in Section we have K = M +. In terms of K, the coefficients c (n) in (33) and (5) are given by the same formula, namely, ( ) K/ + n c (n) =. (57) n A causal maximally flat lowpass differentiator with linear phase is given by ( ) ( ) z + z K L H(z) = z L c(n) ( ) z + z n (58) 9
10 Table : Numerical values of the weights c(n) for 0 K 3, 0 n 0. c(n) n K = 0 K = K = K = with n c(n) = c (k) c (n k) (59) k=0 where c (n) is given by (57) and c (n) is given by (35). Using the algorithms given in [5] we can find the following two-term recursive formula for computing the coefficients c(n). c(n) = (8n + Kn 0n K + 3) c(n ) (n + K 3) c(n ) n (n + ) (60) for n with c(0) =, and c() = K + /3. The formula in (58) with the recursion for c(n) in (60) generates both Type III and Type IV solutions. When K is even we obtain a Type IV transfer function, when K is odd we obtain a Type III transfer function. In either case, the length of the impulse response is N = K + L +. Table gives numerical values for c(n) for several values of K and n. Notice that c(n) does not depend on L; rather L determines how many values of c(n) are needed in (58). As noted in [], this means that extra blocks can be added in a modular fashion to improve the frequency response. See [] for a discussion of implementation issues and block diagrams. 0
11 N = 30; K = 0::, L = (N K)/ 0.8 K=0 H(e jω ) /π K= K= ω/π Figure : Type IV maximally flat lowpass differentiators. To illustrate the maximally flat lowpass differentiator, Figure shows the frequency response for a family of Type IV lowpass differentiators of length N = 30, where K is varied from 0 to in increments of, and where L = (N K)/. The bandwidth of the frequency response depends on the relative values of K and L. When K = 0 we obtain a full band differentiator. Similarly Figure shows the frequency response for a family of Type III lowpass differentiators of length N = 3, where K is varied from to 5 in increments of, and where L = (N K)/. 5 Conclusion This paper describes the design of lowpass linear-phase FIR digital differentiators according to the maximally flat criterion. The solutions can not be obtained from a lowpass filter as in the case of a full-band differentiator. The algorithms for automatic sum simplification described in [5] were used to obtain a simple two-term recurrence relation for computing the coefficients of the impulse response. The equations (58) and (60) contain the main result. There are several possible extensions to the problem described in this paper. For example, the extension of the recursive formulas to the case where the maximally flat approximation to the ideal differentiator is performed not at ω = 0 but at another frequency ω o. For the full-band
12 N = 3; K = ::5, L = (N K)/ 0.8 H(e jω ) /π K= K= ω/π Figure : Type III maximally flat lowpass differentiators. differentiator, solutions are given in [,, ]. This type of solution is relevant when the signal is centered around a known frequency (as in radar using Doppler tracking []). Another remaining question is the existence of low-complexity structures for maximally flat differentiators, of the kind described in [, ] by Samadi, Nishihara, and Iwakura for maximally flat lowpass filters. Those structures are multiplierless and have a regular structure. Another extension of the approach described in this paper is the design of second (and higher) order differentiators as in [0] where the desired frequency response is ω, etc. For applications where a variable fractional sample delay is required, an extension along the lines of [5] will be of interest.
13 A Derivation of c (n) To find c (n) in (3), note that arccos( x) A (x) x x = x x A (x) = x (6) (6) so A (x) satisfies the differential equation A (x) x + A (x) = Substituting the Taylor series for A (x) into this differential equation gives Now note that where c (n) (n + ) xn = x = x. (63) x. (6) g(n) x n (65) g(n) = n! ( ) n and where (a) n is the Pochhammer function, or rising factorial, defined as (66) Then we can write (a) n = (a) (a + ) (a + n ). (67) }{{} n terms c (n) (n + ) xn = g(n) x n. (68) Therefore, matching like powers in (68), we get c (n) = g(n) n + / (69) or c (n) = n + n! ( ). (70) n 3
14 References [] J. Le Bihan. Maximally linear FIR digital differentiators. J. Circuits, Syst., Signal Processing, (5): , 995. [] J. Le Bihan. Coefficients of FIR digital differentiators and Hilbert transformers for midband frequencies. IEEE Trans. on Circuits and Systems II, 3(3):7 7, March 996. [3] B. Carlsson. Maximum flat digital differentiator. Electron. Lett., 7: , April 99. [] IEEE DSP Committee, editor. Selected Papers In Digital Signal Processing,II. IEEE Press, 976. [5] E. Hermanowicz and M. Rojewski. Design of FIR first order digital differentiators of variable fractional sample delay using maximally flat error criterion. Electron. Lett., 30():7 8, 6 January 99. [6] O. Herrmann. On the approximation problem in nonrecursive digital filter design. IEEE Trans. on Circuit Theory, 8(3): 3, May 97. Also in [9]. [7] I. R. Khan and R. Ohba. New design of full band differentiators based on Taylor series. IEE Proc.-Vis. Image Signal Process, 6():85 89, August 999. [8] B. Kumar and S. C. Dutta Roy. Coefficients of maximally linear, FIR digital differentiators for low frequencies. Electron. Lett., (9): , 8 April 988. [9] B. Kumar and S. C. Dutta Roy. Design of digital differentiators for low frequencies. Proc. IEEE, 76(3):87 89, March 988. [0] B. Kumar and S. C. Dutta Roy. Design of efficient FIR digital differentiators and Hilbert transformers for midband frequency ranges. Int. J. Circuit Theory Appl., 7:83 88, 989. [] B. Kumar and S. C. Dutta Roy. Maximally linear FIR digital differentiators for midband frequencies. Int. J. Circuit Theory Appl., 7: 7, 989. [] B. Kumar, S. C. Dutta Roy, and H. Shah. On the design of FIR digital differentiators which are maximally linear at the frequency π/p, p {positive integers}. IEEE Trans. on Signal Processing, 0(9):33 338, September 99. [3] T. W. Parks and C. S. Burrus. Digital Filter Design. John Wiley and Sons, 987.
15 [] S.-C. Pei and P.-H. Wang. Closed-form design of maximally flat FIR Hilbert transformers, differentiators, and fractional delayers by power series expansion. IEEE Trans. on Circuits and Systems I, 8(): , April 00. [5] M. Petkovšek, H. S. Wilf, and D. Zeilberger. A=B. A K Peters, 996. [6] W. Preuss. Symbol Technologies, NY. Personal communication, 000. [7] L. R. Rabiner and B. Gold. Theory and Application of Digital Signal Processing. Prentice Hall, Englewood Cliffs, NJ, 975. [8] L. R. Rabiner, J. H. McClellan, and T. W. Parks. FIR digital filter design techniques using weighted Chebyshev approximation. Proc. IEEE, 63():595 60, April 975. Also in []. [9] L. R. Rabiner and C. M. Rader, editors. Digital Signal Processing. IEEE Press, 97. [0] M. R. R. Reddy, S. C. Dutta Roy, and B. Kumar. Design of efficient second and higher degree FIR digital differentiators for midband frequencies. IEE Proceedings, Part G, 38():9 33, February 99. [] S. Samadi, T. Cooklev, A. Nishihara, and N. Fujii. Multiplierless structure for maximally flat linear phase FIR filters. Electron. Lett., 9():8 85, Jan 993. [] S. Samadi, A. Nishihara, and H. Iwakura. Universal maximally flat lowpass FIR systems. IEEE Trans. on Signal Processing, 8(7):956 96, July 000. [3] H. W. Schüssler and P. Steffen. Some advanced topics in filter design. In J. S. Lim and A. V. Oppenheim, editors, Advanced Topics in Signal Processing, chapter 8, pages 6 9. Prentice Hall, 988. [] I. W. Selesnick and C. S. Burrus. Exchange algorithms for the design of linear phase FIR filters and differentiators having flat monotonic passbands and equiripple stopbands. IEEE Trans. on Circuits and Systems II, 3(9):67 675, September
FIR BAND-PASS DIGITAL DIFFERENTIATORS WITH FLAT PASSBAND AND EQUIRIPPLE STOPBAND CHARACTERISTICS. T. Yoshida, Y. Sugiura, N.
FIR BAND-PASS DIGITAL DIFFERENTIATORS WITH FLAT PASSBAND AND EQUIRIPPLE STOPBAND CHARACTERISTICS T. Yoshida, Y. Sugiura, N. Aikawa Tokyo University of Science Faculty of Industrial Science and Technology
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationClosed-Form Design of Maximally Flat IIR Half-Band Filters
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 49, NO. 6, JUNE 2002 409 Closed-Form Design of Maximally Flat IIR Half-B Filters Xi Zhang, Senior Member, IEEE,
More informationButterworth Filter Properties
OpenStax-CNX module: m693 Butterworth Filter Properties C. Sidney Burrus This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 3. This section develops the properties
More informationEFFICIENT REMEZ ALGORITHMS FOR THE DESIGN OF NONRECURSIVE FILTERS
EFFICIENT REMEZ ALGORITHMS FOR THE DESIGN OF NONRECURSIVE FILTERS Copyright 2003- Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou EFFICIENT
More informationFuliband Differentiators
Design of FIR LS Hibert Transformers Through Fuliband Differentiators Guergana Mollova* Institute of Communications and Radio-Frequency Engineering Vienna University of Technology Gusshausstrasse 25/389,
More informationOn the Frequency-Domain Properties of Savitzky-Golay Filters
On the Frequency-Domain Properties of Savitzky-Golay Filters Ronald W Schafer HP Laboratories HPL-2-9 Keyword(s): Savitzky-Golay filter, least-squares polynomial approximation, smoothing Abstract: This
More informationDesign of Stable IIR filters with prescribed flatness and approximately linear phase
Design of Stable IIR filters with prescribed flatness and approximately linear phase YASUNORI SUGITA Nagaoka University of Technology Dept. of Electrical Engineering Nagaoka city, Niigata-pref., JAPAN
More information4214 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 54, NO. 11, NOVEMBER 2006
4214 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 54, NO. 11, NOVEMBER 2006 Closed-Form Design of Generalized Maxflat R-Regular FIR M th-band Filters Using Waveform Moments Xi Zhang, Senior Member, IEEE,
More informationUsing fractional delay to control the magnitudes and phases of integrators and differentiators
Using fractional delay to control the magnitudes and phases of integrators and differentiators M.A. Al-Alaoui Abstract: The use of fractional delay to control the magnitudes and phases of integrators and
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #24 Tuesday, November 4, 2003 6.8 IIR Filter Design Properties of IIR Filters: IIR filters may be unstable Causal IIR filters with rational system
More informationFilter-Generating Systems
24 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL. 47, NO. 3, MARCH 2000 Filter-Generating Systems Saed Samadi, Akinori Nishihara, and Hiroshi Iwakura Abstract
More informationDigital Wideband Integrators with Matching Phase and Arbitrarily Accurate Magnitude Response (Extended Version)
Digital Wideband Integrators with Matching Phase and Arbitrarily Accurate Magnitude Response (Extended Version) Ça gatay Candan Department of Electrical Engineering, METU, Ankara, Turkey ccandan@metu.edu.tr
More informationFilter Design Problem
Filter Design Problem Design of frequency-selective filters usually starts with a specification of their frequency response function. Practical filters have passband and stopband ripples, while exhibiting
More informationSHAPE OF IMPULSE RESPONSE CHARACTERISTICS OF LINEAR-PHASE NONRECURSIVE 2D FIR FILTER FUNCTION *
FACTA UNIVERSITATIS Series: Electronics and Energetics Vol. 6, N o, August 013, pp. 133-143 DOI: 10.98/FUEE130133P SHAPE OF IMPULSE RESPONSE CHARACTERISTICS OF LINEAR-PHASE NONRECURSIVE D FIR FILTER FUNCTION
More informationDesign of FIR Nyquist Filters with Low Group Delay
454 IEEE TRASACTIOS O SIGAL PROCESSIG, VOL. 47, O. 5, MAY 999 Design of FIR yquist Filters with Low Group Delay Xi Zhang and Toshinori Yoshikawa Abstract A new method is proposed for designing FIR yquist
More informationLINEAR-PHASE FIR FILTER DESIGN BY LINEAR PROGRAMMING
LINEAR-PHASE FIR FILTER DESIGN BY LINEAR PROGRAMMING Often it is desirable that an FIR filter be designed to minimize the Chebyshev error subject to linear constraints that the Parks- McClellan algorithm
More informationA Higher-Density Discrete Wavelet Transform
A Higher-Density Discrete Wavelet Transform Ivan W. Selesnick Abstract In this paper, we describe a new set of dyadic wavelet frames with three generators, ψ i (t), i =,, 3. The construction is simple,
More informationCOMPARISON OF CLASSICAL CIC AND A NEW CLASS OF STOPBAND-IMPROVED CIC FILTERS FORMED BY CASCADING NON-IDENTICAL COMB SECTIONS
FACTA UIVERSITATIS Series: Electronics and Energetics Vol. 29, o 1, March 216, pp. 61-76 DOI: 1.2298/FUEE16161M COMPARISO OF CLASSICAL CIC AD A EW CLASS OF STOPBAD-IMPROVED CIC FILTERS FORMED BY CASCADIG
More informationMinimax Design of Complex-Coefficient FIR Filters with Low Group Delay
Minimax Design of Complex-Coefficient FIR Filters with Low Group Delay Wu-Sheng Lu Takao Hinamoto Dept. of Elec. and Comp. Engineering Graduate School of Engineering University of Victoria Hiroshima University
More informationISSN (Print) Research Article. *Corresponding author Nitin Rawal
Scholars Journal of Engineering and Technology (SJET) Sch. J. Eng. Tech., 016; 4(1):89-94 Scholars Academic and Scientific Publisher (An International Publisher for Academic and Scientific Resources) www.saspublisher.com
More informationTwo Dimensional Linear Phase Multiband Chebyshev FIR Filter
Two Dimensional Linear Phase Multiband Chebyshev FIR Filter Vinay Kumar, Bhooshan Sunil To cite this version: Vinay Kumar, Bhooshan Sunil. Two Dimensional Linear Phase Multiband Chebyshev FIR Filter. Acta
More informationDesign of Orthonormal Wavelet Filter Banks Using the Remez Exchange Algorithm
Electronics and Communications in Japan, Part 3, Vol. 81, No. 6, 1998 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J80-A, No. 9, September 1997, pp. 1396 1402 Design of Orthonormal Wavelet
More informationA DESIGN OF FIR FILTERS WITH VARIABLE NOTCHES CONSIDERING REDUCTION METHOD OF POLYNOMIAL COEFFICIENTS FOR REAL-TIME SIGNAL PROCESSING
International Journal of Innovative Computing, Information and Control ICIC International c 23 ISSN 349-498 Volume 9, Number 9, September 23 pp. 3527 3536 A DESIGN OF FIR FILTERS WITH VARIABLE NOTCHES
More informationModified Pole Re-position Technique for Optimal IIR Multiple Notch Filter Design
Modified Pole Re-position Technique for Optimal IIR Multiple Notch Filter Design 7 Modified Pole Re-position Technique for Optimal IIR Multiple Notch Filter Design Amnart Thamrongmas and Chalie Charoenlarpnopparut,
More informationPower-Efficient Linear Phase FIR Notch Filter Design Using the LARS Scheme
Wei Xu, Jiaxiang Zhao, Hongie Wang, Chao Gu Power-Efficient Linear Phase FIR Notch Filter Design Using the LARS Scheme WEI XU Tianin Polytechnic University Institute of Electronics and Information Engineering
More informationDigital Signal Processing Lecture 9 - Design of Digital Filters - FIR
Digital Signal Processing - Design of Digital Filters - FIR Electrical Engineering and Computer Science University of Tennessee, Knoxville November 3, 2015 Overview 1 2 3 4 Roadmap Introduction Discrete-time
More informationINF3440/INF4440. Design of digital filters
Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 INF3440/INF4440. Design of digital filters October 2004 Last week lecture Today s lecture: Chapter 8.1-8.3, 8.4.2, 8.5.3 Last lectures:
More informationECE 410 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter 12
. ECE 40 DIGITAL SIGNAL PROCESSING D. Munson University of Illinois Chapter IIR Filter Design ) Based on Analog Prototype a) Impulse invariant design b) Bilinear transformation ( ) ~ widely used ) Computer-Aided
More informationLike bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform.
Inversion of the z-transform Focus on rational z-transform of z 1. Apply partial fraction expansion. Like bilateral Laplace transforms, ROC must be used to determine a unique inverse z-transform. Let X(z)
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationResearch Article Design of One-Dimensional Linear Phase Digital IIR Filters Using Orthogonal Polynomials
International Scholarly Research Network ISRN Signal Processing Volume, Article ID 8776, 7 pages doi:.54//8776 Research Article Design of One-Dimensional Linear Phase Digital IIR Filters Using Orthogonal
More informationOptimal Design of Real and Complex Minimum Phase Digital FIR Filters
Optimal Design of Real and Complex Minimum Phase Digital FIR Filters Niranjan Damera-Venkata and Brian L. Evans Embedded Signal Processing Laboratory Dept. of Electrical and Computer Engineering The University
More informationUNIT - III PART A. 2. Mention any two techniques for digitizing the transfer function of an analog filter?
UNIT - III PART A. Mention the important features of the IIR filters? i) The physically realizable IIR filters does not have linear phase. ii) The IIR filter specification includes the desired characteristics
More informationDesign of Biorthogonal FIR Linear Phase Filter Banks with Structurally Perfect Reconstruction
Electronics and Communications in Japan, Part 3, Vol. 82, No. 1, 1999 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J81-A, No. 1, January 1998, pp. 17 23 Design of Biorthogonal FIR Linear
More informationReview of Fundamentals of Digital Signal Processing
Chapter 2 Review of Fundamentals of Digital Signal Processing 2.1 (a) This system is not linear (the constant term makes it non linear) but is shift-invariant (b) This system is linear but not shift-invariant
More informationFourier Series Representation of
Fourier Series Representation of Periodic Signals Rui Wang, Assistant professor Dept. of Information and Communication Tongji University it Email: ruiwang@tongji.edu.cn Outline The response of LIT system
More informationFilter Analysis and Design
Filter Analysis and Design Butterworth Filters Butterworth filters have a transfer function whose squared magnitude has the form H a ( jω ) 2 = 1 ( ) 2n. 1+ ω / ω c * M. J. Roberts - All Rights Reserved
More informationSYNTHESIS OF BIRECIPROCAL WAVE DIGITAL FILTERS WITH EQUIRIPPLE AMPLITUDE AND PHASE
SYNTHESIS OF BIRECIPROCAL WAVE DIGITAL FILTERS WITH EQUIRIPPLE AMPLITUDE AND PHASE M. Yaseen Dept. of Electrical and Electronic Eng., University of Assiut Assiut, Egypt Tel: 088-336488 Fax: 088-33553 E-Mail
More informationDesign of Narrow Stopband Recursive Digital Filter
FACTA UNIVERSITATIS (NIŠ) SER.: ELEC. ENERG. vol. 24, no. 1, April 211, 119-13 Design of Narrow Stopband Recursive Digital Filter Goran Stančić and Saša Nikolić Abstract: The procedure for design of narrow
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #21 Friday, October 24, 2003 Types of causal FIR (generalized) linear-phase filters: Type I: Symmetric impulse response: with order M an even
More informationChapter 7: Filter Design 7.1 Practical Filter Terminology
hapter 7: Filter Design 7. Practical Filter Terminology Analog and digital filters and their designs constitute one of the major emphasis areas in signal processing and communication systems. This is due
More informationAPPLIED SIGNAL PROCESSING
APPLIED SIGNAL PROCESSING DIGITAL FILTERS Digital filters are discrete-time linear systems { x[n] } G { y[n] } Impulse response: y[n] = h[0]x[n] + h[1]x[n 1] + 2 DIGITAL FILTER TYPES FIR (Finite Impulse
More informationReview of Fundamentals of Digital Signal Processing
Solution Manual for Theory and Applications of Digital Speech Processing by Lawrence Rabiner and Ronald Schafer Click here to Purchase full Solution Manual at http://solutionmanuals.info Link download
More informationECSE 512 Digital Signal Processing I Fall 2010 FINAL EXAMINATION
FINAL EXAMINATION 9:00 am 12:00 pm, December 20, 2010 Duration: 180 minutes Examiner: Prof. M. Vu Assoc. Examiner: Prof. B. Champagne There are 6 questions for a total of 120 points. This is a closed book
More informationDigital Signal Processing Lecture 8 - Filter Design - IIR
Digital Signal Processing - Filter Design - IIR Electrical Engineering and Computer Science University of Tennessee, Knoxville October 20, 2015 Overview 1 2 3 4 5 6 Roadmap Discrete-time signals and systems
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science Discrete-Time Signal Processing Fall 2005
1 MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.341 Discrete-Time Signal Processing Fall 2005 FINAL EXAM Friday, December 16, 2005 Walker (50-340) 1:30pm
More informationStability Condition in Terms of the Pole Locations
Stability Condition in Terms of the Pole Locations A causal LTI digital filter is BIBO stable if and only if its impulse response h[n] is absolutely summable, i.e., 1 = S h [ n] < n= We now develop a stability
More informationOn the Designing of Fractional Order FIR Differentiator Using Radial Basis Function and Window
On the Designing of Fractional Order FIR Differentiator Using Radial Basis Function and Window Manjeet Kumar Digital Signal Processing Laboratory, Room No. 135, ECE Division, Netaji Subhas Institute of
More informationT ulation and demodulation. It has applications not only
IEEE TRANSACTIONS ON INSTRUMENTATION AND MEASUREMENT. VOL 39. NO 6. DECEMBER 1990 847 Optimal FIR and IIR Hilbert Transformer Design Via LS and Minimax Fitting ISTVAN KOLLAR, MEMBER, IEEE, RIK PINTELON,
More informationDSP. Chapter-3 : Filter Design. Marc Moonen. Dept. E.E./ESAT-STADIUS, KU Leuven
DSP Chapter-3 : Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ Filter Design Process Step-1 : Define filter specs Pass-band, stop-band,
More informationDigital Filters Ying Sun
Digital Filters Ying Sun Digital filters Finite impulse response (FIR filter: h[n] has a finite numbers of terms. Infinite impulse response (IIR filter: h[n] has infinite numbers of terms. Causal filter:
More informationON FINITE-PRECISION EFFECTS IN LATTICE WAVE-DIGITAL FILTERS
ON FINITE-PRECISION EFFECTS IN LATTICE WAVE-DIGITAL FILTERS Isabel M. Álvarez-Gómez, J.D. Osés and Antonio Álvarez-Vellisco EUIT de Telecomunicación, Camino de la Arboleda s/n, Universidad Politécnica
More informationAll-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials
RADIOENGINEERING, VOL. 3, NO. 3, SEPTEMBER 4 949 All-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials Nikola STOJANOVIĆ, Negovan STAMENKOVIĆ, Vidosav STOJANOVIĆ University of Niš,
More informationFilter structures ELEC-E5410
Filter structures ELEC-E5410 Contents FIR filter basics Ideal impulse responses Polyphase decomposition Fractional delay by polyphase structure Nyquist filters Half-band filters Gibbs phenomenon Discrete-time
More informationON THE REALIZATION OF 2D LATTICE-LADDER DISCRETE FILTERS
Journal of Circuits Systems and Computers Vol. 3 No. 5 (2004) 5 c World Scientific Publishing Company ON THE REALIZATION OF 2D LATTICE-LADDER DISCRETE FILTERS GEORGE E. ANTONIOU Department of Computer
More informationComputer-Aided Design of Digital Filters. Digital Filters. Digital Filters. Digital Filters. Design of Equiripple Linear-Phase FIR Filters
Computer-Aided Design of Digital Filters The FIR filter design techniques discussed so far can be easily implemented on a computer In addition, there are a number of FIR filter design algorithms that rely
More information1. FIR Filter Design
ELEN E4810: Digital Signal Processing Topic 9: Filter Design: FIR 1. Windowed Impulse Response 2. Window Shapes 3. Design by Iterative Optimization 1 1. FIR Filter Design! FIR filters! no poles (just zeros)!
More informationTowards Global Design of Orthogonal Filter Banks and Wavelets
Towards Global Design of Orthogonal Filter Banks and Wavelets Jie Yan and Wu-Sheng Lu Department of Electrical and Computer Engineering University of Victoria Victoria, BC, Canada V8W 3P6 jyan@ece.uvic.ca,
More informationAnalog and Digital Filter Design
Analog and Digital Filter Design by Jens Hee http://jenshee.dk October 208 Change log 28. september 208. Document started.. october 208. Figures added. 6. october 208. Bilinear transform chapter extended.
More informationVARIOUS types of wavelet transform are available for
IEEE TRANSACTIONS ON SIGNAL PROCESSING A Higher-Density Discrete Wavelet Transform Ivan W. Selesnick, Member, IEEE Abstract This paper describes a new set of dyadic wavelet frames with two generators.
More informationA Simple Proof of the Alternation Theorem
A Simple Proof of the Alternation Theorem P. P. Vaidyanathan Dept. of Electrical Engineering California Institute of Technology Pasadena, CA 91125 Email: ppvnath@systems.caltech.edu T. Q. Nguyen Dept.
More informationDesign of Coprime DFT Arrays and Filter Banks
Design o Coprime DFT Arrays and Filter Banks Chun-Lin Liu clliu@caltechedu P P Vaidyanathan ppvnath@systemscaltechedu Digital Signal Processing Group Electrical Engineering Caliornia Institute o Technology
More informationDSP-CIS. Chapter-4: FIR & IIR Filter Design. Marc Moonen
DSP-CIS Chapter-4: FIR & IIR Filter Design Marc Moonen Dept. E.E./ESAT-STADIUS, KU Leuven marc.moonen@esat.kuleuven.be www.esat.kuleuven.be/stadius/ PART-II : Filter Design/Realization Step-1 : Define
More informationPerfect Reconstruction Two- Channel FIR Filter Banks
Perfect Reconstruction Two- Channel FIR Filter Banks A perfect reconstruction two-channel FIR filter bank with linear-phase FIR filters can be designed if the power-complementary requirement e jω + e jω
More informationDIGITAL SIGNAL PROCESSING. Chapter 6 IIR Filter Design
DIGITAL SIGNAL PROCESSING Chapter 6 IIR Filter Design OER Digital Signal Processing by Dr. Norizam Sulaiman work is under licensed Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International
More informationDigital Control & Digital Filters. Lectures 21 & 22
Digital Controls & Digital Filters Lectures 2 & 22, Professor Department of Electrical and Computer Engineering Colorado State University Spring 205 Review of Analog Filters-Cont. Types of Analog Filters:
More informationON THE USE OF GEGENBAUER PROTOTYPES IN THE SYNTHESIS OF WAVEGUIDE FILTERS
Progress In Electromagnetics Research C, Vol. 18, 185 195, 2011 ON THE USE OF GEGENBAUER PROTOTYPES IN THE SYNTHESIS OF WAVEGUIDE FILTERS L. Cifola, A. Morini, and G. Venanzoni Dipartimento di Ingegneria
More informationThe DFT as Convolution or Filtering
Connexions module: m16328 1 The DFT as Convolution or Filtering C. Sidney Burrus This work is produced by The Connexions Project and licensed under the Creative Commons Attribution License A major application
More informationUNIVERSITY OF OSLO. Faculty of mathematics and natural sciences. Forslag til fasit, versjon-01: Problem 1 Signals and systems.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 1th, 016 Examination hours: 14:30 18.30 This problem set
More informationIterative reweighted l 1 design of sparse FIR filters
Iterative reweighted l 1 design of sparse FIR filters Cristian Rusu, Bogdan Dumitrescu Abstract Designing sparse 1D and 2D filters has been the object of research in recent years due mainly to the developments
More informationExercises in Digital Signal Processing
Exercises in Digital Signal Processing Ivan W. Selesnick September, 5 Contents The Discrete Fourier Transform The Fast Fourier Transform 8 3 Filters and Review 4 Linear-Phase FIR Digital Filters 5 5 Windows
More informationSignals and Systems. Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI
Signals and Systems Lecture 11 Wednesday 22 nd November 2017 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Effect on poles and zeros on frequency response
More informationEE 521: Instrumentation and Measurements
Aly El-Osery Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA November 1, 2009 1 / 27 1 The z-transform 2 Linear Time-Invariant System 3 Filter Design IIR Filters FIR Filters
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing 6. DSP - Recursive (IIR) Digital Filters Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher: John Wiley
More informationSignal Processing. Lecture 10: FIR Filter Design. Ahmet Taha Koru, Ph. D. Yildiz Technical University Fall
Signal Processing Lecture 10: FIR Filter Design Ahmet Taha Koru, Ph. D. Yildiz Technical University 2017-2018 Fall ATK (YTU) Signal Processing 2017-2018 Fall 1 / 47 Introduction Introduction ATK (YTU)
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationECE 8440 Unit 17. Parks- McClellan Algorithm
Parks- McClellan Algorithm ECE 8440 Unit 17 The Parks- McClellan Algorithm is a computer method to find the unit sample response h(n for an op=mum FIR filter that sa=sfies the condi=ons of the Alterna=on
More informationEfficient 2D Linear-Phase IIR Filter Design and Application in Image Filtering
Efficient D Linear-Phase IIR Filter Design and Application in Image Filtering Chi-Un Lei, Chung-Man Cheung, and Ngai Wong Abstract We present an efficient and novel procedure to design two-dimensional
More information1 1.27z z 2. 1 z H 2
E481 Digital Signal Processing Exam Date: Thursday -1-1 16:15 18:45 Final Exam - Solutions Dan Ellis 1. (a) In this direct-form II second-order-section filter, the first stage has
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationOn the Dual-Tree Complex Wavelet Packet and. Abstract. Index Terms I. INTRODUCTION
On the Dual-Tree Complex Wavelet Packet and M-Band Transforms İlker Bayram, Student Member, IEEE, and Ivan W. Selesnick, Member, IEEE Abstract The -band discrete wavelet transform (DWT) provides an octave-band
More informationECE503: Digital Signal Processing Lecture 5
ECE53: Digital Signal Processing Lecture 5 D. Richard Brown III WPI 3-February-22 WPI D. Richard Brown III 3-February-22 / 32 Lecture 5 Topics. Magnitude and phase characterization of transfer functions
More informationComputing running DCTs and DSTs based on their second-order shift properties
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering Information Sciences 000 Computing running DCTs DSTs based on their second-order shift properties
More informationONE-DIMENSIONAL (1-D) two-channel FIR perfect-reconstruction
3542 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL. 55, NO. 11, DECEMBER 2008 Eigenfilter Approach to the Design of One-Dimensional and Multidimensional Two-Channel Linear-Phase FIR
More informationAdaptive MMSE Equalizer with Optimum Tap-length and Decision Delay
Adaptive MMSE Equalizer with Optimum Tap-length and Decision Delay Yu Gong, Xia Hong and Khalid F. Abu-Salim School of Systems Engineering The University of Reading, Reading RG6 6AY, UK E-mail: {y.gong,x.hong,k.f.abusalem}@reading.ac.uk
More informationDESIGN OF QUANTIZED FIR FILTER USING COMPENSATING ZEROS
DESIGN OF QUANTIZED FIR FILTER USING COMPENSATING ZEROS Nivedita Yadav, O.P. Singh, Ashish Dixit Department of Electronics and Communication Engineering, Amity University, Lucknow Campus, Lucknow, (India)
More informationDesign of Optimum High Order Finite Wordlength Digital FIR Filters with Linear Phase
EURASIP SIGNAL PROCESSING 1 Design of Optimum High Order Finite Wordlength Digital FIR Filters with Linear Phase Gennaro Evangelista Abstract A novel iterative quantization procedure for the design of
More informationGENERALIZED DEFLATION ALGORITHMS FOR THE BLIND SOURCE-FACTOR SEPARATION OF MIMO-FIR CHANNELS. Mitsuru Kawamoto 1,2 and Yujiro Inouye 1
GENERALIZED DEFLATION ALGORITHMS FOR THE BLIND SOURCE-FACTOR SEPARATION OF MIMO-FIR CHANNELS Mitsuru Kawamoto,2 and Yuiro Inouye. Dept. of Electronic and Control Systems Engineering, Shimane University,
More informationThe Selection of Weighting Functions For Linear Arrays Using Different Techniques
The Selection of Weighting Functions For Linear Arrays Using Different Techniques Dr. S. Ravishankar H. V. Kumaraswamy B. D. Satish 17 Apr 2007 Rajarshi Shahu College of Engineering National Conference
More informationChapter 7: IIR Filter Design Techniques
IUST-EE Chapter 7: IIR Filter Design Techniques Contents Performance Specifications Pole-Zero Placement Method Impulse Invariant Method Bilinear Transformation Classical Analog Filters DSP-Shokouhi Advantages
More informationDIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS. 3.6 Design of Digital Filter using Digital to Digital
DIGITAL SIGNAL PROCESSING UNIT III INFINITE IMPULSE RESPONSE DIGITAL FILTERS Contents: 3.1 Introduction IIR Filters 3.2 Transformation Function Derivation 3.3 Review of Analog IIR Filters 3.3.1 Butterworth
More informationCHEBYSHEV polynomials have been an essential mathematical
IEEE TRANSACTIONS ON COMPUTERS, VOL X, NO X, MONTH YYYY 1 A Karatsuba-based Algorithm for Polynomial Multiplication in Chebyshev Form Juliano B Lima, Student Member, IEEE, Daniel Panario, Member, IEEE,
More informationFIR Filter Design via Semidenite Programming and Spectral Factorization Shao-Po Wu, Stephen Boyd, Lieven Vandenberghe Information Systems Laboratory Stanford University, Stanford, CA 9435 clive@isl.stanford.edu,
More informationINFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters 4. THE BUTTERWORTH ANALOG FILTER
INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters. INTRODUCTION 2. IIR FILTER DESIGN 3. ANALOG FILTERS 4. THE BUTTERWORTH ANALOG FILTER 5. THE CHEBYSHEV-I
More informationPOCHHAMMER SYMBOL WITH NEGATIVE INDICES. A NEW RULE FOR THE METHOD OF BRACKETS.
POCHHAMMER SYMBOL WITH NEGATIVE INDICES A NEW RULE FOR THE METHOD OF BRACKETS IVAN GONZALEZ, LIN JIU, AND VICTOR H MOLL Abstract The method of brackets is a method of integration based upon a small number
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to decrease the sampling rate by an integer
More informationElectronic Circuits EE359A
Electronic Circuits EE359A Bruce McNair B26 bmcnair@stevens.edu 21-216-5549 Lecture 22 578 Second order LCR resonator-poles V o I 1 1 = = Y 1 1 + sc + sl R s = C 2 s 1 s + + CR LC s = C 2 sω 2 s + + ω
More informationLecture 7 Discrete Systems
Lecture 7 Discrete Systems EE 52: Instrumentation and Measurements Lecture Notes Update on November, 29 Aly El-Osery, Electrical Engineering Dept., New Mexico Tech 7. Contents The z-transform 2 Linear
More information