The design of a FIR filter starts with its specifications in either discrete-time domain or DTFT frequency domain, or both. In the time domain, the design objective is the impulse response. In the frequency domain, the requirement is on various parameters of the magnitude response, shown in Fig.5.8 for a lowpass filter. Important parameters are band edge frequencies ωp size 12{ω rSub { size 8{p} } } {} and ωs size 12{ω rSub { size 8{s} } } {}, passband

Typical normalized magnitude response specfications of a lowpass filter [0, ω p ] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqaaeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaacUfacaaIWaGaaiilaiabeM8a3naaBaaaleaacaWGWbaabeaakiaac2faaaa@3C03@ , transition band Δω=ωs−ωp size 12{Δω=ω rSub { size 8{s} } - ω rSub { size 8{p} } } {}, stopband , cutoff frequency ωc=(ωp+ωc)/2 size 12{ω rSub { size 8{c} } = \( ω rSub { size 8{p} } +ω rSub { size 8{c} } \) /2} {}, passband ripple (or deviation) δp size 12{δ rSub { size 8{p} } } {} and stopband ripple δs size 12{δ rSub { size 8{s} } } {}. The two ripples are usually assumed equal . Notice that for ideal (desired) filter, the passband frequency magnitude is normalized to 1 and the stopband to 0 , and the frequency response of the designed filter oscillates between the high amplitude of 1 or the low amplitude of 0. illustrates the specifications of bandpass filter. Here there are two sets of edge frequencies, lower and upper. The rising bandwith Δ ω l MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipDI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqadeaabiGaaiaacaqaaeaadaqaaqaaaOqaaiabfs5aejabeM8a3naaBaaaleaacaWGSbaabeaaaaa@3A52@ and the falling bandwidth Δ ω u MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipDI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqadeaabiGaaiaacaqaaeaadaqaaqaaaOqaaiabfs5aejabeM8a3naaBaaaleaacaWG1baabeaaaaa@3A5B@ are usually assumed equal. Even the initial requirement is imposed on magnitude response but we must design the filter having linear phase or generalized linear phase. Overall filter design From the filter specifications, the first step is to choose between FIR and IIR filters based on their advantages and disadvantages. This chapter only concerns FIR filters. The next step is to select the proper linear phase FIR filters (section 5.2.3). A rather complete procedure for the design of FIR filters are as follows. Specifications of the filter Choosing an appropriate linear phase filter type (section 5.2.3) Choosing the mothod of design such as window, optimal, frequency sampling… Calculation the filter coefficients (impulse response) Finding suitable structure Analysis of the finite wordlength effects (chapter 7) Implementation of the filter in hardware and/or software (chapter 7).Fig.5.9: Bandpass filter specifications

Bandpass filter specifications

Fixed windows The impulse response of ideal filters is infinite duration (IIR). We cannot evaluate the corresponding frequency response and, especially, implement the filter by hardware and/or software. Thus we must truncate the impulse response at both ends with respect to the central. Even we truncate the impulse response when it is small enough, but such a sudden cutoff will cause some undesired effects. The window method will reduce them. In the time domain, windowing means that the we multiply the desired (usually ideal) infinite duration impulse response hd(n) size 12{h rSub { size 8{d} } \( n \) } {} by a finite duration window (or window function) w(n) size 12{w \( n \) } {} to get a soft truncation. The resulted impulse response h(n) size 12{h \( n \) } {} of the designed filter is the product h(n)= h d (n)w(n) , 0≤n≤M MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeqabiGaciaacaqaaeaadaqaaqaaaOqaaiaadIgacaGGOaGaamOBaiaacMcacqGH9aqpcaWGObWaaSbaaSqaaiaadsgaaeqaaOGaaiikaiaad6gacaGGPaGaam4DaiaacIcacaWGUbGaaiykaGqabiaa=bcacaGGSaGaa8hiaiaa=bcacaWFGaGaa8hiaiaa=bcacaWFGaGaaGimaiabgsMiJkaad6gacqGHKjYOcaWGnbaaaa@4CCC@ Here we assume that all desired impulse responses and windows are causal, from n=0 size 12{n=0} {} to n=M size 12{n=M} {}, i.e. window length is M+1 size 12{M+1} {} samples (time indices). Many authors start the design with hd(n) size 12{h rSub { size 8{d} } \( n \) } {} and w(n) size 12{w \( n \) } {} bi-sided (noncausal), i.e. defined in the interval then shift to the right M/2 size 12{M/2} {} indices to make them causal. Multiplying in time domain corresponds to convolution in frequency domain. Thus the frequency response of the designed filter (corresponding to the windowed impulse response h(n) size 12{h \( n \) } {}) is H(ω)=Hd(ω)∗W(ω)=12π∫−ππHd(ω')W(ω−ω')dω' size 12{H \( ω \) =H rSub { size 8{d} } \( ω \) * W \( ω \) = { {1} over {2π} } Int rSub { size 8{ - π} } rSup { size 8{π} } {H rSub { size 8{d} } \( ω' \) W \( ω - ω' \) dω'} } {} where W(ω) size 12{W \( ω \) } {} is the Fourier transform (DTFT) of the window w(n) size 12{w \( n \) } {}. In the window design method , we first evaluate the desired filter impulse response h d (n) MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqadeaabiGaaiaacaqaaeaadaqaaqaaaOqaaiaadIgadaWgaaWcbaGaamizaaqabaGccaGGOaGaamOBaiaacMcaaaa@3A3B@ from the given desired frequency response H d (ω) MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipDI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeqabiGaciaacaqaaeaadaqaaqaaaOqaaiaadIeadaWgaaWcbaGaamizaaqabaGccaGGOaGaeqyYdCNaaiykaaaa@3B16@ , and then apply an appropriate window . Thus the method should be called the Fourier – window method. Retangular window Actually, the sudden truncation mentioned earlier is the simplest window called rectangular window, defined as ()

