发明名称 在无线传感器网络中基于压缩字典学的数据收集方法
摘要 本发明提出了在无线传感器网络中基于压缩字典学的数据收集方法,包含步骤:一、字典的构造:把要学的字典D表示为一个基字典不同原子的线性组合;二、网络数据传输:节点感知周围数据信息,采用链式的压缩感知收集方法进行数据收集;Y;三、字典更新:把字典的自相干程度作为约束项进行优化,并交替进行稀疏编码和字典更新,得到最优化的稀疏表示字典;四、数据重构:得到稀疏表示字典后可以精确的重构原始信号,实现数据的收集。通过以上步骤,基于压缩字典学的数据收集方法,在同样的恢复精度要求下,减少了数据传输次数,延长了网络寿命,同时还能对环境噪声有一定的抑制作用,有广泛的适用性和自适应性。
申请公布号 CN105811993B 申请公布日期 2017.02.15
申请号 CN201610139692.7 申请日期 2016.03.11
申请人 北京航空航天大学 发明人 万江文;王东豪;张强;陈俊英
分类号 H03M7/30(2006.01)I 主分类号 H03M7/30(2006.01)I
代理机构 代理人
主权项 一种在无线传感器网络中基于压缩字典学习的数据收集方法,其特征在于:它包含以下步骤:步骤一,字典的构造:实际信号包含固定的结构信息,能进行稀疏表示;把要学习的字典D表示为一个基字典不同原子的线性组合,基字典是正交稀疏表示矩阵DCT;DCT矩阵能够对信号进行稀疏表示,结合信号的结构性,采用基字典原子的线性组合来学习出新的字典,表示如下:D=ΨΘ,其中<img file="FDA0001185176610000011.GIF" wi="169" he="55" />为基字典,<img file="FDA0001185176610000012.GIF" wi="190" he="54" />ρ≥k,约束为一个稀疏矩阵,k维的向量阵列来描述不同结构信息;这里Θ有不同的形状,当Θ为k×k的单位矩阵时,矩阵D就是一个普通的字典;当Θ=[I<sub>k×k</sub>|Σ<sub>k×σ</sub>],矩阵D为一个长方形的冗余字典;根据汇聚节点是否已知信号的结构信息,得到基字典的不同表达形式;如果汇聚节点已知信号的结构信息,直接进入步骤1.3,否则要首先执行步骤1.1,1.2;步骤1.1:初始化,历史数据收集阶段,网络传输原始信号X=[x<sup>1</sup>,x<sup>2</sup>,...,x<sup>L</sup>]到汇聚节点,其中<img file="FDA0001185176610000013.GIF" wi="149" he="53" />表示一个N维传感数据向量,L为传输数据的数量;步骤1.2:基字典的学习,用矩阵<img file="FDA0001185176610000014.GIF" wi="174" he="54" />来表示这些训练数据X=[x<sup>1</sup>,x<sup>2</sup>,...,x<sup>L</sup>];字典学习问题的表示形式为<img file="FDA0001185176610000015.GIF" wi="620" he="94" />其中||·||<sub>F</sub>为矩阵Frobenius范数,<img file="FDA0001185176610000016.GIF" wi="159" he="54" />是稀疏表示矩阵,S表示信号的稀疏度,使用K‑SVD算法求解上述字典学习问题,得到基字典Ψ;步骤1.3:由上述步骤学习到的基字典Ψ,对字典D加上稀疏矩阵约束,表示为D=ΨΘ,矩阵Θ是稀疏的,每一列都是一个稀疏向量;步骤二,网络数据传输:确定无线传感器网络中汇聚节点的位置,节点感知周围数据信息,按照压缩感知收集方法进行数据收集,最后传递给汇聚节点,获得感知数据Y;步骤2.1:利用步骤一中的字典D对原始信号进行稀疏化,得到稀疏表示矩阵X<sup>*</sup>,D=ΨΘ,X=DX<sup>*</sup>=ΨΘX<sup>*</sup>;步骤2.2:汇聚节点生成M×N维高斯测量矩阵Φ,其中<img file="FDA0001185176610000017.GIF" wi="182" he="46" />中的元素都独立服从均值为0,方差为<img file="FDA0001185176610000018.GIF" wi="59" he="103" />的高斯分布,即<img file="FDA0001185176610000019.GIF" wi="293" he="102" />步骤2.3:通过测量矩阵Φ对稀疏表示矩阵X<sup>*</sup>进行压缩采样,每个节点i在传输时,该节点把所有接收到的值相加,然后再加上x<sub>i</sub>Φ<sub>M,i</sub>作为新的值传出给该节点的下一跳节点,最后得到压缩感知数据Y;步骤三,字典学习:把压缩感知数据Y直接用于字典的学习,对信号的稀疏表示矩阵实时更新,则得到以下约束优化问题,<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>arg</mi><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>D</mi></munder><mo>{</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><mi>&Phi;</mi><mi>D</mi><mo>&CenterDot;</mo><msup><mi>X</mi><mo>*</mo></msup><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>}</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mo>|</mo><mo>|</mo><msubsup><mi>x</mi><mi>i</mi><mo>*</mo></msubsup><mo>|</mo><msub><mo>|</mo><mn>0</mn></msub><mo>&le;</mo><mi>S</mi></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0001185176610000021.GIF" wi="853" he="118" /></maths>式中,||·||<sub>F</sub>表示矩阵Frobenius范数,<img file="FDA0001185176610000022.GIF" wi="117" he="77" />表示向量<img file="FDA0001185176610000023.GIF" wi="41" he="60" />的l<sub>0</sub>范数;上述约束优化问题是一个欠定系统方程,没有唯一解;为了得到唯一解,对上述约束优化问题引入新的约束条件,结合约束条件Θ的稀疏性和字典D的自相干性约束,得到最终的优化问题,<maths num="0002"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>&Theta;</mi><mo>,</mo><msup><mi>X</mi><mo>*</mo></msup></mrow></munder><mo>{</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><msup><mi>&Phi;&Psi;&Theta;X</mi><mo>*</mo></msup><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><msup><mrow><mo>(</mo><mi>&Psi;</mi><mi>&Theta;</mi><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><mi>&Psi;</mi><mi>&Theta;</mi><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>&lambda;</mi><mi>A</mi></msub><mo>|</mo><mo>|</mo><msup><mi>X</mi><mo>*</mo></msup><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>+</mo><msub><mi>&lambda;</mi><mi>&Theta;</mi></msub><mo>|</mo><mo>|</mo><mi>&Theta;</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>}</mo></mrow>]]></math><img file="FDA0001185176610000024.GIF" wi="1110" he="117" /></maths>式中,||·||<sub>F</sub>表示矩阵Frobenius范数,||·||<sub>1</sub>表示向量的l<sub>1</sub>范数,提出使用迭代法来解决上述最终的优化问题,每次迭代过程包含两个步骤:步骤3.1:稀疏逼近:在稀疏逼近阶段,固定字典D,求解压缩感知数据Y在字典D上的稀疏表示矩阵X<sup>*</sup>,即<maths num="0003"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><msup><mi>X</mi><mo>*</mo></msup></munder><mo>{</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><msup><mi>&Phi;&Psi;&Theta;X</mi><mo>*</mo></msup><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>&lambda;</mi><mi>A</mi></msub><mo>|</mo><mo>|</mo><msup><mi>X</mi><mo>*</mo></msup><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>}</mo></mrow>]]></math><img file="FDA0001185176610000025.GIF" wi="610" he="111" /></maths>式中,||·||<sub>F</sub>表示矩阵Frobenius范数,||·||<sub>1</sub>表示向量的l<sub>1</sub>范数,使用正交匹配追踪算法(OMP)求解上述稀疏表示矩阵X<sup>*</sup>;步骤3.2:字典更新:在字典更新阶段,保持矩阵X<sup>*</sup>的值不变,得到下面的优化问题<maths num="0004"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>&Theta;</mi></munder><mo>{</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><msup><mi>&Phi;&Psi;&Theta;X</mi><mo>*</mo></msup><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>+</mo><mo>|</mo><mo>|</mo><msup><mrow><mo>(</mo><mi>&Psi;</mi><mi>&Theta;</mi><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><mi>&Psi;</mi><mi>&Theta;</mi><mo>)</mo></mrow><mo>-</mo><mi>I</mi><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>&lambda;</mi><mi>&Theta;</mi></msub><mo>|</mo><mo>|</mo><mi>&Theta;</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0001185176610000026.GIF" wi="966" he="119" /></maths>式中,||·||<sub>F</sub>表示矩阵Frobenius范数,||·||<sub>1</sub>表示向量的l<sub>1</sub>范数,使用增广拉格朗日乘子法(ALM)来解决上述字典更新的优化问题,令<maths num="0005"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>Y</mi><mo>-</mo><msup><mi>&Phi;&Psi;&Theta;X</mi><mo>*</mo></msup><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>&lambda;</mi><mi>&Theta;</mi></msub><mo>|</mo><mo>|</mo><mi>&Theta;</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>,</mo><mi>h</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mi>&Psi;</mi><mi>&Theta;</mi><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><mi>&Psi;</mi><mi>&Theta;</mi><mo>)</mo></mrow><mo>-</mo><mi>I</mi></mrow>]]></math><img file="FDA0001185176610000027.GIF" wi="1222" he="119" /></maths>式中,||·||<sub>F</sub>表示矩阵Frobenius范数,||·||<sub>1</sub>表示向量的l<sub>1</sub>范数,定义拉格朗日增广项为<maths num="0006"><math><![CDATA[<mrow><mi>L</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>,</mo><mi>Y</mi><mo>,</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>+</mo><mo>&lt;</mo><mi>Y</mi><mo>,</mo><mi>h</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>&gt;</mo><mo>+</mo><mfrac><mi>&mu;</mi><mn>2</mn></mfrac><mo>|</mo><mo>|</mo><mi>h</mi><mrow><mo>(</mo><mi>&Theta;</mi><mo>)</mo></mrow><mo>|</mo><msubsup><mo>|</mo><mi>F</mi><mn>2</mn></msubsup></mrow>]]></math><img file="FDA0001185176610000031.GIF" wi="844" he="111" /></maths>求解上式得到稀疏矩阵Θ,进一步得到更新字典D=ΨΘ;步骤四,数据重构:由步骤三得到了字典D,求解<img file="FDA0001185176610000032.GIF" wi="803" he="111" />得到<img file="FDA0001185176610000033.GIF" wi="189" he="53" />式中,||·||<sub>F</sub>表示矩阵Frobenius范数,||·||<sub>1</sub>表示向量的l<sub>1</sub>范数。
地址 100191 北京市海淀区学院路37号