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

uk=cos(2*pi*f1(i)*k);
% yk=mag_true(i)*cos(2*pi*f1(i)*k+phse_true(i)*pi/180);
yk=mag(i)*cos(2*pi*f1(i)*k+phase(i)*pi/180);
temp(4-k)=yk;
temp(7-k)=uk;
Uk=cos(2*pi*f1(i)*K);
Yk=mag(i)*cos(2*pi*f1(i)*K+phase(i)*pi/180);
end
% subplot (211); stem (Uk);
% subplot(212);stem(Yk);
temp(2:3)=-temp(2:3);
Y(i)=temp(1); A(i,:)=temp(2:6);
% pause;
end
z=Y';
p(1:2)=a(2:3); p(3:5)=b;
phat=inv(A'*A)*A'*z;
A.4.4 Program f 4_6_2
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
%Generates figures 4.15,4.16
clear;clf;
w=20; T=100; f=0.3;records=6;
x=ones([1 T]);
x(w+1:T)=0*x(w+1:T);
for i=1:records
s=(i-1)*T+1; e=s+T-1;
z(s:e)=x;
end
t=1:length(z); noise=1.0*randn(size(t));
y=sin(2*pi*f*t)+noise;
puls=z.*y;
subplot(211);plot(t,puls);grid;
%title('A');
s1=1;e1=s1+w-1;
for i=1:records
s=(i-1)*T+1; e=s+w-1;
xnew(s1:e1)=puls(s:e);
s1=e1+1;e1=s1+w-1;
end
[sxnew,fx]=spk(xnew,1);