发明名称 基于可信度的协同频谱感知方法
摘要 本发明的基于可信度的协同频谱感知方法,包括初始化、未授权用户本地感知、本地判决、感知信息量化、可信度融合、判决融合等步骤。本发明在先验信息未知的条件下给出了未授权用户本地感知的可信度,采用基站与未授权用户联合的混合门限判决策略,并通过量化本地感知结果来减小协同感知的网络开销,具有很好的感知性能。
申请公布号 CN101815305A 申请公布日期 2010.08.25
申请号 CN201010109711.4 申请日期 2010.02.10
申请人 中国人民解放军理工大学 发明人 王金龙;吴启晖;郑学强;沈良;王呈贵;徐以涛
分类号 H04W16/14(2009.01)I;H04W24/00(2009.01)I;H04B17/00(2006.01)I 主分类号 H04W16/14(2009.01)I
代理机构 南京天华专利代理有限责任公司 32218 代理人 夏平
主权项 1.一种基于可信度的协同频谱感知方法,其特征是它包括以下步骤:步骤一:未授权用户在本地使用基于能量检测的频谱感知方法进行感知得到感知结果<maths num="0001"><![CDATA[<math><mrow><mi>&theta;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn><mi>u</mi></mrow></munderover><msup><mrow><mo>|</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow></math>]]></maths>其中y是未授权用户的接收信号,θ是未授权用户的本地感知结果,u是时间带宽积;步骤二:(a).设定本方法的本地虚警概率Q<sub>f</sub><sup>C</sup>,选择软信息融合方法的虚警概率Q<sub>f</sub><sup>S</sup>,根据<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>Q</mi><mi>f</mi><mi>s</mi></msubsup><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>Nu</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>Nu</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>计算软信息融合方法的判决门限λ<sub>0</sub>,其中<maths num="0003"><![CDATA[<math><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mo>&infin;</mo></msubsup><msup><mi>e</mi><mrow><mo>-</mo><mi>t</mi></mrow></msup><msup><mi>t</mi><mrow><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>dt</mi><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow></mrow></mfrac><mo>[</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>x</mi></msubsup><msup><mi>e</mi><mrow><mo>-</mo><mi>t</mi></mrow></msup><msup><mi>t</mi><mrow><mi>a</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>dt</mi><mo>]</mo><mo>,</mo></mrow></math>]]></maths>N是系统中参与协同频谱感知的未授权用户数,然后检验Q<sub>f</sub><sup>S</sup>是否满足以下两个条件;<maths num="0005"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>Q</mi><mi>f</mi><mi>S</mi></msubsup><mo>&lt;</mo><msubsup><mi>Q</mi><mi>f</mi><mi>C</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>Q</mi><mi>f</mi><mi>OR</mi></msubsup><mrow><mo>(</mo><msup><mi>&lambda;</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>></mo><msubsup><mi>Q</mi><mi>f</mi><mi>C</mi></msubsup></mtd></mtr></mtable></mfenced></math>]]></maths>其中Q<sub>f</sub><sup>OR</sup>(λ′)是硬信息融合方法的虚警概率,<maths num="0006"><![CDATA[<math><mrow><msubsup><mi>Q</mi><mi>f</mi><mi>OR</mi></msubsup><mrow><mo>(</mo><msup><mi>&lambda;</mi><mo>&prime;</mo></msup><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>[</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><msup><mi>&lambda;</mi><mo>&prime;</mo></msup><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mrow><mi>N</mi></msup><mo>,</mo></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msup><mi>&lambda;</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac></mrow></math>]]></maths>是硬信息融合方法的本地判决门限;若Q<sub>f</sub><sup>S</sup>满足检验条件,转步骤(b);否则,转步骤(a)重新选择选择软信息融合方法的虚警概率Q<sub>f</sub><sup>S</sup>;(b).根据设定的本地虚警概率Q<sub>f</sub><sup>C</sup>和软信息融合方法的判决门限λ<sub>0</sub>计算未授权用户本地感知判决门限的均值为<img file="FSA00000029686800018.GIF" wi="55" he="108" />和基站可信度融合门限<img file="FSA00000029686800019.GIF" wi="116" he="109" />其中L是参与可信度融合的未授权用户数;得到未授权用户的两个本地判决门限分别为<maths num="0008"><![CDATA[<math><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mo>-</mo><mi>&Delta;&lambda;</mi><mo>,</mo></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mo>+</mo><mi>&Delta;&lambda;</mi><mo>,</mo></mrow></math>]]></maths>Δλ为两个本地判决门限λ<sub>1</sub>,λ<sub>2</sub>与均值<img file="FSA000000296868000112.GIF" wi="56" he="108" />之差,它满足:<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>Q</mi><mi>f</mi><mi>C</mi></msubsup><mo>=</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>[</mo><mi>F</mi><mrow><mo>(</mo><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>]</mo></mrow><mi>N</mi></msup><mo>+</mo><mfrac><mn>1</mn><msup><mrow><mo>[</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>N</mi></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>K</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mfenced open='(' close=')'><mtable><mtr><mtd><mi>N</mi></mtd></mtr><mtr><mtd><mi>K</mi></mtd></mtr></mtable></mfenced><msup><mrow><mo>[</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><msub><mi>&lambda;</mi><mn>1</mn></msub><mn>2</mn></mfrac><mo>)</mo></mrow><mo>-</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><msub><mi>&lambda;</mi><mn>2</mn></msub><mn>2</mn></mfrac><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></msup><msup><mrow><mo>[</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>-</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><msub><mi>&lambda;</mi><mn>1</mn></msub><mn>2</mn></mfrac><mo>)</mo></mrow><mo>]</mo></mrow><mi>K</mi></msup><mi>P</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0011"><![CDATA[<math><mrow><mo>=</mo><mn>1</mn><mo>-</mo><msup><mrow><mo>[</mo><mi>F</mi><mrow><mo>(</mo><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mo>+</mo><mi>&Delta;&lambda;</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>N</mi></msup><mo>+</mo><mfrac><mn>1</mn><msup><mrow><mo>[</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>]</mo></mrow><mi>N</mi></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>K</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mfenced open='(' close=')'><mtable><mtr><mtd><mi>N</mi></mtd></mtr><mtr><mtd><mi>K</mi></mtd></mtr></mtable></mfenced><msup><mrow><mo>[</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><mrow><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mo>-</mo><mi>&Delta;&lambda;</mi></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>-</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><mrow><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mo>+</mo><mi>&Delta;&lambda;</mi></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></msup><msup><mrow><mo>[</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>-</mo><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><mrow><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mo>-</mo><mi>&Delta;&lambda;</mi></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>]</mo></mrow><mi>K</mi></msup><mi>P</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中<maths num="0012"><![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><msup><mi>&sigma;</mi><mrow><mn>2</mn><mi>u</mi></mrow></msup><msup><mn>2</mn><mi>u</mi></msup><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></msup><munder><mrow><mo>&Integral;</mo><mo>&Integral;</mo><mo>&Integral;</mo></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></munderover><msub><mi>&theta;</mi><mi>i</mi></msub><mo>></mo><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>K</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>&lt;</mo><msub><mi>&lambda;</mi><mn>2</mn></msub></mrow></munder><msup><mrow><mo>(</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></munderover><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mi>u</mi><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>e</mi><mrow><mfrac><mn>1</mn><msup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mn>2</mn></msup></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></munderover><msub><mi>&theta;</mi><mi>i</mi></msub></mrow></msup><mi>d</mi><msub><mi>&theta;</mi><mn>1</mn></msub><mi>d</mi><msub><mi>&theta;</mi><mn>2</mn></msub><mo>.</mo><mo>.</mo><mo>.</mo><mi>d</mi><msub><mi>&theta;</mi><mrow><mi>N</mi><mo>-</mo><mi>K</mi></mrow></msub><mo>,</mo></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><mi>F</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>,</mo><mfrac><mi>x</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>σ<sup>2</sup>为未授权用户的本地噪声方差;步骤三:未授权用户进行本地判决x=Φ(θ),判决规则如下:<img file="FSA00000029686800024.GIF" wi="864" he="212" />其中H<sub>1</sub>代表授权用户存在,H<sub>0</sub>代表授权用户不存在,<maths num="0014"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></mfrac></mrow></math>]]></maths>是在本地感知结果θ落在本地判决门限λ<sub>1</sub>、λ<sub>2</sub>之间区域时对信号存在的可信度;步骤四:未授权用户将本地感知和判决的处理结果通过报告信道传给基站,(a).若本地判决为H<sub>1</sub>或H<sub>0</sub>,则直接将1比特判决结果即0或1发送给基站;(b).若本地判决得到的是本地感知可信度α<sub>i</sub>,则将可信度进行M比特量化后再上传;其中,可信度M比特量化即<maths num="0015"><![CDATA[<math><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>=</mo><mi>quantize</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow><mrow><msub><mi>&lambda;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&lambda;</mi><mn>1</mn></msub></mrow></mfrac><mo>,</mo><mi>M</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>quantize(X,M)表示对X进行M比特量化;步骤五:基站对收到的未授权用户本地判决结果x<sub>i</sub>按照长度进行分类并进行可信度融合,即基站将接收到的结果按照数据长度分为两类,即1比特数据b和M比特数据α,然后对所有M比特的可信度进行合并得到可信度融合结果ω,<maths num="0016"><![CDATA[<math><mrow><mi>&omega;</mi><mo>=</mo><mi>&gamma;</mi><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mfrac><mrow><msup><mn>1</mn><mi>T</mi></msup><msup><mi>&Sigma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>&theta;</mi><mo>&prime;</mo></msup></mrow><mrow><msup><mn>1</mn><mi>T</mi></msup><msup><mi>&Sigma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mn>1</mn></mrow></mfrac><mo>&lt;</mo><mi>L</mi><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mfrac><mrow><msup><mn>1</mn><mi>T</mi></msup><msup><mi>&Sigma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><msup><mi>&theta;</mi><mo>&prime;</mo></msup></mrow><mrow><msup><mn>1</mn><mi>T</mi></msup><msup><mi>&Sigma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mn>1</mn></mrow></mfrac><mo>&GreaterEqual;</mo><mi>L</mi><mfrac><msub><mi>&lambda;</mi><mn>0</mn></msub><mi>N</mi></mfrac></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中θ′=λ<sub>2</sub>-(λ<sub>2</sub>-λ<sub>1</sub>)*α,∑为噪声协方差矩阵,即<maths num="0017"><![CDATA[<math><mrow><mi>&Sigma;</mi><mo>=</mo><mi>diag</mi><mrow><mo>(</mo><msubsup><mi>&sigma;</mi><mn>1</mn><mn>2</mn></msubsup><mo>,</mo><msubsup><mi>&sigma;</mi><mn>2</mn><mn>2</mn></msubsup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>&sigma;</mi><mi>N</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>σ<sub>1</sub><sup>2</sup>,σ<sub>2</sub><sup>2</sup>,...,σ<sub>N</sub><sup>2</sup>为未授权用户的本地噪声方差,<img file="FSA00000029686800032.GIF" wi="91" he="109" />是可信度融合门限;步骤六:基站将可信度融合判决结果ω与其他所有1比特信息b使用OR准则进行判决融合得到融合结果x<sub>0</sub>,即x<sub>0</sub>=γ(ω,b),<img file="FSA00000029686800033.GIF" wi="807" he="178" />,若多个未授权用户的感知结果中有任何一个宣布授权用户存在即融合结果x<sub>0</sub>为1,则融合结果判决为授权用户存在。
地址 210007 江苏省南京市白下区御道街标营2号