sobolgenerator.m

来自「very good matlab of Monte carlo」· M 代码 · 共 25 行

%part a
function x=sobolgenerator(n,m1)

%any big N
N=30;
x=zeros(n+1 ,1); %0;
%x= 0;
for i=1:n
   %[dimension i]
   %to find right hand 0 we can look for the first 
   c=1;
   dived=i-1;
   remain=mod(dived,2);
   while remain==1
      dived=floor(dived/2);
      remain=mod(dived,2);
      c=c+1;
   end;
   x(i+1,1)=bitxor(floor(2^N*x(i,1)),floor(2^(N-c)*m1(c)))/2^N;
   %save temporary x
end;