发明名称 用于绿色认知无线电安全能效最佳功率的分配方法
摘要 本发明公开了一种用于绿色认知无线电安全能效最佳功率的分配方法,其主要解决了现有技术无法保证安全能效最大化问题。其实现步骤是:设置功率分配参数;初始化功率分配参数;计算最佳发送功率的真实值;计算的拉格朗日乘子;计算最佳发送功率的迭代真实值;计算拉格朗日乘子迭代误差;判断拉格朗日乘子是否满足终止条件;计算安全能效函数;计算次级用户安全能效;判断是否满足终止条件;得到次级用户安全能效最佳发送功率和次级用户最大安全能效。本发明具有可应用于绿色认知无线电网络中实现安全能效最大化的功率分配,所需迭代次数少,适用于各种衰落信道,易于实现的优点。
申请公布号 CN105307181A 申请公布日期 2016.02.03
申请号 CN201510640929.5 申请日期 2015.09.30
申请人 西安电子科技大学 发明人 李赞;周福辉;黄海燕;周力存;郝本建;司江勃;刘向丽;关磊;齐佩汉;熊天意;胡伟龙
分类号 H04W16/14(2009.01)I;H04W72/04(2009.01)I;H04W72/08(2009.01)I 主分类号 H04W16/14(2009.01)I
代理机构 陕西电子工业专利中心 61205 代理人 田文英;王品华
主权项 一种用于绿色认知无线电安全能效最佳功率的分配方法,包括如下步骤:(1)设置功率分配参数:(1a)设置容错误差ζ,其取值范围为(0,1);(1b)设置次级用户安全能效的最大迭代次数N,其取值为所选用计算设备的容许计数范围内大于10的正整数;(2)初始化功率分配参数:(2a)将次级用户安全能效的迭代次数初始化为1;(2b)将subgradient迭代算法的迭代次数初始化为1;(3)计算次级用户安全能效最佳发送功率的真实值:(3a)按照下式,计算功率比值系数:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>a</mi><mn>0</mn></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>(</mo><msup><mi>&mu;</mi><mn>0</mn></msup><msub><mi>g</mi><mrow><mi>s</mi><mi>p</mi></mrow></msub><mo>+</mo><msup><mi>&lambda;</mi><mn>0</mn></msup><mo>)</mo><mi>l</mi><mi>n</mi><mn>2</mn></mrow></mtd><mtd><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>(</mo><msup><mi>&eta;</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>&beta;</mi><mo>+</mo><msup><mi>&mu;</mi><mn>0</mn></msup><msub><mi>g</mi><mrow><mi>s</mi><mi>p</mi></mrow></msub><mo>+</mo><msup><mi>&lambda;</mi><mn>0</mn></msup><mo>)</mo><mi>l</mi><mi>n</mi><mn>2</mn></mrow></mtd><mtd><mrow><mi>n</mi><mo>&GreaterEqual;</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000815702300000011.GIF" wi="773" he="199" /></maths>其中,a<sub>0</sub>表示初始化的功率比值系数,μ<sup>0</sup>表示初始化的平均干扰功率约束对应的拉格朗日乘子,其值为0.1,λ<sup>0</sup>表示初始化的平均发送功率约束对应的拉格朗日乘子,其值为0.1,g<sub>sp</sub>表示次级用户发送端到主用户接收端信道功率增益,η<sup>n‑1</sup>表示次级用户第n‑1次迭代时次级用户安全能效,n表示次级用户安全能效的迭代次数,n=1,2,3,...,N,N表示次级用户安全能效的最大迭代次数,β表示次级用户发送机的功率放大因子,ln2表示2的自然对数;(3b)按照下式,计算次级用户发送端到接收端的功率增益比例系数:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>g</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msubsup><mi>&sigma;</mi><mi>r</mi><mn>2</mn></msubsup></mfrac></mrow>]]></math><img file="FDA0000815702300000012.GIF" wi="236" he="159" /></maths>其中,γ<sub>sr</sub>表示次级用户发送端s到接收端r的功率增益比例系数,g<sub>sr</sub>表示次级用户发送端s到接收端r的信道功率增益,<img file="FDA0000815702300000013.GIF" wi="65" he="70" />表示次级用户接收端r受到主用户干扰和加性噪声之和的噪声方差;(3c)按照下式,计算次级用户发送端到监听用户接收端的功率增益比例系数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mo>=</mo><mfrac><msub><mi>g</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><msubsup><mi>&sigma;</mi><mi>e</mi><mn>2</mn></msubsup></mfrac></mrow>]]></math><img file="FDA0000815702300000021.GIF" wi="212" he="142" /></maths>其中,γ<sub>se</sub>表示次级用户发送端s到监听用户接收端e的功率增益比例系数,g<sub>se</sub>表示次级用户发送端s到监听用户接收端e的信道功率增益,<img file="FDA0000815702300000022.GIF" wi="69" he="70" />表示监听用户接收端e受到主用户干扰和加性噪声之和的噪声方差;(3d)按照下式,计算次级用户安全能效最佳发送功率的模糊值:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>F</mi><mn>0</mn><mi>n</mi></msubsup><mo>=</mo><mfrac><mrow><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>a</mi><mn>0</mn></msub><mo>(</mo><mrow><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msub><mi>a</mi><mn>0</mn></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>-</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mo>-</mo><msub><mi>a</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></msqrt><mo>-</mo><msub><mi>a</mi><mn>0</mn></msub><mrow><mo>(</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msub><mi>a</mi><mn>0</mn></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000815702300000023.GIF" wi="1268" he="183" /></maths>其中,<img file="FDA0000815702300000024.GIF" wi="69" he="69" />表示第n次迭代的次级用户安全能效最佳发送功率的模糊值,n表示次级用户安全能效的迭代次数,n=1,2,3,...,N,N表示次级用户安全能效的最大迭代次数,a<sub>0</sub>表示初始化的功率比值系数,γ<sub>sr</sub>表示次级用户发送端s到接收端r的功率增益比例系数,γ<sub>se</sub>表示次级用户发送端s到监听用户接收端e的功率增益比例系数,<img file="FDA0000815702300000025.GIF" wi="77" he="76" />表示平方根操作;(3e)按照下式,计算次级用户安全能效最佳发送功率的真实值:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>P</mi><mn>0</mn><mi>n</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>F</mi><mn>0</mn><mi>n</mi></msubsup><mo>+</mo><mo>|</mo><msubsup><mi>F</mi><mn>0</mn><mi>n</mi></msubsup><mo>|</mo></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000815702300000026.GIF" wi="310" he="150" /></maths>其中,<img file="FDA0000815702300000027.GIF" wi="68" he="71" />表示第n次迭代时次级用户安全能效最佳发送功率的真实值,<img file="FDA0000815702300000028.GIF" wi="69" he="68" />表示第n次迭代时次级用户安全能效最佳发送功率的模糊值,n表示次级用户安全能效的迭代次数,n=1,2,3,...,N,N表示次级用户安全能效的最大迭代次数,|·|表示绝对值操作;(4)计算拉格朗日乘子:利用subgradient迭代算法,计算平均干扰功率约束值对应的拉格朗日乘子μ<sup>k</sup>和平均发送功率约束值对应的拉格朗日乘子λ<sup>k</sup>;(5)计算次级用户安全能效最佳发送功率的迭代真实值:(5a)按照下式,计算功率比值系数:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>a</mi><mi>k</mi></msub><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mo>(</mo><msup><mi>&mu;</mi><mi>k</mi></msup><msub><mi>g</mi><mrow><mi>s</mi><mi>p</mi></mrow></msub><mo>+</mo><msup><mi>&lambda;</mi><mi>k</mi></msup><mo>)</mo><mi>l</mi><mi>n</mi><mn>2</mn></mrow></mtd><mtd><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mo>(</mo><msup><mi>&eta;</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>&beta;</mi><mo>+</mo><msup><mi>&mu;</mi><mi>k</mi></msup><msub><mi>g</mi><mrow><mi>s</mi><mi>p</mi></mrow></msub><mo>+</mo><msup><mi>&lambda;</mi><mi>k</mi></msup><mo>)</mo><mi>l</mi><mi>n</mi><mn>2</mn></mrow></mtd><mtd><mrow><mi>n</mi><mo>&GreaterEqual;</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000815702300000031.GIF" wi="783" he="198" /></maths>其中,a<sub>k</sub>表示第k次迭代时功率比值系数,μ<sup>k</sup>表示使用subgradient迭代算法第k次迭代时平均干扰功率约束值对应的拉格朗日乘子,λ<sup>k</sup>表示使用subgradient迭代算法第k次迭代时平均发送功率约束值对应的拉格朗日乘子,k表示subgradient迭代算法的迭代次数,k=1,2,3,....,K,K表示所选用计算设备的容许计数范围内的最大正整数,g<sub>sp</sub>表示次级用户发送端到主用户接收端信道功率增益,η<sup>n‑1</sup>次级用户第n‑1次迭代时获得的次级用户安全能效,n表示次级用户安全能效的迭代次数,n=1,2,3,...,N,N表示次级用户安全能效的最大迭代次数,β表示次级用户发送机的功率放大因子,ln2表示2的自然对数;(5b)按照下式,计算次级用户安全能效最佳发送功率的迭代模糊值:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>F</mi><mi>k</mi><mi>n</mi></msubsup><mo>=</mo><mfrac><mrow><msqrt><mrow><msubsup><mi>a</mi><mi>k</mi><mn>2</mn></msubsup><msup><mrow><mo>(</mo><mrow><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>4</mn><msub><mi>a</mi><mi>k</mi></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>-</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mo>-</mo><msub><mi>a</mi><mi>k</mi></msub><mo>)</mo></mrow></mrow></msqrt><mo>-</mo><msub><mi>a</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><mn>2</mn><msub><mi>a</mi><mi>k</mi></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000815702300000032.GIF" wi="1234" he="174" /></maths>其中,<img file="FDA0000815702300000033.GIF" wi="70" he="63" />表示第k次迭代的次级用户安全能效最佳发送功率的迭代模糊值,a<sub>k</sub>表示第k次迭代时功率比值系数,k表示subgradient迭代算法的迭代次数,k=1,2,3,....,K,K表示所选用计算设备的容许计数范围内的最大正整数,γ<sub>sr</sub>表示次级用户发送端s到接收端r的功率增益比例系数,γ<sub>se</sub>表示次级用户发送端s到监听用户接收端e的功率增益比例系数,<img file="FDA0000815702300000034.GIF" wi="75" he="77" />表示平方根操作;(5c)按照下式,计算次级用户安全能效最佳发送功率的迭代真实值:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>F</mi><mi>k</mi><mi>n</mi></msubsup><mo>+</mo><mo>|</mo><msubsup><mi>F</mi><mi>k</mi><mi>n</mi></msubsup><mo>|</mo></mrow><mn>2</mn></mfrac></mrow>]]></math><img file="FDA0000815702300000035.GIF" wi="309" he="149" /></maths>其中,<img file="FDA0000815702300000036.GIF" wi="67" he="71" />表示第k次迭代时次级用户安全能效最佳发送功率的迭代真实值,<img file="FDA0000815702300000037.GIF" wi="67" he="70" />表示第k次迭代时次级用户安全能效最佳发送功率的迭代模糊值,k表示subgradient迭代算法的迭代次数,k=1,2,3,....,K,K表示所选用计算设备的容许计数范围内的最大正整数,|·|表示绝对值操作;(6)计算拉格朗日乘子迭代误差:利用功率约束算法,计算平均干扰功率约束值对应的拉格朗日乘子迭代误差ζ<sub>μ</sub>和平均发送功率约束值对应的拉格朗日乘子迭代误差ζ<sub>λ</sub>;(7)判断平均干扰功率约束值对应的拉格朗日乘子迭代误差和平均发送功率约束值对应的拉格朗日乘子迭代误差是否满足拉格朗日乘子迭代终止条件,若是,则执行步骤(8),否则,将subgradient迭代算法的迭代次数加1后执行步骤(4);(8)按照下式,计算安全能效函数:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>f</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>&eta;</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mi>E</mi><mo>{</mo><mo>&lsqb;</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow></mtd><mtd><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr><mtr><mtd><mrow><mi>E</mi><mo>{</mo><mo>&lsqb;</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo><mo>-</mo><msup><mi>&eta;</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msup><mi>E</mi><mo>{</mo><msubsup><mi>&beta;P</mi><mi>k</mi><mi>n</mi></msubsup><mo>+</mo><msub><mi>P</mi><mi>C</mi></msub><mo>}</mo></mrow></mtd><mtd><mrow><mi>n</mi><mo>&GreaterEqual;</mo><mn>2</mn></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000815702300000041.GIF" wi="1533" he="231" /></maths>其中,f<sub>n</sub>(η)表示第n次迭代时安全能效函数,n表示次级用户安全能效的迭代次数,n=1,2,3,...,N,N表示次级用户安全能效的最大迭代次数,η次级用户安全能效,γ<sub>sr</sub>表示次级用户发送端s到接收端r的功率增益比例系数,γ<sub>se</sub>表示次级用户发送端s到监听用户接收端e的功率增益比例系数,<img file="FDA0000815702300000044.GIF" wi="59" he="57" />表示第k次迭代时次级用户安全能效最佳发送功率的迭代真实值,η<sup>n‑1</sup>次级用户第n‑1次迭代时次级用户安全能效,β表示次级用户发送机的功率放大因子,P<sub>C</sub>表示固定电路C消耗的功率,log<sub>2</sub>(·)表示以2为底的对数操作,E{·}表示求数学期望操作;(9)按照下式,计算次级用户安全能效:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msup><mi>&eta;</mi><mi>n</mi></msup><mo>=</mo><mfrac><mrow><mi>E</mi><mo>{</mo><mo>&lsqb;</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>r</mi></mrow></msub><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>&gamma;</mi><mrow><mi>s</mi><mi>e</mi></mrow></msub><msubsup><mi>P</mi><mi>k</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>&rsqb;</mo><mo>}</mo></mrow><mrow><mi>E</mi><mo>{</mo><msubsup><mi>&beta;P</mi><mi>k</mi><mi>n</mi></msubsup><mo>+</mo><msub><mi>P</mi><mi>C</mi></msub><mo>}</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000815702300000042.GIF" wi="885" he="193" /></maths>其中,η<sup>n</sup>表示第n次迭代时次级用户安全能效,n表示次级用户安全能效的迭代次数,n=1,2,3,...,N,N表示次级用户安全能效的最大迭代次数,γ<sub>sr</sub>表示次级用户发送端s到接收端r的功率增益比例系数,γ<sub>se</sub>表示次级用户发送端s到监听用户接收端e的功率增益比例系数,<img file="FDA0000815702300000043.GIF" wi="69" he="63" />表示第k次迭代时次级用户安全能效最佳发送功率的迭代真实值,k表示subgradient迭代算法的迭代次数,k=1,2,3,....,K,K表示所选用计算设备的容许计数范围内的最大正整数,β表示次级用户发送机的功率放大因子,P<sub>C</sub>表示固定电路C消耗的功率,log<sub>2</sub>(·)表示以2为底的对数操作,E{·}表示求数学期望操作;(10)判断是否满足安全能效终止条件,若是,则执行步骤(11),否则,将次级用户安全能效的迭代次数加1后执行步骤(3);(11)得到次级用户安全能效最佳发送功率和次级用户最大安全能效。
地址 710071 陕西省西安市太白南路2号