本发明涉及一种基于压缩感知的无线传感器网络动态数据融合树方法，属于无线传感器网络技术领域。该方法包括以下步骤：将压缩感知与动态路由结合；以稀疏随机投影机制为基础；根据投影切换时，节点加入和离开的情况，依据最小代价动态的生成数据融合树；传感器节点将采集到的数据沿数据融合树传输。本发明的技术方案采用稀疏投影机制，可以有效的降低网络数据传输量；将压缩感知与路由协议相结合，可以有效的减少网络能耗；投影切换时，考虑了节点的离开和加入，动态生成数据融合树，方法的计算复杂度很小。因此，本发明的方法在有效地降低了网络能耗的同时，减少了网络节点的计算负担，能够解决网络生存周期短和降低网络计算负载的问题。
1.一种基于压缩感知的无线传感器网络动态数据融合树方法，其特征在于，该方法包括以下步骤：
步骤1：初始化：设置节点集合Si←{}，路径树集合Ti←{}；
步骤2：将网络中的节点一维化Node{N(1),N(2),...,N(n)}，其中N为节点标识，n为节点数目；
步骤3：根据以下公式，获得一个m*n的稀疏随机投影矩阵：
其中Φi,j为第i行j列的投影系数，参数s确定投影矩阵的稀疏度k，p为概率；当s＝n/(lg n)时，O(k lg n)次投影的重构精度等效于最大K系数算法，k为投影矩阵的稀疏度，n为网络的节点数；
矢量φi中非零值所对应的节点即为参与第i次投影的点，将其加入源节点集Si＝{Si(1),Si(2),...,Si(logn)}，1≤i≤m；
步骤4：在第一次投影即i＝1时，依据Dijkstra算法，求出各源节点形成最短路径树Ti＝{Ti(1),Ti(2),...,Ti(logn)}；
步骤5：从第二次投影即i>＝2起，参与上次投影而本次投影不参与的节点，进行删除处理；没有参与上次投影而参与本次投影的节点，进行加入处理；
步骤6：节点采集数据后，按照生成的数据融合树进行传输；
从原生成树中删除一个多播结点s的方法是：如果该节点是多播树中的叶结点，则直接删除s及与s相连的其他节点；若该节点是中间结点，删除s后，将s的所有子孙结点中的端结点组成一个集合Md，将删除了集合Md的子树记为Tt，采用如下算法将Md中的结点加入生成树：
1)初始化：计算Md中每一个结点到子树Tt的最短路径，记为T_DIS，并令集合Q为空；
2)在Md不为空的条件下，选择Md中到Tt距离最短或在原生成树中除接应节点外的节点到接应结点距离最短的端结点min_Md，将min_Md和min_Md到生成树最小路径上所有的Steiner结点和min_Md在原生成树中所有的子孙结点一起加入生成树，从Md中删除所有己加入生成树的端点；若Md不为空转下一步，否则程序结束；
3)若目前的接应结点与在原生成树中的接应结点不同，则在生成树中每加入一个接应结点u，判断其是否是原生成树中的结点；若不是，则将该结点加入Q；每在Q中加入一个结点，考察新加入结点到Md中每一结点v的距离，若该距离小于结点v到子树Tt的距离，则将该距离作为结点v到子树Tt的距离；转2)。
2.如权利要求1所述的方法，其特征在于，在本方法中，在进行投影切换时，把每次投影过程看作一次路由选择，以相同概率选择部分节点参与投影，若每次参与投影的节点能形成最短路径树，就通过该最短路径树将参与投影节点的数据加权和传输到汇聚节点。
3.如权利要求1所述的方法，其特征在于，在步骤5中，投影切换时，加入节点的方法如下：
1)初始化：令k＝1，从源点s开始，将该原点s作为T1；此时Tk＝T1,集合M1＝Mk＝{s}并将所有要加入的结点组成集合N，将N中每一个结点到s的距离作为结点到Tk的距离记为y_DIS[i]；并令集合Q为空，令Mk中所有结点到集合Q的距离为无穷大；
2)静态计算最小代价多播生成树时，各端点加入的顺序：在Mk和N不为空的条件下，依次比较dist[k]与N到Tk的最短路径y_DIS_min；若dist[k]不大于y_DIS_min，则将相应的端点和该端点到生成树最小路径上所有的Steiner结点一起加入生成树(即该端点的路径结点)，从Mk中删除该端点；比较N中每一个结点到每一个新加入结点的距离,若该距离小于y_DIS_min，则将该距离作为y_DIS_min；每加入一个结点，k的值加1，重复本步骤；若dist[k]大于y_DIS_min或Mk为空，则将y_DIS_min对应的端点y和y到z的路径上所有的Steiner结点一起加入生成树，其中z为y的接应结点；若Mk不为空转下一步，否则程序结束；
3)生成树每加入一个结点u，判断新加入是否是原生成树中的结点，若不是，则将该结点加入Q；
4)每在Q中加入一个结点，考察新加入结点到Mk中每一结点v的距离，若该距离小于结点v到集合Q的距离，则将该距离作为结点v到集合Q的距离，并在从中选出到集合Q距离最近的结点；
5)设Mk中到集合Q距离最近的结点到集合的Q距离为D_t，比较y_DIS_min、dist[k]、D_t三者的大小，若dist[k]最小，则将相应的端点vk和该端点到生成树最小路径上所有的Steiner结点一起加入生成树，其中生成树每加入一个结点u，考察u到集合N中每一个元素i的距离，若小于y_DIS[i]，则令y_DIS[i]为u到元素i的距离，令u为元素i的接应结点；并从Mk中删除该端点，若Mk为空，程序结束；若vk就是Mk中到集合Q距离最近的结点，则从Mk中重新选择一个到集合Q距离最近的结点；令k＝k+1直到vk∈Mk，重复步骤5)；若D_t最小，则将Mk中到集合Q距离最近的端点和该端点到生成树最小路径上所有的Steiner结点一起加入生成树从Mk中删除该端点，若Mk为空，程序结束，转3)；若y_DIS_min最小，将y_DIS_min对应的端点y和y到z的路径上所有的Steiner结点一起加入生成树和集合Q，转4)，其中z为y的接应结点。
一种基于压缩感知的无线传感器网络动态数据融合树方法
技术领域
[0001] 本发明属于无线传感器网络技术领域，涉及一种基于压缩感知的无线传感器网络动态数据融合树方法。
背景技术
[0002] 无线传感器网络WSNs是由大量普通节点和一个汇聚节点组成的多跳自组织网络。具有自组织、部署迅捷、高容错性和强隐蔽性等技术优势。经过几十年的发展，WSNs在科学研究和工程领域得到了更广泛应用。但是由于成本、体积等问题，传感器节点的能量和计算能力都是有限的。我们可以采用低复杂度、而且有效的方法来对无线传感器网络需要传输的数据进行融合，来减少传输数据，降低能耗，提高网络寿命。压缩感知是近些年来兴起的信号处理理论，在该理论框架下，网络采集部分信息就可以重构出完整信息。所以如何将压缩感知理论应用到无线传感器网的数据融合中，是一个十分具有研究意义的课题。
[0003] 近来，基于压缩感知理论的数据采集机制得到了较多关注。但是，传统的压缩感知方法也带来了密集投影的问题。由于在一般情况下M＜＜n/2(m和n分别为观测次数和节点个数)，采用压缩感知的方式可以大幅节省能量；但当投影次数m增大时，压缩感知的能耗可能会接近传统方式甚至更高。因此我们需要一种更高效的投影机制，在保证信号能以高概率重构的情况下进一步节省能量。
发明内容
[0004] 有鉴于此，本发明的目的在于提供一种基于压缩感知的无线传感器网络动态数据融合树方法，能够解决一般算法计算复杂度高和网络能耗高的问题。
[0005] 为达到上述目的，本发明提供如下技术方案：
[0006] 一种基于压缩感知的无线传感器网络动态数据融合树方法，该方法包括以下步骤：
[0007] 步骤1：初始化：设置节点集合Si←{}，路径树集合Ti←{}；
[0008] 步骤2：将网络中的节点一维化Node{N(1),N(2),...,N(n)}；
[0009] 步骤3：根据以下公式，获得稀疏随机投影矩阵：
[0010]
[0011] 其中Φi,j为第i行j列的投影系数，参数s确定投影矩阵的稀疏度，p为概率；当s＝n/(lg n)时，O(k lg n)次投影的重构精度等效于最大K系数算法；
[0012] 矢量φi(1≤i≤m)中非零值所对应的节点即为参与第i次投影的点，将其加入源节点集Si＝{Si(1),Si(2),...,Si(logn)}；
[0013] 步骤4：在第一次投影(i＝1)时，依据Dijkstra算法，求出各源节点形成最短路径树Ti＝{Ti(1),Ti(2),...,Ti(logn)}；
[0014] 步骤5：从第二次投影(i>＝2)起，参与上次投影而本次投影不参与的节点，进行删除处理；没有参与上次投影而参与本次投影的节点，进行加入处理；
[0015] 步骤6：节点采集数据后，按照生成的数据融合树进行传输。
[0016] 进一步，在本方法中，在进行投影切换时，把每次投影过程看作一次路由选择，以相同概率选择部分节点参与投影，若每次参与投影的节点能形成最短路径树，就通过该最短路径树将参与投影节点的数据加权和传输到汇聚节点。
[0017] 进一步，在步骤5中，投影切换时，删除节点的方法包括：
[0018] 从原生成树中删除一个多播结点s的方法是：如果该节点是多播树中的叶结点，则直接删除s及与s相连的其他节点；若该节点是中间结点，删除s后，将s的所有子孙结点中的端结点组成一个集合Md，将删除了集合Md的子树记为Tt，采用如下算法将Md中的结点加入生成树：
[0019] 1)初始化：计算Md中每一个结点到子树Tt的最短路径，记为T_DIS，并令集合Q为空；
[0020] 2)在Md不为空的条件下，选择Md中到Tt距离最短或在原生成树中到接应结点距离最短(且接应结点在不中)的端结点min_Md，将min_Md和min_Md到生成树最小路径上所有的Steiner结点和min_Md在原生成树中所有的子孙结点一起加入生成树，从Md中删除所有己加入生成树的端点；若Md不为空转下一步，否则程序结束；
[0021] 3)若x目前的接应结点与在原生成树中的接应结点不同，则在生成树中每加入一个接应结点u，判断其是否是原生成树中的结点；若不是，则将该结点加入Q；每在Q中加入一个结点，考察新加入结点到Md中每一结点v的距离，若该距离小于结点v到子树Tt的距离，则将该距离作为结点v到子树Tt的距离；转2)。
[0022] 进一步，在步骤5中，投影切换时，加入节点的方法如下：
[0023] 1)初始化：令k＝1，从源点s开始，将单结点s作为T1；此时Tk＝T1,V1＝Vk＝{s}并将所有要加入的结点组成集合N，将N中每一个结点到s的距离作为结点到Tk的距离记为y_DIS[i]；并令集合Q为空，令Mk中所有结点到集合Q的距离为无穷大；
[0024] 2)静态计算最小代价多播生成树时，各端点加入的顺序：在Mk和N不为空的条件下，依次比较dist[k]与N到Tk的最短路径y_DIS_min；若dist[k]不大于y_DIS_min，则将相应的端点和该端点到生成树最小路径上所有的Steiner结点一起加入生成树(即该端点的路径结点)，从Mk中删除该端点；比较N中每一个结点到每一个新加入结点的距离,若该距离小于y_DIS_min，则将该距离作为y_DIS_min；每加入一个结点，k的值加1，重复本步骤；若dist[k]大于y_DIS_min或Mk为空，则将y_DIS_min对应的端点y和y到z(z为y的接应结点)的路径上所有的Steiner结点一起加入生成树；若Mk不为空转下一步,否则程序结束；
[0025] 3)生成树每加入一个结点u，判断新加入是否是原生成树中的结点，若不是，则将该结点加入Q；
[0026] 4)每在Q中加入一个结点，考察新加入结点到Mk中每一结点v的距离，若该距离小于结点v到集合Q的距离，则将该距离作为结点v到集合Q的距离，并在从中选出到集合Q距离最近的结点；
[0027] 5)设Mk中到集合Q距离最近的结点到集合的Q距离为D_t,比较y_DIS_min、dist[k]、D_t三者的大小，若dist[k]最小，则将相应的端点vk和该端点到生成树最小路径上所有的Steiner结点一起加入生成树(生成树每加入一个结点u,考察u到集合N中每一个元素i的距离,若小于y_DIS[i]，则令y_DIS[i]为u到元素i的距离，令u为元素i的接应结点)；从Mk中删除该端点(若Mk为空,程序结束)；若vk就是Mk中到集合Q距离最近的结点，则从Mk中重新选择一个到集合Q距离最近的结点；令k＝k+1直到vk∈Mk，重复步骤5)；若D_t最小，则将Mk中到集合Q距离最近的端点和该端点到生成树最小路径上所有的Steiner结点一起加入生成树从Mk中删除该端点(若Mk为空,程序结束)，转3)；若y_DIS_min最小，将y_DIS_min对应的端点y和y到z(z为y的接应结点)的路径上所有的Steiner结点一起加入生成树和集合Q，转4)。
[0028] 本发明的有益效果在于：
[0029] 1、采用稀疏投影机制，可以有效的降低网络数据传输量；
[0030] 2、将压缩感知与路由协议相结合，可以有效的减少网络能耗；
[0031] 3、投影切换时，考虑了节点的离开和加入，动态生成数据融合树，方法的计算复杂度很小。因此，本发明实施例在有效地降低了网络能耗的同时，减少了网络节点的计算负担。
附图说明
[0032] 为了使本发明的目的、技术方案和有益效果更加清楚，本发明提供如下附图进行说明：
[0033] 图1为本发明的稀疏投影示意图；
[0034] 图2为本发明方法的流程图；
具体实施方式
[0035] 下面将结合附图，对本发明的优选实施例进行详细的描述。
[0036] 图2为本发明所述方法的流程图，如图所示，本发明提供的基于压缩感知的无线传感器网络动态数据融合树方法，包括：
[0037] 步骤1：初始化：设置节点集合Si←{}，路径树集合Ti←{}；
[0038] 步骤2：将网络中的节点一维化Node{N(1),N(2),...,N(n)}；
[0039] 步骤3：根据以下公式，获得稀疏随机投影矩阵：
[0040]
[0041] 其中Φi,j为第i行j列的投影系数，参数s确定投影矩阵的稀疏度，p为概率；当s＝n/(lg n)时，O(k lg n)次投影的重构精度等效于最大K系数算法；
[0042] 矢量φi(1≤i≤m)中非零值所对应的节点即为参与第i次投影的点，将其加入源节点集Si＝{Si(1),Si(2),...,Si(logn)}；
[0043] 步骤4：在第一次投影(i＝1)时，依据Dijkstra算法，求出各源节点形成最短路径树Ti＝{Ti(1),Ti(2),...,Ti(logn)}；
[0044] 步骤5：从第二次投影(i>＝2)起，参与上次投影而本次投影不参与的节点，进行删除处理；没有参与上次投影而参与本次投影的节点，进行加入处理；
[0045] 步骤6：节点采集数据后，按照生成的数据融合树进行传输。
[0046] 在本方法中，在进行投影切换时，把每次投影过程看作一次路由选择，以相同概率选择部分节点参与投影，若每次参与投影的节点能形成最短路径树，就通过该最短路径树将参与投影节点的数据加权和传输到汇聚节点。图1为本发明的稀疏投影示意图。
[0047] 在步骤5中，投影切换时，删除节点的方法包括：
[0048] 从原生成树中删除一个多播结点s的方法是：如果该节点是多播树中的叶结点，则直接删除s及与s相连的其他节点；若该节点是中间结点，删除s后，将s的所有子孙结点中的端结点组成一个集合Md，将删除了集合Md的子树记为Tt，采用如下算法(即最小化代价原则)将Md中的结点加入生成树：
[0049] 1)初始化：计算Md中每一个结点到子树Tt的最短路径，记为T_DIS，并令集合Q为空；
[0050] 2)在Md不为空的条件下，选择Md中到Tt距离最短或在原生成树中到接应结点距离最短(且接应结点在不中)的端结点min_Md，将min_Md和min_Md到生成树最小路径上所有的Steiner结点和min_Md在原生成树中所有的子孙结点一起加入生成树，从Md中删除所有己加入生成树的端点；若Md不为空转下一步，否则程序结束；
[0051] 3)若x目前的接应结点与在原生成树中的接应结点不同，则在生成树中每加入一个接应结点u，判断其是否是原生成树中的结点；若不是，则将该结点加入Q；每在Q中加入一个结点，考察新加入结点到Md中每一结点v的距离，若该距离小于结点v到子树Tt的距离，则将该距离作为结点v到子树Tt的距离；转2)。
[0052] 在步骤5中，投影切换时，加入节点的方法如下：
[0053] 1)初始化：令k＝1，从源点s开始，将单结点s作为T1；此时Tk＝T1,V1＝Vk＝{s}并将所有要加入的结点组成集合N，将N中每一个结点到s的距离作为结点到Tk的距离记为y_DIS[i]；并令集合Q为空，令Mk中所有结点到集合Q的距离为无穷大；
[0054] 2)静态计算最小代价多播生成树时，各端点加入的顺序：在Mk和N不为空的条件下，依次比较dist[k]与N到Tk的最短路径y_DIS_min；若dist[k]不大于y_DIS_min，则将相应的端点和该端点到生成树最小路径上所有的Steiner结点一起加入生成树(即该端点的路径结点)，从Mk中删除该端点；比较N中每一个结点到每一个新加入结点的距离,若该距离小于y_DIS_min，则将该距离作为y_DIS_min；每加入一个结点，k的值加1，重复本步骤；若dist[k]大于y_DIS_min或Mk为空，则将y_DIS_min对应的端点y和y到z(z为y的接应结点)的路径上所有的Steiner结点一起加入生成树；若Mk不为空转下一步,否则程序结束；
[0055] 3)生成树每加入一个结点u，判断新加入是否是原生成树中的结点，若不是，则将该结点加入Q；
[0056] 4)每在Q中加入一个结点，考察新加入结点到Mk中每一结点v的距离，若该距离小于结点v到集合Q的距离，则将该距离作为结点v到集合Q的距离，并在从中选出到集合Q距离最近的结点；
[0057] 5)设Mk中到集合Q距离最近的结点到集合的Q距离为D_t,比较y_DIS_min、dist[k]、D_t三者的大小，若dist[k]最小，则将相应的端点vk和该端点到生成树最小路径上所有的Steiner结点一起加入生成树(生成树每加入一个结点u,考察u到集合N中每一个元素i的距离,若小于y_DIS[i]，则令y_DIS[i]为u到元素i的距离，令u为元素i的接应结点)；从Mk中删除该端点(若Mk为空,程序结束)；若vk就是Mk中到集合Q距离最近的结点，则从Mk中重新选择一个到集合Q距离最近的结点；令k＝k+1直到vk∈Mk，重复步骤5)；若D_t最小，则将Mk中到集合Q距离最近的端点和该端点到生成树最小路径上所有的Steiner结点一起加入生成树从Mk中删除该端点(若Mk为空,程序结束)，转3)；若y_DIS_min最小，将y_DIS_min对应的端点y和y到z(z为y的接应结点)的路径上所有的Steiner结点一起加入生成树和集合Q，转4)。
[0058] 最后说明的是，以上优选实施例仅用以说明本发明的技术方案而非限制，尽管通过上述优选实施例已经对本发明进行了详细的描述，但本领域技术人员应当理解，可以在形式上和细节上对其作出各种各样的改变，而不偏离本发明权利要求书所限定的范围。
