三门问题的MonteCarlo仿真方法

一.三门问题

​ 一个人参加某个电视节目，有三个门，其中只有一个门后面有大奖汽车。当该挑战者选好了一个门之后，主持人打开另一个门发现后面没有汽车，问挑战者是否更换门，请问此时这个人是否应该坚持自己最初的想法那个门？

二.问题分析

​ 我们用MonteCarlo法用生成随机数来进行仿真实验。

​ Step1： 首先生成一个随机变量$x_i=j$ ,$j$只能随机取1,2,3之中的一个，代表第$j$个门后有大奖汽车。

​ Step2： 此时再生成一个随机变量$y_i=k$,$k$只能随机取1,2,3之中的一个，代表挑战者心里的想法是$k$门。

​ Step3： 再定义一个随机变量$choic{e}_i=0or1$,代表在主持人询问之后挑战者的决策，$choic{e}_i=0$代表坚持自己 的想法，不改变最初的计划，$choic{e}_i=1$表示改变自己最初的想法。

​ Step4： 现在来判断：如果$x_i=y_i$,即第一次挑战者就选到了正确的门，那么如果此时$choic{e}_i=0$,即坚持自己的 想法，那么此时他坚持自己的想法获胜的次数加一。如果此时$choic{e}_i=1$，即改变了自己的想法，此时 他将挑战失败，那么他失败的次数加一。

​ Step5： 紧接上一步，如果$x_i\ne y_i$,即第一次挑战者没有选到正确的门，那么此时如果$choic{e}_i=0$ ,即坚持自己 的想法，那么他将失败，失败的次数加一。若改变自己的想法，即$choic{e}_i=1$,此时他将成功，那么他不 坚持他自己的想法获胜的次数加一。

​ Step6： 此时$i=i+1$,转向Step1，直到$i\ge n$为止。($n$为总的实验次数。)

​ Step7： 输出：
$$\begin{gathered} P_{赢得奖品|坚持自己的想法}=\frac{n_{坚持他自己的想法获胜的次数}}{n} \\ {P}_{赢得奖品|不坚持自己的想法}=\frac{n_{没有坚持他自己的想法获胜的次数}}{n} \\ {P}_{失败}=\frac{n_{失败的次数}}{n} \end{gathered}$$

三.实现代码

#三门问题
import random,math
numbers = 20000
successfulSelectProp = 0
succUnselectProp = 0
failProp = 0
for i in range(numbers):
    myChoice = random.randint(1,3)
    rightChoice = random.randint(1,3)
    choice = random.randint(0,1)#1代表选择改变，0代表不改变
    if myChoice == rightChoice:
        if choice == 1: #改变了我的意图
            failProp += 1
        else:
            successfulSelectProp += 1
    else:
        if choice == 1:
            succUnselectProp += 1
        else:
            failProp += 1
print('坚持自己的想法最后成功的概率为{:.4f},不坚持想法最后成功的概率{:.4f},最后失败的概率{:.4f}'.format(successfulSelectProp/numbers,succUnselectProp/numbers,failProp/numbers))

四.结果

#结果
坚持自己的想法最后成功的概率为0.1623,不坚持想法最后成功的概率0.3337,最后失败的概率0.5040

Process finished with exit code 0

​ 可以知道在主持人这么问了之后改变自己最初的想法最后成功的概率比较大，几乎是坚持自己的想法获胜的概率的两倍。