CN114422314A - 多载波鉴频方法及系统 - Google Patents

多载波鉴频方法及系统 Download PDF

Info

Publication number
CN114422314A
CN114422314A CN202111583131.3A CN202111583131A CN114422314A CN 114422314 A CN114422314 A CN 114422314A CN 202111583131 A CN202111583131 A CN 202111583131A CN 114422314 A CN114422314 A CN 114422314A
Authority
CN
China
Prior art keywords
carrier
frequency
fourier
band
pass filter
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.)
Granted
Application number
CN202111583131.3A
Other languages
English (en)
Other versions
CN114422314B (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.)
Jiangsu Keda Hengxin Semiconductor Technology Co ltd
Original Assignee
Jiangsu Keda Hengxin Semiconductor Technology Co 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 Jiangsu Keda Hengxin Semiconductor Technology Co ltd filed Critical Jiangsu Keda Hengxin Semiconductor Technology Co ltd
Priority to CN202111583131.3A priority Critical patent/CN114422314B/zh
Publication of CN114422314A publication Critical patent/CN114422314A/zh
Application granted granted Critical
Publication of CN114422314B publication Critical patent/CN114422314B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2649—Demodulators
    • H04L27/265—Fourier transform demodulators, e.g. fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators
    • 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
    • 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
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2657—Carrier synchronisation
    • 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
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2668—Details of algorithms
    • H04L27/2669—Details of algorithms characterised by the domain of operation
    • H04L27/2672—Frequency domain
    • 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/0024—Carrier regulation at the receiver end
    • H04L2027/0026—Correction of carrier offset

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Discrete Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Electromagnetism (AREA)
  • Circuits Of Receivers In General (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开了一种多载波鉴频方法，包括：S1、通过带通滤波器接收包含多载波的调顶信号，并选定其中一个载波的载波频率作为带通滤波器的中心频率；S2、计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数；S3、在连续多个傅里叶积分区间执行步骤S2，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数；S4、利用得到的多个傅里叶系数计算选定载波的频偏；S5、利用选定载波的频偏矫正选定载波的幅度。本发明多载波鉴频方法及系统使用滤波器，通过把选定的载波频率作为中心频率，将其它载波信号滤除以消除其影响，在选定的载波频率上计算傅里叶系数并估计载波的频偏，利用得到的频偏对幅度进行矫正，最终得到误差小的幅度值。

Description

多载波鉴频方法及系统
技术领域
本发明涉及光通信技术领域，特别涉及一种多载波鉴频方法及系统。
背景技术
在光纤通信中，为了有效地跟踪和管理通信网络中所使用光波的波长，让每个光波自带波长信息无疑会带来很大的便利。调顶技术是在几乎不影响主信号通路通信质量的前提下，通过在主信号上叠加调顶信号来传送运行维护管理(OAM)信息。例如在5G前传网络中，主信号是非归零(NRZ)格式的高速数字信号，调顶信号是曼彻斯特(Manchester)编码的低速调幅信号或者多载波调幅信号。
目前发布的调顶技术是为每个波长叠加一个小幅度的低频正弦信号，正弦信号的有无可以传递低速的数字信息，而低频正弦信号的频率被用来对该波长进行标注，也就是说该光波通过调顶信号携带了自身波长信息。调顶后的光波信号如图1所示。其中的载波信号就是叠加的低频正弦信号，它的有无传递的是低速的数字信息。
多载波接收机作为运维管理的一部分，作用是把调顶信号中载波的频率和幅度提取出来。频率用来确定所涉及到的光波波长，幅度用来确定光波的有无或者光波的强度。
在半有源设备里面，一个光感应器会把所有调制在光波上的低频正弦波收下来，经过鉴频后估算低频正弦波的幅度，来确定光纤中的光波强度是否满足要求。
常用的鉴频方法是把接收下来的包含所有载波的信号进行快速傅里叶变换，在变换后的频谱中提取鉴频结果。在特定频率处谱线的幅度就是载波的幅度，而谱线处的频率就是载波的频率。所有载波的频率和幅度可以在一次傅里叶变换后得到。
这种方法的缺点是从客户端发来的载波的频率和幅度的鉴频结果由于自身载波的幅度较小而受其它载波的影响很大，进而造成鉴频结果误差很大。当载波的频率出现频偏的时候，鉴频的误差会更大。
另外，快速傅里叶变换需要等到待处理的信号完全输入后才能开始处理，不但需要大的存储空间，信号处理引起的延迟也很大。
发明内容
本发明要解决的技术问题是提供一种可以准确估计频偏，并利用得到的频偏对幅度进行矫正以得到误差小的幅度值的多载波鉴频方法。
为了解决上述问题，本发明提供了一种多载波鉴频方法，所述多载波鉴频方法包括以下步骤：
S1、通过带通滤波器接收包含多载波的调顶信号，并选定其中一个载波的载波频率作为带通滤波器的中心频率；
S2、计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数；
S3、在连续多个傅里叶积分区间执行步骤S2，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数；
S4、利用得到的多个傅里叶系数计算选定载波的频偏；
S5、利用选定载波的频偏矫正选定载波的幅度。
作为本发明的进一步改进，步骤S2包括：
将调顶信号中包含的多载波表示为：
Figure BDA0003426795090000021
其中，n是载波编号；N是载波个数；An是幅度，为正的实数；φn0是载波的初始相位；fn是载波频率；Δfn是载波对于额定值的频偏；
假设fm是第m个载波频率额定值，则f(t)在fm处产生的傅里叶系数为：
Figure BDA0003426795090000022
其中，T是载波的周期，k是自然数。
作为本发明的进一步改进，步骤S4包括：
所述傅里叶系数的相位与傅里叶积分区间的顺次编号之间的线性关系，表示为：
φm(i)＝πkT·Δfm·i+φn0
其中，φm为傅里叶系数的相位，i为傅里叶积分区间的顺次编号，Δfm为第m个载波的频偏，n是载波编号，T是载波的周期，k是自然数，φn0是载波的初始相位；
上述公式表示的直线的斜率best可以通过最小二乘法得到，记为：
best＝πkT·Δfm
得到第m个载波的频偏为：
Figure BDA0003426795090000035
作为本发明的进一步改进，步骤S5包括：
将公式best＝πkT·Δfm代入以下公式：
Figure BDA0003426795090000031
令m＝n得到：
Figure BDA0003426795090000032
将公式
Figure BDA0003426795090000033
代入以下公式：
Fm＝Bm·ΔFmm
取绝对值，变换后得到：
Figure BDA0003426795090000034
其中，Am是利用载波的频偏矫正得到的载波的幅度。
作为本发明的进一步改进，还包括以下步骤：
S6、依次选定不同载波的载波频率作为带通滤波器的中心频率，并执行步骤S2-S5，得到所有载波的频偏并矫正所有载波的幅度。
为了解决上述问题，本发明还提供了一种多载波鉴频系统，所述多载波鉴频系统包括：
带通滤波器，用于接收包含多载波的调顶信号；
频率选定单元，用于选定其中一个载波的载波频率作为带通滤波器的中心频率；
傅里叶系数计算单元，用于计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数；并在连续多个傅里叶积分区间计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数；
频偏计算单元，用于利用得到的多个傅里叶系数计算选定载波的频偏；
幅度矫正单元，用于利用选定载波的频偏矫正选定载波的幅度。
作为本发明的进一步改进，所述计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数，包括：
将调顶信号中包含的多载波表示为：
Figure BDA0003426795090000041
其中，n是载波编号；N是载波个数；An是幅度，为正的实数；φn0是载波的初始相位；fn是载波频率；Δfn是载波对于额定值的频偏；
假设fm是第m个载波频率额定值，则f(t)在fm处产生的傅里叶系数为：
Figure BDA0003426795090000042
其中，T是载波的周期，k是自然数。
作为本发明的进一步改进，所述利用得到的多个傅里叶系数计算选定载波的频偏，包括：
所述傅里叶系数的相位与傅里叶积分区间的顺次编号之间的线性关系，表示为：
φm(i)＝πkT·Δfm·i+φn0
其中，φm为傅里叶系数的相位，i为傅里叶积分区间的顺次编号，Δfm为第m个载波的频偏，n是载波编号，T是载波的周期，k是自然数，φn0是载波的初始相位；
上述公式表示的直线的斜率best可以通过最小二乘法得到，记为：
best＝πkT·Δfm
得到第m个载波的频偏为：
Figure BDA0003426795090000051
作为本发明的进一步改进，所述利用选定载波的频偏矫正选定载波的幅度，包括：
将公式best＝πkT·Δfm代入以下公式：
Figure BDA0003426795090000052
令m＝n得到：
Figure BDA0003426795090000053
将公式
Figure BDA0003426795090000054
代入以下公式：
Fm＝Bm·ΔFmm
取绝对值，变换后得到：
Figure BDA0003426795090000055
其中，Am是利用载波的频偏矫正得到的载波的幅度。
作为本发明的进一步改进，所述幅度矫正单元还用于：
依次选定不同载波的载波频率作为带通滤波器的中心频率，得到所有载波的频偏并矫正所有载波的幅度。
本发明的有益效果：
本发明多载波鉴频方法及系统使用滤波器，通过把选定的载波频率作为中心频率，将其它载波信号滤除以消除其影响，在选定的载波频率上计算傅里叶系数并估计载波的频偏，利用得到的频偏对幅度进行矫正，最终得到误差小的幅度值。
其中，从信号输入开始就进行傅里叶系数的计算，不需要存储整个积分区间的信号，在积分结束后可以几乎不延迟地给出载波的频偏和矫正后的幅度。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其他目的、特征和优点能够更明显易懂，以下特举较佳实施例，并配合附图，详细说明如下。
附图说明
图1是调顶后的光波信号示意图；
图2是本发明优选实施例中多载波鉴频方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明，以使本领域的技术人员可以更好地理解本发明并能予以实施，但所举实施例不作为对本发明的限定。
实施例一
如图2所示，本实施例公开了一种多载波鉴频方法，其包括以下步骤：
S1、通过带通滤波器接收包含多载波的调顶信号，并选定其中一个载波的载波频率作为带通滤波器的中心频率。其中，带通滤波器把其它载波信号滤除以消除其影响。
S2、计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数。
具体地，将调顶信号中包含的多载波表示为：
Figure BDA0003426795090000071
其中，n是载波编号；N是载波个数；An是幅度，为正的实数；φn0是载波的初始相位；fn是载波频率；Δfn是载波对于额定值的频偏；
假设fm是第m个载波频率额定值，则f(t)在fm处产生的傅里叶系数为：
Figure BDA0003426795090000072
其中，T是载波的周期，k是自然数。公式(2)的计算可以从输入信号的开始就进行，在积分区间结束后，可以几乎无延迟地得到傅里叶系数。把公式(1)带入(2)得到：
Figure BDA0003426795090000073
其中，
Figure BDA0003426795090000074
Figure BDA0003426795090000075
取公式(5)的绝对值得到：
Figure BDA0003426795090000076
从公式(6)可知，当下面公式(7)中的条件满足时，|ΔFnm|＝0；
kT(fn-fm+Δfn)＝p，p＝±1，±2，… (7)
当下面公式(8)中的条件满足时，|ΔFnm|＝1；
fn-fm+Δfn＝0 (8)
S3、在连续多个傅里叶积分区间执行步骤S2，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数。
为了细分公式(3)中fm处的傅里叶系数的组成，把公式(3)改写成公式(9)，如下：
Figure BDA0003426795090000081
其中，等式右边第一项是接收到的第m个载波对fm处的傅里叶系数的贡献，第一个因子Bm是所需要的鉴频结果，第二个因子ΔFmm是由于载波频偏引起的，可以用本发明的方法来估计。等式右边第二项合式是其它载波对对fm处的傅里叶系数计算的干扰。在本发明中，应确保第二项合式足够小，不影响到第m个载波的傅里叶系数的计算。
公式(9)中第二项和式中的第一个因子Bn是第n个载波的幅度和初始相位。可以通过带通滤波器来滤除其它载波，从而使得这个因子很小，进而使得第二项合式对计算fm处的傅里叶系数的贡献可以忽略不记。如果带通滤波器对其它载波的衰减为Afilter，第二项合式中的第一个因子Bn在带通滤波器的输出端为：
Bn，out＝Afilter·Bn (10)
公式(9)中第二项合式中的第二个因子ΔFnm是由其它载波的频率距离第m个载波的频率的远近引起的。为了使得第二项合式对计算fm处的傅里叶系数的贡献可以忽略不记，第二个因子ΔFnm应该尽可能小。因为Δfn是一个随机的量，在确定ΔFnm尽可能小的条件的时候可以忽略。因而，公式(7)变成：
kT(fn-fm)＝p，p＝±1，±2，… (11)
通常与第m个载波相邻的第m-1个载波和第m+1个载波的影响最大，公式(11)变为：
Figure BDA0003426795090000082
公式(11)和(12)可以作为选择傅里叶系数积分区间的参考。
设当其它载波的阈值为Aint时，其它载波对计算fm处的傅里叶系数的贡献可以忽略不记，把对第二项合式中两个因子Bn和ΔFnm的要求合并得到：
|Afilter·Bn·ΔFnm|<Aint (13)
公式(13)可以看成是在确定了计算傅里叶系数的积分区间后，对带通滤波器带外衰减的要求。
S4、利用得到的多个傅里叶系数计算选定载波的频偏。
具体地，利用带通滤波器把其它的载波滤除后，就剩下fm处的载波，公(9)简化为：
Fm＝Bm·ΔFmm (14)
把公(4)和(5)代入到公式(14)得到：
Figure BDA0003426795090000091
由公式(15)得到Fm的相位为：
φm＝πkT·Δfm+φm0 (16)
可以看出，由载波频偏引起的相位变化是πkT·Δfm。如果对接收到的载波以积分区间[0，kT]为单位连续地进行傅里叶系数的计算，那么相邻两次计算得到的相位变化就是πkT·Δfm。计算得到Fm的相位φm与傅里叶系数积分区间的顺次编号i是线性关系，可以表述为：
φm(i)＝πkT·Δfm·i+φn0 (17)
公式(17)是一个直线方程，为了得到准确的频偏估计，需要在连续多个傅里叶积分区间进行多次的傅里叶系数计算。
公式(17)表示的直线的斜率best可以通过最小二乘法得到，记为：
best＝πkT·Δfm (18)
这样第m个载波的频偏就是：
Figure BDA0003426795090000092
S5、利用选定载波的频偏矫正选定载波的幅度。
具体地，得到了第m个载波的频偏后，把公式(18)代入到公式(6)，并令m＝n得到：
Figure BDA0003426795090000101
把公式(4)带入公式(14)并取绝对值，变换后得到：
Figure BDA0003426795090000102
其中，|Fm|根据公式(2)由傅里叶系数的计算得到；|ΔFmm|由公式(20)得到。公式(21)中的Am就是利用载波的频偏矫正得到的载波的幅度。
进一步地，本发明的多载波鉴频方法还包括以下步骤：
S6、依次选定不同载波的载波频率作为带通滤波器的中心频率，并执行步骤S2-S5，得到所有载波的频偏并矫正所有载波的幅度。
本发明多载波鉴频方法使用滤波器，通过把选定的载波频率作为中心频率，将其它载波信号滤除以消除其影响，在选定的载波频率上计算傅里叶系数并估计载波的频偏，利用得到的频偏对幅度进行矫正，最终得到误差小的幅度值。
其中，从信号输入开始就进行傅里叶系数的计算，不需要存储整个积分区间的信号，在积分结束后可以几乎不延迟地给出载波的频偏和矫正后的幅度。
实施例二
本实施例公开了一种多载波鉴频系统，所述多载波鉴频系统包括：
带通滤波器，用于接收包含多载波的调顶信号；
频率选定单元，用于选定其中一个载波的载波频率作为带通滤波器的中心频率；
傅里叶系数计算单元，用于计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数；并在连续多个傅里叶积分区间计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数；
频偏计算单元，用于利用得到的多个傅里叶系数计算选定载波的频偏；
幅度矫正单元，用于利用选定载波的频偏矫正选定载波的幅度。
具体地，所述傅里叶系数计算单元用于执行以下步骤：
将调顶信号中包含的多载波表示为：
Figure BDA0003426795090000111
其中，n是载波编号；N是载波个数；An是幅度，为正的实数；φn0是载波的初始相位；fn是载波频率；Δfn是载波对于额定值的频偏；
假设fm是第m个载波频率额定值，则f(t)在fm处产生的傅里叶系数为：
Figure BDA0003426795090000112
其中，T是载波的周期，k是自然数。公式(2)的计算可以从输入信号的开始就进行，在积分区间结束后，可以几乎无延迟地得到傅里叶系数。把公式(1)带入(2)得到：
Figure BDA0003426795090000113
其中，
Figure BDA0003426795090000114
Figure BDA0003426795090000115
取公式(5)的绝对值得到：
Figure BDA0003426795090000116
从公式(6)可知，当下面公式(7)中的条件满足时，|ΔFnm|＝0；
kT(fn-fm+Δfn)＝p，p＝±1，±2，… (7)
当下面公式(8)中的条件满足时，|ΔFnm|＝1；
fn-fm+Δfn＝0 (8)
具体地，所述频偏计算单元用于执行以下步骤：
利用带通滤波器把其它的载波滤除后，就剩下fm处的载波，公(9)简化为：
Fm＝Bm·ΔFmm (14)
把公(4)和(5)代入到公式(14)得到：
Figure BDA0003426795090000121
由公式(15)得到Fm的相位为：
φm＝πkT·Δfm+φm0 (16)
可以看出，由载波频偏引起的相位变化是πkT·Δfm。如果对接收到的载波以积分区间[0，kT]为单位连续地进行傅里叶系数的计算，那么相邻两次计算得到的相位变化就是πkT·Δfm。计算得到Fm的相位φm与傅里叶系数积分区间的顺次编号i是线性关系，可以表述为：
φm(i)＝πkT·Δfm·i+φn0 (17)
公式(17)是一个直线方程，为了得到准确的频偏估计，需要在连续多个傅里叶积分区间进行多次的傅里叶系数计算。
公式(17)表示的直线的斜率best可以通过最小二乘法得到，记为：
best＝πkT·Δfm (18)
这样第m个载波的频偏就是：
Figure BDA0003426795090000122
具体地，所述幅度矫正单元用于执行以下步骤：
得到了第m个载波的频偏后，把公式(18)代入到公式(6)，并令m＝n得到：
Figure BDA0003426795090000123
把公式(4)带入公式(14)并取绝对值，变换后得到：
Figure BDA0003426795090000124
其中，|Fm|根据公式(2)由傅里叶系数的计算得到；|ΔFmm|由公式(20)得到。公式(21)中的Am就是利用载波的频偏矫正得到的载波的幅度。
进一步的，所述幅度矫正单元还用于执行以下步骤：
依次选定不同载波的载波频率作为带通滤波器的中心频率，得到所有载波的频偏并矫正所有载波的幅度。
本实施例中多载波鉴频系统用于实现前述的多载波鉴频方法，因此该系统的具体实施方式可见前文中的多载波鉴频方法的实施例部分，所以，其具体实施方式可以参照相应的各个部分实施例的描述，在此不再展开介绍。
另外，由于本实施例的多载波鉴频系统用于实现前述的多载波鉴频方法，因此其作用与上述方法的作用相对应，这里不再赘述。
以上实施例仅是为充分说明本发明而所举的较佳的实施例，本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换，均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。

Claims (10)

1.多载波鉴频方法，其特征在于，包括以下步骤：
S1、通过带通滤波器接收包含多载波的调顶信号，并选定其中一个载波的载波频率作为带通滤波器的中心频率；
S2、计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数；
S3、在连续多个傅里叶积分区间执行步骤S2，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数；
S4、利用得到的多个傅里叶系数计算选定载波的频偏；
S5、利用选定载波的频偏矫正选定载波的幅度。
2.如权利要求1所述的多载波鉴频方法，其特征在于，步骤S2包括：
将调顶信号中包含的多载波表示为：
Figure FDA0003426795080000011
其中，n是载波编号；N是载波个数；An是幅度，为正的实数；φn0是载波的初始相位；fn是载波频率；Δfn是载波对于额定值的频偏；
假设fm是第m个载波频率额定值，则f(t)在fm处产生的傅里叶系数为：
Figure FDA0003426795080000012
其中，T是载波的周期，k是自然数。
3.如权利要求1所述的多载波鉴频方法，其特征在于，步骤S4包括：
所述傅里叶系数的相位与傅里叶积分区间的顺次编号之间的线性关系，表示为：
φm(i)＝πkT·Δfm·i+φn0
其中，φm为傅里叶系数的相位，i为傅里叶积分区间的顺次编号，Δfm为第m个载波的频偏，n是载波编号，T是载波的周期，k是自然数，φn0是载波的初始相位；
上述公式表示的直线的斜率best可以通过最小二乘法得到，记为：
best＝πkT·Δfm
得到第m个载波的频偏为：
Figure FDA0003426795080000021
4.如权利要求3所述的多载波鉴频方法，其特征在于，步骤S5包括：
将公式best＝πkT·Δfm代入以下公式：
Figure FDA0003426795080000022
令m＝n得到：
Figure FDA0003426795080000023
将公式
Figure FDA0003426795080000024
代入以下公式：
Fm＝Bm·ΔFmm
取绝对值，变换后得到：
Figure FDA0003426795080000025
其中，Am是利用载波的频偏矫正得到的载波的幅度。
5.如权利要求1所述的多载波鉴频方法，其特征在于，还包括以下步骤：
S6、依次选定不同载波的载波频率作为带通滤波器的中心频率，并执行步骤S2-S5，得到所有载波的频偏并矫正所有载波的幅度。
6.多载波鉴频系统，其特征在于，包括：
带通滤波器，用于接收包含多载波的调顶信号；
频率选定单元，用于选定其中一个载波的载波频率作为带通滤波器的中心频率；
傅里叶系数计算单元，用于计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数；并在连续多个傅里叶积分区间计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数，得到带通滤波器的输出信号在选定的载波频率处的多个傅里叶系数；
频偏计算单元，用于利用得到的多个傅里叶系数计算选定载波的频偏；
幅度矫正单元，用于利用选定载波的频偏矫正选定载波的幅度。
7.如权利要求6所述的多载波鉴频系统，其特征在于，所述计算带通滤波器的输出信号在选定的载波频率处的傅里叶系数，包括：
将调顶信号中包含的多载波表示为：
Figure FDA0003426795080000031
其中，n是载波编号；N是载波个数；An是幅度，为正的实数；φn0是载波的初始相位；fn是载波频率；Δfn是载波对于额定值的频偏；
假设fm是第m个载波频率额定值，则f(t)在fm处产生的傅里叶系数为：
Figure FDA0003426795080000032
其中，T是载波的周期，k是自然数。
8.如权利要求6所述的多载波鉴频系统，其特征在于，所述利用得到的多个傅里叶系数计算选定载波的频偏，包括：
所述傅里叶系数的相位与傅里叶积分区间的顺次编号之间的线性关系，表示为：
φm(i)＝πkT·Δfm·i+φn0
其中，φm为傅里叶系数的相位，i为傅里叶积分区间的顺次编号，Δfm为第m个载波的频偏，n是载波编号，T是载波的周期，k是自然数，φn0是载波的初始相位；
上述公式表示的直线的斜率best可以通过最小二乘法得到，记为：
best＝πkT·Δfm
得到第m个载波的频偏为：
Figure FDA0003426795080000041
9.如权利要求8所述的多载波鉴频系统，其特征在于，所述利用选定载波的频偏矫正选定载波的幅度，包括：
将公式best＝πkT·Δfm代入以下公式：
Figure FDA0003426795080000042
令m＝n得到：
Figure FDA0003426795080000043
将公式
Figure FDA0003426795080000044
代入以下公式：
Fm＝Bm·ΔFmm
取绝对值，变换后得到：
Figure FDA0003426795080000045
其中，Am是利用载波的频偏矫正得到的载波的幅度。
10.如权利要求6所述的多载波鉴频系统，其特征在于，所述幅度矫正单元还用于：
依次选定不同载波的载波频率作为带通滤波器的中心频率，得到所有载波的频偏并矫正所有载波的幅度。
CN202111583131.3A 2021-12-22 2021-12-22 多载波鉴频方法及系统 Active CN114422314B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111583131.3A CN114422314B (zh) 2021-12-22 2021-12-22 多载波鉴频方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111583131.3A CN114422314B (zh) 2021-12-22 2021-12-22 多载波鉴频方法及系统

Publications (2)

Publication Number Publication Date
CN114422314A true CN114422314A (zh) 2022-04-29
CN114422314B CN114422314B (zh) 2024-05-24

Family

ID=81267358

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111583131.3A Active CN114422314B (zh) 2021-12-22 2021-12-22 多载波鉴频方法及系统

Country Status (1)

Country Link
CN (1) CN114422314B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5771443A (en) * 1994-09-19 1998-06-23 Advantest Corporation Method and apparatus for measuring FM frequency deviation
JP2008175552A (ja) * 2007-01-16 2008-07-31 Japan Radio Co Ltd 圧縮係数生成器
CN102308546A (zh) * 2011-07-26 2012-01-04 华为技术有限公司 一种多载波光信号的接收方法和装置
US20130070874A1 (en) * 2010-06-17 2013-03-21 Nippon Telegraph And Telephone Corporation Frequency offset estimation apparatus, frequency offset estimation method, and reception method
CN103312655A (zh) * 2013-05-28 2013-09-18 京信通信技术(广州)有限公司 一种多载波相位优化方法及装置
CN104184471A (zh) * 2014-08-18 2014-12-03 中国航空无线电电子研究所 一种航空无线通信鉴频器及其鉴频方法
CN105158740A (zh) * 2015-08-24 2015-12-16 西安电子科技大学 基于高精度频率估计的噪声调幅干扰抑制方法
CN106230762A (zh) * 2016-07-26 2016-12-14 广州海格通信集团股份有限公司 多普勒频偏估计方法和系统
CN106464372A (zh) * 2014-05-14 2017-02-22 华为技术有限公司 发送器、接收器和频偏校正方法
US20180209851A1 (en) * 2014-06-12 2018-07-26 Colorado State University Research Foundation Apparatus and Method for Measurement of Optical Frequency Shifts
CN109067680A (zh) * 2018-09-19 2018-12-21 深圳市鼎阳科技有限公司 一种基带信号的载波频偏估计方法及其装置

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5771443A (en) * 1994-09-19 1998-06-23 Advantest Corporation Method and apparatus for measuring FM frequency deviation
JP2008175552A (ja) * 2007-01-16 2008-07-31 Japan Radio Co Ltd 圧縮係数生成器
US20130070874A1 (en) * 2010-06-17 2013-03-21 Nippon Telegraph And Telephone Corporation Frequency offset estimation apparatus, frequency offset estimation method, and reception method
CN102308546A (zh) * 2011-07-26 2012-01-04 华为技术有限公司 一种多载波光信号的接收方法和装置
CN103312655A (zh) * 2013-05-28 2013-09-18 京信通信技术(广州)有限公司 一种多载波相位优化方法及装置
CN106464372A (zh) * 2014-05-14 2017-02-22 华为技术有限公司 发送器、接收器和频偏校正方法
US20180209851A1 (en) * 2014-06-12 2018-07-26 Colorado State University Research Foundation Apparatus and Method for Measurement of Optical Frequency Shifts
CN104184471A (zh) * 2014-08-18 2014-12-03 中国航空无线电电子研究所 一种航空无线通信鉴频器及其鉴频方法
CN105158740A (zh) * 2015-08-24 2015-12-16 西安电子科技大学 基于高精度频率估计的噪声调幅干扰抑制方法
CN106230762A (zh) * 2016-07-26 2016-12-14 广州海格通信集团股份有限公司 多普勒频偏估计方法和系统
CN109067680A (zh) * 2018-09-19 2018-12-21 深圳市鼎阳科技有限公司 一种基带信号的载波频偏估计方法及其装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
谢锘;王晓君;葛建华;: "DFT测频的加窗插值算法与实现", 电子科技大学学报, no. 02 *

Also Published As

Publication number Publication date
CN114422314B (zh) 2024-05-24

Similar Documents

Publication Publication Date Title
US9698905B2 (en) Optical transmission system, optical transmission apparatus and wavelength spacing measurement apparatus
CN106160849B (zh) 功率估计装置、频谱特征监测装置和光接收机
EP0608024B1 (fr) Système de transmission avec récupération de rythme
US7539255B2 (en) Using multiple pilot signals for timing error estimation in digital subscriber line communications
CN100473063C (zh) 非线性失真均衡电路及非线性失真均衡方法
US10594406B2 (en) Pilot-aided digital signal processing for reception of coherent optical communication
CN106533573B (zh) 一种非对称孪生单边带调制、解调方法和系统
US9002215B2 (en) Spectral analysis for coherent optical receivers
JP2022000993A (ja) ディジタル信号処理装置
KR101855019B1 (ko) 신호를 송신/수신하기 위한 방법, 이에 대응하는 디바이스와 시스템
EP0762703A1 (fr) Démodulation d&#39;un signal multiporteur diminuant une distorsion blanche en fréquence
CN104365039B (zh) 在相干光通信中减少周跳
FR2721778A1 (fr) Procédé d&#39;estimation d&#39;une erreur de phase résiduelle sur les échantillons d&#39;un signal numérique démodulé, et procédé de correction correspondant.
CN113676254B (zh) 色散估计方法及装置
EP3288195A1 (en) Process for detecting a detuning at a wavelength splitter in an optical fibre transmission network
CN114422314A (zh) 多载波鉴频方法及系统
EP0696392B1 (fr) Procede et dispositif de controle automatique de gain de circuits de reception
US9548812B2 (en) Method and apparatus for sub-carrier frequency control
JP2008211760A (ja) 変調方式推定装置
CN101136707B (zh) 色散补偿装置及色散补偿控制方法
CN107070555A (zh) 处理多个光信号样本的方法、电路、光模块和光通信系统
CN111212006B (zh) 一种信号解调系统及解调方法
CN114978833A (zh) 一种基于组合调制波形的qpsk调制信号偏移补偿方法
JP4149302B2 (ja) 伝送路特性推定装置、ｏｆｄｍ信号復調装置及び回り込みキャンセラ
Yu et al. OSNR monitoring by using single sampling channel generated 2-D phase portrait

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant