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

% make figure 7.8 
Z3=zeros(4,491); 
t=66; 
for i=1:4, 
   eval(['load run' num2str(i)]); 
   z=sum(real(E).^2,2); 
   [y ind]=sort(z); 
   newE=E(ind,:); 
   Z3(i,:)=sum((real(newE(1:400-t,:)).^2))/(400-t);  
z1(i,:)=sum(real(E).^2,1); 
end 
figure(3); semilogy(Z3') 
hold on
title('figure 7.8') 
print -dpsc 7_8