斐波那契等随机序列

1、斐波那契拟序列

function fibonacci
%(1-10维)k=23-37均匀性好,
%(1-20维)k=33-37均匀性好,
%(1维-60维)k=71-73均匀性好,73最好
k=73;%斐波那契参数:33
n=1000;%拟随机数个数
s=60;%维数
%产生斐波那契拟随机数系列
for i=1:s-1     
      quad=fibonacci_qrn(i,n,k);         
      x=quad';
      quad=fibonacci_qrn(i+1,n,k);           
      y=quad';      
      for j=1:n
          xx(j)=x(j,i);
          yy(j)=y(j,i+1);
      end 
      disp(sprintf('%f ',xx)); 
      disp(sprintf('%f ',yy)); 
      %output_filename=strcat('faure_',num2str(m),'_',num2str(n),'.txt');      
      %str0=strcat('第',num2str(i),'维斐波那契拟随机数序列一维分布均匀性'); 
      str1=strcat('第',num2str(i),'维'); 
      str2=strcat('第',num2str(i+1),'维');      
      %figure
      %plot(1:n,xx,'.k','MarkerSize',10)      
      %xlabel('拟随机数个数');
      %ylabel('拟随机数数值');