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 |
|