CN101959244B - 一种适用于无线传感网的层次型路由控制方法 - Google Patents

Abstract

一种适用于无线传感网的层次型路由控制方法，包括以下步骤：1)对无线传感节点随机部署形成的网络拓扑进行距离侦测存储；2)基于N阶近邻理论，在传感节点位置及网络拓扑动态变化时，通过对节点分布离散度及曲率变化的分析确定当前最优的聚簇数量k _opt ；3)根据最优簇数k _opt 计算最佳簇头占有比率p _opt ，并结合节点剩余能量选取簇头集；4)引入Hausdorff距离调整初始快速形成的簇结构；5)数据传输，将各传感器节点收集的数据信息通过簇头发送至基站。本发明提供准确的最优簇数计算，簇头分布更均匀，簇结构合理，节省了各节点能耗，实现了无线传感网的网络负载均衡。

Description

一种适用于无线传感网的层次型路由控制方法

技术领域

本发明涉及无线传感网领域，尤其涉及的是一种适用于无线传感网层次型路由控制方法。

背景技术

无线传感网是一类新兴的科学领域，也是现今物联网的重要组成部分。无线传感网本身是由大规模低成本、低功耗的多功能微型传感器节点以无线方式组成的多跳自组织网络，并隶属于无线网络的范畴。但区别于一般的无线网络，通常无线传感网在节点部署后即保持静止状态，除非节点损坏或者死亡，否则网络拓扑不会改变；网络自身只关心数据的收发，而发送与接收者并不是网络关心的重点。因此，传统基于IP协议的路由机制以及Ad hoc网络中针对移动性较强的多跳路由并不适合于无线传感网。

针对无线传感网的特性，现有的路由按照网络结构可大致被分为三类：平面路由、基于地理位置的路由和层次型路由，而层次型路由被认为是最适合解决无线传感网的路由问题，其中尤以Heinzelman等提出的低能耗自适应等级分簇路由（LEACH）为代表。LEACH以轮转机制周期性地执行簇建立和稳定数据通信过程，等概率随机选举簇头并履行数据融合策略，在一定程度上降低了数据的通信量，有效减少能量消耗，延长了网络的生存时间。

LEACH引入的簇头数据融合机制降低了数据传输的冗余度，改变了传统短距离多跳路由存在的问题，但也带来了新的问题：1）未将节点的剩余能量纳入考虑范畴，导致部分节点很快耗尽能量而缩短网络生存时间；2）簇头选举分布不均，未考虑节点的具体分布，从而导致网络能耗的不均；3）随机选举簇头及成簇方式容易带来大量不确定因素，导致当前的通信网络并非最优，进而造成生存时间大幅缩减；4）伴随网络规模的扩大，容易出现能耗不均的现象（离簇头越远消耗的能量越多，反之越则少），从而进一步影响了网络的整体覆盖和生存时间。

LEACH-C、LEACH-F两种涵盖集中式簇形成算法的路由解决了上述的部分问题，但是LEACH-C过分依赖于GPS与查找算法，效率比LEACH有所降低，而LEACH-F则依赖簇数量的设置以及簇大小的能量成本，一旦网络发生改变（新节点的加入或者节点耗尽能量死亡），其性能指标就迅速下降。为此，Lindsey等提出的PEGASIS协议只保留一个簇头，利用贪婪策略按照离基站的远近及相邻关系成链，通过令牌控制方式沿链以最小功率传输融合后的数据，并最终由簇头把融合后的数据发送给基站。PEGASIS避免了LEACH协议因频繁选举簇头、重构簇带来的通信开销，较之LEACH协议提高了网络的生存时间，但也存在诸如单簇遇到网络失败的机率大大提高，传输时延增大等问题。Manjeshwar等提出适用于响应型WSNs的TEEN协议，该协议通过设立硬、软两阀值来调节数据的发送量，达到精度与系统能耗之间的平衡，实现对突发事件及热点区域的监视，但也存在阀值设置敏感，容易产生噪声干扰，不适用于周期性的数据上报传输。Younis等提出的HEED协议将节点剩余能量纳入考虑范畴，并依据主从关系引入多个约束条件作用于簇头选举，同时首次提出将簇内平均可达能量（AMRP）作为衡量簇内通信成本标准的理念，较之LEACH成簇速度快、网路拓扑更合理，能耗更均匀，但却无法保证所选簇头是当前最优（节点能量最多），即簇头选举的可靠性相对不高。鉴于以上分析，迫切需要一种更适用于无线传感网的路由方法来解决一些实际的节点数据路由传输问题。

