CN107302507A - 基于聚类算法的相干光通信盲均衡方法 - Google Patents
基于聚类算法的相干光通信盲均衡方法 Download PDFInfo
- Publication number
- CN107302507A CN107302507A CN201710358012.5A CN201710358012A CN107302507A CN 107302507 A CN107302507 A CN 107302507A CN 201710358012 A CN201710358012 A CN 201710358012A CN 107302507 A CN107302507 A CN 107302507A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- msub
- signal
- polarization
- 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
- 230000001427 coherent effect Effects 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000004891 communication Methods 0.000 title claims abstract description 18
- 230000010287 polarization Effects 0.000 claims abstract description 44
- 230000003287 optical effect Effects 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims description 20
- 239000000654 additive Substances 0.000 claims description 10
- 230000000996 additive effect Effects 0.000 claims description 10
- 238000013461 design Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000007493 shaping process Methods 0.000 claims description 3
- 230000003111 delayed effect Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 230000006735 deficit Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 102100029764 DNA-directed DNA/RNA polymerase mu Human genes 0.000 description 1
- 101150095211 POLM gene Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03891—Spatial equalizers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J14/00—Optical multiplex systems
- H04J14/06—Polarisation multiplex systems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/02—Arrangements for detecting or preventing errors in the information received by diversity reception
- H04L1/06—Arrangements for detecting or preventing errors in the information received by diversity reception using space diversity
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03082—Theoretical aspects of adaptive time domain methods
- H04L25/03089—Theory of blind algorithms, recursive or not
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
- H04L25/03082—Theoretical aspects of adaptive time domain methods
- H04L25/03095—Theory of fractional equalisers, recursive or not
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L2025/0335—Arrangements for removing intersymbol interference characterised by the type of transmission
- H04L2025/03375—Passband transmission
- H04L2025/0342—QAM
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Theoretical Computer Science (AREA)
- Optical Communication System (AREA)
Abstract
本发明公开了一种基于聚类算法的相干光通信盲均衡方法,通过构造基于聚类算法的代价函数,将相干光通信系统的盲均衡问题转化为一个最优化求解问题,通过将代价函数的概率密度函数对数最大化求解该优化问题。与传统的基于聚类的算法相比较,该算法不仅极大地降低了计算复杂度,而且在算法收敛过程中,步长和方差能够自适应调整。该方法适用于偏振多路复用(PolMux)相干光通信系统盲均衡问题。
Description
技术领域
本发明涉及偏振多路复用相干光通信的信号处理技术领域,特别是一种基于聚类算法的相干光通信盲均衡方法。该方法适用于高阶正交幅度调制(QAM)系统,特别是在算法收敛过程中无法选择合适的方差和步长的场合,是一种高鲁棒性和低算法复杂度的光通信系统盲均衡方法。
背景技术
目前,结合了可消除传播损伤的数字信号处理(DSP)技术的数字相干接收器,已彻底改变了光通信系统的设计方法。而且随着集成技术的快速发展,DSP算法在超高数据速率相干光通信系统中扮演着越来越重要的角色。盲均衡算法因其无需传输导频(训练)数据,大大节约了有限的通信系统带宽,得到了广大研究人员的重视。另一方面,偏振复用(PolMux)相干系统因其能使系统频谱效率加倍,而被广泛应用于相干光通信系统中;基于此技术近年也出现了许多结合了正交相移键控(QPSK)或QAM技术的PolMu光相干系统应用实例。
聚类技术是一种无监督的分类技术,己广泛应用到了信号模式识别、信息检索、图像分割等多个领域,并取得了显著的成果。聚类分析通过定义相似性度量将样本划分成多个不相交的类,其目标是使得同一类中的样本尽可能紧凑,不同类中的样本尽可能分散,以便于发现数据之间隐含的结构信息。将聚类技术应用于均衡算法中,对信号星座图进行处理,降低了均衡算法的计算复杂度,从而整体提高接收端均衡器的性能。然而,当传统的基于聚类的盲均衡算法被应用于高阶QAM系统时,由于其计算复杂度过高,而无法在实际场合应用;此外,传统聚类均衡方法的收敛速度通常很慢。为了克服这些问题,业内研究人员提出了多层聚类算法(Multi-stage Clustering Algorithm,MSA),并将该分层聚类方法运用到盲均衡中。对于高阶QAM系统来说,多层是必须的,所需的聚类层级数多少取决于相关调制的阶数大小。但是该算法仍然存在两个问题:其一,选择两层之间转换的优化时间的选择非常困难。如果选择的时间太短,不能保证收敛,反之,如果时间太长了,收敛速度则变得很慢。其二,该算法对于高斯聚类噪声方差ρ和算法步长μ非常敏感;所选择的ρ值必须确保各类之间具有一个合适的间隔。如果ρ值太大,则无法实现期望的分离度;如果ρ值太小,则不能严格控制各类的大小,导致算法失效。步长μ和ρ情况相同,因为算法需要较快的收敛速度和低输出误差水平。
发明内容
本发明所要解决的技术问题是,针对现有技术不足,提供一种基于聚类算法的相干光通信盲均衡方法。
为解决上述技术问题,本发明所采用的技术方案是:一种基于聚类算法的相干光通信盲均衡方法,包括以下步骤:
1)已调符号序列{s(k)}通过一个成型滤波器g(t)后的基带信号在光信道中传播,受到加性噪声的影响输出连续时间矩阵信号y(t):
其中:y(t)是二进制连续时间矩阵信号;y(t)=[yp(t),yq(t)]T、x(t)=[xp(t),xq(t)]T、n(t)=[np(t),nq(t)]T分别表示接收信号、源信号和加性高斯噪声;p和q分别表示X轴方向偏振和Y轴方向偏振;yp(t),yq(t)分别表示X轴方向偏振接收信号和Y轴方向偏振接收信号;xp(t),xq(t)分别表示X轴方向偏振源信号和Y轴方向偏振源信号;np(t),nq(t)分别表示X轴方向偏振加性高斯噪声和Y轴方向偏振加性高斯噪声;上标T表示矩阵或向量的转置运算;c(t)表示信道脉冲响应;t0为任意时间延迟;Ts是符号周期;
2)分别在X轴方向偏振p和Y轴方向偏振q上,对信号y(t)以Ts/2空间的时间间隔进行采样,得到二进制离散信号 k表示第k个码元,i表示均衡器第i个抽头;
3)利用所述二进制离散信号,在X轴偏振方向,得到抽头间隔为Ts/2的分数间隔均衡器FSE的输出信号为:
其中:L是分数间隔均衡器FSE的长度;Wp,q是输入偏振p和输出偏振q之间的FSE,Wp,p(i)是Wp,q表示均衡器第i个抽头系数;在FSE输出样本中进行2倍抽取,创建TS-间隔输出序列,输出序列矩阵为:z(n*)=wTyL(n*);其中:z(n*)=[zp(n*),zq(q)]T;yL(n*)=[y(n*)T,y(n*-1)T,...y(n*-L+1)T]T是离散时间矩阵信号,y(n*)是y(t)的离散时间信号表达式,n*≥0的整数;w=[w0,w1,...wL-1]T是均衡器向量,w0,w1,...wL1为该均衡器向量的元素,即抽头系数;
4)将高斯白噪声等效为v(n*),则输出序列矩阵z(n*)表示为:表示延迟了d个单元的原始发送信号;
5)设计以下代价函数Jcluster(w,z(n*)):Jcluster(w,z(n*))=E[ρlog(p(w,z(n*)))];其中,pql为yql的先验概率,yql为q方向上分数间隔均衡器FSE输出的第l个值,即rp(k)的第l个值;ρ为高斯聚类噪声方差;E[·]表示期望运算,M为调制信号的阶数;K为一个大于1的正整数;
6)利用下式更新分数间隔均衡器FSE的权值:
则第n*次迭代后的分数间隔均衡器FSE权值为:w(n*)=wR(n)*+jwI(n);*
其中:
其中:μ为步长,为求偏导数运算符;下标R和I分别表示实部和虚部,ρR表示高斯聚类噪声方差ρ的实部;ρI表示高斯聚类噪声方差ρ的虚部; |.|表示绝对值运算;j为虚数单位;
参数ρ和μ的更新方程如下:
其中判决误差e(n*)为:
e(n*)=d(n*)-z(n*)
其中:d(n*)为判决值;sgn(·)为符号函数;和分别代表实部和虚部;利用泰勒公式将ρ和μ的更新方程改写为:
μ(n*+1)=μ(n*)(1+β+β2+…)≈μ(n*)(1+β)
ρ(n*+1)=ρ(n*)(1+β+β2+…)≈ρ(n*)(1+β)
这里β=-γμ(n*)e(n*)2,n*表示第n*次迭代;
7)将经步骤6)更新后的权值代入步骤5))中的代价函数Jcluster(w,z(n*)),直至该代价函数值不再变化,退出。M为4,8,16,32,64…。
与现有技术相比,本发明所具有的有益效果为:本发明在聚类算法收敛的过程中能够自适应调整的全新的盲均衡方法,增强了算法鲁棒性同时提高了收敛速度,是一种高鲁棒性和低算法复杂度的光通信系统盲均衡方法。
附图说明
图1是本发明的信号模型图;
图2是16QAM的判决区域;
图3是本发明所采用的PolMux-16QAM相干系统仿真模型;
图4是在PolMux-16QAM相干系统中,进行均衡前信号以及分别应用了CMA、MMA、基于聚类的多层级算法(M-CLUSTER)和本发明所提出的基于聚类的算法(P-CLUSTER)进行均衡后的信号星座图;
图5是在PolMux-16QAM相干系统中,CMA、MMA、M-CLUSTER和P-CLUSTER的比特误码率(BER)与光信噪比(OSNR)的关系曲线图;
图6是在PolMux-16QAM相干系统中,CMA、MMA、M-CLUSTER和P-CLUSTER的BER与迭代次数的关系曲线图;
图7是在PolMux-16QAM相干系统中,CMA、MMA、M-CLUSTER和P-CLUSTER的BER与旋转速度的关系曲线图;
图8是本发明在X轴偏振方向,判决误差e、步长μ和方差ρ与迭代次数的关系图。
具体实施方式
本发明只涉及线性传播损伤,如色度色散(CD)和偏振膜色散(PMD)。同时,本发明的讨论仅限于光纤产生的线性损伤,相位恢复不属本发明的讨论范围。对于传输速率超过100Gps的系统,符号周期远小于信道的相干时间,因此可以假设光传播信道的所有过程都是线性时不变的。
本发明提供了一种基于聚类算法的偏振复用相干光通信系统盲均衡方法,包括以下步骤:
(1)本发明的信号模型如图1所示,首先,已调符号序列{s(k)}(k表示第k个码元)通过一个成型滤波器g(t)(t表示时间)后的基带信号在光信道中传播,受到加性噪声的影响输出连续时间矩阵信号表示为:
其中:y(t)=[yp(t),yq(t)]T,x(t)=[xp(t),xq(t)]T,n(t)=[np(t),nq(t)]T分别表示接收信号,源信号和加性高斯噪声,p和q分别表示X轴和Y轴方向偏振;上标T表示矩阵或向量的转置运算;c(t)表示信道脉冲响应;t0为任意时间延迟;Ts是符号周期。
(2)分别在偏振p和q上,对接收信号y(t)以Ts/2空间的时间间隔进行采样,得到二进制离散信号:
其中:L是蝶形FSE的长度;Wp,q是输入偏振p和输出偏振q之间FSE,yp和yq分别表示接收信号在偏振p和q上的接收序列的单个接收符号,k表示第k个码元。
(4)最后在FSE输出样本中进行2倍抽取,创建TS-间隔输出序列。为表示方便,本发明假设FSE具有偶数个抽头,并假设在2倍采样中仅保留奇数输出样本。输出序列矩阵为:
其中:L是蝶形FSE的长度;Wp,q是输入偏振p和输出偏振q之间FSE,Wp,p是输入偏振p和输出偏振q之间FSE,k表示第k个码元,i表示第i个均衡器抽头。
(4)最后在FSE输出样本中进行2倍抽取,创建TS-间隔输出序列。为表示方便,本发明假设FSE具有偶数个抽头,并假设在2倍采样中仅保留奇数输出样本。X偏振方向上的输出序列矩阵为:
这里wp,p,wp,q分别表示蝶形均衡器p和p之间的权值以及p和q之间的权值,Wp,p(i)是Wp,q表示第i个均衡器抽头的系数;
记:蝶形均衡器权值系数
wp,p=[wp,p(2i),wp,p(2i+1)]
wp,q=[wp,q(2i),wp,q(2i+1)]
则输出离散信号
可得:
将上式改写为矩阵形式:
其中:均衡器权值系数
wp=[wp,p(0),...,wp,p(L-1),wp,q(0),...,wp,q(L-1)]T
离散输出时间矩阵信号
yL(n)=[yp(n)T,...,yp(n-L+1)T,yq(n)T,...,yq(n-L+1)T]T
(5)将高斯白噪声等效为v(n),均衡器的输出可以表示为:
z(n)=s(n-d)+v(n)
其中d为整数。
假设定义M-QAM符号星座图为:
sil=(2i-Q-1)+j(2l-Q-1),1≤i,l≤Q
其中:M为调制信号的阶数,一般为4,8,16,32,64…,K为一个大于1的正整数。接着将均衡器输出用M-“高斯簇”进行近似建模:zil=sil,1≤i,l≤Q,且近似协方差矩阵为
其中:E[·]表示期望运算;ρ为高斯聚类噪声方差,下标R和I分别表示实部和虚部。
(6)为了降低计算复杂度,如图2所示,将复平面分成M/4个常规区域,每个区域包含4个符号点。如果均衡器的输出z(n)在其中一个区域内,将yql的先验概率记为Pql,z(n)的后验概率密度函数为:
其中:为p的估计值。
为了进一步降低计算复杂度,假设两个一维高斯分布分别对应于MMA的实部与虚部,这两个一维高斯分布可用一个二维高斯分布来代替。那么,两个独立的代价函数定义如下:
JR(wR,zR(n))=E[ρRlog(pR(wR,zR(n)))]
JI(wI,zI(n))=E[ρIlog(pI(wI,zI(n)))]
这里:下标R和I分别表示实部和虚部;pR(wR,zR(n))和pI(wI,zI(n))分别为:
其中|.|表示绝对值运算。
均衡器的权值策略调整如下:
其中:μ为步长,为求偏导数运算符。
那么,第n次迭代后的复均衡器权值为:
w(n)=wR(n)+jwI(n)
式中:j为虚数单位。
参数ρ和μ的更新方程如下:
其中判决误差e(n)为:
e(n)=d(n)-z(n)
其中:d(n)为判决值;sgn(·)为符号函数;和分别代表实部和虚部。
利用泰勒公式
将p和μ的更新方程改写为:
μ(n+1)=μ(n)(1+β+β2+…)≈μ(n)(1+β)
ρ(n+1)=ρ(n)(1+β+β2+…)≈ρ(n)(1+β)
这里β=-γμ(n)e(n)2,n表示第n次迭代。
本发明效果实例:
本发明用仿真结果说明和验证理论的发展。图3是本发明在符号速率为14GBaud的PolMux-16QAM相干系统中所采用的仿真模型。伪随机序列的长度为16384。由于本发明不讨论相位恢复,激光器差拍线宽设置为0。均衡器长度L=9。
仿真中的参数设置如下:CMA和MMA的步长为1×10-5。对于M-CLUSTER,两阶段的步长都采用1×10-5。第一阶段使用3000个样本,ρ1=1.6,第二级中ρ2=0.6。对于P-CLUSTER算法,步长μ和协方差ρ的初值分别为4×10-5和3。
图4为进行均衡前信号以及分别应用了CMA、MMA、M-CLUSTER和P-CLUSTER进行均衡后的信号星座图。传播信道设置如下:CD=1000ps/nm(皮秒/纳米),DGD延迟τDGD=50ps,偏振旋转角度θ=π/4,OSNR=20dB(分贝)。可见,P-CLUSTER的信号星座图比CMA和MMA的更集中,更清晰,其性能与M-CLUSTER相仿。图5是CMA、MMA、M-CLUSTER和P-CLUSTER的比特误码率(BER)与光信噪比(OSNR)的关系曲线图。由图可知,P-CLUSTER比CMA和MMA更高效。例如,当BER=10-3时,P-CLUSTER的OSNR性能比MMA和CMA分别少3dB和5dB。随着OSNR的增加,P-CLUSTER的BER比M-CLUSTER的更低。图6是CMA、MMA、M-CLUSTER和P-CLUSTER的BER与迭代次数的关系曲线图。由图可知,与其它三种算法比较,P-CLUSTER算法的收敛速度最快,并且在稳定阶段有较小的稳定BER。注意,从2000到3000随着迭代增加,M-CLUSTER算法的BER没有变化,这正如前面所述,选择第一级所需要的最优迭代次数是困难的,因此,第一级需要更多次的迭代来保证M-CLUSTER算法收敛。为了了解P-CLUSTER算法的动态特性,假设剩余CD为0,采用琼斯矩阵模拟一个无限偏振旋转:
其中:ω是旋转角速率;cos(·)表示余弦运算;sin(·)表示正弦运算。图7展示了不同均衡器的跟踪性能。由图可知三个均衡器的性能是相似的。均衡器能够跟踪角频率的范围为8×105rad/s(弧度/秒)。当ω<8×105rad/s,P-CLUSTER的性能远优于CMA和MMA,比M-CLUSTER稍好。图8是本发明在X轴偏振方向,判决误差e、步长μ和方差ρ与迭代次数的关系图。图8显示步长μ,方差ρ和判决误差e的一致性,说明了P-CLUSTER算法的合理性。
Claims (2)
1.一种基于聚类算法的相干光通信盲均衡方法,其特征在于,包括以下步骤:
1)已调符号序列{s(k)}通过一个成型滤波器g(t)后的基带信号在光信道中传播,受到加性噪声的影响输出连续时间矩阵信号y(t):
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mo>-</mo>
<mi>&infin;</mi>
</mrow>
<mi>&infin;</mi>
</munderover>
<mi>c</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<msub>
<mi>nT</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中:y(t)是二进制连续时间矩阵信号;y(t)=[yp(t),yq(t)]T、x(t)=[xp(t),xq(t)]T、n(t)=[np(t),nq(t)]T分别表示接收信号、源信号和加性高斯噪声;p和q分别表示X轴方向偏振和Y轴方向偏振;yp(t),yq(t)分别表示X轴方向偏振接收信号和Y轴方向偏振接收信号;xp(t),xq(t)分别表示X轴方向偏振源信号和Y轴方向偏振源信号;np(t),nq(t)分别表示X轴方向偏振加性高斯噪声和Y轴方向偏振加性高斯噪声;上标T表示矩阵或向量的转置运算;c(t)表示信道脉冲响应;t0为任意时间延迟;Ts是符号周期;
2)分别在X轴方向偏振p和Y轴方向偏振q上,对信号y(t)以Ts/2空间的时间间隔进行采样,得到二进制离散信号 k表示第k个码元,i表示第i个抽头;
3)利用所述二进制离散信号,在X轴偏振方向,得到抽头间隔为Ts/2的分数间隔均衡器FSE的输出信号为:
<mrow>
<msub>
<mi>r</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mn>2</mn>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>&lsqb;</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mo>)</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>q</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>q</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mi>i</mi>
</mrow>
<mo>)</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>s</mi>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>;</mo>
</mrow>
其中:L是分数间隔均衡器FSE的长度;Wp,q是输入偏振p和输出偏振q之间的FSE,Wp,p(i)是Wp,q的第i个抽头系数;
在FSE输出样本中进行2倍抽取,创建Ts-间隔输出序列,输出序列矩阵为:z(n*)=wTyL(n*);其中:z(n*)=[zp(n*),zq(q)]T;yL(n*)=[y(n*)T,y(n*-1)T,...y(n*-L+1)T]T是离散时间矩阵信号,y(n*)是y(t)的离散时间信号表达式,n*是大于等于零的整数;w=[w0,w1,...wL-1]T是均衡器向量,w0,w1,...wL1为该均衡器向量的元素,即抽头系数;
4)将高斯白噪声等效为v(n*),则输出序列矩阵z(n*)表示为:s(n*-d)表示延迟了d个单元的原始发送信号;
5)设计以下代价函数Jcluster(w,z(n*)):Jcluster(w,z(n*))=E[ρlog(p(w,z(n*)))];其中,pql为yql的先验概率,yql为q方向上分数间隔均衡器FSE输出的第l个值,即rp(k)的第l个值;ρ为高斯聚类噪声方差;E[·]表示期望运算,M为调制信号的阶数;K为一个大于1的正整数;
6)利用下式更新分数间隔均衡器FSE的权值:
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>w</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&mu;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>J</mi>
<mrow>
<mi>c</mi>
<mi>l</mi>
<mi>u</mi>
<mi>s</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
<mo>,</mo>
<mi>z</mi>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>w</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
则第n*次迭代后的分数间隔均衡器FSE权值为:w(n*)=wR(n*)+jwI(n)*,j为虚数单位;
其中:
<mrow>
<msub>
<mi>w</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&mu;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>J</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>w</mi>
<mi>R</mi>
</msub>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>R</mi>
</msub>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>w</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>w</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>w</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&mu;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>J</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>w</mi>
<mi>I</mi>
</msub>
<mo>,</mo>
<msub>
<mi>z</mi>
<mi>I</mi>
</msub>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>w</mi>
<mi>I</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中:μ为步长,为求偏导数运算符;下标R和I分别表示实部和虚部,ρR表示高斯聚类噪声方差ρ的实部;ρI表示高斯聚类噪声方差ρ的虚部; |.|表示绝对值运算;j为虚数单位;
参数ρ和μ的更新方程如下:
<mrow>
<mfenced open='' close=''>
<mtable>
<mtr>
<mtd>
<mi>&mu;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&mu;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&gamma;&mu;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>|</mo>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>&rho;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<mi>&gamma;&mu;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>|</mo>
<mi>e</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>n</mi>
<mo>*</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中判决误差e(n*)为:
e(n*)=d(n*)-z(n*)
其中:d(n*)为判决值;sgn(·)为符号函数;和分别代表实部和虚部;利用泰勒公式将ρ和μ的更新方程改写为:
μ(n*+1)=μ(n*)(1+β+β2+…)≈μ(n*)(1+β)
ρ(n*+1)=ρ(n*)(1+β+β2+…)≈ρ(n*)(1+β)
这里β=-γμ(n*)|e(n*)|2,n*表示第n*次迭代;
7)将经步骤6)更新后的权值代入步骤5))中的代价函数Jcluster(w,z(n*)),直至该代价函数值不再变化,退出。
2.根据权利要求1所述的基于聚类算法的相干光通信盲均衡方法,其特征在于,M为4,8,16,32,64…。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710358012.5A CN107302507B (zh) | 2017-05-19 | 2017-05-19 | 基于聚类算法的相干光通信盲均衡方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710358012.5A CN107302507B (zh) | 2017-05-19 | 2017-05-19 | 基于聚类算法的相干光通信盲均衡方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107302507A true CN107302507A (zh) | 2017-10-27 |
CN107302507B CN107302507B (zh) | 2020-03-20 |
Family
ID=60137352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710358012.5A Expired - Fee Related CN107302507B (zh) | 2017-05-19 | 2017-05-19 | 基于聚类算法的相干光通信盲均衡方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107302507B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108199777A (zh) * | 2018-01-12 | 2018-06-22 | 温州大学 | 基于概率密度函数拟合和模糊逻辑的相干光通信系统盲均衡方法 |
CN110933005A (zh) * | 2019-12-09 | 2020-03-27 | 北京理工大学 | 一种密度聚类的调制格式识别与osnr估计的联合方法 |
CN114337843A (zh) * | 2021-12-23 | 2022-04-12 | 武汉邮电科学研究院有限公司 | 一种用于相干光通信的自适应信道均衡方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8155244B2 (en) * | 2008-07-08 | 2012-04-10 | The Boeing Company | Demodulation using blind constellation identification for communication signals |
CN103023840A (zh) * | 2012-12-04 | 2013-04-03 | 温州大学 | 多输入多输出泛函网络实现无线激光通信电域信号盲均衡方法 |
-
2017
- 2017-05-19 CN CN201710358012.5A patent/CN107302507B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8155244B2 (en) * | 2008-07-08 | 2012-04-10 | The Boeing Company | Demodulation using blind constellation identification for communication signals |
CN103023840A (zh) * | 2012-12-04 | 2013-04-03 | 温州大学 | 多输入多输出泛函网络实现无线激光通信电域信号盲均衡方法 |
Non-Patent Citations (3)
Title |
---|
MOHAMED L. AMMARI等: "A clustering technique for blind equalization of MIMO frequency selective channels", 《EUROPEAN WIRELESS 2011》 * |
RANGEET MITRA等: "Improved Multi-Stage Clustering Based Blind Equalisation in Distributed Environments", 《IEEE COMPUTER SOCIETY》 * |
季怿: "基于聚类技术的盲均衡算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108199777A (zh) * | 2018-01-12 | 2018-06-22 | 温州大学 | 基于概率密度函数拟合和模糊逻辑的相干光通信系统盲均衡方法 |
CN110933005A (zh) * | 2019-12-09 | 2020-03-27 | 北京理工大学 | 一种密度聚类的调制格式识别与osnr估计的联合方法 |
CN114337843A (zh) * | 2021-12-23 | 2022-04-12 | 武汉邮电科学研究院有限公司 | 一种用于相干光通信的自适应信道均衡方法及系统 |
CN114337843B (zh) * | 2021-12-23 | 2023-11-28 | 武汉邮电科学研究院有限公司 | 一种用于相干光通信的自适应信道均衡方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107302507B (zh) | 2020-03-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108880692B (zh) | 一种相干光通信系统调制格式识别和光信噪比监控方法 | |
US9515742B2 (en) | System and methods for adaptive equalization for optical modulation formats | |
Li et al. | Convergence analysis of finite length blind adaptive equalizers | |
CN107302507B (zh) | 基于聚类算法的相干光通信盲均衡方法 | |
KR100289169B1 (ko) | 전환알고리즘을이용한적응등화기의블라인드수렴을향상시키기위한방법및장치 | |
JP5898374B2 (ja) | 非線形交差偏波軽減アルゴリズム | |
US8995834B2 (en) | Blind equalization for polarization-switched QPSK optical communications | |
CN101599801A (zh) | 滤波器系数调整装置和方法 | |
CN102461035A (zh) | 用于对偏振分集复用信号进行盲解复用的方法和装置 | |
CN108199777A (zh) | 基于概率密度函数拟合和模糊逻辑的相干光通信系统盲均衡方法 | |
KR20190121022A (ko) | 광전송 분산 보상 등화기 장치 및 방법 | |
CN1301440A (zh) | 判决指导的相位探测器 | |
Freire et al. | Domain adaptation: The key enabler of neural network equalizers in coherent optical systems | |
CN113938198B (zh) | 光纤传输系统、基于lda简化非线性均衡器的方法及模块 | |
US6697424B1 (en) | Fast convergent pipelined adaptive decision feedback equalizer using post-cursor processing filter | |
CN112887237B (zh) | 光纤通信系统的复信道均衡器设计方法 | |
CA2310999A1 (en) | The blind equalization algorithm with joint use of the constant r and the sliced symbols | |
Bosco et al. | A novel update algorithm in stokes space for adaptive equalization in coherent receivers | |
CN1312867C (zh) | 偏振模式色散补偿器及其自适应调节方法 | |
Luo et al. | Training-based symbol detection with temporal convolutional neural network in single-polarized optical communication system | |
Wu et al. | Research of Blind Equalization Technology of Coherent Optical Communication Based on Decidion Optimization | |
Yan et al. | Demultiplexing in mode-division multiplexing system using multichannel blind deconvolution | |
US9225433B1 (en) | Timing recovery for optical coherent receivers in the presence of polarization mode dispersion | |
Zhang et al. | Machine learning for long-haul optical systems | |
Jiang et al. | Blind Equalization Algorithms Applied to Dispersion Compensation in POLMUX Coherent Optical Communication System |
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 | ||
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: 20200320 |