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