Leach协议性能分析(2)

最新推荐文章于 2021-04-19 15:28:31 发布

转载最新推荐文章于 2021-04-19 15:28:31 发布 · 1.5k 阅读

文章标签：

#Leach #NS2

网络模拟器(NS2) 专栏收录该内容

89 篇文章

订阅专栏

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

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

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

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

[plain] view plain copy print ?

#计算不同时间剩余节点的数量
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]);
}

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

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 plain copy print ?

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]);
}

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 plain copy print ?

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]);
}

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 plain copy print ?

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

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

得到如下图：