random.m

来自「generate a normal sequence from rayleigh」· M 代码 · 共 36 行

%Generate a random sequence of such length definded in the arguement..
%Programmed by Benjamin.Frank, 05-Mar-2009
%

function y = random(n)
% y = ceil(10*rand(1,n))
% n = 5;
% n = input('Enter the number of the random sequence:');
% n = 10;
if 1,
    y(1) = ceil(n*rand);
    for i=2:n,
        y(i) = ceil(n*rand);
%         for j=1:i-1,
%             if y(j)~=y(i)
%                 ;
%             else 
%                 y(i) = ceil(n*rand);
%                 break;
%             end
%              j = 1;
             while length(find(y([1:i-1])-y(i)))<length(y([1:i-1])),
                 y(i) = ceil(n*rand);
%                  if j<i-1,
%                      j = j+1;
%                  else
%                      break;
%                  end
             end
%         end
    end
end