CN104980379B - 非线性失真的估计装置、方法以及接收机 - Google Patents
非线性失真的估计装置、方法以及接收机 Download PDFInfo
- Publication number
- CN104980379B CN104980379B CN201410145340.3A CN201410145340A CN104980379B CN 104980379 B CN104980379 B CN 104980379B CN 201410145340 A CN201410145340 A CN 201410145340A CN 104980379 B CN104980379 B CN 104980379B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- disturbance
- nonlinear
- nonlinear parameteric
- parameteric
- 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
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/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
- H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
- H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
- H04B10/0795—Performance monitoring; Measurement of transmission parameters
-
- 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/60—Receivers
- H04B10/61—Coherent receivers
- H04B10/616—Details of the electronic signal processing in coherent optical receivers
- H04B10/6163—Compensation of non-linear effects in the fiber optic link, e.g. self-phase modulation [SPM], cross-phase modulation [XPM], four wave mixing [FWM]
Abstract
本发明实施例提供一种非线性失真的估计装置、方法以及接收机。所述估计方法包括:对带限模拟信号进行采样以获得采样序列;基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;以及利用所述非线性微扰项计算非线性失真波形。通过本发明实施例,不仅可以兼容任意调制格式,而且具有精度高,普适性好等优点。
Description
技术领域
本发明涉及长距离光纤通信系统,尤其涉及一种非线性失真的估计装置、方法以及接收机。
背景技术
基于慢变包络近似和恒定偏振态假设,光纤内脉冲演化的传输方程可由非线性薛定谔方程来描述(例如随机偏振下用Manakov方程描述)。此传输方程用于描述光脉冲信号在色散和克尔(Kerr)效应联合作用下的波形演化。但由于非线性薛定谔方程在考虑非线性和色散效应共同作用下没有解析解,故针对光纤非线性损伤的定量研究以及相关的理论模型都是针对非线性薛定谔方程的近似解法发展和建立的。
由于近似解析方法有望显著减小非线性分析的计算复杂度,因而受到了学术界的广泛关注并且在近些年得到了迅速的发展。Volterra级数展开方法作为求解非线性薛定谔方程的一种普适方法,使传统通信系统的分析框架可以被借用到光纤通信系统,并且对不同的脉冲形状和链路类型具有较好的通用性。
Paolo Serena基于Volterra展开方法发展得到了常规微扰法(RP)并赋予各阶微扰较明确的物理意义,从而使微扰求解薛定谔的方法得到了迅速的发展,衍生出了多种理论框架用于在时域或频域定量非线性失真。普遍结果表明,对于典型的长距离光纤传输系统,非线性作用主要由三阶以下的Volterra级数(一阶微扰)充分描述,故目前流行的非线性分析均接受低阶Volterra级数展开的分析框架,即准线性近似。在准线性近似下,用于求解非线性传输方程的一阶微扰框架可概括为求解经过色散(线性)作用的脉冲在传播路径上各点所受非线性失真的矢量和,解析表达为以发送脉冲的时域三项乘积为被积函数的三重积分。
理论分析表明,一阶微扰的解析表达在某些条件下可以进行简化,从而减小微扰方法的计算复杂度。目前,最典型且成功的理论近似为无损大色散链路的解析解,此方法假设光纤传输链路无损耗且积累色散足够大,同时保证发送数字序列的承载脉冲为高斯形状。在以上近似下,一阶微扰的三重积分严格可积,可表达为特殊函数的闭解形式。此方法虽可大幅降低计算复杂度,但由于存在关键的高斯脉冲近似,在应用比较普遍的非高斯脉冲传输系统中的计算精度受限,从而限制了此方法的应用范围。
随着高速数字信号处理(DSP,Digital Signal Process)以及窄带光滤波技术的成熟,具有高频谱利用率的光正交频分复用技术(OOFDM,Optical Orthogonal FrequencyDivision Multiplexing)以及奈奎斯特波分复用(Nyquist-WDM)技术逐渐受到重视。在高频谱利用率传输系统中,由于信号谱密度进一步加大,非线性损伤加剧,对系统的功率预算和传输距离造成可观的负面影响。在这种背景下,为了实现更准确的相干传输系统性能估计和寻找更优的系统设计准则,精确的非线性理论模型研究具有显著意义。
应该注意,上面对技术背景的介绍只是为了方便对本发明的技术方案进行清楚、完整的说明,并方便本领域技术人员的理解而阐述的。不能仅仅因为这些方案在本发明的背景技术部分进行了阐述而认为上述技术方案为本领域技术人员所公知。
下面列出了对于理解本发明和常规技术有益的文献,通过引用将它们并入本文中,如同在本文中完全阐明了一样。
[非专利文献1]:K.V.Peddanarappagari et.al.,IEEE JLT Vol.15,pp.2232-2241,1997
[非专利文献2]:IEEE JLT Vol.16,pp.2046-1055,1998
[非专利文献3]:A.Vannucci et.al.,IEEE JLT Vol.20,No.7,pp.1102-1111,2002
[非专利文献4]:Z.Tao et al.,IEEE JLT Vol.29,pp.2570-2576,2011
[非专利文献5]:A.Carena et.al.,IEEE JLT Vol.30,No.10,pp.1524-1539,2012
[非专利文献6]:A.Mecozzi et.al.,IEEE PTL Vol.12,No.4,pp.392-394,2000
[非专利文献7]:S.Kumar et.al.,Optics Express,Vol.20,No.25,pp.27740-27754,2012
[非专利文献8]:Y.Zhao et al.,ECOC2013,P.4.15.
[非专利文献9]:Y.Fan et al.,ECOC2012,We.2.C.3
发明内容
目前,对非线性理论模型的需求不再局限于高斯脉冲形状,而是向Nyquist脉冲形状甚至任意波形方向展开。基于奈奎斯特(Nyquist)采样定理,任意带限信号均可以由Nyquist采样脉冲表示,所以适用于Nyquist脉冲形状的非线性估计方法是展开具有普适性的非线性模型研究的重要基础。另一方面,从非线性补偿和提升系统性能角度,更高精度的非线性估计模型同时是非线性补偿方法的基础,从而有利于减小非线性损伤,提高系统性能或减小非线性补偿系统复杂度。
本发明实施例提供一种非线性失真的估计装置、方法以及接收机。不仅可以兼容任意调制格式,而且具有精度高,普适性好等优点。
根据本发明实施例的一个方面,提供一种非线性失真的估计装置,所述估计装置包括:
信号采样单元,对带限模拟信号进行采样以获得采样序列;
系数计算单元,基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;
微扰项计算单元,利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;
波形估计单元,利用所述非线性微扰项计算非线性失真波形。
根据本发明实施例的另一个方面,提供一种非线性失真的估计方法,所述估计方法包括:
对带限模拟信号进行采样以获得采样序列;
基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;
利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;以及
利用所述非线性微扰项计算非线性失真波形。
根据本发明实施例的另一个方面,提供一种接收机,其中,所述接收机包括上所述的非线性失真的估计装置。
本发明的有益效果在于:对带限模拟信号进行采样,并基于奈奎斯特脉冲计算非线性微扰系数来估计非线性失真波形;不仅可以兼容任意调制格式,而且具有精度高,普适性好等优点。
参照后文的说明和附图,详细公开了本发明的特定实施方式,指明了本发明的原理可以被采用的方式。应该理解,本发明的实施方式在范围上并不因而受到限制。在所附权利要求的精神和条款的范围内,本发明的实施方式包括许多改变、修改和等同。
针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用,与其它实施方式中的特征相组合,或替代其它实施方式中的特征。
应该强调,术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在,但并不排除一个或更多个其它特征、整件、步骤或组件的存在或附加。
附图说明
所包括的附图用来提供对本发明实施例的进一步的理解,其构成了说明书的一部分,用于例示本发明的实施方式,并与文字描述一起来阐释本发明的原理。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1是长距离光纤传输系统的一结构示意图;
图2是本发明实施例1的估计方法的一流程示意图;
图3是本发明实施例1的任意波形非线性估计模型的一结构示意图;
图4是本发明实施例1的计算Nyquist脉冲形状的非线性微扰系数的一示意图;
图5是本发明实施例1的Nyquist影子脉冲和高斯影子脉冲对比的一示意图;
图6是本发明实施例2的估计装置的一构成示意图;
图7是本发明实施例2的系数计算单元的一构成示意图;
图8是本发明实施例3的接收机的一构成示意图。
具体实施方式
参照附图,通过下面的说明书,本发明的前述以及其它特征将变得明显。在说明书和附图中,具体公开了本发明的特定实施方式,其表明了其中可以采用本发明的原则的部分实施方式,应了解的是,本发明不限于所描述的实施方式,相反,本发明包括落入所附权利要求的范围内的全部修改、变型以及等同物。
本发明实施例适用于长距离光纤通信系统,对数据信号在传输过程中受到的非线性失真进行定量估计。图1是长距离光纤传输系统的一结构示意图,如图1所示,发射机发射的信号经过传输链路中不同的器件(例如光纤、光放大器、色散补偿光纤等)到达接收机。在图1所示的系统中,需要在接收机端进行非线性失真估计。
本发明实施例在一阶微扰框架下,提出一种适用于任意带限波形的非线性失真估计方法以及装置。该方法具有精度高,普适性好等优点,可应用于无色散补偿链路(NDM)与色散管理链路(DM),并且兼容任意调制格式;同时也适用于单偏振或偏振复用系统中。该估计方法以Nyquist脉冲的非线性微扰系数(也可称为加权系数)计算为基础,给出采样序列的非线性微扰计算方法,形成一套基于采样、系数计算、微扰及非线性失真计算的任意波形非线性失真估计模型。
实施例1
本发明实施例提供一种非线性失真的估计方法。图2是本发明实施例的估计方法的一流程示意图,如图2所示,所述估计方法包括:
步骤201,对带限模拟信号进行采样以获得采样序列;
步骤202,基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;
步骤203,利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;以及
步骤204,利用所述非线性微扰项计算非线性失真波形。
在本实施例中,任意带限信号的非线性估计模型以信号波形的采样值以及链路参数为输入,以非线性失真波形为输出。值得注意的是,图2中各个步骤的执行顺序不限于此,例如可以根据实际需要调整步骤201和步骤202的顺序,或者同时执行步骤201和202等。
图3是本发明实施例的任意波形非线性估计模型的一结构示意图,如图3所示,本发明实施例主要由四部分功能实现,分别完成:模拟信号的采样,非线性微扰系数计算,非线性微扰项计算及失真波形的计算。以下结合图2和图3对各部分进行详细说明。
在步骤201中,可以对任意带限模拟信号进行采样。Nyquist采样定理指出,带限模拟信号u(t)(例如频谱占据-B~+B)可以用抽样间隔不大于1/2B的抽样值唯一地表示。在满足抽样定理的前提下,u(t)可表示成正交抽样函数(sinc函数)的无穷级数,如式(1)所示,其中T为采样周期。该无穷级数的基为彼此正交的sinc函数(Nyquist脉冲),系数为u(t)的抽样值u(kT)。式(1)表明带限信号可由其抽样值序列经过矩形滤波器无失真恢复。
在不满足抽样定理前提条件下,即采样频率fs之2B情况下,u(t)无法由采样序列无失真地表示。此时用采样序列来表示模拟信号会引入不可避免的频域混叠失真,失真的大小与信号频谱和采样率有关。
综上所述,无论是否引入混叠,带限模拟信号都可用Nyquist脉冲序列来表示。从非线性模型的适用性角度考虑,若理论模型能够有效描述Nyquist脉冲序列的非线性演化特性,则可基于采样定理应用于任意波形的非线性失真分析。本发明实施例不对脉冲序列的采样率进行限制,理论分析在有无混叠失真的情况下均适用。
在步骤202中,可以基于Nyquist脉冲计算非线性失真估计中的非线性微扰系数,其中整个传输链路的光纤可以被划分为多个光纤跨段,对于每一光纤跨段,可以进行如下处理:基于链路参数对每个光纤跨段的非线性微扰系数进行积分处理;对积分处理后的所述非线性微扰系数进行大色散近似;利用所述奈奎斯特脉冲对大色散近似后的所述非线性微扰系数进行简化,以获得每个光纤跨段的一重积分形式的非线性微扰系数。
在本实施例中,可以对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的一重积分形式的非线性微扰系数。
此外,还可以利用有理函数对信道内非线性失真估计中的链路损耗或增益函数进行近似处理,通过近似后的所述链路损耗或增益函数对所述一重积分形式的非线性微扰系数进行计算,以获得每个光纤跨段的解析形式的非线性微扰系数,并且,对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的解析形式的非线性微扰系数。
此外,还可以忽略所述链路损耗或增益函数后对所述一重积分形式的非线性微扰系数进行计算,以获得每个光纤跨段的解析形式的非线性微扰系数;并且,对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的解析形式的非线性微扰系数。
图4是本发明实施例的计算Nyquist脉冲形状的非线性微扰系数的一示意图,以下基于图4,对本发明实施例的非线性微扰系数的计算进行详细说明。
在已知发送序列脉冲形状及链路参数时,不考虑任何近似的情况下,对于第i个光纤跨段,一阶非线性微扰加权系数可表达为
其中,α,γ为光纤的衰减系数与非线性系数;m,n表示非线性作用的三个脉冲位置分别为t=mT,nT,(m+n)T;g(t,z)与分别为符号脉冲传输至zi处的时域和频域表达;zi表示信号在第i个跨段的传输距离;Ls为光纤跨段长度;β2为色散系数,ζ为残余色散率,C(zi)=β2(i-1)ζLs+β2zi为第i个跨段的积累色散;F(·)为傅里叶变换;F-1(·)为反傅里叶变换;“*”表示卷积操作。
从(2)式可以看出,无近似情况下,非线性微扰系数可以统一表达为四重积分形式(其中,F-1(·)为一重积分,频域二维卷积为二重积分),从而通过数值方法进行计算。
由于微扰系数的计算是非线性失真估计的核心,目前研究趋势统一向如何简化式(2)所示的四重积分展开,以便能以低复杂度硬件实现高精度的非线性估计。以下介绍使式(2)在特定情况下近似成一重积分或解析解的过程。
首先引入大色散近似,在只有色散作用下符号脉冲的演化可用远场近似来描述
远场近似的物理解释为在色散足够大的情况下,脉冲的形状演化为发射脉冲的Fourier变换,以及长距离传输把脉冲频域形状映射到时域,其原理与Fraunhofer远场衍射积分类似,故称为远场近似。
将(3)式带入(2)式,四重积分可简化为二重积分表达式
式(4)表示二维卷积在大色散近似下可用乘法近似。
当考虑采样点的Nyquist成型脉冲(滚降系数为0)时,式(4)可以被进一步简化。考虑Nyquist脉冲的频谱
其中fs=f/T为采样率。将(5)式带入(4)式中,积分因子成为
其中重叠频谱表示为
lm,n与rm,n决定重叠频谱的上下限
由于需要保证lm,n≤rm,n,对于确定的m和n,积分范围别进一步约束为
将(6)式带入(4)式,并对矩形重叠脉冲做反Fourier变换,第i个跨段的时域非线性微扰系数可表达为
(10)式表示Nyquist成型脉冲的非线性微扰系数可以表示对传输距离的一维积分形式,此一维积分可以通过数值方法计算,相对于(2)式所示的4维积分,计算复杂度可显著降低。在数值计算形式下,整个传输链路的微扰系数可以对各光纤跨段直接求和,得到一维积分形式的微扰系数表达。
对于(11)式,一维积分不存在解析形式的原因为积分因子中存在指数衰减因子。
为进一步得到微扰系数的解析解,可以考虑将指数衰减忽略,即类似于高斯脉冲形状下Mecozzi模型对衰减因子的处理方法,令在此近似下可以得到(12)式的解析表达
其中Leff=1/α定义为有效传输距离,min,max分别表示min(m,n,m+n)和max(m,n,m+n)。对(12)式的跨段指数求和即可得到总微扰系数的解析表达。
对衰减因子的另外一种处理方式是:利用特定的有理函数对指数衰减函数进行替换,从而达到令(11)式可积的目的。这种处理方法相对忽略衰减的方法显示出更好的估计精度,是处理指数衰减的有效手段之一。
以上给出了Nyquist成型脉冲的微扰系数计算方法,该方法适用于描述任意波形以任意调制格式在光纤链路中传输的非线性失真演化过程,且不对传输链路施加限制,是一种通用的非线性失真估计方法。
在步骤203中,可以利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项。其中,在计算所述非线性微扰项时加入非线性符号间串扰的影响。
在本实施例中,Cm,n(t)的物理意义为坐落在t=mT,nT,(m+n)T三个位置的Nyquist脉冲在t=0处产生的影子脉冲的时域波形,整个脉冲序列的非线性失真波形为各时间的影子脉冲的叠加。故影子脉冲的性质一定程度上决定非线性失真的特性,对影子脉冲的进一步讨论和非线性失真项的计算方法有必要进一步讨论。
非线性微扰项的传统表示为脉冲序列三项乘积加权和的形式,在偏振复用系统中表达为
其中,ΔuH/V(kT)表示kT时刻H或V偏振态上的非线性失真,表示发送序列在mT时刻H或V偏振态上的数字信息。该表达式只使用非线性微扰系数在t=0时刻的信息,隐含的假设是非线性影子脉冲的宽度足够窄,使影子脉冲对其它样点的影响可以忽略不计。但上述假设只适用于高斯脉冲序列,对于任意波形,理想采样脉冲为Nyquist脉冲,其非线性影子脉冲的性质与高斯影子脉冲性质不同,表现为非线性微扰系数的表达式不同。
图5是本发明实施例的Nyquist影子脉冲和高斯影子脉冲对比的一示意图。不失一般性地,图5所示为m=20T,n=40T的两对高斯脉冲或Nyquist脉冲在t=0出产生的影子脉冲形状。仿真条件为脉冲周期T=31.25ps,脉冲平均功率为-3dBm,传输链路为500×17ps/nm预色散的100km单模光纤。
如图5所示,Nyquist影子脉冲宽度明显大于高斯脉冲宽度,在t=±T时刻,高斯影子脉冲的影响可以忽略,但Nyquist影子脉冲影响显著。所以t=(m+k)T,t=(n+k)T以及t=(m+n+k)T的Nyquist脉冲不仅在t=kT处产生非线性微扰,还会在t=(k±1)T,(k±2)T,…处产生微扰。
综上,可以称坐落在tt(mmk)T,(nnk)T,(m+n+k)T的Nyquist脉冲在t≠kT时刻的非线性影子脉冲贡献为非线性符号串扰(NL-ISI)。由于考虑由Nyquist采样脉冲构成的任意波形非线性失真,所以本发明实施例提出在计算非线性微扰项时要考虑非线性符号串扰项。此情况下,非线性微扰项表达式变为
即在双偏振态下,
其中,l为整数,T为符号周期;和分别表示水平偏振态和垂直偏振态上第m+k+l时刻的脉冲的符号信息;和分别表示水平偏振态和垂直偏振态上第n+k+l时刻的脉冲的符号信息;和分别表示水平偏振态和垂直偏振态上第m+n+k+l时刻的脉冲的符号信息的共轭;
或者,在单偏振态下,所述非线性微扰项的表达式为:
其中,l为整数,T为符号周期;Am+k+l表示单偏振态上第m+k+l时刻的脉冲的符号信息;An+k+l表示单偏振态上第n+k+l时刻的脉冲的符号信息;(Am+n+k+l)*表示单偏振态上第m+n+k+l时刻的脉冲的符号信息的共轭。
在本实施例中,考虑非线性符号串扰的微扰系数解析表达式可由(12)式对光纤跨段直接求和得到
(14)式利用了非线性微扰系数在除0外其它采样点的信息来计算非线性微扰项,应用此式计算非线性微扰适用性更强,结果准确度更高。在利用(14)式计算非线性微扰时,由于l的取值范围有限所以需要对l进行截断操作,由图5所示,非线性符号间串扰在影子脉冲临近的±T,±2T处影响较为显著,其它位置影响较小,故l的截断范围一般取-3≤l≤3。
以下介绍步骤204中非线性失真波形的计算,实际上,一阶非线性微扰理论把(14)式所示的非线性微扰项完全视为加性微扰,从而在线性解的基础上直接叠加非线性微扰项即可视为失真波形的一种计算方法,即
其中,为非线性薛定谔方程的线性解。
(15)式是传统的非线性失真波形计算表达式。虽然(15)式是一阶微扰法的理论解,但把非线性微扰项统一视为加性项存在一定的不合理性,这种不合理性实际上源于一阶微扰法的准线性近似在非线性作用较强时的理论偏差。
进一步研究表明,适当引入高阶非线性项,可以提高估计方法的精度、增加估计方法的功率适用范围,这方面比较典型的结果是之前提出的非线性加乘(AM)模型。加乘模型从非线性脉冲相互作用产生的结果出发,讨论了非线性微扰项在涉及三脉冲相互作用与涉及两脉冲相互作用的区别,指出涉及两脉冲相互作用的微扰项实际表现为非线性相位噪声,从而得到了更具应用价值的非线性失真波形计算方法,以uH(kT)为例,表达为
其中
ΔIFWM(kT)为信道内四波混频项,φIXPM(kT)为信道内交叉相位调制项,分别对应三脉冲相互作用和两脉冲相互作用引入的非线性失真,其中三脉冲相互作用失真仍视作加性项,而两脉冲失真则视作相位失真,表现为与线性项的乘积形式。(16)式在针对高斯脉冲的非线性失真波形计算中表现出良好的计算精确度。
在Nyquist脉冲形状情况下,由于需要考虑非线性符号间串扰,(16a)式与(16b)式所涉及的非线性三项乘积项数增加,在接受加乘模型描述的物理意义的条件下,需要对加乘模型在Nyquist脉冲下的表达式进行修正,修正的理论依据为没有本符号参与的三项乘积项归纳为ΔIFWM(kT)项,而本符号参与的三项乘积项归纳为φIXPM(kT)项。
基于上述准则,(16)式的ΔIFWM(kT)与φIXPM(kT)表达修正为
上两式中的最后一项为考虑非线性符号间串扰所增加的三项乘积项。所以,在Nyquist脉冲形状下,基于加乘模型的非线性失真波形由式(16)和(17a)式(17b)式共同描述。当不考虑非线性符号间串扰时,(17a)式和(17b)式将退化为(16a)式与(16b)式。
由于利用(11)式或(14a)式计算微扰系数比较复杂,在一般系统中(色散补偿或色散管理链路),用高斯脉冲假设下的C0,0(0)项的计算方法(一重积分表达)带入(17b)式并不会引入很大的误差,所以在实际微扰计算中C0,0(0)项常利用高斯脉冲假设进行简化。
由上述实施例可知,对带限模拟信号进行采样,并基于奈奎斯特脉冲计算非线性微扰系数来估计非线性失真波形;不仅可以兼容任意调制格式,而且具有精度高,普适性好等优点。
实施例2
本发明实施例提供一种非线性失真的估计装置,对应于实施例1中的非线性失真的估计方法,相同的内容不再赘述。
图6是本发明实施例的估计装置的一构成示意图,如图6所示,估计装置600包括:信号采样单元601、系数计算单元602、微扰项计算单元603和波形估计单元604;
其中,信号采样单元601对带限模拟信号进行采样以获得采样序列;系数计算单元602基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;微扰项计算单元603利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;波形估计单元604利用所述非线性微扰项计算非线性失真波形。
图7是本发明实施例的系数计算单元的一构成示意图,如图7所示,系数计算单元602可以包括:积分处理单元701、第一近似单元702、第二近似单元703和系数求和单元704;
其中,积分处理单元701基于链路参数对每个光纤跨段的非线性微扰系数进行积分处理;第一近似单元702对积分处理后的所述非线性微扰系数进行大色散近似;第二近似单元703利用所述奈奎斯特脉冲对大色散近似后的所述非线性微扰系数进行简化,以获得每个光纤跨段的一重积分形式的非线性微扰系数;系数求和单元704对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的一重积分形式的非线性微扰系数。
如图7所示,所述系数计算单元602还可以包括:第三近似单元705和第四近似单元706;其中第三近似单元705利用有理函数对信道内非线性失真估计中的链路损耗或增益函数进行近似处理;第四近似单元706通过近似后的所述链路损耗或增益函数对所述一重积分形式的非线性微扰系数进行计算,或者忽略所述链路损耗或增益函数后对所述一重积分形式的非线性微扰系数进行计算,以获得每个光纤跨段的解析形式的非线性微扰系数;所述系数求和单元704还用于对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的解析形式的非线性微扰系数。
在本实施例中,所述微扰项计算单元603在计算所述非线性微扰项时加入非线性符号间串扰的影响。其中,考虑非线性符号间串扰的非线性微扰系数的解析表达式可以为:
其中,l为整数,T为符号周期;Leff=1/α为有效传输距离,α为光纤的衰减系数;m,n表示非线性作用的三个脉冲位置分别为t=mT,nT,(m+n)T,min和max分别表示min(m,n,m+n)和maxxm,n,m+nn;i表示信号在第i个跨段的传输距离;Ls为光纤跨段长度;β2为色散系数,ζ为残余色散率,为第i个跨段的积累色散;Ei表示指数积分函数。
由上述实施例可知,对带限模拟信号进行采样,并基于奈奎斯特脉冲计算非线性微扰系数来估计非线性失真波形;不仅可以兼容任意调制格式,而且具有精度高,普适性好等优点。
实施例3
本发明实施例提供一种接收机,该电子设备包括如实施例2所述的估计装置600。
图8是本发明实施例的接收机的一构成示意图。如图8所示,接收机800可以包括:中央处理器(CPU)200和存储器210;存储器210耦合到中央处理器200。其中该存储器210可存储各种数据;此外还存储信息处理的程序,并且在中央处理器200的控制下执行该程序。
在一个实施方式中,估计装置600的功能可以被集成到中央处理器200中。其中,中央处理器200可以被配置为实现如实施例1所述的估计方法。
在另一个实施方式中,估计装置600可以与中央处理器200分开配置,例如可以将估计装置600配置为与中央处理器200连接的芯片,通过中央处理器200的控制来实现估计装置600的功能。
此外,如图8所示,接收机800还可以包括:收发机220和天线230等;其中,上述部件的功能与现有技术类似,此处不再赘述。值得注意的是,接收机800也并不是必须要包括图8中所示的所有部件;此外,接收机800还可以包括图8中没有示出的部件,可以参考现有技术。
本发明实施例还提供一种计算机可读程序,其中当在电子设备中执行所述程序时,所述程序使得计算机在所述接收机中执行实施例1所述的非线性失真的估计方法。
本发明实施例还提供一种存储有计算机可读程序的存储介质,其中所述计算机可读程序使得计算机在接收机中执行实施例1所述的非线性失真的估计方法。
本发明以上的装置和方法可以由硬件实现,也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序,当该程序被逻辑部件所执行时,能够使该逻辑部件实现上文所述的装置或构成部件,或使该逻辑部件实现上文所述的各种方法或步骤。本发明还涉及用于存储以上程序的存储介质,如硬盘、磁盘、光盘、DVD、flash存储器等。
以上结合具体的实施方式对本发明进行了描述,但本领域技术人员应该清楚,这些描述都是示例性的,并不是对本发明保护范围的限制。本领域技术人员可以根据本发明的精神和原理对本发明做出各种变型和修改,这些变型和修改也在本发明的范围内。
Claims (10)
1.一种非线性失真的估计装置,其特征在于,所述估计装置包括:
信号采样单元,对带限模拟信号进行采样以获得采样序列;
系数计算单元,基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;
微扰项计算单元,利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;
波形估计单元,利用所述非线性微扰项计算非线性失真波形。
2.根据权利要求1所述的估计装置,其中,所述系数计算单元包括:
积分处理单元,基于链路参数对每个光纤跨段的非线性微扰系数进行积分处理;
第一近似单元,对积分处理后的所述非线性微扰系数进行大色散近似;
第二近似单元,利用所述奈奎斯特脉冲对大色散近似后的所述非线性微扰系数进行简化,以获得每个光纤跨段的一重积分形式的非线性微扰系数;
系数求和单元,对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的一重积分形式的非线性微扰系数。
3.根据权利要求2所述的估计装置,其中,所述系数计算单元还包括:
第三近似单元,利用有理函数对信道内非线性失真估计中的链路损耗或增益函数进行近似处理;
第四近似单元,通过近似后的所述链路损耗或增益函数对所述一重积分形式的非线性微扰系数进行计算,或者忽略所述链路损耗或增益函数后对所述一重积分形式的非线性微扰系数进行计算,以获得每个光纤跨段的解析形式的非线性微扰系数;
所述系数求和单元还用于对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的解析形式的非线性微扰系数。
4.根据权利要求1所述的估计装置,其中,所述微扰项计算单元在计算所述非线性微扰项时加入非线性符号间串扰的影响。
5.根据权利要求4所述的估计装置,其中,考虑非线性符号间串扰的非线性微扰系数的解析表达式为:
其中,l为整数,T为符号周期;Leff=1/α为有效传输距离,α为光纤的衰减系数;m,n表示非线性作用的三个脉冲位置分别为t=mT,nT,(m+n)T时的参数;min和max分别表示min(m,n,m+n)和max(m,n,m+n);zi表示信号在第i个跨段的传输距离;Ls为光纤跨段长度;β2为色散系数,ζ为残余色散率,C(zi)=β2(i-1)ζLs+β2zi为第i个跨段的积累色散;Ei表示指数积分函数;fs为采样频率。
6.一种非线性失真的估计方法,其特征在于,所述估计方法包括:
对带限模拟信号进行采样以获得采样序列;
基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数;
利用所述非线性微扰系数以及所述采样序列计算叠加在信号上的非线性微扰项;以及
利用所述非线性微扰项计算非线性失真波形。
7.根据权利要求6所述的估计方法,其中,基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数包括:
基于链路参数对每个光纤跨段的非线性微扰系数进行积分处理;
对积分处理后的所述非线性微扰系数进行大色散近似;
利用所述奈奎斯特脉冲对大色散近似后的所述非线性微扰系数进行简化,以获得每个光纤跨段的一重积分形式的非线性微扰系数;以及
对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的一重积分形式的非线性微扰系数。
8.根据权利要求7所述的估计方法,其中,基于奈奎斯特脉冲计算非线性失真估计中的非线性微扰系数还包括:
利用有理函数对信道内非线性失真估计中的链路损耗或增益函数进行近似处理;
通过近似后的所述链路损耗或增益函数对所述一重积分形式的非线性微扰系数进行计算,或者忽略所述链路损耗或增益函数后对所述一重积分形式的非线性微扰系数进行计算,以获得每个光纤跨段的解析形式的非线性微扰系数;
并且,对分别获得的不同光纤跨段的所述非线性微扰系数进行求和,以得到整个传输链路的解析形式的非线性微扰系数。
9.根据权利要求6所述的估计方法,其中,所述方法还包括:在计算所述非线性微扰项时加入非线性符号间串扰的影响。
10.一种接收机,其中,所述接收机包括如权利要求1至5任一项所述的非线性失真的估计装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410145340.3A CN104980379B (zh) | 2014-04-11 | 2014-04-11 | 非线性失真的估计装置、方法以及接收机 |
JP2015079289A JP2015204622A (ja) | 2014-04-11 | 2015-04-08 | 非線形ディストーションの推定装置、方法及び受信機 |
US14/683,523 US9654211B2 (en) | 2014-04-11 | 2015-04-10 | Estimation apparatus and method for nonlinear distortion and receiver |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410145340.3A CN104980379B (zh) | 2014-04-11 | 2014-04-11 | 非线性失真的估计装置、方法以及接收机 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104980379A CN104980379A (zh) | 2015-10-14 |
CN104980379B true CN104980379B (zh) | 2018-06-08 |
Family
ID=54265949
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410145340.3A Active CN104980379B (zh) | 2014-04-11 | 2014-04-11 | 非线性失真的估计装置、方法以及接收机 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9654211B2 (zh) |
JP (1) | JP2015204622A (zh) |
CN (1) | CN104980379B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9564966B1 (en) * | 2015-09-30 | 2017-02-07 | Osram Sylvania Inc. | Reconstructing light-based communication signals using an alias frequency |
JP6759742B2 (ja) * | 2016-06-16 | 2020-09-23 | 富士通株式会社 | 受信装置及び設定方法 |
US10110263B2 (en) * | 2017-02-01 | 2018-10-23 | Roshmere, Inc. | Crosstalk-correction in digitizers using coupling coefficients |
WO2018144579A1 (en) * | 2017-02-01 | 2018-08-09 | Roshmere, Inc. | Crosstalk correction using pre-compensation |
JP7069578B2 (ja) * | 2017-07-05 | 2022-05-18 | 富士通株式会社 | 光伝送装置及び方法 |
US10236982B1 (en) * | 2017-12-21 | 2019-03-19 | Ciena Corporation | Fiber parameter identification |
CN109714100B (zh) * | 2017-12-26 | 2021-08-06 | 电子科技大学 | 一种多波长信道的光纤非线性串扰计算方法 |
US10985951B2 (en) | 2019-03-15 | 2021-04-20 | The Research Foundation for the State University | Integrating Volterra series model and deep neural networks to equalize nonlinear power amplifiers |
JP7322196B2 (ja) * | 2019-06-21 | 2023-08-07 | フラウンホッファー-ゲゼルシャフト ツァ フェルダールング デァ アンゲヴァンテン フォアシュンク エー.ファオ | 周期的周波数領域における摂動ベースファイバ非線形性補償を用いた光ファイバチャンネルを介したデータを送受信するための装置および方法 |
CN111800194B (zh) * | 2020-06-22 | 2021-06-18 | 北京理工大学 | 针对少模多芯oam光纤传输概率分布的非线性补偿方法 |
US11190272B1 (en) | 2020-12-18 | 2021-11-30 | Huawei Technologies Co., Ltd. | Method and apparatus for evaluating nonlinear impairment of an optical fiber link |
CN114422035B (zh) * | 2022-01-21 | 2023-08-25 | 上海交通大学 | 相干光通信系统中光纤非线性损伤补偿方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102420660A (zh) * | 2010-09-28 | 2012-04-18 | 富士通株式会社 | 非线性补偿装置和发射机 |
CN103259597A (zh) * | 2012-02-20 | 2013-08-21 | 富士通株式会社 | 非线性补偿装置、方法和发射机 |
CN103312414A (zh) * | 2012-03-16 | 2013-09-18 | 富士通株式会社 | 一种逆信道装置和包含该装置的发射机、接收机及系统 |
CN103457606A (zh) * | 2012-05-30 | 2013-12-18 | 英飞凌科技奥地利有限公司 | 用于对输入信号进行采样的方法和装置 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060245766A1 (en) * | 2005-04-29 | 2006-11-02 | Taylor Michael G | Phase estimation for coherent optical detection |
US7570182B2 (en) * | 2006-09-15 | 2009-08-04 | Texas Instruments Incorporated | Adaptive spectral noise shaping to improve time to digital converter quantization resolution using dithering |
US8204389B2 (en) * | 2008-02-27 | 2012-06-19 | University Of Central Florida Research Foundation, Inc. | Electronic post-compensation of optical transmission impairments using digital backward propagation |
EP2146448B1 (en) * | 2008-07-16 | 2010-11-17 | Alcatel Lucent | Adaptive non-linearity compensation in coherent receiver |
WO2010094339A1 (en) * | 2009-02-20 | 2010-08-26 | Telefonaktiebolaget Lm Ericsson (Publ) | Equaliser for an optical transmission system |
US8306418B2 (en) * | 2009-03-10 | 2012-11-06 | Tyco Electronics Subsea Communications Llc | Data pattern dependent distortion compensation in a coherent optical signal receiver |
KR101301147B1 (ko) * | 2009-12-04 | 2013-09-10 | 한국전자통신연구원 | 광 오에프디엠 수신기 및 그것의 광신호 수신 방법 |
WO2011106826A2 (en) * | 2010-03-05 | 2011-09-09 | Ofidium Pty Ltd | Method and system for non-linearity compensation in optical transmission systems |
US8494368B2 (en) * | 2010-04-16 | 2013-07-23 | Alcatel Lucent | Electronic nonlinearity compensation for optical transmission systems |
US8693890B2 (en) * | 2010-08-20 | 2014-04-08 | Nec Laboratories America, Inc. | Look-up table and digital transmitter based architecture for fiber nonlinearity compensation |
AU2011349107B2 (en) * | 2010-12-22 | 2016-08-25 | Ofidium Pty Ltd | Electronic compensation of cross-phase modulation |
JP5327397B2 (ja) * | 2011-02-17 | 2013-10-30 | 日本電気株式会社 | 信号処理回路、信号処理方法、光受信機及び光通信システム |
CN102655432B (zh) * | 2011-03-04 | 2015-08-26 | 富士通株式会社 | 非线性损伤补偿方法和装置 |
CN103023828B (zh) * | 2011-09-26 | 2016-02-17 | 富士通株式会社 | 非线性估计装置、方法和接收机 |
EP2660984B1 (en) * | 2011-12-01 | 2016-03-23 | Huawei Technologies Co., Ltd. | Optical signal compensation device |
US8548097B1 (en) * | 2012-06-20 | 2013-10-01 | MagnaCom Ltd. | Coarse phase estimation for highly-spectrally-efficient communications |
US9184854B2 (en) * | 2012-10-09 | 2015-11-10 | Nec Laboratories America, Inc. | Inter-band cross-phase modulation compensation for the mitigation of intra-channel nonlinear impairments in optical fiber transmission |
JP6287866B2 (ja) * | 2013-02-07 | 2018-03-07 | 日本電気株式会社 | 信号処理装置及び信号処理方法 |
US9853765B2 (en) * | 2013-04-09 | 2017-12-26 | Nec Corporation | Signal processing device and signal processing method for optical polarization multiplexed signal |
JP6135415B2 (ja) * | 2013-09-11 | 2017-05-31 | 富士通株式会社 | 非線形歪み補償装置及び方法並びに光受信器 |
-
2014
- 2014-04-11 CN CN201410145340.3A patent/CN104980379B/zh active Active
-
2015
- 2015-04-08 JP JP2015079289A patent/JP2015204622A/ja active Pending
- 2015-04-10 US US14/683,523 patent/US9654211B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102420660A (zh) * | 2010-09-28 | 2012-04-18 | 富士通株式会社 | 非线性补偿装置和发射机 |
CN103259597A (zh) * | 2012-02-20 | 2013-08-21 | 富士通株式会社 | 非线性补偿装置、方法和发射机 |
CN103312414A (zh) * | 2012-03-16 | 2013-09-18 | 富士通株式会社 | 一种逆信道装置和包含该装置的发射机、接收机及系统 |
CN103457606A (zh) * | 2012-05-30 | 2013-12-18 | 英飞凌科技奥地利有限公司 | 用于对输入信号进行采样的方法和装置 |
Non-Patent Citations (2)
Title |
---|
Aggressive Quantization on Perturbation Coefficients for Nonlinear Pre-Distortion;Qunbi Zhuge 等;《OSA》;20140313;全文 * |
Experimental Demonstration of 448-Gbps+DMT Transmission over 30-km SMF;Toshiki Tanaka 等;《OSA》;20140313;全文 * |
Also Published As
Publication number | Publication date |
---|---|
US20150295643A1 (en) | 2015-10-15 |
US9654211B2 (en) | 2017-05-16 |
CN104980379A (zh) | 2015-10-14 |
JP2015204622A (ja) | 2015-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104980379B (zh) | 非线性失真的估计装置、方法以及接收机 | |
CN104104445B (zh) | 非线性加权系数的计算装置以及方法 | |
US7643752B2 (en) | Testing of transmitters for communication links by software simulation of reference channel and/or reference receiver | |
US8265134B2 (en) | Multi-channel equalization to compensate for impairments introduced by interleaved devices | |
US8639112B2 (en) | Testing of transmitters for communication links by software simulation of reference channel and/or reference receiver | |
CN103259597B (zh) | 非线性补偿装置、方法和发射机 | |
JP2013074625A (ja) | 非線形推定装置、方法及び受信機 | |
WO2020238349A1 (zh) | 群时延波动校准方法、存储介质及电子装置 | |
Biondini | An introduction to rare event simulation and importance sampling | |
Vasylchenkova et al. | Fixed-point realisation of fast nonlinear Fourier transform algorithm for FPGA implementation of optical data processing | |
Span et al. | Successive eigenvalue removal for multi-soliton spectral amplitude estimation | |
Ye et al. | SPICE-compatible equivalent circuit models for accurate time-domain simulations of passive photonic integrated circuits | |
Kremp et al. | Fast split-step wavelet collocation method for WDM system parameter optimization | |
Fougstedt et al. | Power-efficient time-domain dispersion compensation using optimized FIR filter implementation | |
US9941963B2 (en) | Non-linear propagation impairment equalization | |
Arık et al. | Low-complexity implementation of convex optimization-based phase retrieval | |
Moore et al. | A method to compute statistics of large, noise-induced perturbations of nonlinear Schrödinger solitons | |
Huleihel et al. | Analysis of mismatched estimation errors using gradients of partition functions | |
US20110200090A1 (en) | Method for an equalizer computation in a media system using a data set separator sequence | |
Li et al. | Noise-induced perturbations of dispersion-managed solitons | |
Allen et al. | Unbiased Spectral Estimation and System Identification Using Short‐Time Spectral Analysis Methods | |
US7428116B2 (en) | Method for a cyclic dibit response estimation in a media system using a data set separator sequence | |
Medvedev et al. | Fast computation of the direct scattering transform by fourth order conservative multi-exponential scheme | |
Yang et al. | An integral split-step Fourier method for digital back propagation to compensate fiber nonlinearity | |
Premaratne | Split-step spline method for modeling optical fiber communications systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |