一种移动传感网的Sink节点移动路径规划方法技术

本发明专利技术涉及一种移动传感网的Sink节点移动路径规划方法，包括如下步骤：初始化各种参数；将Sink节点的数据收集范围分解成一个半径为r的内圆和厚度均为r的n‑1个圆环；监测区域分成NP个子区域，并取得初始化后的粒子群；计算每一个网格的单位节点能耗；更新个体极值粒子和全局极值粒子；粒子m2分别与个体极值粒子和全局极值粒子进行交叉和变异操作；当该粒子的移动路径长度超过允许的最大移动路程LP时，执行路径修正操作；m1＝m1+1，如果m1大于M1，则结束算法，输出最优移动路径和最大网络生存时间，否则m2＝1，重复计算单位节点能耗。本发明专利技术可寻找到Sink节点的最优移动路径，尽量让传感节点的能量消耗相同，从而平衡传感节点的能耗，提高网络生存时间。

【技术实现步骤摘要】
一种移动传感网的Sink节点移动路径规划方法
本专利技术涉及通信
，尤其涉及无线传感网
，具体是指一种移动传感网的Sink节点移动路径规划方法。
技术介绍
无线传感网(wirelesssensornetworks,WSNs)由空间分布的自主传感节点和Sink节点组成，可协同监测温度、声音、振动、压力、运动等物理和环境参数，在战场监控、战场损伤评估、工业过程监控、机器运行状况监控、家庭自动化、交通监控等军事、工业和民用领域上具有较大的应用价值和市场潜力。在无线传感网中，能量效率是无线传感网的最重要问题之一。由于传感节点的电池能量有限而且其电池更换需要花费较多的时间和精力，因此无线传感网需要在无人干预下自行运行足够长的时间。但是在静态Sink节点的无线传感网中，由于多跳路由具有热点且数据流量集中，靠近Sink节点的传感节点消耗较多的能量且失效较早。这个问题被称为无线传感网的热点问题。热点问题可恶化无线传感网的网络连接和生存时间，甚至使无线传感网不能正常工作。为克服热点问题，很多学者研究Sink节点的移动，让其周围的热点发生变化且出现在网络的各个区域，从而实现统一的能量消耗和延长网络生存时间。如郜帅等人建立优先级高且移动距离不小于阈值的虚拟点集合，采用TSP算法求解能遍历该集合内所有虚拟点的最短路径。郭剑等人将监测区域分成多个大小相同的圆盘。计算每一个圆盘内Sink节点的停留位置，采用量子遗传算法获得最短移动路径。汪林云等人利用分簇技术将传感节点分成通信半径相同的簇，选择剩余能量相对充足的传感节点作为簇头。根据簇头的位置，利用TSP算法构建一条尽可能短的移动路径。SalarianH.等人选择权值较高的传感节点作为RP点(rendezvouspoint)，建立点集合，寻找能访问所有RP点且不超过最大数据传输时延的路径。LeeK.等人考虑传感节点网格分布和邻居传感节点的剩余能量变化，提出一种Sink节点移动的分布式启发算法，寻找其停留位置。但是这些算法都假设所有传感节点都能获知自身位置坐标。为保证定位精度，需要安装GPS或北斗定位模块，这增加了节点能耗和系统硬件成本。而且算法所获得的移动路径是否具有已经达到最大提高网络生存时间的效果，仍需要进一步研究，而且很多移动路径的确定需要传感节点的位置、能量等信息。
技术实现思路
为寻找传感节点均匀分布时Sink节点的最优移动路径和最大网络生存时间，本专利技术考虑传感节点均匀分布的情况，提供了一种移动传感网的Sink节点移动路径规划方法，该方法可寻找到Sink节点的最优移动路径，从而平衡网络能耗，提高网络生存时间。为了实现上述目的，本专利技术具有如下构成：该移动传感网的Sink节点移动路径规划方法，所述的方法包括以下步骤：(1)初始化当前迭代次数m1＝1、当前粒子m2＝1、预设圆环厚度值r、Sink节点的最大通信距离d、预设子区域个数NP、最大迭代次数M1、预设粒子个数M2和允许的最大移动路程LP；(2)将Sink节点的数据收集范围分解成一个半径为r的内圆和厚度均为r的n-1个圆环，其中n＝d/r；(3)将监测区域分成NP个子区域，并对每个子区域依次进行编号，并重复执行M2次粒子的初始操作，得到初始粒子群；(4)依次选择各个粒子m2中Sink节点的停留位置和计算每一个子区域的节点能耗，获得每一个子区域的单位节点能耗；(5)根据每一个子区域的单位节点能耗计算网络生存时间，更新个体极值粒子和全局极值粒子，m2＝m2+1，并判断m2是否大于M2，如果是，则获得各个粒子的适应度值、个体极值和全局极值，令m2＝1，然后继续步骤(6)；否则继续步骤(4)；(6)粒子m2分别与所述个体极值粒子和全局极值粒子进行交叉和变异操作；(7)当该粒子的移动路径长度超过允许的最大移动路程LP时，执行路径修正操作，m2＝m2+1，并判断m2是否大于M2，如果是，则继续步骤(8)，否则继续步骤(4)；(8)判断当前迭代次数m1是否小于预设最大迭代次数M1，如果是，则m2＝1，m1＝m1+1，则继续步骤(4)，否则结束算法，并输出最优移动路径和最大网络生存时间。可选地，所述步骤(2)中，将监测区域分成大小相同的NP个矩形网格，各个矩形网格即为一个子区域，并按照从左到右、从上到下的原则对各个所述矩形网格进行编码。可选地，所述步骤(3)中，重复执行M2次粒子的初始操作，得到初始粒子群，包括如下步骤：在各个所述子区域内随机生成一个Sink节点的停留位置，构成具有NP个停留位置的粒子。可选地，所述步骤(4)中，计算所述单位节点能耗，包括如下步骤：(4-1)Sink节点的当前停留位置为(xs,ys),该子区域中心i的位置为(xi,yi)，监测区域的长和宽为(Lx,Ly)，经过Sink节点和该子区域中心i的直线方程为：Sink节点到该子区域中心i的有向射线与监测区域的边界交点为(xc,yc)：当xs＝xi且ys＝yi时，则两点重合，令yi＝yi+0.1，则xc＝xi，yc＝Ly；当xs＝xi且ys＞yi时，则xc＝xi，yc＝0；当xs＝xi且ys<yi时，则xc＝xi，yc＝Ly；当xs＞xi且ys＝yi时，则xc＝0，yc＝yi；当xs<xi且ys＝yi时，则xc＝Lx，yc＝yi；当xs≠xi且ys≠yi时，如果xs<xi，yi≤ys(xi-Lx)/(xs-Lx)或xs＞xi，yi≤ysxi/xs时，则xc＝-yi(xs-xi)/(ys-yi)+xi，yc＝0；如果xs<xi，ys(xi-Lx)/(xs-Lx)≤yi≤(ys-Ly)(xi-Lx)/(xs-Lx)+Ly时，则xc＝Lx，yc＝(ys-yi)(Lx-xi)/(xs-xi)+yi；如果xs<xi,yi＞(ys-Ly)(xi-Lx)/(xs-Lx)+L或xs＞xi,yi＞(ys-Ly)xi/xs+Ly时，则xc＝(Ly-yi)(xs-xi)/(ys-yi)+xi，yc＝Ly；如果xs＞xi，ys(xi)/xs≤yi≤(ys-Ly)xi/xs+Ly时，则xc＝0，yc＝-xi(ys-yi)/(xs-xi)+yi；(4-2)根据Sink节点到该子区域中心i的有向射线与监测区域的边界交点和Sink节点的多跳数据收集方式，计算经过该子区域中心i的多跳通信收集的最大跳数Ni,max为：其中，表示大于或等于x的最小整数，计算该子区域中心i所在位置到Sink节点的跳数Ni为：(4-3)该子区域中心i的单位节点能耗包括接收和发送该位置到边界交点区域内传感节点的数据能耗和该位置上发送自身数据的能耗，根据如下公式计算该子区域中心i的单位节点能耗其中，Si表示Sink节点的感知数据速率，Eelec表示Sink节点接收或发送单位bit数据的电路能耗，ε表示Sink节点发送单位bit数据的放大器能耗；(4-4)根据如下公式计算该子区域的单位节点能耗Ei：其中，P表示该粒子的移动路径，为一个m×2的矩阵，每一行存放Sink节点停留位置的横坐标和纵坐标，m表示移动路径P中Sink节点停留位置个数，p表示移动路径P中Sink节点的停留位置。可选地，所述步骤(5)中，根据如下公式，计算所述网络生存时间：其中，T表示网络生存时间，即T表示初始能量本文档来自技高网...

