Matlab代码如何用最基本的随机数函数rand产生二项分布B(n, p)的一个随机数呢?先考虑Bernoulli试验,为此产生一个(0,1)上均匀分布随机数,若这个数小于p, 则试验结果记为1,否则记为0,那么试验结果服从0-1分布, n个独立0-1分布随机数的和便是一个二项分布随机数。试根据这样的思路编写B(n, p) 随机数生成函数。
时间: 2024-02-19 09:58:40 浏览: 67
matlab中使用rand函数产生随机数
3星 · 编辑精心推荐
以下是Matlab代码实现B(n, p)随机数生成:
```matlab
function y = binomial(n, p)
count = 0;
for i = 1:n
if rand() < p
count = count + 1;
end
end
y = count;
end
```
其中,`n`是二项分布的参数,即试验次数,`p`是成功的概率。函数通过循环n次,每次生成一个0-1之间的随机数,如果小于p,即认为这次试验成功,计数器加一。最后将计数器的值返回即可。由于每次随机数生成是独立的,因此计数器的值就是n个独立0-1分布随机数的和,满足二项分布B(n, p)的要求。
阅读全文