CN109347562B - 一种co-ofdm***相位噪声优化补偿方法 - Google Patents
一种co-ofdm***相位噪声优化补偿方法 Download PDFInfo
- Publication number
- CN109347562B CN109347562B CN201811167870.2A CN201811167870A CN109347562B CN 109347562 B CN109347562 B CN 109347562B CN 201811167870 A CN201811167870 A CN 201811167870A CN 109347562 B CN109347562 B CN 109347562B
- Authority
- CN
- China
- Prior art keywords
- phase noise
- symbol
- matrix
- ofdm
- channel
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000005457 optimization Methods 0.000 title claims abstract description 34
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 45
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 230000009467 reduction Effects 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 50
- 239000013598 vector Substances 0.000 claims description 35
- 230000006870 function Effects 0.000 claims description 23
- 239000006185 dispersion Substances 0.000 claims description 22
- 230000003595 spectral effect Effects 0.000 claims description 20
- 230000004044 response Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 12
- 125000004122 cyclic group Chemical group 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 230000001427 coherent effect Effects 0.000 claims description 10
- 239000000835 fiber Substances 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 230000009977 dual effect Effects 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 8
- IAZDPXIOMUYVGZ-UHFFFAOYSA-N Dimethylsulphoxide Chemical compound CS(C)=O IAZDPXIOMUYVGZ-UHFFFAOYSA-N 0.000 claims description 6
- 238000001514 detection method Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 6
- 239000000969 carrier Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 description 30
- 230000003287 optical effect Effects 0.000 description 27
- 239000013307 optical fiber Substances 0.000 description 19
- 238000012937 correction Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 2
- 230000006735 deficit Effects 0.000 description 2
- 230000005764 inhibitory process Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 229910052691 Erbium Inorganic materials 0.000 description 1
- HSRJKNPTNIJEKV-UHFFFAOYSA-N Guaifenesin Chemical compound COC1=CC=CC=C1OCC(O)CO HSRJKNPTNIJEKV-UHFFFAOYSA-N 0.000 description 1
- 229930003270 Vitamin B Natural products 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000011719 vitamin A Substances 0.000 description 1
- 235000019156 vitamin B Nutrition 0.000 description 1
- 239000011720 vitamin B Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/25—Arrangements specific to fibre transmission
- H04B10/2507—Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion
- H04B10/2513—Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to chromatic dispersion
- H04B10/2525—Arrangements specific to fibre transmission for the reduction or elimination of distortion or dispersion due to chromatic dispersion using dispersion-compensating fibres
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03159—Arrangements for removing intersymbol interference operating in the frequency domain
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Optical Communication System (AREA)
Abstract
一种CO‑OFDM***相位噪声优化补偿方法,通过减少相位噪声谱结构的约束条件,不管其信道估计精确与否,其优化算法均可有效补偿信道均衡后的残余幅度噪声,因此最终的相位噪声估计精度明显优于其他相位噪声估计算法,并取得了突破性的相位噪声补偿效果。本发明中运用了相位噪声降维模型,对于CO‑OFDM***来说该相位噪声优化方法计算复杂度完全可以容忍。本发明该发明能极大促进CO‑OFDM***在长距离接入网和城域网中的应用。
Description
技术领域
本发明属于光通信网络技术领域,特别涉及一种CO-OFDM***的相位噪声优化补偿方法。
背景技术
在相干光通信***中结合电域正交频分复用(OFDM)调制技术,形成相干光正交频分复用(Coherent Optical Orthogonal Frequency Division Multiplexing,CO-OFDM)传输技术,具有对光纤色散和偏振模色散具有良好的抑制作用、用数字信号处理灵活地补偿***损伤的能力、高频谱利用率等优点,已成为长距离高速通信***和光接入网等领域备受关注的技术之一。
CO-OFDM***结构如图1所示,按其功能可以分为5个模块:CO-OFDM***发射端模块101、光调制模块102、光纤传输模块103、光电检测模块104以及CO-OFDM***接收端模块105,CO-OFDM发射端模块产生的电域信号经过电光调制的上变频变成光域的CO-OFDM信号,CO-OFDM信号经光纤传输、平衡探测器后经光电转换成电域的信号,CO-OFDM接收端再对接收到的电信号进行信号处理以期恢复原始的发送段数据。结合图1,对整个***的工作过程进行详细表述。CO-OFDM***串行输入的数据106经过串并转换模块107,变为并行的N路数据;按照不同的调制格式将串并转换后的信号进行数字调制108;快速傅里叶逆变换IFFT模块109实现信号从频域到时域的转换;加入循环前缀CP110;将得到的电域信号进行并串转换111。上述信号的同相分量和正交分量信号分别通过数模转换器112、113变换为模拟信号并通过低通滤波器114、115;采用放大器将信号的同相分量116和正交分量117放大并注入到I/Q调制器中实现同相分量I和正交分量Q对光信号的正交调制;I/Q调制器由3个双臂的马赫增德尔MZM调制器120、121和122组成,其中两个调制器实现对信号的调制,第三个调制器122控制光调制的同相分量I和正交分量Q的相位差;分别调节两个调制器120、121的直流偏置保证实现信号调制的调制器工作在最小功率点,而第三个控制相位差的调制器工作在正交点以保证两路信号存在90°的相位差;118表示CO-OFDM***的发射激光器,通过分路器119分成两束同样的激光,用于驱动二个光调制器120和121。二个光调制器输出的信号通过合束器123,变成单路的光信号,接着输入到光纤信道进行传输。产生的CO-OFDM信号在光纤124中经过长距离的传输后,经过直接的光-光放大器-掺铒光纤放大器(EDFA)125补偿光纤损耗后再进行传输,表示长距离的光纤,126表示光带通滤波器。经过长距离的光纤传输后,光电检测模块将光域信号变换成电域的信号。127表示CO-OFDM***接收端的本地激光器,通过分路器分成两束同样的激光,128表示一个90°的相移器;129和130表示两个耦合器,驱动4个光电二极管(PD)131、132、133和134。135和136表示两个减法器,分别对应输出接收信号的同相分量I和正交分量Q。得到的同相分量I和正交分量Q经过低通滤波器137、138和模数转换器139、140转换后进入CO-OFDM接收端。CO-OFDM接收端进行数字信号处理141,进行CO-OFDM发送端的逆过程,进行串并转换142,移除循环前缀CP143,然后进行FFT变换144,对CO-OFDM信号进行数字解调145,最后经过并串转换146恢复得到原始的发送端串行数据输出147。
作为多载波传输***,CO-OFDM与单载波相干光通信***相比(如相干光正交幅度调制(QAM)***),由于其符号周期长、子载波间隔小,对激光器相位噪声和光纤的非线性更敏感,所以CO-OFDM传输技术的发展高度依赖于激光器相位噪声和光纤非线性损伤的补偿精度。目前机器学习已被提出用于缓和CO-OFDM***中非线性的影响。虽然在CO-OFDM***中非线性损伤和线性损伤可以通过神经网络等非线性均衡器同时被补偿,但是这样的方法是用大量的训练符号和巨大的计算复杂度作为代价的。本发明中,我们致力于在大激光器线宽和高阶调制的CO-OFDM***中实现高精度的相位噪声补偿方法,使其能应用于一些使用大线宽、低成本激光器的城域接入网。
在CO-OFDM***中有许多用于相位噪声补偿的数字信号处理算法。传统的算法可以分为三种类型:盲相位噪声补偿算法(文献1,Son T.Le,Paul A.Haigh,Andrew D.Ellis,Sergei K.Turitsyn.Blind phase noise compensation for CO-OFDMtransmission.Journal of Lightwave Technology,2016,34(2):745-753.即Son T.Le,Paul A.Haigh,Andrew D.Ellis,Sergei K.Turitsyn.CO-OFDM传输***的盲相位噪声补偿,光波技术学报,2016,34(2):745-753.);判决反馈算法(文献2,Xuezhi Hong,XiaojianHong,Sailing He.Linearly interpolated sub-symbol optical phase noisesuppression in CO-OFDM system.Optics Express,2015,23(4):4691-4702.即XuezhiHong,Xiaojian Hong,Sailing He,CO-OFDM***中线性插值与子符号光相位噪声抑制,光学快报,2015,23(4):4691-4702.);基于导频的算法,包括运用分散的导频结构和射频导频。当CO-OFDM***采用大线宽激光器和高阶调制方式时,利用分散的导频子载波是现实应用中最通用和实用的相位噪声估计方法。在众多基于导频的算法中,线性插值与子符号光相位噪声补偿(LI-SCPEC)是最具代表性的算法之一(文献2,Xuezhi Hong,Xiaojian Hong,Sailing He.Linearly interpolated sub-symbol optical phase noise suppressionin CO-OFDM system.Optics Express,2015,23(4):4691-4702.即Xuezhi Hong,XiaojianHong,Sailing He,CO-OFDM***中线性插值与子符号光相位噪声抑制,光学快报,2015,23(4):4691-4702.),该算法先用最小二乘法估计公共相位噪声(common phase error,CPE):
φn=angle(Yn/Xn),i=0,1,...,Ns
其中Xn、Yn为第n个发送、接收符号。
再线性插值得到各个子载波上的相位噪声,Ncp为循环前缀的长度:
随后进行相位纠正。经过信道均衡和相位纠正的信号为:
再对上述处理过的信号进行判决操作Q(·),并做IFFT变换到时域:
bn=[bn(il+0),bn(il+1),...,bn(il+l-1)]T
an=[an(il+0),an(il+1),...,an(il+l-1)]T
这里i∈[0,NB-1],l=N/NB表示每个子符号的长度。
此时相位噪声矩阵这种方法(linearly interpolated sub-symbol optical PNC scheme,LI-SCPEC)通过以上步骤实现相位噪声估计,但其相位噪声补偿性能在***激光器线宽大和调制阶数高时效果不佳。
最近,有学者在无线OFDM***中,把相位噪声谱结构与降维模型结合,将相位噪声估计问题作为优化问题提出,在相位噪声估计精度上有很大突破(文献3,PramodMathecken,Taneli Riihonen,Stefan Werner,Risto Wichman.Phase noise estimationin OFDM:utilizing its associated spectral geometry.IEEE Transactions onSignal Processing,2016,64(8):1999-2012.即Pramod Mathecken,Taneli Riihonen,Stefan Werner,Risto Wichman,OFDM***中的相位噪声补偿:利用相位噪声谱结构,IEEE信号处理汇刊,2016,64(8):1999-2012.文献4,Pramod Mathecken,Taneli Riihonen,Stefan Werner,Risto Wichman.Constrained phase noise estimation in OFDM usingscattered pilots without decision feedback.IEEE Transactions on SignalProcessing,2017,65(9):2348-2362.即Pramod Mathecken,Taneli Riihonen,StefanWerner,Risto Wichman,OFDM***中无判决反馈的利用分散导频及带约束条件的相位噪声估计,IEEE信号处理汇刊,2017,65(9):2348-2362.)。文献3将相位噪声频谱结构特性(phase noise spectral geometry,PSNG)进行了详细的描述,即每个OFDM符号的相位噪声谱向量δn满足:
Λl是Kronecker delta函数(Λ0=1,Λl=0,l=1,2,...,Nc-1),Pl=(P1)l为置换矩阵由Nc×Nc维的P1构成,P1的第一列是Nc×1的向量[0,1,0,...,0]T,第j列由它向下循环移位j-1次得到。l=0时,P0为Nc×Nc的单位矩阵。使用降维模型可将Nc维相位噪声谱向量δn降至N维γn,故PNSG变为相应的N维等式:该等式可拆写作两个约束条件:和文献4介绍了两种相位噪声优化算法,一种为不带PNSG约束的ULS估计,另一种为带PNSG两个约束条件的GLS估计。本发明中定义矩阵K为K×Nc矩阵,每一列是都单位向量,元素1的位置对应于在发送端记录导频的位置,故发送导频符号为:
Xn,p=KXn
这里是发送端频域符号。因此接收端导频符号的估计值为信道频响矩阵Hn是对角元素为的对角矩阵,矩阵Υn是按列循环矩阵,它的第一列是对应符号的接收数据随后第j列用第一列向下循环位移j-1次,是δn的估计值。
则上述ULS算法的相位噪声估计代价函数为:
GLS算法的代价函数则为
受这种基于导频的相位噪声优化方法启发,本发明将提出一种基于导频的相位噪声优化补偿算法用于CO-OFDM***在大激光器线宽和高阶调制时的高精度相位噪声补偿。在CO-OFDM***中,相位噪声优化补偿算法需要对其信道进行估计和均衡,故信道估计一般采用传统最小二乘估计(least square channel equalization,LS-CE)。即:
其中即为第n个符号,第k个子载波的信道频率响应,Xn,k,Yn,k分别为发送和接收的第n个符号,第k个子载波上的数据,V′n,k为噪声项。LS估计算法中,不考虑噪声项的影响。若想要提高信道估计精度,加入标量卡尔曼滤波减小噪声和载波间干扰(intercarrier interference,ICI)的影响可略微改善算法性能。具体实施如下:
状态预测方程
最小均方误差(MSE)预测方程
卡尔曼增益方程
K(n)=P(n|n-1)/(1+P(n|n-1))
状态更新方程
最小均方误差更新方程
P(n|n)=(1-K(n))P(n|n-1)
这里Kgain,Ka和Kb都是常数。
发明内容
本发明的主要目的在于克服现有算法在大线宽、长距离CO-OFDM***中相位噪声估计精确度不高导致相位噪声补偿效果较差的问题,提出一种相位噪声优化补偿方法(modified geometry constrained least square estimation,M-GLS)用于大线宽CO-OFDM***的相位噪声高精度补偿。
为了解决上述技术问题,本发明提供如下的技术方案:
一种CO-OFDM***相位噪声优化补偿方法,包括以下步骤:
(1)接收端初始信号处理,过程如下:
1-1、接收端对接收到的CO-OFDM信号进行相干探测接收,然后进行模数转换,得到电域的信号;随后进行电域光纤色散补偿,将光纤信道频域传递函数的解析形式经傅立叶变换到时域,设计时域有限长单位冲激响应(FIR)滤波器来实现,该滤波器的阶数随色散累积而增加;
1-2、移除循环前缀CP;
1-3、采用快速傅里叶变换(FFT)将信号从时域变为频域;
(2)最小二乘信道估计(LS-CE)与均衡,过程如下:
2-1、LS粗略信道估计。在发送端,设置前Np个符号为训练符号,在接收端,根据LS得:
这里t为参与信道估计的相邻子载波信道数;
2-3、信道均衡,在每个OFDM帧中,对接收端的Ns个OFDM数据符号进行信道均衡后,第n个OFDM符号第k个频域数据Y′n,k为,
(3)M-GLS相位噪声估计:通过减少相位噪声谱结构的约束条件,不管其信道估计精确与否,其优化算法均可有效补偿信道均衡后的残余幅度噪声,通过构造代价函数和选择约束条件,写出对偶问题最后求出相位噪声精确估计值;
(4)最终相位噪声补偿:对其频域信号进行如下的相位噪声补偿,
进一步,所述步骤(3)包括以下步骤:
3-1、相位噪声谱向量降维,假设CO-OFDM***中,在接收端其他所需的数字信号处理算法都没有误差,则接收的时域第n个OFDM符号、第m个采样点yn,m表示为:
这里xn,m,hn,m和vn,m分别为时域发送信号,信道冲击响应和噪声项‘’
由于在CO-OFDM***中相位噪声可被认为是低通信号,因此仅需估计位于相位噪声谱向量的顶部和底部的元素即低频成分。故通过实现使相位噪声结构特性保持变换(PNSG preserving transformations,PPT)矩阵,即可实现降维,使δn的维度由Nc降为N,计算复杂度也由此大大降低。
相位噪声频域的结构特性PNSG如下所示:
这里是列向量,Λl是Kronecker delta函数(Λ0=1,Λl=0,l=1,2,...,Nc-1),Pl=(P1)l为置换矩阵由Nc×Nc维的P1构成,P1的第一列是Nc×1的向量[0,1,0,...,0]T,第j列由它向下循环移位j-1次得到。l=0时,P0为Nc×Nc的单位矩阵。式中表示共轭转置。T是Nc×N的PPT矩阵,实现相位噪声谱向量的降维,满足δn=Tγn。降维后的相位噪声谱由N维向量γn表示,满足N维频谱结构特性:
上述特性包括l=0和l≠0两个条件:
Xn,p=KXn
这里是发送端频域符号。根据***模型可知,接收端导频符号的估计值则为信道频响矩阵Hn是对角元素为的对角矩阵,矩阵Υn是按列循环矩阵,它的第一列是对应符号的接收数据随后第j列用第一列向下循环位移j-1次,为δn的估计值;
于是相位噪声估计的代价函数定义为:
3-3、M-GLS的算法及其约束条件
于是M-GLS的优化问题为以下带有PNSG部分约束条件的形式:
3-4、写出对偶问题,利用S-procedure写出对应的对偶问题:
(D):Maximizeτ
这里τ,αl,βl都是优化变量;
3-5、求其最优解为:
其中(·)+定义为伪逆矩阵,αl ◇,βl ◇是通过解(D)获得的最小值。
本发明的技术构思为:相位噪声优化补偿方法在发射端每个OFDM帧***若干训练符号,在每个OFDM数据符号中按一定间隔***若干导频子载波作为开销。该相位噪声优化补偿方法在接收端数据进行去循环前缀、FFT、色散补偿、卡尔曼滤波信道估计与均衡的操作后,首先利用相位噪声谱向量的低通特性实现降维,极大地减小计算复杂度。之后根据发送与接收的导频符号间最小二乘误差最小的原则,写出相位噪声优化估计的代价函数。此时,以相位噪声结构特性的部分条件作为约束,再利用S-procedure写出该优化问题的对偶问题,在MATLAB中用CVX工具箱求得最优解。最后将此最优值进行相位噪声补偿。该方法较导频算法LI-SCPEC及其他优化算法相比,取得了较好的补偿效果。在40Gb/s、传输距离为400km的CO-OFDM***中对本发明的算法进行了验证,调制方式为圆的16阶正交振幅调制(C-16QAM)、联合激光器线宽为1.8MHz时,其误码率性能才达硬判决前向纠错上限(HD-FEC),该发明能极大促进CO-OFDM***在长距离接入网和城域网中的应用。
本发明与现有技术相比,具有如下优点和有益效果:
1.对高阶数字调制和大线宽激光器的CO-OFDM***,本发明相位噪声估计的修正办法获得了较好的PNC效果,如对C-16QAM、400km的M-GLS-ICE,联合激光器线宽上限可达1.8MHz。本发明所用频域导频子载波间隔较大,比同类相关算法使用更少的导频子载波,较大提高了***的频谱利用率。
2.本发明提出的M-GLS与传统PNC方法相比,虽然其复杂度略有增加,但以此为代价获得的是误码率性能的显著提高,在大线宽激光器、高阶调制和长距离传输的CO-OFDM***中,该方法补偿效果具有很大的优越性,如对16QAM、0.6MHz联合激光器线宽的M-GLS-ICE,传输距离可达600km,适用于实际中长距离的光纤通信网络。
附图说明
图1是现有技术中的CO-OFDM***的示意图。
图2是本发明提出的相位噪声优化补偿算法M-GLS的原理图。
图3是本发明实施例1中在C-16QAM、传输距离为0(背靠背***)、信道估计采用传统的最小二乘估计(LS-CE)时,ULS(不带约束的优化算法)、GLS(带全部PNSG约束条件的优化算法)、M-GLS、LI-SCPEC的误码率性能随联合激光器线宽变化时的关系曲线。
图4是本发明实施例1中在C-16QAM、传输距离为400km、联合激光器线宽为0.6MHz时,M-GLS-ICE算法误码率性能随卡尔曼滤波中协方差初值与参数Kb变化时的三维曲线。
图5是本发明实施例1中在C-16QAM、传输距离为0、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随联合激光器线宽变化时的关系曲线。
图6是本发明实施例1中在C-16QAM、传输距离为400km、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随联合激光器线宽变化时的关系曲线。
图7是本发明实施例1中在C-16QAM、传输距离为400km、联合激光器线宽为1.4MHz时,ULS-ICE、GLS-ICE、M-GLS-ICE、LI-SCPEC-ICE四种相位噪声补偿算法得到的星座图。
图8是本发明实施例1中在C-16QAM、联合激光器线宽为0.6MHz、光信噪比为32dB、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随传输距离变化时的关系曲线,同时也显示了联合激光器线宽为0MHz,即不需要PNC时误码率性能随传输距离变化的曲线,此为理论上算法性能的上限。
图9是本发明实施例1中在C-16QAM、传输距离为400km、联合激光器线宽为0.6MHz、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随光信噪比(OSNR)变化时的关系曲线,同时也显示了联合激光器线宽为0MHz,即不需要PNC时误码率性能随OSNR变化的曲线,此为理论上算法性能的上限。
具体实施方式
下面结合实施例及附图对本发明作进一步详细地描述,但本发明的实施方式不限于此。
参照图2~图9,一种CO-OFDM***相位噪声优化补偿方法,主要涉及相干光正交频分复用CO-OFDM***接收端的信号处理问题,参考背景技术中对CO-OFDM***结构的详细描述。
如图1所示,CO-OFDM***包括CO-OFDM***发射端模块101、CO-OFDM光调制模块102、光纤传输模块103、光电检测模块104以及CO-OFDM***接收端模块105,***发射端产生的信号经过了光调制的上变频变成光域的CO-OFDM信号,CO-OFDM信号经光纤传输、平衡探测器后经光电转换成电域的信号,***接收端再对接收到的电域信号进行信号处理以期恢复原始的发送端数据。初始40Gb/s伪随机码二进制数据流用高阶QAM调制(16QAM)映射到512个子载波上,FFT或者IFFT的点数为1024。每个OFDM数据符号中的循环前缀CP长度为128点。每50km单模光纤后接一个掺铒光纤放大器EDFA,该放大器增益为10dB,噪声系数为4dB。整个光纤链路共有8段50km单模光纤加放大器EDFA构成(传输距离为400km时)。该单模光纤的色散系数为17ps/nm·km,色散斜率为0.075ps/(nm2·km),非线性系数为1.3W-1·km-1,PMD系数为损耗系数为0.2dB/km。OFDM调制前先对二进制伪随机码进行C-16QAM映射。发射端激光器与相干接收端激光器具有相同的线宽和波长,其波长为1550nm。激光器最优发射功率为0dBm。光带通滤波器的带宽为40GHz。将每个OFDM帧首先4个OFDM符号为训练符号,每个OFDM符号导频序列间隔为16。LI-SCPEC方法中每个OFDM符号划分亚符号数目为NB=4。
下面结合图2,对本发明的一种CO-OFDM***相位噪声优化补偿方法步骤进行详细说明,包括以下步骤:
S201:接收端初始信号处理,过程如下:
S201-1、接收端对接收到的CO-OFDM信号进行相干探测接收,然后进行模数转换,得到电域的信号。随后进行电域光纤色散补偿,将光纤信道频域传递函数的解析形式经傅立叶变换到时域,设计时域有限长单位冲激响应(FIR)滤波器来实现,该滤波器的阶数随色散累积而增加。
S201-2、移除循环前缀CP。
S201-3、采用快速傅里叶变换(FFT)将信号从时域变为频域。
S202:最小二乘信道估计(LS-CE)。假定一个OFDM帧在时域包含Ns=128个OFDM符号,前Np=4个为训练符号,每个OFDM符号在频域包含Nc=512个子载波,过程如下:
S202-1、LS粗略信道估计。在发送端,设置前Np个符号为训练符号,在接收端,根据LS得:
这里t为参与信道估计的相邻子载波信道数。
S203、信道均衡。在每个OFDM帧中,对接收端的Ns个OFDM数据符号进行信道均衡后,第n个OFDM符号第k个频域数据Y′n,k为,
S204:M-GLS相位噪声优化估计与补偿,过程如下:
S204-1、相位噪声谱向量降维。假设CO-OFDM***中,在接收端其他所需的数字信号处理算法都没有误差,则接收的时域第n个OFDM符号、第m个采样点yn,m可以表示为:
这里xn,m,hn,m和vn,m分别为时域发送信号,信道冲击响应和噪声项。
由于在CO-OFDM***中相位噪声可被认为是低通信号,因此仅需估计位于相位噪声谱向量的顶部和底部的元素即低频成分。故通过实现使相位噪声结构特性保持变换(PNSG preserving transformations,PPT)矩阵,即可实现降维,使δn的维度由Nc降为N,计算复杂度也由此大大降低。
相位噪声频域的结构特性PNSG如下所示:
这里是列向量,Λl是Kronecker delta函数(Λ0=1,Λl=0,l=1,2,...,Nc-1),Pl=(P1)l为置换矩阵由Nc×Nc维的P1构成,P1的第一列是Nc×1的向量[0,1,0,...,0]T,第j列由它向下循环移位j-1次得到。l=0时,P0为Nc×Nc的单位矩阵。式中表示共轭转置。T是Nc×N的PPT矩阵,实现相位噪声谱向量的降维,满足δn=Tγn。降维后的相位噪声谱由N维向量γn表示,满足N维频谱结构特性:
上述特性包括l=0和l≠0两个条件:
这里PPT取如下形式:
Xn,p=KXn
这里是发送端频域符号。因此接收端导频符号的估计值为信道频响矩阵Hn是对角元素为的对角矩阵,矩阵Υn是按列循环矩阵,它的第一列是对应符号的接收数据随后第j列用第一列向下循环位移j-1次,是δn的估计值。
于是相位噪声估计的代价函数定义为:
S204-3、M-GLS约束条件。
M-GLS的优化问题为以下带有PNSG部分约束条件的形式:
S204-4、对偶问题。利用S-Procedure写出其对偶问题,如下所示:
(D):Maximizeτ
这里τ,αl,βl都是优化变量。
S204-5、求其最优解为:
其中(·)+定义为伪逆矩阵,αl ◇,βl ◇是通过解(D)获得的最小值。
S205:对其频域信号进行如下的相位噪声补偿,
对该发明提出的相位噪声优化补偿方法进行仿真数值验证。
图3显示了在C-16QAM、传输距离为0(背靠背***)、信道估计采用LS-CE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随联合激光器线宽变化时的关系曲线。此时导频间隔为16,降维维度优化为N=16。由图可知,当联合激光器线宽在0.8MHz~3.2MHz时,M-GLS的优越性一目了然,在3.2MHz时M-GLS才达到硬判决前向纠错上限(HD-FEC),即BER=3.8e-03,且在0.8MHz以下时,用直接对比错误比特法计算M-GLS的误码率显示为0。而带有全部PNSG约束的GLS在线宽不断增大的过程中其性能远差于ULS和M-GLS,且线宽越大相差越多,在1.4MHz时即达到了纠错上限。
图4显示了卡尔曼滤波信道估计(ICE)中,在固定Kgain和Ka都取1时,协方差初值P(0|0)和Kb两个参数的取值对误码性能的影响。Z轴为1/BER,三维曲线图的最高点对应的参数取值即为最优。
激光器最优发射功率0dBm时,图5和6显示了在C-16QAM、传输距离为0(背靠背***)和400km、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随联合激光器线宽变化时的关系曲线。由图可以看出,传输距离为0或400km曲线呈现出的规律大体一致,图3的说明已分析过背靠背***的情况,这里主要分析400km的长距离传输。同样M-GLS-ICE的性能优于其他所有算法,在联合激光器线宽为1.8MHz时才达到HD-FEC上限。对于所有的相位噪声补偿算法,运用ICE比LS-CE可以获得更好地性能,尤其是GLS的性能得到了明显的改善,这说明了在CO-OFDM***中,带有PNSG全部约束的GLS对信道估计的精度十分敏感。对于ULS、M-GLS和LI-SCPEC,运用ICE后的性能相比原来的LS-CE有轻微的改善,表示它们对信道估计的精度都不是特别敏感。
图7显示了在激光器最优发射功率为0dBm以及传输距离为400km、联合激光器线宽为1.4MHz时,ULS-ICE、GLS-ICE、M-GLS-ICE、LI-SCPEC-ICE四种相位噪声补偿算法得到的接收端C-16QAM星座图。接收端OFDM解调后的16QAM信号点受到激光器相位噪声及光纤色散的影响,发生了严重的旋转和发散。在接收端首先根据光纤信道的时域特性设计有限长单位冲激响应(FIR)滤波进行电域色散补偿,接下来进行相位噪声补偿。图7显示的信号点已经均衡为16块数据点,由星座图的发散程度即可看出ICI相位噪声是否得到了有效地抑制,从而降低预判决中判决错误发生的概率。显然M-GLS-ICE的抑制效果最佳。
为了探索本发明提出的M-GLS在传输距离上可达到的上限,图8显示了在C-16QAM、联合激光器线宽为0.6MHz、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随传输距离变化时的关系曲线。随着传输距离的增加,在CO-OFDM***中光纤的非线性损伤变得越来越显著,但从图中可以看出,M-GLS-ICE可以较好地抑制光纤非线性的影响,在达到HD-FEC上限时可比导频算法LI-SCPEC-ICE多传输300km,适用于长距离的城域接入网。
图9对发送端光信噪比(OSNR)对***性能的影响进行了讨论,显示了在C-16QAM、传输距离为400km、联合激光器线宽为0.6MHz、分别采用LS-CE和ICE时,ULS、GLS、M-GLS、LI-SCPEC的误码率性能随OSNR变化时的关系曲线。由图9可以看出,本发明提出的M-GLS即使在低OSNR的情况下依然优于其他所有算法。
本发明中优化算法的降维维度N=16,算法复杂度为O(N4.5),因此在复杂度可以容忍的情况下,本发明提供了一种适用于大线宽、长距离、高阶调制的CO-OFDM***,且PNC性能远远优于传统基于导频算法的算法。
以上对本发明所述的相干光正交频分复用CO-OFDM***中的相位噪声优化补偿算法进行了详细地的介绍,以上的实例说明只是用于帮助理解本发明的方法及其核心思想而非对其进行限制,其他的任何未背离本发明的精神实质与原理下所作改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (2)
1.一种CO-OFDM***相位噪声优化补偿方法,其特征在于,所述方法包括以下步骤:
(1)接收端初始信号处理,过程如下:
1-1、接收端对接收到的CO-OFDM信号进行相干探测接收,然后进行模数转换,得到电域的信号;随后进行电域光纤色散补偿,将光纤信道频域传递函数的解析形式经傅立叶变换到时域,设计时域有限长单位冲激响应(FIR)滤波器来实现,该滤波器的阶数随色散累积而增加;
1-2、移除循环前缀CP;
1-3、采用快速傅里叶变换将信号从时域变为频域;
(2)最小二乘信道估计与均衡,过程如下:
2-1、LS粗略信道估计,在发送端,设置前Np个符号为训练符号,在接收端,根据LS得:
这里t为参与信道估计的相邻子载波信道数;
2-3、信道均衡,在每个OFDM帧中,对接收端的Ns个OFDM数据符号进行信道均衡后,第n个OFDM符号第k个频域数据Y′n,k为,
(3)M-GLS相位噪声估计:通过减少相位噪声谱结构的约束条件,不管其信道估计精确与否,其优化算法均可有效补偿信道均衡后的残余幅度噪声,通过构造代价函数和选择约束条件,写出对偶问题最后求出相位噪声精确估计值;
(4)最终相位噪声补偿:对其频域信号进行如下的相位噪声补偿,
2.如权利要求1所述的一种CO-OFDM***相位噪声优化补偿方法,其特征在于,所述步骤(3)包括以下步骤:
3-1、相位噪声谱向量降维,假设CO-OFDM***中,在接收端其他所需的数字信号处理算法都没有误差,则接收的时域第n个OFDM符号、第m个采样点yn,m表示为:
这里xn,m,hn,m和vn,m分别为时域发送信号,信道冲击响应和噪声项;
相位噪声频域的结构特性PNSG如下所示:
这里是列向量,Λl是Kronecker delta函数(Λ0=1,Λl=0,l=1,2,...,Nc-1),Pl=(P1)l为置换矩阵由Nc×Nc维的P1构成,P1的第一列是Nc×1的向量[0,1,0,...,0]T,第j列由它向下循环移位j-1次得到;l=0时,P0为Nc×Nc的单位矩阵,式中表示共轭转置;T是Nc×N的PPT矩阵,实现相位噪声谱向量的降维,满足δn=Tγn,降维后的相位噪声谱由N维向量γn表示,满足N维频谱结构特性:
上述特性包括l=0和l≠0两个条件:
Xn,p=KXn
这里是发送端频域符号,根据***模型可知,接收端导频符号的估计值则为信道频响矩阵Hn是对角元素为的对角矩阵,矩阵Υn是按列循环矩阵,它的第一列是对应符号的接收数据随后第j列用第一列向下循环位移j-1次,为δn的估计值;
于是相位噪声估计的代价函数定义为:
3-3、M-GLS的算法及其约束条件
于是M-GLS的优化问题为以下带有PNSG部分约束条件的形式:
3-4、写出对偶问题,利用S-procedure写出对应的对偶问题:
(D):Maximizeτ
这里τ,αl,βl都是优化变量;
3-5、求其最优解为:
其中(·)+定义为伪逆矩阵,αl ◇,βl ◇是通过解(D)获得的最小值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811167870.2A CN109347562B (zh) | 2018-10-08 | 2018-10-08 | 一种co-ofdm***相位噪声优化补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811167870.2A CN109347562B (zh) | 2018-10-08 | 2018-10-08 | 一种co-ofdm***相位噪声优化补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109347562A CN109347562A (zh) | 2019-02-15 |
CN109347562B true CN109347562B (zh) | 2020-05-05 |
Family
ID=65308372
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811167870.2A Active CN109347562B (zh) | 2018-10-08 | 2018-10-08 | 一种co-ofdm***相位噪声优化补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109347562B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110190906B (zh) * | 2019-04-15 | 2020-08-04 | 浙江工业大学 | 一种针对高阶qam相干光***的非线性均衡方法 |
CN110011734B (zh) * | 2019-04-18 | 2020-06-30 | 杭州电子科技大学 | 基于导频和二维投影直方图的co-ofdm***中的cpe补偿方法 |
CN110233812B (zh) * | 2019-04-28 | 2022-05-24 | 国网安徽省电力有限公司信息通信分公司 | Co-ofdm***及其接收端和非线性均衡方法、装置 |
CN110830398B (zh) * | 2019-11-13 | 2020-09-25 | 湖南师范大学 | 应用于光纤dmt***中的符号内频域平均信道估计方法 |
CN111917487B (zh) * | 2020-08-12 | 2021-08-03 | 中国科学院长春光学精密机械与物理研究所 | 激光通信干扰消除方法、装置及计算机可读存储介质 |
CN112769497B (zh) * | 2021-02-23 | 2022-02-22 | 苏州大学 | 对大容量高阶qam相干光进行非线性补偿的方法 |
CN113114602B (zh) * | 2021-04-09 | 2022-06-21 | 重庆邮电大学 | Co-ofdm***中基于cp插值与ckf的相位噪声补偿方法 |
CN114095322B (zh) * | 2021-09-28 | 2023-07-14 | 超讯通信股份有限公司 | 相位噪声抑制方法、装置及存储介质 |
CN114465864B (zh) * | 2022-02-17 | 2024-01-12 | 哲库科技(北京)有限公司 | 相位噪声补偿方法、装置、终端、存储介质及程序产品 |
CN115778319B (zh) * | 2022-11-09 | 2024-05-07 | 山东大学 | 基于双光谱仪进行光源噪声补偿的可见光oct*** |
CN117834034B (zh) * | 2024-03-04 | 2024-05-14 | 深圳市光为光通信科技有限公司 | 一种用于100g工业级单波光模块的自适应调制与解调方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8760201B1 (en) * | 2013-03-11 | 2014-06-24 | Analog Devices Technology | Digitally programmed capacitance multiplication with one charge pump |
CN105187345A (zh) * | 2015-09-23 | 2015-12-23 | 浙江工业大学 | 一种适用于co-ofdm***的无迭代盲相位噪声补偿方法 |
CN105187344A (zh) * | 2015-09-23 | 2015-12-23 | 浙江工业大学 | 适用于co-ofdm***的判决辅助时域平均近似ici相位噪声补偿方法 |
CN107171735A (zh) * | 2017-05-16 | 2017-09-15 | 浙江工业大学 | 一种时频域卡尔曼滤波大线宽co‑ofdm相位噪声补偿方法 |
CN107395282A (zh) * | 2017-07-20 | 2017-11-24 | 浙江工业大学 | 一种时域无迹卡尔曼滤波的大线宽co‑ofdm***相位噪声补偿方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104601239B (zh) * | 2015-01-12 | 2017-05-17 | 西南交通大学 | 一种基于强度噪声方差以及低通滤波器的光纤自适应非线性补偿方法 |
-
2018
- 2018-10-08 CN CN201811167870.2A patent/CN109347562B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8760201B1 (en) * | 2013-03-11 | 2014-06-24 | Analog Devices Technology | Digitally programmed capacitance multiplication with one charge pump |
CN105187345A (zh) * | 2015-09-23 | 2015-12-23 | 浙江工业大学 | 一种适用于co-ofdm***的无迭代盲相位噪声补偿方法 |
CN105187344A (zh) * | 2015-09-23 | 2015-12-23 | 浙江工业大学 | 适用于co-ofdm***的判决辅助时域平均近似ici相位噪声补偿方法 |
CN107171735A (zh) * | 2017-05-16 | 2017-09-15 | 浙江工业大学 | 一种时频域卡尔曼滤波大线宽co‑ofdm相位噪声补偿方法 |
CN107395282A (zh) * | 2017-07-20 | 2017-11-24 | 浙江工业大学 | 一种时域无迹卡尔曼滤波的大线宽co‑ofdm***相位噪声补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109347562A (zh) | 2019-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109347562B (zh) | 一种co-ofdm***相位噪声优化补偿方法 | |
CN111010239B (zh) | 一种相干光纤通信***中非线性相位噪声补偿方法及*** | |
JP4872003B2 (ja) | 光ofdm受信器および光伝送システムおよびサブキャリア分離回路およびサブキャリア分離方法 | |
CN107171735B (zh) | 一种时频域卡尔曼滤波大线宽co-ofdm相位噪声补偿方法 | |
US8218979B2 (en) | System, method and apparatus for coherent optical OFDM | |
Spinnler | Equalizer design and complexity for digital coherent receivers | |
US8913901B2 (en) | System and method for blind equalization and carrier phase recovery in a quadrature amplitude modulated system | |
US20050196176A1 (en) | Equalization strategy for dual-polarization optical transport system | |
JP5390607B2 (ja) | 予等化伝送システム | |
CN113300779B (zh) | 一种导频辅助的co-fbmc/oqam***相位噪声补偿方法 | |
CN111865432B (zh) | 一种co-fbmc/oqam***中的盲相位噪声补偿方法 | |
Nhan et al. | Sparse preamble design for polarization division multiplexed CO-OFDM/OQAM channel estimation | |
US8989324B2 (en) | Receiver and receiving method | |
Zhu et al. | Complexity reduction with a simplified MIMO volterra filter for PDM-twin-SSB PAM-4 transmission | |
US8615171B2 (en) | Computationally-efficient equalization in polmux OFDM transmission with direct detection via joint transmitter and receiver processing | |
Zhu et al. | Nyquist-WDM with low-complexity joint matched filtering and adaptive equalization | |
WO2014060031A1 (en) | Method and apparatus for estimating channel coefficients of a mimo communications channel | |
Hamaoka et al. | Super high density multi-carrier transmission system by MIMO processing | |
Zhang et al. | Cost-effective and hardware-efficient coherent scheme for 400G short-reach transmission | |
Cao et al. | Decision-aided joint compensation of channel distortion and transmitter IQ imbalance for coherent optical OFDM | |
Ishihara et al. | Frequency-domain equalization for coherent optical single-carrier transmission systems | |
Fatadin et al. | DSP techniques for 16-QAM coherent optical systems | |
Sato et al. | Mitigation of Inter-Subcarrier Linear Crosstalk with MIMO equalizer | |
Ahmad et al. | Normalized LMS filter for chromatic dispersion compensation in single polarization 40 Gbps CO-OFDM system in presence of timing and frequency offsets | |
WO2023248285A1 (ja) | マルチキャリア信号波形等化回路及びマルチキャリア信号波形等化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |