% Ce programme permet de simuler le modèle alpha de turbulence (Schertzer
% et Lovejoy, 1983, 1985).
% J.-C. POGGIALE -- Mars 2007

% Dans cet exemple, lambda^c=1/2.

clear;
clc;
gammamoins=-1;
%gammaplus=0.5;
lambda=2;
gammaplus=log(2-lambda^gammamoins)/log(lambda);
eps=1;
n=10;
for k=1:n
    eps2=zeros(lambda^k,1);
    for j=1:lambda^(k-1);
        pr=(rand(1,1)<0.5);
        eps2(lambda*j-pr,1)=eps(j,1)*2^(-gammamoins);
        eps2(lambda*j-(1-pr),1)=eps(j,1)*2^(-gammaplus);
    end
    eps=eps2;
end
plot(eps)
