发明名称 一种认知无线电并行合作频谱感知方法
摘要 一种认知无线电并行合作频谱感知方法,该方法首先根据定义的感知开销和传输收益,计算系统效益矩阵,运用Kuhn-Munkres算法,为每一个认知用户分配一个子信道进行感知,如果存在剩余的认知用户,则计算系统效益增量矩阵,然后迭代使用Kuhn-Munkres算法最优地为子信道继续分配认知用户与之前分配的认知用户一同对子信道进行合作感知。该并行合作频谱感知方法充分考虑感知开销、系统效益以及不同认知用户对不同子信道感知性能的差异,克服了现有并行合作频谱感知方法仅考虑检测准确性和需要较多用户进行感知的不足,能够最优的确定感知策略,有效的提高系统效益,减少参加合作感知的认知用户数。
申请公布号 CN103237310B 申请公布日期 2015.09.09
申请号 CN201310118195.5 申请日期 2013.04.08
申请人 黑龙江八一农垦大学 发明人 富爽;杜红;贾美娟;孟艳君
分类号 H04W16/14(2009.01)I;H04W24/00(2009.01)I 主分类号 H04W16/14(2009.01)I
代理机构 哈尔滨市文洋专利代理事务所(普通合伙) 23210 代理人 孙淑荣
主权项 一种认知无线电并行合作频谱感知方法,包括下列系统模型:在一个由一个融合中心和N个认知用户构成的集中式认知无线电网络中,感知对象为M个子信道;定义感知开销为认知用户因频谱感知任务所减少传输的数据量,设R<sub>cr</sub><sup>n</sup>为认知用户n的传输速率,其中n=1,2,...,N,则认知用户n的感知开销O<sup>n</sup>为O<sup>n</sup>=R<sub>cr</sub><sup>n</sup>T<sub>s</sub>      (1)其中,T<sub>s</sub>为感知时隙;定义传输收益为子信道用于认知用户传输时单位时间内所能传输的平均有效数据量,子信道用于认知用户传输时分为两种情况:第一种情况是信道空闲且被认知用户正确感知时,第二种情况是信道被主用户占用,但被认知用户错误感知而误认为是空闲时,对于第二种情况,可以看成是对传输收益的惩罚函数,设R<sub>ch</sub><sup>m</sup>为子信道m的传输速率,P<sub>d_ch</sub>=[P<sub>d_ch</sub><sup>m</sup>]<sub>M</sub>为各子信道的检测概率矢量,P<sub>f_ch</sub>=[P<sub>f_ch</sub><sup>m</sup>]<sub>M</sub>为各子信道的虚警概率矢量,其中P<sub>d_ch</sub><sup>m</sup>和P<sub>f_ch</sub><sup>m</sup>分别为认知网络对子信道m的检测概率和虚警概率,信道m=1,2,...,M,则信道m的传输收益G<sup>m</sup>为G<sup>m</sup>=R<sub>ch</sub><sup>m</sup>T<sub>r</sub>P<sub>on</sub>[1‑P<sub>f_ch</sub><sup>m</sup>]‑R<sub>ch</sub><sup>m</sup>(1‑P<sub>on</sub>)[1‑P<sub>d_ch</sub><sup>m</sup>]T<sub>r</sub>δ      (2)其中:T<sub>r</sub>为传输时隙,P<sub>on</sub>为子信道空闲的概率,δ为惩罚因子;设X=[x<sup>mn</sup>]<sub>M×N</sub>为分配矩阵,当认知用户n被分配感知子信道m时,x<sup>mn</sup>=1,否则x<sup>mn</sup>=0;则分配X时的总系统效益U<sub>all</sub>(X)可以定义为总传输收益和总感知开销之差<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>U</mi><mi>all</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mi>G</mi><mi>m</mi></msup><mo>-</mo><munder><mi>&Sigma;</mi><mrow><msup><mi>x</mi><mi>mn</mi></msup><mo>=</mo><mn>1</mn></mrow></munder><msup><mi>O</mi><mi>n</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000138696220000011.GIF" wi="1211" he="130" /></maths>根据以上描述的系统模型,并行合作频谱感知的感知策略可以等效成如下分配问题:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><munder><mi>max</mi><mi>x</mi></munder></mtd><mtd><msub><mi>U</mi><mi>all</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>;</mo></mtd></mtr></mtable></mfenced>]]></math><img file="FSB0000138696220000024.GIF" wi="453" he="90" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mi>x</mi><mi>mn</mi></msup><mo>=</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>;</mo></mrow>]]></math><img file="FSB0000138696220000026.GIF" wi="836" he="141" /></maths>                                   (4)<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>&beta;</mi><mo>&le;</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mi>x</mi><mi>mn</mi></msup><mo>&le;</mo><msub><mi>N</mi><mrow><mi>co</mi><mo>_</mo><mi>max</mi></mrow></msub><mo>,</mo><mi>m</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>M</mi><mo>;</mo></mrow>]]></math><img file="FSB0000138696220000027.GIF" wi="889" he="139" /></maths>x<sup>mn</sup>∈{0,1}.其中,β表示是否能在一个感知时隙内检测完所有子信道,如果N≥M,即认知用户数大于子信道数,则能够在一个感知时隙内检测完所有子信道,令β=1,否则,置β=0,N为认知用户数,N<sub>co_max</sub>为每个子信道允许的最大的参加合作感知的认知用户数;多个认知用户共同感知一个子信道时使用OR准则;上述的并行合作频谱感知方法是按照如下步骤得到的:步骤1初始化认知用户数N,子信道数M,系统所允许的最大虚警概率Q<sub>fa</sub>,系统要求的检测概率P<sub>d_ch</sub><sup>req</sup>,各认知用户感知的虚警概率p<sub>f</sub>,传输时隙T<sub>r</sub>,感知时隙T<sub>s</sub>,惩罚因子δ,认知用户传输速率R<sub>cr</sub><sup>n</sup>,n=1,2,...,N,子信道传输速率R<sub>ch</sub><sup>m</sup>,m=1,2,...,M,子信道空闲的概率P<sub>on</sub>,检测概率矩阵P<sub>d</sub>=[P<sub>d</sub><sup>mn</sup>]<sub>M×N</sub>,其中P<sub>d</sub><sup>mn</sup>为认知用户n感知子信道m的检测概率,置分配矩阵X=0,当前迭代最优分配矩阵X<sub>0</sub>=0,各子信道的检测概率矢量P<sub>d_ch</sub>=0以及各子信道的虚警概率矢量P<sub>f_ch</sub>=0,利用以下公式计算每个子信道允许的最大的参加合作感知的认知用户数N<sub>co_max</sub>;<img file="FSB0000138696220000025.GIF" wi="1275" he="108" />步骤2根据以下公式,计算信道m分配给认知用户n感知时所能获得的传输收益G<sup>mn</sup>、感知开销O<sup>n</sup>以及系统效益U<sup>mn</sup>,其中m=1,2,...,M,n=1,2,...,N;执行KM算法得到当前迭代最优分配矩阵X<sub>0</sub>=[X<sub>0</sub><sup>mn</sup>]<sub>M×N</sub>使总系统效益U<sub>all</sub>最大,其中如果子信道m分配给认知用户n感知,则X<sub>0</sub><sup>mn</sup>=1,否则为0,更新X=X<sub>0</sub>;G<sup>mn</sup>=R<sub>ch</sub><sup>m</sup>T<sub>r</sub>P<sub>on</sub>[1‑p<sub>f</sub>]‑R<sub>ch</sub><sup>m</sup>(1‑P<sub>on</sub>)[1‑P<sub>d</sub><sup>mn</sup>]T<sub>r</sub>δ     (6)O<sup>n</sup>=R<sub>cr</sub><sup>n</sup>T<sub>s</sub>       (1)U<sup>mn</sup>=G<sup>mn</sup>‑O<sup>n</sup>       (7)步骤3根据以下公式计算分配X下的总系统效益U<sub>all</sub>(X);根据当前迭代最优分配矩阵X<sub>0</sub>,对于X<sub>0</sub><sup>mn</sup>=1的指派,置P<sub>d_ch</sub><sup>m</sup>=P<sub>d</sub><sup>mn</sup>,P<sub>f_ch</sub><sup>m</sup>=p<sub>f</sub>;如果N≤M,则算法终止,得到最优分配矩阵X;<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>U</mi><mi>all</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><msup><mi>x</mi><mi>mn</mi></msup><mo>=</mo><mn>1</mn></mrow></munder><msup><mi>U</mi><mi>mn</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSB0000138696220000031.GIF" wi="1246" he="102" /></maths>步骤4对于每一个未分配子信道的认知用户,以及每一个分配认知用户数小于N<sub>co_max</sub>的子信道,根据以下公式计算检测概率增量ΔP<sub>d_ch</sub><sup>mn</sup>、虚警概率增量ΔP<sub>f_ch</sub><sup>m</sup>、传输收益的增量ΔG<sup>mn</sup>、感知开销的增量ΔO<sup>n</sup>以及系统效益的增量ΔU<sup>mn</sup>;ΔP<sub>d_ch</sub><sup>mn</sup>=1‑[1‑P<sub>d_ch</sub><sup>m</sup>][1‑P<sub>d</sub><sup>mn</sup>]‑P<sub>d_ch</sub><sup>m</sup>      (9)ΔP<sub>f_ch</sub><sup>m</sup>=1‑[1‑P<sub>f_ch</sub><sup>m</sup>][1‑p<sub>f</sub>]‑p<sub>f_ch</sub><sup>m</sup>     (10)ΔG<sup>mn</sup>=R<sub>ch</sub><sup>m</sup>(1‑P<sub>on</sub>)ΔP<sub>d_ch</sub><sup>mn</sup>T<sub>r</sub>‑R<sub>ch</sub><sup>m</sup>T<sub>r</sub>P<sub>on</sub>ΔP<sub>f_ch</sub><sup>m</sup>      (11)ΔO<sup>n</sup>=R<sub>cr</sub><sup>n</sup>T<sub>s</sub>       (12)ΔU<sup>mn</sup>=ΔG<sup>mn</sup>‑ΔO<sup>n</sup>        (13)步骤5在所有未分配子信道的认知用户和分配认知用户数小于N<sub>co_max</sub>的子信道中,执行KM算法,得到使系统效益总增量<img file="FSB0000138696220000032.GIF" wi="378" he="110" />最大的当前迭代最优分配X<sub>0</sub>;对于任意X<sub>0</sub><sup>mn</sup>=1的分配,如果ΔU<sup>mn</sup>=0,或者P<sub>d_ch</sub><sup>m</sup>≥P<sub>d_ch</sub><sup>req</sup>且ΔU<sup>mn</sup>≤0,则置X<sub>0</sub><sup>mn</sup>=0;步骤6根据以下公式,更新P<sub>d_ch</sub><sup>m</sup>,P<sub>f_ch</sub><sup>m</sup>,X以及U<sub>all</sub>;P<sub>d_ch</sub><sup>m</sup>=P<sub>d_ch</sub><sup>m</sup>+ΔP<sub>d_ch</sub><sup>mn</sup>       (14)P<sub>f_ch</sub><sup>m</sup>=P<sub>f_ch</sub><sup>m</sup>+ΔP<sub>f_ch</sub><sup>m</sup>       (15)X=X+X<sub>0</sub>         (16)U<sub>all</sub>=U<sub>all</sub>+ΔU<sub>all</sub>         (17)步骤7当以下两个终止条件满足其一时,迭代终止,得到分配矩阵X;一是所有的认知用户都被分配了子信道进行感知或每一个子信道都被分配了N<sub>co_max</sub>个认知用户进行感知;二是,所有子信道m=1,2,...,M均满足P<sub>d_ch</sub><sup>m</sup>≥P<sub>d_ch</sub><sup>req</sup>,且对于所有m=1,2,...,M,n=1,2,...,N,都有ΔU<sup>mn</sup>≤0,否则,跳到步骤4继续执行迭代。
地址 163319 黑龙江省大庆市开发区黑龙江八一农垦大学信息技术学院