发明名称 一种自适应调整的协同频谱感知方法
摘要 本发明涉及自适应调整的协同频谱感知方法,依次包括认知用户获取到能量检测判决的最优门限值后,分别将能量检测结果及自身信噪比发送给频谱感知融合中心;频谱感知融合中心分别统计感知到授权用户频谱为占用状态、空闲状态的认知用户数目,计算认知用户的感知结果及诚信系数,计算授权用户频谱为占用状态的平均检测概率、全局检测概率、全局漏检概率,以及授权用户频谱为空闲状态的平均检测概率、全局检测概率、全局漏检概率和全局虚警概率;建立基于认知用户数目的频谱感知误差函数后,以频谱感知误差函数最小值对应的数值作为最佳协同认知用户数目,降低频谱感知融合中心融合计算量,并在获取到最佳协同认知用户后完成协同频谱感知。
申请公布号 CN105471528A 申请公布日期 2016.04.06
申请号 CN201510833459.4 申请日期 2015.11.25
申请人 宁波大学 发明人 郑紫微;秦闯
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 宁波诚源专利事务所有限公司 33102 代理人 刘凤钦
主权项 一种自适应调整的协同频谱感知方法,用于授权用户、频谱感知融合中心以及N个认知用户组成的认知无线网络中,其特征在于,依次包括如下步骤:(1)N个认知用户对授权用户的频谱占用情况进行本地能量检测,并将频谱感知结果以及认知用户自身信噪比分别发送给频谱感知融合中心;其中,认知用户标记为CR<sub>i</sub>(i=1,2,…,N),认知用户CR<sub>i</sub>自身信噪比为snr<sub>i</sub>,授权用户标记为PU,频谱感知融合中心标记为FC;(2)频谱感知融合中心FC根据N个认知用户发送的频谱感知结果,统计N个认知用户中感知到授权用户PU频谱为占用状态的认知用户数目为m(1≤m≤N)、感知到授权用户PU频谱为空闲状态的认知用户数目为N‑m;其中,授权用户PU频谱为占用状态记为H<sub>1</sub>,授权用户PU频谱为空闲状态记为H<sub>0</sub>;(3)频谱感知融合中心FC根据N个认知用户发送的信噪比,计算m个感知到授权用户PU频谱为占用状态H<sub>1</sub>的认知用户诚信系数κ<sub>1,j</sub>以及N‑m个感知到授权用户PU频谱为空闲状态H<sub>0</sub>的认知用户诚信系数κ<sub>2,t</sub>;其中,诚信系数κ<sub>1,j</sub>以及κ<sub>2,t</sub>的计算公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&kappa;</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><msubsup><mi>snr</mi><mi>j</mi><mn>2</mn></msubsup></mrow><mrow><mfrac><mn>1</mn><mi>m</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>snr</mi><mi>j</mi><mn>2</mn></msubsup></mrow></mfrac><mo>,</mo><msub><mi>&kappa;</mi><mrow><mn>2</mn><mo>,</mo><mi>t</mi></mrow></msub><mo>=</mo><mfrac><mrow><msup><msub><mi>snr</mi><mi>t</mi></msub><mn>2</mn></msup></mrow><mrow><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mi>m</mi></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>m</mi></mrow></munderover><msup><msub><mi>snr</mi><mi>t</mi></msub><mn>2</mn></msup></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000011.GIF" wi="853" he="223" /></maths>(4)频谱感知融合中心FC根据m个认知用户的各自感知结果以及诚信系数κ<sub>1,j</sub>,分别计算授权用户PU频谱为占用状态H<sub>1</sub>的平均检测概率<img file="FDA0000857652750000012.GIF" wi="151" he="61" />全局检测概率<img file="FDA0000857652750000013.GIF" wi="134" he="55" />和此占用状态H<sub>1</sub>对应的全局漏检概率<img file="FDA0000857652750000014.GIF" wi="183" he="63" />以及授权用户PU的频谱为空闲状态H<sub>0</sub>的平均检测概率<img file="FDA0000857652750000015.GIF" wi="159" he="61" />全局检测概率<img file="FDA0000857652750000016.GIF" wi="165" he="61" />此空闲状态H<sub>0</sub>对应的全局漏检概率<img file="FDA0000857652750000017.GIF" wi="160" he="61" />和全局虚警概率<img file="FDA0000857652750000018.GIF" wi="173" he="62" />其中,该过程包括如下步骤(4‑1)至步骤(4‑3):(4‑1)建立m个认知用户协同感知的全局错误检测概率P<sub>e</sub>,获取关于决策门限的能量检测优化函数γ<sup>*</sup>以及能量检测的最优门限值γ<sub>opt</sub>,并计算授权用户PU频谱为占用状态H<sub>1</sub>的平均检测概率<img file="FDA0000857652750000019.GIF" wi="143" he="55" />其中,m个认知用户协同感知的全局错误检测概率P<sub>e</sub>计算公式如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mi>e</mi></msub><mo>=</mo><msub><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub></msub><msub><mi>P</mi><mi>f</mi></msub><mo>+</mo><msub><mi>P</mi><msub><mi>H</mi><mn>1</mn></msub></msub><msub><mi>P</mi><mi>m</mi></msub><mo>,</mo><msub><mi>P</mi><msub><mi>H</mi><mn>1</mn></msub></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub></msub><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000027.GIF" wi="739" he="62" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>f</mi></msub><mo>=</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><mi>&gamma;</mi><mo>-</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mrow><mfrac><mn>2</mn><mi>m</mi></mfrac><msubsup><mi>&sigma;</mi><mi>n</mi><mn>4</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msub><mi>P</mi><mi>d</mi></msub><mo>=</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><mi>&gamma;</mi><mo>-</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover></mrow><mo>)</mo></mrow><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mrow><mfrac><mn>2</mn><mi>m</mi></mfrac><mrow><mo>(</mo><mrow><mn>2</mn><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><msubsup><mi>&sigma;</mi><mi>n</mi><mn>4</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>;</mo><msub><mi>P</mi><mi>m</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>d</mi></msub><mo>;</mo><mi>Q</mi><mrow><mo>(</mo><mi>z</mi><mo>)</mo></mrow><mo>=</mo><mrow><msubsup><mo>&Integral;</mo><mi>z</mi><mi>&infin;</mi></msubsup><mrow><mfrac><mn>1</mn><msqrt><mrow><mn>2</mn><mi>&pi;</mi></mrow></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>x</mi><mn>2</mn></msup></mrow></msup><mi>d</mi><mi>x</mi></mrow></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000021.GIF" wi="1726" he="263" /></maths>其中,<img file="FDA0000857652750000028.GIF" wi="89" he="56" />为授权用户PU频谱处于空闲状态H<sub>0</sub>的概率,<img file="FDA0000857652750000029.GIF" wi="84" he="55" />为授权用户PU频谱处于占用状态H<sub>1</sub>的概率;P<sub>f</sub>为全局虚警概率,P<sub>d</sub>为全局检测概率,P<sub>m</sub>为全局漏检概率;<img file="FDA0000857652750000022.GIF" wi="85" he="71" />为对应处于授权用户PU频谱处于占用状态H<sub>1</sub>的m个认知用户的平均信噪比,Q(z)表示正态高斯互补积分函数;其中,<img file="FDA0000857652750000023.GIF" wi="349" he="135" />snr<sub>i</sub>为认知用户CR<sub>i</sub>自身的信噪比;关于决策门限的能量检测优化函数γ<sup>*</sup>定义为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mi>&gamma;</mi><mo>*</mo></msup><mo>=</mo><mi>arg</mi><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mi>&gamma;</mi></munder><msub><mi>P</mi><mi>e</mi></msub><mo>=</mo><msub><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub></msub><mo>&CenterDot;</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><mi>&gamma;</mi><mo>-</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mrow><mfrac><mn>2</mn><mi>m</mi></mfrac><msubsup><mi>&sigma;</mi><mi>n</mi><mn>4</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mi>P</mi><msub><mi>H</mi><mn>1</mn></msub></msub><mo>&CenterDot;</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><mi>&gamma;</mi><mo>-</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>+</mo><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover></mrow><mo>)</mo></mrow><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup></mrow><mrow><mfrac><mn>2</mn><mi>m</mi></mfrac><mrow><mo>(</mo><mrow><mn>2</mn><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><msubsup><mi>&sigma;</mi><mi>n</mi><mn>4</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000024.GIF" wi="1237" he="254" /></maths>能量检测的最优门限值γ<sub>opt</sub>为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&gamma;</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo>=</mo><mi>&gamma;</mi><msub><mo>|</mo><mrow><mfrac><mrow><mo>&part;</mo><msub><mi>P</mi><mi>e</mi></msub></mrow><mrow><mo>&part;</mo><mi>&gamma;</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mfrac><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><mn>2</mn></mfrac><mo>+</mo><msubsup><mi>&sigma;</mi><mi>n</mi><mn>2</mn></msubsup><msqrt><mrow><mfrac><mn>1</mn><mn>4</mn></mfrac><mo>+</mo><mfrac><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mn>2</mn></mfrac><mo>+</mo><mfrac><mrow><mn>4</mn><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mo>+</mo><mn>2</mn></mrow><mrow><mi>m</mi><mo>&CenterDot;</mo><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover></mrow></mfrac><mi>ln</mi><mrow><mo>(</mo><mrow><mfrac><msub><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub></msub><msub><mi>P</mi><msub><mi>H</mi><mn>1</mn></msub></msub></mfrac><msqrt><mrow><mn>2</mn><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mo>+</mo><mn>1</mn></mrow></msqrt></mrow><mo>)</mo></mrow></mrow></msqrt></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000025.GIF" wi="1078" he="342" /></maths>授权用户PU频谱为占用状态H<sub>1</sub>的平均检测概率<img file="FDA00008576527500000210.GIF" wi="134" he="60" />计算公式如下:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&gamma;</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><msqrt><mrow><mo>(</mo><mn>2</mn><mo>/</mo><msup><mi>N</mi><mo>&prime;</mo></msup><mo>)</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mover><mrow><mi>s</mi><mi>n</mi><mi>r</mi></mrow><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000026.GIF" wi="661" he="252" /></maths>(4‑2)根据所得权用户PU频谱为占用状态H<sub>1</sub>的平均检测概率<img file="FDA00008576527500000211.GIF" wi="129" he="61" />以及m个认知用户的诚信系数κ<sub>1,j</sub>,计算授权用户PU频谱为占用状态H<sub>1</sub>的全局检测概率<img file="FDA00008576527500000212.GIF" wi="125" he="55" />和此占用状态H<sub>1</sub>对应的全局漏检概率<img file="FDA00008576527500000213.GIF" wi="183" he="55" />其中,全局检测概率<img file="FDA00008576527500000214.GIF" wi="125" he="55" />和全局漏检概率<img file="FDA00008576527500000215.GIF" wi="162" he="57" />计算公式如下:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mroot><mrow><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&kappa;</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub></mrow><mi>m</mi></mroot><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mi>m</mi></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>)</mo></mrow><mi>l</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>l</mi></mrow></msup><mo>;</mo><msub><mi>D</mi><mrow><mi>u</mi><mi>n</mi><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>D</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000031.GIF" wi="1324" he="165" /></maths>(4‑3)根据所得授权用户PU频谱为空闲状态H<sub>0</sub>的平均检测概率<img file="FDA0000857652750000035.GIF" wi="126" he="55" />以及N‑m个认知用户的诚信系数κ<sub>2,t</sub>,计算授权用户PU频谱为空闲状态H<sub>0</sub>的全局检测概率<img file="FDA0000857652750000036.GIF" wi="132" he="62" />和此空闲状态H<sub>0</sub>对应的全局漏检概率<img file="FDA0000857652750000037.GIF" wi="190" he="61" />全局虚警概率<img file="FDA0000857652750000038.GIF" wi="174" he="62" />其中,平均检测概率<img file="FDA0000857652750000039.GIF" wi="155" he="55" />全局检测概率<img file="FDA00008576527500000310.GIF" wi="157" he="55" />全局漏检概率<img file="FDA00008576527500000311.GIF" wi="159" he="55" />和全局虚警概率<img file="FDA00008576527500000312.GIF" wi="147" he="55" />的计算公式分别如下:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>=</mo><mn>1</mn><mo>-</mo><mi>Q</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&gamma;</mi><mrow><mi>o</mi><mi>p</mi><mi>t</mi></mrow></msub><mo>-</mo><mn>1</mn></mrow><msqrt><mrow><mo>(</mo><mn>2</mn><mo>/</mo><msup><mi>N</mi><mo>&prime;</mo></msup><mo>)</mo></mrow></msqrt></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000032.GIF" wi="533" he="222" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>=</mo><mroot><mrow><munderover><mo>&Pi;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>m</mi></mrow></munderover><msub><mi>&kappa;</mi><mrow><mn>2</mn><mo>,</mo><mi>t</mi></mrow></msub></mrow><mrow><mi>N</mi><mo>-</mo><mi>m</mi></mrow></mroot><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mi>N</mi><mo>-</mo><mi>m</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>)</mo></mrow><mi>l</mi></msup><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>l</mi></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000033.GIF" wi="1013" he="158" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>D</mi><mrow><mi>F</mi><mi>a</mi><mi>i</mi><mi>l</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>=</mo><mn>1</mn><mo>-</mo><msub><mi>D</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>;</mo></mrow>]]></math><img file="FDA00008576527500000313.GIF" wi="405" he="62" /></maths>(5)频谱感知融合中心FC根据授权用户PU频谱为占用状态H<sub>1</sub>对应的全局漏检概率<img file="FDA00008576527500000314.GIF" wi="157" he="62" />以及授权用户PU频谱为空闲状态H<sub>0</sub>对应的全局虚警概率<img file="FDA00008576527500000315.GIF" wi="173" he="62" />建立基于认知用户数目的频谱感知误差函数Fun(m);其中,频谱感知误差函数Fun(m)计算公式如下:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>F</mi><mi>u</mi><mi>n</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>P</mi><mrow><mi>p</mi><mi>u</mi></mrow></msub><mo>&CenterDot;</mo><msub><mi>D</mi><mrow><mi>u</mi><mi>n</mi><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>+</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>p</mi><mi>u</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>D</mi><mrow><mi>F</mi><mi>a</mi><mi>i</mi><mi>l</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>P</mi><mrow><mi>p</mi><mi>u</mi></mrow></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>D</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>p</mi><mi>u</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>D</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msub><mi>P</mi><mrow><mi>p</mi><mi>u</mi></mrow></msub><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mroot><mrow><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&kappa;</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub></mrow><mi>m</mi></mroot><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mi>m</mi></mrow><mi>N</mi></munderover><mrow><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub><mo>)</mo></mrow><mi>l</mi></msup><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></msub></mrow><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>l</mi></mrow></msup></mrow></mrow><mo>)</mo></mrow><mo>+</mo></mrow></mtd></mtr><mtr><mtd><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>p</mi><mi>u</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mroot><mrow><munderover><mo>&Pi;</mo><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mi>m</mi></mrow></munderover><msub><mi>&kappa;</mi><mrow><mn>2</mn><mo>,</mo><mi>t</mi></mrow></msub></mrow><mrow><mi>N</mi><mo>-</mo><mi>m</mi></mrow></mroot><mo>&CenterDot;</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mi>N</mi><mo>-</mo><mi>m</mi><mo>+</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub><mo>)</mo></mrow><mi>l</mi></msup><msup><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>det</mi><mo>,</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></msub></mrow><mo>)</mo></mrow><mrow><mi>N</mi><mo>-</mo><mi>l</mi></mrow></msup></mrow></mrow><mo>)</mo></mrow></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000034.GIF" wi="1285" he="477" /></maths>其中,P<sub>pu</sub>表示授权用户PU信号在其授权频谱出现的概率;(6)计算频谱感知误差函数Fun(m)的频谱感知误差最小值Fun(m<sub>0</sub>),并以该频谱感知误差函数最小值Fun(m<sub>0</sub>)对应的数值m<sub>0</sub>(m<sub>0</sub>£m)作为参与协同感知的最佳协同认知用户数目,并对m个认知用户按照其对应的信噪比值snr<sub>i</sub>进行降序排列,得到m个认知用户的降序排列组;(7)选取认知用户降序排列组中的前m<sub>0</sub>个认知用户作为参与协同感知的最佳协同认知用户;其中,分别标记选取的最佳协同认知用户为CR'<sub>r</sub>,其中,r=1,2,…,m<sub>0</sub>;(8)根据步骤(7)中m<sub>0</sub>个最佳协同认知用户CR'<sub>r</sub>的检测概率,频谱感知融合中心FC以加权的OR准则协同后的全局检测概率作为N个认知用户协同感知的最终检测结果;其中,加权的OR准则如下:<maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>d</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mn>0</mn></msub></munderover><msub><mi>&omega;</mi><mi>r</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>Q</mi><mrow><mi>f</mi><mi>a</mi></mrow></msub><mo>=</mo><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mn>0</mn></msub></munderover><msub><mi>&omega;</mi><mi>r</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>&omega;</mi><mi>r</mi></msub><mo>=</mo><mfrac><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>r</mi></mrow></msub><mrow><munderover><mo>&Sigma;</mo><mrow><mi>r</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>m</mi><mn>0</mn></msub></munderover><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>r</mi></mrow></msub></mrow></mfrac><mo>,</mo><mi>r</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>m</mi><mn>0</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000857652750000041.GIF" wi="1652" he="198" /></maths>其中,P<sub>d,r</sub>为最佳协同认知用户CR'<sub>r</sub>的检测概率,P<sub>fa,r</sub>为最佳协同认知用户CR'<sub>r</sub>的虚警概率;Q<sub>d</sub>为协同感知后的全局检测概率,Q<sub>fa</sub>为协同感知后的全局虚警概率;m<sub>0</sub>为最佳协同认知用户的数目;ω<sub>r</sub>为最佳协同认知用户CR'<sub>r</sub>的加权系数。
地址 315211 浙江省宁波市江北区风华路818号