三个城市向一个城市运输资源,途径24个城市,用蚁群算法解决路径的选择的问题,大家谁知道》》谁感兴趣》》谁帮帮小弟,程序有一个,没研究透彻呢》》哎》》欢迎讨论哈!!!
#include <math.h>
#include <time.h>
#include <stdlib.h>
#include <stdio.h>
#define N 6 /*城市的数目*/
#define M 50 /*蚂蚁的数目*/
double inittao=1; /*初始信息素量*/
double tao[N][N]; /*每条路径的信息素量*/
double detatao[N][N]; /*相应路径的信息素增量*/
double distance[N][N]; /*保存城市距离的数组*/
double yita[N][N]; /*启发函数,一般初始其值为*/
int tabu[M][N]; /*tabu[i][j]=1表示蚂蚁i已经走过了j城市*/
int route[M][N]; /*保存蚂蚁k的路径的数组为route[k][N] */
int i,j,k; /*循环变量*/
double solution[M]; /*蚂蚁k访问的路径的长度为solution[k]*/
int BestRoute[N]; /*当前最优的访问路径序列*/
double BestSolution=10000000; /*当前最优的访问路的长度*/
int aa=1;
int NcMax;
double alfa,beta,rou,Q;
/* alfa是信息启发式因子,表示轨迹的相对重要性,其值越大,则该蚂蚁越倾向于选择其他,蚂蚁经过的路径
蚂蚁之间的协作性越强。beta是期望启发式因子,表示能见度?其值越大,则该状态转移概率越接近于贪心规则?
rou是信息残留因子,Q为信息素强度,用于计算蚂蚁留在路径上的信息量 */
int coordinate[N][2]; /*城市的坐标,由用户输入*/
int jianyan(int a[],int n) /*用于检验用户输入的路径序列数组是否合法*/
{
for(i=0;i<n-1;i++){ /*判断数组的元素是否有重复*/
for(j=i+1;j<n;j++){
if(a[i] == a[j]){
k=0;
break;}
else
k=1;
}
if(j<n)break;
}
j=1;
for(i=0;i<n;i++) /*判断数组的元素是否大于等于城市数*/
if(a[i]>=n){
j=0;
break;}
return k*j;
}
int zjianyan(int a[][2],int n) /*判断用户输入的城市坐标是否有相同*/
{
for(i=0;i<n-1;i++)
{
for(j=i+1;j<n;j++)
{
if(a[i][0] == a[j][0]&&a[i][1] == a[j][1])
{
k=0;
break;}
else
k=1;
}
if(j<n)break;
}
return k;
}
void Intial()
{
do{
printf("please input the coordinate of cities:/n");
for(i = 0;i < N;i++)
for(j = 0;j < 2;j++)
scanf("%d",&coordinate[i][j]);
aa=zjianyan(coordinate,N);
if(aa==0)
printf("Error!Please input again!/n");
}while(aa==0);
}
void initparameter(void) /*初始化各参数*/
{
alfa=1.0; beta=5.0; rou=0.9; Q=100;
NcMax=250;
}
double EvalueSolution(int a[])/*评估路径,测试路径长度*/
{
double dist=0;
for(i=0;i<N-1;i++)
dist+=distance[a[i]][a[i+1]];
dist+=distance[a[i]][a[0]]; /*循环到最后要返回起点,形成回路*/
return dist;
}
void dotest(void) /*测试结果准确性*/
{
int rp;
int test[N];
do{
rp = 0;
printf("make a test: /n");
do{
for(i=0;i<N;i++)
scanf("%d",&test[i]);
aa=jianyan(test,N);
if(aa==0)
printf("Error!Please input again!/n");
}
while(aa==0);
printf("the result of test is %4f /n",EvalueSolution(test));
printf("would you want to continue?(1 for yes and 0 for no):");
scanf("%d",&rp);
}while(rp==1);
}
void main()
{
int NC=0;
Intial();
initparameter();
for(i=0;i<N;i++){
for(j=i+1;j<N;j++)
{
distance[j][i]=sqrt(pow(coordinate[i][0]-coordinate[j][0],2)+pow(coordinate[i][1]-coordinate[j][1],2));
distance[i][j]=distance[j][i];
}
}
for(i=0;i<N;i++){
for(j=0;j<N;j++)
{
tao[i][j]=inittao;
if(j!=i)
yita[i][j]=1/distance[i][j];
}
}
for(k=0;k<M;k++)
for(i=0;i<N;i++)
route[k][i]=-1;
srand(time(NULL));
for(k=0;k<M;k++)
{
route[k][0]=k%N; /*设置初始城市*/
tabu[k][route[k][0]]=1; /*设置初始城市被访问标记*/
}
do {
int s=1;
double partsum; /*为方便计算转移概率而设的变量,本身没有实际意义*/
double pper; /*转移概率*/
double drand;
while(s<N)
{
for(k=0;k<M;k++)
{
int jrand=rand()%3000;
drand=jrand/3001.0;
partsum=0;
pper=0;
for(j=0;j<N;j++) /*根据概率函数计算蚂蚁的转移概率 */
{
if(tabu[k][j]==0)
partsum+=pow(tao[route[k][s-1]][j],alfa)*pow(yita[route[k][s-1]][j],beta);
/*route[k][s-1]表示蚂蚁k第s-1步走到的城市,s==1表示蚂蚁从初始城市出发
tao[route[k][s-1]][j]表示蚂蚁k经过的前一个城市到下一个没被访问的城市j的信息量*/
}
for(j=0;j<N;j++)
{
if(tabu[k][j]==0)
pper+=pow(tao[route[k][s-1]][j],alfa)*pow(yita[route[k][s-1]][j],beta)/partsum;
if(pper>drand) /*当drand落在第j个城市上时,选择j城市 */
break;
}
tabu[k][j]=1; /*禁忌表置访问标志 */
route[k][s]=j; /*保存蚂蚁k的第s步经过的城市 */
}
s++;
}
for(i=0;i<N;i++)
for(j=0;j<N;j++)
detatao[i][j]=0;
for(k=0;k<M;k++) /*计算每个蚂蚁经过路径的长度,保存到当前代为止最好的解路径及其长度 */
{
solution[k]=EvalueSolution(route[k]);
if(solution[k]<BestSolution)
{
BestSolution=solution[k];
for(s=0;s<N;s++)
BestRoute[s]=route[k][s];
}
}
for(k=0;k<M;k++) /*计算各个路径上的信息素增量*/
{
for(s=0;s<N-1;s++)
detatao[route[k][s]][route[k][s+1]]+=Q/solution[k];
detatao[route[k][N-1]][route[k][0]]+=Q/solution[k];
}
for(i=0;i<N;i++) /*更新路径上的信息素 */
for(j=0;j<N;j++)
{
tao[i][j]=rou*tao[i][j]+detatao[i][j];
if(tao[i][j]<0.00001)
tao[i][j]=0.00001;
if(tao[i][j]>20)
tao[i][j]=20;
}
for(k=0;k<M;k++) /*将蚂蚁的路径再重新置空,为下一次循环做准备*/
for(j=1;j<N;j++) /*注意起始城市,即j=0没有被清空 */
{
tabu[k][route[k][j]]=0;
route[k][j]=-1;
}
NC++;
} while(NC<NcMax);
printf("*-------------------------------------------------------------------------*/n");
printf("the initialized parameters of ACA are as follows:/n");
printf("alfa=%.1f, beta=%.1f, rou=%.1f, Q=%.1f /n",alfa,beta,rou,Q);
printf("the best route is:");
for(i=0;i<N;i++)
printf("%d ",BestRoute[i]);
printf("/n");
printf("the shortest length of the path is:%2f /n",BestSolution);
printf("*-------------------------------------------------------------------------*/n");
system("PAUSE");
dotest();
}