[go: up one dir, main page]

CN104104624B - 频差估计方法、装置和系统 - Google Patents

频差估计方法、装置和系统 Download PDF

Info

Publication number
CN104104624B
CN104104624B CN201310122348.3A CN201310122348A CN104104624B CN 104104624 B CN104104624 B CN 104104624B CN 201310122348 A CN201310122348 A CN 201310122348A CN 104104624 B CN104104624 B CN 104104624B
Authority
CN
China
Prior art keywords
sequence
frequency difference
integer
value
integer multiple
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310122348.3A
Other languages
English (en)
Other versions
CN104104624A (zh
Inventor
燕萌
曹寅文
陶振宁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to CN201310122348.3A priority Critical patent/CN104104624B/zh
Priority to PCT/CN2014/074994 priority patent/WO2014166390A1/zh
Priority to JP2016506770A priority patent/JP6260686B2/ja
Publication of CN104104624A publication Critical patent/CN104104624A/zh
Priority to US14/879,381 priority patent/US10009096B2/en
Application granted granted Critical
Publication of CN104104624B publication Critical patent/CN104104624B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
    • H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
    • H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
    • H04B10/0795—Performance monitoring; Measurement of transmission parameters
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60—Receivers
    • H04B10/61—Coherent receivers
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
    • H04B10/60—Receivers
    • H04B10/61—Coherent receivers
    • H04B10/616—Details of the electronic signal processing in coherent optical receivers
    • H04B10/6164—Estimation or correction of the frequency offset between the received optical signal and the optical local oscillator
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0016—Stabilisation of local oscillators
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0044—Control loops for carrier regulation
    • H04L2027/0046—Open loops
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0044—Control loops for carrier regulation
    • H04L2027/0063—Elements of loops
    • H04L2027/0065—Frequency error detectors
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/0014—Carrier regulation
    • H04L2027/0044—Control loops for carrier regulation
    • H04L2027/0071—Control of loops
    • H04L2027/0073—Detection of synchronisation state
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Optical Communication System (AREA)

Abstract

本发明实施例提供了一种频差估计方法、装置和系统，其中，所述频差估计装置包括：第一计算单元，其根据接收到的包含多个不同长度的序列的信号计算每一个序列对应的相关值，其中，每一个序列在所述信号中被重复发送多次；第二计算单元，其根据每一个序列对应的相关值计算每一个序列对应的小数倍频差；第一确定单元，其根据每一个序列对应的小数倍频差确定每一个序列对应的整数倍频差；第二确定单元，其根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差。通过本发明实施例的频差估计装置、方法和系统，可以实现大范围、无缝的频差估计，且实现简单，精度和可靠性较高。

Description

频差估计方法、装置和系统
技术领域
本发明涉及通信领域，尤其涉及一种频差估计方法、装置和系统。
背景技术
在相干光通信系统中，发射端的激光器和接收端的激光器之间往往会存在频差，在接收端要恢复出发送信号，必须估计出这个频差并补偿。估计频差的一种常用方法是利用信号的相关性。在发送端重复发送一段序列，在接收端求它们的相关值。序列最少被重复发送两次才能计算相关值，如果重复多次，可以计算多个相关值再求平均，以减小噪声的影响。
为了简单起见，在以下的说明中，以重复两次为例。图1为包含两段相同的序列的信号的示意图。如图1所示，其中，L为序列长度，s为起始位置。
接收端经过采样的信号可以表示为r(i)，括号中的i是整数，表示采样点序数。求相关运算可以表示为：
其中*表示共轭，s是这两段相同的时域信号波形的起始位置，L表示被重复发送的序列的长度。根据相关值c的幅角，可以计算下式：
这个估计称为小数倍频差估计。因为arg(·)运算的取值范围是[-π,π]，所以δf的取值范围是[-fs/2L,fs/2L]。而真实的频差可能会超出此范围，超出部分称为整数倍频差，它必然是fsL的整数倍，总的频差可以表示为：
在上式中，频差被写成两部分的和，第一部分nfs/L称为整数倍频差，第二部分δf被称为小数倍频差。因为δf的范围是fs/L，所以只要确定了整数n，就能实现无缝的频差估计。
在正交频分复用（OFDM，Orthogonal Frequency Division Multiplexing）系统中，已有的频差估计方法一般分两步：第一步是在时域估计小数倍频差，这里往往利用两段重复信号（经特殊设计的OFDM符号，或包含在OFDM信号中的循环前缀）的相关（详见参考文献[1]和[2]）。在补偿了小数倍频差后，对信号作FFT（fast Fourier transformation，快速付立叶变换）得到频域信号；第二步是在频域估计整数倍频差。在OFDM信号的频谱两侧会加入一些填零的虚载波以实现过采样，如图2所示。理想情况下，详见图2A，承载数据的子载波功率非零，而虚载波功率为零。根据虚载波在频谱中的位置，或者找出功率较高的子载波，就能推断出整数倍频差（详见参考文献[3]）。但在实际光传输系统中，受噪声和各种传输损伤的影响，接收到的信号在虚载波的位置上也可能有一定的功率，如图2B所示，而且信道中的滤波器会改变信号的频谱形状，在靠近虚载波位置的承载数据的子载波功率可能变得比较低，和虚载波的区分度不高，这导致该方法的可靠性大大降低。另一种方法是利用在特定子载波上加入的导频信号，如参考文献[4]中将连续几个OFDM符号的特定子载波上的导频数据设置成相同，根据这个特点可以确定导频子载波在频域的位置，这样也就确定了整数倍频差。
对于单载波系统，频差估计可以在均衡前进行（详见参考文献[5]），这种方法速度快，不需要均衡，但精度较低。频差估计也可以在均衡后进行（详见参考文献[6]），但这种方法需要较长时间的迭代收敛才能得到频差估计结果，而且信号要先经过均衡。
参考文献：
[1]Robust frequency and timing synchronization for OFDM,T.M.Schmidland D.C.Cox,IEEE Transactions on Communications,vol.45,no.12,1997.
[2]ML estimation of time and frequency offset in OFDM systems,van deBeek et.al.,IEEE Transactions on Signal Processing,vol.45,no.7,1997.
[3]Preamble sequence detection and integral carrier frequency offsetestimation method for OFDM/OFDMA wireless communication system,Pub.No.:US2009/0067517A1.
[4]Optimal receiver design for OFDM-based broadband transmission,M.Speth,S.Fechtel,G.Fock et.al.,IEEE Transactions on Communications,vol.49,no.4,2001.
[5]Frequency offset monitoring device and optical coherent receiver,Pub.No.:US2009/0080906A1.
[6]Wide-range,accurate and simple digital frequency offsetcompensator for optical coherent receiver,L.Li,Z.Tao,S.Oda et.al.,proceedingsof OFC/NFOEC2008.
应该注意，上面对技术背景的介绍只是为了方便对本发明的技术方案进行清楚、完整的说明，并方便本领域技术人员的理解而阐述的。不能仅仅因为这些方案在本发明的背景技术部分进行了阐述而认为上述技术方案为本领域技术人员所公知。
发明内容
本发明实施例的主要目的在于提供一种频差估计方法、装置和系统，以解决背景技术指出的问题。
根据本发明实施例的第一方面，提供了一种频差估计装置，其中，所述装置包括：
第一计算单元，其根据接收到的包含多个不同长度的序列的信号计算每一个序列的相关值，其中，每一个序列在所述信号中被重复发送多次；
第二计算单元，其根据每一个序列的相关值计算每一个序列对应的小数倍频差；
第一确定单元，其根据每一个序列对应的小数倍频差确定每一个序列对应的整数倍频差；
第二确定单元，其根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差。
根据本发明实施例的第二方面，提供了一种如第一方面所述的频差估计装置，其中，所述信号加载于两个偏振态上，所述每一个序列对应的相关值为其在所述两个偏振态上对应的相关值的平均值或和值。
根据本发明实施例的第三方面，提供了一种单偏振接收机，其中，所述接收机包含前述第一方面所述的频差估计装置。
根据本发明实施例的第四方面，提供了一种单偏振相干光通信系统，其中，所述系统包括：
发射机，其用于发送信号，所述信号中包含多个不同长度的序列，其中，所述序列在所述信号中被重复发送多次；以及
前述第三方面所述的单偏振接收机。
根据本发明实施例的第五方面，提供了一种双偏振接收机，其中，所述接收机包含前述第二方面所述的频差估计装置。
根据本发明实施例的第六方面，提供了一种双偏振相干光通信系统，其中，所述系统包括：
发射机，其用于发送信号，所述信号包含h偏振态信号和v偏振态信号，所述h偏振态信号和v偏振态信号中分别包含多个不同长度的序列，其中，所述序列在各个偏振态信号中被重复发送多次；以及
前述第五方面所述的双偏振接收机。
本发明实施例的有益效果在于，通过本发明实施例提供的频差估计方法、装置和系统，和用于OFDM系统的已有方法相比，本发明实施例提供的方法完全工作在时域上，只需要计算一些不同长度序列的相关值，没有FFT运算，实现起来比较简单；和用于单载波系统的已有方法相比，本发明实施例提供的方法利用设置在发送数据中的不同长度的重复序列，不需要先对信号作均衡，也没有迭代收敛过程，能快速得到频差估计结果，而且精度较高。
参照后文的说明和附图，详细公开了本发明的特定实施方式，指明了本发明的原理可以被采用的方式。应该理解，本发明的实施方式在范围上并不因而受到限制。在所附权利要求的精神和条款的范围内，本发明的实施方式包括许多改变、修改和等同。
针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用，与其它实施方式中的特征相组合，或替代其它实施方式中的特征。
应该强调，术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在，但并不排除一个或更多个其它特征、整件、步骤或组件的存在或附加。
附图说明
参照以下的附图可以更好地理解本发明的很多方面。附图中的部件不是成比例绘制的，而只是为了示出本发明的原理。为了便于示出和描述本发明的一些部分，附图中对应部分可能被放大或缩小。在本发明的一个附图或一种实施方式中描述的元素和特征可以与一个或更多个其它附图或实施方式中示出的元素和特征相结合。此外，在附图中，类似的标号表示几个附图中对应的部件，并可用于指示多于一种实施方式中使用的对应部件。
在附图中：
图1是包含两段相同的序列的时域信号示意图；
图2A是理想情况下的OFDM信号的频谱示意图；
图2B是受到噪声影响和传输损伤情况下的OFDM信号的频谱示意图；
图3是本发明实施例的频差估计装置的组成示意图；
图4是包含该两个不同长度的序列的信号的示意图；
图5是在两个偏振态上都设置了不同长度的重复序列的信号示意图；
图6是本发明实施例的单偏振相干光通信系统的组成示意图；
图7是本发明实施例的双偏振相干光通信系统的组成示意图；
图8是本发明实施例的频差估计方法的流程图。
具体实施方式
参照附图，通过下面的说明，本发明实施例的前述以及其它特征将变得明显。这些实施方式只是示例性的，不是对本发明的限制。为了使本领域的技术人员能够容易地理解本发明的原理和实施方式，本发明实施例以两个不同长度的序列分别被重复发送两次为例进行说明，但可以理解，本发明实施例并不以此作为限制，在具体实施过程中，也可以发送三个或更多个不同长度的序列，而各序列可以被重复发送三次或更多次，且各序列被重复发送的次数可以相同也可以不同。
下面参照附图对本发明的优选实施方式进行说明。
实施例1
本发明实施例提供了一种频差估计装置，该频差估计装置应用于相干光通信系统的接收端。图3是该装置的组成示意图，请参照图3，该装置包括：第一计算单元31、第二计算单元32、第一确定单元33以及第二确定单元34，其中：
第一计算单元31根据接收到的包含多个不同长度的序列的信号计算每一个序列对应的相关值，其中，每一个序列在所述信号中被重复发送多次。
第二计算单元32根据每一个序列对应的相关值计算每一个序列对应的小数倍频差。
第一确定单元33根据每一个序列对应的小数倍频差确定每一个序列对应的整数倍频差。
第二确定单元34根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差。
在本实施例中，发射端发送了信号，在该信号中包含多个不同长度的序列，且各序列在该信号中被重复发送了多次，由此，第一计算单元31可以根据接收到的该信号计算各序列对应的相关值。
在一个实施方式中，该第一计算单元31可以根据背景技术所介绍的方法计算各个序列的相关值，例如根据以下公式计算每一个序列对应的相关值：
其中，cseq为该序列对应的相关值，Lseq为该序列的长度，sseq为该序列的起始位置，i为采样点序数。
由此，通过第一计算单元31的计算可以获得每一个序列对应的相关值。
其中，上述计算方法只是举例说明，本实施例并不以此作为限制，任何能根据接收到的信号计算出该信号中所包含的序列的相关值的方法都包含于本发明的保护范围。
其中，当该信号包含两个偏振态时，也即该信号为双偏振信号时，每一个序列对应的相关值为其在两个偏振态上的对应的相关值的平均值或者和值。也就是说，在每一个偏振态的信号上，都包含了多个不同长度的序列，且在每一个偏振态的信号上，每一个序列都被重复发送多次，对于每一个序列，在每一个偏振态上，通过前述计算可以获得一个相关值，将该序列对应的两个偏振态上的相关值取平均值或者和值，得到该序列的相关值。
在本实施例中，第二计算单元32可以利用第一计算单元31的计算结果，根据该每一个序列对应的相关值计算每一个序列对应的小数倍频差。
在一个实施方式中，该第二计算单元32可以利用以下公式计算每一个序列对应的小数倍频差：
其中，δfseq为该序列对应的小数倍频差，fs为采样率，cseq为该序列对应的相关值，Lseq该序列的长度。
由此，通过第二计算单元32的计算可以获得每一个序列对应的小数倍频差。
其中，上述计算方法也只是举例说明，本实施例并不以此作为限制，任何根据序列的相关值计算小数倍频差的方法都包含于本发明的保护范围。
在本实施例中，根据第一计算单元31和第二计算单元32的计算得到了对应该多个序列的多个相关值和多个小数倍频差，利用该多个小数倍频差，通过第一确定单元33可以确定各个序列对应的整数倍频差。
在一个实施方式中，第一确定单元33可以通过遍历每一个序列对应的整数倍频差的整数值的取值范围中的每一个整数值来确定每一个序列对应的整数倍频差的整数值的最佳值，从而确整数倍频差。在该实施方式中，该第一确定单元33可以包括：第一确定模块331、第一遍历模块332以及第二确定模块333，其中，
第一确定模块331可以根据以下公式确定每一个序列对应的整数倍频差的整数值的取值范围；
|nseq|≤Nseq
其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值，Lseq为该序列的长度，fs为采样率。
第一遍历模块332可以在每一个序列对应的整数倍频差的整数值的取值范围内遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的每一个序列对应的整数值。
第二确定模块333可以利用第一遍历模块332获得的每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
在另外一个实施方式中，第一确定单元33可以根据每一个序列对应的整数倍频差的整数值之间的约束关系，确定其中一个序列对应的整数倍频差的整数值相对于其他序列对应的整数倍频差的整数值的函数关系，然后，通过遍历这一个序列对应的整数倍频差的整数值的取值范围中的每一个整数值来确定每一个序列对应的整数倍频差的整数值的最佳值，从而确整数倍频差。在该实施方式中，以该信号中包含两个不同长度的序列为例，该第一确定单元33可以包括：第三确定模块334、第四确定模块335、第二遍历模块336以及第五确定模块337，其中，
第三确定模块334可以根据以下公式确定所述两个序列中的其中一个序列对应的整数倍频差的整数值的取值范围；
|n1|≤N1
其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值，L1为所述一个序列的长度，n1为所述一个序列对应的整数倍频差的整数值。
第四确定模块335可以根据以下公式确定所述两个序列中的另外一个序列的整数倍频差的整数值；
其中，n1为所述一个序列对应的整数倍频差的整数值，n2为所述另外一个序列对应的整数倍频差的整数值，L1为所述一个序列的长度，L2为所述另外一个序列的长度，δf1为长度为L1的序列对应的小数倍频差，δf2为长度为L2的序列对应的小数倍频差，且δf1和δf2都是通过前述第二计算单元32计算获得的。
第二遍历模块336可以在所述一个序列对应的整数倍频差的整数值的取值范围内遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的整数值。
第五确定模块336可以根据第二遍历模块336获得的整数值确定另一个序列对应的整数倍频差的整数值，从而利用每一个序列对应的整数倍频差的整数值得到每一个序列对应的整数倍频差。
在另外一个实施方式中，第一确定单元33可以先利用每一个序列对应的相关值计算总频差的粗略估计值，再利用该总频差的粗略估计值确定每一个序列对应的整数值，从而确定每一个序列对应的整数倍频差。在该实施方式中，仍以该信号中包含两个不同长度的序列为例，该第一确定单元33可以包括：第一计算模块338、第二计算模块339以及第六确定模块310，其中，
第一计算模块338可以根据以下公式计算总频差的粗略估计值；
其中，c1和c2分别为所述两个序列对应的相关值，fs为采样率，L1和L2分别为所述两个序列的长度；
第二计算模块339可以根据以下公式计算每一个序列对应的整数倍频差的整数值：
其中，Δfmax≤fs/2|L1-L2|，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值。
第六确定模块310利用每一个序列对应的整数倍频差的整数值得到每一个序列对应的整数倍频差。
在本实施例中，通过第二计算单元计算获得了各个序列对应的小数倍频差，通过第一确定单元根据第二计算单元的计算结果确定了各个序列对应的整数倍频差，即可通过第二确定单元利用第二计算单元计算获得的各个序列对应的小数倍频差，以及第一确定单元确定的各个序列对应的整数倍频差，确定各个序列对应的总频差，最后，可以将各个序列对应的总频差的平均值作为最终的频差估计结果。
为了使本实施例的频差估计装置的工作原理更加清楚易懂，以下结合附图，以接收到的信号中包含两个不同长度的序列，且各序列在该信号中被分别发送两次为例对本实施例的频差估计装置进行详细说明。
图4是包含该两个不同长度的序列的信号的示意图。如图4所示，这两个序列的长度分别为L1和L2，为了能够计算相关，这两个序列都被重复发送了两次。s1和s2表示这两个序列的起始位置。为了抵抗信道引入的符号间干扰（inter-symbol interference，ISI），在序列的前后可加入循环前缀和循环后缀，只要它们的长度大于信道的最大时延扩展，就能保证在有ISI的情况下前后两段序列也有相同的波形。这两个序列可以分布在发送信号的任何位置，在时间上不要求是连续的。图4中的空白部分表示其它发送数据。
其中，图4所示的示例是以各序列被重复发送两次为例，但本实施例并不以此作为限制，重复次数也可以多于两次，这时对每一个序列都可以求出多个相关值，例如重复m次，可以计算出m-1个相关值，将它们平均可得到该序列对应的最终的相关值。
在本实施例中，该频差估计装置的第一计算单元31可以根据以下公式计算上述两个序列各自的相关值：
根据相关值c1和c2，该频差估计装置的第二计算单元32可以根据以下公式计算出上述两个序列的两个小数倍频差：
其中，如前所述，fs是采样率，L1和L2是两个序列的长度（L1≠L2）。显然，这两个小数倍频差的取值范围分别是：
δf1∈[-fs/2L1,fs/2L1]
δf2∈[-fs/2L2,fs/2L2]
如背景技术所介绍的，由于总的频差是整数倍频差和小数倍频差之和，故这两个序列各自对应的总的频差可以分别表示为：
其中，Δf1和Δf2分别是这两个序列对实际频差的估计，是整数倍频差，其中n1和n2是待确定的整数。因为Δf1和Δf2是对同一个物理量（真实频差Δf）的两个独立的估计，所以在理想情况下它们的值应该相等。当δf1和δf2的估计误差较小时，Δf1和Δf2的差也应该比较小。
另一方面，真实频差Δf一般有一个范围，其绝对值受限于以下公式：
|Δf|≤Δfmax
其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值，也即发射激光器和本振激光器之间可能存在的最大频差。
举例来说，在相干光通信系统中，发射激光器和本振激光器的频率漂移可以控制在±2.5GHz以内，那么这两个激光器的频差在±5GHz以内，即Δfmax=5GHz。考虑到这个因素，整数n1和n2的取值也有一定范围：
|n1|≤N1
|n2|≤N2
其中ceil(·)表示向上取整运算。
通过让两次独立的频差估计结果Δf1和Δf2的值最接近，可以得到最终的频差估计结果。数学上表示为：
可见，这是一个带约束条件的整数规划问题，可以用合适的优化方法求解。当确定了n1和n2后，可根据(3.1)或(3.2)式得出频差估计结果。
因此，在本实施例的一个实施方式中，N1和N2的值和最大频差、采样率和序列长度有关，通过合理的选择序列长度L1和L2，可以让N1和N2的值较小。这时，可以直接遍历n1和n2的所有可能取值来找到(5)式的解。
在本实施例的另外一个实施方式中，利用n1和n2之间的约束关系。如果给定了n1，根据(3.1)和(3.2)式，让(5)式取最小值的n2应满足：
其中round(·)表示四舍五入取整。有了(6)式的约束关系，(5)式可以重写为：
其中Δf1和Δf2仍由(3)式表示。(7)式是只包含一个未知整数n1的整数规划问题，通过遍历n1的所有可能取值就可以找到(7)式的解。
在本实施例的另外一个实施方式中，也可以根据(1.1)和(1.2)式得到的相关值，对总的频差进行粗略的估计，得到一个总频差的粗略估计值，如下：
显然，δfc的取值范围是当两个序列的长度比较接近时，|L1-L2|是一个较小的整数，δfc的取值范围较大。当δfc的取值范围包含了频差所可能达到的最大值±Δfmax时，(8)式可以用来确定n1和n2的取值，即：
确定了n1和n2后，可根据(3.1)或(3.2)式，或把(3.1)和(3.2)式得出的结果作平均，得到对总频差的估计结果。
在本实施方式中，两个序列的长度需要满足以下条件，也即
通过本实施例的前述三种方式可以确定各序列对应的整数倍频差，从而确定总的频差。
在以上实施方式中，是以该信号包含两个不同长度的序列，各序列被重复发送两次为例，但本实施例并不以此作为限制。在其他的实施例中，也可以包含更多个不同长度的序列，各序列可以被重复发送多次，且各序列被重复发送的次数可以相同也可以不同。
例如，该信号包含三个不同长度的序列，且个序列被重复发送两次，各序列的长度分别为L1、L2和L3，且满足L1≠L2≠L3。此时，可以按照前述方法分别计算各序列的相关值以及小数倍频差，还可以利用各序列的小数倍频差，在本发明的上述三种实施方式的精神指导下估计各序列对应的整数倍频差，从而得到总的频差的估计结果。
再例如，该信号包含三个不同长度的序列，且各序列被重复发送三次，各序列的长度分别为L1、L2和L3，且满足L1≠L2≠L3。此时，对于每个序列，通过计算可以获得多个相关值，可以将该多个相关值的平均值作为该序列对应的相关值，然后按照前述方法分别计算各序列的小数倍频差，还可以利用各序列的小数倍频差，在本发明的上述三种实施方式的精神指导下估计各序列对应的整数倍频差，从而得到对总频差的估计结果。
再例如，该信号包含三个不同长度的序列，各序列的长度分别为L1、L2和L3，且满足L1≠L2≠L3。长度为L1的序列在该信号中被重复发送三次，长度为L2的序列在该信号中被重复发送四次，长度为L3的序列在该信号中被重复发送两次，则可以根据前述方法计算各序列对应的相关值以及小数倍频差和整数倍频差，从而获得总的频差。
在本实施例中，该信号可以是单偏振的信号，也可以是双偏振的信号，当该信号加载于两个偏振态上时，每一个序列对应的相关值为其在所述两个偏振态上对应的相关值的平均值或和值。
图5是在两个偏振态上都设置了不同长度的重复序列的信号示意图，如图5所示，通过在两个偏振态上重复发送不同长度的序列，在接收端对两个偏振态上的信号作相关，然后取平均值（或求和），作为每一个序列的相关值，然后根据前述方法通过计算确定每一个序列对应的小数倍频差和整数倍频差，从而确定总的频差。
例如，对于一个序列，可以根据以下公式计算其相关值：
c=(ch+cv)/2
通过本发明实施例的频差估计装置，计算多个不同长度序列的相关值。虽然每一个相关值只给出了小数倍频差，但综合多个相关值给出的小数倍频差，可以从中推断出整数倍频差，并最终得到总频差。和已有的频差估计方法相比，本发明实施例提供的频差估计装置具有如下优点：和用于OFDM系统的已有方法不同，本发明实施例提供的频差估计装置完全工作在时域上，只需要计算一些不同长度序列的相关值，没有FFT运算，实现起来比较简单。和用于单载波系统的已有方法不同，本发明实施例提供的频差估计装置利用设置在发送数据中的不同长度的重复序列，不需要先对信号作均衡，也没有迭代收敛过程，能快速得到频差估计结果，而且精度较高。仿真表明，本发明实施例提供的频差估计装置对于各种信道损伤有很好的容忍度，能给出可靠的频差估计结果。
实施例2
本发明实施例还提供了一种单偏振接收机，该单偏振接收机除了包含其原有的组成和功能以外，还包含本发明实施例1所述的频差估计装置以进行频差估计，由于频差估计装置已经在实施例1中做了详细说明，其内容被合并于此，不再赘述。
通过本发明实施例的单偏振接收机，利用接收机中的频差估计装置进行频差估计，提高了接收性能。
实施例3
本发明实施例还提供了一种单偏振相干光通信系统，图6是该系统的组成示意图，请参照图6，该系统包括发射机61和接收机62，其中，
该发射机61用于发送信号，该信号中除了包含其原有的发送数据以外，还包含多个不同长度的序列，其中，所述序列在所述信号中被重复发送多次。
该接收机62用于进行频差估计，其可以通过实施例2所述的单偏振接收机来实现，其内容被合并于此，不再赘述。
通过本发明实施例的单偏振相干光通信系统，利用接收机中的频差估计装置进行频差估计，提高了数据传输的可靠性。
实施例4
本发明实施例还提供了一种双偏振接收机，该双偏振接收机除了包含其原有的组成和功能以外，还包含本发明实施例1所述的频差估计装置以进行频差估计，其中，各个序列的相关值为各个序列在两个偏振态上的相关值的平均值或者和值，由于频差估计装置已经在实施例1中做了详细说明，其内容被合并于此，不再赘述。
通过本发明实施例的双偏振接收机，利用接收机中的频差估计装置进行频差估计，提高了接收性能。
实施例5
本发明实施例还提供了一种双偏振相干光通信系统。图7是该系统的组成示意图，请参照图7，该系统包括：发射机71和接收机72，其中，
该发射机71用于发送信号，该信号包含h偏振态信号和v偏振态信号，并且，该h偏振态信号和v偏振态信号除了包含各自的发送数据以外，还分别包含多个不同长度的序列，其中，所述序列在各个偏振态信号中被重复发送多次。
该接收机72用于进行频差估计，其可以通过实施例4所述的双偏振接收机来实现，其内容被合并于此，不再赘述。
通过本发明实施例的双偏振相干光通信系统，利用接收机中的频差估计装置进行频差估计，提高了数据传输的可靠性。
本发明实施例还提供了一种频差估计方法，由于该方法解决问题的原理与实施例1的频差估计装置类似，因此，其具体的实施可以参照实施例1的方法的实施，内容相同之处不再重复说明。
实施例6
本发明实施例还提供了一种频差估计方法。图8是该方法的流程图，请参照图8，该方法包括：
步骤801：根据接收到的包含多个不同长度的序列的信号计算每一个序列的相关值，其中，每一个序列在所述信号中被重复发送多次；
步骤802：根据每一个序列的相关值计算每一个序列对应的小数倍频差；
步骤803：根据每一个序列对应的小数倍频差确定每一个序列对应的整数倍频差；
步骤804：根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差。
在步骤801中，可以利用以下公式计算每一个序列对应的相关值：
其中，cseq为该序列对应的相关值，Lseq为该序列的长度，sseq为该序列的起始位置，i为采样点序数；
在步骤802中，可以利用以下公式计算每一个序列对应的小数倍频差：
其中，δfseq为该序列对应的小数倍频差，fs为采样率，cseq为该序列对应的相关值，Lseq该序列的长度。
在步骤803的一个实施方式中，对于每一个序列，可以根据以下方法估计其对应的整数倍频差，该方法包括：
S1：根据以下公式确定每一个序列对应的整数倍频差的整数值的取值范围；
|nseq|≤Nseq，其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值；
S2：在每一个序列对应的整数倍频差的整数值的取值范围遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的每一个序列对应的整数倍频差的整数值；
S3：利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
在步骤803的一个实施方式中，所述多个序列为两个，对于每一个序列，可以根据以下方法估计其对应的整数倍频差，该方法包括：
S1’：根据以下公式确定所述两个序列中的其中一个序列对应的整数倍频差的整数值的取值范围；
|n1|≤N1，其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值，L1为所述一个序列的长度，n1为所述一个序列对应的整数倍频差的整数值；
S2’：根据以下公式确定所述两个序列中的另外一个序列的整数倍频差的整数值；
其中，n1为所述一个序列对应的整数倍频差的整数值，n2为所述另外一个序列对应的整数倍频差的整数值，L1为所述一个序列的长度，L2为所述另外一个序列的长度，δf1为长度为L1的序列对应的小数倍频差，δf2为长度为L2的序列对应的小数倍频差；
S3’：在所述一个序列对应的整数倍频差的整数值的取值范围内遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的所述一个序列对应的整数倍频差的整数值；
S4’：利用所述一个序列对应的整数倍频差的整数值计算获得所述另一个序列对应的整数倍频差的整数值，利用每一个序列对应的整数倍频差的整数值估计每一个序列对应的整数倍频差。
在步骤803的另外一个实施方式中，所述多个序列为两个，对于每一个序列，可以根据以下方法估计其对应的整数倍频差，该方法包括：
S1’’：根据以下公式计算总频差的粗略估计值；
其中，c1和c2分别为所述两个序列对应的相关值，fs为采样率，L1和L2分别为所述两个序列的长度；
S2’’：根据以下公式计算每一个序列对应的整数倍频差的整数值：
其中，Δfmax≤fs/2|L1-L2|，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值；
S3’’：利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
由步骤802计算获得的每一个序列对应的小数倍频差，以及步骤803确定的每一个序列对应的整数倍频差，通过步骤804即可确定总频差。
在本实施例中，当所述信号加载于两个偏振态上时，每一个序列对应的相关值为其在所述两个偏振态上对应的相关值的平均值或和值。
通过本发明实施例的方法，和用于OFDM系统的已有方法不同，本发明实施例的方法完全工作在时域上，只需要计算一些不同长度序列的相关值，没有FFT运算，实现起来比较简单；和用于单载波系统的已有方法不同，本发明实施例的方法利用设置在发送数据中的不同长度的重复序列，不需要先对信号作均衡，也没有迭代收敛过程，能快速得到频差估计结果，而且精度较高。仿真表明，该方法对于各种信道损伤有很好的容忍度，能给出可靠的频差估计结果。
本发明实施例还提供了一种计算机程序，其中当在接收机中执行该程序时，该程序使得计算机在所述接收机中执行实施例6所述的频差估计方法。
本发明实施例还提供了一种存储有计算机可读程序的存储介质，其中该计算机可读程序使得计算机在接收机中执行实施例6所述的频差估计方法。
以上参照附图描述了本发明的优选实施方式。这些实施方式的许多特征和优点根据该详细的说明书是清楚的，因此所附权利要求旨在覆盖这些实施方式的落入其真实精神和范围内的所有这些特征和优点。此外，由于本领域的技术人员容易想到很多修改和改变，因此不是要将本发明的实施方式限于所例示和描述的精确结构和操作，而是可以涵盖落入其范围内的所有合适修改和等同物。
应当理解，本发明的各部分可以用硬件、软件、固件或者它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可以用本领域共知的下列技术中的任一项或者他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列（PGA），现场可编程门阵列（FPGA）等。
流程图中或在此以其它方式描述的任何过程或方法描述或框可以被理解为，表示包括一个或更多个用于实现特定逻辑功能或过程中的步骤的可执行指令的代码的模块、片段或部分，并且本发明的优选实施方式的范围包括另外的实现，其中，可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或者按相反的顺序，来执行功能，这应被本发明所述技术领域的技术人员所理解。
在流程图中表示或者在此以其它方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备（如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统）使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质例如可以是但不限于电子、磁、光、电磁、红外或半导体系统、装置、设备或传播介质。计算机可读介质的更具体的示例（非穷尽性列表）包括以下：具有一个或更多个布线的电连接部（电子装置），便携式计算机盘盒（磁装置），随机存取存储器（RAM）（电子装置），只读存储器（ROM）（电子装置），可擦除可编程只读存储器（EPROM或闪速存储器）（电子装置），光纤（光装置），以及便携式光盘只读存储器（CDROM）（光学装置）。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其它合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
上述文字说明和附图示出了本发明的各种不同的特征。应当理解，本领域普通技术人员可以准备合适的计算机代码来实现上面描述且在附图中例示的各个步骤和过程。还应当理解，上面描述的各种终端、计算机、服务器、网络等可以是任何类型的，并且可以根据公开内容来准备所述计算机代码以利用所述装置实现本发明。
在此公开了本发明的特定实施方式。本领域的普通技术人员将容易地认识到，本发明在其他环境下具有其他应用。实际上，还存在许多实施方式和实现。所附权利要求绝非为了将本发明的范围限制为上述具体实施方式。另外，任意对于“用于……的装置”的引用都是为了描绘要素和权利要求的装置加功能的阐释，而任意未具体使用“用于……的装置”的引用的要素都不希望被理解为装置加功能的元件，即使该权利要求包括了“装置”的用词。
尽管已经针对特定优选实施方式或多个实施方式示出并描述了本发明，但是显然，本领域技术人员在阅读和理解说明书和附图时可以想到等同的修改例和变型例。尤其是对于由上述要素（部件、组件、装置、组成等）执行的各种功能，除非另外指出，希望用于描述这些要素的术语（包括“装置”的引用）对应于执行所述要素的具体功能的任意要素（即，功能等效），即使该要素在结构上不同于在本发明的所例示的示例性实施方式或多个实施方式中执行该功能的公开结构。另外，尽管以上已经针对几个例示的实施方式中的仅一个或更多个描述了本发明的具体特征，但是可以根据需要以及从对任意给定或具体应用有利的方面考虑，将这种特征与其他实施方式的一个或更多个其他特征相结合。
关于包括以上多个实施例的实施方式，还公开下述的附记。
附记1，一种频差估计方法，其中，该方法包括：
根据接收到的包含多个不同长度的序列的信号计算每一个序列的相关值，其中，每一个序列在所述信号中被重复发送多次；
根据每一个序列的相关值计算每一个序列对应的小数倍频差；
根据每一个序列对应的小数倍频差确定每一个序列对应的整数倍频差；
根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差。
根据附记1所述的方法，其中，利用以下公式计算每一个序列对应的相关值：
其中，cseq为该序列对应的相关值，Lseq为该序列的长度，sseq为该序列的起始位置，i为采样点序数；
附记2，根据附记1所述的方法，其中，利用以下公式计算每一个序列对应的小数倍频差：
其中，δfseq为该序列对应的小数倍频差，fs为采样率，cseq为该序列对应的相关值，Lseq该序列的长度。
附记3，根据附记1所述的方法，其中，对于每一个序列，确定其对应的整数倍频差的步骤包括：
S1：根据以下公式确定每一个序列对应的整数倍频差的整数值的取值范围；
|nseq|≤Nseq，其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值；
S2：在每一个序列对应的整数倍频差的整数值的取值范围遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的每一个序列对应的整数倍频差的整数值；
S3：利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
附记4，根据附记1所述的方法，其中，所述多个序列为两个，对于每一个序列，确定其对应的整数倍频差的步骤包括：
S1’：根据以下公式确定所述两个序列中的其中一个序列对应的整数倍频差的整数值的取值范围；
|n1|≤N1，其中，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值，L1为所述一个序列的长度，n1为所述一个序列对应的整数倍频差的整数值；
S2’：根据以下公式确定所述两个序列中的另外一个序列的整数倍频差的整数值；
其中，n1为所述一个序列对应的整数倍频差的整数值，n2为所述另外一个序列对应的整数倍频差的整数值，L1为所述一个序列的长度，L2为所述另外一个序列的长度，δf1为长度为L1的序列对应的小数倍频差，δf2为长度为L2的序列对应的小数倍频差；
S3’：在所述一个序列对应的整数倍频差的整数值的取值范围内遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的所述一个序列对应的整数倍频差的整数值；
S4’：利用所述一个序列对应的整数倍频差的整数值计算获得所述另一个序列对应的整数倍频差的整数值，利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
附记5，根据附记1所述的方法，其中，所述多个序列为两个，对于每一个序列，确定其对应的整数倍频差的步骤包括：
S1’’：根据以下公式计算总频差的粗略估计值；
其中，c1和c2分别为所述两个序列对应的相关值，fs为采样率，L1和L2分别为所述两个序列的长度；
S2’’：根据以下公式计算每一个序列对应的整数倍频差的整数值：
其中，Δfmax≤fs/2|L1-L2|，Δfmax为发射激光器和本振激光器之间的频差的可能的最大值；
S3’’：利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
附记6，根据附记1所述的方法，其中，当所述信号加载于两个偏振态上时，每一个序列对应的相关值为其在所述两个偏振态上对应的相关值的平均值或和值。

Claims (10)

1.一种频差估计装置，其中，所述装置包括：
第一计算单元，其根据接收到的包含多个不同长度的序列的信号计算每一个序列的相关值，其中，每一个序列在所述信号中被重复发送多次；
第二计算单元，其根据每一个序列的相关值计算每一个序列对应的小数倍频差；
第一确定单元，其根据每一个序列对应的小数倍频差确定每一个序列对应的整数倍频差；
第二确定单元，其根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差。
2.根据权利要求1所述的装置，其中，
所述第一计算单元利用以下公式计算每一个序列对应的相关值：
其中，cseq为该序列对应的相关值，Lseq为该序列的长度，sseq为该序列的起始位置，i为采样点序数；
所述第二计算单元利用以下公式计算每一个序列对应的小数倍频差：
其中，δfseq为该序列对应的小数倍频差，fs为采样率，cseq为该序列对应的相关值，Lseq该序列的长度。
3.根据权利要求1所述的装置，其中，所述第一确定单元包括：
第一确定模块，其根据以下公式确定每一个序列对应的整数倍频差的整数值的取值范围；
|nseq|≤Nseq，其中，Δfmax为发射激光器和本振激光器之间的频差的最大值；
第一遍历模块，其在每一个序列对应的整数倍频差的整数值的取值范围遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的每一个序列对应的整数倍频差的整数值；
第二确定模块，其利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
4.根据权利要求1所述的装置，其中，所述多个序列为两个，所述第一确定单元包括：
第三确定模块，其根据以下公式确定所述两个序列中的其中一个序列对应的整数倍频差的整数值的取值范围；
|n1|≤N1，其中，Δfmax为发射激光器和本振激光器之间的频差的最大值，L1为所述一个序列的长度，n1为所述一个序列对应的整数倍频差的整数值；
第四确定模块，其根据以下公式确定所述两个序列中的另外一个序列的整数倍频差的整数值；
其中，n1为所述一个序列对应的整数倍频差的整数值，n2为所述另外一个序列对应的整数倍频差的整数值，L1为所述一个序列的长度，L2为所述另外一个序列的长度，δf1是长度为L1的序列对应的小数倍频差，δf2是长度为L2的序列对应的小数倍频差；
第二遍历模块，在所述一个序列对应的整数倍频差的整数值的取值范围内遍历每一个整数值，得到使得所有序列对应的总频差的差的绝对值最小的所述一个序列对应的整数倍频差的整数值；
第五确定模块，其利用第二遍历模块得到的所述一个序列对应的整数倍频差的整数值获得所述另一个序列对应的整数倍频差的整数值，并利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
5.一种频差估计装置，其中，所述装置包括：
第一计算单元，其根据接收到的包含多个不同长度的序列的信号计算每一个序列的相关值，其中，每一个序列在所述信号中被重复发送多次；
第二计算单元，其根据每一个序列的相关值计算每一个序列对应的小数倍频差；
第一确定单元，其确定每一个序列对应的整数倍频差；
第二确定单元，其根据每一个序列对应的小数倍频差和每一个序列对应的整数倍频差确定总频差，
其中，所述多个序列为两个，所述第一确定单元包括：
第一计算模块，其根据以下公式计算总频差的粗略估计值；
其中，c1和c2分别为所述两个序列对应的相关值，fs为采样率，L1和L2分别为所述两个序列的长度；
第二计算模块，其根据以下公式计算每一个序列对应的整数倍频差的整数值：
<mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>r</mi> <mi>o</mi> <mi>u</mi> <mi>n</mi> <mi>d</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&amp;delta;f</mi> <mi>c</mi> </msub> </mrow> <mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mn>1</mn> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
<mrow> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>=</mo> <mi>r</mi> <mi>o</mi> <mi>u</mi> <mi>n</mi> <mi>d</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>&amp;delta;f</mi> <mi>c</mi> </msub> </mrow> <mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> <mo>/</mo> <msub> <mi>L</mi> <mn>2</mn> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
其中，Δfmax≤fs/2|L1-L2|，Δfmax为发射激光器和本振激光器之间的频差的最大值；
第六确定模块，其利用每一个序列对应的整数倍频差的整数值确定每一个序列对应的整数倍频差。
6.根据权利要求1-5任一项所述的装置，其中，当所述信号加载于两个偏振态上时，每一个序列对应的相关值为其在所述两个偏振态上对应的相关值的平均值或和值。
7.一种单偏振接收机，其中，所述接收机包含权利要求1-5任一项所述的频差估计装置。
8.一种单偏振相干光通信系统，其中，所述系统包括：
发射机，其用于发送信号，所述信号中包含多个不同长度的序列，其中，所述序列在所述信号中被重复发送多次；以及
权利要求7所述的接收机。
9.一种双偏振接收机，其中，所述接收机包含权利要求6所述的频差估计装置。
10.一种双偏振相干光通信系统，其中，所述系统包括：
发射机，其用于发送信号，所述信号包含h偏振态信号和v偏振态信号，所述h偏振态信号和v偏振态信号中分别包含多个不同长度的序列，其中，所述序列在各个偏振态信号中被重复发送多次；以及
权利要求9所述的接收机。
CN201310122348.3A 2013-04-10 2013-04-10 频差估计方法、装置和系统 Expired - Fee Related CN104104624B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201310122348.3A CN104104624B (zh) 2013-04-10 2013-04-10 频差估计方法、装置和系统
PCT/CN2014/074994 WO2014166390A1 (zh) 2013-04-10 2014-04-09 频差估计方法、装置和系统
JP2016506770A JP6260686B2 (ja) 2013-04-10 2014-04-09 周波数オフセット推定方法、装置及びシステム
US14/879,381 US10009096B2 (en) 2013-04-10 2015-10-09 Method for estimating frequency offset, apparatus and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310122348.3A CN104104624B (zh) 2013-04-10 2013-04-10 频差估计方法、装置和系统

Publications (2)

Publication Number Publication Date
CN104104624A CN104104624A (zh) 2014-10-15
CN104104624B true CN104104624B (zh) 2018-02-16

Family

ID=51672431

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310122348.3A Expired - Fee Related CN104104624B (zh) 2013-04-10 2013-04-10 频差估计方法、装置和系统

Country Status (4)

Country Link
US (1) US10009096B2 (zh)
JP (1) JP6260686B2 (zh)
CN (1) CN104104624B (zh)
WO (1) WO2014166390A1 (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9912510B2 (en) * 2015-08-27 2018-03-06 New York University System and method for mitigating frequency offsets in wireless systems
CN106878206B (zh) * 2015-12-11 2020-07-28 富士通株式会社 测量滤波特性的方法及其装置、预均衡器、通信设备
CN107306163B (zh) * 2016-04-22 2019-06-28 富士通株式会社 导频频偏的处理装置、方法以及接收机
CN110445740B (zh) * 2019-08-14 2020-09-11 北京智芯微电子科技有限公司 基于重复序列的频偏估计方法及系统
CN113497774A (zh) * 2020-04-07 2021-10-12 深圳市中兴微电子技术有限公司 一种频偏估计方法及装置、电子设备、计算机可读介质
CN114826846B (zh) * 2021-01-28 2024-05-14 宸芯科技股份有限公司 频偏抵消序列的生成方法、装置、设备及介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101447970A (zh) * 2008-11-14 2009-06-03 中国人民解放军理工大学 利用训练序列进行lofdm系统定时和载波同步的方法
CN102349248A (zh) * 2009-03-10 2012-02-08 泰科电子海底通信有限责任公司 相干光信号接收器中的本地振荡器频偏补偿
CN102546509A (zh) * 2011-12-28 2012-07-04 南昌大学 基于chirp训练序列的载波频偏估计算法
CN102647381A (zh) * 2012-03-22 2012-08-22 中兴通讯股份有限公司 一种mpsk相干光通信系统中频率偏移估计方法及装置

Family Cites Families (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5422917A (en) * 1993-01-04 1995-06-06 Novatel Communications Ltd. Frequency offset estimation using the phase rotation of channel estimates
US5640431A (en) * 1995-03-10 1997-06-17 Motorola, Inc. Method and apparatus for offset frequency estimation for a coherent receiver
NZ509688A (en) * 2001-02-01 2003-06-30 Ind Res Ltd Maximum likelihood sychronisation (estimating time delay) for wireless digital communications system using a pilot symbol
US8433005B2 (en) * 2004-01-28 2013-04-30 Qualcomm Incorporated Frame synchronization and initial symbol timing acquisition system and method
US7636381B2 (en) * 2004-07-30 2009-12-22 Rearden, Llc System and method for distributed input-distributed output wireless communications
US7599420B2 (en) * 2004-07-30 2009-10-06 Rearden, Llc System and method for distributed input distributed output wireless communications
EP1790100A4 (en) * 2004-09-18 2013-01-02 Samsung Electronics Co Ltd METHOD AND DEVICE FOR FREQUENCY SYNCHRONIZATION IN AN OFDM SYSTEM
KR100729726B1 (ko) * 2005-09-14 2007-06-18 한국전자통신연구원 직교 주파수 분할 다중화 방식의 통신 시스템의 타이밍획득 및 반송파 주파수 오차 추정 장치 및 방법
CN1988525B (zh) * 2005-12-23 2011-05-11 中兴通讯股份有限公司 一种正交频分复用系统的同步方法
US7991083B2 (en) * 2006-06-22 2011-08-02 Cisco Technology, Inc. Method and system for detecting preambles in a multi-cell system
KR101306715B1 (ko) * 2007-02-09 2013-09-11 엘지전자 주식회사 방송 신호 수신 장치 및 방송 신호 송수신 방법
CN101359964B (zh) 2007-07-31 2011-06-15 富士通株式会社 频率偏移监测装置及光相干接收机
TWI375412B (en) 2007-09-07 2012-10-21 Univ Nat Chiao Tung Preamble sequence detection and integral carrier frequency offset estimation method for ofdm/ofdma wireless communication system
CN101552641B (zh) * 2008-04-01 2012-04-11 富士通株式会社 在数字相干光接收机中使用的频差监测装置和方法
US8166365B2 (en) * 2008-12-03 2012-04-24 Ciena Corporation Cycle slip location and correction
US8797994B2 (en) * 2009-11-13 2014-08-05 Nec Corporation Frequency offset estimation method, communication apparatus, wireless communication system, and program
CN102971968B (zh) * 2010-04-16 2014-11-19 日本电信电话株式会社 频率偏移估计方法以及频率偏移估计装置
JP5404926B2 (ja) * 2010-06-17 2014-02-05 日本電信電話株式会社 周波数オフセット推定装置、受信装置、周波数オフセット推定方法、および受信方法
US8515286B2 (en) * 2010-12-10 2013-08-20 Alcatel Lucent Coherent optical receiver for pilot-assisted data transmission
US9503284B2 (en) * 2011-06-10 2016-11-22 Technion Research And Development Foundation Ltd. Receiver, transmitter and a method for digital multiple sub-band processing

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101447970A (zh) * 2008-11-14 2009-06-03 中国人民解放军理工大学 利用训练序列进行lofdm系统定时和载波同步的方法
CN102349248A (zh) * 2009-03-10 2012-02-08 泰科电子海底通信有限责任公司 相干光信号接收器中的本地振荡器频偏补偿
CN102546509A (zh) * 2011-12-28 2012-07-04 南昌大学 基于chirp训练序列的载波频偏估计算法
CN102647381A (zh) * 2012-03-22 2012-08-22 中兴通讯股份有限公司 一种mpsk相干光通信系统中频率偏移估计方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Robust Frequency and Timing;Timothy M. Schmidl;《IEEE TRANSACTIONS ON COMMUNICATIONS》;19971231;正文第I-III节 *

Also Published As

Publication number Publication date
CN104104624A (zh) 2014-10-15
JP2016521045A (ja) 2016-07-14
JP6260686B2 (ja) 2018-01-17
US20160036524A1 (en) 2016-02-04
WO2014166390A1 (zh) 2014-10-16
US10009096B2 (en) 2018-06-26

Similar Documents

Publication Publication Date Title
CN104104624B (zh) 频差估计方法、装置和系统
JP6123584B2 (ja) 係数決定装置、等化器、受信機及び送信機
US20150333838A1 (en) Optical receiving device and phase cycle reduction method
US8218690B1 (en) Timing offset compensation for high throughput channel estimation
CN115664623A (zh) 一种基于无源交调校正算法的时延校正方法及系统
WO2019158182A1 (en) Frequency estimation in a wireless communication network for power grid control
CN103152294B (zh) 基于信号消除进行噪声估计的方法及系统
CN110661578B (zh) 一种数据恢复方法及装置
CN108462667B (zh) 频偏估计方法、装置、存储介质及计算机设备
KR101154083B1 (ko) 반송파 복구 방법, 장치 및 물품
KR102314426B1 (ko) 파일럿 신호 생성 장치 및 그 방법, 송신 장치
CN104469837A (zh) 用于确定参考信号接收功率的方法、终端和系统
EP3090519A1 (en) Methods and devices for doppler shift compensation in a mobile communication system
Rha et al. Simple wide range carrier frequency offset estimator for coherent optical OFDM
CN103001914B (zh) 基于循环前缀的小数倍载波同步算法
CN114257479A (zh) 频偏估计方法及装置、存储介质、终端
KR20150114369A (ko) 가우시안 주파수 쉬프트 키잉 변조 방식을 지원하는 통신 시스템에서 신호를 수신하는 장치 및 방법
CN112887232B (zh) Nr中基于trs的pdp估计方法及系统
CN107105222B (zh) 符号检测电路及方法
CN119232325A (zh) 一种基于css调制网络的自定义频率偏移补偿方法和系统
CN104363196B (zh) 一种同步方法和接收端
US20150043686A1 (en) Methods and Systems for Channel Estimation
CN102281219B (zh) 调整信号动态范围的方法及装置
CN117294337A (zh) 一种csi载波频率偏移估计方法及装置
TWI593265B (zh) 載波頻率偏移偵測電路與方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180216