资源描述：
《31无线传感网络布局优化的虚拟力导向粒子群算法matlab源代码》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、欢迎访问GreenSim团队主页→http://blog.sina.com.cn/greensim邮箱：greensim@163.com无线传感网络布局优化的虚拟力导向粒子群算法MATLAB源代码本源代码主要参考了下面的文献：王雪,王晟,马俊杰.无线传感器网络布局的虚拟力导向微粒群优化策略[J].电子学报,2007,11:(35)2038-2042.%%无线传感器网络布局优化的粒子群导向虚拟力算法仿真程序%GreenSim团队——专业级算法设计&代写程序%欢迎访问GreenSim团队主页→http://blog.sina.com.cn/greensimclcclearcloseall%
2、%ZawFlag=[0,0,0,0,0];%障碍物的启用设置，为0表示不启用BorderLength=1000;%正方形区域的边界长度DL=4;%离散化质点的间距，要求DL能够整除BorderLength[Sx,Sy]=InitSK(BorderLength,DL,ZawFlag);%调用子函数，产生正方形区域，并用离散化质点存储该区域%%rLower=70;%传感器节点覆盖半径取值范围的下界rUpper=70;%传感器节点覆盖半径取值范围的上界NodeNum0=50;%设置固定节点个数[Cx0,Cy0,r0,R0]=InitCovRand(NodeNum0,rLower,rUpper
3、,ZawFlag,BorderLength);%调用子函数，产生固定节点NodeNum1=50;%设置移动节点个数[Cx1,Cy1,r1,R1]=InitCovRand(NodeNum1,rLower,rUpper,ZawFlag,BorderLength);%调用子函数，产生移动节点figure(1)DrawCovering0(Cx0,Cy0,r0,ZawFlag,BorderLength);title('初始化网络布局（固定节点）')figure(2)DrawCovering0(Cx0,Cy0,r0,ZawFlag,BorderLength);holdonDrawCovering1
4、(Cx1,Cy1,r1,ZawFlag,BorderLength);title('初始化网络布局（固定节点和可移动节点）')%%Cx=[Cx0,Cx1];Cy=[Cy0,Cy1];r=[r0,r1];R=[R0,R1];%如需人为指定个别节点的初始位置，请参考下面的语句进行设置（注意指定的坐标不要落入障碍物内）%Cx(1)=400;Cy(1)=300;%人工指定标号为1的节点的坐标Move=ones(1,NodeNum0+NodeNum1);%节点移动控制参数，等于0表示节点不参与移动第5页欢迎访问GreenSim团队主页→http://blog.sina.com.cn/greensi
5、m邮箱：greensim@163.comMove(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:M[CX,CY]=Gravitation(Cx,Cy,r,R,0.2,Sx,Sy,IfCovered);[CX,CY]=DoNotMove(Cx,Cy,CX,CY,Move);[Cx,Cy]=Coulomb(CX,CY,r,R,0.00005);[Cx,Cy]=DoNotMove(CX,CY,Cx,
6、Cy,Move);ifi>0.05*M[r,R]=ProsperityForAll(Cx,Cy,r,R,0.03,0.1);end[C,IfCovered]=CovRates(Cx,Cy,r,Sx,Sy);CC(i+1)=C;ifi==2figure(3)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,Bord
7、erLength);title('粒子群优化网络布局')endifi==50figure(4)DrawCovering0(Cx(1:NodeNum0),Cy(1:NodeNum0),r(1:NodeNum0),ZawFlag,BorderLength);holdonDrawCovering1(0.25*BorderLength+0.5*Cx(NodeNum0+1:end),0.25*BorderLength+0.5*Cy(NodeNum0+1:
