- % 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