发明名称 基于OFDM的PLC系统的低复杂度比特位加载方法
摘要 本发明公开了一种基于OFDM的PLC系统的低复杂度比特位加载方法,包括以下步骤:首先发送端确定所使用子载波集合、总传输功率约束值以及每个子载波允许的最大传输功率约束值,引入权重,将连续型吞吐量优化问题等价为加权最小化均方误差问题,并结合块坐标下降算法和二分法迭代求解该问题,然后对所求得的非整数型子载波比特数进行整数化,以此为基础,再采用贪婪思想迭代求得整数比特位和功率分配方案,发送端最后根据分配方案将信号比特流映射到每个子载波上并设定子载波传输功率,实现PLC系统收发两端的业务传输。本发明方法不仅能够优化系统性能还能使得比特位加载问题的迭代次数大幅减少,从而降低系统的计算复杂度。
申请公布号 CN105162566A 申请公布日期 2015.12.16
申请号 CN201510570072.4 申请日期 2015.09.09
申请人 嘉兴国电通新能源科技有限公司;浙江理工大学;北京国电通网络技术有限公司 发明人 史清江;彭成;齐世强;潘博;秦鹏;苏广波
分类号 H04L5/00(2006.01)I;H04L27/26(2006.01)I;H04B3/54(2006.01)I 主分类号 H04L5/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 邱启旺
主权项 一种基于OFDM的PLC系统的低复杂度比特位加载方法,其特征在于,包括以下步骤:步骤1:发送端确定所使用子载波数目L、每个子载波允许最大传输功率约束值P<sub>max</sub>(m),m∈A<sup>use</sup>以及总传输功率约束值P<sub>total</sub>,其中A<sup>use</sup>为发送端所使用的子载波集合,且L=|A<sup>use</sup>|,|A<sup>use</sup>|表示集合A<sup>use</sup>的元素个数;步骤2:引入权重,将连续型吞吐量优化问题等价为加权最小化均方误差问题,然后结合块坐标下降算法和二分法迭代求解得到一个非整数型子载波比特数集合{C<sub>InI</sub>(m)},其中C<sub>InI</sub>(m),m∈A<sup>use</sup>表示求解连续型吞吐量优化问题后得到的第m个子载波上的非整数比特位数;步骤3:根据式子<img file="FDA00007990148100000112.GIF" wi="1260" he="91" />m∈A<sup>use</sup>对{C<sub>InI</sub>(m)}整数化得到一个整数比特位分配集合{b<sub>InI</sub>(m)},其中Λ<sub>d</sub>(·)表示整数化函数,b<sub>InI</sub>(m)表示第m个子载波最初的比特位数,<img file="FDA0000799014810000012.GIF" wi="58" he="61" />表示子载波采用的M进制‑QAM调制方式所对应的整数比特位数有序集合,<img file="FDA0000799014810000013.GIF" wi="487" he="83" />特别地,M=0表示在<img file="FDA0000799014810000014.GIF" wi="53" he="60" />集合中所对应的比特位数为0,即子载波上不携带比特位数;步骤4:根据贪婪思想,在{b<sub>InI</sub>(m)}基础上进行比特位加载,即初始化迭代次数k=1,令{b<sup>(k)</sup>(m)=b<sub>InI</sub>(m)},则有<img file="FDA0000799014810000015.GIF" wi="739" he="101" />表示子载波比特位分配向量,并求出对应的子载波功率分配向量<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mover><mo>=</mo><mi>&Delta;</mi></mover><msup><mrow><mo>&lsqb;</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>L</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mi>T</mi></msup><mo>=</mo><msup><mi>M</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mrow><mo>(</mo><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mi>N</mi><mo>,</mo></mrow>]]></math><img file="FDA0000799014810000016.GIF" wi="1225" he="115" /></maths>其中:(·)<sup>(k)</sup>表示第k轮迭代所求的对应值,b(m)表示第m个子载波上的比特位数,P(m)表示第m个子载波的传输功率,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>M</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mover><mo>=</mo><mi>&Delta;</mi></mover><msup><mrow><mo>(</mo><mi>I</mi><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mrow><mo>(</mo><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mi>W</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mrow><mo>(</mo><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>)</mo></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mover><mo>=</mo><mi>&Delta;</mi></mover><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mrow><mo>(</mo><msubsup><mi>&gamma;</mi><mn>1</mn><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>&gamma;</mi><mi>m</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mo>...</mo><mo>,</mo><msubsup><mi>&gamma;</mi><mi>L</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000799014810000017.GIF" wi="1396" he="172" /></maths><img file="FDA0000799014810000018.GIF" wi="436" he="210" />Γ和α(m)分别表示信噪比差额和信道增益,<img file="FDA0000799014810000019.GIF" wi="472" he="99" />表示噪声功率向量,其中<img file="FDA00007990148100000110.GIF" wi="315" he="90" />表示第m个子载波的噪声功率,W表示干扰矩阵;步骤5:更新迭代次数k=k+1,L个子载波根据<img file="FDA00007990148100000111.GIF" wi="52" he="60" />依次跳到下一级来更新自己的比特位数,并记录下每次更新后的比特位分配向量B<sup>(k)</sup>(m),m∈A<sup>use</sup>,即B<sup>(k‑1)</sup>中的第m个子载波根据<img file="FDA0000799014810000021.GIF" wi="46" he="63" />更新比特位后得到的比特位分配向量<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mover><mo>=</mo><mi>&Delta;</mi></mover><msup><mrow><mo>&lsqb;</mo><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><mo>,</mo><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>,</mo><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mo>...</mo><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>L</mi><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mi>T</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000799014810000022.GIF" wi="1515" he="103" /></maths>m∈A<sup>use</sup>,由此可以得到更新后的比特位分配向量集合{B<sup>(k)</sup>(m)},并迭代求出对应矩阵集合<img file="FDA0000799014810000023.GIF" wi="332" he="167" />和对应的子载波功率分配向量集合{P<sup>(k)</sup>(m)},即:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mrow><mo>(</mo><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>(</mo><mi>m</mi><mo>)</mo><mo>)</mo></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mover><mo>=</mo><mi>&Delta;</mi></mover><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mrow><mo>(</mo><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>)</mo></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>+</mo><msubsup><mi>&Delta;&gamma;</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><msub><mi>e</mi><mi>m</mi></msub><msubsup><mi>e</mi><mi>m</mi><mi>T</mi></msubsup></mrow>]]></math><img file="FDA0000799014810000024.GIF" wi="841" he="149" /></maths>和<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>M</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><msup><mrow><mo>(</mo><msub><mi>&Lambda;</mi><mrow><mo>(</mo><msup><mi>B</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>(</mo><mi>m</mi><mo>)</mo><mo>)</mo></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mi>N</mi><mo>,</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="FDA0000799014810000025.GIF" wi="869" he="148" /></maths>其中:<img file="FDA0000799014810000026.GIF" wi="775" he="217" />e<sub>m</sub>为列向量,且只有第m行为1,其余行均为0,<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>W</mi><mi>m</mi></msub><mover><mo>=</mo><mi>&Delta;</mi></mover><msubsup><mi>e</mi><mi>m</mi><mi>T</mi></msubsup><mi>W</mi><mo>,</mo><msubsup><mi>&Delta;&gamma;</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mover><mo>=</mo><mi>&Delta;</mi></mover><mfrac><mrow><mo>(</mo><msup><mn>2</mn><mrow><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></msup><mo>-</mo><msup><mn>2</mn><mrow><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></msup><mo>)</mo><mi>&Gamma;</mi></mrow><mrow><mo>|</mo><mi>&alpha;</mi><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><msup><mo>|</mo><mn>2</mn></msup></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000799014810000027.GIF" wi="894" he="206" /></maths>步骤6:执行下列选择方案,即:<maths num="0007" id="cmaths0007"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msup><mi>m</mi><mo>*</mo></msup><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mi>m</mi></munder><mi>&Delta;</mi><mi>f</mi><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mi>m</mi></munder><mfrac><mrow><mi>S</mi><mrow><mo>(</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>)</mo></mrow><mrow><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>b</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mi>S</mi><mrow><mo>(</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mo>(</mo><msup><mi>m</mi><mo>*</mo></msup><mo>)</mo></mrow><mo>)</mo><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></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>&le;</mo><msup><mi>P</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msup><mrow><mo>(</mo><msup><mi>m</mi><mo>*</mo></msup><mo>)</mo></mrow><mo>&le;</mo><msub><mi>P</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mrow><mo>(</mo><msup><mi>m</mi><mo>*</mo></msup><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000799014810000028.GIF" wi="950" he="435" /></maths>其中:S(x)表示对向量x中所有元素的求和;根据上式如果能够找到满足条件的子载波m<sup>*</sup>,那么更新B<sup>(k)</sup>=B<sup>(k)</sup>(m<sup>*</sup>),P<sup>(k)</sup>=P<sup>(k)</sup>(m<sup>*</sup>),<img file="FDA0000799014810000029.GIF" wi="571" he="156" /><img file="FDA00007990148100000210.GIF" wi="769" he="245" />并重复步骤5;如果找不到,则得到最后的比特位分配向量和功率分配向量分别为B<sup>*</sup>=B<sup>(k‑1)</sup>和P<sup>*</sup>=P<sup>(k‑1)</sup>;步骤7:发送端按照最后得到的比特位分配向量B<sup>*</sup>确定每个子载波对应的调制方式,同时根据P<sup>*</sup>设定每个子载波的传输功率,从而实现PLC系统收发两端的业务传输。
地址 314031 浙江省嘉兴市秀洲区秀园路966号智富中心21幢212室