- figure(1)
- for i=1:4,
- eval(['load run' num2str(i)])
- P(1:488,i)=mean(A.^2)';
- end
- semilogy(P)
- set(gca,'xlim',[1 200])
- title('figure 13.6')
- print -dpsc 13_6
- figure(2)
- load run5
- semilogy(mean(A.^2)');
- eval(['load ../9_7_lms/run15'])
- Z=sum((E.^2),2)/200;
- hold on
- semilogy(Z)
- hold off
- set(gca,'xlim',[0 500])
- set(gca,'ylim',[1e-1 1e1])
- title('figure 13.7')
- print -dpsc 13_7