jichu1.rar_ant system_分簇 MATLAB_分簇算法_蚁群 无线_蚁群分簇

clear all clear clc M=82; A=[0;0]; B=[10;10]; rand('state',1) G=rand(2,80)*10; U=[A G B]; C=U'; K=10;Alpha=1;Beta=5;Rho=0.5;Q=100; N=size(C,1);D=ones(N,N); Length_best=ones(K,1); Routh_best=cell(K,1); Length_average=ones(K,1); for i=1:N for j=1:N if i < j D(i,j)=sqrt((C(i,1)-C(j,1))^2+(C(i,2)-C(j,2))^2); end D(j,i)=D(i,j); end end Eta=ones(N,N); D_average=mean(mean(D)); R_max=4*D_average; E=2*ones(N,1); E1=zeros(N,1); E2=zeros(N,1); E3=2*ones(N,1); kd=2000; Eelec=50*10.^-9; Eamp=100*10.^-12; for i=1:N for j=1:N Eta(i,j)=1/D(i,j); end end Tau = ones(N,N); Routes=cell(K,M); PL=zeros(K,M); for k=1:K disp(k); for m=1:M W=1; Path=1; PLkm=0; Tabukm=ones(1,N); Tabukm(1)=0; DD=D; DW=DD(W,:); while W~=N for j=1:N if Tabukm(j)==0 DD(W,j)=inf; end end LJD=find(DD(W,:)<R_max); Len_LJD=length(LJD); if Len_LJD>=1 PP=zeros(1,Len_LJD); for i=1:Len_LJD PP(i)=(Tau(W,LJD(i))^Alpha)*(Eta(W,LJD(i))^Beta); end PP=PP/(sum(PP)); Pcum=cumsum(PP); Select=find(Pcum>=rand); to_visit=LJD(Select(1)); Path =[Path ,to_visit]; PLkm=PLkm+DD(W,to_visit); E1(W)=E1(W)+Eelec*kd+Eamp*kd*(DD(W,to_visit)/10)^2; E2(to_visit)=E2(to_visit)+Eelec*kd; E3=E-E1-E2; Tau(W,to_visit)=(1-Rho)*Tau(W,to_visit)+Q/D(W,to_visit); W=to_visit; Tabukm(W)=0; else break ; end end Routes{k,m}=Path; if Path(end)==N PL(k,m)=PLkm; else PL(k,m)=0; 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); end end if k>1 PL(k,1)=Length_best(k-1); end Length_best(k)=min(PL(k,:)); Position=find(PL(k,:)==Length_best(k)); Routh_best{k}=Routes{k,Position(1)}; Length_average(k)=mean(PL(k,:)); for i=1:length(Position) p=Position(i); Route=Routes{k,p}; for j=2:length(Route) Delta_Tau(Route(j),Route(j-1))=Q/Length_best(k); end end Tau=Tau+Delta_Tau; end figure(1) set(gcf,'Name','Figure of length_best','Color','r') plot(Length_best ,'k'); hold on; figure(2) set(gcf,'Name','Figure of length_average','Color','k') plot(Length_average,'r'); position=find(Length_best==min(Length_best)); shortest_path=Routh_best{position(1)}; shortest_length=Length_best(position(1)); figure(3) set(gcf,'Name','Figure of shortest_path','Color','r') N=length(shortest_path); scatter(C(:,1),C(:,2),10,'filled') hold on for i=2:N plot([C(shortest_path(i-1),1),C(shortest_path(i),1)],[C(shortest_path(i-1),2),C(shortest_path(i),2)]); hold on end Eaverage=mean(mean(E1+E2))