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

p=2*cos(theta);
a3=r^3*(2*p-p^3); a4=r^4*(p*p-1);
b2=r*p;b3=r*r*(p*p-1);
a1=-r*p; a2=r*r;
a11=[1 0 0 a3 a4];
a11=1;
a22=1;
b22=[1 b2 b3];
a33=[1 a1 a2];
b33=1;
[h1 w]=freqz(b33,a33,180); w=w/(2*pi);
[h2 w]=freqz(b11,a11,180); w=w/(2*pi);
m1=abs (h1); m2=abs(h2);
rr=roots(a11); R=abs(rr);Theta=angle(rr);
A(:,k)=R;B(:,k)=Theta;
%polar(Theta,R,'o');grid;
%pause
k=k+1;
end
polar(B,A,'o');grid;
%plot(w,h1,'.-',w,h2);grid
A.8.4 Program f 3_rate
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
close;clear;
N=3;
i=1:N;
a=[sum(i.*i) sum(i)
sum(i) N ];
b=inv(a);
t=1:N; u = ones(size(t));
x= [4 5 6];
c=[t*x' u*x'] ';
mc=b*c
a1=1;
Nmax = 20;
b1 = 2*(1:Nmax);
b2 = (Nmax+1)*(ones(size(b1)));
t1=1:200;
noise=randn(size(t1))*6;
x1=sin(2*pi*0.01*t1);