发明内容

为了克服现有无线传感节点的路由控制方法的簇头分布不均、簇间距不均、簇头选举不合适、能耗负载不均衡的不足，本发明提供了一种簇头分布均匀、簇头选举更合适、簇间距更合理、能耗负载均衡的适用于无线传感器网络的层次型路由控制方法。

本发明解决其技术问题所采用的技术方案是：

一种适用于无线传感网的层次型路由控制方法，包括以下步骤：

1) 节点距离侦测存储：基站首先向各节点广播侦测请求包，节点在接收请求数据包后向基站传回自身的信息；基站收到回复信息后计算各自的距离并形成路由距离矩阵表，同时存储并转发相应的路由节点距离矩阵表给各自的节点，即节点只存储与自身相关的距离信息；

2) 最优聚簇分析：设定在指定区域内部署的N个无线传感节点组成的集合表示为S={x _i | x _i  =(a _i ,b _i )且x _i ∈R ², i = 1,2,…,N }，令                                                

为节点x _i 的一阶最短近邻距离，

为节点x _i 的二阶最短近邻距离，以此类推，

为节点x _i 的N-1阶最短近邻距离，则按照节点x _i 的最短近邻距离从小到大的排列有：

                    (1)

进一步，分别计算j阶最短近邻距离的平均值

和平方均值：

                      (2)

                     (3)

求解方差

：

                (4)

假定y _i 为原传感节点x _i 的次序排列，且满足y ₁ ≤ y ₂ … ≤ y _N-1的关系，同时令

为其对应独立的且具有相同概率的概率密度函数，且必须满足；设定概率密度函数f(y) = 1，运用Beta函数的特性将式(4)推导转换为式(19)：

                 (19)

系数

、分别对应于：

 ,     

          (20)

设

，在传感网络节点曲率变化时，如果在某一节点上曲率变化满足

                 (22)

则认为在该传感网络节点上发生了跃变，该传感网络节点也被称为跃变节点，记为m ^*，其中阀值

是在(0,1)区间中的足够小的常数；

    假定网络初始簇数为k=1，当计算某一节点曲率时满足式(22)，簇数在原来的基础上加1，迭代遍历节点集，统计所有的跃变节点m ^*，得到最优的簇数量k _opt ；

3) 簇头选举：阀值T(n)决定当前传感节点成为簇头节点的概率，而簇头在所有节点中所占百分比p则决定着T(n)大小，通过所述对传感节点的最优聚簇分析得到当前最优的簇数量k _opt ，进而计算当前最优的簇头占用比p _opt ：

                            (23)

则T(n)表示为：

                (24)

其中，r为执行的轮数，N为传感节点集的个数，G为本周期内未当选为簇头的节点集合，当r到达一个周期执行时间，未当选为簇头的节点集合G将被重置；

4) 簇形成：各簇头节点及时向整个网络通告自己当选的消息，该消息包内容只包含簇头节点的ID信息，同时簇头节点等待其他非簇头节点请求加入，其他非簇头节点直接根据初始状态下基站测定并转发的近似距离存储表迅速选择加入最近的簇头，并向簇头发回加入请求消息，该消息包含当前的非簇头节点ID以及欲加入的请求指令；

5) 数据传输：将各传感节点收集的数据信息通过簇头进一步发送至基站。

作为优选的一种方案：在所述步骤3)中，将节点剩余能量纳入簇头选取评价标准，则公式（24）转变为式（25）：

           (25)

其中E _residual 为节点当前评估的剩余能量，E _max 则为节点的最大能量，即为节点初始能量。

作为优选的另一种方案：在所述步骤4)中，假设N个传感节点围绕k个簇头初步形成的簇集为Clusters={C ₁,C ₂,…,C _k }，节点v ₁

C _i 、v ₂

C _j ，且满足v ₁

v ₂、C _i

C _j ，1 i

k,1

j

k，则簇C _i 到C _j 的Hausdorff距离为：

                  (26)

其中的为簇C _i 中的节点到C _j 中节点的范数，如L ₂范数；

最终两簇之间的Hausdorff距离用下式来表示：

 

             (27)

