发明名称 基于OFDM的PLC系统的吞吐量优化方法
摘要 本发明公开了一种基于OFDM的PLC系统的吞吐量优化方法,包括以下步骤:首先发送端设定所使用子载波集合、总传输功率约束值以及每个子载波允许的最大传输功率约束值,然后引入权重,将功率约束下的吞吐量最大化问题等价为加权最小化均方误差问题,再利用块坐标下降算法和二分法迭代求解得到最终功率分配方案,最后发送端按照功率分配方案设定每个子载波的传输功率,从而实现PLC系统的业务传输。由于功率约束下的吞吐量优化问题是一个无法直接的问题,因此本发明通过巧妙构思将该问题等价为加权最小化均方误差问题,使得问题得以采用简单方法迭代求解,从而最大化系统的吞吐量。
申请公布号 CN105282061A 申请公布日期 2016.01.27
申请号 CN201510570058.4 申请日期 2015.09.09
申请人 浙江理工大学 发明人 彭成;史清江;徐伟强
分类号 H04L25/02(2006.01)I;H04L25/03(2006.01)I;H04L27/38(2006.01)I 主分类号 H04L25/02(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种基于OFDM的PLC系统的吞吐量优化方法,其特征在于,包括以下步骤:步骤1:发送端确定所使用的子载波集合A<sup>use</sup>,并求出L=|A<sup>use</sup>|,L表示集合A<sup>use</sup>的元素个数,同时设定每个子载波允许最大传输功率约束值P<sub>max</sub>(m),m∈A<sup>use</sup>以及总传输功率约束值P<sub>total</sub>;步骤2:引入权重值{t<sub>m</sub>},将功率约束下的连续型吞吐量优化问题等价为加权最小化均方差问题,即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><munder><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mrow><mo>{</mo><msub><mi>q</mi><mi>m</mi></msub><mo>,</mo><msub><mi>t</mi><mi>m</mi></msub><mo>,</mo><msub><mi>e</mi><mi>m</mi></msub><mo>}</mo></mrow></munder><msub><mi>F</mi><mn>0</mn></msub><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mrow><mo>(</mo><mi>l</mi><mi>o</mi><mi>g</mi><mo>(</mo><msub><mi>t</mi><mi>m</mi></msub><mo>)</mo><mo>-</mo><msub><mi>t</mi><mi>m</mi></msub><msub><mi>e</mi><mi>m</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000798907790000011.GIF" wi="646" he="135" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi> </mi><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msubsup><mi>q</mi><mi>m</mi><mn>2</mn></msubsup><mo>&le;</mo><msub><mi>P</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub></mrow>]]></math><img file="FDA0000798907790000012.GIF" wi="319" he="138" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mn>0</mn><mo>&le;</mo><msubsup><mi>q</mi><mi>m</mi><mn>2</mn></msubsup><mo>&le;</mo><msub><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>,</mo><mo>&ForAll;</mo><mi>m</mi><mo>&Element;</mo><msup><mi>A</mi><mrow><mi>u</mi><mi>s</mi><mi>e</mi></mrow></msup></mrow>]]></math><img file="FDA0000798907790000013.GIF" wi="582" he="79" /></maths>其中F<sub>0</sub>表示子载波间隔,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>e</mi><mi>m</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><msup><mrow><mo>(</mo><msub><mi>u</mi><mi>m</mi></msub><msub><mi>h</mi><mi>m</mi></msub><msub><mi>q</mi><mi>m</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>u</mi><mi>m</mi><mn>2</mn></msubsup><mi>&Gamma;</mi><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mi>W</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>q</mi><mi>k</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>u</mi><mi>m</mi><mn>2</mn></msubsup><mi>&Gamma;</mi><mi>N</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000798907790000014.GIF" wi="1006" he="134" /></maths>表示MMSE误差估计值,<img file="FDA0000798907790000015.GIF" wi="586" he="100" />α(m)表示第m个子载波的信道增益,P(m)表示第m个子载波的发射功率,u<sub>m</sub>表示第m个子载波的MMSE均衡器,Γ表示信噪比差额,W表示干扰矩阵,其中W(m<sub>0</sub>,m)表示第m个子载波对第m<sub>0</sub>个子载波的干扰,<img file="FDA0000798907790000016.GIF" wi="302" he="71" /><img file="FDA0000798907790000017.GIF" wi="505" he="95" />为噪声功率向量,其中<img file="FDA0000798907790000018.GIF" wi="318" he="87" />表示第m个子载波的噪声功率;步骤3:利用块坐标下降算法和二分法迭代求解该问题,即:初始化:迭代次数n=1、MMSE均衡器<img file="FDA0000798907790000019.GIF" wi="122" he="94" />和相应的<img file="FDA00007989077900000110.GIF" wi="150" he="93" />权重<img file="FDA00007989077900000111.GIF" wi="138" he="94" />其中:<img file="FDA00007989077900000112.GIF" wi="216" he="78" />和<img file="FDA00007989077900000113.GIF" wi="70" he="78" />分别表示第m个子载波的第n次迭代所求的MMSE均衡器、MMSE误差估计值和权重值,然后计算加权最小化均方误差问题的目标值<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mover><mo>=</mo><mi>&Delta;</mi></mover><msub><mi>F</mi><mn>0</mn></msub><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><mrow><mo>(</mo><mi>l</mi><mi>o</mi><mi>g</mi><mo>(</mo><msubsup><mi>t</mi><mi>m</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>)</mo><mo>-</mo><msubsup><mi>t</mi><mi>m</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msubsup><msubsup><mi>e</mi><mi>m</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00007989077900000114.GIF" wi="703" he="130" /></maths>步骤4:更新迭代次数n=n+1,利用二分法求解子载波传输功率集合<img file="FDA00007989077900000115.GIF" wi="230" he="95" />P<sup>(n)</sup>(m)表示第m个子载波第n次迭代所求传输功率值;步骤5:首先求解MMSE均衡器<img file="FDA0000798907790000021.GIF" wi="776" he="207" />和相应的估计误差<img file="FDA0000798907790000022.GIF" wi="397" he="79" />然后求出权重<img file="FDA0000798907790000023.GIF" wi="493" he="135" />从而求出相应的加权最小化均方差问题的目标值<img file="FDA0000798907790000024.GIF" wi="127" he="77" />其中<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mi>q</mi><mi>m</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msubsup><mover><mo>=</mo><mi>&Delta;</mi></mover><msqrt><mrow><msup><mi>P</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000798907790000025.GIF" wi="349" he="93" /></maths>步骤6:判断<img file="FDA0000798907790000026.GIF" wi="351" he="165" />是否成立,其中ε为判定阈值,满足则输出{P<sup>(n)</sup>(m)},即为原问题的最终解,即<img file="FDA0000798907790000027.GIF" wi="402" he="94" />其中P<sup>*</sup>(m)表示第m个子载波最终传输功率,从而求出最终吞吐量<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>R</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow><mo>*</mo></msubsup><mo>=</mo><msub><mi>F</mi><mn>0</mn></msub><munderover><mo>&Sigma;</mo><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><mi>S</mi><mi>I</mi><mi>N</mi><mi>R</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><msup><mi>P</mi><mo>*</mo></msup><mo>(</mo><mi>m</mi><mo>)</mo><mo>)</mo></mrow></mrow><mi>&Gamma;</mi></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000798907790000028.GIF" wi="702" he="175" /></maths>其中<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mi>S</mi><mi>I</mi><mi>N</mi><mi>R</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><msup><mi>P</mi><mo>*</mo></msup><mo>(</mo><mi>m</mi><mo>)</mo><mo>)</mo></mrow><mover><mo>=</mo><mi>&Delta;</mi></mover><mfrac><mrow><msup><mrow><mo>|</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><msup><mi>P</mi><mo>*</mo></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow><mrow><mo>&lsqb;</mo><msup><mi>WP</mi><mo>*</mo></msup><mo>&rsqb;</mo><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&sigma;</mi><mi>b</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000798907790000029.GIF" wi="838" he="175" /></maths><img file="FDA00007989077900000210.GIF" wi="694" he="101" />表示最终传输功率分配向量;否则,重复步骤4到步骤6;步骤7:发送端按照最终功率分配方案设定每个子载波传输功率,从而实现PLC系统的业务传输。
地址 310018 浙江省杭州市江干经济开发区白杨街道2号大街928号