发明名称 一种认知无线电系统中的合作频谱感知方法
摘要 本发明公开了一种认知无线电中的合作频谱感知方法,基于提出的一种改进的能量检测算法,根据接收到的主用户发送的信号yi,将其幅值进行p次幂的运算,从而得到一个本地感知量。然后根据推导得出的本地感知量的分布特性,计算出本地的感知结果。每个认知用户将其本地感知结果发送给融合中心,融合中心使用证据理论进行融合,得出主用户是否出现的总的判决结果。实验表明:本发明的检测概率明显高于基于传统能量检测算法的“或融合”、“与融合”感知方法,体现出了本发明的优越性能。与此同时,本发明的总体的错误概率(漏检概率与虚警概率之和)明显低于“或融合”、“与融合”感知方法,这也体现了本发明具有较高的可靠性。
申请公布号 CN102291188B 申请公布日期 2013.09.25
申请号 CN201110222782.X 申请日期 2011.08.04
申请人 电子科技大学 发明人 刘健;隆克平;陈红宇
分类号 H04B17/00(2006.01)I;H04W24/00(2009.01)I 主分类号 H04B17/00(2006.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 1.一种认知无线电系统中的合作频谱感知方法,其特征在于,包括以下步骤:(1)、认知用户根据频段内接收到的信号y<sub>i</sub>,将其幅值进行p次幂的运算:z<sub>i</sub>=|y<sub>i</sub>|<sup>p</sup>其中,i=0,1,…,N,N为认知用户数量,p为一个大于零的常数;(2)、计算每个用户的本地感知结果:主用户没有出现的感知结果为m<sub>i</sub>(H<sub>0</sub>)、主用户出现的感知结果为m<sub>i</sub>(H<sub>1</sub>)、对于主用户是否出现不确定的结果为m<sub>i</sub>(Ω):<maths num="0001"><![CDATA[<math><mrow><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>z</mi><mi>i</mi></msub><mrow><mo>+</mo><mo>&infin;</mo></mrow></munderover><mfrac><mrow><msqrt><mn>2</mn></msqrt><msup><mi>x</mi><mfrac><mrow><mn>1</mn><mo>-</mo><mi>p</mi></mrow><mi>p</mi></mfrac></msup><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>x</mi><mfrac><mn>2</mn><mi>p</mi></mfrac></msup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><mi>p</mi><msqrt><mi>&pi;</mi></msqrt><msub><mi>&sigma;</mi><mi>n</mi></msub></mrow></mfrac><mi>dx</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><msub><mi>z</mi><mi>i</mi></msub></munderover><mfrac><mrow><msqrt><mn>2</mn></msqrt><msup><mi>x</mi><mfrac><mrow><mn>1</mn><mo>-</mo><mi>p</mi></mrow><mi>p</mi></mfrac></msup><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mi>x</mi><mfrac><mn>2</mn><mi>p</mi></mfrac></msup><mrow><mn>2</mn><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>s</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><mi>p</mi><msqrt><mi>&pi;</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>s</mi><mn>2</mn></msubsup><mo>)</mo></mrow></msqrt></mrow></mfrac><mi>dx</mi></mrow></math>]]></maths>m<sub>i</sub>(Ω)=1-m<sub>i</sub>(H<sub>0</sub>)-m<sub>i</sub>(H<sub>1</sub>)通过换元变换,上述式子中m<sub>i</sub>(H<sub>0</sub>)和m<sub>i</sub>(H<sub>1</sub>)的计算公式可以变换为:<maths num="0003"><![CDATA[<math><mrow><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>&pi;</mi></msqrt></mfrac><munderover><mo>&Integral;</mo><mfrac><msubsup><mi>z</mi><mi>i</mi><mrow><mn>2</mn><mo>/</mo><mi>p</mi></mrow></msubsup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow></mfrac><mrow><mo>+</mo><mo>&infin;</mo></mrow></munderover><msqrt><mi>t</mi></msqrt><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>&pi;</mi></msqrt></mfrac><mi>&Gamma;</mi><mrow><mo>(</mo><mfrac><msubsup><mi>z</mi><mi>i</mi><mrow><mn>2</mn><mo>/</mo><mi>p</mi></mrow></msubsup><mrow><mn>2</mn><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow></mfrac><mo>,</mo><mn>3</mn><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>&pi;</mi></msqrt></mfrac><munderover><mo>&Integral;</mo><mn>0</mn><mfrac><msubsup><mi>z</mi><mi>i</mi><mrow><mn>2</mn><mo>/</mo><mi>p</mi></mrow></msubsup><mrow><mn>2</mn><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>s</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac></munderover><msqrt><mi>v</mi></msqrt><mo>&CenterDot;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>v</mi><mo>)</mo></mrow><mi>dv</mi><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>&pi;</mi></msqrt></mfrac><mo>[</mo><mn>1</mn><mo>-</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mfrac><msubsup><mi>z</mi><mi>i</mi><mrow><mn>2</mn><mo>/</mo><mi>p</mi></mrow></msubsup><mrow><mn>2</mn><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>&sigma;</mi><mi>s</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mn>3</mn><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>其中<img file="FDA00003416824400015.GIF" wi="534" he="111" />表示不完全gamma函数,H<sub>0</sub>表示主用户没有出现,H<sub>1</sub>表示主用户出现,σ<sub>n</sub>表示频段中的零均值加性高斯白噪声的标准差,σ<sub>s</sub>表示主用户出现的时候所发送的信号的标准差;(3)、每个认知用户将其感知计算得到的感知结果m<sub>i</sub>(H<sub>0</sub>)、m<sub>i</sub>(H<sub>1</sub>)以及m<sub>i</sub>(Ω)发送给融合中心;融合中心接收所有认知用户发送过来的感知结果,使用证据理论进行融合,得到一个总的感知结果:<maths num="0005"><![CDATA[<math><mrow><mi>m</mi><mrow><mo>(</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>A</mi><mi>N</mi></msub><mo>=</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&Omega;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>上式中,<img file="FDA00003416824400022.GIF" wi="432" he="151" />表示A<sub>i</sub>取值为H<sub>0</sub>或者Ω的各种组合下的N个m<sub>i</sub>(A<sub>i</sub>),i=0,1,…,N的乘积之和;m(H<sub>0</sub>)表示了主用户没有出现的融合感知结果;<maths num="0006"><![CDATA[<math><mrow><mi>m</mi><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mrow><msub><mi>B</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>B</mi><mn>2</mn></msub><mo>&cap;</mo><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>B</mi><mi>N</mi></msub><mo>=</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>B</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&Omega;</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>上式中,<img file="FDA00003416824400024.GIF" wi="431" he="151" />表示B<sub>i</sub>取值为H<sub>1</sub>或者Ω的各种组合下的N个m<sub>i</sub>(B<sub>i</sub>),i=0,1,…,N的乘积之和;m(H<sub>1</sub>)表示了主用户出现的融合感知结果;(4)、融合中心比较m(H<sub>0</sub>)和m(H<sub>1</sub>),如果m(H<sub>1</sub>)&gt;m(H<sub>0</sub>),则认为主用户存在,如果m(H<sub>1</sub>)≤m(H<sub>0</sub>),则认为主用户没有出现。
地址 611731 四川省成都市高新区(西区)西源大道2006号