WO2014183699A1

WO2014183699A1 - 光相干通信中色散估计方法及装置

Info

Publication number: WO2014183699A1
Application number: PCT/CN2014/078327
Authority: WO
Inventors: 姚扬中; 李强; 李运鹏; 顾国华; 任薇
Original assignee: 中兴通讯股份有限公司
Priority date: 2013-10-22
Filing date: 2014-05-23
Publication date: 2014-11-20
Also published as: EP3048746B1; CN104579476A; CN104579476B; EP3048746A4; EP3048746A1; US20160269108A1; US9729232B2

Abstract

一种光相干通信中色散估计方法及装置，该方法包括：对IQ不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据；对所述频域数据求自相关序列，将所述自相关序列的值分别进行快速逆傅立叶变换，对快速逆傅立叶变换结果求模平方，并将两个偏振方向结果相加，得^Ρ ^[n]；确定多组数据的^P ^[n]平均值^Ρ ^̃[n]，求^Ρ ^̃[n]的最大值的索引n_o，根据最大值的索引n_o估算光纤链路色散值。上述技术方案能够非常精确地估算色散值，且估算速度快。

Description

光相干通信中色散估计方法及装置

技术领域

本发明涉及光相干通信中色散估计技术领域，尤其涉及一种光相干通信中色散估计方法及装置。

背景技术

随着互联网流量的增加，在干线系统的光通信系统中需要更大的容量。当每波长比特率增加时，在传输路径上的色度色散，偏振模色散以及各种非线性效应的波形失真而导致信息质量的退化变得艮严重。

和非相干技术相比，数字相干接收技术有如下优点：大约 3dB的光信噪比 (OSNR)增益；可以方便地釆用电均衡技术来应对信道变化，降低成本等；可以釆用更高效的调制技术以及偏振复用来提高传输容量。因此数字相干技术被认为是高速光通信系统的关键技术。

光相干接收机中，通过将信号光和本振光进行混频，信号光的幅度和相位信息被搬移到基带信号中，因而光相干检测保留了光场的所有信息，可以发挥数字信号处理技术的功能和性能上的优势。釆用电均衡技术能够近乎完全地补偿光信号的线性失真，如补偿色度色散 (CD), 偏振模色散 (PMD)等。

图 1 ( a )为典型的数字相干接收机框图。接收光信号被偏振分光器 (PBS) 分成相互正交的两个偏振态信号。 PBS输出偏振光信号通过 90。光混频器与本振光信号进行混频。混频后的光信号通过平衡光电检测器转换为基带电信号。光电转换后的电信号为每个偏振态有两路信号，但这 4路信号与原始的 4路信号并不对应，因为经过传输信道后，两个偏振态之间有串扰、偏振态也有旋转，所以接收端此处的两个偏振态、每个偏振态有两路正交信号与发射信号没有对应关系。光电转换后的电信号通过 ADC转换为数字信号。可以通用数字信号处理技术对 ADC转换为数字信号进行处理。

数字信号处理部分，如图 1 ( b )所示，包括 IQ Skew补偿， IQ不平衡补偿，色散 (CD)补偿，色散估计，时钟恢复，偏振解复用，载波恢复，判决检测等。 IQ Skew调整 4路信号的延迟，解决光电转换前端 4路信号延迟不一恢复调整釆样相位，为后面的偏振解复用模块提供稳定釆样相位的数据。偏振解复用补偿残余色散和偏振模色散 (PMD)。载波恢复，补偿接收端本振光与发送端载波光频率不一致问题。其中，色度色散 (CD)的值一般较大，色度色散和偏振模色散 (PMD)的均衡一般分两部分完成，首先，补偿静态色散，这里的均衡器通常不能使用标准的自适应算法进行更新系数，如补偿 40000ps/nm色度色散，则滤波器抽头数要达到大几百甚至上千，通常利用快速傅立叶变换技术进行频域快速卷积。由色散估计模块为色散补偿模块提供要补偿的色散值。

