发明名称 基于时空相关性分簇的无线传感器网络压缩数据收集方法
摘要 一种基于时空相关性分簇的无线传感器网络压缩数据收集方法,它包含以下步骤:一、初始化阶段;二、时空相关性度量;三、网络成簇:汇聚节点根据节点间的时空相关性计算值,得到分簇信息后向全网广播信息,簇内节点再动态的选出簇头节点;四、压缩数据收集;五、数据重构。通过以上步骤,基于时空相关性分簇的无线传感器网络压缩数据收集方法,在相同重构误差的情况下,减少了数据传输次数,均衡了网络能耗,避免了某个节点因能耗问题过早死亡,延长了网络生存时间。
申请公布号 CN105682171B 申请公布日期 2017.03.22
申请号 CN201610132885.X 申请日期 2016.03.09
申请人 北京航空航天大学 发明人 万江文;王东豪;聂志朋;姚雷
分类号 H04W40/10(2009.01)I;H04W40/32(2009.01)I;H04W84/20(2009.01)I;H04L12/715(2013.01)I;H04L12/721(2013.01)I 主分类号 H04W40/10(2009.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种基于时空相关性分簇的无线传感器网络压缩数据收集方法,其特征在于:它包含以下步骤:步骤一,初始化阶段:在这个阶段,每个传感器采集数据,并把采集到的原始数据以多跳的方式发送给汇聚节点;汇聚节点保存每个传感器节点在各个时间槽的采样数据;步骤二,时空相关性度量:当汇聚节点收到数据后,进行数据的时空相关性度量阶段;汇聚节点根据步骤一中收集到的数据,按照基于时空相关性分簇的无线传感器网络压缩数据收集方法提出的区域数据相关性度量方法,计算任意两个传感器节点间的数据时空相关性;令x<sup>(k)</sup>,x<sub>i</sub>分别为网络第k次采样的数据向量和第i个传感器在T个时间槽的采样向量,N为网络中节点的数目;定义Gmax_dst为邻域内两个节点间的最大距离,即最大跳数,网内任意两个节点间的距离小于Gmax_dst时,归为一个邻域N<sub>c</sub>;只考虑在邻域N<sub>c</sub>中节点对数据的相关性计算,以符合实际情况;所述基于时空相关性分簇的无线传感器网络压缩数据收集方法提出的区域数据相关性度量方法,具体如下,步骤2.1:计算任意邻域内两个节点的时间相关性,如果邻域内的任意节点对(i,j)存在时间相关性,ρ<sub>m</sub>(i,j)&gt;ρ<sub>m_min</sub>,ρ<sub>m</sub>(i,j)表示邻域内两个节点的时间相关性;其中ρ<sub>m_min</sub>为自定义的时间相关性阈值,<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>&rho;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&rho;</mi><mi>m</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>&rho;</mi><mi>m</mi></msub><mo>(</mo><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>T</mi></mrow></mfrac><mrow><mo>(</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mrow><mo>(</mo><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mi>E</mi><mo>&lsqb;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&rsqb;</mo><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mi>E</mi><mo>&lsqb;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&rsqb;</mo><mo>)</mo></mrow></mrow><msubsup><mi>&sigma;</mi><msub><mi>x</mi><mi>i</mi></msub><mn>2</mn></msubsup></mfrac><mo>+</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mrow><mo>(</mo><msubsup><mi>x</mi><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mi>E</mi><mo>&lsqb;</mo><msub><mi>x</mi><mi>j</mi></msub><mo>&rsqb;</mo><mo>)</mo></mrow><mrow><mo>(</mo><msubsup><mi>x</mi><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>m</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mi>E</mi><mo>&lsqb;</mo><msub><mi>x</mi><mi>j</mi></msub><mo>&rsqb;</mo><mo>)</mo></mrow></mrow><msubsup><mi>&sigma;</mi><msub><mi>x</mi><mi>j</mi></msub><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>,</mo><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>&Element;</mo><msub><mi>N</mi><mi>c</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001205690700000011.GIF" wi="1534" he="327" /></maths>式中,N<sub>c</sub>表示上述中定义的邻域,<img file="FDA0001205690700000012.GIF" wi="66" he="62" />和<img file="FDA0001205690700000013.GIF" wi="107" he="63" />表示第i个节点分别在第k次和第k+m次的采样值,x<sub>i</sub>表示第i个传感器在T个时间槽的采样向量,E[x<sub>i</sub>]和<img file="FDA0001205690700000014.GIF" wi="59" he="63" />分别表示第i个传感器在T个时间槽的采样均值和采样方差,同理可知<img file="FDA0001205690700000015.GIF" wi="266" he="64" />E[x<sub>j</sub>]和<img file="FDA0001205690700000016.GIF" wi="54" he="70" />的表示含义;步骤2.2:计算任意邻域内两个节点的空间相关性,如果邻域内的任意节点对(i,j)存在空间相关性,则需满足ρ<sub>s</sub>(i,j)&gt;ρ<sub>s_min</sub>,ρ<sub>s</sub>(i,j)表示邻域内两个节点的空间相关性;其中ρ<sub>s_min</sub>为自定义的空间相关性阈值,<maths num="0002"><math><![CDATA[<mrow><msub><mi>&rho;</mi><mi>S</mi></msub><mrow><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>T</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>T</mi></munderover><mfrac><mrow><mrow><mo>(</mo><mrow><msubsup><mi>x</mi><mi>i</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mi>E</mi><mrow><mo>&lsqb;</mo><msub><mi>x</mi><mi>i</mi></msub><mo>&rsqb;</mo></mrow></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><msubsup><mi>x</mi><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>-</mo><mi>E</mi><mrow><mo>&lsqb;</mo><msub><mi>x</mi><mi>j</mi></msub><mo>&rsqb;</mo></mrow></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>&sigma;</mi><msub><mi>x</mi><mi>i</mi></msub></msub><msub><mi>&sigma;</mi><msub><mi>x</mi><mi>j</mi></msub></msub></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo></mrow><mo>&Element;</mo><msub><mi>N</mi><mi>c</mi></msub></mrow>]]></math><img file="FDA0001205690700000021.GIF" wi="997" he="150" /></maths>式中,N<sub>c</sub>表示上述中定义的邻域,<img file="FDA0001205690700000022.GIF" wi="68" he="60" />表示第i个节点在第k次的采样值,x<sub>i</sub>表示第i个传感器在T个时间槽的采样向量,E[x<sub>i</sub>]和<img file="FDA0001205690700000023.GIF" wi="57" he="45" />分别表示第i个传感器在T个时间槽的采样均值和采样标准差,同理可知<img file="FDA0001205690700000024.GIF" wi="98" he="62" />E[x<sub>j</sub>]和<img file="FDA0001205690700000025.GIF" wi="56" he="54" />的表示含义;步骤2.3:计算任意邻域内两个节点的时空相关性,如果邻域内的任意节点对(i,j)存在时空相关性,则需满足ρ<sub>ms</sub>(i,j)&gt;ρ<sub>msmin</sub>,ρ<sub>ms</sub>(i,j)表示邻域内两个节点的时空相关性;其中ρ<sub>msmin</sub>为自定义的空间相关性阈值,ρ<sub>ms</sub>(i,j)=ρ<sub>m</sub>(i,j)*ρ<sub>s</sub>(i,j)式中,ρ<sub>m</sub>(i,j)和ρ<sub>s</sub>(i,j)的含义分别在步骤2.1和步骤2.2中给出了;步骤三,网络成簇:汇聚节点根据任意两个传感器节点间的数据时空相关性,得到分簇信息后向全网广播信息,簇内节点再动态的选出簇头节点;具体如下,步骤3.1:根据设定的自定义的空间相关性阈值ρ<sub>msmin</sub>和步骤二中得到的时空相关性数据,若邻域内节点数据相关性大于阈值门限,就认为这两个节点之间存在一条连边(u,v),构造图模型G,图的顶点就是网络中传感器节点,图的边定义为上述连边(u,v);步骤3.2:根据图模型G中节点度值和节点剩余能量确定分簇和簇头节点;根据步骤3.1得到的图模型G中的连边关系,若区域内所有节点都相连,则分为一个簇,簇内以概率P<sub>ch</sub>=E<sub>residual</sub>/E<sub>max</sub>选出簇头节点,其中E<sub>residual</sub>表示节点的剩余能量,E<sub>max</sub>表示节点的初始能量,P<sub>ch</sub>表示节点当选为簇头的概率;步骤3.3:汇聚节点确认簇内节点数目和簇头节点的位置,并由压缩感知理论确定簇内压缩感知采样数目不少于C kln(n/k),其中C为采样因子,k为簇内信号的稀疏度,n为簇内节点数目;C kln(n/k)示以自然数e为底的指数;步骤四,压缩数据收集:由于汇聚节点保存有每个节点的历史数据,当数据的相关性变化小于设定的阈值时,分簇信息不需要改变,直接进行压缩感知数据收集;所述压缩数据收集方法具体如下:步骤4.1:簇内节点和簇头节点使用相同的伪随机数产生器,生成测量系数Φ<sup>i</sup>,所以测量系数不需要在网络中传输;簇内节点采用压缩感知数据收集的方法,将原始数据与测量系数的线性运算结果传递给簇头节点;簇内节点按照步骤4.2进行数据传递;簇内节点采用的压缩感知数据收集方法具体如步骤4.2所述;步骤4.2:测量矩阵Φ<sup>i</sup>为对应第i个簇的子矩阵,CH<sup>i</sup>表示第i个簇的簇头节点,负责收集簇内的所有数据x<sup>i</sup>,则单个簇内数据收集表示为y<sup>i</sup>,y<sup>i</sup>=Φ<sup>i</sup>x<sup>i</sup>;步骤4.3:簇头节点CH<sup>i</sup>采用最小距离平方生成树算法,将测量值y<sup>i</sup>传输到汇聚节点;步骤五,数据重构:汇聚节点接收来自簇头的数据,构成测量值向量<img file="FDA0001205690700000031.GIF" wi="116" he="414" />然后汇聚节点产生相同的随机矩阵Φ,并对原始收集数据进行重构,假设网络中分为5个簇,则重构式满足:<img file="FDA0001205690700000032.GIF" wi="611" he="407" />式中,测量系数Φ<sup>i</sup>为对应第i个簇的子矩阵,x<sup>i</sup>表示第i个簇内的采样向量,y<sup>i</sup>表示对应第i个簇的压缩感知测量值。
地址 100191 北京市海淀区学院路37号