发明名称 一种拥有低峰均功率比的无线光正交多载波通信方法
摘要 本发明公开了一种拥有低峰均功率比(PAPR)的无线光正交多载波通信方法,按如下步骤进行:在发射端,首先调制后的频域符号采用音调注入(TI)算法,用松弛法转化为半正定规划凸优化问题,通过通用凸优化和随机化方法求解;其次,频域信号共轭对称的映射在子载波上,经过反快速傅里叶变换后添加循环前缀;最后,时域发射信号加入直流偏置驱动发光二极管发射。在接收端,光电二极管实现光电转换;信号经放大滤波,模数转换,快速傅里叶变换后移除循环前缀和共轭对称部分;对信号取模将其恢复在原星座图中;最后解调得到接收符号。本发明能够有效降低无线光通信OFDM系统的PAPR,降低对功率放大器和LED线性度的要求,减小非线性失真,提高接收性能。
申请公布号 CN103618687A 申请公布日期 2014.03.05
申请号 CN201310638461.7 申请日期 2013.12.03
申请人 东南大学 发明人 张华;袁颖;许威
分类号 H04L27/26(2006.01)I;H04L27/34(2006.01)I;H04B10/116(2013.01)I 主分类号 H04L27/26(2006.01)I
代理机构 江苏永衡昭辉律师事务所 32250 代理人 王斌
主权项 1.一种拥有低峰均功率比的无线光正交多载波通信方法,包括以下步骤:1)在发射端,无线光通信OFDM系统的子载波数为N,信源产生的二进制数据源经过多进制正交幅度调制M-QAM,形成待发送的频域信号<maths num="0001"><![CDATA[<math><mrow><mi>X</mi><mo>=</mo><mo>[</mo><msub><mi>X</mi><mn>0</mn></msub><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><msub><mi>X</mi><mrow><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>-</mo><mn>1</mn></mrow></msub><mo>]</mo><mo>;</mo></mrow></math>]]></maths>2)根据音调注入方法,将原有星座点扩展,使得频域发送信号变为<img file="FDA0000428558310000012.GIF" wi="316" he="78" />其中C<sub>k</sub>=p<sub>k</sub>D<sub>k</sub>+jq<sub>k</sub>D<sub>k</sub>被称为注入的音调,则p<sub>k</sub>∈{0,-sgn(R<sub>k</sub>)},q<sub>k</sub>∈{0,-sgn(I<sub>k</sub>)},<img file="FDA0000428558310000013.GIF" wi="355" he="90" />R<sub>k</sub>=Re{X<sub>k</sub>},I<sub>k</sub>=Im{X<sub>k</sub>},M<sub>k</sub>表示QAM信号星座的大小,标量因子ρ≥1,Δ<sub>k</sub>表示符号间的最小距离,p<sub>k</sub>,q<sub>k</sub>的值根据经典TI方法选取,sgn()表示取符号,Re{·}、Im{·}分别表示取信号实部和虚部;3)为了满足无线光通信,使得基带信号为实值,无线光通信OFDM中子载波映射后的频域发射信号为Y<sub>k</sub>,k=0,1,...,N-1,映射按照下面公式进行,满足共轭对称性:<maths num="0002"><![CDATA[<math><mrow><msub><mi>Y</mi><mi>k</mi></msub><mo>=</mo><mfenced open='{' close='' separators=' '><mtable><mtr><mtd><mi>Re</mi><mo>{</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>}</mo><mo>,</mo></mtd><mtd><mi>k</mi><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><mo>,</mo></mtd><mtd><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi>Im</mi><mo>{</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>}</mo><mo>,</mo></mtd><mtd><mi>k</mi><mo>=</mo><mfrac><mi>N</mi><mn>2</mn></mfrac></mtd></mtr><mtr><mtd><msubsup><mover><mi>X</mi><mo>&OverBar;</mo></mover><mrow><mi>N</mi><mo>-</mo><mi>k</mi></mrow><mo>*</mo></msubsup><mo>,</mo></mtd><mtd><mi>k</mi><mo>=</mo><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>+</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn></mtd></mtr></mtable><mo>;</mo></mfenced></mrow></math>]]></maths>4)频域信号经过N点的快速傅里叶变换IFFT后转为实值的时域OFDM信号如下:<maths num="0003"><![CDATA[<math><mrow><msub><mi>y</mi><mi>n</mi></msub><mo>=</mo><mfrac><mn>1</mn><msqrt><mi>N</mi></msqrt></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>Y</mi><mi>k</mi></msub><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;kn</mi><mo>/</mo><mi>N</mi></mrow></msup><mo>=</mo><mfrac><mn>2</mn><msqrt><mi>N</mi></msqrt></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mfrac><mi>N</mi><mn>2</mn></mfrac><mo>-</mo><mn>1</mn></mrow></munderover><mi>Re</mi><mo>{</mo><msub><mover><mi>X</mi><mo>&OverBar;</mo></mover><mi>k</mi></msub><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;kn</mi><mo>/</mo><mi>N</mi></mrow></msup><mo>}</mo><mo>+</mo><msub><mi>K</mi><mi>n</mi></msub><mo>,</mo></mrow></math>]]></maths>其中<img file="FDA0000428558310000021.GIF" wi="728" he="135" />通常OFDM系统中直流部分设置<img file="FDA0000428558310000022.GIF" wi="181" he="78" />此时对应K<sub>n</sub>=0;5)重新映射时域信号功率为:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>y</mi><mi>n</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>4</mn><mi>N</mi></mfrac><msup><mrow><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><mi>Re</mi><mo>{</mo><mrow><mo>(</mo><msub><mi>R</mi><mi>k</mi></msub><mo>+</mo><mi>j</mi><msub><mi>I</mi><mi>k</mi></msub><mo>+</mo><msub><mi>p</mi><mi>k</mi></msub><msub><mi>D</mi><mi>k</mi></msub><mo>+</mo><mi>j</mi><msub><mi>q</mi><mi>k</mi></msub><msub><mi>D</mi><mi>k</mi></msub><mo>)</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;kn</mi><mo>/</mo><mi>N</mi></mrow></msup><mo>}</mo><mo>]</mo></mrow><mn>2</mn></msup><mo>=</mo><mfrac><mn>4</mn><mi>N</mi></mfrac><msup><mi>r</mi><mi>T</mi></msup><msub><mi>S</mi><mi>n</mi></msub><mi>r</mi><mo>,</mo></mrow></math>]]></maths>其中S<sub>n</sub>,n=0,…,N-1是与发射信号有关的(N+1)×(N+1)的已知矩阵,r=[1 p<sup>T</sup> q<sup>T</sup>]<sup>T</sup>,p=[p<sub>0</sub> ... p<sub>N/2-1</sub>]<sup>T</sup>,q=[q<sub>0</sub> ... q<sub>N/2-1</sub>]<sup>T</sup>;决定向量元素p<sub>k</sub>,q<sub>k</sub>原本取自{0,±1},为了方便后续步骤的处理,定义新的变量<maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>p</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><mn>2</mn><msub><mi>p</mi><mi>k</mi></msub><mi>sgn</mi><mrow><mo>(</mo><msub><mi>R</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>,</mo><msub><mover><mi>q</mi><mo>~</mo></mover><mi>k</mi></msub><mo>=</mo><mn>2</mn><msub><mi>q</mi><mi>k</mi></msub><mi>sgn</mi><mrow><mo>(</mo><msub><mi>I</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><mn>1</mn><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mover><mi>r</mi><mo>~</mo></mover><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msup><mover><mi>p</mi><mo>~</mo></mover><mi>T</mi></msup></mtd><mtd><msup><mover><mi>q</mi><mo>~</mo></mover><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><msub><mover><mi>p</mi><mo>~</mo></mover><mn>0</mn></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mover><mi>p</mi><mo>~</mo></mover><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mover><mi>q</mi><mo>~</mo></mover><mn>0</mn></msub><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mover><mi>q</mi><mo>~</mo></mover><mrow><mi>N</mi><mo>/</mo><mn>2</mn><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>将TI算法问题等效为<img file="FDA0000428558310000026.GIF" wi="370" he="108" />其中<img file="FDA0000428558310000027.GIF" wi="64" he="83" />是与信号有关的方阵如下:<maths num="0007"><![CDATA[<math><mrow><msub><mover><mi>S</mi><mo>~</mo></mover><mi>n</mi></msub><mo>=</mo><msup><mi>A</mi><mi>T</mi></msup><msub><mi>S</mi><mi>n</mi></msub><mi>A</mi><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>A</mi><mi>T</mi></msup><msub><mi>S</mi><mi>n</mi></msub><mi>b</mi></mtd><mtd><msub><mn>0</mn><mrow><mi>N</mi><mo>+</mo><mn>1</mn><mo>,</mo><mi>N</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>b</mi><mi>T</mi></msup><msub><mi>S</mi><mi>n</mi></msub><mi>A</mi></mtd></mtr><mtr><mtd><msub><mn>0</mn><mrow><mi>N</mi><mo>,</mo><mi>N</mi><mo>+</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>+</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msup><mi>b</mi><mi>T</mi></msup><msub><mi>S</mi><mi>n</mi></msub><mi>b</mi></mtd><mtd><msub><mn>0</mn><mrow><mn>1</mn><mo>,</mo><mi>N</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mn>0</mn><mrow><mi>N</mi><mo>,</mo><mn>1</mn></mrow></msub></mtd><mtd><msub><mn>0</mn><mrow><mi>N</mi><mo>,</mo><mi>N</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mi>A</mi><mo>=</mo><mi>diag</mi><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>sgn</mi><msup><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>sgn</mi><msup><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mo>,</mo><mi>b</mi><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>sgn</mi><msup><mrow><mo>(</mo><mi>R</mi><mo>)</mo></mrow><mi>T</mi></msup></mtd><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>sgn</mi><msup><mrow><mo>(</mo><mi>I</mi><mo>)</mo></mrow><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>,</mo></mrow></math>]]></maths>diag表示构造对角阵;6)令<img file="FDA00004285583100000210.GIF" wi="195" he="80" />则<img file="FDA00004285583100000211.GIF" wi="49" he="73" />是元素仅包含±1并且对角线元素为+1的一个秩为1的方阵,进一步将秩为1这个条件松弛为<img file="FDA00004285583100000212.GIF" wi="49" he="73" />是半正定的,松弛后的问题转化为标准半正定规划SDP凸优化问题:<img file="FDA00004285583100000213.GIF" wi="774" he="233" />其中min表示最小化,s.t.表示优化问题的约束条件;7)使用通用凸优化工具箱CVX解得半正定矩阵<img file="FDA00004285583100000214.GIF" wi="78" he="75" />采用随机化方法将<img file="FDA0000428558310000031.GIF" wi="51" he="73" />特征值分解,选取其最大的特征值λ<sub>max</sub>和对应的特征向量u,则<img file="FDA0000428558310000032.GIF" wi="41" he="58" />中各个元素取值的概率为<img file="FDA0000428558310000033.GIF" wi="519" he="148" /><img file="FDA0000428558310000034.GIF" wi="558" he="149" />Pr{r=r<sub>0</sub>}表示r=r<sub>0</sub>的概率,按照这个概率分布随机产生多组候选向量,选择使得目标函数最小的候选向量作为最优解<img file="FDA0000428558310000035.GIF" wi="68" he="60" />从而计算最终的p,q以及最终发射信号<img file="FDA0000428558310000036.GIF" wi="80" he="71" />8)将经过TI运算后的信号依次经过共轭映射,N点快速傅里叶变换,加上循环前缀得到时域发送信号,将时域信号经过功率放大器以及数模变换器,送往光通信发射模块,加入直流偏置后驱动LED灯发射;9)在接收端,光电二极管将光信号转换为时域电信号,经过放大、滤波处理后,通过模数转换模块,作N点傅里叶变换后去除循环前缀,得到频域接收信号<img file="FDA0000428558310000037.GIF" wi="79" he="74" />10)频域接收信号<img file="FDA0000428558310000038.GIF" wi="44" he="68" />经过对D<sub>k</sub>取模的模块,将接收信号星座点限制在未扩展的M-QAM星座图中,得到<img file="FDA0000428558310000039.GIF" wi="68" he="71" />最后进行解调,获得最终的接收符号。
地址 210096 江苏省南京市四牌楼2号