WO2014019211A1

WO2014019211A1 - 色散估计的方法、装置和系统

Info

Publication number: WO2014019211A1
Application number: PCT/CN2012/079619
Authority: WO
Inventors: 隋琪; 刘伯涛; 吕超; 李良川
Original assignee: Huawei Technologies Co Ltd; Hong Kong Polytechnic University HKPU
Current assignee: Huawei Technologies Co Ltd; Hong Kong Polytechnic University HKPU
Priority date: 2012-08-03
Filing date: 2012-08-03
Publication date: 2014-02-06
Anticipated expiration: 2015-02-03
Also published as: CN103004110A; CN103004110B

Abstract

本发明实施例提供了一种色散估计的方法、装置和系统，该方法包括：计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列；根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的位置相对应；根据该脉沖位置，确定该目标光信号的色散值。根据本发明实施例的色散估计的方法、装置和系统，通过获取目标光信号的功率自相关函数，并确定该功率自相关函数出现脉沖峰值的位置，能够在无需扫描所有可能的色散值的情况下，根据该出现脉沖峰值的位置确定目标光信号的色散值，到达快速估计信号的色散值的目的。

Description

色散估计的方法、装置和系统技术领域

本发明涉及通信领域，并且更具体地，涉及色散估计的方法、装置和系统。背景技术

目前，已知一种色散估计的方法，即，对于未知色散的信号，以一定的步长间隔，扫描所有可能的色散值，并对于每个色散值通过设计代价函数（也可以称为算法设计）计算出对应的函数值，最后通过比较所有计算出的函数值（最大值或最小值），找到实际的色散值，实现色散估计。

但是，由于该方法需要扫描所有的色散值以进行计算，因此需要消耗大量的时间，不能快速估计信号的色散值。发明内容

本发明实施例提供一种，能够快速估计信号的色散值。

第一方面，提供了一种色散估计的方法，该方法包括：计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列；根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的位置相对应；根据该脉沖位置，确定该目标光信号的色散值。

在一种可能的实施方式中，在该目标光信号为非归零码格式信号的情况下，在根据目标光信号的信号功率的自相关函数，确定自相关函数序列前，该方法还包括：对该目标光信号进行第一处理，以消除所述目标光信号的零频干扰。

结合第一方面和第一种可能的实施方式，在第二种可能的实施方式中，在根据目标光信号的信号功率的自相关函数，确定自相关函数序列前，该方法还包括：对该目标光信号进行第二处理，以消除该目标光信号的直流分量干扰。

结合第一方面、第一种可能的实施方式和第二种可能的实施方式，在第三种可能的实施方式中，该根据该自相关函数序列，确定脉沖位置，包括：根据该自相关函数序列，确定第一自相关函数序列，该第一自相关函数序列包括该脉沖峰值；根据该第一自相关函数序列，确定该脉沖位置。

结合第一方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的实施方式，在第四种可能的实施方式中，该根据该第一自相关函数序列，确定该脉沖位置，包括：根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；根据该第一滑动平均处理后的第一自相关函数序歹 ij , 确定该脉沖位置。

结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式和第四种可能的实施方式，在第五种可能的实施方式中，该根据该第一滑动平均处理后的第一自相关函数序列，确定该脉沖位置，包括：根据第一序号和第一预设值，确定第二窗口，该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；根据该第一自相关函数序列和该第二窗口，确定第二自相关函数序列；根据该第二窗口，对该第二自相关函数序列进行第二滑动平均处理；根据该第二滑动平均处理后的第二自相关函数序列，确定该脉沖位置。

结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式和第五种可能的实施方式，在第六种可能的实施方式中，该根据第一序号和第一预设值，确定第二窗口，包括：根据第一序号和第一预设值，从该第一自相关函数序列中确定第一范围，该第一范围以该第一序号为中心，该第一范围的大小与该第一预设值相对应；根据该第一范围中数值最大的第一元素和第二预设值，从该第一自相关函数序列中确定第二范围，该第二范围以该第一元素的序号为中心，该第二范围的大小与该第二预设值相对应；根据该第一元素和该第二范围中数值第二大的第二元素，确定该第二窗口。

结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式，在第七种可能的实施方式中，该根据该脉沖位置，确定该目括：根据以下公式确定该目标光信号的色散值，

其中， CZ)表示色散值， n表示该脉沖位置在该自相关函数序列中的序号, Γ。表示该目标光信号的初始发射机光脉沖宽度， c表示真空光速， Λ表示该目标光信号的波长， SamplingRate表示模数转换器的采样率， SymbolRate表示该目标光信号的符号速率。

第二方面，提供了一种色散估计的装置，该装置包括：功率自相关处理单元，用于计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列，以及向脉沖位置搜索单元传输指示该自相关函数序列的信息；脉沖位置搜索单元，用于从该功率自相关处理单元获取指示该自相关函数序列的信息，并根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的位置相对应，以及向色散值确定单元传输指示该脉沖位置的信息；色散值确定单元，用于从该脉沖位置搜索单元获取指示该脉沖位置的信息，并根据该脉沖位置，确定该目标光信号的色散值。

在一种可能的实施方式中，在该目标光信号为非归零码格式信号的情况下，该装置还包括：第一干扰处理单元，用于对该目标光信号进行第一处理，以消除所述目标光信号的零频干扰。

结合第二方面和第一种可能的实施方式，在第二种可能的实施方式中，该装置还包括：第二干扰处理单元，用于对该目标光信号进行第二处理，以消除该目标光信号的直流分量干扰。

结合第二方面、第一种可能的实施方式和第二种可能的实施方式，在第三种可能的实施方式中，该脉沖位置搜索单元具体用于根据该自相关函数序歹l , 确定第一自相关函数序列，该第一自相关函数序列包括该脉沖峰值；用于根据该第一自相关函数序列，确定该脉沖位置。

结合第二方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的实施方式，在第四种可能的实施方式中，该脉沖位置搜索单元具体用于根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；用于根据该第一滑动平均处理后的第一自相关函数序列，确定该脉沖位置。

结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式和第四种可能的实施方式，在第五种可能的实施方式中，该脉沖位置搜索单元具体用于根据第一序号和第一预设值，确定第二窗口，该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；用于根据该第一自相关函数序列和该第二窗口，确定第二自相关函数序列；用于根据该第二窗口，对该第二自相关函数序列进行第二滑动平均处理；用于根据该第二滑动平均处理后的第二自相关函数序列，确定该脉沖位置。

结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式和第五种可能的实施方式，在第六种可能的实施方式中，该脉沖位置搜索单元具体用于根据第一序号和第一预设值，从该第一自相关函数序列中确定第一范围，该第一范围以该第一序号为中心，该第一范围的大小与该第一预设值相对应；用于根据该第一范围中数值最大的第一元素和第二预设值，从该第一自相关函数序列中确定第二范围，该第二范围以该第一元素的序号为中心，该第二范围的大小与该第二预设值相对应；用于根据该第一元素和该第二范围中数值第二大的第二元素，确定该第二窗口。

结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式，在第七种可能的实施方式中，该色散值确定单元具体用于根的色散值，

其中， CZ)表示色散值， n表示该脉沖位置在该自相关函数序列中的序号， Γ。表示该目标光信号的初始发射机光脉沖宽度， c表示真空光速，表示该目标光信号的波长， Sam ^Rafe表示模数转换器的采样率， SymbolRate表示该目标光信号的符号速率。

第三方面，提供了一种色散估计的系统，该系统包括：总线；与该总线相连的处理器；与该总线相连的存储器；其中，该处理器通过该总线，调用该存储器中存储的程序，以用于计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列；根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的出现位置相对应；根据该脉沖位置，确定该目标光信号的色散值。

在一种可能的实施方式中，在该目标光信号为非归零码格式信号的情况下，该处理器通过该总线，调用该存储器中存储的程序，以用于对该目标光信号进行第一处理，以消除所述目标光信号的零频干扰。结合第三方面和第一种可能的实施方式，在第二种可能的实施方式中，该处理器通过该总线，调用该存储器中存储的程序，以用于对该目标光信号进行第二处理，以消除该目标光信号的直流分量干扰。

结合第三方面、第一种可能的实施方式和第二种可能的实施方式，在第三种可能的实施方式中，该处理器通过该总线，调用该存储器中存储的程序，以具体用于根据该自相关函数序列，确定第一自相关函数序列，该第一自相关函数序列包括该脉沖峰值；根据该第一自相关函数序列，确定该脉沖位置。

结合第三方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的实施方式，在第四种可能的实施方式中，该处理器通过该总线，调用该存储器中存储的程序，以具体用于根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；用于根据该第一滑动平均处理后的第一自相关函数序列，确定该脉沖位置。

结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式和第四种可能的实施方式，在第五种可能的实施方式中，该处理器通过该总线，调用该存储器中存储的程序，以具体用于根据第一序号和第一预设值，确定第二窗口，该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；用于根据该第一自相关函数序列和该第二窗口，确定第二自相关函数序列；用于根据该第二窗口，对该第二自相关函数序列进行第二滑动平均处理；用于根据该第二滑动平均处理后的第二自相关函数序列，确定该脉沖位置。

结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式和第五种可能的实施方式，在第六种可能的实施方式中，该处理器通过该总线，调用该存储器中存储的程序，以具体用于根据第一序号和第一预设值，从该第一自相关函数序列中确定第一范围，该第一范围以该第一序号为中心，该第一范围的大小与该第一预设值相对应；用于根据该第一范围中数值最大的第一元素和第二预设值，从该第一自相关函数序列中确定第二范围，该第二范围以该第一元素的序号为中心，该第二范围的大小与该第二预设值相对应；用于根据该第一元素和该第二范围中数值第二大的第二元素，确定该第二窗口。

结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式，在第七种可能的实施方式中，该处理器通过该总线，调用该存储器中存储的程序，以具体用于根据以下公式确定该目标光信号的色散值

其中， CZ)表示色散值， n表示该色散出现位置在该自相关函数序列中的序号， Γ。表示该目标光信号的初始发射机光脉沖宽度， c表示真空光速，表示该目标光信号的波长， Sam ^Rafe表示模数转换器的采样率， SymbolRate表示该目标光信号的符号速率。

根据本发明实施例的色散估计的方法、装置和系统，通过获取目标光信号的功率自相关函数，并确定该功率自相关函数出现脉沖峰值的位置，能够在无需扫描所有可能的色散值的情况下，根据该出现脉沖峰值的位置确定目标光信号的色散值，到达快速估计信号的色散值的目的。附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作筒单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图 1是根据本发明实施例的色散估计的方法的示意性流程图。

图 2是根据本发明实施例的色散估计的装置的示意性框图。

图 3是根据本发明实施例的色散估计的方法的示意性框图。具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

图 1示出了根据本发明实施例的色散估计的方法 100的示意性流程图，如图 1所示，该方法 100包括： S110, 计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列；

S120, 根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的位置相对应；

S130, 根据该脉沖位置，确定该目标光信号的色散值。

在本发明实施例中，通过试验（例如，获取多个光信号的功率自相关函数，并获取该多个光信号的色散值）发现，在光信号（例如，单偏振的相移键控法（PSK, Phase Shift Keying )信号、偏振分割复用（PDM, Polarization Division Multiplexing ) 的 PSK信号或正交振幅调制（ QAM , Quadrature Amplitude Modulation )信号 )存在色散时，该光信号的功率自相关函数会出现一个脉沖峰值，并且，不同的脉沖峰值出现位置与不同的色散值之间存在映射对应关系，因此，在本发明实施例中，可以根据光信号的信号功率的自相关函数中脉沖峰值的位置，确定该光信号的色散值。

因此，在 S110 中，可以计算目标光信号的信号功率（光功率）的自相关函数。在本发明实施例中，目标光信号为需要估算色散的光信号。在本发明实施例中，可以对该目标光信号进行以下处理，以单偏振系统为例，在相干接收机中，目标光信号和本振光在混频器（Hybrid ) 中混合，通过平衡接收机后转化为同向（in-phase ) 电流和正交电流（quadrature ) 两路电流，经过两路高速模数转换器（ADC, Analog to Digital Converter )采样后，转换为数字信号，得到 I信号（同向信号）和 Q信号（正交信号），从而，可以根据该 I信号和 Q信号来确定目标光信号的信号功率的自相关函数。以下，作为示例而非限定，以该通过该 I信号和 Q信号确定目标光信号的信号功率的自相关函数为例进行说明。

可选地，在本发明实施例中，在该目标光信号为非归零码格式信号的情况下，在根据目标光信号的信号功率的自相关函数，确定自相关函数序列前，该方法还包括：

对该目标光信号进行第一处理，以消除所述目标光信号的零频干扰。具体地说，零频干扰是指，在该目标光信号为非归零（NRZ, No Return

Zero )码格式的情况下，如果直接获取该格式的信号的自相关函数，则获取的自相关函数的零频位置会出现一个较大的峰值，从而对后述搜索脉沖峰值造成影响，因此，作为示例而非限定，在本发明实施例中，可以对 I信号和 Q信号进行延迟处理和相减处理，以实现对该目标光信号的高通滤波，其中，该延迟处理使用的迟延大小可以为例如，一个符号周期。从而，能够避免获取的自相关函数在零频位置出现较大的峰值。

可选地，在本发明实施例中，在根据目标光信号的信号功率的自相关函数，确定自相关函数序列前，该方法还包括：

对该目标光信号进行第二处理，以消除该目标光信号的直流分量干扰。具体地说，作为示例而非限定，可以对 I信号和 Q信号（具体地说，是与该 I信号相对应的 I信号序列以及与该 Q信号相对应的 Q信号序列）进行求平方和处理，然后用经过该求平方和处理后的信号减去信号均值，其中，作为示例而非限定，该信号均值可以通过以下方法确定，即，可以获取与该 I信号相对应的 I信号序列，以及与该 Q信号相对应的 Q信号序列，其后，对 I信号序列和 Q信号序列求和，然后除以 I信号序列和 Q信号序列的长度，从而确定该信号均值。通过以上处理，能够实现对该目标光信号的隔直处理，从而能够去掉直流分量产生的无用的自相关函数部分。

其后，在 S110中，可以对经过上述处理的目标光信号进行傅里叶变换，再做幅度平方，最后进行逆傅里叶变换，得到该目标光信号的信号功率的自相关函数，进而可以根据该自相关函数，生成以该自相关函数的各函数值为元素的自相关函数序列（以下，为了便于说明，记作 [«] )。并且自相关函数以表示为：

其中，表示目标光信号的自相关函数， FFT表示快速傅里叶变换， IFFT表示快速反傅里叶变换， y W表示目标光信号的信号功率（光功率），

R^]为根据该自相关函数（ w )生成的函数序列（自相关函数序列）。

在 S120中，可以从该 W中搜索脉沖峰值的位置（脉沖位置），即，确定该脉沖峰值在 W中的序号。

可选地，在本发明实施例中，可以搜索 W中数值最大的元素（记作元素 a ), 并将该元素 a的位置确定为脉沖峰值的位置，并可以记录该元素 a的序号，以在后述确定色散值时，根据该脉沖峰值在 W中的位置，或者说元素 a在 W中的序号（记作 nj, 确定色散值（随后进行详细说明）。

可选地，在本发明实施例中，该根据该自相关函数序列，确定脉沖位置，包括：根据该自相关函数序列，确定第一自相关函数序列，该第一自相关函数序列包括与该脉沖峰值相对应的元素；

根据该第一自相关函数序列，确定该脉沖位置。

具体地说，在本发明实施例中，可以从该 w中选取一个区间作为搜索范围（即，第一自相关函数序列，以下，为了便于说明，记作 ?₂W ), 并且，在本发明实施例中，可以将 W作为 W₂W, 也可以将 {WJNJ

中从序号为的元素到序号为 N₂的元素）作为 ?₂ [/i] , 其中，表示 W中选取的开始位置（序号）， N₂表示中选取的结束位置（序号）。

可选地，在本发明实施例中，该根据该自相关函数序列，确定第一自相关函数序列，包括：

根据用于获取所述目标光信号的相干接收机的电带宽和用于获取所述目标光信号的高速模数转换器的采样率，从所述自相关函数序列中确定选取所述第一自相关函数序列的开始位置。

具体地说，可以根据以下式 1确定 N 即：

Ni>floor ( 2 x SamplingRate/B_e ) 式 1

其中， floor 为向下取整数函数， Be表示接收机电带宽， SamplingRate 是表示 ADC采样率，通常情况下， Be>10GHz, SamplingRate>20G Sa/s。

并且，在本发明实施例中， N₂可以表示能够进行色散估计的最大范围，在本发明实施例中， N₂>Ni。

可选地，在本发明实施例中，该根据该第一自相关函数序列，确定该脉沖位置，包括：

根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；根据该第一滑动平均处理后的第一自相关函数序列，确定该脉沖位置。具体地说，在确定 ?₂W之后，对 ?₂W进行滑动窗口平均处理（第一滑动平均处理），并且，作为示例而非限定，该滑动窗口平均处理的窗口宽度（记作 widthj可以为例如， 3。在本发明实施例中，滑动窗口平均处理的方法与现有技术相同，这里，为了避免赘述，省略其说明。在对 ?₂[/ι]进行滑动窗口平均处理后（以下，为了便于理解，将处理后生成的序列记作 R₃ [«] ), 可以搜索 R₃W中数值最大的元素（记作元素 b ), 并将该元素 b的位置确定为脉沖峰值的位置，并可以记录该元素 b在 [«]中的序号（记作 n₂ ), 从而，可以才据以下式 2确定即 ni = n₂ + Ni - 1 式 2

通过对该 ?₂W (第一自相关函数序列）进行滑动窗口平均处理，能够减少自相关函数中噪声的干扰，使确定的脉沖峰值的位置（脉沖位置）更加准确。以上列举的作为 width 的数值仅为本发明的一个示例性说明，并不限定于此。

可选地，在本发明实施例中，在确定 ?₂W之后，也可以直接搜索 ?₂W中数值最大的元素（记作元素 C ) ,并将该元素 C的位置确定为脉沖峰值的位置，并可以记录该元素 c在中的序号（记作 n₃ ), 从而，可以根据以下式 3 确定 n 即