在确定簇间的Hausdorff距离后，假设簇C _i 与C _j 的Hausdorff距离H(C _i ,C _j ) = R _h ，那么簇C _i 中每个非簇头节点到C _j 的节点必须满足d(v) < R _h ，v代表簇C _i 中任一非簇头节点，反之亦然。当C _i 中存在一个或者多个节点大于或者等于R _h ，则微调簇C _i 与C _j ，将簇C _i 中离簇C _j 最近的节点划归为C _j ，并且被划分出去的节点同时向C _j 发送请求加入消息和向C _i 发送撤销加入消息。

进一步，在所述步骤5）中，非簇头节点基于不同的TDMA时间片向各自的簇头节点发送数据，簇头节点在接收不同的数据后采取融合策略。

再进一步，在所述步骤5）中，节点会在空闲、发送、接收和睡眠四个状态下流转一旦非簇头节点被授予TDMA时间片，即被唤醒转入空闲状态等待发送数据到相应的簇头节点，而簇头节点也将被唤醒处于准备接收数据的状态。

更进一步，在所述步骤5）中，基站在获取各簇头向全网广播的消息包后及时分配给簇头节点在t时刻的唯一网络ID，其他非簇头节点与簇头节点通信时在数据包头包含此网络ID，簇头节点通过网络ID判定数据包是否发送给自身。

在所述步骤5）中，邻近簇之间的通信采用的 CSMA/CA MAC协议。

本发明的技术构思为：首先进行必要的节点距离侦测，计算存储转发有效的节点距离关系矩阵路由表。其次利用N阶近邻理论实现最优聚簇分析，通过统计计算获取的k _opt 计算最优的簇头占有比p _opt ，并结合节点剩余能量确定簇头选举概率；快速成簇后引入Hausdorff距离进一步优化已形成的簇结构，最终形成全网的稳定通信。

本发明的有益效果主要表现在：最优簇数计算准确，簇头分布更均匀，簇结构合理，节省了各节点能耗，实现了无线传感网的网络负载均衡。

附图说明

图1是无线传感网(WSNs)网络结构图。

图2是无线能耗模型图。

图3是通信执行周期图。

图4是层次型路由方法流程图。

图5是节点状态转换设计图。

具体实施方式

下面结合附图对本发明作进一步描述。

参照图1～图4，一种适用于无线传感网的层次型路由控制方法，包括如下步骤：

1）节点距离侦测存储阶段

无线传感节点在随机部署（如飞机飞行过程中投撒节点）后，即保持固定位置不变，各节点之间形成网络拓扑结构，基站负责统一的配置与通信调度。基站首先向各节点广播侦测请求包，节点在接收请求数据包后向基站传回自身的信息；基站收到回复信息后计算各自的距离并形成路由距离矩阵表，同时存储并转发相应的路由节点距离矩阵表给各自的节点，即节点只存储与自身相关的距离信息。

2）最优聚簇分析阶段

假定在指定区域内部署的N个无线传感节点组成的集合表示为S={x _i | x _i  =(a _i ,b _i )且x _i ∈R ², i = 1,2,…,N }，通过分析节点集的关系，可以找到满足一定条件的k个簇。令

为节点x _i 的一阶最短近邻距离，

为其二阶最短近邻距离，以此类推，则按照节点x _i 的最短近邻距离从小到大的排列有：

                    (1)

进一步，分别计算j阶最短近邻距离的平均值

和平方均值

：

                      (2)

                     (3)

传感节点的分散程度将直接决定最优簇数量，为分析其离散度，求解方差：

                (4)

假定y _i 为原传感节点x _i 的次序排列，且满足y ₁ ≤ y ₂ … ≤ y _N-1的关系，同时令

为其对应独立的且具有相同概率的概率密度函数，且必须满足。计算得到概率密度的K阶统计f _N-1,K(y)：

                     (5)

其中，f _N-1,K(y)为y _i 的累计概率分布函数，P(K(y) ≤ y)为随机试验K阶次序传感节点集中小于等于随机变量y的概率。相应的，上述过程可以看成在N-1次试验中大于等于K次成功的事件概率，由于y _i 的概率密度函数相同，且满足关系y ₁ ≤ y ₂ … ≤ y _N-1，则式(5)可转化成：

          (6)

显然，上式可进一步写成：

               (7)

对上式进行推导得：

          (8)

将式(8)的求和项展开，相同项约掉，便得到下式：

           (9)

