CN105515666A - 一种光信号码元速率识别方法及系统 - Google Patents
一种光信号码元速率识别方法及系统 Download PDFInfo
- Publication number
- CN105515666A CN105515666A CN201510891540.8A CN201510891540A CN105515666A CN 105515666 A CN105515666 A CN 105515666A CN 201510891540 A CN201510891540 A CN 201510891540A CN 105515666 A CN105515666 A CN 105515666A
- Authority
- CN
- China
- Prior art keywords
- correlation function
- signal
- chip rate
- digital
- light signal
- 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
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/25—Arrangements specific to fibre transmission
- H04B10/2575—Radio-over-fibre, e.g. radio frequency signal modulated onto an optical carrier
- H04B10/25752—Optical arrangements for wireless networks
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Optical Communication System (AREA)
Abstract
本发明涉及光纤通信领域的一种光信号码元速率识别方法及系统。系统包括光放大器(OA),光滤波器和数字相干光接收机。通过数字相干光接收机内采样频率为Rs的模数转换器(ADC)和数字信号处理(DSP)获取信号在X、Y偏振方向的N个采样点序列,进而通过DSP进行快速傅立叶变换(FFT)得到rx(m)和ry(m),并计算rx,y(m)的自相关和互相关函数,将两种相关函数取模相加,得到混合相关函数HCF(Ω),该函数在正半轴存在一个明显的脉冲峰值,通过搜索峰值位置ΩCT即可计算信号的码元速率B=Rs(1-ΩCT/N)。本发明适用于多种调制格式信号,对色度色散(CD)、偏振模色散(PMD)和光信噪比(OSNR)劣化等损伤有较大容限,具备精度高和抗干扰强的优点。
Description
技术领域
本发明涉及光纤通信,光信号识别和数字信号处理技术领域,特别涉及光信号码元速率识别及系统。
背景技术
随着光纤通信技术的发展特别是光调制/解调技术的发展,光网络中可能同时传输多种调制格式和速率的光信号。这些光信号具有不同的谱效率,能够容忍的信道损伤和噪声大小也各不相同,因此以满足不同数据传输业务和应用场景的需求。在这种条件下不能够再预设某一个节点或光接收机处光信号的调制格式和速率是一定的。为了有效对光网络进行监测管理,对接收到的光信号进行正确的解调,首先需要识别光信号的码元速率。因此需要研发一种适用信号类型广,对各种信道损伤和噪声容限大,并且准确可靠的光信号码元识别技术与系统。
目前提出的码元速率识别方案中,有一些是基于直接探测的,利用一阶自相关方程曲线的中心宽度对OOK格式信号的码元速率进行识别,原理简单,但结果受信号损伤影响较大。另有一种方案利用二阶自相关方程曲线的周期性实现对OOK信号和NRZ-DPSK信号的码元速率识别,这种方法抗损伤性能有所提高,但是使用范围有限,不能正确识别RZ-DPSK等信号的码元速率。也有研究机构提出利用异步采样强度直方图采用人工神经网络(ANN)等方法识别光信号,但是这种方法需要对ANN进行大量的训练,并只能针对一些确定的码元速率和调制格式信号,对CD和PMD的容限较低,因此这种方法使用不便,应用范围有限,识别精度也不够高。
为解决上述方法需要研究一种适用于多种调制格式信号,对CD、PMD和OSNR劣化等损伤有较大容限,具备精度高和抗干扰强优点的码元速率识别方法。
发明内容
本发明所要解决的技术问题是提出一种基于数字相干光接收机的光信号码元速率识别方法与系统,克服现有技术所存在的识别精度低,对损伤容限小,适用的信号类型少和使用不便等缺陷。
为解决上述技术问题,本发明首先提出一种光信号码元速率识别方法,其特征在于,包括以下步骤:
待识别光信号首先进入光放大器,所述光放大器将光信号放大至一定功率后经光滤波器输入数字相干光接收机;
所述数字相干光接收机通过模数转换以大于信号实际码元速率的采样率Rs分别获得信号X、Y偏振方向的长度为N的采样序列,记为rx(n)和ry(n),并将其输入数字信号处理模块,所述N的取值为512,1024,2048,4096,…,N越大速率估计精度越高,但计算量越大;
所述数字信号处理模块首先通过快速傅立叶变换计算采样序列rx(n)、ry(n)的频谱rx(m)、ry(m),然后计算rx(m)的自相关函数ACFx(Ω),Ω代表频谱位移大小,计算rx(m)与ry(m)的互相关函数XCFx,y(Ω);将自相关函数、互相关函数取模相加得到混合相关函数HCF(Ω),该混合相关函数在正半轴上存在一个明显的脉冲峰值,搜索该混合函数在正半轴的脉冲峰值位置ΩCT,计算得到信号的码元速率B=Rs(1-ΩCT/N)。
作为一种替代方案,所述HCF(Ω)在负半轴也存在同样的一个脉冲峰值,利用这个负半轴的脉冲峰值位置代替码元速率计算公式中的ΩCT,计算得到光信号符号速率。
作为又一种替代方案,计算ry(m)的自相关函数ACFy(Ω),用自相关函数ACFy(Ω)替代所述ACFx(Ω)计算混合相关函数HCF(Ω)。
在偏振模色散PMD损伤不大的情况下,所述HCF(Ω)能够用ACFx(Ω)或ACFy(Ω)的模值取代。
本发明同时还提出了一种实现所述光信号码元速率识别方法的系统,其特征在于,包括光放大器OA,光滤波器,数字相干光接收机和数字信号处理模块;
所述OA,用于放大光信号至一定功率,然后将光信号输入所述光滤波器;
所述光滤波器,用于滤除输入光信号中的信道外噪声;
所述数字相干接收机,用于接收所述光滤波器的输出信号,还包括模数转换器,所述模数转换器用于获得待测光信号在X、Y偏振方向的N个采样点序列rx(n)和ry(n);
所述数字信号处理模块,用于计算所述采样序列rx(n)、ry(n)的FFT频谱rx(m)、ry(m),然后计算rx(m)的自相关函数,rx(m)和ry(m)的互相关函数;将两相关函数取模相加得到混合相关函数,最后通过搜索混合相关函数的脉冲峰值位置计算出信号的码元速率。
优选的,所述数字信号处理模块集成在所述数字相干光接收机内部。
本发明无需任何先验信息,适用于多种调制格式信号,对CD、PMD和OSNR劣化等损伤有较大容限,具备精度高和抗干扰强的优点。
附图说明
下面结合附图和具体实施方式对本发明的技术方案作进一步具体说明。
图1为本发明具体实施的码元速率识别系统结构示意图。
图2为10GBaud偏振复用正交相移键控(DP-QPSK)信号的HCF(Ω)曲线图。
图3分别为10GBaudNRZ(非归零)DP-QPSK信号在PMD造成的差分群速延迟(DGD)和偏振角变化时位于ΩCT处归一化脉冲峰值强度的仿真结果。
图4为10GBaudNRZ-DP-QPSK信号在不同损伤下HCF(Ω)脉冲峰值相对于其均值的峰均值比(PAR)变化图。
图5为10GBaudNRZ-DP-QPSK信号的HCF(Ω)在不同采样速率和参数运算点数N时的ΩCT处归一化脉冲峰值强度的变化图。
图6左图为NRZ、33%-RZ和67%-RZ的10GBaudDP-QPSK信号的码元速率识别误差分布直方图;右图为上述样本的码元速率识别结果误差的标准差随着运算点数N的变化图。
图7左图为NRZ、33%-RZ和67%-RZ的10GBaudDP-16QAM信号的码元速率识别误差分布直方图,右图为码元速率识别结果误差的标准差随着运算点数N的变化图。
具体实施方式
如图1所示的码元速率识别系统包括:光放大器(OA)1,光滤波器2,数字相干光接收机3。数字相干光接收机3内部组成包括光学前端4,模数转换器(ADC)5和数字信号处理器(DSP)6。
下面进一步描述光信号码元速率识别的具体实施方式和系统的工作流程如下:
待识别光信号首先进入光放大器1,光放大器1将光信号放大至一定功率后经光滤波器2输入数字相干光接收机3的光学前端4;
数字相干光接收机通过其内部的采样速率为Rs的ADC模块5以大于信号实际码元速率的采样率分别获得信号X、Y偏振方向的长度为N的采样序列,记为rx(n)和ry(n),并将其输入到其内部的DSP模块6;
DSP模块首先计算两路采样序列rx(n)和ry(n)的FFT频谱rx(m)和ry(m)
其中p=x,y代表X和Y偏振方向信号。然后分别计算rx(m)的自相关函数ACFx,rx(m)和ry(m)的互相关函数XCFx,y
ACFx(Ω)=∑rx *(n)rx(n+Ω),(2)
XCFx,y(Ω)=∑rx *(n)ry(n+Ω)(3)
其中Ω代表频谱位移点数,*表示取复共轭。将两种相关函数取模相加得到混合相关函数HCF
HCF(Ω)=|ACFx(Ω)|+|XCFx,y(Ω)|,(4)
HCF(Ω)存在一个明显的时钟脉冲,脉冲位置记为ΩCT,在相干接收机采样速率不变的情况下,ΩCT与码元速率成正比,有
Ω=ΩCT=±N·(1-B/Rs),(5)
HCF(Ω)在ΩCT的脉冲峰值强度记为HCTM,
HCTM=HCF(ΩCT)=|ACFx(ΩCT)|+|XCFx,y(ΩCT)|.(6)
HCTM结合了|ACFx(ΩCT)|和|XCFx,y(ΩCT)|的优势,在任何PMD损伤影响下,HCTM均可以保证较大的强度,且是HCF(Ω)序列中的最大值的位置,即可根据公式(5)推算出待测光信号的码元速率B=Rs(1-ΩCT/N)。
作为一种等效方法HCF(Ω)在负半轴也存在同样的一个脉冲峰值,利用这个脉冲峰值位置也可以确定光信号符号速率。作为一种等效方法ACFx(Ω)可以用ry(m)自相关函数ACFy(Ω)替代。此外在PMD损伤不大的情况下作为一种等效方法HCF(Ω)可用ACFx(Ω)或ACFy(Ω)的模值取代。
图2为10GBaud偏振复用正交相移键控(DP-QPSK)信号的HCF(Ω)曲线,ADC采样率为2倍码元速率,参与运算点数N为2048,无CD和PMD损伤。
图3分别为10GBaudNRZ(非归零)DP-QPSK信号在PMD造成的差分群速延迟(DGD)和偏振角变化时位于ΩCT处归一化脉冲峰值强度的仿真结果。可以看出PMD对HCTM几乎没有影响。
图4为10GBaudNRZ-DP-QPSK信号在不同损伤下HCF(Ω)脉冲峰值相对于其均值的峰均值比(PAR)变化,可以看出在较强的损伤下PAR仍然维持较高的水平,说些对信号损伤容限较大。
图5为10GBaudNRZ-DP-QPSK信号的HCF(Ω)在不同采样速率和参数运算点数N时的ΩCT处归一化脉冲峰值强度的变化。
图6左图为NRZ、33%-RZ和67%-RZ的10GBaudDP-QPSK信号的码元速率识别误差分布直方图,由32220个样本的计算结果构成。运算点数N设为1024。色散(CD)范围0-100ps/nm,光信噪比(OSNR)范围15-35dB,偏振角θ从0变到90°,DGD从0到100ps,信号采样率为1.3-3倍码元速率。右图为上述样本的码元速率识别结果误差的标准差随着运算点数N的变化,可以看到当N>768时,误差的标准差已经小于信号码元速率的0.02%。
图7左图为NRZ、33%-RZ和67%-RZ的10GBaudDP-16QAM信号的码元速率识别误差分布直方图,由32220个样本的计算结果构成,样本构成与图6相同。右图为码元速率识别结果误差的标准差随着运算点数N的变化,可以看到当N>768时,误差的标准差已经小于信号码元速率的0.02%。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照较佳实施示例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种光信号码元速率识别方法,其特征在于,包括以下步骤:
待识别光信号首先进入光放大器,所述光放大器将光信号放大至一定功率后经光滤波器输入数字相干光接收机;
所述数字相干光接收机通过模数转换以大于信号实际码元速率的采样率Rs分别获得信号X、Y偏振方向的长度为N的采样序列,记为rx(n)和ry(n),并将其输入数字信号处理模块,所述N的取值为512,1024,2048,4096,…,N越大速率估计精度越高,但计算量越大;
所述数字信号处理模块首先通过快速傅立叶变换计算采样序列rx(n)、ry(n)的频谱rx(m)、ry(m),然后计算rx(m)的自相关函数ACFx(Ω),Ω代表频谱位移大小,计算rx(m)与ry(m)的互相关函数XCFx,y(Ω);将自相关函数、互相关函数取模相加得到混合相关函数HCF(Ω),该混合相关函数在正半轴上存在一个明显的脉冲峰值,搜索该混合函数在正半轴的脉冲峰值位置ΩCT,计算得到信号的码元速率B=Rs(1-ΩCT/N)。
2.根据权利要求1所述的光信号码元速率识别方法,其特征在于,所述HCF(Ω)在负半轴也存在同样的一个脉冲峰值,利用这个负半轴的脉冲峰值位置代替码元速率计算公式中的ΩCT,计算得到光信号符号速率。
3.根据权利要求1或2所述的光信号码元速率识别方法,其特征在于,计算ry(m)的自相关函数ACFy(Ω),用自相关函数ACFy(Ω)替代所述ACFx(Ω)计算混合相关函数HCF(Ω)。
4.根据权利要求3所述的光信号码元速率识别方法,其特征在于,在偏振模色散PMD损伤不大的情况下,所述HCF(Ω)能够用ACFx(Ω)或ACFy(Ω)的模值取代。
5.一种实现权利要求1所述光信号码元速率识别方法的系统,其特征在于,包括光放大器OA,光滤波器,数字相干光接收机和数字信号处理模块;
所述OA,用于放大光信号至一定功率,然后将光信号输入所述光滤波器;
所述光滤波器,用于滤除输入光信号中的信道外噪声;
所述数字相干接收机,用于接收所述光滤波器的输出信号,还包括模数转换器,所述模数转换器用于获得待测光信号在X、Y偏振方向的N个采样点序列rx(n)和ry(n);
所述数字信号处理模块,用于计算所述采样序列rx(n)、ry(n)的FFT频谱rx(m)、ry(m),然后计算rx(m)的自相关函数,rx(m)和ry(m)的互相关函数;将两相关函数取模相加得到混合相关函数,最后通过搜索混合相关函数的脉冲峰值位置计算出信号的码元速率。
6.根据权利要求5所述的利用权利要求1或3所述光信号码元速率识别方法的系统,其特征在于,所述数字信号处理模块集成在所述数字相干光接收机内部。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510891540.8A CN105515666B (zh) | 2015-12-04 | 2015-12-04 | 一种光信号码元速率识别方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510891540.8A CN105515666B (zh) | 2015-12-04 | 2015-12-04 | 一种光信号码元速率识别方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105515666A true CN105515666A (zh) | 2016-04-20 |
CN105515666B CN105515666B (zh) | 2018-06-19 |
Family
ID=55723371
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510891540.8A Expired - Fee Related CN105515666B (zh) | 2015-12-04 | 2015-12-04 | 一种光信号码元速率识别方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105515666B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106100930A (zh) * | 2016-06-22 | 2016-11-09 | 中国人民解放军信息工程大学 | 基于精度递进定时估计的波特率细微特征估计方法 |
CN114257303A (zh) * | 2021-12-10 | 2022-03-29 | 武汉邮电科学研究院有限公司 | 一种抑制光纤通信非线性效应的方法和系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103414517A (zh) * | 2013-07-22 | 2013-11-27 | 华中科技大学 | 一种光信号速率识别方法及系统 |
US20140064402A1 (en) * | 2012-09-06 | 2014-03-06 | Electronics And Telecommunications Research Institute | Apparatus and method for modulation classification in wireless communication system |
CN104683027A (zh) * | 2015-02-15 | 2015-06-03 | 华中科技大学 | 调制格式无关光信号速率识别方法及系统 |
-
2015
- 2015-12-04 CN CN201510891540.8A patent/CN105515666B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140064402A1 (en) * | 2012-09-06 | 2014-03-06 | Electronics And Telecommunications Research Institute | Apparatus and method for modulation classification in wireless communication system |
CN103414517A (zh) * | 2013-07-22 | 2013-11-27 | 华中科技大学 | 一种光信号速率识别方法及系统 |
CN104683027A (zh) * | 2015-02-15 | 2015-06-03 | 华中科技大学 | 调制格式无关光信号速率识别方法及系统 |
Non-Patent Citations (3)
Title |
---|
JIN SHANG,ETAL.: "Symbol Rate Identification Using Asynchronous Delayed Sampling", 《LASERS AND ELECTRO-OPTICS (CLEO),2015 CONFERENCE ON》 * |
YUAN ZHOU,ETAL.: "Bit-Rate Identification Using Asynchronous Delayed Sampling", 《IEEE PHOTONICS TECHNOLOGY LETTERS》 * |
彭勃,熊刚: "一种基于二次谱线生成的码速率识别方法", 《电讯技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106100930A (zh) * | 2016-06-22 | 2016-11-09 | 中国人民解放军信息工程大学 | 基于精度递进定时估计的波特率细微特征估计方法 |
CN106100930B (zh) * | 2016-06-22 | 2019-10-29 | 中国人民解放军信息工程大学 | 基于精度递进定时估计的波特率细微特征估计方法 |
CN114257303A (zh) * | 2021-12-10 | 2022-03-29 | 武汉邮电科学研究院有限公司 | 一种抑制光纤通信非线性效应的方法和系统 |
CN114257303B (zh) * | 2021-12-10 | 2023-10-03 | 武汉邮电科学研究院有限公司 | 一种抑制光纤通信非线性效应的方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105515666B (zh) | 2018-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103312645B (zh) | 一种数字相干光通信系统的损伤自适应补偿方法及系统 | |
US9515742B2 (en) | System and methods for adaptive equalization for optical modulation formats | |
CN102231648B (zh) | 基于单载波和多载波的混合传输系统 | |
US11108490B2 (en) | Systems and methods for carrier phase recovery | |
US20200106523A1 (en) | Fiber nonlinear noise monitoring method and apparatus and optical receiver | |
CN103281120B (zh) | 一种光信号调制格式识别方法及系统 | |
CN102983910A (zh) | 相干光通信系统中色散和非线性补偿方法及系统 | |
CN107634814A (zh) | 一种自零差检测模分复用系统中载波路串扰的消除方法 | |
CN110071738A (zh) | 基于多伪码支路接收机的扩频信号解扩跟踪方法 | |
CN107078982A (zh) | 数字相干光接收器的色散估计 | |
CA2929182C (en) | Spreading signal generating method, generating device, receiving method and receiving device | |
CN104683027B (zh) | 调制格式无关的光信号速率识别方法及系统 | |
CN103944638A (zh) | 基于数字非线性处理的光信号调制格式识别方法及系统 | |
CN102377485A (zh) | 一种解调光调顶信号的方法和装置 | |
CN103178900A (zh) | 基于支持向量机的相干光通信系统非线性损伤消除装置 | |
WO2016173324A1 (zh) | 一种光信噪比监测的方法及装置 | |
EP2465219A1 (en) | Filter device | |
CN104868969B (zh) | 一种基于斯托克斯分析的非正交偏振复用信号传输方法 | |
Saif et al. | Separability of histogram based features for optical performance monitoring: An investigation using t-SNE technique | |
CN109347776B (zh) | 一种差分相位幅度比的光通信信号调制格式识别方法 | |
Fougstedt et al. | ASIC implementation of time-domain digital backpropagation with deep-learned chromatic dispersion filters | |
CN105515666A (zh) | 一种光信号码元速率识别方法及系统 | |
US9571224B2 (en) | Transmission device, receiving device, and communication method | |
CN103051375A (zh) | 无线激光通信外差检测系统及其检测方法 | |
CN104297766B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180619 Termination date: 20181204 |