在了解了蒙特卡洛模拟的基本操作后,我们再举三个案例来进一步巩固知识。感兴趣的同学可以尝试自己写出代码再来参考笔记。
第一个案例是关于自然常数e的估计问题,这个问题跟Buffon投针有异曲同工之妙;在投针问题里,我们想方设法把pi弄到了概率表达式中,于是同理可知,只要把e设计到概率中去那么这个估计问题也就不难了。(实际上我们已经给出了背景-- 伯努利错装信封问题;要想自己想出一个实验设计还是蛮难的)在写第二题的代码的过程中,我想了半天也不知道该怎么写出“按概率出现的随机数”,搞了几个for循环最后运行量太大导致无功而返。但是后来一查才知道,原来matlab是有现成的函数randsrc()的【捂脸】;那么有了这个强力的randsrc()就好办了,只需要设置一个while循环,对rank=1的初始武器随机升级直到到达5级为止即可。
最后一题也不难,看看就好