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

clear;
i=1;
for phi=20:1:150
r=0.9; theta=phi*(pi/180); p=2*cos(theta);
a1=-r*p; a2=r*r;
P1=[1 a1 a2]; P2=[1 -a1 a1*a1-a2];
P3=conv(P1,P2);
a3=r^3*(2*p-p^3); a4=r^4*(p*p-1);
P4=[1 0 0 a3 a4];
P4-P3;
R=abs(roots(P2)); r1(i)=R(1);r2(i)=R(2);Theta(i)=phi;
%pause
i=i+1;
end
subplot(211)
plot(Theta,r1,'*',Theta,r2,'o');grid
%pause
%print -depsc f51_1
A.5.3 MATLAB Program f 51_2
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
%Generates figure 5.13
close;clear
kmax=5000; theta=20; f=theta/360;
t=1:kmax;
z=sin(2*pi*f*t)+0.2*randn(size(t));
a1=0.1; a2=0.1;
r=0.98;
v=0.0;
Hk=[1e-9 0 0
0 1e-9 0
0 0 1e-9];
delp_k_old=[0 0 0]';
delp_k= [0 0 0]';
xk=0; xk1=0.1; sk=0; sk1=0; sk2=0; xk_corkt=0;
zk=0; zk1=0;
del_a1=0; del_a2=0;
old_del_a1=0; old_del_a2=0;
alpha1=0.95; delx_k=0;
pk=[a1 a2 v]';