shen.m

神经网络计算双色球程序

function y=shen(Z)

m1=3;
m2=23;
w1=18;
w2=18;
k1=15;
k2=15;
aa=6;




for i=1:640
    for j=1:24
        H(j,i)=Z(i+j-1);
    end
end
save H H

Z=H;
n=length(Z);


for kk=k1:k2

for w=w1:w2;
    
    
for m=m1:m2
    
for i=1:n
   for j=1:24
       X(j,i)=(Z(j,i)-1)/w;
   end
end
    
for i=m:610
    for j=1:6
         for k=1:kk
             XX((j+(k-1)*6),i-m+1)=X(j,(k-1)+i);
         end
      
     end
end

for i=m:610
    for j=1:aa

       YY(j,i-m+1)=X(j,i-1);
     end
end



save XX XX
save YY YY




  for i=2:(m-1)
      for j=1:6
           for k=1:kk
               p((j+(k-1)*6),i-1)=X(j,(k-1)+i);
           end
        
      end
  end
    for i=2:(m-1)
      for j=1:aa

         Fi(j,i-1)=X(j,i-1);
      end
  end
  
    for i=1:1
      for j=1:6
           for k=1:kk
               J((j+(k-1)*6),i)=X(j,(k-1)+i);
           end
      end
    end
  save J J
  save p p 
  save Fi Fi



   XXX(1:(6*kk),1:120)=XX(1:(6*kk),1:120);
   YYY(1:aa,1:120)=YY(1:aa,1:120);



    

    net=newrbe(XXX,YYY,1);
      
     XXX=[];
     YYY=[];   
     XX=[];
     YY=[];         

    save net net 
    S(:,1)=sim(net,p(:,(m-2)));
    M(1,m-m1+1)=Fi(1,m-2);
    M(2,m-m1+1)=S(1,1);

            
        
        

end


save M M

for i=1:(m2-m1+1)
    F(i)=abs(sum(M(1,i)-M(2,i)));
   
    if F(i)<1/w
          Q(i)=0;
       
           else F(i)>1/w;
           Q(i)=1;
       end
    
end

R(kk-k1+1,w-w1+1)=sum(Q);


end

end
save R  R