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

figure(1) 
 
for i=1:4, 
	eval(['load run' num2str(i)]) 
	P(1:488,i)=mean(A.^2)'; 
end 
 
semilogy(P) 
set(gca,'xlim',[1 200]) 
title('figure 13.6') 
print -dpsc 13_6 
 
figure(2) 
 
load run5 
semilogy(mean(A.^2)'); 
 
eval(['load ../9_7_lms/run15']) 
Z=sum((E.^2),2)/200; 
hold on 
semilogy(Z) 
hold off 
 
set(gca,'xlim',[0 500]) 
set(gca,'ylim',[1e-1 1e1]) 
title('figure 13.7') 
print -dpsc 13_7