Matlab计算诈金花豹子概率(蒙特卡洛模拟)

x=1:52                            %一共52张牌
y=reshape(x,13,4)
z=y'                                %52张牌分成4种花色，13种点数
s=0
for i=1:1:10000000;      %蒙特卡洛模拟，试验次数取一百万次
a=randperm(52,3);        %随机发三张牌(序号1~52)
b=sort(a);                       %三张牌从小到大排序(本例也可以不排序，但排序对顺子设置条件更有利)
b1=b(1);                         %序号最小一张牌(1~50中的一个数)
b2=b(2);                         %序号中间一张牌(2~51中的一个数)
b3=b(3);                         %序号最大一张牌(3~52中的一个数)
[m1,n1]=find(z==b1);    %第一张牌花色和点数，m代表花色，n代表点数；本例通过Z矩阵进行转化，
[m2,n2]=find(z==b2);    %也可以用mod(b1,4)，mod(b1,13)求余的方法分别获取花色和点数
[m3,n3]=find(z==b3);
if (n1==n2 && n2==n3) s=s+1;     %点数相同  (同花顺，金花，顺子可类似设定条件)，豹子次数加1
end;
end;
s                                      %豹子出现总次数，除以试验次数1