24小时热门版块排行榜    

查看: 1240  |  回复: 4

wolffporter

新虫 (初入文坛)

[求助] matlab做Monte Carlo 如何多次产生 均匀 且 低重复 的随机数?

用 Matlab 模拟聚合反应,有一些反应的概率比其他反应低4个数量级。Matlab自带的 rand、randi、unidrnd 等伪随机数函数以及 拟随机数列 haltonset、sobolset等都无法满足均匀且低重复性的要求,导致分子量分布曲线在高分子量尾部带有一个极高点。不知道大家有没有建议可以解决这个问题?谢谢。
PS:这是短反应时间时的情况,反应时间足够长之后,这个异常便消失。见附图,左上图为段时间;右边两幅为长时间的结果(纵坐标的表示方法不一致而已),其中右上图的纵坐标为对数坐标(与左上图一致),右下图为普通线性坐标。

matlab做Monte Carlo 如何多次产生 均匀 且 低重复 的随机数?
无标题.png
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

manutdhk

木虫 (小有名气)

看不懂你的图。但产生高接收率的随机数是有办法的。一般采用反正切函数导数的形式来逼近你需要抽样的随机数分布。

发自小木虫IOS客户端
2楼2016-03-11 21:11:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolffporter

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by manutdhk at 2016-03-11 21:11:15
看不懂你的图。但产生高接收率的随机数是有办法的。一般采用反正切函数导数的形式来逼近你需要抽样的随机数分布。

我需要平均分布(uniformly distributed)的随机数来抽样反应类型,附图是聚合物的分子量分布曲线,横坐标是分子量的对数,曲线面积正比于聚合物总质量。
3楼2016-03-11 21:17:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manutdhk

木虫 (小有名气)

引用回帖:
3楼: Originally posted by wolffporter at 2016-03-11 21:17:34
我需要平均分布(uniformly distributed)的随机数来抽样反应类型,附图是聚合物的分子量分布曲线,横坐标是分子量的对数,曲线面积正比于聚合物总质量。...

均匀分布是最理想的啊,怎么会出现低接受率的情况呢?

发自小木虫IOS客户端
4楼2016-03-11 21:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wolffporter

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by manutdhk at 2016-03-11 21:20:43
均匀分布是最理想的啊,怎么会出现低接受率的情况呢?
...

不是低接受率的问题,觉得主要还是分布不够均匀的问题,一些增长链没有分配到低概率的反应,而另一些过度得到。如果将随机数对其序号(也就是增长链标号)作图,很明显有空洞,而且这种空洞很难被重复产生的新随机数组覆盖掉,也就意味着这些空洞对应的增长链很难遍历所有类型的反应。
5楼2016-03-12 00:36:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wolffporter 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 303求调剂 +3 睿08 2026-03-17 5/250 2026-03-17 22:01 by 睿08
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 281求调剂(0805) +3 烟汐忆海 2026-03-16 7/350 2026-03-17 20:16 by peike
[考研] 070300化学319求调剂 +4 锦鲤0909 2026-03-17 4/200 2026-03-17 18:21 by 重科小霸王
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +3 小李0854 2026-03-16 3/150 2026-03-17 13:40 by 热情沙漠
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
信息提示
请填处理意见