Filter Design Problem
|
|
- Claribel Walsh
- 5 years ago
- Views:
Transcription
1 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 a gradual roll-off in the transition band. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 1/40
2 Filter specifications Absolute specifications In the passband: 1 δ p H(e jω ) 1 + δ p, for 0 ω ω p, with δ p 1. In the stopband: H(e jω ) δ s, for ω s ω π, with δ s 1. Relative specifications 1 δ p 1 + δ p H(e jω ) 1, 0 ω ω p, H(e jω ) δ s 1 + δ p, ω s ω π. The passband A p and stopband A s ripples are defined as ( ) 1 + δp A p 20 log 10 1 δ p and A s 20 log 10 ( 1 + δp δ s ). Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 2/40
3 Continuous-time filter specification In practical applications, the passband F p and stopband F s cutoff frequencies are specified in Hz. The values of ω p and ω s are calculated according to ω p = 2π F pass F s, ω s = 2π F stop F s. IIR filters are typically derived from analog filters, which are traditionally specified using the quantities ɛ and A as follows. ( ) 20 log ɛ 2 = A p and 20 log 10 A = A s, which gives ɛ = Ap 1 and A = As. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 3/40
4 Optimality criteria for filter design Let H(z) be a rational system function that approximates some ideal frequency response H d (e jω ). The coefficients of H(z) depends on the approximation criterion. Mean-squared-error approximation: E 2 [ 1/2 1 H(e jω ) H d (e jω ) dω] 2, 2π B with B being a frequency interval of interest. Minimax approximation: E max ω B H(ejω ) H d (e jω ). The solution minimizing this error function is called a minimax or Chebyshev approximation. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 4/40
5 FIR filters with linear phase Without loss of generality, we focus on the ideal lowpass filter (LPF) with its impulse and frequency responses given by { h lp [n] = sin ω c(n α), H lp (e jω e jαω, ω ω c, ) = π(n α) 0, ω c < ω π. Depending on the choice of α, there are three possible cases 1 If α Z, h lp [n] is symmetric about its sample at n = α. 2 If 2α Z, h lp [n] is symmetric about the middle point between the samples at α 1/2 and α + 1/2. 3 If 2α is not an integer, there is no symmetry at all. A causal FIR filter is created by setting h[n] = h lp [n] for n = 0,..., M and zero elsewhere. There are four types of FIR filter with linear phase. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 5/40
6 Type-I FIR linear-phase filters A Type-I FIR system has a symmetric impulse response with even order M, that is, h[n] = h[m n], 0 n M. The frequency response of a Type-I FIR filter can be expressed as M/2 H(e jω ) = a[k] cos ωk e j(m/2)ω = A(e jω )e j(m/2)ω, k=0 } {{ } A(e jω ) where A(e jω ) is a real, even, and periodic function of ω with a[0] = h[m/2], a[k] = 2h[(M/2) k], k = 1,..., M/2. This type of filters is the most general one. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 6/40
7 Type-II FIR linear-phase filters A Type-I FIR system has a symmetric impulse response with odd order M. The frequency response of a Type-II FIR filter can be expressed as (M+1)/2 H(e jω ) = b[k] cos [ω(k 1/2)] e j(m/2)ω = k=1 } {{ } A(e jω ) where the delay M/2 is fractional and = A(e jω )e j(m/2)ω, b[k] = 2h[(M + 1)/2 k], k = 1,..., (M + 1)/2. We note that at ω = π, A(e jω ) = 0, independently of h[k]. This implies that Type-II filters cannot be high-pass. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 7/40
8 Type-III FIR linear-phase filters A Type-III FIR system has an antisymmetric impulse response with even order M, i.e. h[n] = h[m n], 0 n M. In this case, the frequency response is given by M/2 H(e jω ) = c[k] sin ωk j e j(m/2)ω = j A(e jω )e j(m/2)ω, where k=1 } {{ } A(e jω ) c[k] = 2h[(M/2) k], k = 1,..., M/2. The condition h[m/2] = h[m/2] requires that h[m/2] = 0 for every Type-III filter. Note that A(e jω ) = 0, at ω {0, π}. Such filters are most suitable for designing differentiators and Hilbert transformers. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 8/40
9 Type-IV FIR linear-phase filters If the impulse response is antisymmetric and M is odd, then we have (M+1)/2 H(e jω ) = d[k] sin [ω(k 1/2)] j e j(m/2)ω = k=1 } {{ } A(e jω ) where the delay M/2 is fractional and = ja(e jω )e j(m/2)ω, d[k] = 2h[(M + 1)/2 k], k = 1,..., (M + 1)/2. For Type-IV filters, ones has A(e jω ) = 0 at ω = 0, independently of the choice of h[k]. Thus, this class of filters is most suitable for approximating differentiators and Hilbert transformers. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 9/40
10 Frequency responses of the FIR filters The frequency responses of all Type I-IV FIR filters with linear phase can be expressed in the form H(e jω ) = M h[n]e jωn = A(e jω )e jψ(ejω), n=0 where the amplitude response A(e jω ) is real and Ψ(e jω ) αω + β. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 10/40
11 Frequency responses of the FIR filters (cont.) The constant β results from the presence of j in H(e jω ), and it does not affect the group delay τ gd (ω) = dψ(ejω ) dω = α. In contrast to H(e jω, A(e jω ) is analytic, that is, its derivative exists for all ω. As a result, in contrast to H(e jω ), Ψ(e jω ) a continuous function of ω. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 11/40
12 Additional insights Ψ(e jω ) = αω + β is called a generalized linear phase. Recall the following Fourier transform relations: h e [n] A e (e jω ) = h e [n n 0 ] A e (e jω )e jωn0. h o [n] ja o (e jω ) = h o [n n 0 ] A o (e jω )e jωn0+jπ/2. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 12/40
13 Design of FIR filters by windowing Suppose we wish to approximate an ideal frequency response H d (e jω ) = h d [n]e jωn n= with an FIR filter h[n], n = 0,..., M by minimizing the MSE ɛ 2 = 1 2π π π H(e jω ) H d (e jω ) 2 dω. Using Parsevals identity, one can express ɛ 2 as ɛ 2 = M (h[n] h d [n]) 2 + n=0 1 n= (h d [n]) 2 + Thus, the optimal solution is { h d [n], n = 0,..., M h[n] =, 0, otherwise which is nothing else but a truncation of h d [n]. n=m+1 (h d [n]) 2 Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 13/40
14 Truncation in the DTFT domain The truncated impulse response is the product of h d [n] with which suggests w[n] = u[n] u[n M 1], H(e jω ) = 1 2π π π H d (e jθ )W (e j(ω θ) )dθ, where W (e jω ) = sin [ω(m + 1)/2] e jωm/2. sin ω/2 The convolution produces passband and stopband ripples which is a result of the Gibbs phenomenon. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 14/40
15 Truncation in the DTFT domain (cont.) Truncating h d [n] yields passband ripples about 20 log 10 (1.0895) = = 0.75 db, as well as a minimum stopband attenuation of about 20 log 10 (1/0.0895) = 21 db, irrespective of the filter length M + 1. Both are not sufficient for practical applications. The design can be improved by using smoother windows. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 15/40
16 Filter design by windowing Below are the most commonly used windows. Bartlett: Hann: Hamming: w[n] = w[n] = 2n/M, 0 n M/2, w[n] = 2 2n/M, M/2 < n M,. 0, otherwise { cos(2πn/m), n = 0,..., M, 0, otherwise. { cos(2πn/m), n = 0,..., M, 0, otherwise. Blackman: { cos(2πn/m) cos(4πn/m), n = 0,..., M, w[n] = 0, otherwise. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 16/40
17 Fixed windows Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 17/40
18 Fixed windows (cont.) Windows are always symmetric, that is { w[m n], 0 n M, w[n] = 0, otherwise. Therefore, windowing an FIR filter with linear phase does not change its type. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 18/40
19 Kaiser window We see that there is a trade-off between main-lobe width and side-lobe height; we cannot reduce both quantities at the same time. Kaiser window maximizes the ratio of the energy in a frequency band about ω = 0 over the total energy of the window. Formally, [ ] I 0 β 1 [(n α)/α] 2 w[n] = I 0(β), 0 n M, 0, otherwise, where α = M/2 and I 0 is the zeroth-order modifier Bessel function of the first kind. The shape of the Kaiser window is controlled by the shape parameter β (in MATLAB: kaiser(l,beta)), which controls the aforementioned trade-off. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 19/40
20 Kaiser window (cont.) Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 20/40
21 Kaiser window (cont.) A prescribed value A for A s is related to β according to 0, A < 21, β = (A 21) (A 21), 21 A 50, (A 8.7), A > 50, where A = 20 log 10 δ, with δ = max{δ p, δ s } (since the windowing method produces filters with δ p δ s, irrespective of the design specifications.) The range 0 β 8 provides useful windows. The order M required to achieve prescribed values of A and ω is estimated as M = A ω. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 21/40
22 FIR filter design using fixed windows 1 Given the design specifications {ω p, ω s, A p, A s }, determine the ripples δ p and δ s and set δ = min{δ p, δ s }. 2 Determine the cutoff frequency of the ideal LPF as ω c = ωp+ωs 2. 3 Determine the parameters A = 20 log 10 δ and ω = ω s ω p. 4 From the table, choose the window function that provides the smallest stopband attenuation greater than A. 5 For this window function, determine the required value of M by selecting the corresponding value of ω. 6 Determine the impulse response of the ideal lowpass filter by h d [n] = sin [ω c(n M/2)]. π(n M/2) 7 Compute the impulse response h[n] = h d [n] w[n] using the chosen window. 8 Check whether the designed filter satisfies the prescribed specifications; if not, increase M and go back to step 6. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 22/40
23 Example: M = 40, ω c = π/4 Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 23/40
24 Multiband filter design The above approach is straightforward to generalize to the design of highpass, bandpass, bandstop, and other multiband filters. For example, the magnitude response of a 5-band ideal multiband filter can look like In general, the impulse response of a multiband filter is given by h mb [n] = K (A k A k+1 ) sin [ω k(n M/2)], π(n M/2) k=1 where K is the number of bands. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 24/40
25 Design of FIR filters by frequency sampling Suppose we are given L samples of a desired frequency response H d [k] = H d (e j2πk/l ), k = 0, 1,..., L 1. The associated impulse response is given by h d [n] = 1 L which is L periodic. L 1 H d [k]w nk L = h d [n ml], m Z k=0 Multiplying h d with a window of length L results in with its DTFT given by h[n] = h d [n]w[n], H(e jω ) = 1 L 1 H d [k]w (e j(ω 2πk/L) ). L k=0 Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 25/40
26 Design by frequency sampling (cont.) Thus, the frequency response of the designed filter is obtained by interpolating between the samples H d [k] using W (e jω ) as an interpolation function. The windowing method creates an FIR filter by truncating h d [n]; the frequency sampling method creates an FIR filter by aliasing or folding h d [n]. If w[n] is an L-point rectangular window then h[n] is the primary period of h[n]. In this case the approximation error is zero at the sampling frequencies, and finite between them. For nonrectangular windows, the approximation error is not zero at the sampling frequencies. In practice, to minimize the time-domain aliasing distortion, we start with a large value for L, and then we choose a much smaller value for the length of the window w[n]. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 26/40
27 Design by frequency sampling (cont.) To guarantee the resulting filter has linear phase, we use the following formulas H d [k] = A d [k]e jψ d[k], A d [k] = Ψ d [k] = { A d (e j0 ), k = 0, A d (e j2π(l k)/l ), k = 1,..., L, { ± π 2 q L 1 2 π 2 q + L 1 2 where Q = (L 1)/2. 2π L 2π L k, k = 0, 1,..., Q, (L k), k = Q + 1,..., L 1, Note that the parameter q = 0 for type I-II FIR filters, and q = 1 for type III-IV FIR filters. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 27/40
28 Example: w c = 0.45π, M = 19 Note that the samples of the frequency response transition sharply from passband to stopband. Also, the approximation error is larger near the sharp transition and is smaller away from it. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 28/40
29 Better design approaches Since the transition band is typically unspecified, we can improve the quality of the approximation by enforcing a smoother transition band. For many applications, it is sufficient to use either a straight-line roll-off A(e jω ) = (ω s ω)/(ω s ω p ), ω p < ω < ω s, or a raised-cosine roll-off A(e jω ) = cos(π(ω s ω)/(ω s ω p )), ω p < ω < ω s. The main drawback is to create a wider transition bandwidth in the resulting design which can be reduced by increasing the order M. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 29/40
30 Example: w c = 0.45π, M = 19 Notice a substantial reduction in passband and stopband ripples. The ripples obtained using the raised-cosine approach are almost non-noticeable. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 30/40
31 Example: w c = 0.45π, M = 19 Using non-rectangular windows can also be used to reduce passband and stopband ripple. The resulting frequency response is similar to the one obtained using the raised-cosine smooth transition. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 31/40
32 Design procedure To design a lowpass filter with specifications ω p, ω s, A p, and A s, the following steps can be used: 1 Choose the order of the filter M by placing at least two samples in the transition band. 2 For a window design approach obtain samples of the desired frequency response H d [k]. For a smooth transition band approach, add transition band samples. 3 Compute the (M + 1)-point IDFT of H d [k] to obtain h[n]. For a window design approach multiply h[n] by the appropriate w[n]. 4 Compute the log-magnitude response and verify the design over passband and stopband. 5 If the specifications aren t met, increase M and go back to step 1. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 32/40
33 Example Consider the following specifications of a lowpass filter: ω p = 0.25π, ω s = 0.35π, A p = 0.1 db, A s = 40 db. Since the transition bandwidth is 0.1π, we will need more than 40 samples of the ideal frequency response. M = 44; L = M+1; alpha = M/2; Q = floor(alpha); Psi = -alpha*2*pi/l*[(0:q), -(L-(Q+1:M))]; Dw = 2*pi/L; k1 = floor(wp/dw); k2 = ceil(ws/dw); w = ((k2-1):-1:(k1+1))*dw; A = *cos(pi*(ws-w)/(ws-wp)); Ad = [ones(1,k1+1),a,zeros(1,l-2*k2+1),... flip(a),ones(1,k1)]; Hd = Ad.*exp(1j*Psi); h = real(ifft(hd)); %phase delay parameter %desired phase %transition band %trans. band samples Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 33/40
34 Example (cont.) The stopband attenuation (in db) of the designed filter can be measured as: H = freqz(h,1,linspace(0,2*pi,1001)); dw = 2*pi/1000; maxmag = max(abs(h)); Asd = min(-20*log10(abs(h(ceil(ws/dw):501))/maxmag)) Asd = The measured value of the stopband attenuation does not satisfy the given specifications. Clearly, we need to increase M or use another approach. Since a Hamming window provides more than 50 db of attenuation, one can use it instead of the rectangular window. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 34/40
35 Example (cont.) In this case, the script remains largely the same except for: omc = (wp+ws)/2; %cutoff frequency k1 = floor(omc/dw); k2 = ceil(omc/dw); Ad = [ones(1,k1+1),zeros(1,l-2*k2+1),ones(1,k1)]; Hd = Ad.*exp(1j*Psi); h = real(ifft(hd)).*hamming(l)'; %windowed impulse response In this case, the measured value of the stopband attenuation is 37.4 db, which is still insufficient. By increasing the value of M and checking for the resulting stopband attenuation, one can verify that M = 50 satisfies the given specifications for both approaches. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 35/40
36 Example (cont.) The above example shows that the frequency sampling technique is not well suited to standard filter design, because it is difficult, if not impossible, to determine a priori the number of samples needed for the correct design. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 36/40
37 Equiripple filter design Problem formulation Given a desired amplitude response A d (e jω ) and a weighting function W (ω) 0, the actual amplitude response A(e jω ) is characterized by the approximation error defined as E(ω) W (ω) [A d (e jω ) A(e jω )]. The design objective is to find the coefficients of a Type I-IV FIR filter minimizing the weighted Chebyshev error, defined by E = max ω B E(ω), where B is a union of disjoint closed subsets of [0, π]. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 37/40
38 Equiripple filter design (cont.) Since we focus on FIR filters with linear phase, A(e jω ) can be written as A(e jω ) = Q(e jω )P (e jω ) with 1, Type I R P (e jω ) = p[k] cos(ωk) and Q(e jω cos(ω/2), Type II ) = sin(ω), Type III k=0 sin(ω/2), Type IV where R = M/2, if M is even, and R = (M 1)/2, if M is odd. Consequently, the error E(ω) can be redefined as E(ω) = W (ω)q(e jω ) }{{} W (ω) [ Ad (e jω ) Q(e jω ) }{{} Ā d (e jω ) ] P (e jω ) = W ] (ω)[ād (e jω ) P (e jω ). The Chebyshev approximation problem {p[k]} R k=0 = arg min E. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 38/40
39 Alternation theorem for FIR filters The solution to the Chebyshev approximation problem is characterized by the alternation theorem. Theorem (Alternation theorem for FIR filters) A necessary and sufficient condition for P (e jω ) to be the unique solution of the Chebyshev approximation problem is that E(ω) exhibits at least R + 2 alternations in B. That is, there must exist R + 2 extremal frequencies ω 1 < ω 2 <... < ω R+2 such that for every k = 1,..., R + 2 E(ω k ) = E(ω k+1 ), E(ω k ) = max E(ω) δ. ω B The alternation theorem implies that the best Chebyshev approximation must have an equiripple error function. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 39/40
40 Alternation theorem for FIR filters (cont.) The conditions characterizing the optimum filter can be used to develop an efficient algorithm for obtaining filter coefficients. Parks and McClellan (1972) solved this problem by means of the Remez exchange algorithm. Modern methods of filter design are based on the tools of convex optimization. See, for example, FIR Filter Design Factorization and Convex Optimization, by S. P. Wu, S. Boyd, and L. Vandenberghe. Prof. O. Michailovich, Dept of ECE, Winter 2017 ECE 413: Digital Signal Processing / Section 8 40/40
Digital 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 information-Digital Signal Processing- FIR Filter Design. Lecture May-16
-Digital Signal Processing- FIR Filter Design Lecture-17 24-May-16 FIR Filter Design! FIR filters can also be designed from a frequency response specification.! The equivalent sampled impulse response
More informationBasic Design Approaches
(Classic) IIR filter design: Basic Design Approaches. Convert the digital filter specifications into an analog prototype lowpass filter specifications. Determine the analog lowpass filter transfer function
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 informationUNIT - 7: FIR Filter Design
UNIT - 7: FIR Filter Design Dr. Manjunatha. P manjup.jnnce@gmail.com Professor Dept. of ECE J.N.N. College of Engineering, Shimoga October 5, 06 Unit 7 Syllabus Introduction FIR Filter Design:[,, 3, 4]
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 information(Refer Slide Time: 01:28 03:51 min)
Digital Signal Processing Prof. S. C. Dutta Roy Department of Electrical Engineering Indian Institute of Technology, Delhi Lecture 40 FIR Design by Windowing This is the 40 th lecture and our topic for
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 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 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 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 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 informationDesign of IIR filters
Design of IIR filters Standard methods of design of digital infinite impulse response (IIR) filters usually consist of three steps, namely: 1 design of a continuous-time (CT) prototype low-pass filter;
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 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 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 informationTransform Analysis of Linear Time-Invariant Systems
Transform Analysis of Linear Time-Invariant Systems Discrete-Time Signal Processing Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-Sen University Kaohsiung, Taiwan ROC Transform
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 3 Signals & Systems Prof. ark Fowler Note Set #28 D-T Systems: DT Filters Ideal & Practical /4 Ideal D-T Filters Just as in the CT case we can specify filters. We looked at the ideal filter for the
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 informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
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 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 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 informationLABORATORY 3 FINITE IMPULSE RESPONSE FILTERS
LABORATORY 3 FINITE IMPULSE RESPONSE FILTERS 3.. Introduction A digital filter is a discrete system, used with the purpose of changing the amplitude and/or phase spectrum of a signal. The systems (filters)
More informationMultirate signal processing
Multirate signal processing Discrete-time systems with different sampling rates at various parts of the system are called multirate systems. The need for such systems arises in many applications, including
More informationContents. Digital Signal Processing, Part II: Power Spectrum Estimation
Contents Digital Signal Processing, Part II: Power Spectrum Estimation 5. Application of the FFT for 7. Parametric Spectrum Est. Filtering and Spectrum Estimation 7.1 ARMA-Models 5.1 Fast Convolution 7.2
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 informationLINEAR-PHASE FIR FILTERS DESIGN
LINEAR-PHASE FIR FILTERS DESIGN Prof. Siripong Potisuk inimum-phase Filters A digital filter is a minimum-phase filter if and only if all of its zeros lie inside or on the unit circle; otherwise, it is
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 informationThe Discrete-Time Fourier
Chapter 3 The Discrete-Time Fourier Transform 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 3-1-1 Continuous-Time Fourier Transform Definition The CTFT of
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 informationLecture 14: Windowing
Lecture 14: Windowing ECE 401: Signal and Image Analysis University of Illinois 3/29/2017 1 DTFT Review 2 Windowing 3 Practical Windows Outline 1 DTFT Review 2 Windowing 3 Practical Windows DTFT Review
More informationELEG 5173L Digital Signal Processing Ch. 5 Digital Filters
Department of Electrical Engineering University of Aransas ELEG 573L Digital Signal Processing Ch. 5 Digital Filters Dr. Jingxian Wu wuj@uar.edu OUTLINE 2 FIR and IIR Filters Filter Structures Analog Filters
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 informationProblem Set 9 Solutions
Problem Set 9 Solutions EE23: Digital Signal Processing. From Figure below, we see that the DTFT of the windowed sequence approaches the actual DTFT as the window size increases. Gibb s phenomenon is absent
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 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 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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
IT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete all 2008 or information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. assachusetts
More informationMultimedia Signals and Systems - Audio and Video. Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2
Multimedia Signals and Systems - Audio and Video Signal, Image, Video Processing Review-Introduction, MP3 and MPEG2 Kunio Takaya Electrical and Computer Engineering University of Saskatchewan December
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 informationELEG 305: Digital Signal Processing
ELEG 305: Digital Signal Processing Lecture : Design of Digital IIR Filters (Part I) Kenneth E. Barner Department of Electrical and Computer Engineering University of Delaware Fall 008 K. E. Barner (Univ.
More informationDigital Signal Processing:
Digital Signal Processing: Mathematical and algorithmic manipulation of discretized and quantized or naturally digital signals in order to extract the most relevant and pertinent information that is carried
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 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 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 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 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 informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationDHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A
DHANALAKSHMI COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONICS ENGINEERING EC2314- DIGITAL SIGNAL PROCESSING UNIT I INTRODUCTION PART A Classification of systems : Continuous and Discrete
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 informationChirp Transform for FFT
Chirp Transform for FFT Since the FFT is an implementation of the DFT, it provides a frequency resolution of 2π/N, where N is the length of the input sequence. If this resolution is not sufficient in a
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 informationAnalog vs. discrete signals
Analog vs. discrete signals Continuous-time signals are also known as analog signals because their amplitude is analogous (i.e., proportional) to the physical quantity they represent. Discrete-time signals
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 informationMaximally Flat Lowpass Digital Differentiators
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: 78 60-36 fax: 78
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 informationDiscrete-time Fourier transform (DTFT) representation of DT aperiodic signals Section The (DT) Fourier transform (or spectrum) of x[n] is
Discrete-time Fourier transform (DTFT) representation of DT aperiodic signals Section 5. 3 The (DT) Fourier transform (or spectrum) of x[n] is X ( e jω) = n= x[n]e jωn x[n] can be reconstructed from its
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science. Fall Solutions for Problem Set 2
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science Issued: Tuesday, September 5. 6.: Discrete-Time Signal Processing Fall 5 Solutions for Problem Set Problem.
More information( ) John A. Quinn Lecture. ESE 531: Digital Signal Processing. Lecture Outline. Frequency Response of LTI System. Example: Zero on Real Axis
John A. Quinn Lecture ESE 531: Digital Signal Processing Lec 15: March 21, 2017 Review, Generalized Linear Phase Systems Penn ESE 531 Spring 2017 Khanna Lecture Outline!!! 2 Frequency Response of LTI System
More informationFIR 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 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 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 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 informationChapter 5 Frequency Domain Analysis of Systems
Chapter 5 Frequency Domain Analysis of Systems CT, LTI Systems Consider the following CT LTI system: xt () ht () yt () Assumption: the impulse response h(t) is absolutely integrable, i.e., ht ( ) dt< (this
More informationVALLIAMMAI ENGINEERING COLLEGE. SRM Nagar, Kattankulathur DEPARTMENT OF INFORMATION TECHNOLOGY. Academic Year
VALLIAMMAI ENGINEERING COLLEGE SRM Nagar, Kattankulathur- 603 203 DEPARTMENT OF INFORMATION TECHNOLOGY Academic Year 2016-2017 QUESTION BANK-ODD SEMESTER NAME OF THE SUBJECT SUBJECT CODE SEMESTER YEAR
More informationUNIVERSITY OF OSLO. Please make sure that your copy of the problem set is complete before you attempt to answer anything.
UNIVERSITY OF OSLO Faculty of mathematics and natural sciences Examination in INF3470/4470 Digital signal processing Day of examination: December 9th, 011 Examination hours: 14.30 18.30 This problem set
More informationChap 2. Discrete-Time Signals and Systems
Digital Signal Processing Chap 2. Discrete-Time Signals and Systems Chang-Su Kim Discrete-Time Signals CT Signal DT Signal Representation 0 4 1 1 1 2 3 Functional representation 1, n 1,3 x[ n] 4, n 2 0,
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 informationQuestion Bank. UNIT 1 Part-A
FATIMA MICHAEL COLLEGE OF ENGINEERING & TECHNOLOGY Senkottai Village, Madurai Sivagangai Main Road, Madurai -625 020 An ISO 9001:2008 Certified Institution Question Bank DEPARTMENT OF ELECTRONICS AND COMMUNICATION
More informationLecture 3 - Design of Digital Filters
Lecture 3 - Design of Digital Filters 3.1 Simple filters In the previous lecture we considered the polynomial fit as a case example of designing a smoothing filter. The approximation to an ideal LPF can
More informationResponses of Digital Filters Chapter Intended Learning Outcomes:
Responses of Digital Filters Chapter Intended Learning Outcomes: (i) Understanding the relationships between impulse response, frequency response, difference equation and transfer function in characterizing
More informationDiscrete-Time Signals and Systems. Frequency Domain Analysis of LTI Systems. The Frequency Response Function. The Frequency Response Function
Discrete-Time Signals and s Frequency Domain Analysis of LTI s Dr. Deepa Kundur University of Toronto Reference: Sections 5., 5.2-5.5 of John G. Proakis and Dimitris G. Manolakis, Digital Signal Processing:
More informationOptimum Design of Frequency-Response-Masking Filters Using Convex-Concave Procedure. Haiying Chen
Optimum Design of Frequency-Response-Masking Filters Using Convex-Concave Procedure by Haiying Chen B.Eng., University of Electrical Science and Technology of China, 2013 A Dissertation Submitted in partial
More informationDigital Signal Processing
COMP ENG 4TL4: Digital Signal Processing Notes for Lecture #20 Wednesday, October 22, 2003 6.4 The Phase Response and Distortionless Transmission In most filter applications, the magnitude response H(e
More informationDigital Speech Processing Lecture 10. Short-Time Fourier Analysis Methods - Filter Bank Design
Digital Speech Processing Lecture Short-Time Fourier Analysis Methods - Filter Bank Design Review of STFT j j ˆ m ˆ. X e x[ mw ] [ nˆ m] e nˆ function of nˆ looks like a time sequence function of ˆ looks
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 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 informationDesign of Nonuniform Filter Banks with Frequency Domain Criteria
Blekinge Institute of Technology Research Report No 24:3 Design of Nonuniform Filter Banks with Frequency Domain Criteria Jan Mark de Haan Sven Nordholm Ingvar Claesson School of Engineering Blekinge Institute
More informationMultidimensional digital signal processing
PSfrag replacements Two-dimensional discrete signals N 1 A 2-D discrete signal (also N called a sequence or array) is a function 2 defined over thex(n set 1 of, n 2 ordered ) pairs of integers: y(nx 1,
More informationDiscrete-Time Fourier Transform
C H A P T E R 7 Discrete-Time Fourier Transform In Chapter 3 and Appendix C, we showed that interesting continuous-time waveforms x(t) can be synthesized by summing sinusoids, or complex exponential signals,
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 informationECE-700 Review. Phil Schniter. January 5, x c (t)e jωt dt, x[n]z n, Denoting a transform pair by x[n] X(z), some useful properties are
ECE-7 Review Phil Schniter January 5, 7 ransforms Using x c (t) to denote a continuous-time signal at time t R, Laplace ransform: X c (s) x c (t)e st dt, s C Continuous-ime Fourier ransform (CF): ote that:
More informationLTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Phase Transfer Functions Linear Phase Transfer Functions
LTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Pase Transfer Functions Linear Pase Transfer Functions Tania Stataki 811b t.stataki@imperial.ac.uk Types of Transfer Functions Te time-domain
More informationIT DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A
DEPARTMENT OF ELECTRONICS AND COMMUNICATION ENGINEERING IT6502 - DIGITAL SIGNAL PROCESSING (2013 regulation) UNIT-1 SIGNALS AND SYSTEMS PART-A 1. What is a continuous and discrete time signal? Continuous
More informationComputer Engineering 4TL4: Digital Signal Processing
Computer Engineering 4TL4: Digital Signal Processing Day Class Instructor: Dr. I. C. BRUCE Duration of Examination: 3 Hours McMaster University Final Examination December, 2003 This examination paper includes
More informationFrequency-Domain C/S of LTI Systems
Frequency-Domain C/S of LTI Systems x(n) LTI y(n) LTI: Linear Time-Invariant system h(n), the impulse response of an LTI systems describes the time domain c/s. H(ω), the frequency response describes the
More informationDiscrete-Time Signals and Systems
ECE 46 Lec Viewgraph of 35 Discrete-Time Signals and Systems Sequences: x { x[ n] }, < n
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 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 informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
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 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 information2.161 Signal Processing: Continuous and Discrete Fall 2008
IT OpenCourseWare http://ocw.mit.edu 2.161 Signal Processing: Continuous and Discrete all 2008 or information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. assachusetts
More informationVel Tech High Tech Dr.Ranagarajan Dr.Sakunthala Engineering College Department of ECE
Subject Code: EC6502 Course Code:C302 Course Name: PRINCIPLES OF DIGITAL SIGNAL PROCESSING L-3 : T-1 : P-0 : Credits 4 COURSE OBJECTIVES: 1. To learn discrete Fourier transform and its properties 2. To
More informationOptimum Ordering and Pole-Zero Pairing of the Cascade Form IIR. Digital Filter
Optimum Ordering and Pole-Zero Pairing of the Cascade Form IIR Digital Filter There are many possible cascade realiations of a higher order IIR transfer function obtained by different pole-ero pairings
More informationEEL3135: Homework #4
EEL335: Homework #4 Problem : For each of the systems below, determine whether or not the system is () linear, () time-invariant, and (3) causal: (a) (b) (c) xn [ ] cos( 04πn) (d) xn [ ] xn [ ] xn [ 5]
More informationDiscrete Time Signals and Systems Time-frequency Analysis. Gloria Menegaz
Discrete Time Signals and Systems Time-frequency Analysis Gloria Menegaz Time-frequency Analysis Fourier transform (1D and 2D) Reference textbook: Discrete time signal processing, A.W. Oppenheim and R.W.
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 informationDiscrete-Time Fourier Transform (DTFT)
Discrete-Time Fourier Transform (DTFT) 1 Preliminaries Definition: The Discrete-Time Fourier Transform (DTFT) of a signal x[n] is defined to be X(e jω ) x[n]e jωn. (1) In other words, the DTFT of x[n]
More information