发明名称 一种低复杂度多载波PLC中继系统功率分配方法
摘要 本发明公开了一种低复杂度多载波PLC中继系统功率分配方法,包括以下步骤:首先确定系统所使用的子载波,利用导频信道估计方法得到各子载波上的信道系数,设定系统总速率目标值并初始化系统功率分配;然后利用凹凸优化思想将非凸的系统功率最小化问题近似为凸问题,并利用拉格朗日对偶方法迭代地求解近似的系统功率最小化问题得到最后的系统发射功率分配值;最后系统根据计算得到的发射功率值设定源端、中继各子载波上的发射功率,实现多载波PLC系统的信息传输。本发明通过利用凹凸优化方法和拉格朗日对偶方法设计PLC中继系统的功率分配,从而在达到系统总速率要求的同时降低PLC中继系统的总功率。
申请公布号 CN105451324A 申请公布日期 2016.03.30
申请号 CN201510909431.4 申请日期 2015.12.10
申请人 嘉兴国电通新能源科技有限公司;浙江理工大学;北京国电通网络技术有限公司 发明人 史清江;唐剑聪;黄品源;齐世强;潘博;苏广波
分类号 H04W52/34(2009.01)I;H04B3/54(2006.01)I 主分类号 H04W52/34(2009.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种低复杂度多载波PLC中继系统功率分配方法,其特征在于,包括以下步骤:步骤1:系统确定选择使用的子载波集合,总的子载波数为K;利用导频方法进行信道估计得到各子载波上的信道系数<img file="FDA0000873177070000011.GIF" wi="138" he="95" />k=1,2…K,其中<img file="FDA0000873177070000012.GIF" wi="106" he="95" />表示L<sub>1</sub>∈{S,R}到L<sub>2</sub>∈{R,D}之间第k个子载波的信道系数,S表示源节点,R表示中继节点,D表示目的节点;设定总速率设计目标值q;步骤2:初始化迭代次数:n=0,令P<sup>[1]</sup>=P<sup>[2]</sup>=...=P<sup>[K]</sup>,并对等式C({P<sup>[k]</sup>})=q采用二分法求得初始功率分配(P<sup>[k]</sup>)<sup>(n)</sup>,同时令<img file="FDA0000873177070000013.GIF" wi="311" he="110" />并计算<img file="FDA0000873177070000014.GIF" wi="382" he="134" />其中P<sup>[k]</sup>是系统中源、中继节点在第k个子载波上使用的发射功率,<img file="FDA0000873177070000015.GIF" wi="308" he="110" />分别是第n次迭代的系统总传输功率值和所求得的第k个子载波当前发射功率,C({P<sup>[k]</sup>})表示系统总速率函数;步骤3:在<img file="FDA0000873177070000016.GIF" wi="134" he="103" />处,对系统总功率最小化问题中的系统总速率约束进行凸逼近,得到<img file="FDA0000873177070000017.GIF" wi="230" he="109" />从而将原问题近似为如下凸问题:<math><![CDATA[<mrow><munder><mi>min</mi><mrow><mo>{</mo><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>}</mo></mrow></munder><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup></mrow>]]></math><img file="FDA0000873177070000018.GIF" wi="222" he="143" /></maths><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mover><mi>C</mi><mo>~</mo></mover><mrow><mo>(</mo><mrow><mo>{</mo><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>}</mo></mrow><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mi>q</mi></mrow>]]></math><img file="FDA0000873177070000019.GIF" wi="317" he="103" /></maths>P<sup>[k]</sup>≥0,k=1,2…K其中,<math><![CDATA[<mrow><mover><mi>C</mi><mo>~</mo></mover><mrow><mo>(</mo><mo>{</mo><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>}</mo><mo>)</mo></mrow><mover><mo>=</mo><mi>&Delta;</mi></mover><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><msup><mi>a</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>+</mo><msup><mi>b</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><msqrt><mrow><msup><mi>c</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>+</mo><msup><mi>d</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup></mrow></msqrt><mo>+</mo><msup><mi>e</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA00008731770700000110.GIF" wi="1085" he="142" /></maths><math><![CDATA[<mrow><msup><mi>a</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>=</mo><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><mo>+</mo><mn>2</mn><mfrac><msub><mi>x</mi><mn>0</mn></msub><msub><mi>t</mi><mn>0</mn></msub></mfrac><msqrt><mrow><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>R</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><msubsup><mi>&gamma;</mi><mrow><mi>R</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup></mrow></msqrt><mo>-</mo><mfrac><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><msubsup><mi>t</mi><mn>0</mn><mn>2</mn></msubsup></mfrac><mrow><mo>(</mo><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>R</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><mo>+</mo><msubsup><mi>&gamma;</mi><mrow><mi>R</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA00008731770700000111.GIF" wi="862" he="150" /></maths><math><![CDATA[<mrow><msup><mi>b</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>=</mo><mn>2</mn><mfrac><msub><mi>x</mi><mn>0</mn></msub><msub><mi>t</mi><mn>0</mn></msub></mfrac><mo>,</mo></mrow>]]></math><img file="FDA00008731770700000119.GIF" wi="243" he="127" /></maths><math><![CDATA[<mrow><msup><mi>c</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>=</mo><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><mo>+</mo><mn>2</mn><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>R</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><msup><mover><mi>P</mi><mo>~</mo></mover><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA00008731770700000112.GIF" wi="500" he="85" /></maths><math><![CDATA[<mrow><msup><mi>d</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>=</mo><mo>-</mo><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>R</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><msup><mrow><mo>(</mo><msup><mover><mi>P</mi><mo>~</mo></mover><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mrow>]]></math><img file="FDA00008731770700000113.GIF" wi="461" he="111" /></maths><math><![CDATA[<mrow><msup><mi>e</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><msubsup><mi>x</mi><mn>0</mn><mn>2</mn></msubsup><msubsup><mi>t</mi><mn>0</mn><mn>2</mn></msubsup></mfrac><mo>,</mo></mrow>]]></math><img file="FDA00008731770700000114.GIF" wi="270" he="146" /></maths><math><![CDATA[<mrow><msub><mi>x</mi><mn>0</mn></msub><mo>=</mo><msup><mover><mi>P</mi><mo>~</mo></mover><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><msqrt><mrow><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>R</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><msubsup><mi>&gamma;</mi><mrow><mi>R</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup></mrow></msqrt><mo>+</mo><msqrt><mrow><msup><mover><mi>P</mi><mo>~</mo></mover><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>D</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mover><mi>P</mi><mo>~</mo></mover><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><msubsup><mi>&gamma;</mi><mrow><mi>S</mi><mi>R</mi></mrow><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msubsup><mo>)</mo></mrow></mrow></msqrt><mo>,</mo></mrow>]]></math><img file="FDA00008731770700000115.GIF" wi="854" he="118" /></maths><img file="FDA00008731770700000116.GIF" wi="492" he="103" />节点L<sub>1</sub>∈{S,R}到L<sub>2</sub>∈{R,D}之间第k个子载波的信道归一化增益<img file="FDA00008731770700000117.GIF" wi="317" he="190" /><img file="FDA00008731770700000118.GIF" wi="94" he="87" />是L<sub>2</sub>∈{R,D}处第k个子载波上的噪声功率;步骤4:利用对偶方法求解得到上述凸近似问题的最优解<img file="FDA0000873177070000021.GIF" wi="254" he="189" />更新迭代次数:n=n+1,令<math><![CDATA[<mrow><msup><mrow><mo>(</mo><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>=</mo><msup><mrow><mo>(</mo><mover><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>~</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><mo>*</mo><mo>)</mo></mrow></msup><mo>,</mo></mrow>]]></math><img file="FDA0000873177070000022.GIF" wi="394" he="174" /></maths>k=1,2…K,并计算<math><![CDATA[<mrow><msubsup><mi>P</mi><mi>&Sigma;</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msup><mrow><mo>(</mo><msup><mi>P</mi><mrow><mo>&lsqb;</mo><mi>k</mi><mo>&rsqb;</mo></mrow></msup><mo>)</mo></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mo>;</mo></mrow>]]></math><img file="FDA0000873177070000023.GIF" wi="383" he="134" /></maths>步骤5:判断<img file="FDA0000873177070000024.GIF" wi="315" he="92" />是否成立,其中ε<sub>1</sub>表示判定阈值,其值在0.001~0.000001之间,如果成立则令<img file="FDA0000873177070000025.GIF" wi="319" he="111" />k=1,2…K,然后重复步骤3‑5;否则输出问题最后的解(P<sup>[k]</sup>)<sup>(*)</sup>=(P<sup>[k]</sup>)<sup>(n)</sup>,k=1,2…K;步骤6:多载波PLC中继系统中源、中继节点按照{(P<sup>[k]</sup>)<sup>(*)</sup>}设定各个子载波上的发射功率,从而实现PLC系统收发两端的信息传输。
地址 314031 浙江省嘉兴市秀洲区秀园路966号智富中心21幢212室