发明名称 基于Fountain码的功率和中继联合优化方法及其协作通信方法
摘要 本发明公开一种基于Fountain码的功率和中继联合优化方法,其特征在于采用Fountain码对信号进行编码传输,基于Fountain码的信息累积特点和实际信道环境构建功率和中继联合优化模型,从而在保证信息传输的可靠性前提下,最优化其有效性指标,以实现信息传输的有效性和可靠性的良好折中,在此基础上,求解联合优化模型的理论最优解作为协作通信的依据。本发明还公开了基于上述联合优化方法的协作通信方法。本发明基于Fountain码的功率和中继联合优化方法的协作通信方法非常简单而易于实现,具有很好的应用前景。
申请公布号 CN102427586B 申请公布日期 2013.12.04
申请号 CN201110399899.5 申请日期 2011.12.06
申请人 南京邮电大学;东南大学 发明人 张晖;朱洪波;杨龙祥;赵海涛;王再见;吴鹏;赵远林
分类号 H04L1/00(2006.01)I;H04W24/02(2009.01)I;H04W52/00(2009.01)I 主分类号 H04L1/00(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 艾中兰
主权项 1.一种基于Fountain码的功率和中继联合优化方法,其特征在于采用Fountain码对信号进行编码传输,基于Fountain码的信息累积特点和实际信道环境构建功率和中继联合优化模型,从而在保证信息传输的可靠性前提下,最优化其有效性指标,以实现信息传输的有效性和可靠性的良好折中,在此基础上,求解联合优化模型的理论最优解作为协作通信的依据,所述功率和中继联合优化模型如下所示:<maths num="0001"><![CDATA[<math><mrow><mo>{</mo><msup><mi>C</mi><mo>*</mo></msup><mo>,</mo><msubsup><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>&Delta;</mi><mn>1</mn><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>&Delta;</mi><mn>2</mn><mo>*</mo></msubsup><mo>}</mo><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><mo>{</mo><mi>C</mi><mo>/</mo><mi>C</mi><mo>&Element;</mo><mi>N</mi><mo>}</mo><mo>,</mo><mo>{</mo><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>}</mo><mo>,</mo><mo>{</mo><msub><mi>&Delta;</mi><mi>i</mi></msub><mo>/</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mn>2</mn><mo>}</mo></mrow></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msub><mi>&Delta;</mi><mi>i</mi></msub></mrow></math>]]></maths>subject to0≤P<sub>1S</sub>≤P<sub>max</sub>0≤P<sub>2S</sub>≤P<sub>max</sub>        (1)P<sub>1C</sub>=00≤P<sub>2C</sub>≤P<sub>max</sub>Δ<sub>i</sub>≥0,1≤i≤2<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>SC</mi><mn>1</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo>)</mo></mrow><msub><mi>&Delta;</mi><mn>1</mn></msub><mo>&GreaterEqual;</mo><mi>I</mi></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>S</mi><mo>,</mo><mi>C</mi></mrow></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msubsup><mi>f</mi><mi>jD</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mi>ij</mi></msub><mo>)</mo></mrow><msub><mi>&Delta;</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><mi>I</mi><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>k</mi><mo>&le;</mo><mn>2</mn></mrow></math>]]></maths>其中,Δ<sub>1</sub>表示从源节点S发送数据到中继节点C解出数据的时间间隔,即第1个时隙;Δ<sub>2</sub>表示从中继节点C解出数据到目的节点D解出数据的时间间隔,即第2个时隙,P<sub>1S</sub>,P<sub>2S</sub>,P<sub>1C</sub>,P<sub>2C</sub>表示源节点S和中继节点C分别在2个时隙的发射功率,Ν表示候选中继节点集,P<sub>max</sub>表示任意节点的发射功率最大值,C<sup>*</sup>,<img file="FDA0000381293260000014.GIF" wi="344" he="75" />分别表示变量C,P<sub>1S</sub>,P<sub>2S</sub>,P<sub>2C</sub>,Δ<sub>1</sub>,Δ<sub>2</sub>的最优解,优化模型的最后2个约束条件分别规定了节点C和节点D解出数据所需满足的条件,其中,<img file="FDA0000381293260000015.GIF" wi="181" he="77" />表示在时隙1的节点S和节点C之间的信道传输速率,Ι为数据比特数,<img file="FDA0000381293260000016.GIF" wi="182" he="82" />j=S,C,表示在时隙i的节点j和节点D之间的信道传输速率,<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>SC</mi><mn>1</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>W</mi><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>SC</mi></msub><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub></mrow><mrow><msub><mi>N</mi><mi>C</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msubsup><mi>f</mi><mi>SD</mi><mn>1</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>W</mi><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub></mrow><mrow><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>,</mo><msubsup><mi>f</mi><mi>SD</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mi>W</mi><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mi>CD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>CD</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mrow><mi>W</mi><mi>log</mi></mrow><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>CD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>这里,W表示信号带宽,N<sub>C</sub>和N<sub>D</sub>分别表示以节点C和节点D为接收端的信道的加性高斯噪声功率谱密度,h<sub>SC</sub>、h<sub>SD</sub>和h<sub>CD</sub>分别表示信道l<sub>SC</sub>、l<sub>SD</sub>和l<sub>CD</sub>的信道增益;求解联合优化模型的理论最优解的方法是,将联合优化模型分解成2个子问题:中继选择问题和功率分配问题,首先固定中继选择问题,对联合优化模型进行化简求解其功率分配问题,进而求解中继选择问题,从而获得联合优化模型的理论最优解,求解功率分配问题的方法如下:固定中继选择问题,将优化模型化简为:<maths num="0006"><![CDATA[<math><mrow><mo>{</mo><msubsup><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>&Delta;</mi><mn>1</mn><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>&Delta;</mi><mn>2</mn><mo>*</mo></msubsup><mo>}</mo><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><mo>{</mo><msub><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>}</mo><mo>,</mo><mo>{</mo><msub><mi>&Delta;</mi><mi>i</mi></msub><mo>/</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&le;</mo><mn>2</mn><mo>}</mo></mrow></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><msub><mi>&Delta;</mi><mi>i</mi></msub><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中约束条件同于公式(1),将优化过程分为2个时隙,首先通过理论分析方法获得第1个时隙的最优解;<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>=</mo><msub><mi>P</mi><mi>max</mi></msub></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>&Delta;</mi><mn>1</mn><mo>*</mo></msubsup><mo>=</mo><mfrac><mi>I</mi><mrow><msubsup><mi>f</mi><mi>SC</mi><mn>1</mn></msubsup><mrow><mo>(</mo><msubsup><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mi>I</mi><mrow><mi>W</mi><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>SC</mi></msub><msubsup><mi>P</mi><mrow><mn>1</mn><mi>S</mi></mrow><mo>*</mo></msubsup></mrow><mrow><msub><mi>N</mi><mi>C</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>在此基础上,得到第2个时隙的优化模型,将其等价转化为:<maths num="0009"><![CDATA[<math><mrow><mo>{</mo><msubsup><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow><mo>*</mo></msubsup><mo>,</mo><msubsup><mi>&Delta;</mi><mn>2</mn><mo>*</mo></msubsup><mo>}</mo><mo>=</mo><munder><mrow><mi>arg</mi><mi>max</mi></mrow><mrow><mo>{</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>}</mo><mo>,</mo><mo>{</mo><msub><mi>&Delta;</mi><mn>2</mn></msub><mo>}</mo></mrow></munder><msubsup><mi>f</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><mo>=</mo><munder><mrow><mi>arg</mi><mi>max</mi></mrow><mrow><mo>{</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>}</mo><mo>,</mo><mo>{</mo><msub><mi>&Delta;</mi><mn>2</mn></msub><mo>}</mo></mrow></munder><mrow><mo>(</mo><msub><mrow><mi>W</mi><mi>log</mi></mrow><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mi>CD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mrow><mi>W</mi><mi>log</mi></mrow><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>CD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>subject to0≤P<sub>2S</sub>≤P<sub>max</sub>0≤P<sub>2C</sub>≤P<sub>max</sub><maths num="0011"><![CDATA[<math><mrow><msub><mi>&Delta;</mi><mn>2</mn></msub><mo>=</mo><mfrac><msubsup><mi>I</mi><mi>D</mi><mn>2</mn></msubsup><mrow><msubsup><mi>f</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>&GreaterEqual;</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA0000381293260000031.GIF" wi="64" he="78" />表示在第二个时隙节点D解码所需信息量,<img file="FDA0000381293260000032.GIF" wi="641" he="285" /><maths num="0012"><![CDATA[<math><mrow><msubsup><mi>f</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><msub><mrow><mi>W</mi><mi>log</mi></mrow><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mi>CD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow><mo>+</mo><msub><mrow><mi>W</mi><mi>log</mi></mrow><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><mrow><msub><mi>h</mi><mi>CD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mi>SD</mi></msub><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>+</mo><msub><mi>N</mi><mi>D</mi></msub><mi>W</mi></mrow></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>表示在第2个时隙到节点D的总的传输速率,从而求解得到第2个时隙的最优解:<maths num="0013"><![CDATA[<math><mrow><msubsup><mi>&Delta;</mi><mn>2</mn><mo>*</mo></msubsup><mo>=</mo><mfrac><msubsup><mi>I</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mi>max</mi><msubsup><mi>f</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>24</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中<img file="FDA0000381293260000035.GIF" wi="314" he="65" /><img file="FDA0000381293260000036.GIF" wi="1968" he="255" />求解功率分配问题基础上,求解中继选择问题,从而获得联合优化模型的理论最优解,具体方法是:比较候选中继节点集Ν中各中继节点相应的最优传输时延,获得最佳中继节点C<sup>*</sup>:<maths num="0014"><![CDATA[<math><mrow><msup><mi>C</mi><mo>*</mo></msup><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mrow><mi>C</mi><mo>&Element;</mo><mi>N</mi></mrow></munder><msubsup><mi>&Delta;</mi><mi>C</mi><mo>*</mo></msubsup><mo>,</mo></mrow></math>]]></maths><maths num="0015"><![CDATA[<math><mrow><msubsup><mi>&Delta;</mi><mi>C</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>&Delta;</mi><mn>1</mn><mo>*</mo></msubsup><mo>+</mo><msubsup><mi>&Delta;</mi><mn>2</mn><mo>*</mo></msubsup></mrow></math>]]></maths>表示最优传输时延,<maths num="0016"><![CDATA[<math><mrow><msubsup><mi>&Delta;</mi><mn>2</mn><mo>*</mo></msubsup><mo>=</mo><mfrac><msubsup><mi>I</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mi>max</mi><msubsup><mi>f</mi><mi>D</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>S</mi></mrow></msub><mo>,</mo><msub><mi>P</mi><mrow><mn>2</mn><mi>C</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>表示第2个时隙的最优解,因此,最佳中继节点C<sup>*</sup>和相应的最优功率分配即为联合优化模型的理论最优解。
地址 210003 江苏省南京市鼓楼区新模范马路66号