Rectangular window of order M or length M + 1 w(n)=1,0≤n≤M0,otherwisealignl { stack { size 12{w \( n \) = matrix { {} # 1{} } matrix { {} } , matrix { {} # {} } 0 <= n <= M} {} # size 12{ matrix { {} # {} # matrix { {} # 0{} } {} } matrix { {} } , matrix { {} # {} } "otherwise"} {} } } {} The Fourier transform of this window is W(ω)= ∑ n=−∞ ∞ w(n) e −jωn = ∑ n=0 M e −jωn = 1− e −jω(M+1) 1− e −jω = e −j M 2 ω sin⁡ω(M+1)/2 sin⁡ω/2 , ω≠0 M+1 , ω=0 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeGabiGaciaacaqaaeaadaqaaqaaaOabaeqabaGaam4vaiaacIcacqaHjpWDcaGGPaGaeyypa0ZaaabCaeaacaWG3bGaaiikaiaad6gacaGGPaGaamyzamaaCaaaleqabaGaeyOeI0IaamOAaiabeM8a3jaad6gaaaGccqGH9aqpdaaeWbqaaiaadwgadaahaaWcbeqaaiabgkHiTiaadQgacqaHjpWDcaWGUbaaaaqaaiaad6gacqGH9aqpcaaIWaaabaGaamytaaqdcqGHris5aOGaeyypa0ZaaSaaaeaacaaIXaGaeyOeI0IaamyzamaaCaaaleqabaGaeyOeI0IaamOAaiabeM8a3jaacIcacaWGnbGaey4kaSIaaGymaiaacMcaaaaakeaacaaIXaGaeyOeI0IaamyzamaaCaaaleqabaGaeyOeI0IaamOAaiabeM8a3baaaaaabaGaamOBaiabg2da9iabgkHiTiabg6HiLcqaaiabg6HiLcqdcqGHris5aaGcbaqbaeqabeGaaaqaaaqaauaabeqabiaaaeaafaqabeqadaaabaqbaeqabeGaaaqaauaabeqabiaaaeaaaeaaaaaabaaaaaqaaaqaaaaaaeaafaqabeqacaaabaqbaeqabeGaaaqaauaabeqabiaaaeaaaeaaaaaabaaaaaqaaaaaaaaaauaabeqabeaaaeaacqGH9aqpcaWGLbWaaWbaaSqabeaacqGHsislcaWGQbWaaSaaaeaacaWGnbaabaGaaGOmaaaacqaHjpWDaaGcdaWcaaqaaiGacohacaGGPbGaaiOBaiabeM8a3jaacIcacaWGnbGaey4kaSIaaGymaiaacMcacaGGVaGaaGOmaaqaaiGacohacaGGPbGaaiOBaiabeM8a3jaac+cacaaIYaaaaaaafaqabeqabaaabaaaaiaacYcacaaMf8UaeqyYdCNaeyiyIKRaaGimaaqaauaabeqabiaaaeaaaeaafaqabeqacaaabaqbaeqabeWaaaqaauaabeqabiaaaeaafaqabeqacaaabaaabaaaaaqaaaaaaeaaaeaaaaaabaqbaeqabeWaaaqaauaabeqabiaaaeaaaeaaaaaabaaabaaaaaaaaaGaaGzbVlaaywW7caaMf8UaamytaiabgUcaRiaaigdafaqabeqabaaabaaaaiaacYcafaqabeqacaaabaaabaGaaGzbVlaaywW7aaGaaGzbVlaaywW7caaMe8UaeqyYdCNaeyypa0JaaGimaaaaaa@9C45@ This is the same as Equation (3.54) when replacing M size 12{M} {} by 2N size 12{2N} {}. Also, instead of the above expression we can use the sum of cosines (Equations (3.46), (3.53)). W(ω) size 12{W \( ω \) } {} has includes a phase factor showing the time shift of a two-sided symmetric window to a causal window. The magnitude and phase responses are, respectively,

Normalized frequency response of rectangular window for M = 30 and 50 | W(ω) |=| sin⁡ω(M+1)/2 sin⁡ω/2 | , ω≠0 M+1 , ω=0 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqadeaabiGaciaacaqabeaadaqaaqaaaOabaeqabaWaaqWaaeaacaWGxbGaaiikaiabeM8a3jaacMcaaiaawEa7caGLiWoacqGH9aqpdaabdaqaamaalaaabaGaci4CaiaacMgacaGGUbGaeqyYdCNaaiikaiaad2eacqGHRaWkcaaIXaGaaiykaiaac+cacaaIYaaabaGaci4CaiaacMgacaGGUbGaeqyYdCNaai4laiaaikdaaaaacaGLhWUaayjcSdGaaGjbVlaacYcacaaMf8UaaGzbVlabeM8a3jabgcMi5kaaicdaaeaacaaMf8UaaGzbVlaaywW7caaMf8EbaeqabeGaaaqaaiaad2eacqGHRaWkcaaIXaaabaGaaiilaiaaywW7caaMf8UaaGzbVdaacaaMf8UaaGzbVlabeM8a3jabg2da9iaaicdaaaaa@6F96@ Φ(ω)=− M 2 ω+β MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeGabiGaciaacaqaaeaadaqaaqaaaOqaaiabfA6agjaacIcacqaHjpWDcaGGPaGaeyypa0JaeyOeI0YaaSaaaeaacaWGnbaabaGaaGOmaaaacqaHjpWDcqGHRaWkcqaHYoGyaaa@4265@ where β size 12{β} {} is 0 or π size 12{π} {} as discussed in previous section. Thus the window is linear phase. Fig.3.27 illustrates the variation of the amplitude of ∣W(ω)∣ size 12{ lline W \( ω \) rline } {} (remember to replace 2N size 12{2N} {} by M size 12{M} {}). It is maximum and equals to M+1 size 12{M+1} {} at ω=0 size 12{ω=0} {}. The zero-crossing points are multiples of 2π/(M+1) size 12{2π/ \( M+1 \) } {}. The response consists of a main lobe and many sidelobes. As M size 12{M} {} gets larger the main-lobe gets smaller, and there are more sidelobes, also narrower but the amplitude of the first sidelobe remains the same (Fig.5.11). As M→∞ size 12{M rightarrow infinity } {}, W(ω) size 12{W \( ω \) } {} tends to be an unit sample δ(ω) size 12{δ \( ω \) } {}.

Magnitude response of designed lowpass filter when using rectangular window with M = 44 The oscillatory ∣W(ω)∣ size 12{ lline W \( ω \) rline } {} when convolved with the ideal lowpass response Hd(ω) size 12{H rSub { size 8{d} } \( ω \) } {} will result in a response H(ω) size 12{H \( ω \) } {} having a non-zero transition width and ripples is both the passband and stopband. Fig.5.12 shows the designed lowpass filter when using rectangular window having M=44 size 12{M="44"} {} to abruptly truncate the impulse response of an ideal lowpass filter with cutoff frequency ωc=π/2 size 12{w rSub { size 8{c} } =π/2} {}. Fig.5.13 shows the magnitude in linear scale for two case of filter order, M = 10 and M = 18. In order that H(ω) size 12{H \( ω \) } {} approaches Hd(ω) size 12{H rSub { size 8{d} } \( ω \) } {} the rectangular window must be of infinite duration, i.e. we must take into account the whole impulse response without truncation. In frequency domain this means W(ω) size 12{W \( ω \) } {} is an unit sample δ(ω) size 12{δ \( ω \) } {} as said. We can reason in the reverse direction: Since H(ω)=Hd(ω)∗δ(ω)=Hd(ω) size 12{H \( ω \) =H rSub { size 8{d} } \( ω \) * δ \( ω \) =H rSub { size 8{d} } \( ω \) } {} then δ(ω) size 12{δ \( ω \) } {} must be a unit sample and the windows must be of infinite duration. When we simulate H(ω) size 12{H \( ω \) } {}on a computer using Maltlab sofware, or else, with M size 12{M} {} in the hundreds we will still see the non-zero transition width and ripples due to the Gibbs phenomenon (section 3.1.4). An infinitely long windows is not practical, so the idea is to look for finite duration windows which perform better than the rectangular. Actually in examples 5.2.1 and 5.2.2 the truncation of the impulse responses of the ideal filters meant the rectangular window had been applied.

Magnitude response of designed filter for M = 10 and M = 18 Other windows The Gibbs phenomenon can be reduced considerably by the use of a softer (less abrupt) truncation, i.e. by tapering the rectangular smootly to zero at both ends. Unfortunately, the reduction in amplitude of sidelobes of H(ω) size 12{H \( ω \) } {}is accompanied by a widening of its mainlobe. For comparision, the window responses W(ω) size 12{W \( ω \) } {}are usually plotted in dB size 12{ ital "dB"} {} scale.

Triangular window (Bartlett window) Several smooth window functions have been proposed and used. It might come to our mind that the first applicant would be the triangular window, also called Bartlett window , depicted in Fig.5.14. From the figure we can write the window function ( table 5.1). Table 5.1 lists common fixed windows and Fig.5.15 plots their functions. Fig.5.16 plots the dB magnitude for and Fig.5.17 for M = 50.The first sidelobe is at -25dB compared to the -13dB of the corresponding (same length) rectangular window. However the mainlobe width is about twice as large.

Common window plotted as function of continuous time n (whereas actual windows are function of discrete index n) The sidelobe of the triangular window is still high because the tapering is still rather coarse. For smoother tapering, cosinusoid is incorporated into the window function. This observation has led to the three well known windows: Hanning (or von Hann), Hamming, and Blackman, all defined in the interval 0≤n≤M size 12{0 <= n <= M} {}, otherwise zero (table 5.1). We can check that the functions are normalized (peak value of 1 at n=M/2 size 12{n=M/2} {}) and zero at both ends except the rectangular and the Hamming (0.08 instead of 0). Besides, there are many other less-used fixed windows. Notice that all windows mentioned, from rectangular to Blackman, are simple functions easily evaluated, and their frequency responses concentrate around ω=0 size 12{ω=0} {}(the mainlobe) as expected. Also notice that all the windows mentioned are symmetric about the mid-point n=M/2 size 12{n=M/2} {}, this when combined with the symmetry or antisymmetry of the filter’s impulse response will make the corresponding designed filter linear phase or generalized linear phase.

dB magnitude frequency response of fixed windows for M = 20 (a) Bartlett, (b) Hanning, (c) Hamming, (d) Blackman

dB magnitude frequency response of various windows for M = 50 (a) Bartlett, (b) Hanning, (c) Hamming, (d) Blackman Table 5.2 lists various features of common windows for comparison. The passband ripple and the stopband ripple δs size 12{δ rSub { size 8{s} } } {} and the stopband ripple δ s MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqadeaabiGaciaacaqabeaadaqaaqaaaOqaaiabes7aKnaaBaaaleaacaWGZbaabeaaaaa@38AE@ , are taken equal (the smaller of the two): A=−20 log⁡ 10 [min⁡( δ p , δ s )] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaadgeacqGH9aqpcqGHsislcaaIYaGaaGimaiGacYgacaGGVbGaai4zamaaBaaaleaacaaIXaGaaGimaaqabaGccaGGBbGaciyBaiaacMgacaGGUbGaaiikaiabes7aKnaaBaaaleaacaWGWbaabeaakiaacYcacqaH0oazdaWgaaWcbaGaam4CaaqabaGccaGGPaGaaiyxaaaa@4ACA@ Fig.5.16, Fig.5.17 and Table 5.2 show that no window is the best in all aspects but there is a tradeoff of features, and the choice of an appropriate window depending on our requirement. For smallest mainlobe width it’s the rectangular window, for best sidelobe attenuation it’s the Blackman. In between, the Hamming is a good choice. The Bartlett is the transition from the rectangular to the other windows. See section 5.3.3 for design examples.

The Kaiser window The windows discussed so far are fixed windows, only the length (M+1) size 12{ \( M+1 \) } {} is adjustable. The Kaiser window has an additional ripple parameter β size 12{β} {} , enabling the designer to tradeoff the transition and ripple. It is defined in the interval 0≤n≤M size 12{0 <= n <= M} {}, otherwise zero: w(n)= I 0 { 1−( n−M/2 M/2 ) β 2 } I 0 (β) , 0≤n≤M MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaadEhacaGGOaGaamOBaiaacMcacqGH9aqpdaWcaaqaaiaadMeadaWgaaWcbaGaaGimaaqabaGcdaGadaqaamaakeaabaGaaGymaiabgkHiTmaabmaabaWaaSaaaeaacaWGUbGaeyOeI0Iaamytaiaac+cacaaIYaaabaGaamytaiaac+cacaaIYaaaaaGaayjkaiaawMcaaaWcbaqcLbyacqaHYoGycaaMe8oaaOqbaeqabeqaaaqaamaaxacabaaaleqabaGaaGOmaaaaaaaakiaawUhacaGL9baaaeaacaWGjbWaaSbaaSqaaiaaicdaaeqaaOGaaiikaiabek7aIjaacMcaaaqbaeqabeqaaaqaaaaacaGGSaqbaeqabeGaaaqaaaqaaaaacaaIWaGaeyizImQaamOBaiabgsMiJkaad2eaaaa@5897@ where I0(x) size 12{I rSub { size 8{0} } \( x \) } {} is the modified zeroth-order Bessel function, computed by a power series expansion: I0(x)=1+∑k=1∞x(2)kk!2=1+x24+122⋅x24⋅x24+1z2⋅x24⋅x416+⋅ size 12{I rSub { size 8{0} } \( x \) =1+ Sum cSub { size 8{k=1} } cSup { size 8{ infinity } } { left [ { {x \( 2 \) rSup { size 8{k} } } over {k!} } right ]} rSup { size 8{2} } =1+ { {x rSup { size 8{2} } } over {4} } + { {1} over {2 rSup { size 8{2} } } } cdot { {x rSup { size 8{2} } } over {4} } cdot { {x rSup { size 8{2} } } over {4} } + { {1} over {z rSup { size 8{2} } } } cdot { {x rSup { size 8{2} } } over {4} } cdot { {x rSup { size 8{4} } } over {"16"} } + cdot cdot cdot } {} Following is the values of I0(x) size 12{I rSub { size 8{0} } \( x \) } {} for small x size 12{x} {}: X I 0 (n) X I 0 (n) 0.0 1.000 1.2 1.3938 0.1 1.002 2.0 2.2796 0.2 1.010 3.0 4.3306 0.4 1.040 4.0 11.302 0.6 1.092 5.0 27.240 0.8 1.166 6.0 67.235 1.0 1.266 7.0 168.59 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOabaeqabaGaamiwaiaaywW7caaMf8UaaGjbVlaaysW7caWGjbWaaSbaaSqaaiaaicdaaeqaaOGaaiikaiaad6gacaGGPaGaaGzbVlaaywW7caaMf8UaaGzbVlaadIfacaaMf8UaaGzbVlaaywW7caaMe8UaamysamaaBaaaleaacaaIWaaabeaakiaacIcacaWGUbGaaiykaaqaaiaaicdacaGGUaGaaGimaiaaywW7caaMf8UaaGymaiaac6cacaaIWaGaaGimaiaaicdacaaMf8UaaGzbVlaaywW7caaMf8UaaGymaiaac6cacaaIYaGaaGzbVlaaywW7caaMf8UaaGymaiaac6cacaaIZaGaaGyoaiaaiodacaaI4aaabaGaaGimaiaac6cacaaIXaGaaGzbVlaaywW7caaIXaGaaiOlaiaaicdacaaIWaGaaGOmaiaaywW7caaMf8UaaGzbVlaaywW7caaIYaGaaiOlaiaaicdacaaMf8UaaGzbVlaaywW7caaIYaGaaiOlaiaaikdacaaI3aGaaGyoaiaaiAdaaeaacaaIWaGaaiOlaiaaikdacaaMf8UaaGzbVlaaigdacaGGUaGaaGimaiaaigdacaaIWaGaaGzbVlaaywW7caaMf8UaaGzbVlaaiodacaGGUaGaaGimaiaaywW7caaMf8UaaGzbVlaaisdacaGGUaGaaG4maiaaiodacaaIWaGaaGOnaaqaaiaaicdacaGGUaGaaGinaiaaywW7caaMf8UaaGymaiaac6cacaaIWaGaaGinaiaaicdacaaMf8UaaGzbVlaaywW7caaMf8UaaGinaiaac6cacaaIWaGaaGzbVlaaywW7caaMf8UaaGymaiaaigdacaGGUaGaaG4maiaaicdacaaIYaaabaGaaGimaiaac6cacaaI2aGaaGzbVlaaywW7caaIXaGaaiOlaiaaicdacaaI5aGaaGOmaiaaywW7caaMf8UaaGzbVlaaywW7caaI1aGaaiOlaiaaicdacaaMf8UaaGzbVlaaywW7caaIYaGaaG4naiaac6cacaaIYaGaaGinaiaaicdaaeaacaaIWaGaaiOlaiaaiIdacaaMf8UaaGzbVlaaigdacaGGUaGaaGymaiaaiAdacaaI2aGaaGzbVlaaywW7caaMf8UaaGzbVlaaiAdacaGGUaGaaGimaiaaywW7caaMf8UaaGzbVlaaiAdacaaI3aGaaiOlaiaaikdacaaIZaGaaGynaaqaaiaaigdacaGGUaGaaGimaiaaywW7caaMf8UaaGymaiaac6cacaaIYaGaaGOnaiaaiAdacaaMf8UaaGzbVlaaywW7caaMf8UaaG4naiaac6cacaaIWaGaaGzbVlaaywW7caaMf8UaaGymaiaaiAdacaaI4aGaaiOlaiaaiwdacaaI5aaaaaa@0ADE@ The upper limit of the summation does not need to be too large, say about 20 at the most.

(a) Kaiser window for β=0, 3 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiabek7aIjabg2da9iaaicdacaGGSaGaaGjbVlaaiodaaaa@3C3D@ and 6 and M = 20, (b) Magnitude response corresponding to the windows in (a), (c) magnitude response for β=6 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiabek7aIjabg2da9iaaiAdaaaa@3949@ and M = 10, 20,For β=0 size 12{β=0} {} both the numerator and denominator are 1 and the Kaiser window becomes the rectangular window, and for β=5.44 size 12{β=5 "." "44"} {} it is near to the Hamming window. The parameter β size 12{β} {}is so chosen that the magnitude response of the filter lies in the allowed region (see Fig.5.1). Fig.5.18 shows the Kaiser for various values of β size 12{β} {} and M size 12{M} {}. Notice that the Kaiser is also symmetric. It has minimum stopband attenuation from 50dB (for β=4.54 size 12{β=4 "." "54"} {}) to 90dB (for β=8.96 size 12{β=8 "." "96"} {}), so we can choose a Kaiser window to have maximum stopband atternuation equal to the Blackman window (74, Table 5.2) or higher. For ideal filters (Fig.5.1 and Fig.5.2) the Kaiser window has equal ripple δ size 12{δ} {} in both the passband and stopband. For real filters the ripples are slightly different . In design we take the common ripple as δ=min(δp,δs) size 12{δ="min" \( δ rSub { size 8{p} } ,δ rSub { size 8{s} } \) } {} of the desired values δp,δs size 12{δ rSub { size 8{p} } ,δ rSub { size 8{s} } } {}. As said earlier, actually δ size 12{δ} {} is usually specified in term of dBs: A=−20log10δ size 12{A= - "20""log" rSub { size 8{"10"} } δ} {}, for example δ=0.01 size 12{δ=0 "." "01"} {} then A=40dB size 12{A="40""dB"} {}, δ=0.001 size 12{δ=0 "." "001"} {} then A=60dB size 12{A="60""dB"} {}. Since δ size 12{δ} {} is a fraction of 1 then A is always ponsitive. Kaiser has developed empirical formulae and guide for the application of the window. Following is design steps for lowpass filters (similarly for other frequency selective filters). 1/ Specifications of edge frequencies ωp,ωs size 12{ω rSub { size 8{p} } ,ω rSub { size 8{s} } } {}; ripples δp,δs size 12{δ rSub { size 8{p} } ,δ rSub { size 8{s} } } {} 2/Evaluate the minimum ripple A=−20log10[min(δp,δs)] size 12{A= - "20""log" rSub { size 8{"10"} } \[ "min" \( δ rSub { size 8{p} } ,δ rSub { size 8{s} } \) \] } {}, the cutoff frequency ω c =( ω p + ω s )/2 MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiabeM8a3naaBaaaleaacaWGJbaabeaakiabg2da9iaacIcacqaHjpWDdaWgaaWcbaGaamiCaaqabaGccqGHRaWkcqaHjpWDdaWgaaWcbaGaam4CaaqabaGccaGGPaGaai4laiaaikdaaaa@4371@ , the transition width △ω= ω s − ω p MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiabloBjwjabeM8a3jabg2da9iabeM8a3naaBaaaleaacaWGZbaabeaakiabgkHiTiabeM8a3naaBaaaleaacaWGWbaabeaaaaa@414B@ 3/ Determine the filter filter order M M=A−7.952.285Δω,A>21M=5.79Δω,A≤21alignl { stack { size 12{M= { {A - 7 "." "95"} over {2 "." "285"Δω} } matrix { } , matrix { {} # {} # {} } A>"21"} {} # size 12{M= { {5 "." "79"} over {Δω} } matrix { {} } , matrix { matrix { {} # {} } {} # {} # {} } matrix { {} } A <= "21"} {} } } {} 4/ Determine the parameter β size 12{β} {}, depending on A: β=0, A≤21 β=0.5842 (A−21) 0.4 +0.07886(A−21), 21≤A≤50 β=0.1102(A−8.7), A≥50 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOabaeqabaGaeqOSdiMaeyypa0JaaGimaiaacYcacaaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaadgeacqGHKjYOcaaIYaGaaGymaaqaaiabek7aIjabg2da9iaaicdacaGGUaGaaGynaiaaiIdacaaI0aGaaGOmaiaacIcacaWGbbGaeyOeI0IaaGOmaiaaigdacaGGPaWaaWbaaSqabeaacaaIWaGaaiOlaiaaisdaaaGccqGHRaWkcaaIWaGaaiOlaiaaicdacaaI3aGaaGioaiaaiIdacaaI2aGaaiikaiaadgeacqGHsislcaaIYaGaaGymaiaacMcacaGGSaGaaGzbVlaaywW7caaMf8UaaGOmaiaaigdacqGHKjYOcaWGbbGaeyizImQaaGynaiaaicdaaeaacqaHYoGycqGH9aqpcaaIWaGaaiOlaiaaigdacaaIXaGaaGimaiaaikdacaGGOaGaamyqaiabgkHiTiaaiIdacaGGUaGaaG4naiaacMcacaGGSaGaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caaMf8UaaGzbVlaaywW7caWGbbGaeyyzImRaaGynaiaaicdaaaaa@9F66@ 5/ Computer the window coefficients w(n) size 12{w \( n \) } {}from the formula (5.40) 6/ Compute the impulse response hd(n) size 12{h rSub { size 8{d} } \( n \) } {}of the desired filter which is linear phase 7/ Compute the designed impulse response h(n)=hd(n)w(n) size 12{h \( n \) =h rSub { size 8{d} } \( n \) w \( n \) } {}. 8/ Find the frequency response H(ω)=IDTFT[h(n)] MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaadIeacaGGOaGaeqyYdCNaaiykaiabg2da9iaadMeacaWGebGaamivaiaadAeacaWGubGaai4waiaadIgacaGGOaGaamOBaiaacMcacaGGDbaaaa@43E8@ to check with specifications. For bandpass filter the two transition bandwidths are usually equal, if not we take the smaller one.

Design examples Design a linear phase lowpass filter using the Hamming window to meet the specificationsCutoff frequency 2.5 kHz Transition width 1.65 kHz Stopband attenuation >50dB Sampling frequency 10kHz Solution The relation between the analog linear frequency F size 12{F} {} samples/sec (or Hz) and the digital angular frequency ω size 12{ω} {} radians/sample is (Equation 1.40) is ω = 2π F f s size 12{ω=2π { {F} over {f rSub { size 8{s} } } } } {} where fs size 12{f rSub { size 8{s} } } {} is the sampling frequency (samples/sec, or Hz). Thus Cutoff frequency ωc=2π2.510=0.5π size 12{ω rSub { size 8{c} } =2π { {2 "." 5} over {"10"} } =0 "." 5π} {} rad/sample Transition width Δω=2π1.6510=0.33π size 12{Δω=2π { {1 "." "65"} over {"10"} } =0 "." "33"π} {} rad/sample Table 5.2 shows that the Hamming satisfy (so can the Blackman) the maximum stopband attenuation ( −20log10δs) size 12{ - "20""log" rSub { size 8{"10"} } δ rSub { size 8{s} } \) } {} requirement. Table 5.2 includes the relation between transition width Δω size 12{Δω} {} and window order M size 12{M} {}: Δω = 6 . 6π M ⇒ M = 6 . 6π 0 . 33 π = 20 size 12{Δω= { {6 "." 6π} over {M} } matrix { {} # drarrow {} # {} } M= { {6 "." 6π} over {0 "." "33"π} } ="20"} {} Now we compute the causal impulse response hLP(n) size 12{h rSub { size 8{ ital "LP"} } \( n \) } {}of ideal lowpass filter (Equation (5.32)) for 0≤n≤20 size 12{0 <= n <= "20"} {}. Notice that for M size 12{M} {} even and h(n) size 12{h \( n \) } {}symmetric we have linear phase FIR-1 (Fig.5.4a). The result is plotted in Fig.5.19a. Next we evaluate the Hamming window w(n)=0.54−0.46cos⁡ 2πn M , 0≤n≤M=20 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipDI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaadEhacaGGOaGaamOBaiaacMcacqGH9aqpcaaIWaGaaiOlaiaaiwdacaaI0aGaeyOeI0IaaGimaiaac6cacaaI0aGaaGOnaiGacogacaGGVbGaai4CamaalaaabaGaaGOmaiabec8aWjaad6gaaeaacaWGnbaaauaabeqabeaaaeaaaaGaaiilauaabeqabiaaaeaaaeaaaaGaaGimaiabgsMiJkaad6gacqGHKjYOcaWGnbGaeyypa0JaaGOmaiaaicdaaaa@515E@ The result is plotted in Fig.5.19b. The impulse response of the designed filter is the product h(n)= h LP (n)w(n) , 0≤n≤M=20 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipDI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqadeaabiGaaiaacaqaaeaadaqaaqaaaOqaaiaadIgacaGGOaGaamOBaiaacMcacqGH9aqpcaWGObWaaSbaaSqaaiaadYeacaWGqbaabeaakiaacIcacaWGUbGaaiykaiaadEhacaGGOaGaamOBaiaacMcafaqabeqabaaabaaaaiaacYcafaqabeqacaaabaaabaaaaiaaicdacqGHKjYOcaWGUbGaeyizImQaamytaiabg2da9iaaikdacaaIWaaaaa@4BCB@ We can evaluale w(n) and h LP (n) MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaadIgadaWgaaWcbaGaamitaiaadcfaaeqaaOGaaiikaiaad6gacaGGPaaaaa@3AF8@ separately as above then take the product, or evaluate the product directly . The result is plotted in Fig.5.19c. We then take the inverse transform of h(n) size 12{h \( n \) } {} to get the frequency response H(ω) size 12{H \( ω \) } {}: H ( ω ) = ∑ n = − ∞ ∞ h ( n ) e − jωn = ∑ n = 0 M h ( n ) e − jωn , M = 20 size 12{H \( ω \) = Sum cSub { size 8{n= - infinity } } cSup { size 8{ infinity } } {h \( n \) } e rSup { size 8{ - jωn} } = Sum cSub { size 8{n=0} } cSup { size 8{M} } {h \( n \) e rSup { size 8{ - jωn} } } matrix { } , matrix { {} # {} } M="20"} {} Fig.5.19d shows |H(ω)| MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaacYhacaWGibGaaiikaiabeM8a3jaacMcacaGG8baaaa@3BD6@ , and Fig.5.19c shows |H(ω) | dB MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaacYhacaWGibGaaiikaiabeM8a3jaacMcacaGG8bWaaSbaaSqaaiaadsgacaWGcbaabeaaaaa@3DB2@ . We need ro extract a low frequency signal having bandwidth 200Hz and center frequency 300Hz from a blackground noise. Design a bandpass filter having transition width 100Hz, passband ripple 40dB, stopband ripple 60dB. The sampling frequency is 1200Hz. Solution The passband of the signal is from 300 – 200/2 to 300 + 200/2, i.e. from 200Hz to 400Hz. We are to design a bandpass FIR filter to meet the specifications Passband: 200Hz Transition width: 100Hz

(linear phase lowpass filter with ω c =0.5π, Δω=0.33π MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiabeM8a3naaBaaaleaacaWGJbaabeaakiabg2da9iaaicdacaGGUaGaaGynaiabec8aWjaacYcacaaMe8UaeuiLdqKaeqyYdCNaeyypa0JaaGimaiaac6cacaaIZaGaaG4maiabec8aWbaa@48D5@ ) using the Hamming window

(linear phase bandpass filter using Blackman window Passband ripple: (−20log10δp)=40dB(δp=0.01) size 12{ \( - "20""log" rSub { size 8{"10"} } δ rSub { size 8{p} } \) ="40""dB" matrix { } \( δ rSub { size 8{p} } =0 "." "01" \) } {} Stopband ripple: (−20 log⁡ 10 δ p )=60dB⇒ δ s =0.001 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaacIcacqGHsislcaaIYaGaaGimaiGacYgacaGGVbGaai4zamaaBaaaleaacaaIXaGaaGimaaqabaGccqaH0oazdaWgaaWcbaGaamiCaaqabaGccaGGPaGaeyypa0JaaGOnaiaaicdacaGGKbGaaiOqaiabgkDiElaaysW7cqaH0oazdaWgaaWcbaGaam4CaaqabaGccqGH9aqpcaaIWaGaaiOlaiaaicdacaaIWaGaaGymaaaa@5075@ Sampling frequency: 1200Hz From Table 5.2 see that only the Blackman window can satisfy the stopband ripple (so can the Kaiser window, see next subsection). Since the passband and stopband attenuations are different, we take the smaller one in the design, i.e. A size 12{A} {}= 60dB. The transition width and the filter order are, respectively, Δω = 2π 100 1200 = π 6 size 12{Δω=2π { {"100"} over {"1200"} } = { {π} over {6} } } {} M= 11.1π π/6 =66.6≈67 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaad2eacqGH9aqpdaWcaaqaaiaaigdacaaIXaGaaiOlaiaaigdacqaHapaCaeaacqaHapaCcaGGVaGaaGOnaaaacqGH9aqpcaaI2aGaaGOnaiaac6cacaaI2aGaeyisISRaaGOnaiaaiEdaaaa@46C4@ Thus the Blackman window is W(n)=0.42−0.5cos⁡ 2πn M +0.08cos⁡ 4πn M , 0≤n≤M=67 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqipDI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeqabiGaciaacaqaaeaadaqaaqaaaOqaaiaadEfacaGGOaGaamOBaiaacMcacqGH9aqpcaaIWaGaaiOlaiaaisdacaaIYaGaeyOeI0IaaGimaiaac6cacaaI1aGaci4yaiaac+gacaGGZbWaaSaaaeaacaaIYaGaeqiWdaNaamOBaaqaaiaad2eaaaGaey4kaSIaaGimaiaac6cacaaIWaGaaGioaiGacogacaGGVbGaai4CamaalaaabaGaaGinaiabec8aWjaad6gaaeaacaWGnbaaauaabeqabeaaaeaaaaGaaiilauaabeqabiaaaeaaaeaaaaGaaGimaiabgsMiJkaad6gacqGHKjYOcaWGnbGaeyypa0JaaGOnaiaaiEdaaaa@5B75@ Ideal bandpass filter impulse response (Equation 5.10b) h BP (n)= sin⁡ ω u n πn − sin⁡ ω l n πn MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiaadIgadaWgaaWcbaGaamOqaiaadcfaaeqaaOGaaiikaiaad6gacaGGPaGaeyypa0ZaaSaaaeaaciGGZbGaaiyAaiaac6gacqaHjpWDdaWgaaWcbaGaamyDaaqabaGccaWGUbaabaGaeqiWdaNaamOBaaaacqGHsisldaWcaaqaaiGacohacaGGPbGaaiOBaiabeM8a3naaBaaaleaacaWGSbaabeaakiaad6gaaeaacqaHapaCcaWGUbaaaaaa@4FE7@ where ω u = 2π 400 1200 = 2π 3 size 12{ω rSub { size 8{u} } =2π { {"400"} over {"1200"} } = { {2π} over {3} } } {} ω l =2π 200 1200 = π 3 MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeGabiGaciaacaqaaeaadaqaaqaaaOqaaiabeM8a3naaBaaaleaacaWGSbaabeaakiabg2da9iaaikdacqaHapaCdaWcaaqaaiaaikdacaaIWaGaaGimaaqaaiaaigdacaaIYaGaaGimaiaaicdaaaGaeyypa0ZaaSaaaeaacqaHapaCaeaacaaIZaaaaaaa@4512@ Since M is odd we have linear phase type 2 (FIR-2) . The designed impulse response is h ( n ) = h BP ( n ) w ( n ) size 12{h \( n \) =h rSub { size 8{ ital "BP"} } \( n \) w \( n \) } {} We then take the inverse transform of h(n) to get or the frequency response of the designed filter (Fig.5.20). Design a highpass FIR filter using Kaiser window to meet the specifications :Cutoff frequency 2,5 kHz Transition width 0.5 kHz Passband ripple 0.001 Stopband attenuation 40dB Sampling frequency 10kHz Solution First let’s derive the expression of the impule response of ideal highpass filter (Fig.5.1b) with generalized linear phase (Fig.5.4) (Fig.5.1b) with generalized linear phase (Fig.5.4). The frequency response of the causal highpass filter in the interval ω=[0,π] MathType@MTEF@5@5@+=feaagaart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOqaaiabeM8a3jabg2da9iaacUfacaaIWaGaaiilaiabec8aWjaac2faaaa@3D9D@ is H HP (ω)=0, 0≤ω≤ ω c e −j M 2 ω , ω c <ω≤π MathType@MTEF@5@5@+=feaagaart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVCI8FfYJH8YrFfeuY=Hhbbf9v8qqaqFr0xc9pk0xbba9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr0=vqpWqabeaabiGaciaacaqaaeaadaqaaqaaaOabaeqabaGaamisamaaBaaaleaacaWGibGaamiuaaqabaGccaGGOaGaeqyYdCNaaiykaiabg2da9iaaicdacaGGSaqbaeqabeGaaaqaaaqaauaabeqabmaaaeaaaeaaaeaaaaaaaiaaicdacqGHKjYOcqaHjpWDcqGHKjYOcqaHjpWDdaWgaaWcbaGaam4yaaqabaaakeaafaqabeqadaaabaaabaaabaqbaeqabeqaaaqaaaaaaaqbaeqabeqaaaqaaaaacaaMe8UaaGjbVlaadwgadaahaaWcbeqaaiabgkHiTiaadQgadaWcaaqaaiaad2eaaeaacaaIYaaaaiabeM8a3baakiaacYcacaaMf8UaaGzbVlaaysW7caaMe8UaeqyYdC3aaSbaaSqaaiaadogaaeqaaOGaeyipaWJaeqyYdCNaeyizImQaeqiWdahaaaa@60D7@ or related to the frequency response of the lowpass filter HLP(ω) size 12{H rSub { size 8{ ital "LP"} } \( ω \) } {}(Equation (5.32)) as H HP ( ω ) = e − j M 2 ω − H LP ( ω ) size 12{H rSub { size 8{ ital "HP"} } \( ω \) =e rSup { size 8{ - j { {M} o