Fortran Coder

查看: 10917|回复: 4
打印 上一主题 下一主题

[数值问题] 请教,随机数的问题

[复制链接]

8

帖子

4

主题

0

精华

入门

F 币
41 元
贡献
25 点
跳转到指定楼层
楼主
发表于 2014-3-24 12:09:42 | 显示全部楼层 回帖奖励 |倒序浏览 |阅读模式
不用写出程序,只求指出我的程序哪里出错了,谢谢啊
1.我要建一个模型,100元素的一维数组设为S,每个元素非1即-1.随机数判断产生。计算相邻的乘积值的负值首尾也算相邻-S(1)*S(100)设为M

2.再从1到一百取随机一个数J提出S(J),取S(J)=-S(J),计算M2,再随机取一个数X判断,如果X小于EXP(M2-M),就不改S(J)了,反之,S(J)=-S(J)(意思就是再改回来),然后M=M2(与下一个M2进行比较),再输出M的值保存

3.第二步进行10000步

结果每次都是90.00000000什么的,数据里全是这个,我吧第二步进行十次,数据一点也没变,求助大家我程序哪里出错了?谢谢大家,感激不尽啊

Text1.f90

671 Bytes, 下载次数: 9

分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

8

帖子

4

主题

0

精华

入门

F 币
41 元
贡献
25 点
沙发
 楼主| 发表于 2014-3-24 20:09:09 | 显示全部楼层
fcode 发表于 2014-3-24 12:42
1.随机数需要种子。建议你加上
2.再随机取一个数X判断,如果X小于EXP(M2-M)这句话里的 X 是随机数,还是  ...

是的,就是单纯的随机数,谢谢你啊

8

帖子

4

主题

0

精华

入门

F 币
41 元
贡献
25 点
板凳
 楼主| 发表于 2014-3-24 20:10:50 | 显示全部楼层
chuxf 发表于 2014-3-24 13:00
我写了个版本,看是否是你需要的。

[mw_shl_code=fortran,true]Program www_fcode_cn

还是谢谢你啊!我可以参考看看,您真是好心人嘞!
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-5-7 07:18

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表