1.两点分布 clc clear a=rand(1,10); for ii=1:10 if a(ii)<0.2 a(ii)=0; else a(ii)=1; end end a x=0的概率为0.2,x=1的概率为0.8:两点分布 clc clear a=rand(1,10); b=(a>0.2) 循环用向量化表示 2.伯努利分布(二项分布) clc clear N=1000000; r=binornd(19,0.3,1,N);%A事件发生概率0.3,重复19次.事件A发生的次数.仿真1000…
原文地址:<如何做Gibbs采样(how to do gibbs-sampling)> 随机模拟 随机模拟(或者统计模拟)方法最早有数学家乌拉姆提出,又称做蒙特卡洛方法.蒙特卡洛是一个著名的赌场,赌博总是和统计有着密切的关系,所以这个命名风趣而贴切,被广为接受. 随机模拟的一个重要问题就是给定一个概率分布\(p(x)\),如何生成它的样本.均匀分布\(Uniform(0,1)\)的样本可以通过线性同余发生器生成的伪随机数来模拟.常见的概率分布,无论是离散的还是连续的分布,都可以基于\(Unif…