ii! = n₃ + Ni - 1 式 3

从而，在后述确定色散值时，根据该脉沖峰值在 W中的位置，或者说

11!的值，确定色散值 (随后进行详细说明）。

可选地，在本发明实施例中，该根据该第一滑动平均处理后的第一自相关函数序列，确定该脉沖位置，包括：

根据第一序号和第一预设值，确定第二窗口，该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；

根据该第一自相关函数序列和该第二窗口，确定第二自相关函数序列；根据该第二窗口，对该第二自相关函数序列进行第二滑动平均处理；根据该第二滑动平均处理后的第二自相关函数序列，确定该脉沖位置。在本发明实施例中，可以将自相关函数中与产生的色散相对应的脉沖峰值的全宽半高（FWHM, Full Width Half Maximum ) 的一半，即，将脉沖峰值（最大值）的二分之一作为门限（后述第一门限），寻找脉沖宽度，再将脉沖宽度的例如，二分之一作为第二窗口，具体地说，在确定 n₂后，可以从 W₂W中，在以 n₂为中心的预设范围（与第一预设值相对应，例如，可以为 5 个单位，这里， "单位" 可以是单位序号）内，搜索数值最大的元素（记作元素 d ),并根据该元素 d的值，以确定第一门限，这里，作为示例而非限定，将该元素 d的值的二分之一作为该第一门限，其后，以元素 d为中心，向左侧搜索第一个数值小于该第一门限的元素（记作元素 e ), 并记录该元素 e的序号（记作 ^ );向右侧搜索第一个数值小于该第一门限的元素（记作元素 f ), 并记录该元素 f的序号（记作 t₂ ),可以根据以下式 4,确定滑动窗口宽度（第二窗口的一例，记作 width₂ ), 即 width₂ = floor { ( t₂ - ti ) /2} - 1 式 4

其中， floor为向下取整数函数。应理解，以上列举的以脉沖峰值（最大值）的二分之一作为门限，确定第二窗口的方式仅为本发明的一个示例性说明，本发明并不限定于此。

在确定了 width₂ (第二窗口）后，可以从该 ?₂W中选取一个区间作为搜索范围（即，第二自相关函数序列，以下，为了便于说明，记作 ?₄W), 在本发明实施例中，可以将 w₂w作为 w₄w, 也可以将 {R₂[N₃]:R₂[N₄]} (即， w₂w 中从序号为 N₃的元素到序号为 N₄的元素）作为 ?₄ ], 其中， N₃表示中选取的开始位置（序号）， N₄表示 ?₂[/ι]中选取的结束位置（序号）。

具体地说，可以根据以下式 5确定 N₃, 即：

N₃ = ti - width₂ - 1 式 5

可以根据以下式 6确定 N₄, 即：

N₄ = t₂ + width₂ +1 式 6

在确定了 W后，根据 width₂, 对该 ?₄W进行滑动平均处理（第二滑动平均处理），以下，为了便于理解，将处理后生成的序列记作 R₅[«] , 可以搜索 R₅w中数值最大的元素（记作元素 g ),并将该元素 g的位置确定为脉沖峰值的位置，并可以记录该元素 g在 [«]中的序号（记作 n₄), 由于自相关函数第一位对应 τ=0, R [W]对应延迟为 (n— i)'SamplingInterval (修正一位），其中， Samplinglnterval表示采样间隔，从而，可以才艮据以下式 7确定 11 , 即

ii! = n₄ + tj - width₂ + Νι - 3 式 7

其中， width₂ = floor{ ( t₂— ） /2}— 1。

11!的值，确定色散值 (随后进行详细说明）。

通过对该 ?₂W (第一自相关函数序列）进行滑动窗口平均处理，能够减少自相关函数中噪声的干扰，使确定的脉沖峰值的位置（脉沖位置）更加准确。

可选地，在本发明实施例中，该根据第一序号和第一预设值，确定第二窗口，包括：

根据第一序号和第一预设值，从该第一自相关函数序列中确定第一范围，该第一范围以该第一序号为中心，该第一范围的大小与该第一预设值相对应；根据该第一范围中数值最大的第一元素和第二预设值，从该第一自相关函数序列中确定第二范围，该第二范围以该第一元素的序号为中心，该第二范围的大小与该第二预设值相对应；

根据该第一元素和该第二范围中数值第二大的第二元素，确定该第二窗口。

具体地说，如果偏振模色散 ( PMD, Polarization Mode Dispersion )较大，则可能会出现一个自相关函数中的脉沖（峰值）分裂为两个脉沖（峰值）的情况，此情况下，如果仅以脉沖峰值（最大值）的位置为中心寻找脉沖宽度，则寻找到的脉沖宽度值会产生偏差，因此，在本发明实施例中，在确定 n₂ 后，还可以从 W₂W中，在以 n₂为中心的预设范围（即，第一范围，与第一预设值相对应，例如，可以为 5个单位，这里为 "单位" 可以是单位序号）内，搜索数值最大的元素（即，元素 d )。

其后，在以元素 d在 ?₂W中的序号为中心的预设范围（即，第二范围，与第二预设值相对应，例如，可以为 10个单位，这里为 "单位" 可以为单位序号），搜索数值第二大的元素（记作元素 h )。

其后，确定元素 d与元素 h在中的位置关系，这里，为了便于理解，以元素 d在元素 h的右侧为例（即，元素 d在 ?₂ W中的序号大于元素 h在 ?₂ W 中的序号）进行说明。

其后， ^据该元素 d的值，以确定第一门限，这里，作为示例而非限定，将该元素 d的值的二分之一作为该第一门限。

其后，以元素 h为中心，向左侧搜索第一个数值小于该第一门限的元素 (记作元素 i, 并记录该元素 i的序号（记作 t₃ ); 以元素 d为中心，向右侧搜索第一个数值小于该第一门限的元素（记作元素 j ), 并记录该元素 j的序号（记作 t₄ ), 可以 ^据以下式 8, 确定一个中间参量（记作 width₃ ), 即 Width₃ = floor { ( t₄ - 1₃ ) /2} - 1 式 8

在确定了 width₃后，可以从该 W₂W中选取一个区间作为搜索范围（即，第二自相关函数序列，以下，为了便于说明，记作在本发明实施例中，可以将 w₂w作为 R₆[«] , 也可以将 {R₂[N₅] : R₂[N₆]} (即， w₂w中从序号为 N₅的元素到序号为 N₆的元素）作为 R₆[«], 其中， N₅表示?₂W中选取的开始位置（序号）， N₆表示中选取的结束位置（序号）。

具体地说，可以根据以下式 9确定 N₅, 即：

N₅ = t₃ - width₃ - 1 式 9

可以根据以下式 10确定 N₆, 即：

N₆ = t₄ + width₃ + 1 式 10

在确定了 R₆W后，可以根据以下式 11,确定对 R₆[«]进行滑动平均处理时使用的窗口宽度（第二窗口的另一例，记作 width₄), 即

width₄ = 2 X width₃ +1 式 11

其后，对该 R₆[«]进行滑动平均处理（第二滑动平均处理），以下，为了便于理解，将处理后生成的序列记作 R₇[«], 可以从该 R₇[«]中选取一个区间作为搜索范围（以下，为了便于说明，记作 R₈[«])，在本发明实施例中，可以将 R₇W作为 R₈W,也可以将 {R₇[N₇]:R₇[N₈]} (即， R₇W中从序号为 N₇的元素到序号为 N₈的元素）作为 R₈[«], 其中， N₇表示 R₇[_w]中选取的开始位置（序号）， N₈表示 R₇[_w]中选取的结束位置（序号）。

具体地说，可以根据以下式 12确定 N₇, 即：

N₇ = 2 式 12

可以 ^据以下式 13确定 N₈, 即：

N₈ = t₄ - t₃ + 2 X width₃ + 2 式 13

在确定 R₈W后，可以搜索 R₈W中数值最大的元素（记作元素 k),并记录该元素 k在 R₈[«]中的序号（记作 n₅ )。