残余的色度色散和偏振模色散的补偿通过 FIR蝶形均衡器来实现， FIR 蝶形滤波器釆用自适应算法对系数进行更新，以跟踪补偿随时间动态变化的偏振模色散。 FIR蝶形均衡器又称为偏振解复用。 FIR蝶形自适应均衡器，有均衡、匹配滤波和釆样位置调整的作用。当釆样位置变化范围过大，或釆样频偏存在，釆样相位变化范围超过 FIR蝶形自适应均衡器调整范围，会引起 FIR蝶形自适应均衡器无法正常工作。需要在 FIR蝶形均衡器之前放一个时钟恢复模块。

时钟恢复，估计输入符号的釆样时间误差，并对符号的釆样时间进行插值调整，或者通过 VCO调整 ADC釆样频率，以保证提供稳定的符号釆样相位。时钟恢复的鉴相器需要能够容忍一定的信号失真，而传统的时钟恢复通常仅能容忍很小的色散。为了不增加时钟恢复模块的复杂性，且目前缺少容忍大残余色散值的时钟恢复方法，需要精确的色散补偿。

综上所述，相关技术中提出的利用传输质量信息（如差错率、 Q因数等）来进行可变色散补偿器的反馈控制的方法，以一定步长的色散间隔，改变色散补偿滤波器的色散补偿量，直到系统收敛。系统启动时，搜索过程是緩慢的，且精度低。再者，光纤链路受到环境温度变化的影响，链路色散值也会发生緩慢变化。此外，工作过程中，相关技术的基于搜索的色散估计方法，很难判断运行中色散值的緩慢变化。发明内容

本发明实施例要解决的技术问题是提出一种光相干通信中色散估计方法及装置，以提高色散测量精度和效率。

为了解决上述问题，釆用如下技术方案：

一种光相干通信中的色散估计方法，包括：

对 IQ 不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据；

对两个偏振方向的所述频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值；将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果；对两个快速逆傅立叶变换结果分别求模平方，得到两个模平方结果；将这两个模平方结果相加 , 得到模平方结果之和 ^{p w}；针对多个 IQ不平衡补偿后的数据，得到多个模平方结果之和，对所有的

模平方结果之和求平均，得到色散目标函 ; 求所述色散目标函数 ^Ρ[Ώ]的最大值的索引 ,根据该索引估算光纤链路色散值。可选地，对 IQ不平衡补偿后的数据做快速傅立叶变换，得到两个偏振方向的频域数据的步骤包括：

对 IQ不平衡补偿后的两个偏振数据，分别进行快速傅立叶变换后，按照如下方式得到两个偏振方向的频域数据； ί]和 ¾ , 其中， .k = it, ... , N - Ί , 为频率索引，为傅立叶变化点数：按照如下公式取两个偏振方向非正交的信号的频谱，得到两个偏振方向的频域数据：

i Ik] = X k cos 8； + Y[k] sin i - 1,2 其中， =。二^ 可选地，对 IQ不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据的步骤包括：

使用频域速卷积变换或者快速傅立叶变换变换得到两个偏振方向的频域数据。

可选地，两个偏振方向的所述频域数据的自相关序列间隔为波特率；对两个所述偏振方向的频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值的步骤包括：按照如下公式，对两个偏振方向的频域数据中的频语' ^Yl 求自相关序列对频谱 ^½ 求自相关序列

C_t [k] =

[k + fe_¾aud] k - 6_f ... , K i - 1,2 其中，不大于的整数, 为波特率大小对应频率索引

可选地，将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果的步骤包括：

按照如下公式，将两个所述自相关序列的值分别进行快速逆傅立叶变换:

其中 ^{3 Si}为傅立叶变化点数。可选地，将这两个模平方结果相加，得到模平方结果之和 ^[Ώ]的步骤包括: 按照如下公式得到所述模平方结果之和：

= 一 ― Ϊ

,其中 ' ^为傅立叶变化点数。可选地，求所述色散目标函数 ^P[n的最大值的索引^¾后，该方法还包括: 釆用插值公式对所述索引"'进行修正。可选地，釆用插值公式对' ^¾进行修正的步骤包括：按照如下抛物线插值公式对所述索引 ^¾进行修正，得到修正后的结果： r _ — 1】—P【？ ¾ + 1]

H '

· 2(Ρ [η_ϋ - 11 - 2P[a₀| - P[n_c + 1])· 其中， ^为修正后的结果。可选地，根据最大值的索引估算光纤链路色散值的步骤包括：按照如下公式计算所述光纤链路色散值 CD:

CD =

其中， c为真空中光速，为光载波在真空中的波长， ^ 为波特率。

一种光相干通信中色散估计装置，包括色散补偿单元和色散估计单元，其中：

所述色散补偿单元设置成：对 IQ进行不平衡补偿，对 IQ不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据；

所述色散估计单元包括第一模块、第二模块和第三模块，其中：所述第一模块设置成：对两个偏振方向的所述频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值；

所述第二模块设置成：将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果；对两个快速逆傅立叶变换结果分别求模平方，得到两个模平方结果；将这两个模平方结果相加，得到模平方结果

、 . FW

之和；

所述第三模块设置成：对针对多个 IQ不平衡补偿后的数据，得到的模平方结果之和求平均，得到色散目标函数 ; 求所述色散目标函数 ^的最大值的索引 , 根据该索引"'估算光纤链路色散值。

可选地，所述色散补偿单元设置成按照如下方式对 IQ不平衡补偿后的数据做快速傅立叶变换，得到两个偏振方向的频域数据：

对 IQ不平衡补偿后的两个偏振数据，分别进行快速傅立叶变换后，按照如下方式得到两个偏振方向的频域数据 ^. ]和 , 其中， k = 0, ... N_m _ 1 , ^Κ 为频率索引， ^m为傅立叶变化点数：按照如下公式取两个偏振方向非正交的信号的频谱，得到两个偏振方向的频域数据：

『fcl 二 X L'u + Y\k] sin ί - 1,2 其中，

可选地，所述第一模块设置成按照如下方式对两个偏振方向的所述频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值：按照如下公式，对两个偏振方向的频域数据中频谱 ^¾求自相关序列 l^M,对频谱' ^i¾求自相关序列 ^ί ,其中，频域数据的自相关序列间隔为波特率：

Q[¾ 十 k = , ί = 1,2 上式中，为不大于¹ ^— ^{1 _ h i:l£l}的整数，』为波特率大小对应频率索引间隔。可选地，所述第二模块设置成按照如下方式将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果：

按照如下公式，将所述自相关序列的值分别进行快速逆傅立叶变换：

S

Pt M = Xf₌₁ Q = -lV_{i ¾j} ... , _f¾— _ 1.2

其中^ ^为傅立叶变化点数。可选地，所述第二模块设置成按照如下方式将这两个模平方结果相加，得到模平方结果之和^¹¹¹：

按照如下公式对得到所述模平方结果之和

[n] = \Pjnl + ' = - 一 - 1 可选地，所述色散估计单元还包括修正模块，其中：所述修正模块设置成：釆用插值公式对所述索引进行修正。可选地，所述修正模块设置成按照如下抛物线插值公式对所述索引' ^¾进行修正：

P f?lg _ 1】— P【？ ¾ 4- Ϊ]

2( ?ΐ ₃ - 1| - 2Ρ[η₀| - P[j¾十 1]) 其中，为修正后的结果。可选地，所述第三模块设置成按照如下方式根据所述索引' 估算光纤链路色散值：按照如下公式计算所述光纤链路色散值：

其中， c为真空中光速，为光载波在真空中的波长，为波特率。可选地，所述色散补偿单元设置成按照如下方式对 IQ不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据：对 IQ 不平衡补偿后的数据使用频域速卷积变换得到两个偏振方向的频域数据。可选地，所述色散估计单元还包括：傅立叶变换模块，对 IQ不平衡补偿后的数据分别使用快速傅立叶变换得到两个偏振方向的频域数据。

上述技术方案能够直接计算出色散，无需搜索，这样系统启动时，可以快速估算到链路色散值。系统运行中，也可以继续估算链路色散值，跟踪其变化，为色散补偿模块精确补偿色散提供保证。可以实现精确的色散补偿，从而可以减小后面时钟恢复模块的复杂性等。综上所述，上述技术方案与相关技术相比，能够非常精确地估算色散值，且估算速度快。

附图概述

图 1 ( a )是相关技术中典型的数字相干接收机框图；

图 1 ( b )是图 1 ( a )所示数字相干接收机中数字信号处理原理图；图 2为本发明实施例的光相干通信中色散估计流程图；图 3 是本发明实施例的仿真^{13 [n}曲线图；图 4 是本发明实施例的仿真' ^曲线图，其中，横坐标用色散值表示；图 5 是本发明实施例的仿真' Η曲线的局部放大显示图；

图 6为本发明实施例的色散估计装置示意图。

本发明的较佳实施方式

下文将结合附图对本发明技术方案作进一步详细说明。需要说明的是，在不冲突的情况下，本申请的实施例和实施例中的特征可以任意相互组合。实施例 1

本发明实施例提供了一种光相干通信中色散估计方法，如图 2所示，包括如下步骤：

步骤 200,对 IQ不平衡补偿后的数据做快速傅立叶变换 FFT得到两个偏振方向的频域数据。

通常，色散补偿模块使用频域速卷积技术，可以为 CD估计模块提供频域数据。或者， CD估计模块单独对 IQ不平衡补偿后的时域信号做 FFT得频域数据。

对于偏振复用数字相干接收机，有两个正交的偏振态，分别计算得到频域信号

X 和 Υ , 。

τ其中，为频率索引，为傅立叶变化点数，频率间隔为 A^{/ :=} ' ^s/Nfft，其中为时域了数字信号的釆样间隔。

为了消除偏振模色散的影响，本发明实施例按照如下公式取两个偏振方向非正交的信号的频谱 = ] 咽 ^²。其中， ^¾:=

步骤 201, 对两个偏振方向的频域信号分别求自相关序列，得到两个偏振方向的频域信号相应的自相关序列的值，其中，计算自相关序列的间隔可为波特率。

该步骤中__; 按照如下公式对频语 ^[ 求自相关序列对频语求自相关序列 ^[¾]：

¾ Ik] = Xi Ik.4- k_baVi(t] — 0"„ i — 12 其中，为不大于 — ^{1_ aud}的整数， ^aud为波特率大小对应频率索引间隔。

步骤 202, 对步骤 201的两个自相关序列值分别进行快速逆傅立叶变换，得到两个逆傅里叶变换结果。其中' 为傅立叶变化点数。

步骤 203 , 对步骤 202的两个逆傅立叶变换结果分别求模平方，得到两个模平方结果，将这两个模平方结果相加，得到模平方结果之和 ^人

luj = \P, M\^L + l »JI^; ?Ϊ二 — - i 步骤 204, 将多组 IQ不平衡补偿后的数据执行上述步骤 201-203 , 得到多个模平方结果之和，对所有的模平方结果之和求平均，得到色散目标函数图 3所示即为所求得的 ^曲线。仿真条件为，偏振复用 DQPSK相干光通信系数，系统比特率为 112Gbit/s, 波特率为 28Gbaud, 光数字相干接收机的 ADC釆样速率为两倍符合速率。步骤 200的傅立叶变化点数为 4096, 步骤 203的逆傅立叶变换点数也为 4096。图 3中，离线仿真所用数据源的色度色散值为 48000ps/nm。

步骤 205 , 求色散目标函数^的最大值的索引 ^ί¾。

为了提高整数索引^¾ '的分辨率，在步骤 205的操作之后，还可以通过插值技术，对 ^£进行修正，例如，釆用如下抛物线插值公式对 ³进行修正：

, 一 [?½— 1】—P【？ ¾ + 1]

H '

· 2(Ρ [?ι₀ - 11 - 2Ρ[π_Ό] - P[n_c + 1]) 如图 3所示，离线仿真所用数据源的色度色散值为 48000ps/nm得到的曲线，求得整数索引为 605 , 经插值修正精确索引为 605.0695。

步骤 206, 根据色散目标函数' ^的最大值的索引估算光纤链路色散值。该步骤中，可通过下述公式估算出光纤链路的色散值 CD:

U rn ^fba A 其中， c为真空中光速，为光载波在真空中的波长， ^isa5id为波特率。由于的索引 n与色散有对应关系，将图 3 中的横坐标用色散表示，见图 4 所示。由图 4 可以看出，可以估计的色散值范围为 -162ns/nm 到 162ns/nm。将图 4局部放大后，得图 5。图 5中，索引 605对应的色散值为 47994ps/nm, 修正精确索引 605.0695对应的色散值为 48000ps/匪。说明本发明实施例色散估计精确度很高。

经过大量数据仿真，釆用上述技术方案进行色散估计误差为几个 ps/nm。而相关技术的色散估计方法，如通过自适应均衡器收敛的 MSE值为目标函数的色散搜索法，估计误差为几百 ps/nm。本申请技术方案的色散估计误差比相关技术的色散估计方法精度提高上百倍。

另外，上述技术方案中，可以直接计算出色散，无需搜索，故系统启动时，可以快速估算到链路色散值。系统运行中，也可以继续估算链路色散值，跟踪其变化，为色散补偿模块精确补偿色散提供保证。精确的色散补偿，可以减小后面时钟恢复模块的复杂性等。

实施例 2

本发明实施例还提供了一种光相干通信中色散估计装置，如图 6所示，包括至少一个用来执行色散补偿单元 601和色散估计单元 602的处理器，其中，

色散补偿单元 601设置成：对 IQ不平衡补偿后经过快速傅立叶变换得到的两个偏振方向的频域数据；

色散估计单元 602又划分为第一模块 6021、第二模块 6022和第三模块 6023 , 其中：

第一模块 6021设置成：对 IQ不平衡补偿后经过快速傅立叶变换得到的两个偏振方向的频域数据求自相关序列；

上述第一模块 6021 按照如下公式，对两个偏振方向的频域数据中频谱

'^¾'^{ι ί¾]}求自相关序列对频语求自相关序列^ ^],其中，频域数据的自相关序列间隔为波特率：

€dk] = i [k]X; [k + k_b:aiid] = ©,„., A' i = 12 上式中， ^A '为不大于'^¾— 的整数，为波特率大小对应频率索引间隔。

