路径规划

三个城市向一个城市运输资源,途径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();
}
 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值