可选地，在本发明实施例中，可以将该元素 k的位置确定为脉沖峰值的位置，由于自相关函数第一位对应 τ=0 , R [η]对应延迟为 - 1) · Samplinglnterval (修正一位），其中， Sam /^/wfervaZ表示采样间隔，从而，可以根据以下式 14确定即

ii! = n₅ + t₃ - width₃ + Ni - 3 式 14

其中， Width₃ = floor{ ( t₄ - t₃ ) /2} - 1

11!的值，确定色散值 (随后进行详细说明）。

可选地，在本发明实施例中，在确定 n₅后，还可以利用 R₈[«₅-l]、 R,[n₅],

R₈[«₅+l]三点做抛物线拟合，从而，可以根据以下式 15, 确定该抛物线顶点在 R₈W中的位置（序号，记作 n₆), 即 n - „ I R₈["₅- 1]- R₈["₅+l] 式

⁵ 2(R_&[n₅ -i] + R_&[n₅ +i]-2R_&[n₅])

其后，可以根据 n₆和以下式 16,确定该抛物线顶点在 ?₂W中的位置（序号），即

n₇ = n_{5 +} t₃ - width₃ - 2 + R₈["₅ _1]₇R₈["_{5 +}1] 式 ₁₆ 由于自相关函数第一位对应 τ=0 , R [η]对应延迟为 - 1) · Samplinglnterval (修正一位），其中， Sam /^toervaZ表示采样间隔，从而，可以根据以下式 17确定即

II！ =n₅.t₃ - width₃ + ( R,[n₅ -l]-R,[n +\] _{+ N} 3 式 ₁₇ 其中， width₃ = floor{ ( t₄— t₃ ) /2}— 1。

从而，在 S130中，根据在 S120中确定的自相关函数中的脉沖峰值出现位置（即， nj, 确定目标光信号的色散值。

在本发明实施例中，根据该脉沖位置，确定该目标光信号的色散值，包括：

根据离线数据，确定指示脉沖峰值位置与色散值的对应关系的指示信息，

根据该脉沖位置和该指示信息，确定该光信号的色散值。

具体地说，作为示例而非限定，例如，可以分别对多个光信号（离线数据的一例 )进行色散值测试等试验，记录各光信号的色散值和它们各自对应的自相关函数中脉沖峰值出现的位置，作为训练数据，从而，可以采用例如期望最大化算法，通过逐步改进公式的参数，使参数和训练数据的似然概率逐渐增大，最后终止于一个极大点。直观地理解，期望最大化算法也可被看作为一个逐次逼近算法：事先并不知道模型的参数，可以随机的选择一套参数或者事先粗略地给定某个初始参数 λθ , 确定出对应于这组参数的最可能的状态，计算每个训练样本的可能结果的概率，在当前的状态下再由样本对参数修正，重新估计参数 λ , 并在新的参数下重新确定模型的状态，这样，通过多次的迭代，循环直至某个收敛条件满足为止，就可以使得模型的参数逐渐逼近真实参数。因此，该步骤可以以离线的方式自动进行。将离线计算得到的上述各种脉沖峰值位置与各自对应的色散值预先存储起来，例如以表的形式存储起来。应理解，以上列举的各公式以及各公式中的常数的训练方法仅是本发明的一个实施例，本发明并不限定于此。

作为示例而非限定，以下式 18示出了色散值与自相关函数中脉沖峰值出现的位置的对应关系的一例，即

CD = (x + ^x² - 4Τ₀ ⁴ )- π - -^ 式 18

2π · SamplingRate · SymbolRate

其中， CD表示色散值， n表示所述脉沖位置在该自相关函数序列中的序号， Γ。表示所述目标光信号的初始发射机光脉沖宽度， c表示真空光速，表示所述目标光信号的波长， SampZ ^Rafe表示模数转换器的采样率， SymbolRate 表示所述目标光信号的符号速率。

应理解，以上列举的各公式、参数以及数值仅为本发明的示例性说明，本发明并不限定于此，其他能够用于确定光信号的自相关函数的脉沖峰值出现位置的各公式、参数以及数值均落入本发明的保护范围内。

根据本发明实施例的色散估计的方法，通过获取目标光信号的光信号功率自相关函数，并确定该光信号功率自相关函数出现脉沖峰值的位置，能够在无需扫描所有可能的色散值的情况下，根据该出现脉沖峰值的位置确定目标光信号的色散值，到达快速估计信号的色散值的目的。

上文中，结合图 1 , 详细描述了根据本发明实施例的色散估计的方法，下面将结合图 2 , 详细描述根据本发明实施例的色散估计的装置。

图 2示出了根据本发明实施例的色散估计的装置 200的示意性框图。如图 2所示，该装置可以包括：

功率自相关处理单元 210 , 用于计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列，以及向脉沖位置搜索单元 220 传输指示该自相关函数序列的信息；

脉沖位置搜索单元 220 , 用于从该功率自相关处理单元 210获取指示该自相关函数序列的信息，并根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的位置相对应，以及向色散值确定单元 230 传输指示该脉沖位置的信息；

色散值确定单元 230 , 用于从该脉沖位置搜索单元 220指示该脉沖位置的信息，并根据该脉沖位置，确定该目标光信号的色散值。在本发明实施例中，通过试验（例如，获取多个光信号的功率自相关函数，并获取该多个光信号的色散值）发现，在光信号（例如，单偏振的 PSK 信号、 PDM的 PSK信号或 QAM信号存在色散时，该光信号功率的自相关函数会出现一个脉沖峰值，并且，不同的脉沖峰值位置与不同的色散值之间存在映射对应关系，因此，在本发明实施例中，可以根据光信号的信号功率的自相关函数中脉沖峰值的位置，确定该光信号的色散值。

因此，功率自相关处理单元 210可以计算目标光信号的信号功率（光功率）的自相关函数。在本发明实施例中，目标光信号为需要估算色散的光信号。在本发明实施例中，可以对该目标光信号进行以下处理，以单偏振系统为例，在相干接收机中，目标光信号和本振光在 Hybrid 中混合，通过平衡接收机后转化为两路电流（同向电流和正交电流；)，经过两路 ADC采样后，转换为数字信号，得到 I信号和 Q信号，从而，可以根据该 I信号和 Q信号来确定目标光信号的信号功率的自相关函数。以下，作为示例而非限定，以该通过该 I信号和 Q信号确定目标光信号的信号功率的自相关函数为例进行说明。

可选地，在本发明实施例中，在该目标光信号为非归零码格式信号的情况下，该装置还包括：

第一干扰处理单元，用于对该目标光信号进行第一处理，以消除所述目标光信号的零频干扰。

具体地说，零频干扰是指，在该目标光信号为非归零（ NRZ , No Return

可选地，在本发明实施例中，该装置还包括：

第二干扰处理单元，用于对该目标光信号进行第二处理，以消除该目标光信号的直流分量干扰。

具体地说，作为示例而非限定，第二干扰处理单元可以对 I信号和 Q信号（具体地说，是与该 I信号相对应的 I信号序列以及与该 Q信号相对应的 Q信号序列）进行求平方和处理，然后用经过该求平方和处理后的信号减去信号均值，其中，作为示例而非限定，该信号均值可以通过以下方法确定，即，可以获取与该 I信号相对应的 I信号序列，以及与该 Q信号相对应的 Q 信号序列，其后，对 I信号序列和 Q信号序列求和，然后除以 I信号序列和 Q信号序列的长度，从而确定该信号均值。通过以上处理，能够实现对该目标光信号的隔直处理，从而能够去掉直流分量产生的无用的自相关函数部分。傅里叶变换，再做幅度平方，最后进行逆傅里叶变换，得到该目标光信号的信号功率的自相关函数，进而可以根据该自相关函数，生成以该自相关函数的各函数值为元素的自相关函数序列（以下，为了便于说明，记作

脉沖位置搜索单元 220可以从该 W中搜索脉沖峰值的位置，即，确定该脉沖峰值在 w中的序号。

可选地，在本发明实施例中，脉沖位置搜索单元 220可以搜索 W中数值最大的元素（记作元素 a ), 并将该元素 a的位置确定为脉沖峰值的位置，并可以记录该元素 a的序号，以在后述色散值确定单元 230确定色散值时，根据该脉沖峰值在 W中的位置，或者说元素 a在 W中的序号（记作 nj, 确定色散值 (随后进行详细说明 )。

可选地，在本发明实施例中，该脉沖位置搜索单元 220具体用于根据该自相关函数序列，确定第一自相关函数序列，该第一自相关函数序列包括与该脉沖峰值相对应的元素；

用于根据该第一自相关函数序列，确定该脉沖位置。

具体地说，在本发明实施例中，脉沖位置搜索单元 220可以从该 W中选取一个区间作为搜索范围（即，第一自相关函数序列，以下，为了便于说明，记作 ?₂W ),并且，在本发明实施例中，脉沖位置搜索单元 220可以将 W 作为 ] , 也可以将 [A^ WJNJ (即, W中从序号为的元素到序号为 N₂的元素）作为 W₂W , 其中，表示 W中选取的开始位置（序号）， N₂表示 W中选取的结束位置（序号）。

可选地，在本发明实施例中，该脉沖位置搜索单元 220具体用于根据用于获取该目标光信号的相干接收机的电带宽和用于获取该目标光信号的高速模数转换器的采样率，从该自相关函数序列中确定选取该第一自相关函数序列的开始位置。

具体地说，脉沖位置搜索单元 220可以根据式 1确定。

可选地，在本发明实施例中，该脉沖位置搜索单元 220具体用于根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；

用于根据该第一滑动平均处理后的第一自相关函数序列，确定该脉沖位置。

具体地说，脉沖位置搜索单元 220在确定 ?₂W之后，可以对 ?₂W进行滑动窗口平均处理（第一滑动平均处理），并且，作为示例而非限定，该滑动窗口平均处理的窗口宽度（记作 widthj可以为例如， 3。在本发明实施例中，滑动窗口平均处理的方法与现有技术相同，这里，为了避免赘述，省略其说明。在对 ?₂W进行滑动窗口平均处理后（以下，为了便于理解，将处理后生成的序列记作 R₃W ), 可以搜索 R₃ [«]中数值最大的元素（记作元素 b ), 并将该元素 b的位置确定为脉沖峰值的位置，并可以记录该元素 b在 R₃ [«]中的序号（记作 n₂ ), 从而，可以根据式 2确定 _ηι。

