- % remake 5.22
- figure
- for i=4:-1:1,
- eval(['load run', num2str(i)])
- Z(:,i)=sum((E.^2),2)/200;
- j=abs(4-i)+1;
- subplot(4,1,j)
- stem(mean(Wo,2),'k')
- set(gca,'box','off')
- set(gca,'ylim',[-1 1.7])
- % set(gca,'xticklabel',[5 10])
- line([0 11],[0 0])
- end
- clear Z;
- subplot(4,1,1); title('figure 5.22')
- print -dps 5_22