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

% remake figure 5.17 
 
figure 
load run5;plot(mean(Wx'), '--'); hold on 
 
mu=rp.mu; n=1:5000; 
a=-0.99; lambda=abs(a); Ew1=-a+a*(1-mu*lambda).^n; 
a=0.99;  lambda=abs(a); Ew2=-a+a*(1-mu*lambda).^n; 
plot(Ew1); plot(Ew2) 
load run5b;plot(mean(Wx'),'--'); 
line([1 5000], [0 0]);hold off 
text(250,0.62,'a = - 0.99'); text(250,-0.62,'a = + 0.99') 
title(' figure 5.17 \mu=0.001') 
     legend('Experiment', 'Theory') 
print -dpsc 5_17