qos组播路由问题的蚁群算法通用matlab源1.docVIP

QoS组播路由问题的蚁群算法通用Matlab源码(2007-10-29 19:00:40) 知识/探索? function [MRT,EDGES,cost]=ACA_QoS_MR(C,D,S,E,Dmax,K,M,Alpha,Beta,Gamma,Tau,Rho,Q)%% Ant Colony Algorithm for QoS Multicast Routing%? QoS组播路由蚁群算法 %? GreenSim团队原创作品，转载请注明 %? 此源代码完整无删减，请放心使用 %? 欲与原作者技术交流请发邮件aihuacheng@%% 输入参数列表%? C??????????? 费用邻接矩阵（N×N）%? D??????????? 延时邻接矩阵（N×N）%? S??????????? 源节点%? E??????????? 组播目的节点（行向量）%? Dmax???????? 延时约束%? K??????????? 迭代次数（指蚂蚁出动多少波）%? M??????????? 蚂蚁个数（每一波蚂蚁有多少个）%? Alpha??????? 表征信息素重要程度的参数%? Beta???????? 表征启发式因子（费用）重要程度的参数%? Gamma??????? 表征启发式因子（延时）重要程度的参数%? Tau????????? 初始信息素矩阵%? Rho????????? 信息素蒸发系数%? Q??????????? 信息素增加强度系数%% 输出参数列表%? MRT????????? 最优组播树（01邻接矩阵表示）%? EDGES??????? 最优组播树所有的边%? cost???????? 最优组播树的费用%% %% 第一步：变量初始化N=size(C,1);%网络节点个数为NP=length(E);%目的节点个数为MMRT=zeros(N,N);cost=inf;ROUTES=cell(P,K,M);%用细胞结构存储到每一个目的节点的每一代的每一只蚂蚁的爬行路线DELAYS=inf*ones(P,K,M);%用三维数组存储每代每个蚂蚁爬行到各个目的节点的延时COSTS=inf*ones(P,K,M);%用三维数组存储每代每个蚂蚁爬行到各个目的节点的费用%% 第二步：启动到P个目的节点的K轮蚂蚁觅食活动，每轮派出M只蚂蚁for p=1:P??? Tau=ones(N,N);??? for k=1:K??????? for m=1:M%%??????? 第三步：状态初始化??????????? W=S;%当前节点初始化为起始点??????????? Path=S;%爬行路线初始化??????????? PD=0;%爬行路线延时初始化??????????? PC=0;%爬行路线费用初始化??????????? TABU=ones(1,N);%禁忌表初始化??????????? TABU(S)=0;%S已经在初始点了，因此要排除??????????? CC=C;%费用邻接矩阵备份??????????? DD=D;%延时邻接矩阵备份%%??????? 第四步：下一步可以前往的节点??????????? DW=DD(W,:);??????????? DW1=find(DWinf);??????????? for j=1:length(DW1)??????????????? if TABU(DW1(j))==0??????????????????? DW(j)=inf;??????????????? end??????????? end??????????? LJD=find(DWinf);%可选节点集??????????? Len_LJD=length(LJD);%可选节点的个数%%??????? 觅食停止条件：蚂蚁未遇到食物或者陷入死胡同??????????? while (W~=E(p))(Len_LJD=1)%%??????????? 第五步：转轮赌法选择下一步怎么走??????????????? PP=zeros(1,Len_LJD);??????????????? for i=1:Len_LJD??????????????????? PP(i)=(Tau(W,LJD(i))^Alpha)*(C(W,LJD(i))^Beta)*(D(W,LJD(i))^Gamma);??????????????? end??????????????? PP=PP/(sum(PP));%建立概率分布??????????????? Pcum=cumsum(PP);??????????????? Select=find(Pc