由于项

，则K阶统计的概率密度f _N-1,K(y)为：

              (10)

其中的F(y)是对应的累积概率分布函数：

                             (11)

此时，式(2)、式(3)可等价转化成下式：

                         (12)

                           (13)

再将式(12)与式(13)代入式(4)得：

           (14)

同时，再将式(10)与式(11)分别代入可得：

  (15)

从上式可以发现，若想计算最终的

，必须给定相应的概率密度函数f(y)。由于一般部署的无线传感节点满足均匀分布的条件，因此f(y) = 1，上式便可简化成：

               (16)

进一步运用Beta函数及其性质

代入可得：

               (17)

利用求和公式与平方求和公式可将上式转化为：

                (18)

显然，式(18)可看成函数：

                 (19)

系数

、分别对应于：

 ,               (20)

       从上面的一系列推导发现，节点离散度

只与传感节点的个数N与节点m相关，将

做

线性化处理后得到：

                   (21)

设

，在统计传感节点曲率变化时，如果在某一节点上曲率变化满足

                 (22)

则认为在这一节点上发生了跃变，该节点也被称为跃变节点，记为m ^*。其中m为节点，阀值

是在(0,1)区间中的足够小的常数，且上式必须满足

的情况，否则将不成立。通过判定近邻节点曲率变化关系，可确定跃变的发生及获取跃变节点，即可以通过求解曲率变化确定簇的存在和簇数量大小。不失一般性，假定网络初始簇数为k=1，当计算某一节点曲率时满足式(22)，则意味着在该节点发生跃变，簇数在原来的基础上加1，即当前必定存在两个簇，使得在计算统计节点曲率时发生跃变。迭代遍历节点集，统计所有的跃变节点m ^*，就可得到最终的簇数量；而且随着跃变节点的增加，簇数k也必定为单调递增。迭代执行全部过程，即可得到最优的簇数量k _opt ，这也为下一步簇头选举打下基础。

3）簇头选举阶段

簇头的选取在层次型路由中起着关键作用，阀值T(n)决定当前传感节点成为簇头节点的概率，而簇头在所有节点中所占百分比p则决定着T(n)大小，其设定值也直接影响着整个网络的性能。通过前述对传感节点的最优聚簇分析可以得到当前最优的簇数量k _opt ，进而计算当前最优的簇头占用比p _opt ：

                              (23)

则T(n)可表示为：

                (24)

其中，r为执行的轮数，N为传感节点集的个数，G为本周期内未当选为簇头的节点集合。当r到达一个周期执行时间，未当选为簇头的节点集合G将被重置。

考虑传感节点剩余能量对簇头选取乃至整个网络生存周期的影响， 将节点剩余能量纳入簇头选取评价标准，则T(n)进一步转变为：

           (25)

其中E _residual 为节点当前评估的剩余能量，E _max 则为节点的最大能量，即为节点初始能量。

4）簇形成阶段

一旦簇头选举完成，各簇头节点采用基于载波监听多路访问/防止冲突（CSMA/CA）MAC协议及时向整个网络通告自己当选的消息，该消息包内容只包含簇头节点的ID信息，同时簇头节点等待其他非簇头节点请求加入。其他非簇头节点直接根据初始状态下基站测定并转发的近似距离存储表迅速选择加入最近的簇头，仍按照CSMA/CA协议向簇头发回加入请求消息，该消息包含当前的非簇头节点ID以及欲加入的请求指令。

快速成簇过程存在许多缺陷，即使前述已经实行了最优聚簇分析和最优簇头选取策略，也无法保证当前簇就是最优的。因此，为使网络划分的簇更趋均匀合理，引入Hausdorff距离来调整已经初步形成的簇，使得调整后的簇是最优的。

假设N个传感节点围绕k个簇头初步形成的簇集为Clusters={C ₁,C ₂,…,C _k}，节点v ₁

C _i 、v ₂

C _j ，且满足v ₁

v ₂、C _i

C_j，1 i

k,1 j

k，则簇C _i 到C _j 的Hausdorff距离为：

                  (26)

其中的

为簇C _i中的节点到C _j 中节点的范数，如L ₂范数。反过来C _j到C _i的Hausdorff距离不一定与C _i 到C _j 的Hausdorff距离相同，因此，最终两簇之间的Hausdorff距离可以用下式来表示：

 

             (27)

