【分享】matlab各种产生随机数函数的使用方法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1.rand
rand产生的是0到1(不包括1)的随机数.
matlab的rand函数生的是伪随机数,即由种子递推出来的,相同的种子,生成相同的随机数.
matlab刚运行起来时,种子都为初始值,因此每次第一次执行rand得到的随机数都是相同的.
1.多次运行,生成相同的随机数方法:

用rand('state',S)设定种子
S为35阶向量，最简单的设为0就好
例:
rand('state',0);rand(10)
2. 任何生成相同的随机数方法:
试着产生和时间相关的随机数,种子与当前时间有关.
rand('state',sum(100*clock))
即:
rand('state',sum(100*clock)) ;rand(10)
只要执行rand('state',sum(100*clock)) ;的当前计算机时间不现,生成的随机值就不现.
也就是如果时间相同,生成的随机数还是会相同.
在你计算机速度足够快的情况下,试运行一下:
rand('state',sum(100*clock));A=rand(5,5);rand('state',sum(100*clock));B=rand(5,5);
A和B是相同.
所以建议再增加一个随机变量,变成:
rand('state',sum(100*clock)*rand(1));
%
据说matlab 的rand 函数还存在其它的根本性的问题,似乎是非随机性问题

clear
b=round(7*rand(1,10000))+21;
a21=length(find(b==21));
a22=length(find(b==22));
a23=length(find(b==23));
a24=length(find(b==24));
a25=length(find(b==25));
a26=length(find(b==26));
a27=length(find(b==27));
a28=length(find(b==28));
[a21;a22;a23;a24;a25;a26;a27;a28]
ans =
         721
        1442
        1340
        1404
        1393
        1447
        1495
         758
修改:先多产生两个,两头一头多一个数,然后去掉最大和最小的两个数,或多产生两头两个后,直接取中间的数,方法是
a=[3 1 4 1 5 9 2 6 2 7 1 8 2 8 1 8 2 8];
b=a;
b([find(b==1)])=[];
b([find(b==8)])=[];
a,b

或
clear
b=round(9*rand(1,10000))+20;
a21=length(find(b==21));
a22=length(find(b==22));
a23=length(find(b==23));
a24=length(find(b==24));
a25=length(find(b==25));
a26=length(find(b==26));
a27=length(find(b==27));
a28=length(find(b==28));
[a21;a22;a23;a24;a25;a26;a27;a28]
2.exprnd

3.unifrnd
产生21到28之间的均匀分布随机数:
clear
a=unifrnd(20.5+eps,28.5-eps,1,5000);
b=round(a);
a21=length(find(b==21));
a22=length(find(b==22));
a23=length(find(b==23));
a24=length(find(b==24));
a25=length(find(b==25));
a26=length(find(b==26));
a27=length(find(b==27));
a28=length(find(b==28));
[a21;a22;a23;a24;a25;a26;a27;a28]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB atlab matla Mat Lab 计算机 clock matlab