24小时热门版块排行榜    

查看: 3167  |  回复: 9

田山东

捐助贵宾 (著名写手)

[求助] matlab如何生成不同的随机数

程序1.下面这个程序生成的随机数都是相同的
h=zeros(10,10);
t=1;
fid1=fopen('tyh4.dat','a');
for jj=1:4
rand('state',sum(1000000*clock));
x=rand(1,10)-0.5;
%fprintf(fid2,'%f\n',x);
for  j=1:10;
    H(j,j)=x(j)*t;
    fprintf(fid1,'%f\n', H(j,j));
end  
end
把文件打开看下,4次生成的随机数相同。
去掉rand('state',sum(1000000*clock));
h=zeros(10,10);
t=1;
fid1=fopen('tyh4.dat','a');
for jj=1:4
x=rand(1,10)-0.5;
%fprintf(fid2,'%f\n',x);
for  j=1:10;
    H(j,j)=x(j)*t;
    fprintf(fid1,'%f\n', H(j,j));
end  
end
4次循环生成的随机数不同,但是多次运行程序会发现,每次得到的文件内容都是相同的。
目的:每次jj循环得到不同的随机数,每次运行程序得到不同的随机数,怎么能实现呢?
回复此楼

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

everythinghasitsseason.enjoyyourlife.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

送鲜花一朵
这样就可以,大侠解释下吧
h=zeros(10,10);
t=1;
fid1=fopen('tyh66.dat','a');
for jj=1:4
rng('shuffle')
x=rand(1,10)-0.5;
%fprintf(fid2,'%f\n',x);
for  j=1:10;
    H(j,j)=x(j)*t;
    fprintf(fid1,'%f\n', H(j,j));
end  
end
everythinghasitsseason.enjoyyourlife.
2楼2012-12-12 10:14:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

在2011下可以了,我linux下用的是r2010b,为什么有这个提示呢?Undefined function or method 'rng' for input arguments of type 'char'
everythinghasitsseason.enjoyyourlife.
3楼2012-12-12 10:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

最后没有办法只能让系统产生一个大的随机数数组,没循环一次调用这个大数组里面的一行或一列,一直调用下去,反正大数组是大的随机数,不会重复。
everythinghasitsseason.enjoyyourlife.
4楼2012-12-12 11:04:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lurencyj

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你要进行撒种子seed,不然,都是一样的!~~~

» 本帖已获得的红花(最新10朵)

很女子很弓虽大
5楼2012-12-12 12:44:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

送鲜花一朵
引用回帖:
5楼: Originally posted by lurencyj at 2012-12-12 12:44:01
你要进行撒种子seed,不然,都是一样的!~~~

取随机数rand()。取系统时间作为随机数种子,是避免产生伪随机数,不然每次重新打开matlab所取的随机数序列是相同的。

%取0-1之间的随机数
    ctime = datestr(now, 30);%取系统时间
    tseed = str2num(ctime((end - 5) : end)) ;%将时间字符转换为数字
    rand('seed', tseed) ;%设置种子,若不设置种子则可取到伪随机数
%取0-1之间的随机数
    for i = 1:5      
        weight(i) = rand();
    end
    weight
这样可以吧?呵呵,网上搜到了
everythinghasitsseason.enjoyyourlife.
6楼2012-12-12 14:37:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

这样还是不行啊。妈的,。抑郁了!!
这样的话有下面的情况:
第一次运行程序,生成一组随机数1
第二次运行程序,又生成一组随机数2
1和2不同,。
但是!!!!!!程序循环10次,每次循环得到的随机数相同!!我哭了!!!
:cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry::cry:
everythinghasitsseason.enjoyyourlife.
7楼2012-12-12 19:04:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

rand('seed',jj+7);
好像每次改变jj+7的数值就可以,我试下。
everythinghasitsseason.enjoyyourlife.
8楼2012-12-12 19:54:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与,应助指数 +1
整那么复杂干嘛?我用下面的代码run了几次,每次的数也不一样,同一次10次循环得到的也不一样.
CODE:
for i=1:10
    fprintf(1,'%f\n',rand(1));
end

结果:
CODE:
>> temp
0.757740
0.743132
0.392227
0.655478
0.171187
0.706046
0.031833
0.276923
0.046171
0.097132
>> temp
0.823458
0.694829
0.317099
0.950222
0.034446
0.438744
0.381558
0.765517
0.795200
0.186873
>> temp
0.489764
0.445586
0.646313
0.709365
0.754687
0.276025
0.679703
0.655098
0.162612
0.118998
>>

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
9楼2012-12-12 20:21:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

田山东

捐助贵宾 (著名写手)

引用回帖:
9楼: Originally posted by libralibra at 2012-12-12 20:21:48
整那么复杂干嘛?我用下面的代码run了几次,每次的数也不一样,同一次10次循环得到的也不一样.
for i=1:10
    fprintf(1,'%f\n',rand(1));
end

结果:
>> temp
0.757740
0.743132
0.392227
0.65547 ...

超级郁闷,我用的是2011b,linux下装的是2010b,
for i=1:10
    fprintf(1,'%f\n',rand(1));
end
在2011b下完全是正确的,但是在2010下就悲剧了!!

找个最新版的matlab装下吧。没办法了
everythinghasitsseason.enjoyyourlife.
10楼2012-12-13 16:47:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 田山东 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大佬帮着分析今年的评审意见?该课题连续申请了三年了,明年要不要换个方向。 100+3 cai-11138 2026-08-31 11/550 2026-09-06 20:49 by archvillain
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:01 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 17:49 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才!!! +5 SCSIOyxguo 2026-09-03 8/400 2026-09-06 15:06 by 可爱小花朵
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +10 瞬息宇宙 2026-09-01 19/950 2026-09-06 11:15 by jiangxuan2004
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗? +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
信息提示
请填处理意见