31无线传感网络布局优化的虚拟力导向粒子群算法MATLAB源代码

1、无线传感网络布局优化的虚拟力导向粒子群算法MATLAB源代码本源代码主要参考了下面的文献：王雪,王晟,马俊杰.无线传感器网络布局的虚拟力导向微粒群优化策略J.电子学报，2007,11:(35)2038-2042.%无线传感器网络布局优化的粒子群导向虚拟力算法仿真程序%GreenSim团队专业级算法设计&代写程序%欢迎访问GreenSim团队主页/greensimclcclearcloseall%ZawFlag=0,0,0,0,0;%障碍物的启用设置，为0表示不启用BorderLength=1000;%正方形区域的边界长度DL=4;%离散化质点的间距，要求DL能够整除BorderLengthSx

2、,Sy=InitSK(BorderLength,DL,ZawFlag);%调用子函数，产生正方形区域，并用离散化质点存储该区域%rLower=70;%传感器节点覆盖半径取值范围的下界rUpper=70;%传感器节点覆盖半径取值范围的上界NodeNum0=50;%设置固定节点个数Cx0,Cy0,r0,R0=InitCovRand(NodeNum0,rLower,rUpper,ZawFlag,BorderLength);%调用子函数，产生固定节点NodeNum1=50;%设置移动节点个数Cx1,Cy1,r1,R1=InitCovRand(NodeNum1,rLower,rUpper,ZawFlag

3、,BorderLength);%调用子函数，产生移动节点figure(1)DrawCovering0(Cx0,Cy0,r0,ZawFlag,BorderLength);title(初始化网络布局(固定节点)figure(2)DrawCovering0(Cx0,Cy0,r0,ZawFlag,BorderLength);holdon4DrawCovering1(Cx1,Cy1,r1,ZawFlag,BorderLength);title(初始化网络布局(固定节点和可移动节点)%Cx=Cx0,Cx1;Cy=Cy0,Cy1;r=r0,r1;R=R0,R1;%如需人为指定个别节点的初始位置，请参考下面的

4、语句进行设置(注意指定的坐标不要落入障碍物内)%Cx(1)=400;Cy(1)=300;%人工指定标号为1的节点的坐标Move=ones(1,NodeNum0+NodeNum1);%节点移动控制参数，等于0表示节点不参与移动第1页欢迎访问GreenSim团队主页/greensim由卩箱：第 页Move(1:NodeNum0)=0;M=200;%迭代次数CC=zeros(1,M+1);%计算初始覆盖率C,IfCovered=CovRates(Cx,Cy,r,Sx,Sy);%计算初始覆盖率CC(1)=C;counter=2;fori=1:MCX,CY=Gravitation(Cx,Cy,r,R,0

5、.2,Sx,Sy,IfCovered);CX,CY=DoNotMove(Cx,Cy,CX,CY,Move);Cx,Cy=Coulomb(CX,CYr,R,0.00005);Cx,Cy=DoNotMove(CX,CYCx,Cy,Move);ifi0.05*Mr,R=ProsperityForAll(Cx,Cy,r,R,0.03,0.1);endC,IfCovered=CovRates(Cx,Cy,r,Sx,Sy);CC(i+1)=C;ifi=2figure(3)DrawCoveringO(Cx(1:NodeNumO),Cy(1:NodeNumO),r(1:NodeNumO),ZawFlag,Bo

6、rderLength);holdonJIDrawCovering1(Cx(NodeNum0+1:end),Cy(NodeNum0+1:end),r(NodeNum0+1:end),ZawFlag,BorderLength);title(粒子群优化网络布局)*1endifi=50figure(4)|TDrawCoveringO(Cx(1:NodeNumO),Cy(1:NodeNumO),r(1:NodeNumO),ZawFlag,BorderLength);holdon*LrIJDrawCovering1(0.25*BorderLength+0.5*Cx(NodeNum0+1:end),0.25

7、*BorderLength+0.5*Cy(NodeNumO+1:end),r(NodeNumO+1:end),ZawFlag,BorderLength);title(虚拟力优化网络布局)endifi=Mfigure(5)DrawCovering0(Cx(1:NodeNum0),Cy(1:NodeNum0),r(1:NodeNum0),ZawFlag,BorderLength);holdonDrawCovering1(Cx(NodeNum0+1:end),Cy(NodeNum0+1:end),r(NodeNum0+1:end),ZawFlag,BorderLength);title(虚拟力导向粒子群优化网络布局)endend%figure(6)plot(CC)xlabel(迭代次数,FontName,TimesNewRoman,FontSize,10);ylabel(覆盖率(%),FontName,TimesNewRoman,FontSize,10);title(虚拟力导向粒子群优化网络布局覆盖率收敛曲线)；源代码运行结果展示Or初始化网络布局（固定节点和可移动节