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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » Matlab 按概率生成随机数的一个小问题
51/1返回列表
查看: 1657  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ice_huya

金虫 (小有名气)

[求助] Matlab 按概率生成随机数的一个小问题 已有2人参与

问题如下：

已知A=[6 5 4 3 2 1], 各元素值对应的概率为P=[0.0049 0.0978 0.1559 0.2220 0.2617 0.1468]
想要生成n个（比如，n=100）服从上述概率的随机数。

通常的办法是用rand函数，但我的问题是，各个值的概率加起来并不是等于1，且概率并没有单调递增或递减，，而rand随机产生的是0～1之间的数，所以想请教大家，该怎么做？

先跪谢！
回复此楼

» 猜你喜欢
1楼 2016-11-18 03:16:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小圆盾

铜虫 (小有名气)

【答案】应助回帖

A = [6 5 4 3 2 1];
P = [0.0049 0.0978 0.1559 0.2220 0.2617 0.1468];
[p, I] = sort(P, 'ascend');
m = length(I);
B = A;
for i = 1:m
    A(i) = B(I(i));
end
n = 100;
p = p./sum(p);
for i = 1:n
r = rand;
if r < p(1)
x(i) = A(1);
elseif r < p(2)
x(i) = A(2);
elseif r < p(3)
x(i) = A(3);
elseif r < p(4)
x(i) = A(4);
elseif r < p(5)
x(i) = A(5);
else
x(i) = A(6);
end
end
x
赞一下(1人) 回复此楼
9楼2016-12-19 10:47:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

ice_huya

金虫 (小有名气)
好像知道了 是用r=a+(b-a)*rand 吗？

还有一个问题，这样生产的随机数，是大致符合给定的概率的分布，而不是严格服从概率吧？每一次重新生成就会变

那是不是越多越准确些呢？比如生产1000个比100个更符合概率分布？
赞一下 回复此楼
2楼2016-11-18 04:30:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
ice_huya: 金币+10, ★有帮助, 非常感谢回复！有启发！！ 2016-11-18 22:31:43
这属于生成离散分布随机数的问题。不知道MATLAB是否有直接生成离散分布随机数的函数，但有一个变通的办法，将[0,1]区间划分为[0.0049 0.0978 0.1559 0.2220 0.2617 0.1468] 6个区间。用MATLAB生成【0，1】上的均匀随机数，落在哪个区间，就取哪个区间对应的整数(即[6,5,4,3,2,1]其中之一)。

第二个问题：MATLAB生成的是伪随机数，不过一般使用足够了。生成的随机样本越多，对概率分布的描述越好，或者说样本分布更加接近总体分布。
赞一下(2人) 回复此楼
知其然，知其所以然。
3楼2016-11-18 14:54:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChuYanWu

金虫 (正式写手)
function Ans = superRand
p = [0.0049 0.0978 0.1559 0.222 0.2617 0.1468];
pz = sum( p );
Num = round( 1000000*p./pz );
Top = sum( Num );
div = cumsum(Num);
Case = randi( Top );
n = find( sign( Case - div )== -1 );
Ans = n(1);

发自小木虫Android客户端
赞一下 回复此楼
4楼2016-11-19 02:46:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错，点赞！ +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题，请指教！ 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[基金申请] 情人节自我反思：在爱情中有过遗憾吗？ +4 瞬息宇宙 2026-02-15 5/250 2026-02-15 22:28 by baiboxie
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
关闭
请填处理意见
关闭