在确定簇间的Hausdorff距离后，判断两簇之间的非簇头节点是否满足Hausdorff距离条件，如果不满足则对簇进行调整。假设簇C _i 与C _j 的Hausdorff距离H(C _i ,C _j ) = R _h ，那么簇C _i 中每个非簇头节点到C _j 的节点必须满足d(v) < R _h （v代表簇C _i 中任一非簇头节点），反之亦然。当C _i 中存在一个或者多个节点大于或者等于R _h ，则微调簇C _i 与C _j ，将簇C _i 中离簇C _j 最近的节点划归为C _j ，并且被划分出去的节点同时向C _j 发送请求加入消息和向C _i 发送撤销加入消息。

5）数据传输阶段

数据传输阶段主要是将各传感节点收集的数据信息通过簇头进一步发送至基站。非簇头节点基于不同的TDMA时间片向各自的簇头节点发送数据，簇头节点在接收不同的数据后采取融合策略以降低冗余度或还原关键原始数据。通常节点会在空闲、发送、接收和睡眠四个状态下流转，鉴于节省能耗的设计原则，其大部分的时间处于睡眠状态，一旦非簇头节点被授予TDMA时间片，即被唤醒转入空闲状态等待发送数据到相应的簇头节点，而簇头节点也将被唤醒处于准备接收数据的状态。由于节点主要的能耗都在数据通信上，数据融合及睡眠状态时消耗极小部分能量，因此能保证尽可能地延长网络生存时间。

此外，对于邻近簇间的通信干扰，按照特定的网络编码方式，即基站在获取各簇头向全网广播的消息包后及时分配给簇头节点在t时刻的唯一网络ID，其他非簇头节点与簇头节点通信时必须在数据包头包含此网络ID，簇头节点通过网络ID判定数据包是否发送给自身，这样可以过滤掉大量邻簇的数据包。同时，采用的 CSMA/CA MAC协议同样可以避免大量通信冲突干扰的发生，使得全网能够保持稳定的数据通信传输。

如图1，无线传感节点随机部署在一定区域内，该区域可以是规则的几何图形区域，也可以是不规则图形区域，节点与节点之间、节点与基站之间形成多层网络拓扑结构。无线传感节点通过选举机制选取出簇头集，节点采集的各类信息由这些簇头集传输给基站，这种层次型的网络拓扑有利于网络能耗的均衡。

图2是本发明所采用的无线能耗模型，当节点v ₁向距离d外的v ₂发送L字节数据包时的能耗为：

               (28)

E _Tx (L,d)包含发射电路及功率放大损耗两部分，其中E _elec 为节点发射每位(bit)数据消耗的能量，参数

_fs 和

_mp 依赖于具体的传输放大模型，d _threshold 为阀值，当传感节点v ₁与v ₂的距离d

d _threshold 时，则使用Friss自由空间模型，反之当d > d _threshold 时则采用多路径衰减模型。特别的，当d = d _threshold ，则有：

                         (29)

同时，v ₂接收v ₁发送过来的数据包消耗的能量为：

                         (30)

       如图3通信执行周期所示，本发明采用轮机制，主要分为5个步骤，分别为距离侦测存储、最优聚簇分析、簇头选举、簇形成和数据通信。距离侦测阶段主要是节点初始化部署后，网络拓扑自然形成，需要对网络整体进行侦测。通过基站向节点广播询问消息以及节点的回复确立距离关系，并在基站计算存储该距离矩阵路由表，同时基站也会向节点转发相应的距离关系表。距离侦测存储只在节点初始化部署后执行，而且只执行一次。

接着，最优聚簇阶段实行N阶近邻理论分析，通过对节点的离散度及曲率变化分析，统计确定现有网络拓扑条件下可分簇数量。离散度的计算见式(4)及其推导，曲率分析计算主要见式(22)，当计算某一节点曲率时满足式(22)，则意味着在该节点发生了曲率上的跃变，那么当前就必定存在两个簇，迭代运行至计算得出完全的簇数量k _opt 。该过程只在节点初始化部署及网络拓扑动态变化的情况发生。

然后簇头选举阶段将基于最优簇数k _opt 计算最佳簇头占有比率p _opt ，同时考虑节点的剩余能量，选举当前最优的簇头集。该过程主要是计算节点当选簇头的概率T(n)计算，每一轮都将根据实际的情况选举新的簇头，结合节点的剩余能量考虑，能量越多的节点将被优先考虑选为簇头。而簇头占有比率p _opt 和簇的地域划分也决定着簇头的选举。