【技术保护点】
一种移动传感网的Sink节点移动路径规划方法，其特征在于，所述的方法包括以下步骤：(1)初始化当前迭代次数m1＝1、当前粒子m2＝1、预设圆环厚度值r、Sink节点的最大通信距离d、预设子区域个数NP、最大迭代次数M1、预设粒子个数M2和允许的最大移动路程LP；(2)将Sink节点的数据收集范围分解成一个半径为r的内圆和厚度均为r的n‑1个圆环，其中n＝d/r；(3)将监测区域分成NP个子区域，并对每个子区域依次进行编号，并重复执行M2次粒子的初始操作，得到初始粒子群；(4)依次选择各个粒子m2中Sink节点的停留位置和计算每一个子区域的节点能耗，获得每一个子区域的单位节点能耗；(5)根据每一个子区域的单位节点能耗计算网络生存时间，更新个体极值粒子和全局极值粒子，m2＝m2+1，并判断m2是否大于M2，如果是，则获得各个粒子的适应度值、个体极值和全局极值，令m2＝1，然后继续步骤(6)；否则继续步骤(4)；(6)粒子m2分别与所述个体极值粒子和全局极值粒子进行交叉和变异操作；(7)当该粒子的移动路径长度超过允许的最大移动路程LP时，执行路径修正操作，m2＝m2+1，并判断m2是否大于M2，如果是，则继续步骤(8)，否则继续步骤(4)；(8)判断当前迭代次数m1是否小于预设最大迭代次数M1，如果是，则m2＝1，m1＝m1+1，则继续步骤(4)，否则结束算法，并输出最优移动路径和最大网络生存时间。...