第二模块 6022设置成：将所述自相关序列的值分别进行快速逆傅立叶变换，对快速逆傅立叶变换结果求模平方，并将两个偏振方向结果相加，；本发明实施例中，第二模块 6022可以按照如下公式，将所述自相关序列的值分别进行 n快速逆傅立叶变换：

其中' 为傅立叶变化点数。再按照如下公 ^对所述快速逆傅立叶变换结果求模平方，并将两个偏振方向结果相加，得 ^ΡίΏ :

= "- ² + \ ₂ΜΪ¹ = — ¾赚'— - 1 第三模块 6023设置成：确定多组数据的 ^平均值求的最大值的索引 ^η 根据最大值的索引'' '估算光纤链路色散值。

上述第三模块 6023可以按照如下公式计算光纤链路色散值：

其中， c为真空中光速，为光载波在真空中的波长， ^isa5ia为波特率。由于的索引 n与色散有对应关系，将图 3 中的横坐标用色散表示，见图 4 所示。由图 4 可以看出，可以估计的色散值范围为 -162ns/nm 到 162ns/nm。将图 4局部放大后，得图 5。图 5中，索引 605对应的色散值为 47994ps/nm, 修正精确索引 605.0695对应的色散值为 48000ps/匪。说明本发明实施例色散估计精确度很高。

