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

% make figure 7.7 
 
Z2=zeros(4,491); 
load run2 
tail=[0 5 8 66];  
z=sum(real(E).^2,2); 
[y ind]=sort(z); 
newE=E(ind,:); 
for i=1:4, 
   Z2(i,:)=sum((real(newE(1:400-tail(i),:)).^2))/(400-tail(i)); 
end 
figure(1);semilogy(Z2','k') 
title('figure 7.7') 
print -dps 7_7