可选地，在本发明实施例中，在确定 ?₂W之后，脉沖位置搜索单元 220 也可以直接搜索 ?₂W中数值最大的元素（记作元素 c ), 并将该元素 c的位置确定为脉沖峰值的位置，并可以记录该元素 c在 ?₂W中的序号（记作 n₃ ), 从而，可以根据式 3确定 1^。

从而，在后述确定色散值时，色散值确定单元 230可以根据该脉沖峰值在 W中的位置，或者说 _ηι的值，确定色散值（随后进行详细说明）。

可选地，在本发明实施例中，该脉沖位置搜索单元 220具体用于根据第一序号和第一预设值，确定第二窗口，该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；

用于根据该第一自相关函数序列和该第二窗口，确定第二自相关函数序歹' J ;

用于根据该第二窗口，对该第二自相关函数序列进行第二滑动平均处理；用于根据该第二滑动平均处理后的第二自相关函数序列，确定该脉沖位置。

在本发明实施例中，脉沖位置搜索单元 220例如，可以将自相关函数中与产生的色散相对应的脉沖峰值的 FWHM的一半，即，将脉沖峰值（最大值）的二分之一作为门限（后述第一门限），寻找脉沖宽度，再将脉沖宽度的例如，二分之一作为第二窗口，具体地说，在确定 n₂后，可以从?₂W中，在以 n₂的中心的预设范围（与第一预设值相对应，例如，可以为 5个单位，这里 "单位" 可以是单位序号）内，搜索数值最大的元素（记作元素 d), 并才艮据该元素 d的值，以确定第一门限，这里，作为示例而非限定，将该元素 d的值的二分之一作为该第一门限，其后，以元素 d为中心，向左侧搜索第一个数值小于该第一门限的元素（记作元素 e ), 并记录该元素 e的序号（记作；向右侧搜索第一个数值小于该第一门限的元素（记作元素 f), 并记录该元素 f的序号（记作 t₂), 可以根据式 4, 确定滑动窗口宽度（第二窗口的一例，记作 width₂)。

在确定了 width₂ (第二窗口）后，脉沖位置搜索单元 220可以从该 W₂W 中选取一个区间作为搜索范围（即，第二自相关函数序列，以下，为了便于说明，记作 ?₄W ), 在本发明实施例中，可以将 ?₂W作为 ?₄W, 也可以将 {R₂[N₃] :R₂[N₄}}(即， ?₂W中从序号为 N₃的元素到序号为 N₄的元素 M乍为 ?₄[/i], 其中， N₃表示 ?₂W中选取的开始位置（序号）， N₄表示 ?₂W中选取的结束位置（序号）。具体地说，脉沖位置搜索单元 220可以根据式 5确定 N₃, 并可以根据式 6确定 N₄。

在确定了 ?₄W后，脉沖位置搜索单元 220根据 width₂, 对该 ?₄W进行滑动平均处理（第二滑动平均处理），以下，为了便于理解，将处理后生成的序列记作 [«], 可以搜索 R₅[«]中数值最大的元素（记作元素 g), 并将该元素 g的位置确定为脉沖峰值的位置，并可以记录该元素 g在 [«]中的序号（记作 n₄ ) , 由于自相关函数第一位对应 τ=0 , R W对应延迟为 0 - 1) · Samplinglnterval (修正一位 ) , 其中, Samplinglnterval表示采样间 Ρι¾ , 从而，可以 ^据式 7确定。

从而，在后述确定色散值时，色散值确定单元 230可以根据该脉沖峰值在中的位置，或者说 _ηι的值，确定色散值（随后进行详细说明）。

可选地，在本发明实施例中，该脉沖位置搜索单元 220具体用于根据第一序号和第一预设值，从该第一自相关函数序列中确定第一范围，该第一范围以该第一序号为中心，该第一范围的大小与该第一预设值相对应；

用于根据该第一范围中数值最大的第一元素和第二预设值，从该第一自相关函数序列中确定第二范围，该第二范围以该第一元素的序号为中心，该第二范围的大小与该第二预设值相对应；

用于根据该第一元素和该第二范围中数值第二大的第二元素，确定该第二窗口。

具体地说，如果 PMD较大，则可能会出现一个自相关函数中的脉沖（峰值）分裂为两个脉沖（峰值）的情况，此情况下，如果仅以脉沖峰值（最大值）的位置为中心寻找脉沖宽度，则寻找到的脉沖宽度值会产生偏差，因此，在本发明实施例中，脉沖位置搜索单元 220在确定 n₂后，还可以从 W₂W中，在以 n₂为中心的预设范围（即，第一范围，与第一预设值相对应，例如，可以为 5个单位）内，搜索数值最大的元素（即，元素 d )。

其后，在以元素 d在中的序号为中心的预设范围（即，第二范围，与第二预设值相对应，例如，可以为 10个单位）内，搜索数值第二大的元素（记作元素 h )。

其后，以元素 h为中心，向左侧搜索第一个数值小于该第一门限的元素

(记作元素 i, 并记录该元素 i的序号（记作 t₃ ); 以元素 d为中心，向右侧搜索第一个数值小于该第一门限的元素（记作元素 j ), 并记录该元素 j的序号（记作 t₄ ), 可以 ^据式 8 , 确定一个中间参量（记作 width₃ )。

脉沖位置搜索单元 220在确定了 width₃后，可以从该 ?₂W中选取一个区间作为搜索范围（即，第二自相关函数序列，以下，为了便于说明，记作 R₆ W ), 在本发明实施例中，可以将 w₂w作为 R₆[«] ,也可以将 {R₂[N₅ ] : R₂[N₆]} (即， w₂w 中从序号为 N₅的元素到序号为 N₆的元素）作为 R₆[_w], 其中， N₅表示中选取的开始位置（序号）， N₆表示 ?₂[/1]中选取的结束位置（序号）。具体地说，脉沖位置搜索单元 220可以根据式 9确定 N₅, 并可以根据式 10确定 N₆。

脉沖位置搜索单元 220在确定了 R₆[«]后，可以根据式 11,确定对 R₆[«]进行滑动平均处理时使用的窗口宽度（第二窗口的另一例，记作 width₄)。

其后，脉沖位置搜索单元 220对该 R₆[«]进行滑动平均处理（第二滑动平均处理 ),以下，为了便于理解，将处理后生成的序列记作 R₇[«],可以从该 R₇[«] 中选取一个区间作为搜索范围（以下，为了便于说明，记作 ¾[_w])，在本发明实施例中，可以将 R₇[«]作为 R₈[«], 也可以将 {R₇[N₇]:R₇[N₈]} (即， R₇[«]中从序号为 N₇的元素到序号为 N₈的元素）作为 R₈[«], 其中， N₇表示 R₇[«]中选取的开始位置（序号）， N₈表示 R₇[«]中选取的结束位置（序号）。具体地说，脉沖位置搜索单元 220可以根据式 12确定 N₇, 并可以根据式 13确定 N₈。

脉沖位置搜索单元 220在确定 R₈[«]后，可以搜索 R₈[«]中数值最大的元素 (记作元素 k), 并记录该元素 ]^在 [«]中的序号（记作 n₅)。

在本发明实施例中，可以将该元素 k的位置确定为脉沖峰值的位置，由于自相关函数第一位对应 τ=0 , R [η]对应延迟为 - 1) · Samplinglnterval (修正一位），其中， ^mpZ^tom^表示采样间隔，从而，脉沖位置搜索单元 220 可以根据式 14确定 _ηι。

从而，在后述色散值确定单元 230确定色散值时，根据该脉沖峰值在 [n] 中的位置，或者说 _ηι的值，确定色散值（随后进行详细说明）。

可选地，在本发明实施例中，脉沖位置搜索单元 220在确定 n₅后，还可以利用 R₈[«₅-l]、 R,[n₅], R₈[«₅+l]三点做抛物线拟合，从而，可以根据式 15, 确定该抛物线顶点在 R₈[«]中的位置（序号，记作 n₆)。

其后，脉沖位置搜索单元 220可以根据 n₆和式 16, 确定该抛物线顶点在?₂W中的位置（序号）。

由于自相关函数第一位对应 τ=0 , R [η]对应延迟为 - 1) · Samplinglnterval (修正一位），其中， ^m/^Mg/MferraZ表示采样间隔，从而，可以根据式 17 确定 n

从而，色散值确定单元 230可以根据脉沖位置搜索单元 220确定的自相关函数中的脉沖峰值出现位置（即， nj, 确定目标光信号的色散值。

在本发明实施例中，该色散值确定单元 230具体用于根据离线数据，确定指示脉沖峰值位置与色散值的对应关系的指示信息，根据该脉沖位置和该指示信息，确定该光信号的色散值。

具体地说，作为示例而非限定，例如，色散值确定单元 230可以分别对多个光信号（离线数据的一例 )进行色散值测试等试验，记录各光信号的色散值和它们各自对应的自相关函数中脉沖峰值出现的位置，作为训练数据，从而，可以采用例如期望最大化算法，通过逐步改进公式的参数，使参数和训练数据的似然概率逐渐增大，最后终止于一个极大点。直观地理解，期望最大化算法也可被看作为一个逐次逼近算法：事先并不知道模型的参数，可以随机的选择一套参数或者事先粗略地给定某个初始参数 λθ , 确定出对应于这组参数的最可能的状态，计算每个训练样本的可能结果的概率，在当前的状态下再由样本对参数修正，重新估计参数 λ , 并在新的参数下重新确定模型的状态，这样，通过多次的迭代，循环直至某个收敛条件满足为止，就可以使得模型的参数逐渐逼近真实参数。因此，该步骤可以以离线的方式自动进行。将离线计算得到的上述各种脉沖峰值位置与各自对应的色散值预先存储起来，例如以表的形式存储起来。应理解，以上列发明并不限定于此。

