8.2 Programs
141
eta=zeros(N,NR);
for j=1:NR
for i=1:N-1
Pat1at = 0.5*(1. + exp(-i*dt/tau)); %Probability to stay in "-1"
Pat2at1 = Pat1at;
Pat1bt = 0.5*(1. - exp(-i*dt/tau)); %Transition Probability 1->0
Pat2bt1 = Pat1bt;
Pbt2bt1 = 0.5*(1. + exp(-i*dt/tau)); %Probability to stay in "1"
Pbt1bt = Pbt2bt1;
Pbt1at = 0.5*(1. - exp(-i*dt/tau)); %Transition Probability 0->1
Pbt2at1 = Pbt1at;
R=unifrnd(0,1);
if(Pat1at > R)
eta(i+1,j)=-1.;
else
eta(i+1,j)=1.;
end
if(Pat1at < R || Pat1at > R)
R1=unifrnd(0,1);
elseif(Pat2bt1 > R1 || Pat2at1 > R1) % || significa or eta(i+2,j)=-1.;
else
eta(i+2,j)=1.;
end
end
end
%======================================================
tR=tc/ts;
INVTR=1/tR;
A1=6.50618; A2=3.89653; A3=3.01154;
t1=8.33507e-5; D1=A1*t1;
t2=0.01425; D2=A2*t2;
t3=0.3896; D3=A3*t3;
x=zeros(N-1,NR);
v=zeros(N-1,NR);
Dd1U=zeros(N-1,NR);
Dd2U=zeros(N-1,NR);
Dd3U=zeros(N-1,NR);
DDeff=zeros(N-1,NR);
eps=zeros(N-1,NR);
h=zeros(N-1,NR);
eps1=zeros(N-1,NR);
h1=zeros(N-1,NR);
eps2=zeros(N-1,NR);