1. 无线传感器网络路由协议概述
无线传感器网络（WSN）作为物联网的基础架构之一，其路由协议设计直接影响着网络性能和能耗效率。在众多路由协议中，LEACH（Low-Energy Adaptive Clustering Hierarchy）因其创新的分簇思想成为经典方案，后续衍生的LEACH-C和TS-I-LEACH等改进版本进一步优化了性能表现。
注：本文所有Matlab实现代码均基于R2021a版本测试通过，兼容2016b及以上版本
2. 三种协议的核心原理对比
2.1 经典LEACH协议工作机制
LEACH采用随机轮换簇首机制，每轮通信分为两个阶段：
建立阶段（Setup Phase）：
- 节点以概率p决定是否成为簇首
- 计算公式：T(n) = p/(1-p*(r mod (1/p)))
- 非簇首节点选择信号最强的簇首加入
稳定阶段（Steady Phase）：
- 簇内采用TDMA调度
- 簇首进行数据聚合后转发至基站
- 典型参数：p=0.05，轮次时长20s
2.2 LEACH-C的集中式改进
LEACH-C引入基站集中计算能力：
- 基站收集全网节点位置和能量信息
- 使用模拟退火算法优化簇首选择
- 确保簇首均匀分布且剩余能量较高
- 计算复杂度：O(kN)，k为迭代次数
2.3 TS-I-LEACH的双阶段优化
TS-I-LEACH的创新点在于：
- 两阶段簇首选举：
- 初级选举：类似传统LEACH
- 次级选举：考虑节点度和剩余能量
- 引入阈值函数：
- T_new = T(n)*[1 + α(E_res/E_init) + β(D_avg/D_max)]
- 其中α=0.4，β=0.3为权重系数
3. Matlab实现关键代码解析
3.1 基础参数设置
% 网络参数 nodeNum = 100; % 节点数量 areaSize = 100; % 区域尺寸(m) baseStation = [50,175]; % 基站坐标 % 通信参数 E_elec = 50e-9; % 电路能耗(J/bit) E_fs = 10e-12; % 自由空间模型能耗 E_mp = 0.0013e-12; % 多径衰减模型能耗 d0 = sqrt(E_fs/E_mp); % 临界距离3.2 LEACH簇首选举实现
function isCH = electClusterHead(p, r, G) % p: 预设簇首比例 % r: 当前轮次 % G: 最近1/p轮未当选节点集合 T = p/(1-p*mod(r,round(1/p))); for i = 1:length(G) if rand() < T isCH(G(i)) = true; end end end3.3 LEACH-C的模拟退火优化
function CHs = saOptimization(nodes, k) % 初始化温度 T = 1000; coolingRate = 0.95; currentSolution = randomSolution(nodes); currentCost = calculateCost(currentSolution); for i = 1:k newSolution = perturb(currentSolution); newCost = calculateCost(newSolution); if acceptanceProbability(currentCost, newCost, T) > rand() currentSolution = newSolution; currentCost = newCost; end T = T * coolingRate; end CHs = currentSolution; end4. 性能对比实验设计
4.1 实验环境配置
| 参数 | 值 |
|---|---|
| 节点数量 | 50-200可变 |
| 初始能量 | 0.5J/节点 |
| 数据包大小 | 4000bits |
| 仿真轮次 | 5000 |
| 信道模型 | 混合衰减模型 |
4.2 关键性能指标
网络生命周期：
- 首节点死亡轮次(FND)
- 半数节点死亡轮次(HND)
- 末节点死亡轮次(LND)
能量效率：
- 每轮总能耗
- 能量标准差
吞吐量：
- 成功传输数据包数
- 基站接收速率
5. 实验结果分析与优化建议
5.1 生命周期对比
![生命周期对比曲线]
- TS-I-LEACH的FND比LEACH延长约43%
- LEACH-C在HND指标上表现最优
- 密集网络下TS-I-LEACH优势更明显
5.2 能量消耗分布
% 能量标准差计算代码 energy_std = zeros(1,maxRound); for r = 1:maxRound energy_std(r) = std([nodes.energy]); end- LEACH-C的能量均衡性最佳
- 传统LEACH在后期出现能量空洞
5.3 改进建议
- 动态调整簇首比例：
p_adaptive = p_init * (1 + 0.2*sin(r/100)); - 混合路由策略：
- 初期采用LEACH-C
- 后期切换为TS-I-LEACH
- 引入移动基站支持
6. 常见问题排查指南
6.1 Matlab实现问题
节点不进行簇首选举：
- 检查随机数种子设置
- 验证T(n)阈值计算
- 确保G集合更新正确
能量计算异常：
% 典型能耗计算公式 if d < d0 E_tx = E_elec*k + E_fs*k*d^2; else E_tx = E_elec*k + E_mp*k*d^4; end
6.2 协议参数调优
| 参数 | 推荐范围 | 影响分析 |
|---|---|---|
| p | 0.03-0.1 | 值越小簇首负担越重 |
| 轮次时长 | 15-30秒 | 影响数据采集频率 |
| 数据聚合比 | 0.5-0.8 | 影响簇首计算能耗 |
7. 进阶研究方向
跨层优化设计：
- 结合MAC层调度
- 考虑物理层信道状态
机器学习辅助：
% 使用神经网络预测最佳簇首 net = fitnet(10); net = train(net, inputs, targets);三维空间扩展：
- 无人机辅助网络
- 水下传感器网络应用
提示：完整仿真代码建议采用模块化设计，将协议实现、能量模型、可视化分离为不同m文件