【技术特征摘要】
1.一种移动传感网的Sink节点移动路径规划方法，其特征在于，所述的方法包括以下步骤：(1)初始化当前迭代次数m1＝1、当前粒子m2＝1、预设圆环厚度值r、Sink节点的最大通信距离d、预设子区域个数NP、最大迭代次数M1、预设粒子个数M2和允许的最大移动路程LP；(2)将Sink节点的数据收集范围分解成一个半径为r的内圆和厚度均为r的n-1个圆环，其中n＝d/r；(3)将监测区域分成NP个子区域，并对每个子区域依次进行编号，并重复执行M2次粒子的初始操作，得到初始粒子群；(4)依次选择各个粒子m2中Sink节点的停留位置和计算每一个子区域的节点能耗，获得每一个子区域的单位节点能耗；(5)根据每一个子区域的单位节点能耗计算网络生存时间，更新个体极值粒子和全局极值粒子，m2＝m2+1，并判断m2是否大于M2，如果是，则获得各个粒子的适应度值、个体极值和全局极值，令m2＝1，然后继续步骤(6)；否则继续步骤(4)；(6)粒子m2分别与所述个体极值粒子和全局极值粒子进行交叉和变异操作；(7)当该粒子的移动路径长度超过允许的最大移动路程LP时，执行路径修正操作，m2＝m2+1，并判断m2是否大于M2，如果是，则继续步骤(8)，否则继续步骤(4)；(8)判断当前迭代次数m1是否小于预设最大迭代次数M1，如果是，则m2＝1，m1＝m1+1，则继续步骤(4)，否则结束算法，并输出最优移动路径和最大网络生存时间。2.根据权利要求1所述的移动传感网的Sink节点移动路径规划方法，其特征在于，所述步骤(2)中，将监测区域分成大小相同的NP个矩形网格，各个矩形网格即为一个子区域，并按照从左到右、从上到下的原则对各个所述矩形网格进行编码。3.根据权利要求1所述的移动传感网的Sink节点移动路径规划方法，其特征在于，所述步骤(3)中，重复执行M2次粒子的初始操作，得到初始粒子群，包括如下步骤：在各个所述子区域内随机生成一个Sink节点的停留位置，构成具有NP个停留位置的粒子。4.根据权利要求1所述的移动传感网的Sink节点移动路径规划方法，其特征在于，所述步骤(4)中，计算所述单位节点能耗，包括如下步骤：(4-1)Sink节点的当前停留位置为(xs,ys),该子区域中心i的位置为(xi,yi)，监测区域的长和宽为(Lx,Ly)，经过Sink节点和该子区域中心i的直线方程为：Sink节点到该子区域中心i的有向射线与监测区域的边界交点为(xc,yc)：当xs＝xi且ys＝yi时，则两点重合，令yi＝yi+0.1，则xc＝xi，yc＝Ly；当xs＝xi且ys＞yi时，则xc＝xi，yc＝0；当xs＝xi且ys<yi时，则xc＝xi，yc＝Ly；当xs＞xi且ys＝yi时，则xc＝0，yc＝yi；当xs<xi且ys＝yi时，则xc＝Lx，yc＝yi；当xs≠xi且ys≠yi时，如果xs<xi，yi≤ys(xi-Lx)/(xs-Lx)或xs＞xi，yi≤ysxi/xs时，则xc＝-yi(xs-xi)/(ys-yi)+xi，yc＝0；如果xs<xi，ys(xi-Lx)/(xs-Lx)≤yi≤(ys-Ly)(xi-Lx)/(xs-Lx)+Ly时，则xc＝Lx，yc＝(ys-yi)(Lx-xi)/(xs-xi)+yi；如果xs<xi,yi＞(ys-Ly)(xi-Lx)/(xs-Lx)+L或xs＞xi,yi＞(ys-Ly)xi/xs+Ly时，则xc＝(Ly-yi)(xs-xi)/(ys-yi)+xi，yc＝Ly；如果xs＞xi，ys(xi)/xs≤yi≤(ys-Ly)xi/xs+Ly时，则xc＝0，yc＝-xi(ys-yi)/(xs-xi)+yi；(4-2)根据Sink节点到该子区域中心i的有向射线与监测区域的边界交点和Sink...

【专利技术属性】
技术研发人员：陈友荣，万锦昊，陆思一，刘半藤，
申请(专利权)人：浙江树人学院，
类型：发明
国别省市：浙江,33