在簇头选举结束后，簇头节点采用基于载波监听多路访问/防止冲突（CSMA/CA）MAC协议及时向整个网络通告自己当选的消息，同时等待其他非簇头节点请求加入。其他非簇头节点直接根据初始状态下基站测定并转发的近似距离存储表迅速选择加入最近的簇头，仍按照CSMA/CA协议向簇头发回加入请求消息。在快速成簇，本发明又引入了Hausdorff距离，通过对簇间的Hausdorff距离判定，调整已经形成的簇结构，使之更趋合理和均衡。

最后进入数据稳定通信阶段，该阶段主要是在已经形成的WSNs通信网络上进行通信，非簇头节点基于不同的TDMA时间片向各自的簇头节点发送数据，簇头节点在接收不同的数据后采取融合策略以降低冗余度或还原关键原始数据。经过融合的数据信息通过簇头最后传送至基站，基站再对这些数据进行分类整理或二次融合，以达到进一步降低数据量的目的。

详细的流程如图4所示。为了满足能耗节省的设计要求，如图5所示，通常节点被设计成空闲、发送、接收和睡眠四个状态。节点大部分时间将处于睡眠状态，一旦非簇头节点被授予TDMA时间片，即被唤醒转入空闲状态等待发送数据到相应的簇头节点，而簇头节点也将被唤醒处于准备接收数据的状态。这样的设计模式能够尽可能地降低能力消耗，节省有限的电力资源。此外，为了充分保证整个通信网络的工作状态和效率，满足关系T _analysis + T _select + T _form << T _data，其中T _analysis表示最优聚簇分析时间，T _select表示簇头选举时间，以及T _form表示簇形成时间。

Claims (5)

1.一种适用于无线传感网的层次型路由控制方法，包括以下步骤：

1)节点距离侦测存储：基站首先向各节点广播侦测请求包，节点在接收请求数据包后向基站传回自身的信息；基站收到回复信息后计算节点与节点之间的距离并形成路由距离矩阵表，同时存储并转发相应的路由节点距离矩阵表给各自的节点，即节点只存储与自身相关的距离信息；其特征在于：所述层次型路由控制方法还包括以下步骤：

2)最优聚簇分析：设定在指定区域内部署的N个无线传感节点组成的集合表示为S＝{x_i|x_i=(a_i,b_i)且x_i∈R²,i=1,2,…,N}，令

为节点x_i的一阶最短近邻距离，

为节点x_i的二阶最短近邻距离，以此类推，

为节点x_i的N-1阶最短近邻距离，则按照节点x_i的最短近邻距离从小到大的排列有：

$d_1^i<d_2^i<\&CenterDot;\&CenterDot;\&CenterDot;<d_{N-1}^i\&ForAll;i=\mathrm{1,2},...,N-1---\left(1\right)$

进一步，分别计算j阶最短近邻距离的平均值

和平方均值

${\stackrel{\&OverBar;}{d}}_j=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{d}_j^i/N\&ForAll;j=\mathrm{1,2},...N-1---\left(2\right)$

${\stackrel{\&OverBar;}{d}}_j^2=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{\left(d_j^i\right)}^2/N\&ForAll;j=\mathrm{1,2},...N-1---\left(3\right)$

求解方差

${\mathrm{\&sigma;}}_{N,m}^2=\frac{1}{m}\underset{j=1}{\overset{m}{\mathrm{\&Sigma;}}}\left({{\stackrel{\&OverBar;}{d}}_j^2-\left({\stackrel{\&OverBar;}{d}}_j\right)}^2\right)\&ForAll;j,m=\mathrm{1,2}...,N-1---\left(4\right)$

假定y_i为原无线传感节点x_i的次序排列，且满足y₁≤y₂…≤y_N-1的关系，同时令

为其对应独立的且具有相同概率的概率密度函数

且必须满足

设定概率密度函数f(y)=1，运用Beta函数的特性将式（4）推导转换为式（19）：

${\mathrm{\&sigma;}}_{N,m}^2=\mathrm{\&alpha;}\left(N\right){(m+1)}^2+\mathrm{\&beta;}\left(N\right)(m+1)---\left(19\right)$

系数α(N)、β(N)分别对应于：

