轮盘赌算法

轮盘赌长这个样子,每个格子的概率是1/37,我们需要用到的模型如右边这个图,即每个有颜色格子的概率是不同的,整体概率为1。


先撇开遗传算法,觉得上来讲染色体群体的选择 有点不地道。通俗的讲一下我对轮盘赌算法的理解。
右上边饼图不同颜色的区域,面积大小对应着不同的概率,面积越大,代表概率越大。假想把这张图打印到一张纸上,随机扔一把小米,落在3区域的小米相对来说数量最多。好了,现在我一粒一粒的扔,扔了10粒米(意味着只选了10个样本),假如5个落在3区域,3个落在1区域,1个落在4区域,1个落在5区域。


在应用中,比方说,7号米粒利用概率38%(因为落在了3号区域),8号米粒利用概率14%,9号米粒利用概率38%,10号米粒利用概率31%。
这样有什么好处?避免了所有的米粒都选择概率最大的区域3(所谓的最优值问题),换句话,各个概率(各种情况)都相应的被使用到了,避免了陷入局部最优的问题。
问题又来了,实际编程中该怎么用?往下看。


把概率整合到一条线上,然后随机产生数据a,a属于[0,1],比方说



这样,不同的数据对应不同的数据概率,并且整体上还保留了“区域概率越大,对应数据越多”这一分布!

附上代码吧,c语言写的

<span style="font-size:18px;">#include<math.h>
#include<stdio.h>
#include<stdlib.h>

float generate_random()
{
	//randomly generate number between [0,1]
	float rand_num = 0.0;
	rand_num = (float)rand()/RAND_MAX;
	//printf("%f\n", rand_num);
	return rand_num;
}

float *generate_probability_band(float *prob_arr, int length)
{
	//产生概率带
	//传入参数是在main函数中指定的概率分布
	//传出参数就是概率带(积累概率)
	//计算累计概率,保存在sum_array数组中,用malloc写了一下。
	
	printf("一共有%d个概率数据\n",length);

	float *sum_array = NULL;
	float *sum_array_tmp = NULL;
	sum_array = (float *)malloc((length+1)*sizeof(float));//多一个是为了保存概率带最左边那个0
	sum_array_tmp = sum_array;
	int i = 0;
	float sum = 0.0;
	sum_array[0] = 0.0;
	for(i;i<length;i++)
	{
		sum = sum + prob_arr[i];
		printf("sum=%f\n",sum);
		sum_array++;//先执行++,目的是把概率带最左边那个0保存在sum_array[0]上
		*sum_array = sum;
		//printf("%x-->%f\n", sum_array,*sum_array);
		
	}
	//sum_array = sum_array_tmp;
	return sum_array_tmp;
	//free(sum_array);
}

//判断随机数位于概率带的哪个位置
int *judge_random_location(float *sum_array, int length)
{
	
	int i = 0;
	int j = 0;
	float rand_num=0;
	int *count=NULL;
	count = (int *)malloc((length-1)*sizeof(int)); //count数组用来统计 落在某个概率上的数量
	for(i=0;i<length-1;i++)
	{
		//初始化为0
		count[i] = 0;
	}

	for(i=0;i<length;i++)
	{
		printf("sum_array[%d] = %f\n", i, sum_array[i]);
	}

	for(j=0;j<1000;j++) //产生1000个随机数
	{
		rand_num = generate_random();
		if(rand_num>0 && rand_num<1)
		{
			for(i=0;i<length;i++)
			{
				if(rand_num>sum_array[i] && rand_num<=sum_array[i+1])
				{
					//对应的概率带 计数器加1
					count[i] = count[i]+1;
				}

			}
		}
	}
	
	return count;
}
void main()
{
	int i=0;

	//probability array
	float prob_arr[] = {0.1, 0.2, 0.3, 0.4}; //按照概率1:2:3:4
	int length = sizeof(prob_arr)/sizeof(prob_arr[0]);
	float *sum_array = NULL;
	int *count=NULL;

	sum_array = generate_probability_band(prob_arr, length);
	count = judge_random_location(sum_array, length+1);

	for(i=0;i<length;i++)
	{
		printf("%d:", count[i]);
		if(i == (length-1))
			printf("\b");
	}
	printf("\n");
	system("pause");
}</span>


  • 27
    点赞
  • 77
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值