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

ary_scl(temp1,Jk_hat);
ary_eql(temp1,delp_k);
delx_k=sk1*delp_k[0] + sk2*delp_k[1];
xk_corkt=xk+delx_k;
xk=xk_corkt;
alpha1=alpha1*0.9995 + 0.0005;
return(xk);
}
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
close;clear
! f51_2main
load f512.dat;
mrem=f512;
p1 = mrem(:,1);
p2 = mrem(:,2);
p3 = mrem(:,3);
p4 = mrem(:,4);
p5 = mrem(:,5);
k_c_code = length(p1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%
for k=1:k_c_code;
Dz=[1 -p1(k) -p2(k)];
R=roots(Dz);
rho(1,k)=abs(R(1)); theta(1,k)=angle(R(1));
rho(2,k)=abs(R(2)); theta(2,k)=angle(R(2));
x(k)=p4(k); % Output Signal
y(k)=p5(k); % Input Signal
A(k)=p1(k);
B(k)=p2(k);
C(k)=p3(k);
T(k)=k;
%%%
end
s=800; e=850;