加入收藏
大学数学, 研究生数学,大学数学资料下载,免费 大学数学课件,研究生数学课件,免费下载
matlab中生成无重复随机整数的方法
2014-09-05 15:38:56

数学建模中常常需要产生随机数,一般的用rand()就可以解决,但是如果要生成无重复的整数随机数列,可以用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);  

[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

赞一个(151) | 阅读(11347)
上一篇:2009高教社杯全国大学生数学建模竞赛AB题评阅要点
下一篇:Matlab如何给满足条件的元素批量赋值
 

胡桃木屋版权所有@2013 湘ICP备13006789号-1