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

pause
print -depsc f6_9
A.3.17 Program f 3_61_01
%Digital Signal Processing:A Practitioner's Approach
%Dr.Kaluri Venkata Ranga Rao, [e-mail]
%Generates figures 3.26,3.25
%This program needs rbn_sig.m and expand.m files.
clear;
close;
%N=5000;
% generate BFSK data
xx=sign(randn([1 200]));
fsk=rbn_sig(150,xx);
ii=find(fsk==0);
fsk=fsk*0.27;fsk(ii)=fsk(ii)+0.23;N=length(fsk);
kk=1:N; z=sin(2*pi*fsk.*kk);
[sz,ff]=spk(z,1);qq=800;
subplot(221);
plot(kk(1:qq), z(1: qq)*0.2+1.1,kk(1:qq),fsk(1:qq)*3);grid;
title('A');xlabel(' Time subplot(222);plot(ff,sz);grid;
title('B');xlabel(' Frequency ');
pause;
% load fsk_data
% K=1:N;
% z=sin(2*pi*f*K);
noise=0.2236*randn(size(z));
z=z+0.1*noise;
lambda=0.9;
mu =0.002;
r=0.9;
vk=0;
delp =0;
N=5000;
range_p = 4;
delP = range_p/N;
xk1=0; xk=0;
sk1=0; sk=0;
uk1=0; uk=0;
P=0;
for k=1:N
uk2=uk1; uk1=uk;