可选地，色散估计单元还可以包括修正模块 6024, 该修正模块 6024设置成：在求出的最大值的索引后，釆用插值公式对进行修正。例如，可以按照如下抛物线插值公式对所述^¾>进行修正：

还要说明的是，色散估计单元 602中第一模块 6021 , 对 IQ不平衡补偿后经过快速傅立叶变换得到的两个偏振方向的频域数据求自相关序列时， IQ 不平衡补偿后经过快速傅立叶变换得到的两个偏振方向的频域数据可以是色散补偿单元，对 IQ不平衡补偿后的数据使用频域速卷积变换得到的。也可以是由色散估计单元使用 FFT变换得到的。其中，由散估计单元使用 FFT变换得到两个偏振方向的频域数据时，可在散估计单元中增加傅立叶变换模块，以使用 FFT变换得到两个偏振方向的频域数据。

上述光相干通信中色散估计装置的其他细节可参见实施例 1 的相应内容，在此不再赘述。

本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序来指令相关硬件完成，所述程序可以存储于计算机可读存储介质中，如只读存储器、磁盘或光盘等。可选地，上述实施例的全部或部分步骤也可以使用一个或多个集成电路来实现。相应地，上述实施例中的各模块 /单元可以釆用硬件的形式实现，也可以釆用软件功能模块的形式实现。本申请不限制于任何特定形式的硬件和软件的结合。

