- % 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