无线传感器网络DEEC算法.docVIP

无线传感器网络DEEC算法 clear %1.初始参数设定模块 %.传感器节点区域界限(单位 M) xm=100; ym=100; %(1)汇聚节坐标给定 sink.x=0.5*xm; sink.y=0.5*ym; %区域内传器节数 n=100 %簇头优化比例（当选簇头的概率） P=0.1; %能量模型（单位 焦） %初始化能量模型 Eo=0.5; %Eelec=Etx=Erx ETX=50*0.000000001; ERX=50*0.000000001; %Transmit Amplifier types Efs=10*0.000000000001; Emp=0.0013*0.000000000001; %Data Aggregation Energy EDA=5*0.000000001; %高能量节点超出一节点能量的百分比 a=1; %最大循环次数 rmax=5000 %算出参数 do do=sqrt(Efs/Emp); Et=0; %2.无线传感器网络模型产生模块 %构建无线传感器网络,在区域内均匀投放100个节点,并画出图形 for i=1:1:n S(i).xd=rand(1,1)*xm; XR(i)=S(i).xd; S(i).yd=rand(1,1)*ym; YR(i)=S(i).yd; S(i).G=0; S(i).E=Eo*(1+rand*a); E(i)= S(i).E; Et=Et+E(i); %initially there are no cluster heads only nodes S(i).type=N; end d1=0.765*xm/2; K=sqrt(0.5*n*do/pi)*xm/d1^2; d2=xm/sqrt(2*pi*K); Er=4000*(2*n*ETX+n*EDA+K*Emp*d1^4+n*Efs*d2^2); S(n+1).xd=sink.x; S(n+1).yd=sink.y; %3.网络运行模块 %簇头节点数 countCHs=0; cluster=1;%此定义的目的仅仅是给定一个1开始的下标参数，真正的簇头数应该还减去1 flag_first_dead=0; flag_teenth_dead=0; flag_all_dead=0; %死亡节点数 dead=0; first_dead=0; teenth_dead=0; all_dead=0; %活动节点数 allive=n; %counter for bit transmitted to Bases Station and to Cluster Heads packets_TO_BS=0; packets_TO_CH=0; %(1)循环模式设定 for r=0:1:rmax %该 for 循环将下面的所有程序包括在内，直到最后一 end 才结束循环 r %每过一个轮转周期(本程序为10次)使各节点的S（i）.G参数（该参数用于后面的簇选举，在该轮转周期内已当选过簇头的节点不能再当选）恢复为零 if(mod(r, round(1/P) )==0) for i=1:1:n S(i).G=0; S(i).cl=0; end end Ea=Et*(1-r/rmax)/n; %(2)死亡节点检查模块 dead=0; for i=1:1:n %检查有无死亡节点 if (S(i).E=0) dead=dead+1; %(3)第一个死亡节点的产生时间(用轮次表示) %第一个节点死亡时间 if (dead==1) if(flag_first_dead==0) first_dead=r; flag_first_dead=1; end end %10%的节点死亡时间 if(dead==0.1*n) if(flag_teenth_dead==0) teenth_dead=r; flag_teenth_dead=1; end end if(dead==n) if(flag_all_dead==0) all_dead=r; flag_all_dead=1;