发明名称 一种基于量子蛙跳的频谱感知算法
摘要 本发明涉及一种基于量子蛙跳的频谱感知算法,包括:根据目标函数所表示的适应度值确定全局最优的青蛙量子位置以及每个族群内最好、最差适应度值的青蛙量子位置;根据跳跃性方程更新最差适应度值的青蛙量子位置的相关信息以及向全局最优解收敛。本发明通过先将青蛙种群划分为若干个族群;根据目标函数所表示的适应度值来确定全局最优的青蛙量子位置以及每个族群内最好、最差适应度值的青蛙量子位置;再按照提出的跳跃性方程更新其中最差个体的相关信息,快速向全局最优解收敛。
申请公布号 CN104467999B 申请公布日期 2017.02.22
申请号 CN201410659258.2 申请日期 2014.11.18
申请人 北京邮电大学 发明人 张勇;宋梅;魏翼飞;滕颖蕾;郭达;王莉;满毅;成晨;张亚男;鲍叙言;程刚;王东安;李沸乐;方一鸣
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 李相雨
主权项 一种基于量子蛙跳的频谱感知算法,其特征在于,包括以下具体步骤:S1:根据目标函数所表示的适应度值确定全局最优的青蛙量子位置以及每个族群内最好、最差适应度值的青蛙量子位置;S2:根据跳跃性方程更新最差适应度值的青蛙量子位置的相关信息以及向所述全局最优解收敛;所述步骤S2具体包括:S21:根据所述青蛙族群中的最优解进行更新;S22:若没有改进,则根据全局最优解进行更新,其中,每个青蛙的量子位置表示方程的一个解集,每个量子表示该解集中的一个解,每个解代表每个认知用户的权重向量;所述步骤S2还包括:具体地,在QSFLA中,第i只青蛙的量子位置表示为<maths num="0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>d</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>d</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>d</mi><mrow><mi>i</mi><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>x</mi><mrow><mi>i</mi><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&beta;</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>&beta;</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>&beta;</mi><mrow><mi>i</mi><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0001110132420000011.GIF" wi="961" he="150" /></maths>其中,|x<sub>ij</sub>|<sup>2</sup>+|β<sub>ij</sub>|<sup>2</sup>=1,(j=1,2,…,M);定义x<sub>ij</sub>和β<sub>ij</sub>为实数并且0≤x<sub>ij</sub>≤1,0≤β<sub>ij</sub>≤1。量子旋转角为θ<sub>ij</sub>,第i只青蛙的量子位位置d<sub>ij</sub>运用量子旋转门U(θ<sub>ij</sub>)更新;第i只青蛙的量子位置的第j个量子位位置按照如下方式更新<maths num="0002"><math><![CDATA[<mrow><msubsup><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mi>n</mi><mi>e</mi><mi>w</mi></mrow></msubsup><mo>=</mo><mi>a</mi><mi>b</mi><mi>s</mi><mrow><mo>(</mo><mi>U</mi><mo>(</mo><msub><mi>&theta;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>a</mi><mi>b</mi><mi>s</mi><mo>(</mo><mrow><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>cos&theta;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>sin&theta;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>sin&theta;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mtd><mtd><mrow><msub><mi>cos&theta;</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow><mo>)</mo></mrow>]]></math><img file="FDA0001110132420000012.GIF" wi="917" he="159" /></maths>其中,<img file="FDA0001110132420000013.GIF" wi="84" he="71" />为第i只青蛙的量子位置的第j个量子位位置;所述步骤S2还包括:通过公式<maths num="0003"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&theta;</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>=</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>b</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>z</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow><mrow><mi>n</mi><mi>e</mi><mi>w</mi></mrow></msubsup><mo>=</mo><mo>|</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><msub><mi>cos&theta;</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>-</mo><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><msub><mi>sin&theta;</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>|</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001110132420000014.GIF" wi="766" he="214" /></maths>更新具有最差适应度值的青蛙的量子位置x<sub>w</sub>;其中,θ<sub>wj</sub>:每个族群中适应度最差的青蛙w的量子位置的第j个量子位的自旋角;x<sub>w</sub>j:每个族群中适应度最差的青蛙w的量子位置的第j个量子位的位置;x<sub>bj</sub>:每个族群中适应度最优的青蛙b的量子位置的第j个量子位的位置;<img file="FDA0001110132420000021.GIF" wi="102" he="70" />每族群中适应度最差的青蛙w量子位置的第j个量子位的更新后的位置;其中,rand为[0,1]之间的均匀随机数,randn为满足均值为0,方差为1的高斯分布的随机数;z=[z<sub>1</sub>,z<sub>2</sub>,...,z<sub>M</sub>]为族群内全部青蛙的量子位置的平均值,其中<img file="FDA0001110132420000022.GIF" wi="600" he="123" />其中,n为每组的青蛙个数;如果新的量子位置<img file="FDA0001110132420000023.GIF" wi="89" he="62" />的适应度值没有改进,则按照下述方式进行更新:<maths num="0004"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>&theta;</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>=</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>g</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>+</mo><mi>r</mi><mi>a</mi><mi>n</mi><mi>d</mi><mi>n</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>z</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow><mrow><mi>n</mi><mi>e</mi><mi>w</mi></mrow></msubsup><mo>=</mo><mo>|</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><msub><mi>cos&theta;</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>-</mo><msqrt><mrow><mn>1</mn><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><msub><mi>sin&theta;</mi><mrow><mi>w</mi><mi>j</mi></mrow></msub><mo>|</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001110132420000024.GIF" wi="774" he="213" /></maths>其中,x<sub>gj</sub>:为全局最优解g的量子位置的第j个量子位的位置。
地址 100876 北京市海淀区北京邮电大学