打怪升级的monteCarlo仿真方法

该题目来自公众号数学建模交流的清风老师之手。这个公众号有着学习数学建模的优质资源和课程，强烈推荐学习。

一.问题说明

二.问题的分析

Step1: 初始化等级$degree=1$级，金钱数为$gol{d}_i=0$。

Step2: 产生随机数，将等级$degree$按照一定的概率进行转移到其他$degree$,金币加10000。

Step3:如果$degree\ne 5$,转向Step2。

Step4:将$gol{d}_i$相加,转向Step1,除非$i>n$退出。

输出：$\frac{1}{n}{\sum }_{i=1}^{n}gol{d}_i$。

三.代码实现

lc,clear;
successRate = [0.65 0.2 0.1 0.05 0;0.25 0.4 0.2 0.1 0.05;0.1 0.2 0.4 0.2 0.1;0 0.1 0.3 0.4 0.2];
numbers = 10000;
alphet = [1,2,3,4,5];
sumGold = 0;
x = zeros(1,numbers);
for i = 1:numbers
    degree = 1;%一开始的等级是1级
    gold = 0;%一开始一毛钱都没花
    while degree < 5
        degree = randsrc(1,1,[alphet;successRate(degree,:)]); %依据固定的概率产生字母表中的数。
        gold = gold + 10000;
    end
    sumGold = sumGold + gold;
    x(i) = sumGold/i;
end
plot(1:numbers,x,'b');
xlabel('仿真次数');
ylabel('升到5级的平均花费');
title('打怪升级需要的平均金币')
legend('升到5级别的金币');
disp(['升级到15级平均需要的金钱是：',num2str(x(numbers))]);
  

四.结果

结果是：

  % 升级到15级平均需要的金钱是：159991