$\mathrm{\&alpha;}\left(N\right)=-\frac{1}{3N^2(N+1)},$ $\mathrm{\&beta;}\left(N\right)=-\frac{3N+1}{6N^2(N+1)}---\left(20\right)$

设m≥2，在无线传感节点曲率变化时，如果在某一节点上曲率变化满足

$|\frac{{\mathrm{\&sigma;}}_{N,m+1}^2}{m+2}+\frac{{\mathrm{\&sigma;}}_{N,m-1}^2}{m}-2\frac{{\mathrm{\&sigma;}}_{N,m}^2}{m+1}|<\mathrm{\&epsiv;}\left|\frac{{\mathrm{\&sigma;}}_{N,m}^2}{m+1}\right|---\left(22\right)$

则认为在该无线传感节点上发生了跃变，该无线传感节点也被称为跃变节点，记为m^*，其中的阀值ε是在(0,1)区间中的足够小的常数；

假定网络初始簇数为k＝1，当计算某一节点曲率时满足式(22)，簇数在原来的基础上加1，迭代遍历节点集，统计所有的跃变节点m^*，得到最优的簇数量k_opt；3)簇头选举：阀值T(n)决定当前无线传感节点成为簇头节点的概率，而簇头在所有节点中所占百分比p则决定着T(n)大小，通过对无线传感节点的最优聚簇分析得到当前最优的簇数量k_opt，进而计算当前最优的簇头占用比p_opt：

$p_{\mathrm{opt}}=\frac{k_{\mathrm{opt}}}{N}---\left(23\right)$

则T(n)表示为：

$T\left(n\right)=\left\{\begin{array}{cc}\frac{k_{\mathrm{opt}}}{N-k_{\mathrm{opt}}\&times;\left[r\mathrm{mod}(N/k_{\mathrm{opt}})\right]}& \mathrm{ifn}\&Element;G\\ 0& \mathrm{otherwise}\end{array}\right.---\left(24\right)$

其中，r为执行的轮数，N为无线传感节点集的个数，G为本周期内未当选为簇头的节点集合，当r到达一个周期执行时间，未当选为簇头的节点集合G将被重置；

4)簇形成：各簇头节点及时向整个网络通告自己当选的消息，该消息包内容只包含簇头节点的ID信息，同时簇头节点等待其他非簇头节点请求加入，其他非簇头节点直接根据初始状态下基站测定并转发的近似距离存储表迅速选择加入最近的簇头，并向簇头发回加入请求消息，该消息包含当前的非簇头节点ID以及欲加入的请求指令；

5)数据传输：将各无线传感节点收集的数据信息通过簇头进一步发送至基站。

2.如权利要求1所述的一种适用于无线传感网的层次型路由控制方法，其特征在于：在所述步骤3)中，将节点剩余能量纳入簇头选取评价标准，则公式（24）转变为式（25）：

$T\left(n\right)=\left\{\begin{array}{cc}\frac{k_{\mathrm{opt}}}{N-k_{\mathrm{opt}}\&times;\left[r\mathrm{mod}(N/k_{\mathrm{opt}})\right]}\&times;\frac{E_{\mathrm{residual}}}{E_{\max }}& \mathrm{ifn}\&Element;G\\ 0& \mathrm{otherwise}\end{array}\right.---\left(25\right)$

其中E_residual为节点当前评估的剩余能量，E_max则为节点的最大能量，即为节点初始能量。

3.如权利要求1或2所述的一种适用于无线传感网的层次型路由控制方法，其特征在于：在所述步骤5）中，节点会在空闲、发送、接收和睡眠四个状态下流转，一旦非簇头节点被授予TDMA时间片，即被唤醒转入空闲状态等待发送数据到相应的簇头节点，而簇头节点也将被唤醒处于准备接收数据的状态。

4.如权利要求1或2所述的一种适用于无线传感网的层次型路由控制方法，其特征在于：在所述步骤5）中，基站在获取各簇头向全网广播的消息包后及时分配给簇头节点在t时刻的唯一网络ID，其他非簇头节点与簇头节点通信时在数据包头包含此网络ID，簇头节点通过网络ID判定数据包是否发送给自身。

5.如权利要求4所述的一种适用于无线传感网的层次型路由控制方法，其特征在于：在所述步骤5）中，邻近簇之间的通信采用的CSMA/CAMAC协议。

Images