www.pudn.com > aft.rar > plot_5_23.m

% 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