以上所述，仅为本发明的较佳实例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

工业实用性上述技术方案能够直接计算出色散，无需搜索，这样系统启动时，可以快速估算到链路色散值。系统运行中，也可以继续估算链路色散值，跟踪其变化，为色散补偿模块精确补偿色散提供保证。可以实现精确的色散补偿，从而可以减小后面时钟恢复模块的复杂性等。综上所述，上述技术方案与相关技术相比，能够非常精确地估算色散值，且估算速度快。因此本发明具有 4艮强的工业实用性。

Claims

权利要求书

1、一种光相干通信中的色散估计方法，包括：

对两个偏振方向的所述频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值；将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果；对两个快速逆傅立叶变换结果分别求模平方，得到两个模平方结果；将这两个模平方结果相加，得到模平方结果之和 ^[Ώ]; 针对多个 IQ不平衡补偿后的数据，得到多个模平方结果之和，对所有的模平方结果之和求平均，得到色散目标函数 ^PM; 求所述色散目标函数 ^1]的最大值的索引 ^ί¾,根据该索引估算光纤链路色散值。

2、如权利要求 1所述的色散估计方法，其中，对 IQ不平衡补偿后的数据做快速傅立叶变换，得到两个偏振方向的频域数据的步骤包括：

对 IQ不平衡补偿后的两个偏振数据，分别进行快速傅立叶变换后，按照如下方式得到两个偏振方向的频域数据； q i和 ] , 其中， = o，， λ¾ _ 1 , 为频率索引， ^¾为傅立叶变化点数：按照如下公式取两个偏振方向非正交的信号的频谱，得到两个偏振方向的频域数据：

¾ [¾] = X [¾] cu Θ_ι + ¥\k] sin 6, i - 1 2 其中， ^{= ί} , ^{2 := ϊ}

3、如权利要求 2所述的色散估计方法，其中，对 IQ不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据的步骤包括：

4、如权利要求 3所述的色散估计方法，其中：

两个偏振方向的所述频域数据的自相关序列间隔为波特率；对两个所述偏振方向的频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值的步骤包括：按照如下公式，对两个偏振方向的频域数据中的频语' 求自相关序列 ^ClM,对频语 ^ί 求自相关序列 ^[ ：

其中， ^为不大于 ^¾— ^{1 _} ^ 的整数， ^¾5 为波特率大小对应频率索引间隔。

5、如权利要求 4所述的色散估计方法，其中，将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果的步骤包括：按照如下公式，将两个所述自相关序列的值分别进行快速逆傅立叶变换：

其中^ 为傅立叶变化点数。

6、如权利要求 5所述的色散估计方法，其中，将这两个模平方结果相加，得到模平方结果之和的步骤包括：

按照如下公式得到所述模平方结果之和：傅立叶变化点数。

7、如权利要求 6所述的色散估计方法，其中，求所述色散目标函数' ^ 的最大值的索引后，该方法还包括：釆用插值公式对所述索引 ^¾ '进行修正。

