应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » matlab做Monte Carlo 如何多次产生 均匀 且 低重复 的随机数？
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1275  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wolffporter

新虫 (初入文坛)

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

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2016-03-10 18:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

manutdhk

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

发自小木虫IOS客户端
赞一下(1人) 回复此楼
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客户端
赞一下(1人) 回复此楼
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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +13 努力奋斗112 2026-04-06 14/700 2026-04-08 15:07 by screening
[考研] 285求调剂 +7 恶法大二的气味� 2026-04-05 10/500 2026-04-08 14:34 by zhq0425
[考研] 化学调剂求助 +11 LULONG1 2026-04-03 13/650 2026-04-08 13:15 by 努力向上，不断�
[考研] 353求调剂 +6 晴空万里air 2026-04-07 6/300 2026-04-08 10:06 by 不吃魚的貓
[考研] 277求调剂 数一104分 +9 瓶子PZ 2026-04-05 14/700 2026-04-07 17:52 by 蓝云思雨
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +9 Mr. Z 2026-04-05 9/450 2026-04-07 15:30 by 尽舜尧1
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 08600生物与医药-327 +9 18755400796 2026-04-05 9/450 2026-04-06 22:35 by 52305043001
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 10/500 2026-04-06 19:47 by jean5056
[考研] 材料334求调剂 +19 Eecho# 2026-04-03 19/950 2026-04-06 08:37 by 小小树2024
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +13 努力奋斗112 2026-04-04 13/650 2026-04-06 07:13 by jj987
[考研] 085600，320分求调剂 +7 大馋小子 2026-04-01 8/400 2026-04-05 21:19 by 学员8dgXkO
[考研] 工科277分求调剂材料 +8 上了上了上哦 2026-04-05 9/450 2026-04-05 13:05 by wwytracy
[考研] 一志愿江南大学085501机械工程专硕326分，本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 材料专业383求调剂 +8 郭阳阳阳成 2026-04-03 8/400 2026-04-04 10:29 by Rednal.
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
信息提示
关闭
请填处理意见
关闭