发明名称 基于岭回归分析的相干光通信色散补偿方法
摘要 本发明涉及一种基于的基于岭回归分析的相干光通信色散补偿方法,本发明的意义在于为光通信系统,提供了一种仅依靠小数据量且收敛速度快的色散补偿方法。通过构建基于岭回归分析的色散补偿性能函数,将色散补偿问题转化为求解无约束优化问题,最后构造批处理形式迭代算法求解该优化问题。该方法不仅克服了传统LMS/LS色散补偿方法的系数矩阵病态问题,而且无论从收敛速度、收敛后的性能或者是所依赖的数据量要求,均优越于传统的LMS/LS色散补偿方法。该方法适用于低阶和高阶的多值复数域正交幅度调制(QAM)信号以及多进制相移键控(MPSK)信号的色散补偿问题。
申请公布号 CN103346839B 申请公布日期 2015.12.09
申请号 CN201310287833.6 申请日期 2013.07.09
申请人 温州大学 发明人 阮秀凯;陈孝敬;张耀举;谈燕花;吴平;施肖菁;罗素琴
分类号 H04B10/2513(2013.01)I 主分类号 H04B10/2513(2013.01)I
代理机构 温州瓯越专利代理有限公司 33211 代理人 于艳玲
主权项 一种基于岭回归分析的相干光通信色散补偿方法,其特征在于,包括以下步骤:(1)光通讯系统中分集光接收机接收到信号光和本振激光器产生的本振光经过混频后,形成的叠加光场经光探测器转换为光电流信号,该光电流信号经相关处理后输出基带信号;标记接收端实现色散补偿的均衡器的输出信号为y<sub>k</sub>,y<sub>k</sub>=w<sup>T</sup>x<sub>k</sub>,其中w=[w<sub>0</sub>,w<sub>1</sub>,...,w<sub>L‑1</sub>]<sup>T</sup>是均衡器向量;x<sub>k</sub>=[x<sub>k</sub>,x<sub>k‑1</sub>...,x<sub>k‑L+1</sub>]<sup>T</sup>是基带信号接收矩阵,L是均衡器的抽头数,上标T表示矩阵转置;(2)构造岭回归分析方法的无约束最优化问题,其中无约束最优化问题的数学模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>w</mi><mo>^</mo></mover><mo>=</mo><mi>arg</mi><mi> </mi><munder><mi>min</mi><mi>w</mi></munder><msup><mrow><mo>|</mo><mrow><msup><mrow><mo>|</mo><mi>y</mi><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mi>R</mi><mi>p</mi></msub></mrow><mo>|</mo></mrow><mi>q</mi></msup><mo>+</mo><mi>&lambda;</mi><msup><mrow><mo>|</mo><mi>w</mi><mo>|</mo></mrow><mi>q</mi></msup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mi>L</mi><mo>-</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msup><mrow><mo>|</mo><mrow><msubsup><mi>x</mi><mi>k</mi><mi>T</mi></msubsup><mi>w</mi></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mi>R</mi><mi>p</mi></msub><mo>)</mo></mrow><mi>q</mi></msup><mo>+</mo><mi>&lambda;</mi><msup><mrow><mo>|</mo><mi>w</mi><mo>|</mo></mrow><mi>q</mi></msup><mo>;</mo></mrow>]]></math><img file="FDA0000782251730000011.GIF" wi="1252" he="147" /></maths>其中:<img file="FDA0000782251730000012.GIF" wi="53" he="63" />为w的估计值,R<sub>p</sub>是星座特征常数,N为接收数据长度,q为正数,λ为岭参数,λ≥0;y=[y<sub>L‑1</sub>,y<sub>L</sub>,…y<sub>N</sub>]<sup>T</sup>;(3)令q=2,得到性能函数J(w);<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><msup><mrow><mo>|</mo><mi>y</mi><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mi>R</mi><mi>p</mi></msub><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mi>L</mi><mo>-</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>(</mo><msup><mrow><mo>|</mo><mrow><msubsup><mi>x</mi><mi>k</mi><mi>T</mi></msubsup><mi>w</mi></mrow><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msub><mi>R</mi><mi>p</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mi>&lambda;w</mi><mi>H</mi></msup><mi>w</mi><mo>;</mo></mrow>]]></math><img file="FDA0000782251730000013.GIF" wi="1186" he="143" /></maths>(4)计算J(w)的一阶和二阶偏导数;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mfrac><mrow><mo>&part;</mo><mi>J</mi><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>w</mi></mrow></mfrac><mo>=</mo><mrow><mo>(</mo><msup><mi>X</mi><mi>H</mi></msup><mi>D</mi><mi>X</mi><mo>+</mo><mi>&lambda;</mi><mi>I</mi><mo>)</mo></mrow><mi>w</mi><mo>-</mo><msub><mi>R</mi><mi>p</mi></msub><msup><mi>X</mi><mi>H</mi></msup><mi>y</mi><mo>,</mo><mfrac><mrow><msup><mo>&part;</mo><mn>2</mn></msup><mi>J</mi><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow></mrow><mrow><mo>&part;</mo><mi>w</mi><mo>&part;</mo><msup><mi>w</mi><mi>H</mi></msup></mrow></mfrac><mo>=</mo><msup><mi>X</mi><mi>H</mi></msup><mi>D</mi><mi>X</mi><mo>+</mo><mi>&lambda;</mi><mi>I</mi><mo>;</mo></mrow>]]></math><img file="FDA0000782251730000014.GIF" wi="1283" he="148" /></maths>其中I为单位矩阵,<img file="FDA0000782251730000015.GIF" wi="1050" he="316" /><img file="FDA0000782251730000021.GIF" wi="797" he="390" />(5)令<img file="FDA0000782251730000022.GIF" wi="256" he="146" />得到q=2时均衡器向量:w=(X<sup>H</sup>DX+λI)<sup>‑1</sup>R<sub>p</sub>X<sup>H</sup>y;(6)根据q=2时均衡器输出,获得q&lt;2时的均衡器向量:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><msup><mrow><mo>(</mo><msup><mi>X</mi><mi>H</mi></msup><msub><mi>D</mi><msup><mrow><mo>|</mo><mi>y</mi><mo>|</mo></mrow><mn>2</mn></msup></msub><mi>X</mi><mo>+</mo><msub><mi>&lambda;D</mi><mi>d</mi></msub><mo>+</mo><mi>&delta;</mi><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>R</mi><mi>p</mi></msub><msup><mi>X</mi><mi>H</mi></msup><mi>y</mi></mrow>]]></math><img file="FDA0000782251730000023.GIF" wi="804" he="130" /></maths>其中:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>D</mi><mi>d</mi></msub><mo>=</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>|</mo><mi>w</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow><mo>|</mo></mrow><mrow><mn>2</mn><mo>-</mo><mi>q</mi></mrow></msup><mo>+</mo><mi>&epsiv;</mi></mrow></mfrac><mo>,</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>|</mo><mi>w</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo></mrow><mrow><mn>2</mn><mo>-</mo><mi>q</mi></mrow></msup><mo>+</mo><mi>&epsiv;</mi></mrow></mfrac><mo>,</mo><mo>...</mo><mo>,</mo><mfrac><mn>1</mn><mrow><msup><mrow><mo>|</mo><mi>w</mi><mrow><mo>(</mo><mi>L</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>|</mo></mrow><mrow><mn>2</mn><mo>-</mo><mi>q</mi></mrow></msup><mo>+</mo><mi>&epsiv;</mi></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000782251730000024.GIF" wi="1237" he="196" /></maths>diag(·)表示对角矩阵,<img file="FDA0000782251730000025.GIF" wi="92" he="90" />是由主对角元素为|y<sub>j</sub>|<sup>2</sup>组成的对角矩阵,j=1,…,N;ε为一预先设定的常数,ε>0;δ为一预先设定的机器扰动;(7)运用迭代算法获得均衡器向量的估计值<img file="FDA0000782251730000026.GIF" wi="82" he="67" />(8)根据均衡器输入输出关系,获得不同q值时的均衡器输出<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>y</mi><mo>=</mo><mi>X</mi><mover><mi>w</mi><mo>^</mo></mover><mo>.</mo></mrow>]]></math><img file="FDA0000782251730000027.GIF" wi="210" he="75" /></maths>
地址 325000 浙江省温州市茶山高教园区