回复：=== 蚁群算法源代码大全 ===

39．蚂蚁觅食过程模拟
说明：
1) 各个蚂蚁在没有事先告诉他们食物在什么地方的前提下开始寻找食物。
2) 当一只找到食物以后，它会向环境释放一种信息素，吸引其他的蚂蚁过来，这样越来越多的蚂蚁会找到食物。
3) 有些蚂蚁并没有象其它蚂蚁一样总重复同样的路，他们会另辟蹊径，如果令开辟的道路比原来的其他道路更短，那么，渐渐，更多的蚂蚁被吸引到这条较短的路上来。
4) 最后，经过一段时间运行，可能会出现一条最短的路径被大多数蚂蚁重复着。
题目难度：难
设计要求：请使用C语言编程，设计一个有效的算法，模拟蚂蚁觅食的过程。

# include <iostream.h>
# include <fostream.h>
# include <math.h>
# include <conio.h>
# include <stdlib.h>
# include <iomanip.h>
# define N 20 //city size
# define M 10 //ant number
double inittao=1;
double tao[N][N];
double detatao[N][N];
double distance[N][N];
double yita[N][N];
int tabu[M][M];
int route[M][N];
double solution[M];
int bestroute[N];
double Bestsolution=10000000000;
double alfa,beta,rou,Q;
int Ncmax;
void initparameter();
double Evaluesolution(int *a);
void Incityxy(double x[],double y[]);
void main()
{
int NC=0;
initparameter();
double x[];
double y[];
Incityxy(x,y);
for(int i=0;i<N;i++)
for(int j=i+1;j<N;j++)
{
distance[j][i]=sqrt((x[i]-x[j])*(x[i]-x[j])-(y[i]-y[j])*(y[i]-y[j]));
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]=100/(distance[i][j]+50-int(distance[i][j])%50)
}
for(int 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)%N;
tabu[k][route[k][0]]=1;
}
do {
int s=1;
double partsum;
double pper;
double drand;
while(s,N)
{
for(int k=0;k<M;k++)
{
int jrand=rand()%3000;
drand=double(jrand)/3001;
partsum=0;
pper=0;
for(int 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);
if(pper>drand)
break;
}
tabu[k][j]=1;
route[k][s]=j;
}
s++;
}
for(i=0;i<N;i++)
for(int j=0;j<N;j++)
detatao[i][j]=0;
for(int k=0;k<M;k++)
solution[k]=Evaluesolution(route[k]);
for(k=1;k<M;k++)
if(solution[k]<Bestsolution)
{
Bestsolution=solution[k];
for(s=0;s<N;s++)
Bestsolution[s]=route[k][s];
}
for(k=0;k<M;k++)
{
for(s=0;s<N-1;s++)
detatao[route[k][s]][route[k][s]]+=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(int j=1;j<N;j++)
{
tabu[k][route[k][j]]=0;
route[k][j]=-1;
}
NC++;
} while(NC<Ncmax);
fstream result;
result.open("optimal_results.dat",ios::app);
if(!result)
{
cout<<"can't open the 'optimal_results' file\n";
abort();
}
result<<"************************************"<<endl;
result<<"the initialized parameters of aca are as follows:"<<endl;
result<<"alfa="<<alfa<<",beta="<<beta<<endl;
result<<"Q="<<Q<<", the maximum iteration number of aca is:"<<Ncmax<<endl;
result<<"the evaporation parameter of aca is:"<<rou<<endl;
result<<"**************************************"<<endl;
for(i=0;i<N;i++)
result<<bestroute[i]<<".";
result<<endl;
result.close();
}
double Evaluesolution(int *a)
{
double dist=0;
for(int i=0;i<N;i++)
dist+=distance[*(a+(i+N)%N)][*(a+(i+N+1)%N)];
return dist;
}
void Incityxy(double x[],double y[])
{
fstream inxyfile;
inxyfile.open("city_coordinates.dat",ios::in);
if(!inxyfile)
{
cout<<"can't open the city_coordinates' file\n";
abort();
}
char ch1,ch2;
while(!inxyfile.eof())
{
inxyfile.get(ch1);
if(ch1>='0')
break;
int i=0,j=0;
x[0]=0;y[0]=0;
while(!inxyfile.eof())
{
inxyfile.get(ch1);
if(ch1>='0'&&ch1<='9')
{
ch2=ch1;
while(ch2>='0'&&ch2<='9')
{
switch(i)
{
case 0:break;
case 1:x[j]=x[j]*10+(double(ch2)-48);break;
case 2:y[j]=y[j]*10+(double(ch2)-48);break;
}
inxyfile.get(ch2);
}
i=(++i)%3;
if(i==0&&j<N-1)
{
j++;
x[j]=0;
y[j]=0;
}
}
}
}
}

请问为什么在common.h中声明了
extern double g_Trial[N_CITY_COUNT][N_CITY_COUNT]; //两两城市间信息素
extern double g_Distance[N_CITY_COUNT][N_CITY_COUNT]; //两两城市间距离
extern double g_Prob[N_CITY_COUNT][N_CITY_COUNT];
extern double g_DistanceBeta[N_CITY_COUNT][N_CITY_COUNT];
在common.cpp中还要重新声明
double g_Trial[N_CITY_COUNT][N_CITY_COUNT]; //两两城市间信息素
double g_Distance[N_CITY_COUNT][N_CITY_COUNT]; //两两城市间距离
double g_Prob[N_CITY_COUNT][N_CITY_COUNT];
double g_DistanceBeta[N_CITY_COUNT][N_CITY_COUNT];


我正在看你编写的最大最小蚂蚁算法，有一个地方不明白，在InitData() 中有下面代码
//用贪心算法生成最开始的解
m_cBestAnt.Init();
m_cBestAnt.m_blGreedSearch=true;
m_cBestAnt.Search();
在m_cBestAnt.init()函数中把m_cBestAnt.m_blGreedSearch=false；就是不是用贪心算法。
然后又把m_cBestAnt.m_blGreedSearch=true;
然后m_cBestAnt.Search();中再次调用了init().这不是等于又把m_cBestAnt.m_blGreedSearch=false;了。这样最后还是没有使用贪心搜索啊？
我是初学者，希望大虾指点一二。

以前发的两个问题都已经解决了。现在想问问大家，有没研究过连续域蚁群算法？

这个邮箱登不进去了 请大神赐码 471678756@qq.com

靠，那个家伙把密码改了，龌龊。
我正在修复……

已经恢复了。密码还是 antalgorithm

为什么运行总是出现??? Error using ==> load
Unable to read file data.txt: No such file or directory.
这样的错误啊

楼主大爱啊

楼主有没有关于翼型优化的蚁群算法代码发一份给我
agirlxixi@yahoo.cn
非常感谢！

楼主，我下了你的基本蚁群算法，为什么打开中文都变成乱码呢》？我用vc++6.0打开的。

在物流配送中的应用，我想改成多个配送中心，由于刚刚学编程，楼主能否给点思路？不胜感激