作为示例而非限定，色散值确定单元 230可以根据式 18确定色散值。应理解，以上列举的各公式、参数以及数值仅为本发明的示例性说明，本发明并不限定于此，其他能够用于确定光信号的自相关函数的脉沖峰值出现位置的各公式、参数以及数值均落入本发明的保护范围内。

根据本发明实施例的色散估计的装置可对应于本发明实施例的色散估计的方法的执行主体，并且，该色散估计的装置中的各单元和上述其他操作和 /或功能分别为了实现图 1中的方法 100的相应流程，为了筒洁，在此不再赘述。

根据本发明实施例的色散估计的装置，通过获取目标光信号的光信号功率自相关函数，并确定该光信号功率自相关函数出现脉沖峰值的位置，能够在无需扫描所有可能的色散值的情况下，根据该出现脉沖峰值的位置确定目标光信号的色散值，到达快速估计信号的色散值的目的。

上文中，结合图 1和图 2, 详细描述了根据本发明实施例的色散估计的方法和装置，下面将结合图 3 , 详细描述根据本发明实施例的色散估计的系统。图 3示出了根据本发明实施例的色散估计的装置 300的示意性框图。如图 3所示，该系统可以包括：

总线 310;

与该总线相连的处理器 320;

与该总线相连的存储器 330;

其中，该处理器 320通过该总线，调用该存储器 330中存储的程序，以用于计算目标光信号的信号功率的自相关函数，并获取该自相关函数的自相关函数序列；

根据该自相关函数序列，确定脉沖位置，该脉沖位置与该自相关函数的脉沖峰值的出现位置相对应；

根据该脉沖位置，确定该目标光信号的色散值。

在本发明实施例中，通过试验（例如，获取多个光信号的功率自相关函数，并获取该多个光信号的色散值）发现，在光信号（例如，单偏振的 PSK 信号、 PDM的 PSK信号或 QAM信号存在色散时，该光信号功率的自相关函数会出现一个脉沖峰值，并且，不同的脉沖峰值位置与不同的色散值之间存在映射对应关系，因此，在本发明实施例中，可以根据光信号的信号功率的自相关函数中脉沖峰值的位置，确定该光信号的色散值。

因此，处理器 320可以计算目标光信号的信号功率（光功率）的自相关函数。在本发明实施例中，目标光信号为需要估算色散的光信号。在本发明实施例中，可以对该目标光信号进行以下处理，以单偏振系统为例，在相干接收机中，目标光信号和本振光在 Hybrid 中混合，通过平衡接收机后转化为两路电流（同向电流和正交电流），经过两路 ADC采样后，转换为数字信号，得到 I信号和 Q信号，从而，可以根据该 I信号和 Q信号来确定目标光信号的信号功率的自相关函数。以下，作为示例而非限定，以该通过该 I信号和 Q信号确定目标光信号的信号功率的自相关函数为例进行说明。。

可选地，在本发明实施例中，在该目标光信号为非归零码格式信号的情况下，该处理器 320通过该总线 310, 调用该存储器 330中存储的程序，以消除所述目标光信号的零频干扰。

具体地说，零频干扰是指，在该目标光信号为非归零（NRZ, No Return Zero )码格式的情况下，如果直接获取该格式的信号的自相关函数，则获取的自相关函数的零频位置会出现一个较大的峰值，从而对后述搜索脉沖峰值造成影响，因此，作为示例而非限定，在本发明实施例中，可以对 I信号和 Q信号进行延迟处理和相减处理，以实现对该目标光信号的高通滤波，其中，该延迟处理使用的迟延大小可以为例如，一个符号周期。从而，能够避免获取的自相关函数在零频位置出现较大的峰值。

可选地，在本发明实施例中，该处理器 320通过该总线 310, 调用该存储器 330中存储的程序，以用于对该目标光信号进行第二处理，以消除该目标光信号的直流分量干扰。

具体地说，作为示例而非限定，处理器 320可以对 I信号和 Q信号（具体地说，是与该 I信号相对应的 I信号序列以及与该 Q信号相对应的 Q信号序列）进行求平方和处理，然后用经过该求平方和处理后的信号减去信号均值，其中，作为示例而非限定，该信号均值可以通过以下方法确定，即，可以获取与该 I信号相对应的 I信号序列，以及与该 Q信号相对应的 Q信号序歹 ij , 其后，对 I信号序列和 Q信号序列求和，然后除以 I信号序列和 Q信号序列的长度，从而确定该信号均值。通过以上处理，能够实现对该目标光信号的隔直处理，从而能够去掉直流分量产生的无用的自相关函数部分。

其后，处理器 320可以对经过上述处理的目标光信号进行傅里叶变换，再做幅度平方，最后进行逆傅里叶变换，得到该目标光信号的信号功率的自相关函数，进而可以根据该自相关函数，生成以该自相关函数的各函数值为元素的自相关函数序列（以下，为了便于说明，记作 w )

处理器 320可以从该 W中搜索脉沖峰值的位置，即，确定该脉沖峰值在 w中的序号。

可选地，在本发明实施例中处理器 320可以搜索 W中数值最大的元素 (记作元素 a ), 并将该元素 a的位置确定为脉沖峰值的位置，并可以记录该元素 a的序号，以在后述色散值确定单元 230确定色散值时，根据该脉沖峰值在 W中的位置，或者说元素 a在 W中的序号（记作 nj, 确定色散值 (随后进行详细说明）。

可选地，在本发明实施例中，该处理器 320通过该总线 310, 调用该存储器 330中存储的程序，以具体用于根据该自相关函数序列，确定第一自相关函数序列，该第一自相关函数序列包括与该脉沖峰值相对应的元素；

根据该第一自相关函数序列，确定该脉沖位置。

具体地说，在本发明实施例中，处理器 320可以从该 W中选取一个区间作为搜索范围（即，第一自相关函数序列，以下，为了便于说明 ,记作 ?₂ W ), 并且，在本发明实施例中，处理器 320 可以将 W作为也可以将 {^[NJi^fNJIi即， W中从序号为的元素到序号为 N₂的元素 M乍为 ?₂[/i] , 其中，表示 W中选取的开始位置（序号）， N₂表示 W中选取的结束位置（序号）。

可选地，在本发明实施例中，处理器 320具体用于根据用于获取该目标光信号的相干接收机的电带宽和用于获取该目标光信号的高速模数转换器的采样率，从该自相关函数序列中确定选取该第一自相关函数序列的开始位置。

具体地说，处理器 320可以根据式 1确定 Ni。

可选地，在本发明实施例中，该处理器 320通过该总线 310, 调用该存储器 330中存储的程序，以具体用于根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；

具体地说，处理器 320在确定 ?₂W之后，可以对 ?₂W进行滑动窗口平均处理（第一滑动平均处理），并且，作为示例而非限定，该滑动窗口平均处理的窗口宽度（记作 widthj可以为例如， 3。在本发明实施例中，滑动窗口平均处理的方法与现有技术相同，这里，为了避免赘述，省略其说明。在对 ?₂[/ι]进行滑动窗口平均处理后（以下，为了便于理解，将处理后生成的序列记作 R₃W ), 可以搜索 R₃ [«]中数值最大的元素（记作元素 b ), 并将该元素 b的位置确定为脉沖峰值的位置，并可以记录该元素 b在 [«]中的序号（记作 n₂ ), 从而，可以根据式 2确定 _ηι。

可选地，在本发明实施例中，在确定 ?₂W之后，处理器 320也可以直接搜索 ?₂W中数值最大的元素（记作元素 c ), 并将该元素 c的位置确定为脉沖峰值的位置，并可以记录该元素 c在 ?₂W中的序号（记作 n₃ ), 从而，可以根据式 3确定 _ηι。从而，在后述确定色散值时，处理器 320可以根据该脉沖峰值在 W中的位置，或者说的值，确定色散值（随后进行详细说明）。

可选地，在本发明实施例中，该处理器 320通过该总线 310, 调用该存储器 330中存储的程序，以具体用于根据第一序号和第一预设值，确定第二窗口，该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；

用于根据该第二窗口，对该第二自相关函数序列进行第二滑动平均处理；

用于根据该第二滑动平均处理后的第二自相关函数序列，确定该脉沖位置。

