电子科技大学《数字信号处理》Ch10FIR digital Filter Design.pptVIP

% Lowpass Filter Design Using the Kaiser Window fpts = input(Type in the bandedges = ); mag = input(Type in the desired magnitude values = ); dev = input(Type in the ripples in each band = ); [N, Wn, beta, ftype] = kaiserord(fpts,mag,dev) kw = kaiser(N+1,beta); b = fir1(N,Wn, kw); [h,omega] = freqz(b,1,512); plot(omega/pi,20*log10(abs(h)));grid; xlabel(\omega/\pi); ylabel(Gain, dB); Example 10.25 Type in the bandedges = [0.3 0.4] Type in the desired magnitude values = [1 0] Type in the ripples in each band = [0.003162 0.003162] N = 59 Wn = 0.3500 beta = 4.5513 ftype = low Example 10.26 Use the same program as Example 10.25 by modifying b = fir1(N,Wn,high,kw); Type in the bandedges = [0.4 0.55] Type in the desired magnitude values = [0 1] Type in the ripples in each band = [0.02 0.02] N = 26 Wn = 0.4750 beta = 2.6523 ftype = high FIR Digital Filter Design Using MATLAB FIR Digital Filter Design Using MATLAB Design a linear phase FIR LPF: If given the passband edge frequency fp=300Hz, the passband ripple is 1dB，the stopband edge frequency f s=450Hz, the minimum stopband attenuation-25dB. The sampling rate is 3000Hz. Homework: 10.17, 10.20, 10.31, 10.32 M10.8, M10.9, M10.17 Passband and stopband ripples are same. Thus, Observe Distance between the locations of the maximum passband deviation and minimum stopband value ? ?ML Width of transition band: ?w = ?s - ?p ?ML To ensure a fast transition from passband to stopband, window should have a very small main lobe width. To reduce the passband and stopband ripple d, the area under the sidelobes should be very small. Unfortunately, these two requirements are contradictory. Rectangular and Blackman windows In the case of rectangular, Hann, Hamming, and Blackman windows, the value of ripple does not depend on filter length or cutoff frequency ?c, and is essentially constant. In addition, ?w ? c / M where c is a constant for most practical purposes. ?So, we can determine ?w b