24小时热门版块排行榜    

查看: 1244  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wolffporter

新虫 (初入文坛)

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

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

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

» 猜你喜欢

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

已阅   回复此楼   关注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的回帖
查看全部 5 个回答

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的回帖

wolffporter

新虫 (初入文坛)

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

不是低接受率的问题,觉得主要还是分布不够均匀的问题,一些增长链没有分配到低概率的反应,而另一些过度得到。如果将随机数对其序号(也就是增长链标号)作图,很明显有空洞,而且这种空洞很难被重复产生的新随机数组覆盖掉,也就意味着这些空洞对应的增长链很难遍历所有类型的反应。
5楼2016-03-12 00:36:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302求调剂 +6 呼呼呼。。。。 2026-03-17 6/300 2026-03-17 23:56 by 星空星月
[考研] 311求调剂 +9 冬十三 2026-03-15 10/500 2026-03-17 22:53 by lbsjt
[考研] 085601专硕,总分342求调剂,地区不限 +4 share_joy 2026-03-16 4/200 2026-03-17 21:32 by hmn_wj
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 0703化学调剂 ,六级已过,有科研经历 +8 曦熙兮 2026-03-15 8/400 2026-03-17 20:31 by xilongliang
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 293求调剂 +6 世界首富 2026-03-11 6/300 2026-03-17 17:04 by ruiyingmiao
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 材料与化工304求B区调剂 +7 邱gl 2026-03-11 8/400 2026-03-17 09:36 by 努力学习赚彩礼
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
信息提示
请填处理意见