DSP-Mitra-10FIR Digital Filter.pptVIP

284 10.3 CAD of Equiripple Linear-Phase FIR Filters (II) Weighted error function: (10.47) or (10.62) (10.68) 10.3 CAD of Equiripple Linear-Phase FIR Filters (III) Chebyshev or Minimax criterion: equiripple FIR filter Minimize the peak absolute value of Linear-phase FIR filters obtained by the criterion polynomial approximation 10.3 CAD of Equiripple Linear-Phase FIR Filters (IV) Alternation Theorem: Let R be a union of disjoint closed subsets of ? Let a desired function D(x) and weighting function W(x) be continuous on R Define the error function E (x) = W (x) [PL(x)－D(x)] Maximum error 10.3 CAD of Equiripple Linear-Phase FIR Filters (V) necessary and sufficient condition for PL(x) being the unique Lth order polynomial under the Minimax criterion can be expressed by the alternation theorem: E(x) has at least L + 2 alterations on F , i.e. ? xi , i = 1, . . . , L’ ≥L + 2 such that xi xi+1, E(xi) = － E(xi+1), for i = 1, . . . , L’－1 and E(xi) = ±Emax, for i = 1, . . . , L’ 10.3 CAD of Equiripple Linear-Phase FIR Filters (VI) Parks-McClellan Algorithm Iterative method to determine the alternation frequencies ωi and the ripple 1. initialize ωi to some values 2. compute — the ripple corresponding to the alternation frequencies 3. interpolate a polynomial between the alternation points 4. find the maximum/minimum values of the error 5. if |E(ω)| ≤ : stop else compute new ωi’ as extreme of E(ω), and go to 2 (else recursive) 10.5 FIR Digital Filter Design Using Matlab Order Estimation: kaiord() Kaiser’s Formula bellangord() Bellanger’s Formula remezord() Hermann’s Formula kaiserord() filter order for Kaiser window-based design 10.5 FIR Digital Filter Design Using Matlab(II) Equiripple linear-phase FIR filter design: remez() equiripple FIR filter design using Parks-McClellan a