pla算法 matlab,求助Matlab蚁群算法求一般函数极值的算法(函数的维数任意)

function [ROUTES,PL,Tau]=ACASP(G,Tau,K,M,S,E,Alpha,Beta,Rho,Q)

%% ---------------------------------------------------------------

% ACASP.m

%蚁群算法动态寻路算法% ChengAihua,PLA Information Engineering University,ZhengZhou,China

% Email:aihuacheng@gmail.com

% All rights reserved

%% ---------------------------------------------------------------

%输入参数列表% G地形图为01矩阵，如果为1表示障碍物% Tau初始信息素矩阵(认为前面的觅食活动中有残留的信息素)% K迭代次数(指蚂蚁出动多少波)% M蚂蚁个数(每一波蚂蚁有多少个)% S起始点(最短路径的起始点)% E终止点(最短路径的目的点)% Alpha表征信息素重要程度的参数% Beta表征启发式因子重要程度的参数% Rho信息素蒸发系数% Q信息素增加强度系数%

%输出参数列表% ROUTES每一代的每一只蚂蚁的爬行路线% PL每一代的每一只蚂蚁的爬行路线长度% Tau输出动态修正过的信息素%% --------------------变量初始化----------------------------------

%load

D=G2D(G);

N=size(D,1);%N表示问题的规模(象素个数)MM=size(G,1);

a=1;%小方格象素的边长Ex=a*(mod(E,MM)-0.5);%终止点横坐标if Ex==-0.5

Ex=MM-0.5;

end

Ey=a*(MM+0.5-ceil(E/MM));%终止点纵坐标Eta=zeros(1,N);%启发式信息，取为至目标点的直线距离的倒数%下面构造启发式信息矩阵for i=1:N

if ix==-0.5

ix=MM-0.5;

end

iy=a*(MM+0.5-ceil(i/MM));

if i~=E

Eta(1,i)=1/((ix-Ex)^2+(iy-Ey)^2)^0.5;

else

Eta(1,i)=100;

end

end

ROUTES=cell(K,M);%用细胞结构存储每一代的每一只蚂蚁的爬行路线PL=zeros(K,M);%用矩阵存储每一代的每一只蚂蚁的爬行路线长度%% -----------启动K轮蚂蚁觅食活动，每轮派出M只蚂蚁--------------------

for k=1:K

disp(k);

for m=1:M

%%第一步：状态初始化W=S;%当前节点初始化为起始点Path=S;%爬行路线初始化PLkm=0;%爬行路线长度初始化TABUkm=ones(1,N);%禁忌表初始化TABUkm(S)=0;%已经在初始点了，因此要排除DD=D;%邻接矩阵初始化%%第二步：下一步可以前往的节点DW=DD(W,:);

DW1=find(DW

for j=1:length(DW1)

if TABUkm(DW1(j))==0

DW(j)=inf;

end

end

LJD=find(DW

Len_LJD=length(LJD);%可选节点的个数%%觅食停止条件：蚂蚁未遇到食物或者陷入死胡同while W~=E&&Len_LJD>=1

%%第三步：转轮赌法选择下一步怎么走PP=zeros(1,Len_LJD);

for i=1:Len_LJD

PP(i)=(Tau(W,LJD(i))^Alpha)*(Eta(LJD(i))^Beta);

end

PP=PP/(sum(PP));%建立概率分布Pcum=cumsum(PP);

Select=find(Pcum>=rand);

%%第四步：状态更新和记录Path=[Path,to_visit];%路径增加PLkm=PLkm+DD(W,to_visit);%路径长度增加W=to_visit;%蚂蚁移到下一个节点for kk=1:N

if TABUkm(kk)==0

DD(W,kk)=inf;

DD(kk,W)=inf;

end

end

TABUkm(W)=0;%已访问过的节点从禁忌表中删除for j=1:length(DW1)

if TABUkm(DW1(j))==0

DW(j)=inf;

end

end

LJD=find(DW

Len_LJD=length(LJD);%可选节点的个数end

%%第五步：记下每一代每一只蚂蚁的觅食路线和路线长度ROUTES{k,m}=Path;

if Path(end)==E

PL(k,m)=PLkm;

else

PL(k,m)=inf;

end

end

%%第六步：更新信息素Delta_Tau=zeros(N,N);%更新量初始化for m=1:M

if PL(k,m) ROUT=ROUTES{k,m};

TS=length(ROUT)-1;%跳数PL_km=PL(k,m);

for s=1:TS

x=ROUT(s);

Delta_Tau(y,x)=Delta_Tau(y,x)+Q/PL_km;

end

end

end

Tau=(1-Rho).