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

%remake plot 5.19 
 
figure 
load run5 
mu = rp.mu; a  = rp.a; sigu2 = 0.995;  sigv2 = (1-(a^2)) * sigu2; n = 1:5000;  
a2=a^2; p5mu=0.5*mu; 
j = sigu2*(1-a2)*(1+p5mu*sigu2) + sigu2*(a2+p5mu*a2*sigu2-0.5*mu*sigu2)*(1-mu*sigu2).^(2*n); 
semilogy(j,'--'); hold on 
semilogy(mean((E').^2));semilogy(j,'-'); 
title('figure 5.18 \mu=0.001') 
print -dps 5_18 
 
 
figure 
load run6 
mu = rp.mu; a  = rp.a; sigu2 = 0.995;  sigv2 = (1-(a^2)) * sigu2; n = 1:5000;  
a2=a^2; p5mu=0.5*mu; 
j = sigu2*(1-a2)*(1+p5mu*sigu2) + sigu2*(a2+p5mu*a2*sigu2-0.5*mu*sigu2)*(1-mu*sigu2).^(2*n); 
semilogy(j,'r'); hold on 
semilogy(mean((E').^2));semilogy(j,'r'); 
title('figure 5.19 \mu=0.005') 
%print -dpsc 5_19b