% Chapter 07: Example 7.8 Lowpass filter design - Hamming window%wp = 0.2*pi; ws = 0.3*pi;tr_width = ws - wpM = ceil(6.6*pi/tr_width) + 1n=[0:1:M-1];wc = (ws+wp)/2hd = ideal_lp(wc,M);w_ham = (hamming(M))';h = hd .* w_ham;[db,mag,pha,grd,w] = freqz_m(h,[1]);delta_w = 2*pi/1000;Rp = -(min(db(1:1:wp/delta_w+1)))        % Passband RippleAs = -round(max(db(ws/delta_w+1:1:501))) % Min Stopband attenuation% plotssubplot(1,1,1)subplot(2,2,1); stem(n,hd); title('Ideal Impulse Response')axis([0 M-1 -0.1 0.3]); xlabel('n'); ylabel('hd(n)')subplot(2,2,2); stem(n,w_ham);title('Hamming Window')axis([0 M-1 0 1.1]); xlabel('n'); ylabel('w(n)')subplot(2,2,3); stem(n,h);title('Actual Impulse Response')axis([0 M-1 -0.1 0.3]); xlabel('n'); ylabel('h(n)')subplot(2,2,4); plot(w/pi,db);title('Magnitude Response in dB');gridaxis([0 1 -100 10]); xlabel('frequency in pi units'); ylabel('Decibels')set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])set(gca,'YTickMode','manual','YTick',[-50,0])set(gca,'YTickLabelMode','manual','YTickLabels',['50';' 0'])