plot_7_8.m

matlab里面的一些经典的例子程序

% 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