发明名称 无线传感网监测数据的空间相关性建模与定量分析方法及系统
摘要 本发明公开了一种无线传感网监测数据的空间相关性建模与定量分析方法及系统,该方法包括:从小规模无线传感网中获取N个节点采集到的实际数据序列W=(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>)<sup>T</sup>,对实际数据序列W进行稀疏变换,并按节点访问的顺利将最远的非零系数对应的距离确定为参数h;根据公式γ(i),得到h+1个关于参数α<sub>i</sub>,β,σ<sub>i</sub>的非线性方程组;求解该非线性方程组,得到h+1个参数α<sub>i</sub>,β;将得到的2h+1个参数α<sub>i</sub>,β,σ<sub>i</sub>代入模型X(x,y)中,生成与实际数据序列W的空间相关性保持一致的大规模数据集W<sup>*</sup>。本发明不仅可从少量的无线传感网实际监测数据中快速地抽取参数,高精度地生成与之空间相关性匹配的数据集;而且还可以灵活地控制模型的参数,人工合成具有任意空间相关性强弱的大规模数据集。
申请公布号 CN105574290A 申请公布日期 2016.05.11
申请号 CN201610030621.3 申请日期 2016.01.15
申请人 中国地质大学(武汉) 发明人 陈分雄;胡凯;赵天明;凌承昆;唐曜曜;王典洪
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 湖北武汉永嘉专利代理有限公司 42102 代理人 唐万荣
主权项 一种无线传感网监测数据的空间相关性建模与定量分析方法,其特征在于,包括以下步骤:S1、从小规模无线传感网中获取N个节点采集到的实际数据序列W=(x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>)<sup>T</sup>,对实际数据序列W进行稀疏变换,得到k个非零系数,并按节点访问的顺利将最远的非零系数对应的距离确定为参数h,h为整数;S2、根据公式γ(i),得到h个方程,将其与方程<img file="FDA0000907646430000011.GIF" wi="291" he="142" />联立,得到h+1个关于参数α<sub>i</sub>,β,σ<sub>i</sub>的非线性方程组,其中:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>&gamma;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>&beta;</mi></mrow><mrow><mn>2</mn><mi>&beta;</mi></mrow></mfrac><mrow><mo>(</mo><msub><mi>&alpha;</mi><mn>1</mn></msub><msup><msub><mi>&sigma;</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msub><mi>&alpha;</mi><mi>2</mi></msub><msup><msub><mi>&sigma;</mi><mn>2</mn></msub><mn>2</mn></msup><mo>+</mo><mo>...</mo><mo>+</mo><msub><mi>&alpha;</mi><mi>h</mi></msub><msup><msub><mi>&sigma;</mi><mi>h</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></munderover><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>h</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></munderover><mfrac><msub><mi>&alpha;</mi><mi>k</mi></msub><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mi>&gamma;</mi><mrow><mo>(</mo><msub><mi>d</mi><mrow><msub><mi>i</mi><mi>j</mi></msub><msub><mi>k</mi><mi>l</mi></msub></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000907646430000012.GIF" wi="1357" he="143" /></maths>S3、采用最小二乘法计算h个参数σ<sub>i</sub>,并代入非线性方程组;S4、求解该非线性方程组,得到h+1个参数α<sub>i</sub>,β;S5、将得到的2h+1个参数α<sub>i</sub>,β,σ<sub>i</sub>代入模型X(x,y)中,生成与实际数据序列W的空间相关性保持一致的大规模数据集W<sup>*</sup>,模型X(x,y)的公式为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>X</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msubsup><mi>X</mi><mn>1</mn><mn>1</mn></msubsup><mo>+</mo><msub><mi>Z</mi><mn>1</mn></msub></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mfrac><msub><mi>&alpha;</mi><mn>1</mn></msub><mrow><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>X</mi><mn>1</mn><mrow><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msub><mi>Z</mi><mn>1</mn></msub></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mfrac><msub><mi>&alpha;</mi><mn>1</mn></msub><mrow><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mrow><msubsup><mi>X</mi><mn>2</mn><mn>1</mn></msubsup><mo>+</mo><msub><mi>Z</mi><mn>2</mn></msub></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mfrac><msub><mi>&alpha;</mi><mn>2</mn></msub><mrow><mi>N</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>X</mi><mn>2</mn><mrow><mi>N</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msub><mi>Z</mi><mn>2</mn></msub></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mfrac><msub><mi>&alpha;</mi><mn>2</mn></msub><mrow><mi>N</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>X</mi><mi>h</mi><mn>1</mn></msubsup><mo>+</mo><msub><mi>Z</mi><mi>h</mi></msub></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mfrac><msub><mi>&alpha;</mi><mi>h</mi></msub><mrow><mi>N</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mrow></mrow></mtd><mtd><mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>X</mi><mi>h</mi><mrow><mi>N</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msub><mi>Z</mi><mi>h</mi></msub></mrow></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mfrac><msub><mi>&alpha;</mi><mi>h</mi></msub><mrow><mi>N</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mi>Y</mi></mtd><mtd><mrow><mi>w</mi><mo>.</mo><mi>p</mi><mo>.</mo></mrow></mtd><mtd><mi>&beta;</mi></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000907646430000013.GIF" wi="963" he="1295" /></maths>其中,X<sub>i</sub>,Z<sub>i</sub>,Y共同决定了X(x,y),其中Z<sub>i</sub>(i=1,2,…,h)、Y和X为两两相互独立的随机变量,Z<sub>i</sub>表示相邻节点采集数据值间的差异,Y表示与节点X(x,y)距离超出h倍d后,X(x,y)取值仅与Y有关而与相邻节点X<sub>i</sub>无关;w.p.表示为中文的意思“关于”,可用“关于”两个字代替;<img file="FDA0000907646430000021.GIF" wi="126" he="126" />表示随机变量X(x,y)取值为<img file="FDA0000907646430000022.GIF" wi="167" he="78" />的平均概率,β表示随机变量X(x,y)取值为Y的概率。
地址 430074 湖北省武汉市洪山区鲁磨路388号