蒙提霍尔问题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，验证了蒙提霍尔问题的答案。