1.算法描述
LEACH协议，全称是“低功耗自适应集簇分层型协议” (Low Energy Adaptive Clustering Hierarchy)，是一种无线传感器网络路由协议。基于LEACH协议的算法，称为LEACH算法。LEACH算法是一种无线传感器网络路由协议，来源于Wendi Rabiner Heinzelman, Anantha Chandrakasan, 和Hari Balakrishnan三人在2000年Proceedings of the 33rd Hawaii International Conference on System Sciences上的一篇文章Energy-Efficient Communication Protocol for Wireless Microsensor Networks。 该算法基本思想是：以循环的方式随机选择簇头节点，将整个网络的能量负载平均分配到每个传感器节点中，从而达到降低网络能源消耗、提高网络整体生存时间的目的。仿真表明，与一般的平面多跳路由协议和静态分层算法相比，LEACH分簇协议可以将网络生命周期延长15%。
LEACH在运行过程中不断的循环执行簇的重构过程，每个簇重构过程可以用回合的概念来描述。每个回合可以分成两个阶段：簇的建立阶段和传输数据的稳定阶段。为了节省资源开销，稳定阶段的持续时间要大于建立阶段的持续时间。簇的建立过程可分成4个阶段：簇头节点的选择、簇头节点的广播、簇头节点的建立和调度机制的生成。
簇头节点的选择依据网络中所需要的簇头节点总数和迄今为止每个节点已成为簇头节点的次数来决定。具体的选择办法是：每个传感器节点随机选择0-1之间的一个值。如果选定的值小于某一个阈值，那么这个节点成为簇头节点。
选定簇头节点后，通过广播告知整个网络。网络中的其他节点根据接收信息的信号强度决定从属的簇，并通知相应的簇头节点，完成簇的建立。最后，簇头节点采用TDMA方式为簇中每个节点分配向其传递数据的时间点。
稳定阶段中，传感器节点将采集的数据传送到簇头节点。簇头节点对簇中所有节点所采集的数据进行信息融合后再传送给汇聚节点，这是一种较少通信业务量的合理工作模型。稳定阶段持续一段时间后，网络重新进入簇的建立阶段，进行下一回合的簇重构，不断循环，每个簇采用不同的CDMA代码进行通信来减少其他簇内节点的干扰。
LEACH路由协议主要分为两个阶段：即簇建立阶段(setup phase)和稳定运行阶段(ready phase)。簇建立阶段和稳定运行阶段所持续的时间总和为一轮(round)。为减少协议开销，稳定运行阶段的持续时间要长于簇建立阶段。
在簇建立阶段，传感器节点随机生成一个0，1之间的随机数，并且与阈值T(n)做比较，如果小于该阈值，则该节点就会当选为簇头。T(n)按照下列公式计算：式中：P为节点成为簇头节点的百分数，r为当前轮数，G为在最近的1/p轮中未当选簇头的节点集合。簇头节点选定后，广播自己成为簇头的消息，节点根据接收到的消息的强度决定加入哪个簇，并告知相应的簇头，完成簇的建立过程。然后，簇头节点采用TDMA的方式，为簇内成员分配传送数据的时隙。
在稳定阶段，传感器节点将采集的数据传送到簇头节点。簇头节点对采集的数据进行数据融合后再将信息传送给汇聚节点，汇聚节点将数据传送给监控中心来进行数据的处理。稳定阶段持续一段时间后，网络重新进入簇的建立阶段，进行下一轮的簇重建，不断循环。
2.仿真效果预览
matlab2022a仿真结果如下：
3.MATLAB部分代码预览
for i=1:1:n
if Ea>0
p(i)=P*n*S4(i).E*E4(i)/(Et*Ea);
if(S4(i).E>0)
temp_rand=rand;
if ( (S4(i).G)<=0)
%簇头的选举，当选的簇头会把各种相关信存入下面程序所给定的变量中
if(temp_rand<= (p(i)/(1-p(i)*mod(r,round(1/p(i))))))
countCHs4=countCHs4+1;
packets_TO_BS4=packets_TO_BS4+1;
PACKETS_TO_BS4(r+1)=packets_TO_BS4;
S4(i).type='C';
S4(i).G=round(1/p(i))-1;
C4(cluster4).xd=S4(i).xd;
C4(cluster4).yd=S4(i).yd;
distance=sqrt( (S4(i).xd-(S4(n+1).xd) )^2 + (S4(i).yd-(S4(n+1).yd) )^2 );
C4(cluster4).distance=distance;
C4(cluster4).id=i;
X4(cluster4)=S4(i).xd;
Y4(cluster4)=S4(i).yd;
cluster4=cluster4+1;
%计算簇头发送4000bit数据到基站的能量消耗（这里应是所有节点包括簇头每一轮发送4000bit数据）
distance;
if (distance>do)
S4(i).E=S4(i).E- ( (ETX+EDA)*(4000) + Emp*4000*( distance*distance*distance*distance ));
end
if (distance<=do)
S4(i).E=S4(i).E- ( (ETX+EDA)*(4000) + Efs*4000*( distance * distance ));
end
end
end
% S4(i).G=S4(i).G-1;
end
end
end
A40