8、如权利要求 7所述的色散估计方法，其中，釆用插值公式对进行修正的步骤包括：按照如下抛物线插值公式对所述索引 ^¾进行修正，得到修正后的结果： t — ^: ^ -] ；

· 2(Ρ [π_α - 1】 - 2F[¾1 - Ρ[?Ϊ₀ + !]) _ 其中， ^为修正后的结果。

Ί、

9、如权利要求 8所述的色散估计方法，其中，根据最大值的索引 ' ⁰估算光纤链路色散值的步骤包括：按照如下公式计算所述光纤链路色散值 CD:

CD = ~ ^ —

其中， c为真空中光速，为光载波在真空中的波长， ^/¾sa^为波特率。

10、一种光相干通信中色散估计装置，包括色散补偿单元和色散估计单元，其中：

所述色散估计单元包括第一模块、第二模块和第三模块，其中：所述第一模块设置成：对两个偏振方向的所述频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值；所述第二模块设置成：将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果；对两个快速逆傅立叶变换结果分别求得到两个模平方结果；将这两个模平方结果相加，得到模平方结果

所述第三模块设置成：对针对多个 IQ不平衡补偿后的数据，得到的模平方结果之和求平均，得到色散目标函数' 求所述色散目标函数々最大值的索引 ' 根据该索引估算光纤链路色散值。

11、如权利要求 10所述的色散估计装置，其中，所述色散补偿单元设置成按照如下方式对 IQ不平衡补偿后的数据做快速傅立叶变换，得到两个偏振方向的频域数据：

θ, = 0

其中

12、如权利要求 11所述的色散估计装置，其中，所述第一模块设置成按照如下方式对两个偏振方向的所述频域数据分别求自相关序列，得到两个偏振方向的频域数据相应的两个自相关序列的值：按照如下公式，对两个偏振方向的频域数据中频谱求自相关序列

^ClM,对频谱 ^ί 求自相关序列其中，频域数据的自相关序列间隔为波特率：上式中， ^为不大于 ^¾— ^{1 _ ¾aU} 整数， ^¾aud为波特率大小对应频率索引间隔。

13、如权利要求 12所述的色散估计装置，其中，所述第二模块设置成按照如下方式将两个所述自相关序列的值分别进行快速逆傅立叶变换，得到两个逆傅立叶变换结果：

按照如下公式，将所述自相关序列的值分别进行快速逆傅立叶变换： Pi [ϊΐ] =∑f₌₁ Q [feje⁵²⁵¹^^ _?l = -N^ ^ ^— 1 i - 1.2

其中为傅立叶变化点数。

14、如权利要求 13所述的色散估计装置，其中，所述第二模块设置成按照如下方式将这两个模平方结果相加，得到模平方结果之和

按照如下公式对得到所述模平方结果之和

15、如权利要求 14所述的色散估计装置，其中，所述色散估计单元还包括修正模块，其中：所述修正模块设置成：釆用插值公式对所述索引进行修正。

16、如权利要求 15所述的色散估计装置，其中，所述修正模块设置成按照如下抛物线插值公式对所述索引 ' 进行修正： _ n₀—1】— P【3¾ + ί]

其中， ^为修正后的结果。

17、如权利要求 16所述的色散估计装置，其中，所述第三模块设置成按照如下方式根据所述索引 '"估算光纤链路色散值：按照如下公式计算所述光纤链路色散值：

其中， c为真空中光速，为光载波在真空中的波长， ^为波特率。

18、如权利要求 10至 17中任一项所述的装置，其中，所述色散补偿单元设置成按照如下方式对 IQ 不平衡补偿后的数据做快速傅立叶变换得到两个偏振方向的频域数据：对 IQ 不平衡补偿后的数据使用频域速卷积变换得到两个偏振方向的频域数据。

19、如权利要求 10至 17中任一项所述的色散估计装置，其中，所述色散估计单元还包括：傅立叶变换模块，对 IQ不平衡补偿后的数据分别使用快速傅立叶变换得到两个偏振方向的频域数据。