在本发明实施例中，处理器 320例如，可以将自相关函数中与产生的色散相对应的脉沖峰值的 FWHM的一半，即，将脉沖峰值（最大值）的二分之一作为门限（后述第一门限），寻找脉沖宽度，再将脉沖宽度的例如，二分之一作为第二窗口，具体地说，在确定 n₂后，可以从 ?₂W中，在以 n₂为中心的预设范围（即，第一范围，例如，可以为 5个单位，这里 "单位" 可以是单位序号）内，搜索数值最大的元素（记作元素 d ), 并根据该元素 d 的值，以确定第一门限，这里，作为示例而非限定，将该元素 d的值的二分之一作为该第一门限，其后，以元素 d为中心，向左侧搜索第一个数值小于该第一门限的元素（记作元素 e ), 并记录该元素 e的序号（记作；向右侧搜索第一个数值小于该第一门限的元素（记作元素 f), 并记录该元素 f的序号（记作 t₂ ), 可以根据式 4, 确定滑动窗口宽度（第二窗口的一例，记作 width₂ )。

处理器 320在确定了 width₂ (第二窗口）后，可以从该 ?₂W中选取一个区间作为搜索范围（即，第二自相关函数序列，以下，为了便于说明，记作 R₄[n] ), 在本发明实施例中，可以将 W₂W作为 W₄W , 也可以将 {R₂ [N₃ ] : R₂ [N₄]} (即， ?₂W中从序号为 N₃的元素到序号为 N₄的元素）作为 ?₄W , 其中， N₃ 表示 ?₂W中选取的开始位置（序号 ), N₄表示 ?₂W中选取的结束位置（序号）。具体地说，处理器 320可以根据式 5确定 N₃, 并可以根据式 6确定 N₄。

处理器 320在确定了 ?₄W后，可以根据 width₂, 对该 ?₄W进行滑动平均处理（第二滑动平均处理），以下，为了便于理解，将处理后生成的序列记作 R₅w , 可以搜索 R₅w中数值最大的元素（记作元素 g ), 并将该元素 g的位置确定为脉沖峰值的位置，并可以记录该元素 g在 R₅ [«]中的序号（记作 n₄ ), 由于自相关函数第一位对应 τ=0 , R [η]对应延迟为 - 1) · Samplinglnterval (修正一位），其中， ^m/^/^/wferraZ表示采样间隔，从而，处理器 320可以根据式 7确定。

从而，处理器 320在后述确定色散值时，可以根据该脉沖峰值在 W中的位置，或者说的值，确定色散值（随后进行详细说明）。

通过对该 (第一自相关函数序列）进行滑动窗口平均处理，能够减少自相关函数中噪声的干扰，使确定的脉沖峰值的位置（脉沖位置）更加准确。

可选地，在本发明实施例中，该处理器 320通过该总线 310, 调用该存储器 330中存储的程序，以具体用于根据第一序号和第一预设值，从该第一自相关函数序列中确定第一范围，该第一范围以该第一序号为中心，该第一范围的大小与该第一预设值相对应；

具体地说，如果 PMD较大，则可能会出现一个自相关函数中的脉沖（峰值）分裂为两个脉沖（峰值）的情况，此情况下，如果仅以脉沖峰值（最大值）的位置为中心寻找脉沖宽度，则寻找到的脉沖宽度值会产生偏差，因此，在本发明实施例中，处理器 320在确定 n₂后，还可以从 ?₂W中，在以 n₂为中心的预设范围（即，第一范围，与第一预设值相对应，例如，可以为 5个单位）内，搜索数值最大的元素（即，元素 d )。

其后，在以元素 d在中的序号为中心的的预设范围（即，第二范围，与第二预设值相对应，例如，可以为 10个单位，这里为单位序号），搜索数值第二大的元素（记作元素 h )。

其后，确定元素 d与元素 h在 ?₂W中的位置关系，这里，为了便于理解，以元素 d在元素 h的右侧为例（即，元素 d在 ?₂ W中的序号大于元素 h在 ?₂ W 中的序号）进行说明。

其后，以元素 h为中心，向左侧搜索第一个数值小于该第一门限的元素 (记作元素 i, 并记录该元素 i的序号（记作 t₃ ); 以元素 d为中心，向右侧搜索第一个数值小于该第一门限的元素（记作元素 j), 并记录该元素 j的序号（记作 t₄), 可以 ^据式 8, 确定一个中间参量（记作 width₃)。

处理器 320在确定了 width₃后，可以从该中选取一个区间作为搜索范围（即，第二自相关函数序列，以下，为了便于说明，记作在本发明实施例中，可以将 w₂W作为 R₆[«], 也可以将 {R₂[N₅]:R₂[N₆]} (即， w₂W中从序号为 N₅的元素到序号为 N₆的元素）作为 R₆W,其中， N₅表示?₂W中选取的开始位置（序号）， N₆表示 ?₂[/ι]中选取的结束位置（序号）。具体地说，处理器 320可以根据式 9确定 N₅, 并可以根据式 10确定 N₆。

其后，处理器 320对该 R₆W进行滑动平均处理（第二滑动平均处理），以下，为了便于理解，将处理后生成的序列记作 R₇[«], 可以从该 R₇[«]中选取一个区间作为搜索范围（以下，为了便于说明，记作 R₈[«])，在本发明实施例中，可以将 R₇[«]作为 R₈[«], 也可以将 {R₇[N₇]:R₇[N₈]} (即， R₇[«]中从序号为 N₇的元素到序号为 N₈的元素）作为 R₈[«], 其中， N₇表示 R₇[«]中选取的开始位置（序号）， N₈表示 R₇[«]中选取的结束位置（序号）。具体地说，脉沖位置搜索单元 220可以根据式 12确定 N₇, 并可以根据式 13确定 N₈。

处理器 320在确定 R₈[«]后，可以搜索 R₈[«]中数值最大的元素（记作元素 k), 并记录该元素 ]^在 [«]中的序号（记作 n₅)。

在本发明实施例中，处理器 320可以将该元素 k的位置确定为脉沖峰值的位置，由于自相关函数第一位对应 τ=0 , R W对应延迟为 0 - 1) · Samplinglnterval (修正一位 ) , 其中, Samplinglnterval表示采样间 Ρι¾ , 从而，处理器 320可以 ^据式 14确定 _ηι。

从而，在后述处理器 320确定色散值时，根据该脉沖峰值在 W中的位置，或者说 _ηι的值，确定色散值（随后进行详细说明）。

可选地，在本发明实施例中，处理器 320 在确定 n₅后，还可以利用 R_& [n₅ - \] , R_& [n₅ ] , R₈[«₅ + l]三点做抛物线拟合，从而，可以根据式 15 , 确定该抛物线顶点在中的位置（序号，记作 n₆ )。

其后，处理器 320可以根据 n₆和式 16, 确定该抛物线顶点在 ?₂W中的位置（序号）。

由于自相关函数第一位对应 τ=0 , R [η]对应延迟为 - 1) · Samplinglnterval

(修正一位），其中， ^m/^Mg/MferraZ表示采样间隔，从而，可以根据式 17 确定 n

从而，处理器 320可以根据所确定的自相关函数中的脉沖峰值出现位置 (即， nj, 确定目标光信号的色散值。

在本发明实施例中，处理器 320根据离线数据，确定指示脉沖峰值位置与色散值的对应关系的指示信息，

根据该脉沖位置和该指示信息，确定该光信号的色散值。

具体地说，作为示例而非限定，例如，处理器 320可以分别对多个光信号（离线数据的一例 )进行色散值测试等试验，记录各光信号的色散值和它们各自对应的自相关函数中脉沖峰值出现的位置，作为训练数据，从而，可以采用例如期望最大化算法，通过逐步改进公式的参数，使参数和训练数据的似然概率逐渐增大，最后终止于一个极大点。直观地理解，期望最大化算法也可被看作为一个逐次逼近算法：事先并不知道模型的参数，可以随机的选择一套参数或者事先粗略地给定某个初始参数 λθ , 确定出对应于这组参数的最可能的状态，计算每个训练样本的可能结果的概率，在当前的状态下再由样本对参数修正，重新估计参数 λ , 并在新的参数下重新确定模型的状态，这样，通过多次的迭代，循环直至某个收敛条件满足为止，就可以使得模型的参数逐渐逼近真实参数。因此，该步骤可以以离线的方式自动进行。将离线计算得到的上述各种脉沖峰值位置与各自对应的色散值预先存储起来，例如以表的形式存储起来。应理解，以上列举的各不限定于此。

作为示例而非限定，处理器 320可以根据式 18确定色散值。

应理解，以上列举的各公式、参数以及数值仅为本发明的示例性说明，本发明并不限定于此，其他能够用于确定光信号的自相关函数的脉沖峰值出现位置的各公式、参数以及数值均落入本发明的保护范围内。根据本发明实施例的色散估计的系统可对应于本发明实施例的色散估计的方法的执行主体，并且，该色散估计的装置中的各单元和上述其他操作和 /或功能分别为了实现图 1中的方法 100的相应流程，为了筒洁，在此不再赘述。

根据本发明实施例的色散估计的系统，通过获取目标光信号的光信号功率自相关函数，并确定该光信号功率自相关函数出现脉沖峰值的位置，能够在无需扫描所有可能的色散值的情况下，根据该出现脉沖峰值的位置确定目标光信号的色散值，到达快速估计信号的色散值的目的。应理解，在本发明的各种实施例中，上述各过程的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

所属领域的技术人员可以清楚地了解到，为描述的方便和筒洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。

所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括： U盘、移动硬盘、只读存储器（ ROM, Read-Only Memory )、随机存取存储器（RAM, Random Access Memory ), 磁碟或者光盘等各种可以存储程序代码的介质。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims

