发明名称 一种新型的基于能量检测的双门限协作频谱感知算法
摘要 一种新型的基于能量检测的双门限协作频谱感知算法,它属于无线系统的频谱感知技术领域。它的方法步骤一:考虑一个中心单元或者基站和M个认知用户参与合作频谱感知;中心单元负责管理认知无线电网络和所有的M个合作认知用户;步骤二:每个参与合作频谱感知的认知用户首先检查它所接收到的授权用户信号的能量;判决能量统计量Ti是否小于门限λ<sub>0</sub>或者大于等于门限λ<sub>1</sub>,小于门限λ<sub>0</sub>则认为所感知的授权信道当前空闲,判决为“0”;大于等于门限λ<sub>1</sub>,则认为授权用户正在使用当前频段。本发明通过将K秩融合策略和EGC软融合策略的有效结合,有效提升合作频谱感知效率,节省系统开销,相比较单用户自适应双门限能量检测,多融合准则的双门限协作频谱感知算法检测。
申请公布号 CN105227253A 申请公布日期 2016.01.06
申请号 CN201510515201.X 申请日期 2015.08.20
申请人 黑龙江科技大学 发明人 江晓林;苗雨;谢玉鹏;夏洪洋;康辉
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 哈尔滨龙科专利代理有限公司 23206 代理人 高媛
主权项 一种新型的基于能量检测的双门限协作频谱感知算法,其特征在于它的方法步骤为:步骤一:考虑一个中心单元或者基站和M个认知用户参与合作频谱感知;中心单元负责管理认知无线电网络和所有的M个合作认知用户;在此网络中假设每个认知用户进行本地感知的过程是相互独立的,均采用双门限能量频谱感知算法来进行;第i个合作认知用户的感知情况,检测模型表示如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>{</mo><mtable><mtr><mtd><mrow><msub><mi>H</mi><mn>0</mn></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>H</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>s</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>,</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000011.GIF" wi="1308" he="162" /></maths>式中,s(n)为输入的授权用户信息,假设服从均值为零,方差为<img file="FDA0000785084580000012.GIF" wi="66" he="79" />的正态分布,w(n)为高斯白噪声,服从均值为零,方差为<img file="FDA0000785084580000013.GIF" wi="69" he="83" />的正态分布,H<sub>0</sub>表示测试端只有噪声,H<sub>1</sub>表示测试端有授权用户信息和噪声,假设噪声采样值相互独立同分布,并且与信号采样值也相互独立,则在H<sub>0</sub>下<img file="FDA0000785084580000014.GIF" wi="360" he="79" />在H<sub>1</sub>下<img file="FDA0000785084580000015.GIF" wi="483" he="79" />也即在N个抽样求和后,输出统计量T<sub>i</sub>可以表示为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><msup><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000016.GIF" wi="1060" he="149" /></maths>协作频谱检测充分利用各个认知用户感知数据信息,显著提高频谱感知结果的可靠性能,因此我们采用了基于能量检测的双门限的协作频谱感知策略;步骤二:每个参与合作频谱感知的认知用户首先检查它所接收到的授权用户信号的能量;判决能量统计量Ti是否小于门限λ<sub>0</sub>或者大于等于门限λ<sub>1</sub>,小于门限λ<sub>0</sub>则认为所感知的授权信道当前空闲,判决为“0”;大于等于门限λ<sub>1</sub>,则认为授权用户正在使用当前频段,判决为“1”;我们用D<sub>i</sub>表示第i个用户的本地感知结果,D<sub>i</sub>取数据“0”或“1”,判决结果D<sub>i</sub>将送入融合中心,融合中心根据本地的感知结果,采用自适应K秩硬融合策略:下面对对K秩硬融合准则判决进行分析,假设一共有M个认知用户参与协作感知,基于能量感知的双门限算法的认知无线网络中,其k秩融合策略联合检测概率和虚警概率表达式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>D</mi><mo>_</mo><mi>d</mi><mi>o</mi><mi>u</mi></mrow></msub><mo>=</mo><mi>P</mi><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>></mo><mi>k</mi><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mi>k</mi></mrow><mi>K</mi></munderover><mfenced open='(' close=')'><mtable><mtr><mtd><mi>K</mi></mtd></mtr><mtr><mtd><mi>i</mi></mtd></mtr></mtable></mfenced><msup><mrow><mo>(</mo><mi>P</mi><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>></mo><mi>i</mi><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo><mo>)</mo></mrow><mi>i</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>P</mi><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>></mo><mi>i</mi><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo><mo>)</mo></mrow><mrow><mi>K</mi><mo>-</mo><mi>i</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000021.GIF" wi="1591" he="162" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>F</mi><mo>_</mo><mi>d</mi><mi>o</mi><mi>u</mi></mrow></msub><mo>=</mo><mi>P</mi><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>&gt;</mo><mi>k</mi><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mi>k</mi></mrow><mi>K</mi></munderover><mfenced open = '(' close = ')'><mtable><mtr><mtd><mi>K</mi></mtd></mtr><mtr><mtd><mi>i</mi></mtd></mtr></mtable></mfenced><msup><mrow><mo>(</mo><mi>P</mi><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>&gt;</mo><mi>k</mi><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo><mo>)</mo></mrow><mi>i</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>P</mi><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>T</mi><mi>i</mi></msub><mo>&gt;</mo><mi>k</mi><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo><mo>)</mo></mrow><mrow><mi>K</mi><mo>-</mo><mi>i</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000022.GIF" wi="1571" he="162" /></maths>k=1即为“或”融合,k=M即为“与”融合;传统的K秩融合投票门限值固定,一般为<img file="FDA0000785084580000023.GIF" wi="287" he="83" />其中<img file="FDA0000785084580000024.GIF" wi="78" he="82" />为ceiling函数,此时称为Majority准则;对频谱状态作出判决,判决结果用R<sub>1</sub>来表示,表达式可写为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>R</mi><mn>1</mn></msub><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>D</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><mi>K</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>D</mi><mi>i</mi></msub><mo>&lt;</mo><mi>K</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000025.GIF" wi="1221" he="298" /></maths>若R<sub>1</sub>为1,表示融合中心判决授权用户存在,授权用户不可以占用信道;若R<sub>1</sub>为0,判决授权用户不存在,授权用户可以接入信道进行通信;步骤三:如果Ti落入由门限λ<sub>0</sub>和λ<sub>1</sub>构成的不可靠区间,即λ<sub>0</sub>≤Ti&lt;λ<sub>1</sub>,则第i个合作认知用户将自己接收到的信号直接上传到融合中心,在融合中心采用EGC软融合策略数据进行加权求和,并与软融合门限值进行比较,并作出判决,采用EGC准则的融合中心的检测过程可以表示为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&psi;</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><msub><mi>T</mi><mi>i</mi></msub><msqrt><mi>M</mi></msqrt></mfrac><mo>&GreaterEqual;</mo><msub><mi>&lambda;</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&psi;</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><msub><mi>T</mi><mi>i</mi></msub><msqrt><mi>M</mi></msqrt></mfrac><mo>&lt;</mo><msub><mi>&lambda;</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000026.GIF" wi="1315" he="298" /></maths>其中,λ<sub>EGC</sub>为预设门限值;文献[6]指出ψ<sub>EGC</sub>服从正态分布:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&psi;</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>N</mi><mi>o</mi><mi>r</mi><mi>m</mi><mi>a</mi><mi>l</mi><mrow><mo>(</mo><mi>N</mi><msqrt><mi>M</mi></msqrt><mo>,</mo><mn>2</mn><mi>N</mi><mo>)</mo></mrow></mrow></mtd><mtd><msub><mi>H</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mrow><mi>N</mi><mi>o</mi><mi>r</mi><mi>m</mi><mi>a</mi><mi>l</mi><mrow><mo>(</mo><mfrac><mi>N</mi><msqrt><mi>M</mi></msqrt></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd><mtd><msub><mi>H</mi><mn>1</mn></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000031.GIF" wi="1398" he="248" /></maths>其中,γ<sub>i</sub>表示第i个用户的接收信噪比,N表示采样点数,M表示认知用户数;所以联合检测概率和虚警概率为:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>D</mi><mo>_</mo><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>=</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>f</mi><mo>_</mo><mi>e</mi><mi>g</mi><mi>c</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msqrt><mfrac><mi>N</mi><mrow><mn>2</mn><mi>M</mi></mrow></mfrac></msqrt><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&gamma;</mi><mi>i</mi></msub></mrow><msqrt><mrow><mn>1</mn><mo>+</mo><mfrac><mn>2</mn><mi>M</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&gamma;</mi><mi>i</mi></msub></mrow></msqrt></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000032.GIF" wi="1399" he="279" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>Q</mi><mrow><mi>F</mi><mo>_</mo><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>=</mo><mi>Q</mi><mrow><mo>(</mo><msqrt><mrow><mn>1</mn><mo>+</mo><mfrac><mn>2</mn><mi>M</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&gamma;</mi><mi>i</mi></msub></mrow></msqrt><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>(</mo><msub><mi>Q</mi><mrow><mi>D</mi><mo>_</mo><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msqrt><mfrac><mi>N</mi><mrow><mn>2</mn><mi>M</mi></mrow></mfrac></msqrt><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000033.GIF" wi="1433" he="179" /></maths>式中,<img file="FDA0000785084580000034.GIF" wi="193" he="143" />为信息融合中心瞬时信噪比;采用EGC准则,参与的M个认知节点将本地的感知观测值Ti直接传送给融合中心,融合中心再对数据进行加权求和,若各节点权值取<img file="FDA0000785084580000035.GIF" wi="179" he="78" />融合中心接收的信号:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>&psi;</mi><mrow><mi>E</mi><mi>G</mi><mi>C</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mfrac><msub><mi>T</mi><mi>i</mi></msub><msqrt><mi>M</mi></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000036.GIF" wi="1110" he="152" /></maths>将判决结果传送给用户,若判决结果用R<sub>2</sub>来表示,表达式可写为判决规则如下:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>R</mi><mn>2</mn></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>Y</mi><mo>&GreaterEqual;</mo><mi>&xi;</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>Y</mi><mo>&lt;</mo><mi>&xi;</mi></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000785084580000037.GIF" wi="1219" he="156" /></maths>其中,<img file="FDA0000785084580000038.GIF" wi="400" he="142" />若R<sub>2</sub>为1,表示融合中心判决授权用户存在,认知用户不可以占用信道;若R<sub>2</sub>为0,表示融合中心判决授权用户不存在,认知用户可以接入信道进行通信。
地址 150027 黑龙江省哈尔滨市松北区浦源路2468号