x=rand(1,1024);                    
x=wfbm(0.3,1024);
stem(x);                           
xlabel('随机序列')                  
n=size(x,2);                       
s=sum(x);                        
Tt=log2(n);                        
e=2.^(1:Tt);                           
N=n./e;                           
for k=1:Tt
    d=N(k);
    for m=1:e(k)
        A=x(d*(m-1)+[1:d]);
        s0=sum(A);
        pe(m,k)=s0/s;
    end
end         
能解释一下每步程序的意思吗，谢谢