leach和leach-c协议仿真

1.复制leach_test为leach-c_test，修改里面的文件夹和输出文件名。并且手动建立相应的文件夹。

很多教程说只修改文件名，没提到修改文件夹，如果同时运行两个协议会覆盖一部分实验结果的，而且是并行的就更难分解出是哪个协议的实验数据了。所以还是分开较好。

2.添加到test文件最后，后台运行的，需要等待。

3.写脚本分析实验数据。

[plain]
view plaincopy

#计算不同时间剩余节点的数量

BEGIN {

countcyl=0;
totalleft=0;

lasttime=0;
time[0]=0;
node=0;
total[0]=100;
}
{

simtime = $1;
nodeid = $2;
statenode = $3;

if (simtime>lasttime ) {
countcyl++;
lasttime=simtime;
time[countcyl]=simtime;
totalleft=0;
}

if (statenode==1)
totalleft++;
total[countcyl]=totalleft;
}
END {
for(i=0;i<=countcyl;i++)
printf( "%f %d\n",time[i],total[i]);
}

[plain]
view plaincopy

BEGIN {

countcyl=0;
packetsum=0; #当前时间发包且接收的总数

lasttime=0;
time[0]=0;
sum[0]=0;
}
{

simtime = $1;
nodeid = $2;
packet = $3; #读取当前时间，当前节点的发包且被成功接收的数目

if (simtime>lasttime ) { #具体算法
packtsum=0;
countcyl++;
lasttime=simtime;
time[countcyl]=simtime;
}

if (simtime==lasttime ) {
packetsum=packetsum+packet;
sum[countcyl]=packetsum;
}
}
END { #对应时段，打印发包总量。
for(i=0;i<=countcyl;i++)
printf ( "%f %f\n",time[i],sum[i]);
}

[plain]
view plaincopy

BEGIN {

countcyl=0; #计数器，用来记录当前的轮数 或者当前的时段数
energysum=0; #用来暂时储存当前时段的节点消耗能量总和

lasttime=0; #记录当前最后的时间
time[0]=0; #存储各（轮）分段的时间
sum[0]=0; #存储各分段时间对应的消耗的能量
}
{

simtime = $1; #文件中第一字段的值，当前时间
nodeid = $2; #。。。第二字段的值，节点ID
nodeenergy = $3; #。。。第三字段的值，当前节点使用的能量

if (simtime>lasttime ) { #这里就是一个简单的方法来提取每个时间段，节点总共消耗了多少能量
energysum=0;
countcyl++;
lasttime=simtime;
time[countcyl]=simtime;
}

if (simtime==lasttime ) {
if (nodeenergy<2.0) {
energysum=energysum+nodeenergy;
sum[countcyl]=energysum;
}
else if (nodeenergy>=2.0) { #节点初始的总能量为2
energysum=energysum+2.0;
sum[countcyl]=energysum;
}
}
}
END { #对应时间段，得出能量消耗量并打印
for(i=0;i<=countcyl;i++)
printf ( "%f %f\n",time[i],sum[i]);
}

4.gnuplot来画图

[plain]
view plaincopy

set multiplot
set origin 0.0,0.5
set size 0.5,0.5
plot 'leach.alive.rst' with linespoint,'leach-c.alive.rst' with linespoint
set origin 0.5,0.5
set size 0.5,0.5
plot 'leach.data.rst' with linespoint,'leach-c.data.rst' with linespoint
set origin 0.33,0.0
set size 0.5,0.5
plot 'leach.energy.rst' with linespoint,'leach-c.energy.rst' with linespoint

得到如下图：



明显实验结果说明leach-c死亡节点比较早，虽然bs收到数据量大，但是能量消耗是leach协议更关心的问题。

修改bs坐标为(0,0)，再次计算得到如下结果：



从各方面来说c都比leach更优，有人说若bs在节点区域内部，则leach更优，若在外部则leach-c更优！

想知道为什么还需要看懂这两个协议喽。

接下来代码和理论分析。