MATLAB and C Programs - Digital Signal Processing: A Practitioner's Approach - page 148

phi_new=angle(freqz(bnew,1,N))*180/pi;
subplot(221);
stem(bnew);grid;title('A');ylabel('hk');
subplot(222);plot(fn,phi,fn,phi_new);grid;title('B');ylabel('deg');
pause;
%print -depsc f3_7;
m_new=abs(freqz(bnew,1,N));m_new=m_new/(max(m_new));
subplot(211);plot(fn,m,fn,m_new);grid;
pause;
%print -depsc f3_8;pause;
A.3.11 Program f 3_789
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
%Generates figure 3.21,3.22
clear
close
[m1,f]=freq(0.5);
[m2,f]=freq(0.6);
[m3,f]=freq(0.7);
[m4,f]=freq(0.8);
[m5,f]=freq(0.9);
[m6,f]=freq(0.95);
[m7,f]=freq(0.98);
mn=[m1,m2,m3,m4,m5,m6,m7];
[m11,f]=freqt(36);
[m21,f]=freqt(54);
[m31,f]=freqt(72);
[m41,f]=freqt(90);
[m51,f]=freqt(108);
[m61,f]=freqt(126);
[m71,f]=freqt(144);
m12=[m11 m21 m31 m41 m51 m61 m71];
subplot(211)
plot(f,m12);
xlabel ('Normalised Frequency');
ylabel ('Magnitude')
grid;
%gtext('theta = 36');
%gtext('theta = 144');
pause
%print -depsc f31_8