www.pudn.com > adaptivefiltering.rar > plot_5_18.m

%remake plot 5.18 
 
figure 
load run5 
mu = rp.mu; a  = rp.a; sigu2 = 0.995;  sigv2 = (1-(a^2)) * sigu2; n = 1:5000;  
a2=a^2; p5mu=0.5*mu; 
j = sigu2*(1-a2)*(1+p5mu*sigu2) + sigu2*(a2+p5mu*a2*sigu2-0.5*mu*sigu2)*(1-mu*sigu2).^(2*n); 
semilogy(j,'--'); hold on 
semilogy(mean((E').^2));semilogy(j,'-'); 
%legend('Theory', 'Experiment') 
title('figure 5.18 \mu=0.001') 
print -dps 5_18