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

xk=xk1;
step=step/1.5;
end
if change_in_g == -2
xk=xk1;
step=step/1.5;
end
g(i)=change_in_g;
%if g(i) == 2
% xl(i)=xk; xl_prv(i)=xk1;
%end
%if g(i) == -2
% xr(i)=xk; xr_prv(i)=xk1;
%end
end
subplot(211);stem(t,g);grid;
subplot(212);plot(t,x1,t,ref);grid;
%pause;
%subplot(211);plot(t,xl_prv, 'o', t,xl, '*', t,ref);grid
%subplot(212);plot(t,xr_prv, 'o', t,xr, '*', t,ref);grid
A.8.2 Program f 2_lemma
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
clear;close
A=rand([3 3])-0.5; A=fix(A*20);
B=A+A';
c=rand([3 1])-0.5; c=fix(c*20);
delB=c*c';
Z=inv(B+delB);invB=inv(B);
Zhat=invB - invB*c*c'*invB/(1+c'*invB*c);
Z-Zhat
A.8.3 Program f 5_pipe
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
clear;clf ;k=1;
for i=30:1:50
theta=i*(pi/180);r=0.6;