色散估计的方法、装置和系统.pdf

上传人:小** 文档编号:4088480 上传时间:2018-08-27 格式:PDF 页数:29 大小:2.61MB
返回 下载 相关 举报
摘要
申请专利号:

CN201280001461.X

申请日:

2012.08.03

公开号:

CN103004110A

公开日:

2013.03.27

当前法律状态:

授权

有效性:

有权

法律详情:

授权|||实质审查的生效IPC(主分类):H04B 10/077申请日:20120803|||公开

IPC分类号:

H04B10/077(2013.01)I; H04B10/2513(2013.01)I

主分类号:

H04B10/077

申请人:

华为技术有限公司; 香港理工大学

发明人:

隋琪; 刘伯涛; 吕超; 李良川

地址:

518129 广东省深圳市龙岗区坂田华为总部办公楼

优先权:

专利代理机构:

北京龙双利达知识产权代理有限公司 11329

代理人:

毛威;张亮

PDF下载: PDF下载
内容摘要

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

权利要求书

权利要求书一种色散估计的方法,其特征在于,所述方法包括:
计算目标光信号的信号功率的自相关函数,并获取所述自相关函数的自相关函数序列;
根据所述自相关函数序列,确定脉冲位置,所述脉冲位置与所述自相关函数的脉冲峰值的位置相对应;
根据所述脉冲位置,确定所述目标光信号的色散值。
根据权利要求1所述的方法,其特征在于,在所述目标光信号为非归零码格式信号的情况下,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,所述方法还包括:
对所述目标光信号进行第一处理,以消除所述目标光信号的零频干扰。
根据权利要求1或2所述的方法,其特征在于,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,所述方法还包括:
对所述目标光信号进行第二处理,以消除所述目标光信号的直流分量干扰。
根据权利要求1至3中任一项所述的方法,其特征在于,所述根据所述自相关函数序列,确定脉冲位置,包括:
根据所述自相关函数序列,确定第一自相关函数序列,所述第一自相关函数序列包括所述脉冲峰值;
根据所述第一自相关函数序列,确定所述脉冲位置。
根据权利要求4所述的方法,其特征在于,所述根据所述第一自相关函数序列,确定所述脉冲位置,包括:
根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;
根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置。
根据权利要求5所述的方法,其特征在于,所述根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置,包括:
根据第一序号和第一预设值,确定第二窗口,所述第一序号是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;
根据所述第一自相关函数序列和所述第二窗口,确定第二自相关函数序列;
根据所述第二窗口,对所述第二自相关函数序列进行第二滑动平均处理;
根据所述第二滑动平均处理后的第二自相关函数序列,确定所述脉冲位置。
根据权利要求6所述的方法,其特征在于,所述根据第一序号和第一预设值,确定第二窗口,包括:
根据第一序号和第一预设值,从所述第一自相关函数序列中确定第一范围,所述第一范围以所述第一序号为中心,所述第一范围的大小与所述第一预设值相对应;
根据所述第一范围中数值最大的第一元素和第二预设值,从所述第一自相关函数序列中确定第二范围,所述第二范围以所述第一元素的序号为中心,所述第二范围的大小与所述第二预设值相对应;
根据所述第一元素和所述第二范围中数值第二大的第二元素,确定所述第二窗口。
根据权利要求1至7中任一项所述的方法,其特征在于,所述根据所述脉冲位置,确定所述目标光信号的色散值,包括:
根据以下公式确定所述目标光信号的色散值,
<mrow><MI>CD</MI><MO>=</MO><MROW><MO>(</MO><MI>x</MI><MO>+</MO><MSQRT><MSUP><MI>x</MI><MN>2</MN></MSUP><MO>-</MO><MSUBSUP><MROW><MN>4</MN><MI>T</MI></MROW><MN>0</MN><MN>4</MN></MSUBSUP></MSQRT><MO>)</MO></MROW><MO>&amp;CenterDot;</MO><MI>π</MI><MO>&amp;CenterDot;</MO><MFRAC><MI>c</MI><MSUP><MI>λ</MI><MN>2</MN></MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS><BR><MATHS id=cmaths0002 num="0002"><MATH><![CDATA[<mrow><MI>x</MI><MO>=</MO><MFRAC><MI>n</MI><MROW><MN>2</MN><MI>π</MI><MO>&amp;CenterDot;</MO><MI>SamplingRate</MI><MO>&amp;CenterDot;</MO><MI>SymbolRate</MI></MROW></MFRAC></MROW>]]&gt;</MATH></MATHS><BR>其中,CD表示色散值,n表示所述脉冲位置在所述自相关函数序列中的序号,T0表示所述目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示所述目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光信号的符号速率。<BR>一种色散估计的装置,其特征在于,所述装置包括:<BR>功率自相关处理单元,用于计算目标光信号的信号功率的自相关函数,并获取所述自相关函数的自相关函数序列,以及向脉冲位置搜索单元传输指示所述自相关函数序列的信息;<BR>脉冲位置搜索单元,用于从所述功率自相关处理单元获取指示所述自相关函数序列的信息,并根据所述自相关函数序列,确定脉冲位置,所述脉冲位置与所述自相关函数的脉冲峰值的位置相对应,以及向色散值确定单元传输指示所述脉冲位置的信息;<BR>色散值确定单元,用于从所述脉冲位置搜索单元获取指示所述脉冲位置的信息,并根据所述脉冲位置,确定所述目标光信号的色散值。<BR>根据权利要求9所述的装置,其特征在于,在所述目标光信号为非归零码格式信号的情况下,所述装置还包括:<BR>第一干扰处理单元,用于对所述目标光信号进行第一处理,以消除所述目标光信号的零频干扰。<BR>根据权利要求9或10所述的装置,其特征在于,所述装置还包括:<BR>第二干扰处理单元,用于对所述目标光信号进行第二处理,以消除所述目标光信号的直流分量干扰。<BR>根据权利要求9至11中任一项所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据所述自相关函数序列,确定第一自相关函数序列,所述第一自相关函数序列包括所述脉冲峰值;<BR>用于根据所述第一自相关函数序列,确定所述脉冲位置。<BR>根据权利要求12所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;<BR>用于根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置。<BR>根据权利要求13所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据第一序号和第一预设值,确定第二窗口,所述第一序号是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;<BR>用于根据所述第一自相关函数序列和所述第二窗口,确定第二自相关函数序列;<BR>用于根据所述第二窗口,对所述第二自相关函数序列进行第二滑动平均处理;<BR>用于根据所述第二滑动平均处理后的第二自相关函数序列,确定所述脉冲位置。<BR>根据权利要求14所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据第一序号和第一预设值,从所述第一自相关函数序列中确定第一范围,所述第一范围以所述第一序号为中心,所述第一范围的大小与所述第一预设值相对应;<BR>用于根据所述第一范围中数值最大的第一元素和第二预设值,从所述第一自相关函数序列中确定第二范围,所述第二范围以所述第一元素的序号为中心,所述第二范围的大小与所述第二预设值相对应;<BR>用于根据所述第一元素和所述第二范围中数值第二大的第二元素,确定所述第二窗口。<BR>根据权利要求9至15中任一项所述的装置,其特征在于,所述色散值确定单元具体用于根据以下公式确定所述目标光信号的色散值,<BR><MATHS id=cmaths0003 num="0003"><MATH><![CDATA[<mrow><MI>CD</MI><MO>=</MO><MROW><MO>(</MO><MI>x</MI><MO>+</MO><MSQRT><MSUP><MI>x</MI><MN>2</MN></MSUP><MO>-</MO><MSUBSUP><MROW><MN>4</MN><MI>T</MI></MROW><MN>0</MN><MN>4</MN></MSUBSUP></MSQRT><MO>)</MO></MROW><MO>&amp;CenterDot;</MO><MI>π</MI><MO>&amp;CenterDot;</MO><MFRAC><MI>c</MI><MSUP><MI>λ</MI><MN>2</MN></MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS><BR><MATHS id=cmaths0004 num="0004"><MATH><![CDATA[<mrow><MI>x</MI><MO>=</MO><MFRAC><MI>n</MI><MROW><MN>2</MN><MI>π</MI><MO>&amp;CenterDot;</MO><MI>SamplingRate</MI><MO>&amp;CenterDot;</MO><MI>SymbolRate</MI></MROW></MFRAC></MROW>]]&gt;</MATH></MATHS><BR>其中,CD表示色散值,n表示所述脉冲位置在所述自相关函数序列中的序号,T0表示所述目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示所述目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光信号的符号速率。<BR>一种色散估计的系统,其特征在于,所述系统包括:<BR>总线;<BR>与所述总线相连的处理器;<BR>与所述总线相连的存储器;<BR>其中,所述处理器通过所述总线,调用所述存储器中存储的程序,以用于计算目标光信号的信号功率的自相关函数,并获取所述自相关函数的自相关函数序列;<BR>根据所述自相关函数序列,确定脉冲位置,所述脉冲位置与所述自相关函数的脉冲峰值的出现位置相对应;<BR>根据所述脉冲位置,确定所述目标光信号的色散值。<BR>根据权利要求17所述的系统,其特征在于,在所述目标光信号为非归零码格式信号的情况下,所述处理器通过所述总线,调用所述存储器中存储的程序,以用于对所述目标光信号进行第一处理,以消除所述目标光信号的零频干扰。<BR>根据权利要求17或18所述的系统,其特征在于,所述处理器通过所述总线,调用所述存储器中存储的程序,以用于对所述目标光信号进行第二处理,以消除所述目标光信号的直流分量干扰。<BR>根据权利要求17至19中任一项所述的系统,其特征在于,所述处理器通过所述总线,调用所述存储器中存储的程序,以具体用于根据所述自相关函数序列,确定第一自相关函数序列,所述第一自相关函数序列包括所述脉冲峰值;<BR>根据所述第一自相关函数序列,确定所述脉冲位置。<BR>根据权利要求20所述的系统,其特征在于,所述处理器通过所述总线,调用所述存储器中存储的程序,以具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;<BR>用于根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置。<BR>根据权利要求21所述的系统,其特征在于,所述处理器通过所述总线,调用所述存储器中存储的程序,以具体用于根据第一序号和第一预设值,确定第二窗口,所述第一序号是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;<BR>用于根据所述第一自相关函数序列和所述第二窗口,确定第二自相关函数序列;<BR>用于根据所述第二窗口,对所述第二自相关函数序列进行第二滑动平均处理;<BR>用于根据所述第二滑动平均处理后的第二自相关函数序列,确定所述脉冲位置。<BR>根据权利要求21所述的系统,其特征在于,所述处理器通过所述总线,调用所述存储器中存储的程序,以具体用于根据第一序号和第一预设值,从所述第一自相关函数序列中确定第一范围,所述第一范围以所述第一序号为中心,所述第一范围的大小与所述第一预设值相对应;<BR>用于根据所述第一范围中数值最大的第一元素和第二预设值,从所述第一自相关函数序列中确定第二范围,所述第二范围以所述第一元素的序号为中心,所述第二范围的大小与所述第二预设值相对应;<BR>用于根据所述第一元素和所述第二范围中数值第二大的第二元素,确定所述第二窗口。<BR>根据权利要求17至23中任一项所述的系统,其特征在于,所述处理器通过所述总线,调用所述存储器中存储的程序,以具体用于根据以下公式确定所述目标光信号的色散值,<BR><MATHS id=cmaths0005 num="0005"><MATH><![CDATA[<mrow><MI>CD</MI><MO>=</MO><MROW><MO>(</MO><MI>x</MI><MO>+</MO><MSQRT><MSUP><MI>x</MI><MN>2</MN></MSUP><MO>-</MO><MSUBSUP><MROW><MN>4</MN><MI>T</MI></MROW><MN>0</MN><MN>4</MN></MSUBSUP></MSQRT><MO>)</MO></MROW><MO>&amp;CenterDot;</MO><MI>π</MI><MO>&amp;CenterDot;</MO><MFRAC><MI>c</MI><MSUP><MI>λ</MI><MN>2</MN></MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS><BR><MATHS id=cmaths0006 num="0006"><MATH><![CDATA[<mrow><MI>x</MI><MO>=</MO><MFRAC><MI>n</MI><MROW><MN>2</MN><MI>π</MI><MO>&amp;CenterDot;</MO><MI>SamplingRate</MI><MO>&amp;CenterDot;</MO><MI>SymbolRate</MI></MROW></MFRAC></MROW>]]&gt;</MATH></MATHS><BR>其中,CD表示色散值,n表示所述色散出现位置在所述自相关函数序列中的序号,T0表示所述目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示所述目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光信号的符号速率。</p></div> </div> </div> <div class="zlzy"> <div class="zltitle">说明书</div> <div class="gdyy"> <div class="gdyy_show"><p>说明书色散估计的方法、装置和系统 <BR>技术领域 <BR>本发明涉及通信领域,并且更具体地,涉及色散估计的方法、装置和系统。 <BR>背景技术 <BR>目前,已知一种色散估计的方法,即,对于未知色散的信号,以一定的步长间隔,扫描所有可能的色散值,并对于每个色散值通过设计代价函数(也可以称为算法设计)计算出对应的函数值,最后通过比较所有计算出的函数值(最大值或最小值),找到实际的色散值,实现色散估计。 <BR>但是,由于该方法需要扫描所有的色散值以进行计算,因此需要消耗大量的时间,不能快速估计信号的色散值。 <BR>发明内容 <BR>本发明实施例提供一种,能够快速估计信号的色散值。 <BR>第一方面,提供了一种色散估计的方法,该方法包括:计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列;根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的位置相对应;根据该脉冲位置,确定该目标光信号的色散值。 <BR>在一种可能的实施方式中,在该目标光信号为非归零码格式信号的情况下,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括:对该目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 <BR>结合第一方面和第一种可能的实施方式,在第二种可能的实施方式中,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括:对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 <BR>结合第一方面、第一种可能的实施方式和第二种可能的实施方式,在第三种可能的实施方式中,该根据该自相关函数序列,确定脉冲位置,包括:根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括该脉冲峰值;根据该第一自相关函数序列,确定该脉冲位置。 <BR>结合第一方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的实施方式,在第四种可能的实施方式中,该根据该第一自相关函数序列,确定该脉冲位置,包括:根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 <BR>结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式和第四种可能的实施方式,在第五种可能的实施方式中,该根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置,包括:根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列;根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理;根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 <BR>结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式和第五种可能的实施方式,在第六种可能的实施方式中,该根据第一序号和第一预设值,确定第二窗口,包括:根据第一序号和第一预设值,从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应;根据该第一范围中数值最大的第一元素和第二预设值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相对应;根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 <BR>结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式,在第七种可能的实施方式中,该根据该脉冲位置,确定该目标光信号的色散值,包括:根据以下公式确定该目标光信号的色散值, <BR><MATHS num="0001"><MATH><![CDATA[ <mrow><MI>CD</MI> <MO>=</MO> <MROW><MO>(</MO> <MI>x</MI> <MO>+</MO> <MSQRT><MSUP><MI>x</MI> <MN>2</MN> </MSUP><MO>-</MO> <MSUBSUP><MROW><MN>4</MN> <MI>T</MI> </MROW><MN>0</MN> <MN>4</MN> </MSUBSUP></MSQRT><MO>)</MO> </MROW><MO>&amp;CenterDot;</MO> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MFRAC><MI>c</MI> <MSUP><MI>λ</MI> <MN>2</MN> </MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0002"><MATH><![CDATA[ <mrow><MI>x</MI> <MO>=</MO> <MFRAC><MI>n</MI> <MROW><MN>2</MN> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MI>SamplingRate</MI> <MO>&amp;CenterDot;</MO> <MI>SymbolRate</MI> </MROW></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR>其中,CD表示色散值,n表示该脉冲位置在该自相关函数序列中的序号,T0表示该目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示该目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示该目标光信号的符号速率。 <BR>第二方面,提供了一种色散估计的装置,该装置包括:功率自相关处理单元,用于计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列,以及向脉冲位置搜索单元传输指示该自相关函数序列的信息;脉冲位置搜索单元,用于从该功率自相关处理单元获取指示该自相关函数序列的信息,并根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的位置相对应,以及向色散值确定单元传输指示该脉冲位置的信息;色散值确定单元,用于从该脉冲位置搜索单元获取指示该脉冲位置的信息,并根据该脉冲位置,确定该目标光信号的色散值。 <BR>在一种可能的实施方式中,在该目标光信号为非归零码格式信号的情况下,该装置还包括:第一干扰处理单元,用于对该目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 <BR>结合第二方面和第一种可能的实施方式,在第二种可能的实施方式中,该装置还包括:第二干扰处理单元,用于对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 <BR>结合第二方面、第一种可能的实施方式和第二种可能的实施方式,在第三种可能的实施方式中,该脉冲位置搜索单元具体用于根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括该脉冲峰值;用于根据该第一自相关函数序列,确定该脉冲位置。 <BR>结合第二方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的实施方式,在第四种可能的实施方式中,该脉冲位置搜索单元具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;用于根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 <BR>结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式和第四种可能的实施方式,在第五种可能的实施方式中,该脉冲位置搜索单元具体用于根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;用于根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列;用于根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理;用于根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 <BR>结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式和第五种可能的实施方式,在第六种可能的实施方式中,该脉冲位置搜索单元具体用于根据第一序号和第一预设值,从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应;用于根据该第一范围中数值最大的第一元素和第二预设值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相对应;用于根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 <BR>结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式,在第七种可能的实施方式中,该色散值确定单元具体用于根据以下公式确定该目标光信号的色散值, <BR><MATHS num="0003"><MATH><![CDATA[ <mrow><MI>CD</MI> <MO>=</MO> <MROW><MO>(</MO> <MI>x</MI> <MO>+</MO> <MSQRT><MSUP><MI>x</MI> <MN>2</MN> </MSUP><MO>-</MO> <MSUBSUP><MROW><MN>4</MN> <MI>T</MI> </MROW><MN>0</MN> <MN>4</MN> </MSUBSUP></MSQRT><MO>)</MO> </MROW><MO>&amp;CenterDot;</MO> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MFRAC><MI>c</MI> <MSUP><MI>λ</MI> <MN>2</MN> </MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0004"><MATH><![CDATA[ <mrow><MI>x</MI> <MO>=</MO> <MFRAC><MI>n</MI> <MROW><MN>2</MN> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MI>SamplingRate</MI> <MO>&amp;CenterDot;</MO> <MI>SymbolRate</MI> </MROW></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR>其中,CD表示色散值,n表示该脉冲位置在该自相关函数序列中的序号,T0表示该目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示该目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示该目标光信号的符号速率。 <BR>第三方面,提供了一种色散估计的系统,该系统包括:总线;与该总线相连的处理器;与该总线相连的存储器;其中,该处理器通过该总线,调用该存储器中存储的程序,以用于计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列;根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的出现位置相对应;根据该脉冲位置,确定该目标光信号的色散值。 <BR>在一种可能的实施方式中,在该目标光信号为非归零码格式信号的情况下,该处理器通过该总线,调用该存储器中存储的程序,以用于对该目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 <BR>结合第三方面和第一种可能的实施方式,在第二种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以用于对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 <BR>结合第三方面、第一种可能的实施方式和第二种可能的实施方式,在第三种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括该脉冲峰值;根据该第一自相关函数序列,确定该脉冲位置。 <BR>结合第三方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的实施方式,在第四种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;用于根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 <BR>结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式和第四种可能的实施方式,在第五种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;用于根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列;用于根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理;用于根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 <BR>结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式和第五种可能的实施方式,在第六种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据第一序号和第一预设值,从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应;用于根据该第一范围中数值最大的第一元素和第二预设值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相对应;用于根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 <BR>结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式,在第七种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据以下公式确定该目标光信号的色散值, <BR><MATHS num="0005"><MATH><![CDATA[ <mrow><MI>CD</MI> <MO>=</MO> <MROW><MO>(</MO> <MI>x</MI> <MO>+</MO> <MSQRT><MSUP><MI>x</MI> <MN>2</MN> </MSUP><MO>-</MO> <MSUBSUP><MROW><MN>4</MN> <MI>T</MI> </MROW><MN>0</MN> <MN>4</MN> </MSUBSUP></MSQRT><MO>)</MO> </MROW><MO>&amp;CenterDot;</MO> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MFRAC><MI>c</MI> <MSUP><MI>λ</MI> <MN>2</MN> </MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR><MATHS num="0006"><MATH><![CDATA[ <mrow><MI>x</MI> <MO>=</MO> <MFRAC><MI>n</MI> <MROW><MN>2</MN> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MI>SamplingRate</MI> <MO>&amp;CenterDot;</MO> <MI>SymbolRate</MI> </MROW></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR>其中,CD表示色散值,n表示该色散出现位置在该自相关函数序列中的序号,T0表示该目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示该目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示该目标光信号的符号速率。 <BR>根据本发明实施例的色散估计的方法、装置和系统,通过获取目标光信号的功率自相关函数,并确定该功率自相关函数出现脉冲峰值的位置,能够在无需扫描所有可能的色散值的情况下,根据该出现脉冲峰值的位置确定目标光信号的色散值,到达快速估计信号的色散值的目的。 <BR>附图说明 <BR>为了更清楚地说明本发明实施例的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。 <BR>图1是根据本发明实施例的色散估计的方法的示意性流程图。 <BR>图2是根据本发明实施例的色散估计的装置的示意性框图。 <BR>图3是根据本发明实施例的色散估计的方法的示意性框图。 <BR>具体实施方式 <BR>下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。 <BR>图1示出了根据本发明实施例的色散估计的方法100的示意性流程图,如图1所示,该方法100包括: <BR>S110,计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列; <BR>S120,根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的位置相对应; <BR>S130,根据该脉冲位置,确定该目标光信号的色散值。 <BR>在本发明实施例中,通过试验(例如,获取多个光信号的功率自相关函数,并获取该多个光信号的色散值)发现,在光信号(例如,单偏振的相移键控法(PSK,Phase&nbsp;Shift&nbsp;Keying)信号、偏振分割复用(PDM,PolarizationDivision&nbsp;Multiplexing)的PSK信号或正交振幅调制(QAM,QuadratureAmplitude&nbsp;Modulation)信号)存在色散时,该光信号的功率自相关函数会出现一个脉冲峰值,并且,不同的脉冲峰值出现位置与不同的色散值之间存在映射对应关系,因此,在本发明实施例中,可以根据光信号的信号功率的自相关函数中脉冲峰值的位置,确定该光信号的色散值。 <BR>因此,在S110中,可以计算目标光信号的信号功率(光功率)的自相关函数。在本发明实施例中,目标光信号为需要估算色散的光信号。在本发明实施例中,可以对该目标光信号进行以下处理,以单偏振系统为例,在相干接收机中,目标光信号和本振光在混频器(Hybrid)中混合,通过平衡接收机后转化为同向(in‑phase)电流和正交电流(quadrature)两路电流,经过两路高速模数转换器(ADC,Analog&nbsp;to&nbsp;Digital&nbsp;Converter)采样后,转换为数字信号,得到I信号(同向信号)和Q信号(正交信号),从而,可以根据该I信号和Q信号来确定目标光信号的信号功率的自相关函数。以下,作为示例而非限定,以该通过该I信号和Q信号确定目标光信号的信号功率的自相关函数为例进行说明。 <BR>可选地,在本发明实施例中,在该目标光信号为非归零码格式信号的情况下,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括: <BR>对该目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 <BR>具体地说,零频干扰是指,在该目标光信号为非归零(NRZ,No&nbsp;ReturnZero)码格式的情况下,如果直接获取该格式的信号的自相关函数,则获取的自相关函数的零频位置会出现一个较大的峰值,从而对后述搜索脉冲峰值造成影响,因此,作为示例而非限定,在本发明实施例中,可以对I信号和Q信号进行延迟处理和相减处理,以实现对该目标光信号的高通滤波,其中,该延迟处理使用的迟延大小可以为例如,一个符号周期。从而,能够避免获取的自相关函数在零频位置出现较大的峰值。 <BR>可选地,在本发明实施例中,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括: <BR>对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 <BR>具体地说,作为示例而非限定,可以对I信号和Q信号(具体地说,是与该I信号相对应的I信号序列以及与该Q信号相对应的Q信号序列)进行求平方和处理,然后用经过该求平方和处理后的信号减去信号均值,其中,作为示例而非限定,该信号均值可以通过以下方法确定,即,可以获取与该I信号相对应的I信号序列,以及与该Q信号相对应的Q信号序列,其后,对I信号序列和Q信号序列求和,然后除以I信号序列和Q信号序列的长度,从而确定该信号均值。通过以上处理,能够实现对该目标光信号的隔直处理,从而能够去掉直流分量产生的无用的自相关函数部分。 <BR>其后,在S110中,可以对经过上述处理的目标光信号进行傅里叶变换,再做幅度平方,最后进行逆傅里叶变换,得到该目标光信号的信号功率的自相关函数,进而可以根据该自相关函数,生成以该自相关函数的各函数值为元素的自相关函数序列(以下,为了便于说明,记作R1[n])。并且自相关函数以表示为: <BR><MATHS num="0007"><MATH><![CDATA[ <mrow><MSUB><MOVER><MI>R</MI> <MO>&amp;OverBar;</MO> </MOVER><MI>yy</MI> </MSUB><MO>[</MO> <MI>n</MI> <MO>]</MO> <MO>=</MO> <MI>IFFT</MI> <MO>{</MO> <MSUP><MROW><MO>|</MO> <MI>FFT</MI> <MROW><MROW><MO>(</MO> <MI>y</MI> <MO>[</MO> <MI>n</MI> <MO>]</MO> <MO>)</MO> </MROW><MO>|</MO> </MROW></MROW><MN>2</MN> </MSUP><MO>}</MO> <MO>=</MO> <MFRAC><MN>1</MN> <MI>N</MI> </MFRAC><MO>&amp;CenterDot;</MO> <MI>FFT</MI> <MO>{</MO> <MSUP><MROW><MO>|</MO> <MI>FFT</MI> <MROW><MO>(</MO> <MI>y</MI> <MO>[</MO> <MI>n</MI> <MO>]</MO> <MO>)</MO> </MROW><MO>|</MO> </MROW><MN>2</MN> </MSUP><MO>}</MO> </MROW>]]&gt;</MATH></MATHS> <BR>其中,表示目标光信号的自相关函数,FFT表示快速傅里叶变换,IFFT表示快速反傅里叶变换,y[n]表示目标光信号的信号功率(光功率),R1[n]为根据该自相关函数生成的函数序列(自相关函数序列)。 <BR>在S120中,可以从该R1[n]中搜索脉冲峰值的位置(脉冲位置),即,确定该脉冲峰值在R1[n]中的序号。 <BR>可选地,在本发明实施例中,可以搜索R1[n]中数值最大的元素(记作元素a),并将该元素a的位置确定为脉冲峰值的位置,并可以记录该元素a的序号,以在后述确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说元素a在R1[n]中的序号(记作n1),确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,该根据该自相关函数序列,确定脉冲位置,包括: <BR>根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括与该脉冲峰值相对应的元素; <BR>根据该第一自相关函数序列,确定该脉冲位置。 <BR>具体地说,在本发明实施例中,可以从该R1[n]中选取一个区间作为搜索范围(即,第一自相关函数序列,以下,为了便于说明,记作R2[n]),并且,在本发明实施例中,可以将R1[n]作为R2[n],也可以将{R1[N1]:R1[N2]}(即,R1[n]中从序号为N1的元素到序号为N2的元素)作为R2[n],其中,N1表示R1[n]中选取的开始位置(序号),N2表示R1[n]中选取的结束位置(序号)。 <BR>可选地,在本发明实施例中,该根据该自相关函数序列,确定第一自相关函数序列,包括: <BR>根据用于获取所述目标光信号的相干接收机的电带宽和用于获取所述目标光信号的高速模数转换器的采样率,从所述自相关函数序列中确定选取所述第一自相关函数序列的开始位置。 <BR>具体地说,可以根据以下式1确定N1,即: <BR>N1&gt;floor(2×SamplingRate/Be)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式1 <BR>其中,floor为向下取整数函数,Be表示接收机电带宽,SamplingRate是表示ADC采样率,通常情况下,Be&gt;10GHz,SamplingRate&gt;20G&nbsp;Sa/s。 <BR>并且,在本发明实施例中,N2可以表示能够进行色散估计的最大范围,在本发明实施例中,N2&gt;N1。 <BR>可选地,在本发明实施例中,该根据该第一自相关函数序列,确定该脉冲位置,包括: <BR>根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理; <BR>根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 <BR>具体地说,在确定R2[n]之后,对R2[n]进行滑动窗口平均处理(第一滑动平均处理),并且,作为示例而非限定,该滑动窗口平均处理的窗口宽度(记作width1)可以为例如,3。在本发明实施例中,滑动窗口平均处理的方法与现有技术相同,这里,为了避免赘述,省略其说明。在对R2[n]进行滑动窗口平均处理后(以下,为了便于理解,将处理后生成的序列记作R3[n]),可以搜索R3[n]中数值最大的元素(记作元素b),并将该元素b的位置确定为脉冲峰值的位置,并可以记录该元素b在R3[n]中的序号(记作n2),从而,可以根据以下式2确定n1,即 <BR>n1=n2+N1-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式2 <BR>通过对该R2[n](第一自相关函数序列)进行滑动窗口平均处理,能够减少自相关函数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。以上列举的作为width1的数值仅为本发明的一个示例性说明,并不限定于此。 <BR>可选地,在本发明实施例中,在确定R2[n]之后,也可以直接搜索R2[n]中数值最大的元素(记作元素c),并将该元素c的位置确定为脉冲峰值的位置,并可以记录该元素c在R2[n]中的序号(记作n3),从而,可以根据以下式3确定n1,即 <BR>n1=n3+N1-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式3 <BR>从而,在后述确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,该根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置,包括: <BR>根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号; <BR>根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列; <BR>根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理; <BR>根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 <BR>在本发明实施例中,可以将自相关函数中与产生的色散相对应的脉冲峰值的全宽半高(FWHM,Full&nbsp;Width&nbsp;Half&nbsp;Maximum)的一半,即,将脉冲峰值(最大值)的二分之一作为门限(后述第一门限),寻找脉冲宽度,再将脉冲宽度的例如,二分之一作为第二窗口,具体地说,在确定n2后,可以从R2[n]中,在以n2为中心的预设范围(与第一预设值相对应,例如,可以为5个单位,这里,“单位”可以是单位序号)内,搜索数值最大的元素(记作元素d),并根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之一作为该第一门限,其后,以元素d为中心,向左侧搜索第一个数值小于该第一门限的元素(记作元素e),并记录该元素e的序号(记作t1);向右侧搜索第一个数值小于该第一门限的元素(记作元素f),并记录该元素f的序号(记作t2),可以根据以下式4,确定滑动窗口宽度(第二窗口的一例,记作width2),即 <BR>width2=floor{(t2-t1)/2}-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式4 <BR>其中,floor为向下取整数函数。应理解,以上列举的以脉冲峰值(最大值)的二分之一作为门限,确定第二窗口的方式仅为本发明的一个示例性说明,本发明并不限定于此。 <BR>在确定了width2(第二窗口)后,可以从该R2[n]中选取一个区间作为搜索范围(即,第二自相关函数序列,以下,为了便于说明,记作R4[n]),在本发明实施例中,可以将R2[n]作为R4[n],也可以将{R2[N3]:R2[N4]}(即,R2[n]中从序号为N3的元素到序号为N4的元素)作为R4[n],其中,N3表示R2[n]中选取的开始位置(序号),N4表示R2[n]中选取的结束位置(序号)。 <BR>具体地说,可以根据以下式5确定N3,即: <BR>N3=t1-width2-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式5 <BR>可以根据以下式6确定N4,即: <BR>N4=t2+width2+1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式6 <BR>在确定了R4[n]后,根据width2,对该R4[n]进行滑动平均处理(第二滑动平均处理),以下,为了便于理解,将处理后生成的序列记作R5[n],可以搜索R5[n]中数值最大的元素(记作元素g),并将该元素g的位置确定为脉冲峰值的位置,并可以记录该元素g在R5[n]中的序号(记作n4),由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,可以根据以下式7确定n1,即 <BR>n1=n4+t1-width2+N1-3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式7 <BR>其中,width2=floor{(t2-t1)/2}-1。 <BR>从而,在后述确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>通过对该R2[n](第一自相关函数序列)进行滑动窗口平均处理,能够减少自相关函数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。 <BR>可选地,在本发明实施例中,该根据第一序号和第一预设值,确定第二窗口,包括: <BR>根据第一序号和第一预设值,从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应; <BR>根据该第一范围中数值最大的第一元素和第二预设值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相对应; <BR>根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 <BR>具体地说,如果偏振模色散(PMD,Polarization&nbsp;Mode&nbsp;Dispersion)较大,则可能会出现一个自相关函数中的脉冲(峰值)分裂为两个脉冲(峰值)的情况,此情况下,如果仅以脉冲峰值(最大值)的位置为中心寻找脉冲宽度,则寻找到的脉冲宽度值会产生偏差,因此,在本发明实施例中,在确定n2后,还可以从R2[n]中,在以n2为中心的预设范围(即,第一范围,与第一预设值相对应,例如,可以为5个单位,这里为“单位”可以是单位序号)内,搜索数值最大的元素(即,元素d)。 <BR>其后,在以元素d在R2[n]中的序号为中心的预设范围(即,第二范围,与第二预设值相对应,例如,可以为10个单位,这里为“单位”可以为单位序号),搜索数值第二大的元素(记作元素h)。 <BR>其后,确定元素d与元素h在R2[n]中的位置关系,这里,为了便于理解,以元素d在元素h的右侧为例(即,元素d在R2[n]中的序号大于元素h在R2[n]中的序号)进行说明。 <BR>其后,根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之一作为该第一门限。 <BR>其后,以元素h为中心,向左侧搜索第一个数值小于该第一门限的元素(记作元素i,并记录该元素i的序号(记作t3);以元素d为中心,向右侧搜索第一个数值小于该第一门限的元素(记作元素j),并记录该元素j的序号(记作t4),可以根据以下式8,确定一个中间参量(记作width3),即 <BR>Width3=floor{(t4-t3)/2}-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式8 <BR>其中,floor为向下取整数函数。应理解,以上列举的以脉冲峰值(最大值)的二分之一作为门限,确定第二窗口的方式仅为本发明的一个示例性说明,本发明并不限定于此。 <BR>在确定了width3后,可以从该R2[n]中选取一个区间作为搜索范围(即,第二自相关函数序列,以下,为了便于说明,记作R6[n]),在本发明实施例中,可以将R2[n]作为R6[n],也可以将{R2[N5]:R2[N6]}(即,R2[n]中从序号为N5的元素到序号为N6的元素)作为R6[n],其中,N5表示R2[n]中选取的开始位置(序号),N6表示R2[n]中选取的结束位置(序号)。 <BR>具体地说,可以根据以下式9确定N5,即: <BR>N5=t3-width3-1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式9 <BR>可以根据以下式10确定N6,即: <BR>N6=t4+width3+1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式10 <BR>在确定了R6[n]后,可以根据以下式11,确定对R6[n]进行滑动平均处理时使用的窗口宽度(第二窗口的另一例,记作width4),即 <BR>width4=2×width3+1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式11 <BR>其后,对该R6[n]进行滑动平均处理(第二滑动平均处理),以下,为了便于理解,将处理后生成的序列记作R7[n],可以从该R7[n]中选取一个区间作为搜索范围(以下,为了便于说明,记作R8[n]),在本发明实施例中,可以将R7[n]作为R8[n],也可以将{R7[N7]:R7[N8]}(即,R7[n]中从序号为N7的元素到序号为N8的元素)作为R8[n],其中,N7表示R7[n]中选取的开始位置(序号),N8表示R7[n]中选取的结束位置(序号)。 <BR>具体地说,可以根据以下式12确定N7,即: <BR>N7=2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式12 <BR>可以根据以下式13确定N8,即: <BR>N8=t4-t3+2×width3+2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式13 <BR>在确定R8[n]后,可以搜索R8[n]中数值最大的元素(记作元素k),并记录该元素k在R8[n]中的序号(记作n5)。 <BR>可选地,在本发明实施例中,可以将该元素k的位置确定为脉冲峰值的位置,由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,可以根据以下式14确定n1,即 <BR>n1=n5+t3‑width3+N1‑3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;式14 <BR>其中,Width3=floor{(t4‑t3)/2}‑1 <BR>从而,在后述确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,在确定n5后,还可以利用R8[n5‑1]、R8[n5]、R8[n5+1]三点做抛物线拟合,从而,可以根据以下式15,确定该抛物线顶点在R8[n]中的位置(序号,记作n6),即<MATHS num="0008"><MATH><![CDATA[ <mrow> <MSUB><MI>n</MI> <MN>6</MN> </MSUB><MO>=</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MFRAC><MROW><MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>-</MO> <MN>1</MN> <MO>]</MO> <MO>-</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MN>1</MN> <MO>]</MO> </MROW><MROW><MN>2</MN> <MROW><MO>(</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>-</MO> <MN>1</MN> <MO>]</MO> <MO>+</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MN>1</MN> <MO>]</MO> <MO>-</MO> <MSUB><MROW><MN>2</MN> <MI>R</MI> </MROW><MN>5</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>]</MO> <MO>)</MO> </MROW></MROW></MFRAC></MROW>]]&gt;</MATH></MATHS>式15 <BR>其后,可以根据n6和以下式16,确定该抛物线顶点在R2[n]中的位置(序号),即<MATHS num="0009"><MATH><![CDATA[ <mrow> <MSUB><MI>n</MI> <MN>7</MN> </MSUB><MO>=</MO> <MSUB><MI>n</MI> <MROW><MN>5</MN> <MO>+</MO> </MROW></MSUB><MSUB><MI>t</MI> <MN>3</MN> </MSUB><MO>-</MO> <MSUB><MI>width</MI> <MN>3</MN> </MSUB><MO>-</MO> <MN>2</MN> <MO>+</MO> <MFRAC><MROW><MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>-</MO> <MN>1</MN> <MO>]</MO> <MO>-</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MN>1</MN> <MO>]</MO> </MROW><MROW><MN>2</MN> <MROW><MO>(</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>-</MO> <MN>1</MN> <MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MN>1</MN> <MO>]</MO> <MO>-</MO> <MSUB><MROW><MN>2</MN> <MI>R</MI> </MROW><MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>]</MO> </MROW></MFRAC></MROW>]]&gt;</MATH></MATHS>式16 <BR>由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,可以根据以下式17确定n1,即<MATHS num="0010"><MATH><![CDATA[ <mrow> <MSUB><MI>n</MI> <MN>1</MN> </MSUB><MO>=</MO> <MSUB><MI>n</MI> <MROW><MN>5</MN> <MO>+</MO> </MROW></MSUB><MSUB><MI>t</MI> <MN>3</MN> </MSUB><MO>-</MO> <MSUB><MI>width</MI> <MN>3</MN> </MSUB><MO>+</MO> <MFRAC><MROW><MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>-</MO> <MN>1</MN> <MO>]</MO> <MO>-</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MN>1</MN> <MO>]</MO> </MROW><MROW><MN>2</MN> <MROW><MO>(</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>-</MO> <MN>1</MN> <MO>)</MO> </MROW><MO>+</MO> <MSUB><MI>R</MI> <MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>+</MO> <MN>1</MN> <MO>]</MO> <MO>-</MO> <MSUB><MROW><MN>2</MN> <MI>R</MI> </MROW><MN>8</MN> </MSUB><MO>[</MO> <MSUB><MI>n</MI> <MN>5</MN> </MSUB><MO>]</MO> </MROW></MFRAC><MO>+</MO> <MSUB><MI>N</MI> <MN>1</MN> </MSUB><MO>-</MO> <MN>3</MN> </MROW>]]&gt;</MATH></MATHS>式17 <BR>其中,width3=floor{(t4-t3)/2}-1。 <BR>从而,在S130中,根据在S120中确定的自相关函数中的脉冲峰值出现位置(即,n1),确定目标光信号的色散值。 <BR>在本发明实施例中,根据该脉冲位置,确定该目标光信号的色散值,包括: <BR>根据离线数据,确定指示脉冲峰值位置与色散值的对应关系的指示信息, <BR>根据该脉冲位置和该指示信息,确定该光信号的色散值。 <BR>具体地说,作为示例而非限定,例如,可以分别对多个光信号(离线数据的一例)进行色散值测试等试验,记录各光信号的色散值和它们各自对应的自相关函数中脉冲峰值出现的位置,作为训练数据,从而,可以采用例如期望最大化算法,通过逐步改进公式的参数,使参数和训练数据的似然概率逐渐增大,最后终止于一个极大点。直观地理解,期望最大化算法也可被看作为一个逐次逼近算法:事先并不知道模型的参数,可以随机的选择一套参数或者事先粗略地给定某个初始参数λ0,确定出对应于这组参数的最可能的状态,计算每个训练样本的可能结果的概率,在当前的状态下再由样本对参数修正,重新估计参数λ,并在新的参数下重新确定模型的状态,这样,通过多次的迭代,循环直至某个收敛条件满足为止,就可以使得模型的参数逐渐逼近真实参数。因此,该步骤可以以离线的方式自动进行。将离线计算得到的上述各种脉冲峰值位置与各自对应的色散值预先存储起来,例如以表的形式存储起来。应理解,以上列举的各公式以及各公式中的常数的训练方法仅是本发明的一个实施例,本发明并不限定于此。 <BR>作为示例而非限定,以下式18示出了色散值与自相关函数中脉冲峰值出现的位置的对应关系的一例,即<MATHS num="0011"><MATH><![CDATA[ <mrow> <MI>CD</MI> <MO>=</MO> <MROW><MO>(</MO> <MI>x</MI> <MO>+</MO> <MSQRT><MSUP><MI>x</MI> <MN>2</MN> </MSUP><MO>-</MO> <MSUBSUP><MROW><MN>4</MN> <MI>T</MI> </MROW><MN>0</MN> <MN>4</MN> </MSUBSUP></MSQRT><MO>)</MO> </MROW><MO>&amp;CenterDot;</MO> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MFRAC><MI>c</MI> <MSUP><MI>λ</MI> <MN>2</MN> </MSUP></MFRAC></MROW>]]&gt;</MATH></MATHS>式18 <BR><MATHS num="0012"><MATH><![CDATA[ <mrow><MI>x</MI> <MO>=</MO> <MFRAC><MI>n</MI> <MROW><MN>2</MN> <MI>π</MI> <MO>&amp;CenterDot;</MO> <MI>SamplingRate</MI> <MO>&amp;CenterDot;</MO> <MI>SymbolRate</MI> </MROW></MFRAC></MROW>]]&gt;</MATH></MATHS> <BR>其中,CD表示色散值,n表示所述脉冲位置在该自相关函数序列中的序号,T0表示所述目标光信号的初始发射机光脉冲宽度,c表示真空光速,λ表示所述目标光信号的波长,SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光信号的符号速率。 <BR>应理解,以上列举的各公式、参数以及数值仅为本发明的示例性说明,本发明并不限定于此,其他能够用于确定光信号的自相关函数的脉冲峰值出现位置的各公式、参数以及数值均落入本发明的保护范围内。 <BR>根据本发明实施例的色散估计的方法,通过获取目标光信号的光信号功率自相关函数,并确定该光信号功率自相关函数出现脉冲峰值的位置,能够在无需扫描所有可能的色散值的情况下,根据该出现脉冲峰值的位置确定目标光信号的色散值,到达快速估计信号的色散值的目的。 <BR>上文中,结合图1,详细描述了根据本发明实施例的色散估计的方法,下面将结合图2,详细描述根据本发明实施例的色散估计的装置。 <BR>图2示出了根据本发明实施例的色散估计的装置200的示意性框图。如图2所示,该装置可以包括: <BR>功率自相关处理单元210,用于计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列,以及向脉冲位置搜索单元220传输指示该自相关函数序列的信息; <BR>脉冲位置搜索单元220,用于从该功率自相关处理单元210获取指示该自相关函数序列的信息,并根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的位置相对应,以及向色散值确定单元230传输指示该脉冲位置的信息; <BR>色散值确定单元230,用于从该脉冲位置搜索单元220指示该脉冲位置的信息,并根据该脉冲位置,确定该目标光信号的色散值。 <BR>在本发明实施例中,通过试验(例如,获取多个光信号的功率自相关函数,并获取该多个光信号的色散值)发现,在光信号(例如,单偏振的PSK信号、PDM的PSK信号或QAM信号存在色散时,该光信号功率的自相关函数会出现一个脉冲峰值,并且,不同的脉冲峰值位置与不同的色散值之间存在映射对应关系,因此,在本发明实施例中,可以根据光信号的信号功率的自相关函数中脉冲峰值的位置,确定该光信号的色散值。 <BR>因此,功率自相关处理单元210可以计算目标光信号的信号功率(光功率)的自相关函数。在本发明实施例中,目标光信号为需要估算色散的光信号。在本发明实施例中,可以对该目标光信号进行以下处理,以单偏振系统为例,在相干接收机中,目标光信号和本振光在Hybrid中混合,通过平衡接收机后转化为两路电流(同向电流和正交电流),经过两路ADC采样后,转换为数字信号,得到I信号和Q信号,从而,可以根据该I信号和Q信号来确定目标光信号的信号功率的自相关函数。以下,作为示例而非限定,以该通过该I信号和Q信号确定目标光信号的信号功率的自相关函数为例进行说明。 <BR>可选地,在本发明实施例中,在该目标光信号为非归零码格式信号的情况下,该装置还包括: <BR>第一干扰处理单元,用于对该目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 <BR>具体地说,零频干扰是指,在该目标光信号为非归零(NRZ,No&nbsp;ReturnZero)码格式的情况下,如果直接获取该格式的信号的自相关函数,则获取的自相关函数的零频位置会出现一个较大的峰值,从而对后述搜索脉冲峰值造成影响,因此,作为示例而非限定,在本发明实施例中,可以对I信号和Q信号进行延迟处理和相减处理,以实现对该目标光信号的高通滤波,其中,该延迟处理使用的迟延大小可以为例如,一个符号周期。从而,能够避免获取的自相关函数在零频位置出现较大的峰值。 <BR>可选地,在本发明实施例中,该装置还包括: <BR>第二干扰处理单元,用于对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 <BR>具体地说,作为示例而非限定,第二干扰处理单元可以对I信号和Q信号(具体地说,是与该I信号相对应的I信号序列以及与该Q信号相对应的Q信号序列)进行求平方和处理,然后用经过该求平方和处理后的信号减去信号均值,其中,作为示例而非限定,该信号均值可以通过以下方法确定,即,可以获取与该I信号相对应的I信号序列,以及与该Q信号相对应的Q信号序列,其后,对I信号序列和Q信号序列求和,然后除以I信号序列和Q信号序列的长度,从而确定该信号均值。通过以上处理,能够实现对该目标光信号的隔直处理,从而能够去掉直流分量产生的无用的自相关函数部分。 <BR>其后,功率自相关处理单元210可以对经过上述处理的目标光信号进行傅里叶变换,再做幅度平方,最后进行逆傅里叶变换,得到该目标光信号的信号功率的自相关函数,进而可以根据该自相关函数,生成以该自相关函数的各函数值为元素的自相关函数序列(以下,为了便于说明,记作R1[n]) <BR>脉冲位置搜索单元220可以从该R1[n]中搜索脉冲峰值的位置,即,确定该脉冲峰值在R1[n]中的序号。 <BR>可选地,在本发明实施例中,脉冲位置搜索单元220可以搜索R1[n]中数值最大的元素(记作元素a),并将该元素a的位置确定为脉冲峰值的位置,并可以记录该元素a的序号,以在后述色散值确定单元230确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说元素a在R1[n]中的序号(记作n1),确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,该脉冲位置搜索单元220具体用于根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括与该脉冲峰值相对应的元素; <BR>用于根据该第一自相关函数序列,确定该脉冲位置。 <BR>具体地说,在本发明实施例中,脉冲位置搜索单元220可以从该R1[n]中选取一个区间作为搜索范围(即,第一自相关函数序列,以下,为了便于说明,记作R2[n]),并且,在本发明实施例中,脉冲位置搜索单元220可以将R1[n]作为R2[n],也可以将{R1[N1]:R1[N2]}(即,R1[n]中从序号为N1的元素到序号为N2的元素)作为R2[n],其中,N1表示R1[n]中选取的开始位置(序号),N2表示R1[n]中选取的结束位置(序号)。 <BR>可选地,在本发明实施例中,该脉冲位置搜索单元220具体用于根据用于获取该目标光信号的相干接收机的电带宽和用于获取该目标光信号的高速模数转换器的采样率,从该自相关函数序列中确定选取该第一自相关函数序列的开始位置。 <BR>具体地说,脉冲位置搜索单元220可以根据式1确定N1。 <BR>可选地,在本发明实施例中,该脉冲位置搜索单元220具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理; <BR>用于根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 <BR>具体地说,脉冲位置搜索单元220在确定R2[n]之后,可以对R2[n]进行滑动窗口平均处理(第一滑动平均处理),并且,作为示例而非限定,该滑动窗口平均处理的窗口宽度(记作width1)可以为例如,3。在本发明实施例中,滑动窗口平均处理的方法与现有技术相同,这里,为了避免赘述,省略其说明。在对R2[n]进行滑动窗口平均处理后(以下,为了便于理解,将处理后生成的序列记作R3[n]),可以搜索R3[n]中数值最大的元素(记作元素b),并将该元素b的位置确定为脉冲峰值的位置,并可以记录该元素b在R3[n]中的序号(记作n2),从而,可以根据式2确定n1。 <BR>通过对该R2[n](第一自相关函数序列)进行滑动窗口平均处理,能够减少自相关函数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。以上列举的作为width1的数值仅为本发明的一个示例性说明,并不限定于此。 <BR>可选地,在本发明实施例中,在确定R2[n]之后,脉冲位置搜索单元220也可以直接搜索R2[n]中数值最大的元素(记作元素c),并将该元素c的位置确定为脉冲峰值的位置,并可以记录该元素c在R2[n]中的序号(记作n3),从而,可以根据式3确定n1。 <BR>从而,在后述确定色散值时,色散值确定单元230可以根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,该脉冲位置搜索单元220具体用于根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号; <BR>用于根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列; <BR>用于根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理; <BR>用于根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 <BR>在本发明实施例中,脉冲位置搜索单元220例如,可以将自相关函数中与产生的色散相对应的脉冲峰值的FWHM的一半,即,将脉冲峰值(最大值)的二分之一作为门限(后述第一门限),寻找脉冲宽度,再将脉冲宽度的例如,二分之一作为第二窗口,具体地说,在确定n2后,可以从R2[n]中,在以n2的中心的预设范围(与第一预设值相对应,例如,可以为5个单位,这里“单位”可以是单位序号)内,搜索数值最大的元素(记作元素d),并根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之一作为该第一门限,其后,以元素d为中心,向左侧搜索第一个数值小于该第一门限的元素(记作元素e),并记录该元素e的序号(记作t1);向右侧搜索第一个数值小于该第一门限的元素(记作元素f),并记录该元素f的序号(记作t2),可以根据式4,确定滑动窗口宽度(第二窗口的一例,记作width2)。 <BR>在确定了width2(第二窗口)后,脉冲位置搜索单元220可以从该R2[n]中选取一个区间作为搜索范围(即,第二自相关函数序列,以下,为了便于说明,记作R4[n]),在本发明实施例中,可以将R2[n]作为R4[n],也可以将{R2[N3]:R2[N4]}(即,R2[n]中从序号为N3的元素到序号为N4的元素)作为R4[n],其中,N3表示R2[n]中选取的开始位置(序号),N4表示R2[n]中选取的结束位置(序号)。具体地说,脉冲位置搜索单元220可以根据式5确定N3,并可以根据式6确定N4。 <BR>在确定了R4[n]后,脉冲位置搜索单元220根据width2,对该R4[n]进行滑动平均处理(第二滑动平均处理),以下,为了便于理解,将处理后生成的序列记作R5[n],可以搜索R5[n]中数值最大的元素(记作元素g),并将该元素g的位置确定为脉冲峰值的位置,并可以记录该元素g在R5[n]中的序号(记作n4),由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,可以根据式7确定n1。 <BR>从而,在后述确定色散值时,色散值确定单元230可以根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>通过对该R2[n](第一自相关函数序列)进行滑动窗口平均处理,能够减少自相关函数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。 <BR>可选地,在本发明实施例中,该脉冲位置搜索单元220具体用于根据第一序号和第一预设值,从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应; <BR>用于根据该第一范围中数值最大的第一元素和第二预设值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相对应; <BR>用于根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 <BR>具体地说,如果PMD较大,则可能会出现一个自相关函数中的脉冲(峰值)分裂为两个脉冲(峰值)的情况,此情况下,如果仅以脉冲峰值(最大值)的位置为中心寻找脉冲宽度,则寻找到的脉冲宽度值会产生偏差,因此,在本发明实施例中,脉冲位置搜索单元220在确定n2后,还可以从R2[n]中,在以n2为中心的预设范围(即,第一范围,与第一预设值相对应,例如,可以为5个单位)内,搜索数值最大的元素(即,元素d)。 <BR>其后,在以元素d在R2[n]中的序号为中心的预设范围(即,第二范围,与第二预设值相对应,例如,可以为10个单位)内,搜索数值第二大的元素(记作元素h)。 <BR>其后,确定元素d与元素h在R2[n]中的位置关系,这里,为了便于理解,以元素d在元素h的右侧为例(即,元素d在R2[n]中的序号大于元素h在R2[n]中的序号)进行说明。 <BR>其后,根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之一作为该第一门限。 <BR>其后,以元素h为中心,向左侧搜索第一个数值小于该第一门限的元素(记作元素i,并记录该元素i的序号(记作t3);以元素d为中心,向右侧搜索第一个数值小于该第一门限的元素(记作元素j),并记录该元素j的序号(记作t4),可以根据式8,确定一个中间参量(记作width3)。 <BR>脉冲位置搜索单元220在确定了width3后,可以从该R2[n]中选取一个区间作为搜索范围(即,第二自相关函数序列,以下,为了便于说明,记作R6[n]),在本发明实施例中,可以将R2[n]作为R6[n],也可以将{R2[N5]:R2[N6]}(即,R2[n]中从序号为N5的元素到序号为N6的元素)作为R6[n],其中,N5表示R2[n]中选取的开始位置(序号),N6表示R2[n]中选取的结束位置(序号)。具体地说,脉冲位置搜索单元220可以根据式9确定N5,并可以根据式10确定N6。 <BR>脉冲位置搜索单元220在确定了R6[n]后,可以根据式11,确定对R6[n]进行滑动平均处理时使用的窗口宽度(第二窗口的另一例,记作width4)。 <BR>其后,脉冲位置搜索单元220对该R6[n]进行滑动平均处理(第二滑动平均处理),以下,为了便于理解,将处理后生成的序列记作R7[n],可以从该R7[n]中选取一个区间作为搜索范围(以下,为了便于说明,记作R8[n]),在本发明实施例中,可以将R7[n]作为R8[n],也可以将{R7[N7]:R7[N8]}(即,R7[n]中从序号为N7的元素到序号为N8的元素)作为R8[n],其中,N7表示R7[n]中选取的开始位置(序号),N8表示R7[n]中选取的结束位置(序号)。具体地说,脉冲位置搜索单元220可以根据式12确定N7,并可以根据式13确定N8。 <BR>脉冲位置搜索单元220在确定R8[n]后,可以搜索R8[n]中数值最大的元素(记作元素k),并记录该元素k在R8[n]中的序号(记作n5)。 <BR>在本发明实施例中,可以将该元素k的位置确定为脉冲峰值的位置,由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,脉冲位置搜索单元220可以根据式14确定n1。 <BR>从而,在后述色散值确定单元230确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,脉冲位置搜索单元220在确定n5后,还可以利用R8[n5‑1]、R8[n5]、R8[n5+1]三点做抛物线拟合,从而,可以根据式15,确定该抛物线顶点在R8[n]中的位置(序号,记作n6)。 <BR>其后,脉冲位置搜索单元220可以根据n6和式16,确定该抛物线顶点在R2[n]中的位置(序号)。 <BR>由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,可以根据式17确定n1。 <BR>从而,色散值确定单元230可以根据脉冲位置搜索单元220确定的自相关函数中的脉冲峰值出现位置(即,n1),确定目标光信号的色散值。 <BR>在本发明实施例中,该色散值确定单元230具体用于根据离线数据,确定指示脉冲峰值位置与色散值的对应关系的指示信息, <BR>根据该脉冲位置和该指示信息,确定该光信号的色散值。 <BR>具体地说,作为示例而非限定,例如,色散值确定单元230可以分别对多个光信号(离线数据的一例)进行色散值测试等试验,记录各光信号的色散值和它们各自对应的自相关函数中脉冲峰值出现的位置,作为训练数据,从而,可以采用例如期望最大化算法,通过逐步改进公式的参数,使参数和训练数据的似然概率逐渐增大,最后终止于一个极大点。直观地理解,期望最大化算法也可被看作为一个逐次逼近算法:事先并不知道模型的参数,可以随机的选择一套参数或者事先粗略地给定某个初始参数λ0,确定出对应于这组参数的最可能的状态,计算每个训练样本的可能结果的概率,在当前的状态下再由样本对参数修正,重新估计参数λ,并在新的参数下重新确定模型的状态,这样,通过多次的迭代,循环直至某个收敛条件满足为止,就可以使得模型的参数逐渐逼近真实参数。因此,该步骤可以以离线的方式自动进行。将离线计算得到的上述各种脉冲峰值位置与各自对应的色散值预先存储起来,例如以表的形式存储起来。应理解,以上列举的各公式以及各公式中的常数的训练方法仅是本发明的一个实施例,本发明并不限定于此。 <BR>作为示例而非限定,色散值确定单元230可以根据式18确定色散值。 <BR>应理解,以上列举的各公式、参数以及数值仅为本发明的示例性说明,本发明并不限定于此,其他能够用于确定光信号的自相关函数的脉冲峰值出现位置的各公式、参数以及数值均落入本发明的保护范围内。 <BR>根据本发明实施例的色散估计的装置可对应于本发明实施例的色散估计的方法的执行主体,并且,该色散估计的装置中的各单元和上述其他操作和/或功能分别为了实现图1中的方法100的相应流程,为了简洁,在此不再赘述。 <BR>根据本发明实施例的色散估计的装置,通过获取目标光信号的光信号功率自相关函数,并确定该光信号功率自相关函数出现脉冲峰值的位置,能够在无需扫描所有可能的色散值的情况下,根据该出现脉冲峰值的位置确定目标光信号的色散值,到达快速估计信号的色散值的目的。 <BR>上文中,结合图1和图2,详细描述了根据本发明实施例的色散估计的方法和装置,下面将结合图3,详细描述根据本发明实施例的色散估计的系统。 <BR>图3示出了根据本发明实施例的色散估计的装置300的示意性框图。如图3所示,该系统可以包括: <BR>总线310; <BR>与该总线相连的处理器320; <BR>与该总线相连的存储器330; <BR>其中,该处理器320通过该总线,调用该存储器330中存储的程序,以用于计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列; <BR>根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的出现位置相对应; <BR>根据该脉冲位置,确定该目标光信号的色散值。 <BR>在本发明实施例中,通过试验(例如,获取多个光信号的功率自相关函数,并获取该多个光信号的色散值)发现,在光信号(例如,单偏振的PSK信号、PDM的PSK信号或QAM信号存在色散时,该光信号功率的自相关函数会出现一个脉冲峰值,并且,不同的脉冲峰值位置与不同的色散值之间存在映射对应关系,因此,在本发明实施例中,可以根据光信号的信号功率的自相关函数中脉冲峰值的位置,确定该光信号的色散值。 <BR>因此,处理器320可以计算目标光信号的信号功率(光功率)的自相关函数。在本发明实施例中,目标光信号为需要估算色散的光信号。在本发明实施例中,可以对该目标光信号进行以下处理,以单偏振系统为例,在相干接收机中,目标光信号和本振光在Hybrid中混合,通过平衡接收机后转化为两路电流(同向电流和正交电流),经过两路ADC采样后,转换为数字信号,得到I信号和Q信号,从而,可以根据该I信号和Q信号来确定目标光信号的信号功率的自相关函数。以下,作为示例而非限定,以该通过该I信号和Q信号确定目标光信号的信号功率的自相关函数为例进行说明。。 <BR>可选地,在本发明实施例中,在该目标光信号为非归零码格式信号的情况下,该处理器320通过该总线310,调用该存储器330中存储的程序,以消除所述目标光信号的零频干扰。 <BR>具体地说,零频干扰是指,在该目标光信号为非归零(NRZ,No&nbsp;ReturnZero)码格式的情况下,如果直接获取该格式的信号的自相关函数,则获取的自相关函数的零频位置会出现一个较大的峰值,从而对后述搜索脉冲峰值造成影响,因此,作为示例而非限定,在本发明实施例中,可以对I信号和Q信号进行延迟处理和相减处理,以实现对该目标光信号的高通滤波,其中,该延迟处理使用的迟延大小可以为例如,一个符号周期。从而,能够避免获取的自相关函数在零频位置出现较大的峰值。 <BR>可选地,在本发明实施例中,该处理器320通过该总线310,调用该存储器330中存储的程序,以用于对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 <BR>具体地说,作为示例而非限定,处理器320可以对I信号和Q信号(具体地说,是与该I信号相对应的I信号序列以及与该Q信号相对应的Q信号序列)进行求平方和处理,然后用经过该求平方和处理后的信号减去信号均值,其中,作为示例而非限定,该信号均值可以通过以下方法确定,即,可以获取与该I信号相对应的I信号序列,以及与该Q信号相对应的Q信号序列,其后,对I信号序列和Q信号序列求和,然后除以I信号序列和Q信号序列的长度,从而确定该信号均值。通过以上处理,能够实现对该目标光信号的隔直处理,从而能够去掉直流分量产生的无用的自相关函数部分。 <BR>其后,处理器320可以对经过上述处理的目标光信号进行傅里叶变换,再做幅度平方,最后进行逆傅里叶变换,得到该目标光信号的信号功率的自相关函数,进而可以根据该自相关函数,生成以该自相关函数的各函数值为元素的自相关函数序列(以下,为了便于说明,记作R1[n]) <BR>处理器320可以从该R1[n]中搜索脉冲峰值的位置,即,确定该脉冲峰值在R1[n]中的序号。 <BR>可选地,在本发明实施例中处理器320可以搜索R1[n]中数值最大的元素(记作元素a),并将该元素a的位置确定为脉冲峰值的位置,并可以记录该元素a的序号,以在后述色散值确定单元230确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说元素a在R1[n]中的序号(记作n1),确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,该处理器320通过该总线310,调用该存储器330中存储的程序,以具体用于根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括与该脉冲峰值相对应的元素; <BR>根据该第一自相关函数序列,确定该脉冲位置。 <BR>具体地说,在本发明实施例中,处理器320可以从该R1[n]中选取一个区间作为搜索范围(即,第一自相关函数序列,以下,为了便于说明,记作R2[n]),并且,在本发明实施例中,处理器320可以将R1[n]作为R2[n],也可以将{R1[N1]:R1[N2]}(即,R1[n]中从序号为N1的元素到序号为N2的元素)作为R2[n],其中,N1表示R1[n]中选取的开始位置(序号),N2表示R1[n]中选取的结束位置(序号)。 <BR>可选地,在本发明实施例中,处理器320具体用于根据用于获取该目标光信号的相干接收机的电带宽和用于获取该目标光信号的高速模数转换器的采样率,从该自相关函数序列中确定选取该第一自相关函数序列的开始位置。 <BR>具体地说,处理器320可以根据式1确定N1。 <BR>可选地,在本发明实施例中,该处理器320通过该总线310,调用该存储器330中存储的程序,以具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理; <BR>用于根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 <BR>具体地说,处理器320在确定R2[n]之后,可以对R2[n]进行滑动窗口平均处理(第一滑动平均处理),并且,作为示例而非限定,该滑动窗口平均处理的窗口宽度(记作width1)可以为例如,3。在本发明实施例中,滑动窗口平均处理的方法与现有技术相同,这里,为了避免赘述,省略其说明。在对R2[n]进行滑动窗口平均处理后(以下,为了便于理解,将处理后生成的序列记作R3[n]),可以搜索R3[n]中数值最大的元素(记作元素b),并将该元素b的位置确定为脉冲峰值的位置,并可以记录该元素b在R3[n]中的序号(记作n2),从而,可以根据式2确定n1。 <BR>通过对该R2[n](第一自相关函数序列)进行滑动窗口平均处理,能够减少自相关函数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。以上列举的作为width1的数值仅为本发明的一个示例性说明,并不限定于此。 <BR>可选地,在本发明实施例中,在确定R2[n]之后,处理器320也可以直接搜索R2[n]中数值最大的元素(记作元素c),并将该元素c的位置确定为脉冲峰值的位置,并可以记录该元素c在R2[n]中的序号(记作n3),从而,可以根据式3确定n1。 <BR>从而,在后述确定色散值时,处理器320可以根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,该处理器320通过该总线310,调用该存储器330中存储的程序,以具体用于根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号; <BR>用于根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列; <BR>用于根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理; <BR>用于根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 <BR>在本发明实施例中,处理器320例如,可以将自相关函数中与产生的色散相对应的脉冲峰值的FWHM的一半,即,将脉冲峰值(最大值)的二分之一作为门限(后述第一门限),寻找脉冲宽度,再将脉冲宽度的例如,二分之一作为第二窗口,具体地说,在确定n2后,可以从R2[n]中,在以n2为中心的预设范围(即,第一范围,例如,可以为5个单位,这里“单位”可以是单位序号)内,搜索数值最大的元素(记作元素d),并根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之一作为该第一门限,其后,以元素d为中心,向左侧搜索第一个数值小于该第一门限的元素(记作元素e),并记录该元素e的序号(记作t1);向右侧搜索第一个数值小于该第一门限的元素(记作元素f),并记录该元素f的序号(记作t2),可以根据式4,确定滑动窗口宽度(第二窗口的一例,记作width2)。 <BR>处理器320在确定了width2(第二窗口)后,可以从该R2[n]中选取一个区间作为搜索范围(即,第二自相关函数序列,以下,为了便于说明,记作R4[n]),在本发明实施例中,可以将R2[n]作为R4[n],也可以将{R2[N3]:R2[N4]}(即,R2[n]中从序号为N3的元素到序号为N4的元素)作为R4[n],其中,N3表示R2[n]中选取的开始位置(序号),N4表示R2[n]中选取的结束位置(序号)。具体地说,处理器320可以根据式5确定N3,并可以根据式6确定N4。 <BR>处理器320在确定了R4[n]后,可以根据width2,对该R4[n]进行滑动平均处理(第二滑动平均处理),以下,为了便于理解,将处理后生成的序列记作R5[n],可以搜索R5[n]中数值最大的元素(记作元素g),并将该元素g的位置确定为脉冲峰值的位置,并可以记录该元素g在R5[n]中的序号(记作n4),由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,处理器320可以根据式7确定n1。 <BR>从而,处理器320在后述确定色散值时,可以根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>通过对该R2[n](第一自相关函数序列)进行滑动窗口平均处理,能够减少自相关函数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。 <BR>可选地,在本发明实施例中,该处理器320通过该总线310,调用该存储器330中存储的程序,以具体用于根据第一序号和第一预设值,从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应; <BR>用于根据该第一范围中数值最大的第一元素和第二预设值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相对应; <BR>用于根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 <BR>具体地说,如果PMD较大,则可能会出现一个自相关函数中的脉冲(峰值)分裂为两个脉冲(峰值)的情况,此情况下,如果仅以脉冲峰值(最大值)的位置为中心寻找脉冲宽度,则寻找到的脉冲宽度值会产生偏差,因此,在本发明实施例中,处理器320在确定n2后,还可以从R2[n]中,在以n2为中心的预设范围(即,第一范围,与第一预设值相对应,例如,可以为5个单位)内,搜索数值最大的元素(即,元素d)。 <BR>其后,在以元素d在R2[n]中的序号为中心的的预设范围(即,第二范围,与第二预设值相对应,例如,可以为10个单位,这里为单位序号),搜索数值第二大的元素(记作元素h)。 <BR>其后,确定元素d与元素h在R2[n]中的位置关系,这里,为了便于理解,以元素d在元素h的右侧为例(即,元素d在R2[n]中的序号大于元素h在R2[n]中的序号)进行说明。 <BR>其后,根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之一作为该第一门限。 <BR>其后,以元素h为中心,向左侧搜索第一个数值小于该第一门限的元素(记作元素i,并记录该元素i的序号(记作t3);以元素d为中心,向右侧搜索第一个数值小于该第一门限的元素(记作元素j),并记录该元素j的序号(记作t4),可以根据式8,确定一个中间参量(记作width3)。 <BR>处理器320在确定了width3后,可以从该R2[n]中选取一个区间作为搜索范围(即,第二自相关函数序列,以下,为了便于说明,记作R6[n]),在本发明实施例中,可以将R2[n]作为R6[n],也可以将{R2[N5]:R2[N6]}(即,R2[n]中从序号为N5的元素到序号为N6的元素)作为R6[n],其中,N5表示R2[n]中选取的开始位置(序号),N6表示R2[n]中选取的结束位置(序号)。具体地说,处理器320可以根据式9确定N5,并可以根据式10确定N6。 <BR>脉冲位置搜索单元220在确定了R6[n]后,可以根据式11,确定对R6[n]进行滑动平均处理时使用的窗口宽度(第二窗口的另一例,记作width4)。 <BR>其后,处理器320对该R6[n]进行滑动平均处理(第二滑动平均处理),以下,为了便于理解,将处理后生成的序列记作R7[n],可以从该R7[n]中选取一个区间作为搜索范围(以下,为了便于说明,记作R8[n]),在本发明实施例中,可以将R7[n]作为R8[n],也可以将{R7[N7]:R7[N8]}(即,R7[n]中从序号为N7的元素到序号为N8的元素)作为R8[n],其中,N7表示R7[n]中选取的开始位置(序号),N8表示R7[n]中选取的结束位置(序号)。具体地说,脉冲位置搜索单元220可以根据式12确定N7,并可以根据式13确定N8。 <BR>处理器320在确定R8[n]后,可以搜索R8[n]中数值最大的元素(记作元素k),并记录该元素k在R8[n]中的序号(记作n5)。 <BR>在本发明实施例中,处理器320可以将该元素k的位置确定为脉冲峰值的位置,由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,处理器320可以根据式14确定n1。 <BR>从而,在后述处理器320确定色散值时,根据该脉冲峰值在R1[n]中的位置,或者说n1的值,确定色散值(随后进行详细说明)。 <BR>可选地,在本发明实施例中,处理器320在确定n5后,还可以利用R8[n5‑1]、R8[n5]、R8[n5+1]三点做抛物线拟合,从而,可以根据式15,确定该抛物线顶点在R8[n]中的位置(序号,记作n6)。 <BR>其后,处理器320可以根据n6和式16,确定该抛物线顶点在R2[n]中的位置(序号)。 <BR>由于自相关函数第一位对应τ=0,R[n]对应延迟为(n‑1)·SamplingInterval(修正一位),其中,SamplingInterval表示采样间隔,从而,可以根据式17确定n1。 <BR>从而,处理器320可以根据所确定的自相关函数中的脉冲峰值出现位置(即,n1),确定目标光信号的色散值。 <BR>在本发明实施例中,处理器320根据离线数据,确定指示脉冲峰值位置与色散值的对应关系的指示信息, <BR>根据该脉冲位置和该指示信息,确定该光信号的色散值。 <BR>具体地说,作为示例而非限定,例如,处理器320可以分别对多个光信号(离线数据的一例)进行色散值测试等试验,记录各光信号的色散值和它们各自对应的自相关函数中脉冲峰值出现的位置,作为训练数据,从而,可以采用例如期望最大化算法,通过逐步改进公式的参数,使参数和训练数据的似然概率逐渐增大,最后终止于一个极大点。直观地理解,期望最大化算法也可被看作为一个逐次逼近算法:事先并不知道模型的参数,可以随机的选择一套参数或者事先粗略地给定某个初始参数λ0,确定出对应于这组参数的最可能的状态,计算每个训练样本的可能结果的概率,在当前的状态下再由样本对参数修正,重新估计参数λ,并在新的参数下重新确定模型的状态,这样,通过多次的迭代,循环直至某个收敛条件满足为止,就可以使得模型的参数逐渐逼近真实参数。因此,该步骤可以以离线的方式自动进行。将离线计算得到的上述各种脉冲峰值位置与各自对应的色散值预先存储起来,例如以表的形式存储起来。应理解,以上列举的各公式以及各公式中的常数的训练方法仅是本发明的一个实施例,本发明并不限定于此。 <BR>作为示例而非限定,处理器320可以根据式18确定色散值。 <BR>应理解,以上列举的各公式、参数以及数值仅为本发明的示例性说明,本发明并不限定于此,其他能够用于确定光信号的自相关函数的脉冲峰值出现位置的各公式、参数以及数值均落入本发明的保护范围内。 <BR>根据本发明实施例的色散估计的系统可对应于本发明实施例的色散估计的方法的执行主体,并且,该色散估计的装置中的各单元和上述其他操作和/或功能分别为了实现图1中的方法100的相应流程,为了简洁,在此不再赘述。 <BR>根据本发明实施例的色散估计的系统,通过获取目标光信号的光信号功率自相关函数,并确定该光信号功率自相关函数出现脉冲峰值的位置,能够在无需扫描所有可能的色散值的情况下,根据该出现脉冲峰值的位置确定目标光信号的色散值,到达快速估计信号的色散值的目的。应理解,在本发明的各种实施例中,上述各过程的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。 <BR>本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。 <BR>所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统、装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。 <BR>在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。 <BR>所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。 <BR>另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。 <BR>所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read‑Only&nbsp;Memory)、随机存取存储器(RAM,Random&nbsp;Access&nbsp;Memory)、磁碟或者光盘等各种可以存储程序代码的介质。 <BR>以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。</p></div> </div> </div> </div> <div class="tempdiv cssnone" style="line-height:0px;height:0px; overflow:hidden;"> </div> <div id="page"> <div class="page"><img src='https://img.zhuanlichaxun.net/fileroot2/2018-8/27/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb95/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb951.gif' alt="色散估计的方法、装置和系统.pdf_第1页" width='100%'/></div><div class="pageSize">第1页 / 共29页</div> <div class="page"><img src='https://img.zhuanlichaxun.net/fileroot2/2018-8/27/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb95/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb952.gif' alt="色散估计的方法、装置和系统.pdf_第2页" width='100%'/></div><div class="pageSize">第2页 / 共29页</div> <div class="page"><img src='https://img.zhuanlichaxun.net/fileroot2/2018-8/27/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb95/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb953.gif' alt="色散估计的方法、装置和系统.pdf_第3页" width='100%'/></div><div class="pageSize">第3页 / 共29页</div> </div> <div id="pageMore" class="btnmore" onclick="ShowSvg();">点击查看更多>></div> <div style="margin-top:20px; line-height:0px; height:0px; overflow:hidden;"> <div style=" font-size: 16px; background-color:#e5f0f7; font-weight: bold; text-indent:10px; line-height: 40px; height:40px; padding-bottom: 0px; margin-bottom:10px;">资源描述</div> <div class="detail-article prolistshowimg"> <p>《色散估计的方法、装置和系统.pdf》由会员分享,可在线阅读,更多相关《色散估计的方法、装置和系统.pdf(29页珍藏版)》请在专利查询网上搜索。</p> <p >1、(10)申请公布号 CN 103004110 A (43)申请公布日 2013.03.27 C N 1 0 3 0 0 4 1 1 0 A *CN103004110A* (21)申请号 201280001461.X (22)申请日 2012.08.03 H04B 10/077(2013.01) H04B 10/2513(2013.01) (71)申请人华为技术有限公司 地址 518129 广东省深圳市龙岗区坂田华为 总部办公楼 申请人香港理工大学 (72)发明人隋琪 刘伯涛 吕超 李良川 (74)专利代理机构北京龙双利达知识产权代理 有限公司 11329 代理人毛威 张亮 (54) 发明名称 。</p> <p >2、色散估计的方法、装置和系统 (57) 摘要 本发明实施例提供了一种色散估计的方法、 装置和系统,该方法包括:计算目标光信号的信 号功率的自相关函数,并获取该自相关函数的自 相关函数序列;根据该自相关函数序列,确定脉 冲位置,该脉冲位置与该自相关函数的脉冲峰值 的位置相对应;根据该脉冲位置,确定该目标光 信号的色散值。根据本发明实施例的色散估计的 方法、装置和系统,通过获取目标光信号的功率自 相关函数,并确定该功率自相关函数出现脉冲峰 值的位置,能够在无需扫描所有可能的色散值的 情况下,根据该出现脉冲峰值的位置确定目标光 信号的色散值,到达快速估计信号的色散值的目 的。 (85)PCT申请进入国。</p> <p >3、家阶段日 2012.12.10 (86)PCT申请的申请数据 PCT/CN2012/079619 2012.08.03 (51)Int.Cl. 权利要求书4页 说明书22页 附图2页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 4 页 说明书 22 页 附图 2 页 1/4页 2 1.一种色散估计的方法,其特征在于,所述方法包括: 计算目标光信号的信号功率的自相关函数,并获取所述自相关函数的自相关函数序 列; 根据所述自相关函数序列,确定脉冲位置,所述脉冲位置与所述自相关函数的脉冲峰 值的位置相对应; 根据所述脉冲位置,确定所述目标光信号的色散值。 2.根据权利要求。</p> <p >4、1所述的方法,其特征在于,在所述目标光信号为非归零码格式信号 的情况下,在根据目标光信号的信号功率的自相关函数,确定自相关函数序列前,所述方法 还包括: 对所述目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 3.根据权利要求1或2所述的方法,其特征在于,在根据目标光信号的信号功率的自相 关函数,确定自相关函数序列前,所述方法还包括: 对所述目标光信号进行第二处理,以消除所述目标光信号的直流分量干扰。 4.根据权利要求1至3中任一项所述的方法,其特征在于,所述根据所述自相关函数序 列,确定脉冲位置,包括: 根据所述自相关函数序列,确定第一自相关函数序列,所述第一自相关函数序列包括 所述。</p> <p >5、脉冲峰值; 根据所述第一自相关函数序列,确定所述脉冲位置。 5.根据权利要求4所述的方法,其特征在于,所述根据所述第一自相关函数序列,确定 所述脉冲位置,包括: 根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理; 根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置。 6.根据权利要求5所述的方法,其特征在于,所述根据所述第一滑动平均处理后的第 一自相关函数序列,确定所述脉冲位置,包括: 根据第一序号和第一预设值,确定第二窗口,所述第一序号是所述第一滑动平均处理 后的第一自相关函数序列中数值最大的元素的序号; 根据所述第一自相关函数序列和所述第二窗口,确定第二自相关函数。</p> <p >6、序列; 根据所述第二窗口,对所述第二自相关函数序列进行第二滑动平均处理; 根据所述第二滑动平均处理后的第二自相关函数序列,确定所述脉冲位置。 7.根据权利要求6所述的方法,其特征在于,所述根据第一序号和第一预设值,确定第 二窗口,包括: 根据第一序号和第一预设值,从所述第一自相关函数序列中确定第一范围,所述第一 范围以所述第一序号为中心,所述第一范围的大小与所述第一预设值相对应; 根据所述第一范围中数值最大的第一元素和第二预设值,从所述第一自相关函数序列 中确定第二范围,所述第二范围以所述第一元素的序号为中心,所述第二范围的大小与所 述第二预设值相对应; 根据所述第一元素和所述第二范围中数值第。</p> <p >7、二大的第二元素,确定所述第二窗口。 8.根据权利要求1至7中任一项所述的方法,其特征在于,所述根据所述脉冲位置,确 定所述目标光信号的色散值,包括: 权 利 要 求 书CN 103004110 A 2/4页 3 根据以下公式确定所述目标光信号的色散值, 其中,CD表示色散值,n表示所述脉冲位置在所述自相关函数序列中的序号,T 0 表示所 述目标光信号的初始发射机光脉冲宽度,c表示真空光速,表示所述目标光信号的波长, SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光信号的符号速率。 9.一种色散估计的装置,其特征在于,所述装置包括: 功率自相关处理单元,用于计算。</p> <p >8、目标光信号的信号功率的自相关函数,并获取所述自 相关函数的自相关函数序列,以及向脉冲位置搜索单元传输指示所述自相关函数序列的信 息; 脉冲位置搜索单元,用于从所述功率自相关处理单元获取指示所述自相关函数序列的 信息,并根据所述自相关函数序列,确定脉冲位置,所述脉冲位置与所述自相关函数的脉冲 峰值的位置相对应,以及向色散值确定单元传输指示所述脉冲位置的信息; 色散值确定单元,用于从所述脉冲位置搜索单元获取指示所述脉冲位置的信息,并根 据所述脉冲位置,确定所述目标光信号的色散值。 10.根据权利要求9所述的装置,其特征在于,在所述目标光信号为非归零码格式信号 的情况下,所述装置还包括: 第一干扰处。</p> <p >9、理单元,用于对所述目标光信号进行第一处理,以消除所述目标光信号的 零频干扰。 11.根据权利要求9或10所述的装置,其特征在于,所述装置还包括: 第二干扰处理单元,用于对所述目标光信号进行第二处理,以消除所述目标光信号的 直流分量干扰。 12.根据权利要求9至11中任一项所述的装置,其特征在于,所述脉冲位置搜索单元具 体用于根据所述自相关函数序列,确定第一自相关函数序列,所述第一自相关函数序列包 括所述脉冲峰值; 用于根据所述第一自相关函数序列,确定所述脉冲位置。 13.根据权利要求12所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据 预设的第一窗口,对第一自相关函数序列进行第一滑动平。</p> <p >10、均处理; 用于根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置。 14.根据权利要求13所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据 第一序号和第一预设值,确定第二窗口,所述第一序号是所述第一滑动平均处理后的第一 自相关函数序列中数值最大的元素的序号; 用于根据所述第一自相关函数序列和所述第二窗口,确定第二自相关函数序列; 用于根据所述第二窗口,对所述第二自相关函数序列进行第二滑动平均处理; 用于根据所述第二滑动平均处理后的第二自相关函数序列,确定所述脉冲位置。 15.根据权利要求14所述的装置,其特征在于,所述脉冲位置搜索单元具体用于根据 第一序号和第一预设值,从。</p> <p >11、所述第一自相关函数序列中确定第一范围,所述第一范围以所 述第一序号为中心,所述第一范围的大小与所述第一预设值相对应; 权 利 要 求 书CN 103004110 A 3/4页 4 用于根据所述第一范围中数值最大的第一元素和第二预设值,从所述第一自相关函数 序列中确定第二范围,所述第二范围以所述第一元素的序号为中心,所述第二范围的大小 与所述第二预设值相对应; 用于根据所述第一元素和所述第二范围中数值第二大的第二元素,确定所述第二窗 口。 16.根据权利要求9至15中任一项所述的装置,其特征在于,所述色散值确定单元具体 用于根据以下公式确定所述目标光信号的色散值, 其中,CD表示色散值,n表示所。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>12、述脉冲位置在所述自相关函数序列中的序号,T 0 表示所 述目标光信号的初始发射机光脉冲宽度,c表示真空光速,表示所述目标光信号的波长, SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光信号的符号速率。 17.一种色散估计的系统,其特征在于,所述系统包括: 总线; 与所述总线相连的处理器; 与所述总线相连的存储器; 其中,所述处理器通过所述总线,调用所述存储器中存储的程序,以用于计算目标光信 号的信号功率的自相关函数,并获取所述自相关函数的自相关函数序列; 根据所述自相关函数序列,确定脉冲位置,所述脉冲位置与所述自相关函数的脉冲峰 值的出现位置相对应; 根据所述。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>13、脉冲位置,确定所述目标光信号的色散值。 18.根据权利要求17所述的系统,其特征在于,在所述目标光信号为非归零码格式信 号的情况下,所述处理器通过所述总线,调用所述存储器中存储的程序,以用于对所述目标 光信号进行第一处理,以消除所述目标光信号的零频干扰。 19.根据权利要求17或18所述的系统,其特征在于,所述处理器通过所述总线,调用所 述存储器中存储的程序,以用于对所述目标光信号进行第二处理,以消除所述目标光信号 的直流分量干扰。 20.根据权利要求17至19中任一项所述的系统,其特征在于,所述处理器通过所述总 线,调用所述存储器中存储的程序,以具体用于根据所述自相关函数序列,确定第一自相关。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>14、 函数序列,所述第一自相关函数序列包括所述脉冲峰值; 根据所述第一自相关函数序列,确定所述脉冲位置。 21.根据权利要求20所述的系统,其特征在于,所述处理器通过所述总线,调用所述存 储器中存储的程序,以具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑 动平均处理; 用于根据所述第一滑动平均处理后的第一自相关函数序列,确定所述脉冲位置。 22.根据权利要求21所述的系统,其特征在于,所述处理器通过所述总线,调用所述存 储器中存储的程序,以具体用于根据第一序号和第一预设值,确定第二窗口,所述第一序号 是所述第一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号; 权 利 要 求 。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>15、书CN 103004110 A 4/4页 5 用于根据所述第一自相关函数序列和所述第二窗口,确定第二自相关函数序列; 用于根据所述第二窗口,对所述第二自相关函数序列进行第二滑动平均处理; 用于根据所述第二滑动平均处理后的第二自相关函数序列,确定所述脉冲位置。 23.根据权利要求21所述的系统,其特征在于,所述处理器通过所述总线,调用所述存 储器中存储的程序,以具体用于根据第一序号和第一预设值,从所述第一自相关函数序列 中确定第一范围,所述第一范围以所述第一序号为中心,所述第一范围的大小与所述第一 预设值相对应; 用于根据所述第一范围中数值最大的第一元素和第二预设值,从所述第一自相关函数 序列中。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>16、确定第二范围,所述第二范围以所述第一元素的序号为中心,所述第二范围的大小 与所述第二预设值相对应; 用于根据所述第一元素和所述第二范围中数值第二大的第二元素,确定所述第二窗 口。 24.根据权利要求17至23中任一项所述的系统,其特征在于,所述处理器通过所述总 线,调用所述存储器中存储的程序,以具体用于根据以下公式确定所述目标光信号的色散 值, 其中,CD表示色散值,n表示所述色散出现位置在所述自相关函数序列中的序号,T 0 表 示所述目标光信号的初始发射机光脉冲宽度,c表示真空光速,表示所述目标光信号的 波长,SamplingRate表示模数转换器的采样率,SymbolRate表示所述目标光。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>17、信号的符号速 率。 权 利 要 求 书CN 103004110 A 1/22页 6 色散估计的方法、 装置和系统 技术领域 0001 本发明涉及通信领域,并且更具体地,涉及色散估计的方法、装置和系统。 背景技术 0002 目前,已知一种色散估计的方法,即,对于未知色散的信号,以一定的步长间隔,扫 描所有可能的色散值,并对于每个色散值通过设计代价函数(也可以称为算法设计)计算出 对应的函数值,最后通过比较所有计算出的函数值(最大值或最小值),找到实际的色散值, 实现色散估计。 0003 但是,由于该方法需要扫描所有的色散值以进行计算,因此需要消耗大量的时间, 不能快速估计信号的色散值。 发明内容。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>18、 0004 本发明实施例提供一种,能够快速估计信号的色散值。 0005 第一方面,提供了一种色散估计的方法,该方法包括:计算目标光信号的信号功率 的自相关函数,并获取该自相关函数的自相关函数序列;根据该自相关函数序列,确定脉冲 位置,该脉冲位置与该自相关函数的脉冲峰值的位置相对应;根据该脉冲位置,确定该目标 光信号的色散值。 0006 在一种可能的实施方式中,在该目标光信号为非归零码格式信号的情况下,在根 据目标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括:对该目标 光信号进行第一处理,以消除所述目标光信号的零频干扰。 0007 结合第一方面和第一种可能的实施方式,在第二种。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>19、可能的实施方式中,在根据目 标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括:对该目标光信 号进行第二处理,以消除该目标光信号的直流分量干扰。 0008 结合第一方面、第一种可能的实施方式和第二种可能的实施方式,在第三种可能 的实施方式中,该根据该自相关函数序列,确定脉冲位置,包括:根据该自相关函数序列,确 定第一自相关函数序列,该第一自相关函数序列包括该脉冲峰值;根据该第一自相关函数 序列,确定该脉冲位置。 0009 结合第一方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的 实施方式,在第四种可能的实施方式中,该根据该第一自相关函数序列,确定该脉冲位置, 包括。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>20、:根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;根据该第一 滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 0010 结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式和第四种可能的实施方式,在第五种可能的实施方式中,该根据该第一滑动平均处 理后的第一自相关函数序列,确定该脉冲位置,包括:根据第一序号和第一预设值,确定第 二窗口,该第一序号是该第一滑动平均处理后的第一自相关函数序列中数值最大的元素的 序号;根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列;根据该第二 说 明 书CN 103004110 A 2/22页 7 窗口,对。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>21、该第二自相关函数序列进行第二滑动平均处理;根据该第二滑动平均处理后的第 二自相关函数序列,确定该脉冲位置。 0011 结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式、第四种可能的实施方式和第五种可能的实施方式,在第六种可能的实施方式中,该 根据第一序号和第一预设值,确定第二窗口,包括:根据第一序号和第一预设值,从该第一 自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围的大小与 该第一预设值相对应;根据该第一范围中数值最大的第一元素和第二预设值,从该第一自 相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大 小与。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>22、该第二预设值相对应;根据该第一元素和该第二范围中数值第二大的第二元素,确定 该第二窗口。 0012 结合第一方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式,在第七 种可能的实施方式中,该根据该脉冲位置,确定该目标光信号的色散值,包括:根据以下公 式确定该目标光信号的色散值, 0013 0014 0015 其中,CD表示色散值,n表示该脉冲位置在该自相关函数序列中的序号,T 0 表示 该目标光信号的初始发射机光脉冲宽度,c表示真空光速,表示该目标光信号的波长, SamplingRate表示模数转换器的采样。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>23、率,SymbolRate表示该目标光信号的符号速率。 0016 第二方面,提供了一种色散估计的装置,该装置包括:功率自相关处理单元,用于 计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列,以及 向脉冲位置搜索单元传输指示该自相关函数序列的信息;脉冲位置搜索单元,用于从该功 率自相关处理单元获取指示该自相关函数序列的信息,并根据该自相关函数序列,确定脉 冲位置,该脉冲位置与该自相关函数的脉冲峰值的位置相对应,以及向色散值确定单元传 输指示该脉冲位置的信息;色散值确定单元,用于从该脉冲位置搜索单元获取指示该脉冲 位置的信息,并根据该脉冲位置,确定该目标光信号的色散值。 00。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>24、17 在一种可能的实施方式中,在该目标光信号为非归零码格式信号的情况下,该装 置还包括:第一干扰处理单元,用于对该目标光信号进行第一处理,以消除所述目标光信号 的零频干扰。 0018 结合第二方面和第一种可能的实施方式,在第二种可能的实施方式中,该装置还 包括:第二干扰处理单元,用于对该目标光信号进行第二处理,以消除该目标光信号的直流 分量干扰。 0019 结合第二方面、第一种可能的实施方式和第二种可能的实施方式,在第三种可能 的实施方式中,该脉冲位置搜索单元具体用于根据该自相关函数序列,确定第一自相关函 数序列,该第一自相关函数序列包括该脉冲峰值;用于根据该第一自相关函数序列,确定该 脉冲位。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>25、置。 0020 结合第二方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的 实施方式,在第四种可能的实施方式中,该脉冲位置搜索单元具体用于根据预设的第一窗 说 明 书CN 103004110 A 3/22页 8 口,对第一自相关函数序列进行第一滑动平均处理;用于根据该第一滑动平均处理后的第 一自相关函数序列,确定该脉冲位置。 0021 结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式和第四种可能的实施方式,在第五种可能的实施方式中,该脉冲位置搜索单元具体 用于根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理后的 第一自相关函数序。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>26、列中数值最大的元素的序号;用于根据该第一自相关函数序列和该第二 窗口,确定第二自相关函数序列;用于根据该第二窗口,对该第二自相关函数序列进行第二 滑动平均处理;用于根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 0022 结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式、第四种可能的实施方式和第五种可能的实施方式,在第六种可能的实施方式中,该 脉冲位置搜索单元具体用于根据第一序号和第一预设值,从该第一自相关函数序列中确定 第一范围,该第一范围以该第一序号为中心,该第一范围的大小与该第一预设值相对应;用 于根据该第一范围中数值最大的第一元素和第二预。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>27、设值,从该第一自相关函数序列中确定 第二范围,该第二范围以该第一元素的序号为中心,该第二范围的大小与该第二预设值相 对应;用于根据该第一元素和该第二范围中数值第二大的第二元素,确定该第二窗口。 0023 结合第二方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式,在第七 种可能的实施方式中,该色散值确定单元具体用于根据以下公式确定该目标光信号的色散 值, 0024 0025 0026 其中,CD表示色散值,n表示该脉冲位置在该自相关函数序列中的序号,T 0 表示 该目标光信号的初始发射机光脉冲宽度,c表示真。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>28、空光速,表示该目标光信号的波长, SamplingRate表示模数转换器的采样率,SymbolRate表示该目标光信号的符号速率。 0027 第三方面,提供了一种色散估计的系统,该系统包括:总线;与该总线相连的处理 器;与该总线相连的存储器;其中,该处理器通过该总线,调用该存储器中存储的程序,以 用于计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关函数序列; 根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲峰值的出现位 置相对应;根据该脉冲位置,确定该目标光信号的色散值。 0028 在一种可能的实施方式中,在该目标光信号为非归零码格式信号的情况下,该处 理器通。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>29、过该总线,调用该存储器中存储的程序,以用于对该目标光信号进行第一处理,以消 除所述目标光信号的零频干扰。 0029 结合第三方面和第一种可能的实施方式,在第二种可能的实施方式中,该处理器 通过该总线,调用该存储器中存储的程序,以用于对该目标光信号进行第二处理,以消除该 目标光信号的直流分量干扰。 0030 结合第三方面、第一种可能的实施方式和第二种可能的实施方式,在第三种可能 的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据该自相 关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括该脉冲峰值;根据该第 说 明 书CN 103004110 A 4/22页 9 。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>30、一自相关函数序列,确定该脉冲位置。 0031 结合第三方面、第一种可能的实施方式、第二种可能的实施方式和第三种可能的 实施方式,在第四种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程 序,以具体用于根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理;用于 根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 0032 结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式和第四种可能的实施方式,在第五种可能的实施方式中,该处理器通过该总线,调用 该存储器中存储的程序,以具体用于根据第一序号和第一预设值,确定第二窗口,该第一序 号是该第。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>31、一滑动平均处理后的第一自相关函数序列中数值最大的元素的序号;用于根据该 第一自相关函数序列和该第二窗口,确定第二自相关函数序列;用于根据该第二窗口,对该 第二自相关函数序列进行第二滑动平均处理;用于根据该第二滑动平均处理后的第二自相 关函数序列,确定该脉冲位置。 0033 结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式、第四种可能的实施方式和第五种可能的实施方式,在第六种可能的实施方式中,该 处理器通过该总线,调用该存储器中存储的程序,以具体用于根据第一序号和第一预设值, 从该第一自相关函数序列中确定第一范围,该第一范围以该第一序号为中心,该第一范围 的大小与。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>32、该第一预设值相对应;用于根据该第一范围中数值最大的第一元素和第二预设 值,从该第一自相关函数序列中确定第二范围,该第二范围以该第一元素的序号为中心,该 第二范围的大小与该第二预设值相对应;用于根据该第一元素和该第二范围中数值第二大 的第二元素,确定该第二窗口。 0034 结合第三方面、第一种可能的实施方式、第二种可能的实施方式、第三种可能的实 施方式、第四种可能的实施方式、第五种可能的实施方式和第六种可能的实施方式,在第七 种可能的实施方式中,该处理器通过该总线,调用该存储器中存储的程序,以具体用于根据 以下公式确定该目标光信号的色散值, 0035 0036 0037 其中,CD表示色散值,n。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>33、表示该色散出现位置在该自相关函数序列中的序号,T 0 表 示该目标光信号的初始发射机光脉冲宽度,c表示真空光速,表示该目标光信号的波长, SamplingRate表示模数转换器的采样率,SymbolRate表示该目标光信号的符号速率。 0038 根据本发明实施例的色散估计的方法、装置和系统,通过获取目标光信号的功率 自相关函数,并确定该功率自相关函数出现脉冲峰值的位置,能够在无需扫描所有可能的 色散值的情况下,根据该出现脉冲峰值的位置确定目标光信号的色散值,到达快速估计信 号的色散值的目的。 附图说明 0039 为了更清楚地说明本发明实施例的技术方案,下面将对实施例或现有技术描述中 所需要使用。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>34、的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实 施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图 说 明 书CN 103004110 A 5/22页 10 获得其他的附图。 0040 图1是根据本发明实施例的色散估计的方法的示意性流程图。 0041 图2是根据本发明实施例的色散估计的装置的示意性框图。 0042 图3是根据本发明实施例的色散估计的方法的示意性框图。 具体实施方式 0043 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完 整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发 。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>35、明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施 例,都属于本发明保护的范围。 0044 图1示出了根据本发明实施例的色散估计的方法100的示意性流程图,如图1所 示,该方法100包括: 0045 S110,计算目标光信号的信号功率的自相关函数,并获取该自相关函数的自相关 函数序列; 0046 S120,根据该自相关函数序列,确定脉冲位置,该脉冲位置与该自相关函数的脉冲 峰值的位置相对应; 0047 S130,根据该脉冲位置,确定该目标光信号的色散值。 0048 在本发明实施例中,通过试验(例如,获取多个光信号的功率自相关函数,并获取 该多个光信号的色散值)发现,。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>36、在光信号(例如,单偏振的相移键控法(PSK,Phase Shift Keying)信号、偏振分割复用(PDM,PolarizationDivision Multiplexing)的PSK信号或正 交振幅调制(QAM,QuadratureAmplitude Modulation)信号)存在色散时,该光信号的功率 自相关函数会出现一个脉冲峰值,并且,不同的脉冲峰值出现位置与不同的色散值之间存 在映射对应关系,因此,在本发明实施例中,可以根据光信号的信号功率的自相关函数中脉 冲峰值的位置,确定该光信号的色散值。 0049 因此,在S110中,可以计算目标光信号的信号功率(光功率)的自相关函数。在 本。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>37、发明实施例中,目标光信号为需要估算色散的光信号。在本发明实施例中,可以对该 目标光信号进行以下处理,以单偏振系统为例,在相干接收机中,目标光信号和本振光 在混频器(Hybrid)中混合,通过平衡接收机后转化为同向(in-phase)电流和正交电流 (quadrature)两路电流,经过两路高速模数转换器(ADC,Analog to Digital Converter) 采样后,转换为数字信号,得到I信号(同向信号)和Q信号(正交信号),从而,可以根据该I 信号和Q信号来确定目标光信号的信号功率的自相关函数。以下,作为示例而非限定,以该 通过该I信号和Q信号确定目标光信号的信号功率的自相关函数为。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>38、例进行说明。 0050 可选地,在本发明实施例中,在该目标光信号为非归零码格式信号的情况下,在根 据目标光信号的信号功率的自相关函数,确定自相关函数序列前,该方法还包括: 0051 对该目标光信号进行第一处理,以消除所述目标光信号的零频干扰。 0052 具体地说,零频干扰是指,在该目标光信号为非归零(NRZ,No ReturnZero)码格式 的情况下,如果直接获取该格式的信号的自相关函数,则获取的自相关函数的零频位置会 出现一个较大的峰值,从而对后述搜索脉冲峰值造成影响,因此,作为示例而非限定,在本 发明实施例中,可以对I信号和Q信号进行延迟处理和相减处理,以实现对该目标光信号的 说 明 书。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>39、CN 103004110 A 10 6/22页 11 高通滤波,其中,该延迟处理使用的迟延大小可以为例如,一个符号周期。从而,能够避免获 取的自相关函数在零频位置出现较大的峰值。 0053 可选地,在本发明实施例中,在根据目标光信号的信号功率的自相关函数,确定自 相关函数序列前,该方法还包括: 0054 对该目标光信号进行第二处理,以消除该目标光信号的直流分量干扰。 0055 具体地说,作为示例而非限定,可以对I信号和Q信号(具体地说,是与该I信号相 对应的I信号序列以及与该Q信号相对应的Q信号序列)进行求平方和处理,然后用经过该 求平方和处理后的信号减去信号均值,其中,作为示例而非限定,该信。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>40、号均值可以通过以下 方法确定,即,可以获取与该I信号相对应的I信号序列,以及与该Q信号相对应的Q信号 序列,其后,对I信号序列和Q信号序列求和,然后除以I信号序列和Q信号序列的长度,从 而确定该信号均值。通过以上处理,能够实现对该目标光信号的隔直处理,从而能够去掉直 流分量产生的无用的自相关函数部分。 0056 其后,在S110中,可以对经过上述处理的目标光信号进行傅里叶变换,再做幅度 平方,最后进行逆傅里叶变换,得到该目标光信号的信号功率的自相关函数,进而可以根据 该自相关函数,生成以该自相关函数的各函数值为元素的自相关函数序列(以下,为了便于 说明,记作R 1 n)。并且自相关函数以表示为。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>41、: 0057 0058 其中,表示目标光信号的自相关函数,FFT表示快速傅里叶变换,IFFT表示 快速反傅里叶变换,yn表示目标光信号的信号功率(光功率),R 1 n为根据该自相关函数 生成的函数序列(自相关函数序列)。 0059 在S120中,可以从该R 1 n中搜索脉冲峰值的位置(脉冲位置),即,确定该脉冲峰 值在R 1 n中的序号。 0060 可选地,在本发明实施例中,可以搜索R 1 n中数值最大的元素(记作元素a),并将 该元素a的位置确定为脉冲峰值的位置,并可以记录该元素a的序号,以在后述确定色散值 时,根据该脉冲峰值在R 1 n中的位置,或者说元素a在R 1 n中的序号(记作n 1。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>42、 ),确定色散 值(随后进行详细说明)。 0061 可选地,在本发明实施例中,该根据该自相关函数序列,确定脉冲位置,包括: 0062 根据该自相关函数序列,确定第一自相关函数序列,该第一自相关函数序列包括 与该脉冲峰值相对应的元素; 0063 根据该第一自相关函数序列,确定该脉冲位置。 0064 具体地说,在本发明实施例中,可以从该R 1 n中选取一个区间作为搜索范围(即, 第一自相关函数序列,以下,为了便于说明,记作R 2 n),并且,在本发明实施例中,可以将 R 1 n作为R 2 n,也可以将R 1 N 1 :R 1 N 2 (即,R 1 n中从序号为N 1 的元素到序号为N 2 的 元素。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>43、)作为R 2 n,其中,N 1 表示R 1 n中选取的开始位置(序号),N 2 表示R 1 n中选取的结 束位置(序号)。 0065 可选地,在本发明实施例中,该根据该自相关函数序列,确定第一自相关函数序 列,包括: 0066 根据用于获取所述目标光信号的相干接收机的电带宽和用于获取所述目标光信 说 明 书CN 103004110 A 11 7/22页 12 号的高速模数转换器的采样率,从所述自相关函数序列中确定选取所述第一自相关函数序 列的开始位置。 0067 具体地说,可以根据以下式1确定N 1 ,即: 0068 N 1 floor(2SamplingRate/B e ) 式1 0069 。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>44、其中,floor为向下取整数函数,Be表示接收机电带宽,SamplingRate是表示ADC 采样率,通常情况下,Be10GHz,SamplingRate20G Sa/s。 0070 并且,在本发明实施例中,N 2 可以表示能够进行色散估计的最大范围,在本发明实 施例中,N 2 N 1 。 0071 可选地,在本发明实施例中,该根据该第一自相关函数序列,确定该脉冲位置,包 括: 0072 根据预设的第一窗口,对第一自相关函数序列进行第一滑动平均处理; 0073 根据该第一滑动平均处理后的第一自相关函数序列,确定该脉冲位置。 0074 具体地说,在确定R 2 n之后,对R 2 n进行滑动窗口平均。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>45、处理(第一滑动平均处 理),并且,作为示例而非限定,该滑动窗口平均处理的窗口宽度(记作width 1 )可以为例如, 3。在本发明实施例中,滑动窗口平均处理的方法与现有技术相同,这里,为了避免赘述,省 略其说明。在对R 2 n进行滑动窗口平均处理后(以下,为了便于理解,将处理后生成的序列 记作R 3 n),可以搜索R 3 n中数值最大的元素(记作元素b),并将该元素b的位置确定为 脉冲峰值的位置,并可以记录该元素b在R 3 n中的序号(记作n 2 ),从而,可以根据以下式2 确定n 1 ,即 0075 n 1 n 2 N 1 1 式2 0076 通过对该R 2 n(第一自相关函数序列)进行滑动。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>46、窗口平均处理,能够减少自相关函 数中噪声的干扰,使确定的脉冲峰值的位置(脉冲位置)更加准确。以上列举的作为width 1 的数值仅为本发明的一个示例性说明,并不限定于此。 0077 可选地,在本发明实施例中,在确定R 2 n之后,也可以直接搜索R 2 n中数值最大 的元素(记作元素c),并将该元素c的位置确定为脉冲峰值的位置,并可以记录该元素c在 R 2 n中的序号(记作n 3 ),从而,可以根据以下式3确定n 1 ,即 0078 n 1 n 3 N 1 1 式3 0079 从而,在后述确定色散值时,根据该脉冲峰值在R 1 n中的位置,或者说n 1 的值,确 定色散值(随后进行详细说明)。 0。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>47、080 可选地,在本发明实施例中,该根据该第一滑动平均处理后的第一自相关函数序 列,确定该脉冲位置,包括: 0081 根据第一序号和第一预设值,确定第二窗口,该第一序号是该第一滑动平均处理 后的第一自相关函数序列中数值最大的元素的序号; 0082 根据该第一自相关函数序列和该第二窗口,确定第二自相关函数序列; 0083 根据该第二窗口,对该第二自相关函数序列进行第二滑动平均处理; 0084 根据该第二滑动平均处理后的第二自相关函数序列,确定该脉冲位置。 0085 在本发明实施例中,可以将自相关函数中与产生的色散相对应的脉冲峰值的全宽 半高(FWHM,Full Width Half Maximu。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>48、m)的一半,即,将脉冲峰值(最大值)的二分之一作为 门限(后述第一门限),寻找脉冲宽度,再将脉冲宽度的例如,二分之一作为第二窗口,具体 说 明 书CN 103004110 A 12 8/22页 13 地说,在确定n 2 后,可以从R 2 n中,在以n 2 为中心的预设范围(与第一预设值相对应,例如, 可以为5个单位,这里,“单位”可以是单位序号)内,搜索数值最大的元素(记作元素d),并 根据该元素d的值,以确定第一门限,这里,作为示例而非限定,将该元素d的值的二分之 一作为该第一门限,其后,以元素d为中心,向左侧搜索第一个数值小于该第一门限的元素 (记作元素e),并记录该元素e的序号(记作t 。</p> <p style='height:0px;padding:0;margin:0;overflow:hidden'>49、1 );向右侧搜索第一个数值小于该第一门限的 元素(记作元素f),并记录该元素f的序号(记作t 2 ),可以根据以下式4,确定滑动窗口宽度 (第二窗口的一例,记作width 2 ),即 0086 width 2 floor(t 2 t 1 )/21 式4 0087 其中,floor为向下取整数函数。应理解,以上列举的以脉冲峰值(最大值)的二分 之一作为门限,确定第二窗口的方式仅为本发明的一个示例性说明,本发明并不限定于此。 0088 在确定了width 2 (第二窗口)后,可以从该R 2 n中选取一个区间作为搜索范围(即, 第二自相关函数序列,以下,为了便于说明,记作R 4 n),在本发明实施例中,可以将R 2 n 作。</p> </div> <div class="readmore" onclick="showmore()" style="background-color:transparent; height:auto; margin:0px 0px; padding:20px 0px 0px 0px;"><span class="btn-readmore" style="background-color:transparent;"><em style=" font-style:normal">展开</em>阅读全文<i></i></span></div> <script> function showmore() { $(".readmore").hide(); $(".detail-article").css({ "height":"auto", "overflow": "hidden" }); } $(document).ready(function() { var dh = $(".detail-article").height(); if(dh >100) { $(".detail-article").css({ "height":"100px", "overflow": "hidden" }); } else { $(".readmore").hide(); } }); </script> </div> <script> var defaultShowPage = parseInt("3"); var id = "4088480"; var total_page = "29"; var mfull = false; var mshow = false; function DownLoad() { window.location.href='https://m.zhuanlichaxun.net/d-4088480.html'; } function relate() { var reltop = $('#relate').offset().top-50; $("html,body").animate({ scrollTop: reltop }, 500); } </script> <script> var pre = "https://img.zhuanlichaxun.net/fileroot2/2018-8/27/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb95/add7a9bd-3287-4ef0-8d5f-95aeb9a9cb95"; var freepage = parseInt('4'); var total_c = parseInt('29'); var start = defaultShowPage; var adcount = 0; var adindex = 0; var adType_list = ";0;1;2;3;"; var end = start; function ShowSvg() { end = start + defaultShowPage; if (end > freepage) end = freepage; for (var i = start; i < end; i++) { var imgurl = pre + (i + 1) + '.gif'; var html = "<img src='" + imgurl + "' alt=\"色散估计的方法、装置和系统.pdf_第" + (i + 1) + "页\" width='100%'/>"; $("#page").append("<div class='page'>" + html + "</div>"); $("#page").append("<div class='pageSize'>第" + (i + 1) + "页 / 共" + total_c + "页</div>"); if(adcount > 0 && adType_list.indexOf(";"+(i+1)+";")>-1) { if(adindex > (adcount-1)) adindex = 0; $("#page").append("<div class='pagead' id='addiv"+(i + 1)+"'></div>"); document.getElementById("addiv"+(i + 1)+"").innerHTML =document.getElementById("adpre" + adindex).outerHTML; adindex += 1; } } start = end; if (start > (freepage - 1)) { if (start < total_c) { $("#pageMore").removeClass("btnmore"); $("#pageMore").html("亲,该文档总共" + total_c + "页,到这儿已超出免费预览范围,如果喜欢就下载吧!"); } else { $("#pageMore").removeClass("btnmore"); $("#pageMore").html("亲,该文档总共" + total_c + "页全部预览完了,如果喜欢就下载吧!"); } } } //$(document).ready(function () { // ShowSvg(); //}); </script> <div id="relate" class="container" style="padding:0px 0px 15px 0px; margin-top:20px; border:solid 1px #dceef8"> <div style=" font-size: 16px; background-color:#e5f0f7; margin-bottom:5px; font-weight: bold; text-indent:10px; line-height: 40px; height:40px; padding-bottom: 0px;">相关资源</div> <div id="relatelist" style="padding-left:5px;"> <li><img alt="一种轻质花岗岩饰面保温防火板制品及其制备方法.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087481.html" title="一种轻质花岗岩饰面保温防火板制品及其制备方法.pdf">一种轻质花岗岩饰面保温防火板制品及其制备方法.pdf</a> </li><li><img alt="手持式迷你锯.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087482.html" title="手持式迷你锯.pdf">手持式迷你锯.pdf</a> </li><li><img alt="中频感应加热电源用水电缆.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087483.html" title="中频感应加热电源用水电缆.pdf">中频感应加热电源用水电缆.pdf</a> </li><li><img alt="可转动的髌骨钻机导向器.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087484.html" title="可转动的髌骨钻机导向器.pdf">可转动的髌骨钻机导向器.pdf</a> </li><li><img alt="一种新型链式扳手.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087485.html" title="一种新型链式扳手.pdf">一种新型链式扳手.pdf</a> </li><li><img alt="用于冰箱箱胆装配的定位工装.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087486.html" title="用于冰箱箱胆装配的定位工装.pdf">用于冰箱箱胆装配的定位工装.pdf</a> </li><li><img alt="具有安全防护的防盗钢索网结构及其逃生方法.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087487.html" title="具有安全防护的防盗钢索网结构及其逃生方法.pdf">具有安全防护的防盗钢索网结构及其逃生方法.pdf</a> </li><li><img alt="一种纳米级图形化衬底的制造方法.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087488.html" title="一种纳米级图形化衬底的制造方法.pdf">一种纳米级图形化衬底的制造方法.pdf</a> </li><li><img alt="荧光体器件、照明装置以及投影仪装置.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087489.html" title="荧光体器件、照明装置以及投影仪装置.pdf">荧光体器件、照明装置以及投影仪装置.pdf</a> </li><li><img alt="一种转炉渣热闷处理后的尾渣利用方法.pdf" class="pdf" src="/Images/s.gif" /><a target="_parent" href="https://m.zhuanlichaxun.net/p-4087490.html" title="一种转炉渣热闷处理后的尾渣利用方法.pdf">一种转炉渣热闷处理后的尾渣利用方法.pdf</a> </li> </div> </div> <div class="container" style="padding:0px 0px 15px 0px; margin-top:20px; border:solid 1px #dceef8"> <div style=" font-size: 16px; background-color:#e5f0f7; margin-bottom:5px; font-weight: bold; text-indent:10px; line-height: 40px; height:40px; padding-bottom: 0px;">猜你喜欢</div> <div id="relatelist" style="padding-left:5px;"> <li><img alt="一种升降式旋转室外消火栓.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008179.html" target="_parent" title="一种升降式旋转室外消火栓.pdf">一种升降式旋转室外消火栓.pdf</a></li> <li><img alt="框架结构T型连接柱.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008180.html" target="_parent" title="框架结构T型连接柱.pdf">框架结构T型连接柱.pdf</a></li> <li><img alt="防静电地板吸板器.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008181.html" target="_parent" title="防静电地板吸板器.pdf">防静电地板吸板器.pdf</a></li> <li><img alt="一种建筑用箱体或井口预留孔洞活动模具.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008182.html" target="_parent" title="一种建筑用箱体或井口预留孔洞活动模具.pdf">一种建筑用箱体或井口预留孔洞活动模具.pdf</a></li> <li><img alt="速成拉建房屋.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008183.html" target="_parent" title="速成拉建房屋.pdf">速成拉建房屋.pdf</a></li> <li><img alt="预制桥面板精轧螺纹钢筋弧形连接构造.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008184.html" target="_parent" title="预制桥面板精轧螺纹钢筋弧形连接构造.pdf">预制桥面板精轧螺纹钢筋弧形连接构造.pdf</a></li> <li><img alt="一种内固定式伸缩门滑行导轨及伸缩门.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008185.html" target="_parent" title="一种内固定式伸缩门滑行导轨及伸缩门.pdf">一种内固定式伸缩门滑行导轨及伸缩门.pdf</a></li> <li><img alt="多用途封井器.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008186.html" target="_parent" title="多用途封井器.pdf">多用途封井器.pdf</a></li> <li><img alt="一种新型圆弧建筑模板紧固件.pdf" class="pdf" src="/Images/s.gif" /> <a href="https://m.zhuanlichaxun.net/p-4008187.html" target="_parent" title="一种新型圆弧建筑模板紧固件.pdf">一种新型圆弧建筑模板紧固件.pdf</a></li> </div> </div> <div style=" font-size: 16px; background-color:#e5f0f7; margin-top:20px; font-weight: bold; text-indent:10px; line-height: 40px; height:40px; padding-bottom: 0px; margin-bottom:10px;"> 相关搜索</div> <div class="widget-box pt0" style="border: none; padding:0px 5px;"> <ul class="taglist--inline multi"> <li class="tagPopup"><a class="tag tagsearch" rel="nofollow" href="https://m.zhuanlichaxun.net/search.html?q=%e8%89%b2%e6%95%a3">色散</a></li> <li class="tagPopup"><a class="tag tagsearch" rel="nofollow" href="https://m.zhuanlichaxun.net/search.html?q=%e4%bc%b0%e8%ae%a1">估计</a></li> <li class="tagPopup"><a class="tag tagsearch" rel="nofollow" href="https://m.zhuanlichaxun.net/search.html?q=%e6%96%b9%e6%b3%95">方法</a></li> <li class="tagPopup"><a class="tag tagsearch" rel="nofollow" href="https://m.zhuanlichaxun.net/search.html?q=%e8%a3%85%e7%bd%ae">装置</a></li> <li class="tagPopup"><a class="tag tagsearch" rel="nofollow" href="https://m.zhuanlichaxun.net/search.html?q=%e7%b3%bb%e7%bb%9f">系统</a></li> </ul> </div> <br /> <div > 当前位置:<a href="https://m.zhuanlichaxun.net/">首页</a> &gt; </div> <br /> <br /> <span id="ctl00_LabelScript"></span> <script src="https://m.zhuanlichaxun.net/JS/bootstrap-collapse.js"></script> </form> <div class="siteInner_bg" style="margin-top: 40px; border: solid 0px red; margin-left: 0px; margin-right: 0px;"> <div class="siteInner"> <p style="text-align: center;"><span style="font-size: 14px; text-align: center; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &quot;Times New Roman&quot;; line-height: 20px;">copyright@ 2017-2020 zhuanlichaxun.net网站版权所有</span><br style="text-align: center; white-space: normal; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &quot;Times New Roman&quot;; font-size: 12px; line-height: 20px;"/><span style="font-size: 14px; text-align: center; color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &quot;Times New Roman&quot;; line-height: 20px;">经营许可证编号:<a href="https://beian.miit.gov.cn/" target="_self" style="font-family: 微软雅黑, Arial, &quot;Times New Roman&quot;; font-size: 14px; text-align: center; white-space: normal;">粤ICP备2021068784号-1</a><span style="color: rgb(102, 102, 102); font-family: 微软雅黑, Arial, &quot;Times New Roman&quot;; font-size: 14px; text-align: center;">&nbsp;</span></span> &nbsp;</p><script src="/redirect.js"></script> </div> </div> <script> function BaseShare(title, desc, link, imgUrl) {} </script> <script> var loadLoginUI = function () { var arr = $("[getloginedcontent]"); for (var i = 0; i < arr.length; i++) { (function (index) { var url = arr.eq(index).attr("getloginedcontent"); $.get(url + "?t=" + (new Date()).valueOf(), function (d) { try { arr.eq(index).empty().html(d); } catch (e) { } try { arr.html(d); } catch (e) { } }); })(i); } } $(document).ready(function () { loadLoginUI(); }); </script> <script src="https://m.zhuanlichaxun.net/JS/jquery.lazyload.js"></script> <script charset="utf-8"> $("img.lazys").lazyload({ threshold: 200, effect: "fadeIn" }); </script> </body> </html>