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

plot(F,amp);
grid;
xlabel p
ylabel J(p)
%title (' Objective function for r=0.95 at 3dB SNR ');
pause
%print -depsc f31_9
subplot(211)
plot(F,slope);
grid;
xlabel p
ylabel ( 'slope of [J(p)]')
pause
%print -depsc f32_0
%zoom
A.3.13 Program f 31_9a
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
%Generates figure 3.29
clear;clf;
b1=[1 0.6]; a1=[1 -1.7 1.53 -0.648];
[h,w]=freqz(b1,a1,180) ;m1=abs(h);f=w/(2*pi);
subplot(221);
plot(f,m1);grid;title('A');xlabel('Frequency');
b2=0.999317; a2=[1 -2.16988 2.41286 -1.47436 0.365412];
[h,w]=freqz(b2,a2,180);m2=abs(h);f=w/(2*pi);
subplot(222);
plot(f,m2,f,m1,'.-');grid;title('B');xlabel('Frequency')
%pause;
%print -depsc f31_9a
A.3.14 Program f 3_6_1a
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
%Generates figure 3.31
clear;close;
vw=15;N=2000;
yw_max=N;yw=0:yw_max;
delta_y=2*pi/yw_max;
xw=500*sin(delta_y*yw);