发明名称 一种认知无线电协作频谱感知方法
摘要 本发明涉及一种认知无线电协作频谱感知方法,该方法包括各认知用户分别独立地发送自身信噪比和检测结果给频谱感知融合中心,由频谱感知融合中心根据预设信噪比筛选值选出参与协作的初选认知用户,计算初选认知用户的信任度和所有初选认知用户的信噪比均方根值,得到各初选认知用户信噪比和信噪比均方根值间的商值,计算信噪比预设阈值、信噪比最佳阈值和初选认知用户的联合筛选参数值,由频谱感知融合中心根据联系筛选参数值,选定参与协作的复选认知用户,并得到复选认知用户调整后的检测概率,最后将重新选择的终选认知用户按照加权的OR准则进行协作检测。该协作频谱感知方法减小了低信噪比认知用户对协作检测影响,提高了协作检测性能。
申请公布号 CN105141384A 申请公布日期 2015.12.09
申请号 CN201510595118.8 申请日期 2015.09.17
申请人 宁波大学 发明人 郑紫微;胡峰
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 宁波诚源专利事务所有限公司 33102 代理人 刘凤钦
主权项 一种认知无线电协作频谱感知方法,其特征在于,依次包括如下步骤:(1)N个认知用户CR<sub>i</sub>分别独立地获取自身的信噪比SNR<sub>i</sub>和频谱检测结果,并分别将获取的信噪比SNR<sub>i</sub>和频谱检测结果发送至频谱感知融合中心FC,其中,频谱检测结果包括认知用户CR<sub>i</sub>的检测概率P<sub>d,i</sub>以及虚警概率P<sub>f,i</sub>,i=1,2,…,N;(2)频谱感知融合中心FC接收各认知用户CR<sub>i</sub>发送来的信噪比SNR<sub>i</sub>和频谱检测结果,并判断信噪比SNR<sub>i</sub>大于预设的信噪比筛选值SNR<sub>chose</sub>时,则选择此信噪比对应的认知用户为参与协作检测的初选认知用户组成员,并执行步骤(3);否则,选择具有最高信噪比的认知用户所对应的频谱检测结果为频谱感知融合中心FC的最终检测结果;(3)根据参与协作检测的初选认知用户组成员的信噪比和检测概率,计算各初选认知用户的信任度y<sub>i</sub>以及初选认知用户组的信噪比均方根值<img file="FDA0000805053190000011.GIF" wi="66" he="78" />并令信噪比SNR'<sub>i</sub>=γ<sub>i</sub>,初选认知用户组成员的数量为N',且N'≤N,其中,初选认知用户的信任度Ψ<sub>i</sub>、信噪比均方根值<img file="FDA0000805053190000012.GIF" wi="39" he="81" />的计算公式如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&psi;</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>i</mi></mrow></msub><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mo>&prime;</mo></msup></munderover><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><msup><mi>N</mi><mo>&prime;</mo></msup></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>N</mi><mo>&prime;</mo></msup></munderover><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mn>2</mn></mrow></msup><mi>i</mi></msub></mrow></msqrt><mo>,</mo></mrow></mtd><mtd><mrow><msup><mi>N</mi><mo>&prime;</mo></msup><mo>&le;</mo><mi>N</mi></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000013.GIF" wi="925" he="247" /></maths>(4)频谱感知融合中心FC分别依次计算步骤(3)中各初选认知用户信噪比SNR'<sub>i</sub>与所有初选认知用户信噪比均方根值<img file="FDA0000805053190000014.GIF" wi="40" he="78" />之间的商值η<sub>i</sub>,其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mtable><mtr><mtd><mrow><msub><mi>&eta;</mi><mi>i</mi></msub><mo>=</mo><mo>|</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>/</mo><mover><mi>&gamma;</mi><mo>&OverBar;</mo></mover><mo>|</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi></mrow></mtd></mtr></mtable><mo>&prime;</mo></msup><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000015.GIF" wi="537" he="149" /></maths>(5)频谱感知融合中心FC计算、获取信噪比预设阈值λ、信噪比最佳阈值λ<sub>optimal</sub>,以及初选认知用户CR'<sub>i</sub>的信任度Ψ<sub>i</sub>,计算各初选认知用户CR'<sub>i</sub>的联合筛选参数值ξ<sub>i</sub>,并选定参与协作的复选认知用户:(a)频谱感知融合中心FC根据接收的N'个初选认知用户CR'<sub>i</sub>对应的信噪比集合{SNR'<sub>i</sub>},获取初选认知用户信噪比集合{SNR'<sub>i</sub>}中的信噪比最大值,其中,记该信噪比最大值为SNR'max;(b)以获取的信噪比最大值SNR'max为参考,并将信噪比最大值SNR'max分别与N'个初选认知用户CR'<sub>i</sub>的信噪比SNR'<sub>i</sub>作商处理,计算得到各初选认知用户信噪比SNR'<sub>i</sub>所对应的初始阈值λ<sub>i</sub>,其中,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&lambda;</mi><mi>i</mi></msub><mo>=</mo><mo>|</mo><msub><msup><mi>SNR</mi><mo>&prime;</mo></msup><mi>max</mi></msub><mo>/</mo><msub><msup><mi>SNR</mi><mo>&prime;</mo></msup><mi>i</mi></msub><mo>|</mo><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>N</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000021.GIF" wi="721" he="129" /></maths>(c)根据各初选认知用户CR'<sub>i</sub>的信任度Ψ<sub>i</sub>和信噪比商值η<sub>i</sub>,计算各初选认知用户CR'<sub>i</sub>的联合筛选参数值ξ<sub>i</sub>,并根据联合筛选参数值ξ<sub>i</sub>,选取参与协作的复选认知用户CR”<sub>j</sub>,其中,复选认知用户CR”<sub>j</sub>的数量为M,ξ<sub>i</sub>=η<sub>i</sub>·Ψ<sub>i</sub>,i=1,2,…,N',j=1,2,…,M,M≤N':若联合筛选参数值ξ<sub>i</sub>位于预设数值区间范围[ξ<sub>a</sub>,ξ<sub>b</sub>]内,即ξ<sub>a</sub>≤ξ<sub>i</sub>≤ξ<sub>b</sub>时,则选取该联合筛选参数值ξ<sub>i</sub>对应的初选认知用户为复选认知用户,并参与协作检测;否则,该初选认知用户不予选取;(d)根据步骤(c)中的信噪比预设阈值λ,获取M个复选认知用户CR”<sub>j</sub>分别在OR准则和AND准则下的协作检测性能曲线,其中,OR准则:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>Q</mi><mi>d</mi></msub><mo>=</mo><mn>1</mn><mo>-</mo><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&omega;</mi><mi>j</mi></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><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>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>&omega;</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>j</mi></msub></mrow><mrow><mn>0.5</mn><mo>&CenterDot;</mo><mrow><mo>(</mo><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>max</mi></msub><mo>+</mo><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>min</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000023.GIF" wi="1703" he="159" /></maths>AND准则:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>Q</mi><mi>d</mi></msub><mo>=</mo><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><msub><mi>Q</mi><mrow><mi>f</mi><mi>a</mi></mrow></msub><mo>=</mo><munderover><mo>&Pi;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi><mo>,</mo><mi>j</mi></mrow></msub></mrow></mtd></mtr></mtable><mo>,</mo></mrow>]]></math><img file="FDA0000805053190000022.GIF" wi="726" he="148" /></maths>j=1,2,…,M,M≤N';其中,P<sub>d,j</sub>为第j个复选认知用户CR”<sub>j</sub>的检测概率,P<sub>fa,j</sub>为第j个复选认知用户CR”<sub>j</sub>的虚警概率;Q<sub>d</sub>为协作检测后的全局检测概率,Q<sub>fa</sub>为协作检测后的全局虚警概率;ω<sub>j</sub>表示信噪比SNR”<sub>j</sub>的权重系数,SNR”<sub>j</sub>是第j个复选认知用户CR”<sub>j</sub>的信噪比,SNR”<sub>max</sub>表示M个复选认知用户的信噪比最大值,SNR”<sub>min</sub>表示M个复选认知用户的信噪比最小值;(e)根据OR准则和AND准则下的协作检测性能曲线,得到在OR准则和AND准则下的最佳检测性能值Q<sub>d‑max</sub>,并以该最佳检测性能值Q<sub>d‑max</sub>所对应的信噪比预设阈值为信噪比最佳阈值,其中记信噪比最佳阈值为λ<sub>optimal</sub>,最佳检测性能值即为全局检测概率Q<sub>d</sub>中的最大值;(6)根据获取的信噪比最佳阈值λ<sub>optimal</sub>,得到该信噪比最佳阈值λ<sub>optimal</sub>对应的复选认知用户CR”,获取该复选认知用户CR”的调整因子α以及其他M‑1个复选认知用户CR”<sub>k</sub>的调整因子α<sub>k</sub>,并分别根据调整因子α、α<sub>k</sub>对应调整复选认知用户CR”、CR”<sub>k</sub>的虚警概率,其中,复选认知用户CR”调整后的虚警概率记为P<sub>fa</sub>,复选认知用户CR”<sub>k</sub>调整后的虚警概率记为P<sub>fa,k</sub>;P<sub>fa,k</sub>=α<sub>k</sub>·P<sub>fa</sub>,k=1,2,…,M‑1;<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>=</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mover><mrow><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></mrow><mo>&OverBar;</mo></mover><mo>-</mo><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>k</mi></msub></mrow><mover><mrow><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></mrow><mo>&OverBar;</mo></mover></mfrac><mo>,</mo></mrow></mtd><mtd><mrow><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi><mo>-</mo><mn>1</mn></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000031.GIF" wi="903" he="159" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mover><mrow><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><msqrt><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo><mn>2</mn></mrow></msup><mi>j</mi></msub></mrow></msqrt><mo>,</mo></mrow></mtd><mtd><mrow><mi>M</mi><mo>&le;</mo><msup><mi>N</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000032.GIF" wi="759" he="169" /></maths>其中,α<sub>k</sub>为复选认知用户CR”<sub>k</sub>的调整因子,用来根据复选认知用户CR”<sub>k</sub>自身的信噪比实现对其虚警概率大小的调整;SNR”<sub>j</sub>为第j个复选认知用户CR”<sub>j</sub>的信噪比;(7)根据步骤(6)中获取的M个复选认知用户的调整因子α<sub>j</sub>以及对应调整后的虚警概率P<sub>fa,j</sub>,计算复选认知用户CR”<sub>j</sub>调整后的判决阈值λ'<sub>i</sub>和检测概率P<sub>d,j</sub>,其中,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><msup><mi>&lambda;</mi><mo>&prime;</mo></msup><mi>j</mi></msub><mo>=</mo><msubsup><mi>&sigma;</mi><mi>w</mi><mn>2</mn></msubsup><mo>&lsqb;</mo><msqrt><mrow><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msqrt><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>N</mi><mn>1</mn></msub><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msubsup><mi>&sigma;</mi><mi>w</mi><mn>2</mn></msubsup><mo>&lsqb;</mo><msqrt><mrow><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msqrt><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mi>&delta;</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>N</mi><mn>1</mn></msub><mo>&rsqb;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><msubsup><mi>&sigma;</mi><mi>w</mi><mn>2</mn></msubsup><mo>&lsqb;</mo><msqrt><mrow><mn>2</mn><msub><mi>N</mi><mn>1</mn></msub></mrow></msqrt><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mover><mrow><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></mrow><mo>&OverBar;</mo></mover><mo>-</mo><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>j</mi></msub></mrow><mover><mrow><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup></mrow><mo>&OverBar;</mo></mover></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi></mrow></msub><mo>)</mo><mo>+</mo><msub><mi>N</mi><mn>1</mn></msub><mo>&rsqb;</mo></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000033.GIF" wi="1028" he="335" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>d</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mi>Q</mi><mo>&lsqb;</mo><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msqrt><mrow><msub><mi>N</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>j</mi></msub></mrow></msqrt><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000034.GIF" wi="731" he="118" /></maths><maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>N</mi><mn>1</mn></msub><mo>=</mo><mn>2</mn><msup><mrow><mo>&lsqb;</mo><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>Q</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi></mrow></msub><mo>)</mo></mrow><msqrt><mrow><mn>1</mn><mo>+</mo><mn>2</mn><msub><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mi>j</mi></msub></mrow></msqrt><mo>&rsqb;</mo></mrow><mn>2</mn></msup><msup><mi>SNR</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><msup><msub><mrow></mrow><mi>j</mi></msub><mrow><mo>-</mo><mn>2</mn></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0000805053190000035.GIF" wi="1064" he="125" /></maths>其中,<img file="FDA0000805053190000036.GIF" wi="528" he="140" />j=1,2,…,M,M≤N';N<sub>1</sub>为采样点数;(8)根据步骤(7)中M个复选认知用户的信噪比SNR”<sub>j</sub>以及得到的调整后的检测概率P<sub>d,j</sub>,返回步骤(5),重新在M个复选认知用户中选择参与协作的终选认知用户,并以加权的OR准则协作后的全局检测概率为频谱感知融合中心FC的最终检测结果。
地址 315211 浙江省宁波市江北区风华路818号