fid1=fopen('z1.bin','r');
[xxx,l1]=fread(fid1,'ubit64');
xxx=reshape(xxx,8,128);
fid2=fopen('x.bin','r');
[xx,l2]=fread(fid2,'ubit64');
xx=xx';
xx=xx-x2;
xx=xx*256;
xx=reshape(xx,8,128);
xxxx=xx+xxx;
xxxx=xxxx';
 M=2^20;
Na=10;
L=0;
time=M/Na;
time=fix(time);
a=round(rand(1,M));
d=[0:Na*8-1];
s=[0,a(1),a(2),a(3),a(4)];
x=5;
 tt=[0:Na*8-1];
hold on
  for i=1:500
      for j=0:9
        s=[s(2:5),a(x)];
        if x==5,
        L=0;
        else a(x-5)=2*a(x-5)-1;       
        L=L+a(x-5); 
        end,
        x=x+1;        
        L=rem(L,4);
        if L<0,L=L+4;end,
        su=s(1)*64+s(2)*32+s(3)*16+s(4)*8+s(5)*4+L+1;
        bb=xxxx(su,(1:8));
        d(1+j*8:8+j*8)=bb;
    end
    plot(tt,d);
end