0% found this document useful (0 votes)
31 views100 pages

Fir Part 1

Uploaded by

jdn6qkm3l1
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF or read online on Scribd
0% found this document useful (0 votes)
31 views100 pages

Fir Part 1

Uploaded by

jdn6qkm3l1
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF or read online on Scribd
You are on page 1/ 100
> ? Finite Impulse Response Filters pve. Introduction —- Inthe previous chapter \ digital processing Applica ‘The following are the ma < Stalled the design of IIR filter and its properties, tn many Nions BIR filters are proferted over their IIR. counterparts Nn advantages of the FIR filter over HR sitter, 1, FIR filters are Always stable, 2, FIR filters with exactly linear phase ¢ 3, FIR filters can be rey 4, FIR filters are tree of Himit eyele word length digital system, 5. Excellent design Methods are n easily be designed lized in both rec Isive and Non-recursive structures, oscillations, when implemented on a finite Available for various kinds of FIR filters, The disadvantages of FIR filter are 1. The implementation of narrow transition band FIR filters are Very costly, as it requires considerably more arithmetic operations and hantware vomponents Such as multipliers, adders and delay elements, % Memory requirement and execution time are very high, 62 Linear phase FIR filter The transfer function of a FIR causal filter is given by H(z) Ss h(n)sr" (1) nad Where A(n) is the impulse response of the filter, Fourier transform of h(n) is Nal x > (nent (6.2 ned H(e!) 6.2 Digital Signal Processing which is periodic in frequency with period 27, H(e™) = £|H(e%)|¢i(w) where |#7(e%)| is magnitude response and 8(.2) is phase respon r We define the phase delay and group delay of a filter ay = 9) ~d9(w) w diy and Tp For FIR filters with linear phase we can define 6(w)=-aw —7 Rin)e=dor N x (Jerson 4 ( MN ya +> aw—1_ m)e~#(N=1~m mo : ~~62) s ~ Drolet en (Nw) tae aX ) th oy eS4(N-1)/2 + DAN ~ 1 — n)e-de(w-1-n) » nye 6%) Fora sy Substig mn tettical impulse respo: mu Neale ‘tituting this relation in Eq (ox) a ‘al n= 2 () = = SY hinjenson +h ( } e012 4. (ner fa(N—1-m) no DY p(njessler-n2- no my +P hlnperFl D-H ( n=0 an Ay =i Neal eg) [E 2h(n) cosw = =n) +h( | ; a0 6.6 Digital Signal Processing Let 431 — n= p,then N=} -y2l< Eel H(e) = e FN -V/2 [x 2h (? i ”) coswp +h p= = ¢S(N-1)/2 an (> =n} eosun ss (SH) ral B 2 Noi = eF4N-1/2S > a(n) coswn (6m, 7=0 where a(o) =» (4) a(n) = 2h (SS -n) We can write Eq. (6.24) as H (4) = eH -D AF ( 54) = H(i) a where os) H(e#) = > a(n) coswn 7=0 and Nod > a Ow) = -aw = - ( } ” ‘ , so of Hl H(e) is a real and even function of w. The magnitude and pr \4(@)| = H(e)| and LHe) = 0(w) = aw when Hel) 2 * ZH(E*) = ~aw + when He”) a and (w) = aw when Hie) < fron i stinguish 3 Hi (©) is called as zero-phase frequency response (0 is my a ai nitude response, The zero-phase frequency response of the INT y oat gl tive and negative values, whereas the magnitude response °°) ry odd Mot! ‘The frequency response of symmetric impulse resPO"S Ut (ei) 8 ‘erence between | (e)) and H (e!) and between’) Fig. 63. 3.83. Relation between magnitude response |H(e*)] and the zero phase ssmse H(e*) and between ZH (el) and 0(). Gell: Symmetric impulse response for N even Ne He) = za (men N-1 = Jenson 4 aM A(n)e" n=0 “SA (n)e-0" 4. Si (NV = 1 = nee =1-m) n=0 Nkoow nin) =hN-+1-n) Hele) as Mea )s eS Anerson 4. = A(njenietN-1-0) an nad on S e-ion-1)/9 [x h(n)ellO¥-D/2-M1 4.7 (njer#—DPn) n=0 n=0 i Sersuwity. |e 1S an(n) contd (% ned tal signal processin8 68 Digi i, Eo, (M a) coe (n-a el * (Ge n N/2 = ten? YH”) cos(n- 1/2) where N 5. b(n) = Qh ( 2 = n) (629) Weean write Eq. (6.28) 2° s H(e) = ese) = He) where a Ale) = Leno G @ al . Ned 3) oe ae l ise respon ‘Te frequency esponse of linear phase fiter with jmmessic impu! _ for N even is shown in Fig. 6.4. te) ot? Fig.6.4 Freq ve response fr linear phase FIR fiter, symmet"° impuls? Case I: Antis . symmetric For this type of asia odd nA) -6 = Nz1 ; H(e) = SS bine” i = ole : Finite Impulse Response Filters 6.9 m Be ‘ N=1 Exe +h () enFa(N-1)/2 se Ae ” nati N-1 SSH Dyer =) n= NH 7 », SO a = 1 = njeel=1-n) n=0 Sa pee =yhlnye —1~n), therefore, He) = Y. h(nje™" — 57 h(njenMIN-1-n) n= =e HN-D/2 = h(n)efelN-D)/2—] _ = h(njerdell- n=0 n=0 5 ==) =e tN-V2ei#/2 \ ~ e(n) sinwn (6.2) ni Sem 77 = 2h(n) sinw ( a0 wa = e-d(N-1)/2,5n/2 SD Pr eC : St (6.33) (6.34) ; - (* (6.35) Yay Ue : : Vea J Tesponse of linear phase FIR filter for antisymmetric sequence ‘wn in Fig, 65. jun a et = F pment ¥ h(nje n=N/2 —jun Xs a See ee 2 os wenave s(p) Set Cieatean eo H(&*) fs h(n)je Je" = Sm hi oa” ] n=0 ‘| = en fetN-1/2) = ninjeielln VAM - H a(nye nl on Nan n=0 NI: = enHulN~1)/2gin/2 ba » (F-n) sinw(n - va = e-bo(N-U)/2gin/2 [pan Bn é . 4] ¢ 2 where d(n) = 2h (F-») eee: =e IHN-1)/2¢5n/2 (eH) janes amcor ‘gs Frequency response of linear phase fiter for antisymmetric impulse re- od se with NV even. np Ale) = So d(n)sinu(n= 1/2) (6.40) nal 6(w) ; “ao - (6.41) The frequency response of linear phase filter for antisymmetric sequence with Nv enis shown in Fig. 6.6. The impulse response of symmetric with odd number of samples can be used to ésign all types of filters. = The symmetric impulse response having even number of samples cannot be used indesign highpass filters. The frequency response of antisymmetric impulse response is imaginary and this {prof filters are most suitable for such filters as Hilbert transformers and differen- tiaors, ample 6.1 Determine the frequency response of FIR filter defined by y(n) = Mie(n) + 2(n — 1) + 0.252(n — 2). Calculate the phase delay and group delay. Suton Given, E y(n) = 0.252(n) + a(n ~ 1) + 0.25e(n = 2) vg Fourier transform on both sides ¥(e*) = 0,25 (el) + e-J4X (elt) 0.25e72 X (el) Hid» _ ¥ (ei) X (ei = 0:25 + 7 4 0,25e7 24 e4* (0,253 + 1 + 0.25e%) = e™(1 + 0,508) eH (3%) : (41a) "BE. (6.41a) with Eq, (6.25) we get 0(w) = —wr Sedelay 7, — ~O(w) _ w os e: Te "OUD delay A causal TI system have linear phase iis impulse response h(n) We orantisymmetrical about n = Y—* satisfying the condition h(n) = ACN 1. rte peeavely The Given inpulsc reqeee ofl Therefore for linear phase the impulse response is either symmetrical or antisymnex rical about n = 2. That is, for symmetrical impulse response ag = 4630, = a5;07 = a4 For antisymmetrical impulse response Qo = —06; a) = —45; a2 = —a43a3 = 0 Example 6.3 The impulse response of a non-causal LTI-FIR discrete time systems given by : A(n) = ap5(n + 2) + 014(n + 1) + a26(n) + a35(n ~ 1) + ax5(n-2) Find the values of the impulse response samples for which its frequency "570% H(e*) have a zero phase. Solution 4 i ‘A non-causal LTL-FIR system have zerophase if its impulse response) ‘symmetrical or antisymmetrical about n = 0. the impul™ The given impulse response is of length 5. Therefore for zerophas response is symmetrical or antisymmetrical about n = 0. ‘That is, for symmetrical impulse response ao = 44; 41 “i =a For antisymmetrical impulse response a2 = 0; ag = — 41) 44 = o give ters zie Problem 6.1 The frequency response of a linearphase FIR y 3 bw H(e) = cos 5 + 5008 - Determine the impulse response h(n). x0 ET RS in ier Practice Problem 6.2 If the frequency response of a linear-phase by H(e) = €-34(9,30 + 0.5 cosw + 0.3008) Determine filter coefficients. : Naa H(2)= Do h(n)e- n=0 “zowsanfinite ero of F(z), then . N-1 @)|__ = H(z) = Yo A(n)aq" =0 —- n=0 =A(0) + A(1)2g1 +... 4+ A(N - 1) =0 Fora linear phase filter h(n) = h(N -1=n) . (IN = 1) ACN — 2)251 +... + A(1)zp"— + h(0)25 "9 =0 IND Th(N — 1)2f"! + A(N = 2)z9/-? +... + h(1)29 + A(O)] =O Fig. 6.7 Zero locations of linear phase FIR filter. Not = 2") ninja =0 H(z) =) piayiagty™ = 2 WD H(251) =0 n=0 Ho) = (ze) ar (6.44) AE Ere («- S) =o te o(E-s) m09'E acme? | : uan9 NV “gsuodsaz osindust| : rounouusks| m Sa ‘ See i [ ( (n@)H 234 = 0 Be al fe © 0 < (wr?) 304 0 = 9 204m gee eee fumes (uo T< | aya yee? po gt m— re ay. |’ asuodsau osindutt pomautusxs| (n®)H7 (2) suoneoyddy ‘asuodsas ast asuodsar opmruseyy asuodsax Aouonbasy adh Soy Ula aseyd TwaUy Jo sonsuaoemyD Jo AEUNUNG 1-9 S14BL («- =) = (u)P 0> (m2)H 105 4 = 9 14) Ire 1s taxa O<(oA)H310=02%4N | |(2—u) mus(up<| |(E-") mms ~X ey -pyne- Pear? N ssuodsar asindu 0+ S40 i & * oer, sia 2 (ene 0>(,9)H 10} 4=9 rol oy a o< ecu agus (ue Be SG zht-n)oe—e/2 | y ee Paar a: a Pomeuurtsnuy | (op)H7 Moe) (ode | ] suodsat ose asuodsai apmnudoyy 2suodsai sottonbany oak SOMOS 19 TAR, — 6.16 Digital Signal Processing From Eq. (6.44) we find that if zp is a zero of A(z), then zal From the above result we can find the location of zeros in a linear, as follows (See Fig. 6.7). 1. If z; = —1, then 271 = 2, the zero lies at z one zero on the unit circle. s . If zz is real zero with.|2o| < 1 then 25+ zeros in this group. 3. Ifzs is a complex zero with |z3| = 1 then 251 = 2* and there ae toa this group. . = =1, this group is also a real zero and there ar 4. If zis complex zero with |z4| # 1. This group contains four zeros 23712 23) Example 6.4 If H(z) has zeros at z) = 1/24 5/V2, z2 = 2Determine the long degree H(z) that has a linear phase. Solution If zp = 2 isa zero fora linear Phase filter, then z; = 1/2is alsoazem Har = 5 + 4, whose |2| = lisa zero, then 2, = 2} is also azer0, The total number of zeros are four zy 1 mome-aen(cde)b-()-Cl = ( = Te =?) ( ~ V2 +27) Practice Problem 63 If A(. lowest degree H(z) that has e 2) has zeros at $0-3"/? and 26-0", dee a linear phase. 65 The Fourier Series Method ae? mee Fsponse (©) of a system is periodic in 27. From pc cata Bett funtion can be expressed a ie epee of Designing FIR Filtes é Hae?) = DX haln)e-*" a where the Fourier : or a the filter “SPs han) ae the desired impulse P= é Es han) = = ip Ha)" du she Eg. (6.47) represents a non-causal digital fiter of i ppercanster function, the series can be truncated by: A(n) = ha(n) for |n| < x =0 otherwise 64) HR)= >) - Anz” 6% shi 2) OM 4 4 hs +0) +h Det Qe +4[- SJ)“ a =h0) + > [h(n)z-™ + h(-n)z"] 650) nt For a symmetrical impulse response having symmetry at 1 = 0 h(—n) = h(n) 651 Therefore, Eq. (6.50) can be writen as wet H(z) = h(0) + > a(n" +27) 652) n=l Me shove transfer finction is not physically realizable Realizability can be by multiplying Eg, (6.52) by 2-(—)/ where ~ H(z) = DPA (2) a mse a2 fx + Daye a 653) n=l — 6.18 Digital Signal Processing Example 6.5 Design an ideal lowpass filter with a frequency respon se Hale) = 1 for = 0 for \ Find the values of h(n) for NV = 11. Find H(z). Plot the magnitude reg ‘ Ponse, Solution The frequency response-iowpass filter with w, = § is shown in Fj {Given 1B 68 Ha(e”) = 1 for 7 = Ofor= <|w) qn) (+27")] ol 5 ss BPS aay ("+427") = S05 Neg E881 + 274) — 0.106 (22 42-4) + 0.06986 (2° + 2-8), “inetion ofthe realizable filter is H"(z) = 2-(N-D/2 7/2) =o oe + 0.063662 From the above Eq. (6.57) the fier coefficients of (0) = h(10) = 9.06366; n(1) = h(9) = 05 AQ) = a(3) = A(7) (4) = h(6) = 0.3183; M6) =05 ‘The frequency response is given by om He) = Do a(n) coswn_-where n=0 oo) =n (Se+) =A@)=08 a(n) = 2h (S -n) a(t) = 2h(5 — 1) = 2h(4) = 0.6366 a(2) = 2h(5 — 2) = 2h(8) = a(3) = 2h(5 ~ 3) = 2h(2) = -0.212 (4) : (5) = 2h(5 ~ 4) = 2h(1) = a(5) = 2h(5 — 5) = 2h(0) = 0.127 a4 Fc) = 0.5 + 0.6366 cos w — 0.212 cos 3u + 0.127 60s: ‘The magnitude in dB is calculated by varying w from 0 to The magnitude | (e)|ag = 20 log |H(e”)}. 04 021 -0.26-0.517-0.21 0.42 0.77 9 100 110 120 130 140 150 16 6 -14,56-31,89-206 -26 32 -24.7-3055 ney response plot is shown in Fig. 6:9, pete Fig.6.9 Frequency response of lowpass fier of example 6.5. ae i Fame 66 Design an ideal highpass filter with a frequency response Hya(e*) = 1 for 7
pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy