www.pudn.com > adaptivefiltering.rar > plot_7_9_abcd.m

% make figure 7.9abcd 
base=3; 
strb=['a'; 'b'; 'c'; 'd']; 
for i=1:4, 
   figure; 
   eval(['load run' num2str(i)]) 
   z=sum(real(E).^2,2); 
   [y ind]=sort(z); 
   newE=E(ind,:); 
   t=66; 
   z_dct=sum((real(newE(1:400-t,:)).^2))/(400-t); 
   semilogy(z_dct); 
   hold on 
   eval(['load ../Chap_5_EQ/run' num2str(i)]) 
   semilogy(mean(E'.^2),'r') 
   eval(['load ../Chap_9/run' num2str(i)]) 
   semilogy(mean(A.^2),'m') 
   hold off 
   str=['title(''figure 7.9' strb(i) ''')']; 
   eval(str) 
   eval(['print -dpsc 7_9' strb(i) ]); 
end