分享

matlab中生成无重复随机整数的方法

 成长中辉煌 2015-05-16
  Matlab自带函数randperm(n)产生1到n的整数的无重复的随机排列,利用它就可以得到无重复的随机数。例如:


randperm(n):产生一个1到n的随机顺序。 
>> randperm(10)
ans =
     6     4     8     9     3     5     7    10     2     1


randperm(n)源代码:
function p = randperm(n);
%RANDPERM Random permutation.
% RANDPERM(n) is a random permutation of the integers from 1 to n.
% For example, RANDPERM(6) might be [2 4 5 6 1 3].
%
% Note that RANDPERM calls RAND and therefore changes RAND's state. %
% See also PERMUTE. % Copyright 1984-2002 The MathWorks, Inc.
% $Revision: 5.10 $ $Date: 2002/04/09 00:26:14 $
[ignore,p] = sort(rand(1,n));


randperm(n)生成无重复整数的原理:


1. rand(1, n)产生1行n列的0-1之内的随机数矩阵。
2. sort()把这个矩阵排序,返回的ignore是排序后的序列,p是排序后的序列的各数原来的索引(即位置的编号),这个索引肯定是随机的,而且是在1到n间无重复的整数。


question:randperm(n)只能产生一个长度为n的且每个数都在1到n之间的随机顺序,我想产生一个长度为m(m<=n)的且每个数都在1到n之间的随机顺序,又该如何呢?


很简单,可以令p=randperm(n);然后选p的前m个数即可:p(1,m);


例如:m=3,n=8时
>> p=randperm(8);
>> p(1:3)                                                      

ans = 

 2     1     6

 

当然也可以把上述两步编程为一个子函数randnorepeat(m,n)
 

function d=randnorepeat(m,n)
%生成一列在[1,n]范围内的m个不重复的整数
p=randperm(n);
d=p(1:m);
 


例如:

>> d=randnorepeat(5,10)
 

d =
     3     2     5     6    10


附:《Matlab中的几个随机函数》

rand

rand(n):生成0到1之间的n阶随机数方阵

rand(m,n):生成0到1之间的m×n的随机数矩阵

randint

randint(m,n,[1 N]):生成m×n的在1到N之间的随机整数矩阵,randint(m,n,N+1)生成m×n的在0到N之间的随机整数矩阵,注:高版本的matlab中将randint改为randi,使用方法为randint([1 N],m,n)

 

>> randint(3,4,[1 10])
 

ans =

     3    10     4    10
     3     3     2     9
     1     3     3     2

>> randint(3,4,11)
 

ans =

     4     4     3     2
     9     5     0     9
    10     5     6     6
 


>> randi([1 10],3,4) 

ans =                                                           


     8     1     3     7
    10     1     5     9
     7     1     8     3
 


randperm


randperm(n):产生一个1到n的随机顺序。 
>> randperm(5)


ans =
     4     3     2     5     1
   

    本站是提供个人知识管理的网络存储空间,所有内容均由用户发布,不代表本站观点。请注意甄别内容中的联系方式、诱导购买等信息,谨防诈骗。如发现有害或侵权内容,请点击一键举报。
    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多