发明名称 基于能量有效的多任务贝叶斯压缩感知宽带频谱检测方法
摘要 本发明公开了一种基于能量有效的多任务贝叶斯压缩感知宽带频谱检测方法。本发明方法针对认知无线网络中主用户信号在空频域的稀疏性,基于贝叶斯压缩感知的信号重构通过层次化贝叶斯分析分级先验模型获得稀疏信号估计。利用多认知用户感知信号的时空相关性实现在多用户多任务传输条件下的稀疏信号重构与宽带压缩频谱检测。认知基站通过基于期望最大化算法和相关向量机模型进行多任务贝叶斯压缩感知参数估计。在满足一定检测性能和贝叶斯压缩感知重构均方误差的条件下,实现节点感知能耗最小化。本发明方法有效保障了多节点多任务宽带频谱检测的节点能量有效性。
申请公布号 CN104703216A 申请公布日期 2015.06.10
申请号 CN201510058061.8 申请日期 2015.02.04
申请人 杭州电子科技大学 发明人 许晓荣;王赞;姚英彪;章坚武
分类号 H04W24/06(2009.01)I;H04W24/10(2009.01)I;H04B17/382(2015.01)I 主分类号 H04W24/06(2009.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 基于能量有效的多任务贝叶斯压缩感知宽带频谱检测方法,其特征在于,该方法的具体步骤是:步骤(1)认知用户本地频谱感知:根据主用户PU在授权频谱上接入行为的不同,主用户PU的信号在空频域具有稀疏性,认知用户SU的感知参数随着主用户PU接入的变化而变化;认知用户SU通过低速采样获取观测数据进行本地频谱感知,并通过报告信道向认知基站CBS报告本地感知信息,认知基站利用贝叶斯压缩感知方法重构主用户PU宽带信道上的功率谱,实现在主用户PU宽带信道上的频谱检测;其具体方法是:主用户发射机PUT与主用户接收机PUR利用授权频谱进行通信,K个认知用户SU对主用户PU频谱占用情况进行本地感知,并选择最佳认知用户SU,利用最佳认知用户SU的报告信道向认知基站CBS汇报本地感知信息,认知基站CBS基于多任务贝叶斯压缩感知BCS进行信息融合与感知信号稀疏重构,然后根据信道能量累积进行在主用户PU宽带信道上的频谱检测;在t时刻,第k个认知用户SU进行本地频谱感知信号<img file="FDA0000667308740000011.GIF" wi="180" he="93" />为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>x</mi><msub><mi>CR</mi><mi>k</mi></msub></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>h</mi><mi>k</mi></msub><msub><mi>x</mi><mi>PUT</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><msub><mi>H</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>n</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><msub><mi>H</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000667308740000012.GIF" wi="769" he="165" /></maths>其中,h<sub>k</sub>表示第k个认知用户SU的感知信道增益,x<sub>PUT</sub>(t)表示在<sub>t</sub>时刻主用户发射机PUT的发射信号,n<sub>k</sub>(t)表示感知信道的加性高斯白噪声,H<sub>1</sub>和H<sub>0</sub>分别表示主用户发射机PUT存在和不存在的两种假设;步骤(2)认知基站对感知信号进行统一观测:记L组长度为N的原始信号{x<sub>i</sub>}<sub>i=1,2,…,L</sub>,把这些信号映射到L组M<sub>i</sub>×1维的观测向量{y<sub>i</sub>}<sub>i=1,2,…,L</sub>,映射的观测矩阵为<img file="FDA0000667308740000013.GIF" wi="267" he="86" />其中R为实数集合;原始信号{x<sub>i</sub>}<sub>i=1,2,…,L</sub>在变换基Ψ上稀疏表示为{s<sub>i</sub>}<sub>i=1,2,…,L</sub>,则y<sub>i</sub>=Φ<sub>i</sub>x<sub>i</sub>+E<sub>i</sub>=Φ<sub>i</sub>Ψ<sub>i</sub>s<sub>i</sub>+E<sub>i</sub>=Θ<sub>i</sub>s<sub>i</sub>+E<sub>i</sub>,i=1,2,…,L其中Θ<sub>i</sub>表示第i个重构任务压缩感知矩阵,每一个由观测向量y<sub>i</sub>恢复得到信号x<sub>i</sub>的过程称为第i个重构任务;步骤(3)通过观测向量采用贝叶斯压缩感知方法进行感知参数估计;所述的贝叶斯压缩感知方法为期望最大化方法或相关向量机模型方法;步骤(4)认知基站基于贝叶斯压缩感知BCS进行稀疏重构和宽带频谱检测:在获得稀疏重构估计向量s={s<sub>i</sub>}<sub>i=1,2,…,L</sub>后,得到原始多任务信号x的估计值x<sup>*</sup>=ψs;当节点感知信道为瑞利衰落信道时,宽带压缩频谱检测的判决门限λ为<img file="FDA0000667308740000021.GIF" wi="358" he="83" />对于重构信号x<sup>*</sup>,归一化重构均方误差MSE为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>MSE</mi><mo>=</mo><mn>10</mn><mi>lg</mi><mo>[</mo><mi>E</mi><mrow><mo>(</mo><mfrac><msubsup><mrow><mo>|</mo><mo>|</mo><msup><mi>x</mi><mo>*</mo></msup><mo>-</mo><mi>x</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><msubsup><mrow><mo>|</mo><mo>|</mo><mi>x</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup></mfrac><mo>)</mo></mrow><mo>]</mo><mrow><mo>(</mo><mi>dB</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000667308740000022.GIF" wi="754" he="180" /></maths>认知基站进行宽带压缩频谱检测得到平均检测概率Pr<sub>d</sub>为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Pr</mi><mi>d</mi></msub><mo>=</mo><mi>Pr</mi><mrow><mo>(</mo><mi>Y</mi><mo>></mo><mi>&lambda;</mi><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mi>&lambda;</mi><mn>2</mn></mfrac></mrow></msup><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>u</mi><mo>-</mo><mn>2</mn></mrow></munderover><mfrac><mn>1</mn><mrow><mi>p</mi><mo>!</mo></mrow></mfrac><msup><mrow><mo>(</mo><mfrac><mi>&lambda;</mi><mn>2</mn></mfrac><mo>)</mo></mrow><mi>p</mi></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>+</mo><mi>&gamma;</mi></mrow><mi>&gamma;</mi></mfrac><mo>)</mo></mrow><mrow><mi>u</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>[</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mi>&lambda;</mi><mrow><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&lambda;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></msup><mo>-</mo><mi>e</mi><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>u</mi><mo>-</mo><mn>2</mn></mrow></munderover><mfrac><mn>1</mn><mrow><mi>p</mi><mo>!</mo></mrow></mfrac><msup><mrow><mo>(</mo><mfrac><mi>&lambda;&gamma;</mi><mrow><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mi>&gamma;</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mi>p</mi></msup><mo>]</mo></mrow>]]></math><img file="FDA0000667308740000023.GIF" wi="1571" he="179" /></maths>其中Pr表示概率,Y表示判决统计量,u为时间带宽积,γ为重构信号接收信噪比;认知基站进行宽带频谱检测时的平均能耗<img file="FDA0000667308740000024.GIF" wi="178" he="82" />为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>E</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>i</mi></msub><mo>,</mo><mi>N</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><msup><msub><mi>&beta;</mi><mn>1</mn></msub><mn>2</mn></msup><mrow><mo>(</mo><msup><mrow><mo>|</mo><msub><mi>h</mi><mn>1</mn></msub><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>PUT</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><msub><mi>CR</mi><mn>1</mn></msub><mn>2</mn></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&beta;</mi><mi>i</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msup><mrow><mo>|</mo><msub><mi>h</mi><mi>i</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><msubsup><mi>&sigma;</mi><mi>PUT</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>CRi</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>]</mo><msub><mi>NT</mi><mi>S</mi></msub><mo>/</mo><mn>2</mn></mrow>]]></math><img file="FDA0000667308740000025.GIF" wi="1227" he="91" /></maths>其中,β<sub>1</sub>和β<sub>i</sub>分别为第1个和第i个认知用户SU的放大转发增益,N为贝叶斯压缩感知BCS重构的采样点数,T<sub>s</sub>为采样间隔,<img file="FDA0000667308740000026.GIF" wi="103" he="77" />为主用户发射机PUT的发射功率,<img file="FDA0000667308740000027.GIF" wi="102" he="78" />和<img file="FDA0000667308740000028.GIF" wi="88" he="86" />分别为第1个和第i个认知用户SU的发射功率;通过数值计算求解能耗优化问题得到平均感知能耗、平均检测概率、归一化重构均方误差MSE之间的关系,该优化问题为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><munder><mi>min</mi><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>,</mo><mi>N</mi></mrow></munder></mtd><mtd><mover><mi>E</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mi>&beta;</mi><mi>i</mi></msub><mo>,</mo><mi>N</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000667308740000029.GIF" wi="353" he="124" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='' close='' separators=''><mtable><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mtd><mtd><msub><mi>Pr</mi><mi>d</mi></msub><mo>&GreaterEqual;</mo><mover><msub><mi>Pr</mi><mi>d</mi></msub><mo>&OverBar;</mo></mover></mtd></mtr></mtable><mo>;</mo></mfenced>]]></math><img file="FDA00006673087400000210.GIF" wi="582" he="89" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>MSE</mi><mo>&le;</mo><mover><mi>MSE</mi><mo>&OverBar;</mo></mover></mrow>]]></math><img file="FDA00006673087400000211.GIF" wi="314" he="72" /></maths>其中,平均感知能耗的目标函数<img file="FDA0000667308740000038.GIF" wi="194" he="74" />约束条件均由上式给出,<img file="FDA0000667308740000039.GIF" wi="85" he="85" />和<img file="FDA00006673087400000310.GIF" wi="134" he="84" />分别表示宽带压缩频谱检测的检测概率阈值和归一化重构均方误差MSE阈值;所述的约束条件包括平均检测概率Pr<sub>d</sub>和归一化重构均方误差MSE。
地址 310018 浙江省杭州市下沙高教园区2号大街