权利要求

1、一种色散估计的方法，其特征在于，所述方法包括：

计算目标光信号的信号功率的自相关函数，并获取所述自相关函数的自相关函数序列；

根据所述自相关函数序列，确定脉沖位置，所述脉沖位置与所述自相关函数的脉沖峰值的位置相对应；

根据所述脉沖位置，确定所述目标光信号的色散值。

2、根据权利要求 1 所述的方法，其特征在于，在所述目标光信号为非归零码格式信号的情况下，在根据目标光信号的信号功率的自相关函数，确定自相关函数序列前，所述方法还包括：

对所述目标光信号进行第一处理，以消除所述目标光信号的零频干扰。

3、根据权利要求 1或 2所述的方法，其特征在于，在根据目标光信号的信号功率的自相关函数，确定自相关函数序列前，所述方法还包括：对所述目标光信号进行第二处理，以消除所述目标光信号的直流分量干扰。

4、根据权利要求 1至 3中任一项所述的方法，其特征在于，所述根据所述自相关函数序列，确定脉沖位置，包括：

根据所述自相关函数序列，确定第一自相关函数序列，所述第一自相关函数序列包括所述脉沖峰值；

根据所述第一自相关函数序列，确定所述脉沖位置。

5、根据权利要求 4所述的方法，其特征在于，所述根据所述第一自相关函数序列，确定所述脉沖位置，包括：

根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；根据所述第一滑动平均处理后的第一自相关函数序列，确定所述脉沖位置。

6、根据权利要求 5所述的方法，其特征在于，所述根据所述第一滑动平均处理后的第一自相关函数序列，确定所述脉沖位置，包括：

根据第一序号和第一预设值，确定第二窗口，所述第一序号是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；

根据所述第一自相关函数序列和所述第二窗口，确定第二自相关函数序歹' J ; 根据所述第二窗口，对所述第二自相关函数序列进行第二滑动平均处理；

根据所述第二滑动平均处理后的第二自相关函数序列，确定所述脉沖位置。

7、根据权利要求 6所述的方法，其特征在于，所述根据第一序号和第一预设值，确定第二窗口，包括：

根据第一序号和第一预设值，从所述第一自相关函数序列中确定第一范围，所述第一范围以所述第一序号为中心，所述第一范围的大小与所述第一预设值相对应；

根据所述第一范围中数值最大的第一元素和第二预设值，从所述第一自相关函数序列中确定第二范围，所述第二范围以所述第一元素的序号为中心，所述第二范围的大小与所述第二预设值相对应；

根据所述第一元素和所述第二范围中数值第二大的第二元素，确定所述第二窗口。

8、根据权利要求 1至 7中任一项所述的方法，其特征在于，所述根据所述脉沖位置，确定所述目标光信号的色散值，包括：

根据以下公式确定所述目标光信号的色散值，

η

χ =

2π · SamplingRate · SymbolRate

其中， CD表示色散值， n表示所述脉沖位置在所述自相关函数序列中的序号， Γ。表示所述目标光信号的初始发射机光脉沖宽度， c表示真空光速， λ 表示所述目标光信号的波长， ^m/^^R 表示模数转换器的采样率， ¾™bo/ ¾fe表示所述目标光信号的符号速率。

9、一种色散估计的装置，其特征在于，所述装置包括：

功率自相关处理单元，用于计算目标光信号的信号功率的自相关函数，并获取所述自相关函数的自相关函数序列，以及向脉沖位置搜索单元传输指示所述自相关函数序列的信息；

脉沖位置搜索单元，用于从所述功率自相关处理单元获取指示所述自相关函数序列的信息，并根据所述自相关函数序列，确定脉沖位置，所述脉沖位置与所述自相关函数的脉沖峰值的位置相对应，以及向色散值确定单元传输指示所述脉沖位置的信息；

色散值确定单元，用于从所述脉沖位置搜索单元获取指示所述脉沖位置的信息，并根据所述脉沖位置，确定所述目标光信号的色散值。

10、根据权利要求 9所述的装置，其特征在于，在所述目标光信号为非归零码格式信号的情况下，所述装置还包括：

第一干扰处理单元，用于对所述目标光信号进行第一处理，以消除所述目标光信号的零频干扰。

11、根据权利要求 9或 10所述的装置，其特征在于，所述装置还包括：第二干扰处理单元，用于对所述目标光信号进行第二处理，以消除所述目标光信号的直流分量干扰。

12、根据权利要求 9至 11 中任一项所述的装置，其特征在于，所述脉沖位置搜索单元具体用于根据所述自相关函数序列，确定第一自相关函数序列，所述第一自相关函数序列包括所述脉沖峰值；

用于根据所述第一自相关函数序列，确定所述脉沖位置。

13、根据权利要求 12所述的装置，其特征在于，所述脉沖位置搜索单元具体用于根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；

用于根据所述第一滑动平均处理后的第一自相关函数序列，确定所述脉沖位置。

14、根据权利要求 13所述的装置，其特征在于，所述脉沖位置搜索单元具体用于根据第一序号和第一预设值，确定第二窗口，所述第一序号是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；用于根据所述第一自相关函数序列和所述第二窗口，确定第二自相关函数序列；

用于根据所述第二窗口，对所述第二自相关函数序列进行第二滑动平均处理；

用于根据所述第二滑动平均处理后的第二自相关函数序列，确定所述脉沖位置。

15、根据权利要求 14所述的装置，其特征在于，所述脉沖位置搜索单元具体用于根据第一序号和第一预设值，从所述第一自相关函数序列中确定第一范围，所述第一范围以所述第一序号为中心，所述第一范围的大小与所述第一预设值相对应；

用于根据所述第一范围中数值最大的第一元素和第二预设值，从所述第一自相关函数序列中确定第二范围，所述第二范围以所述第一元素的序号为中心，所述第二范围的大小与所述第二预设值相对应；

用于根据所述第一元素和所述第二范围中数值第二大的第二元素，确定所述第二窗口。

16、根据权利要求 9至 15中任一项所述的装置，其特征在于，所述色

种色散估计的系统，其特征在于，所述系统包括:

总线；

与所述总线相连的处理器；

与所述总线相连的存储器；

其中，所述处理器通过所述总线，调用所述存储器中存储的程序，以用于计算目标光信号的信号功率的自相关函数，并获取所述自相关函数的自相关函数序列；

根据所述自相关函数序列，确定脉沖位置，所述脉沖位置与所述自相关函数的脉沖峰值的出现位置相对应；

根据所述脉沖位置，确定所述目标光信号的色散值。

18、根据权利要求 17所述的系统，其特征在于，在所述目标光信号为非归零码格式信号的情况下，所述处理器通过所述总线，调用所述存储器中存储的程序，以用于对所述目标光信号进行第一处理，以消除所述目标光信号的零频干扰。

19、根据权利要求 17或 18所述的系统，其特征在于，所述处理器通过所述总线，调用所述存储器中存储的程序，以用于对所述目标光信号进行第二处理，以消除所述目标光信号的直流分量干扰。

20、根据权利要求 17至 19中任一项所述的系统，其特征在于，所述处理器通过所述总线，调用所述存储器中存储的程序，以具体用于根据所述自相关函数序列，确定第一自相关函数序列，所述第一自相关函数序列包括所述脉沖峰值；

根据所述第一自相关函数序列，确定所述脉沖位置。

21、根据权利要求 20所述的系统，其特征在于，所述处理器通过所述总线，调用所述存储器中存储的程序，以具体用于根据预设的第一窗口，对第一自相关函数序列进行第一滑动平均处理；

22、根据权利要求 21所述的系统，其特征在于，所述处理器通过所述总线，调用所述存储器中存储的程序，以具体用于根据第一序号和第一预设值，确定第二窗口，所述第一序号是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号；

用于根据所述第一自相关函数序列和所述第二窗口，确定第二自相关函数序列；

23、根据权利要求 21所述的系统，其特征在于，所述处理器通过所述总线，调用所述存储器中存储的程序，以具体用于根据第一序号和第一预设值，从所述第一自相关函数序列中确定第一范围，所述第一范围以所述第一序号为中心，所述第一范围的大小与所述第一预设值相对应；

24、根据权利要求 17至 23中任一项所述的系统，其特征在于，所述处理器通过所述总线，调用所述存储器中存储的程序，以具体用于根据以下公式确定所述目标光信号的色散值，

2π · SamplingRate · SymbolRate

其中， CD表示色散值， n表示所述色散出现位置在所述自相关函数序列中的序号， Γ。表示所述目标光信号的初始发射机光脉沖宽度， c表示真空光速，表示所述目标光信号的波长， ^m ^R 表示模数转换器的采样率， ¾™bo/ ¾fe表示所述目标光信号的符号速率。