【数学建模】轮盘赌算法 |matlab验证

一、轮盘赌算法

轮盘赌算法大量应用于遗传算法，蚁群等算法里。
第一次看的时候硬是没看懂，后面发现是把一个简单的概念弄复杂了。

简单的说，轮盘赌算法就是一个随机算法，就是实现了一些概率问题，其他倒没啥了。核心是存在一个累积概率，这个要理解清楚，我这里就不写了，可以参考论文和其他博客。

看下面这道题就懂了

二、具体实例——TSP旅行商问题选择城市

假设一位商人要去旅游，一共有４座城市可供选择。去城市A的概率为0.27,城市B为0.54.城市C为0.11，城市D为0.08。
用轮盘赌算法验证一下。（虽然感觉像个白痴问题）

%轮盘
list=[0.27 0.54 0.11 0.08];
sum=cumsum(list)
iter=0;iter_max=1000;
a1=0;a2=0;a3=0;a4=0;
while iter<iter_max
    r=rand %生成0-1之间的随机数
    index=find(sum>=rand)%找到比随机数大的那些点，但是只取第一个。index（1）中可以体现
    if list(index(1))== 0.27
        a1=a1+1
    end
    if list(index(1))== 0.54
        a2=a2+1
    end
    if list(index(1))== 0.11
        a3=a3+1
    end 
    if list(index(1))==0.08
        a4=a4+1
    end 
    iter=iter+1;
end
p1=a1./(a1+a2+a3+a4);%  求频率
p2=a2./(a1+a2+a3+a4);
p3=a3./(a1+a2+a3+a4);
p4=a4./(a1+a2+a3+a4);
disp(['0.27出现的次数为：' num2str(a1) ',同时0.27出现的频率为：' num2str(p1)])%空格一定要有
disp(['0.54出现的次数为：' num2str(a2) ',同时0.54出现的频率为：' num2str(p2)])
disp(['0.11出现的次数为：' num2str(a3) ',同时0.11出现的频率为：' num2str(p3)])
disp(['0.08出现的次数为：' num2str(a4) ',同时0.08出现的频率为：' num2str(p4)])

最后结果如下，可以发现，各个城市出现的频率和其本身概率差球不多的，并且也不会因城市的顺序变化而改变。