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

% remake figure 9.14 
figure(1); title('figure 9.14'); hold on 
load run1;plot(mean(Wx'));hold on;plot(Wx(:,69)) 
load run2;plot(mean(Wx'));plot(Wx(:,69)) 
line([1 500], [0 0]) 
ylabel('Tap weight') 
xlabel('Number of iterations') 
text(300,0.6,'a = - 0.99') 
text(300,-0.5,'a = + 0.99') 
set(gca,'ylim',[-1.5 1.5]);hold off 
print -dpsc 9_14 
 
 
% remake figure 9.15 
figure(2); title('figure 9.15');hold on 
load run2 
semilogy(mean((E').^2)); 
plot(E(:,69).^2) 
set(gca,'ylim',[1e-3 1e0]) 
ylabel('Squared error'); xlabel('Number of iterations');hold off 
print -dpsc 9_15 
 
 
% remake figure 9.16 
figure(3); title('figure 9.16'); hold on 
load run3;semilogy(mean((E').^2),'r-.'); 
load run2;semilogy(mean((E').^2),'b-'); 
load run4;semilogy(mean((E').^2),'g--'); 
ylabel('Mean-squared error') 
xlabel('Number of iterations') 
legend('\mu=0.01','\mu=0.05','\mu=0.1') 
hold off 
print -dpsc 9_16 
 
 
% remake figure 9.17 
load run1 
mu=rp.mu; n=1:500; 
a=-0.99; lambda=abs(a); Ew1=-a+a*(1-mu*lambda).^n; 
a=0.99;  lambda=abs(a); Ew2=-a+a*(1-mu*lambda).^n; 
figure(4); title('figure 9.17'); 
plot(Ew1); hold on; plot(Ew2) 
load run1;plot(mean(Wx')) 
load run2;plot(mean(Wx')) 
line([1 500], [0 0]); 
text(250,0.6,'a = - 0.99'); text(250,-0.5,'a = + 0.99') 
hold off 
print -dpsc 9_17 
 
 
% remake figure 9.18 
figure(5); title('figure 9.18');hold on 
load run4 
mu = rp.mu; a  = rp.a; sigu2 = 0.995;  sigv2 = (1-(a^2)) * sigu2; n = 1:500;  
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);semilogy(mean((E').^2)); 
hold off 
print -dpsc 9_18