发明名称 一种基于遗传算法的并行合作频谱感知的方法
摘要 本发明属于认知无线电频谱感知技术领域,具体涉及一种基于遗传算法的并行合作频谱感知的方法。本发明包括:初始群体的产生;个体修正;计算适应度函数;选取适应性最强的个体;变异;生存竞争;优胜者交叉;随机产生个体;个体修正;最适应个体直接进入下一代;截止。本发明能够充分考虑不同认知用户对不同异构信道感知性能的差异性、系统有效吞吐量以及感知时长,最优的确定感知策略,即确定感知任务分配和感知时长,相比随机感知方法,能够获得更大的系统有效吞吐量。
申请公布号 CN104243063B 申请公布日期 2017.01.11
申请号 CN201410431708.2 申请日期 2014.08.28
申请人 哈尔滨工程大学 发明人 李一兵;富爽;叶方;林云;黄诚;李一晨;陈杰
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 代理人
主权项 一种基于遗传算法的并行合作频谱感知的方法,其特征在于:(1)初始群体的产生:算法首先产生一个个体数为L的种群,L取4j,其中j为正整数,每一个个体代表一种感知方案,采用整数编码,将感知信道分配矢量X和感知时长等级k连接得到染色体;(2)个体修正:当信道感知的认知用户大于N<sub>co_max</sub>时,则置多余的认知用户对应的基因值为0,即该认知用户不感知信道;(3)计算适应度函数:适应度函数定义为系统有效吞吐量R<sub>sum</sub>(X,k);(4)选取适应性最强的个体:在当前种群中选取适应性最强的个体,作为步骤(5)的父代;(5)变异:由于最适应个体的基因是最优良的,选用最适应个体作为变异的父代,变异产生L/4个变异子代,直接进入下一代,变异概率为P<sub>l</sub>,P<sub>l</sub>=0.1;(6)生存竞争:L个个体两两组成一组,进行生存竞争,得到L/2个优胜个体;(7)优胜者交叉:由生存竞争产生的L/2个优胜个体以交叉概率P<sub>c</sub>两两进行双点交叉后,产生L/2个子代,直接进入下一代,P<sub>c</sub>=0.5;(8)随机产生个体:为了避免遗传算法早熟的问题,每次进化时产生L/4个个体,加入下一代;(9)个体修正:对不满足要求的个体进行修正;(10)最适应个体直接进入下一代:计算当前种群的适应度,用步骤(4)中选取的最优个体,替换当前群体中的最差个体,直接进入下一代,保证种群能够稳定地进化;(11)截止:当两个截止条件满足其一时,进化停止,得到最优个体,获得最优感知任务分配及感知时长,截止条件一是达到最大进化代数W,二是最优个体和最差个体之间适应度差异足够小,即小于最优个体适应度的V倍的时候,如果不满足截止条件,则跳到步骤(4)进行下一次进化;所述的系统模型为:在一个由一个融合中心和N个认知用户构成的集中式认知无线电网络中,感知对象包括M个异构信道;每一个子信道具有不同的带宽,不同的传输速率,不同的使用概率;分别用H<sub>0</sub>和H<sub>1</sub>指示信道空闲和被主用户占用两种状态;认知用户在每一个时隙T内分为感知时隙T<sub>s</sub>和传输时隙T<sub>r</sub>两部分,T=T<sub>s</sub>+T<sub>r</sub>;在每个时隙T内,T<sub>s</sub>和T<sub>r</sub>可在一定范围内变化;为便于各节点同步及减少传输开销,T<sub>s</sub>只能在K个离散值之中选取,各离散值由小到大排列,记为T<sub>s</sub>(k),其中k=1,2,...,K;每个认知用户在一个感知时隙T<sub>s</sub>(k)内仅感知一个子信道,而每个信道可由多个认知用户合作感知以提高准确性,合作时采用OR准则;设各节点使用能量检测方法进行频谱感知,信号为复PSK信号,噪声为循环对称复高斯噪声;由于每个用户所处的位置不同,不同认知用户处不同信道的信噪比也不同,设已知各认知用户处各信道的信噪比Υ=[γ<sub>mn</sub>]<sub>M×N</sub>,其中γ<sub>mn</sub>表示第n个认知用户处第m个信道的信噪比;每个认知用户的虚警概率<img file="FDA0001003505050000021.GIF" wi="85" he="70" />保持一致,设<img file="FDA0001003505050000022.GIF" wi="203" he="72" />其中n=1,2,…,N,m=1,2,…,M;则第n个认知用户对第m个信道的检测概率P<sub>d</sub><sup>mn</sup>(T<sub>s</sub>(k))可由下式得出<maths num="0001"><math><![CDATA[<mrow><msup><msub><mi>P</mi><mi>d</mi></msub><mrow><mi>m</mi><mi>n</mi></mrow></msup><mrow><mo>(</mo><msub><mi>T</mi><mi>s</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>e</mi><mi>r</mi><mi>f</mi><mi>c</mi><mo>{</mo><mfrac><mn>1</mn><msqrt><mrow><mn>2</mn><mrow><mo>(</mo><mn>2</mn><msub><mi>&gamma;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msqrt></mfrac><mo>&lsqb;</mo><msqrt><mn>2</mn></msqrt><msup><mi>erf</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mn>2</mn><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>f</mi></msub><mo>)</mo></mrow><mo>-</mo><msqrt><mrow><msub><mi>T</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>f</mi><mi>s</mi></msub></mrow></msqrt><msub><mi>&gamma;</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>&rsqb;</mo><mo>}</mo></mrow>]]></math><img file="FDA0001003505050000023.GIF" wi="1370" he="166" /></maths>其中,erfc为互补误差函数,f<sub>s</sub>为采样频率;设X=[x<sub>n</sub>]<sub>N</sub>为分配向量,其中,x<sub>n</sub>代表第n个认知用户感知的信道序号,x<sub>n</sub>=0,1,2,...,M,为0时表示该认知用户不感知任何信道;则信道m的检测概率P<sub>D</sub><sup>m</sup>(X,T<sub>s</sub>(k))和虚警概率P<sub>F</sub><sup>m</sup>(X)可表示为<maths num="0002"><math><![CDATA[<mrow><msup><msub><mi>P</mi><mi>D</mi></msub><mi>m</mi></msup><mrow><mo>(</mo><mi>X</mi><mo>,</mo><msub><mi>T</mi><mi>s</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><munder><mo>&Pi;</mo><mrow><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><mi>m</mi></mrow></munder><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msup><msub><mi>P</mi><mi>d</mi></msub><mrow><mi>m</mi><mi>n</mi></mrow></msup><mrow><mo>(</mo><msub><mi>T</mi><mi>s</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001003505050000024.GIF" wi="773" he="111" /></maths><maths num="0003"><math><![CDATA[<mrow><msup><msub><mi>P</mi><mi>F</mi></msub><mi>m</mi></msup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn><mo>-</mo><munder><mo>&Pi;</mo><mrow><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><mi>m</mi></mrow></munder><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>f</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001003505050000025.GIF" wi="485" he="112" /></maths>设<img file="FDA0001003505050000029.GIF" wi="74" he="71" />为信道m在信道空闲时所能达到的信道速率;当主用户不存在时,且认知用户能够正常的感知到空闲的频谱,该信道能够达到的平均吞吐量<img file="FDA00010035050500000210.GIF" wi="75" he="70" />为<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>R</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mi>T</mi></mfrac><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msubsup><mi>P</mi><mi>F</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msubsup><mi>C</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup></mrow>]]></math><img file="FDA0001003505050000026.GIF" wi="746" he="118" /></maths>设信道m空闲的概率为<img file="FDA00010035050500000212.GIF" wi="97" he="70" />则该信道所能达到的平均吞吐量的近似值<img file="FDA00010035050500000211.GIF" wi="186" he="70" />为<maths num="0005"><math><![CDATA[<mrow><msup><mover><mi>R</mi><mo>~</mo></mover><mi>m</mi></msup><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>&ap;</mo><msubsup><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup><msubsup><mi>R</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mi>T</mi></mfrac><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msubsup><mi>P</mi><mi>F</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msubsup><mi>C</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup></mrow>]]></math><img file="FDA0001003505050000027.GIF" wi="1116" he="117" /></maths>定义系统有效吞吐量R<sub>sum</sub>(X,k)为检测概率大于检测概率要求<img file="FDA00010035050500000213.GIF" wi="57" he="67" />的信道的吞吐量的总和,表达式如下:<maths num="0006"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>R</mi><mrow><mi>s</mi><mi>u</mi><mi>m</mi></mrow></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><munder><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo></mrow></munder><mrow><msubsup><mi>P</mi><mi>D</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>,</mo><msub><mi>T</mi><mi>s</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>D</mi></msub></mrow><mi>M</mi></munderover><msub><mover><mi>R</mi><mo>~</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munderover><munder><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo></mrow></munder><mrow><msubsup><mi>P</mi><mi>D</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>,</mo><msub><mi>T</mi><mi>s</mi></msub><mo>(</mo><mi>k</mi><mo>)</mo><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mover><mi>P</mi><mo>&OverBar;</mo></mover><mi>D</mi></msub></mrow><mi>M</mi></munderover><msubsup><mi>P</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mi>T</mi></mfrac><mo>&lsqb;</mo><mn>1</mn><mo>-</mo><msubsup><mi>P</mi><mi>F</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>&rsqb;</mo><msubsup><mi>C</mi><msub><mi>H</mi><mn>0</mn></msub><mi>m</mi></msubsup></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001003505050000028.GIF" wi="1054" he="355" /></maths>此外,为了充分的利用频谱,设系统所允许的最大虚警概率为<img file="FDA0001003505050000032.GIF" wi="83" he="70" />则每个信道允许的最大认知用户数N<sub>co_max</sub>为<img file="FDA0001003505050000031.GIF" wi="558" he="95" />
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室