激光与光电子学进展, 2018, 55 (12): 122603, 网络出版: 2019-08-01   

基于邦加球轨迹的穆勒矩阵测量法 下载: 1990次

Measurement Method of Müller Matrix Based on Trajectories on Poincare Sphere
作者单位
上海大学特种光纤与光接入网省部共建重点实验室, 上海 200072
摘要
提出了基于邦加球轨迹的穆勒矩阵测量法,将极分解得到的矩阵与其在邦加球上的轨迹联系起来,通过输入光和输出光的偏振态在邦加球上的位置关系,直接获取非消偏元件的双折射矩阵和二向色性矩阵。结果表明:基于所提方法,偏振控制系统只需要一个可旋转起偏器,在测试前后提供两个确定的偏振态,即可直接获得光纤器件的双折射矩阵和二向色性矩阵,且在测量过程中无需改变待测光纤的状态。
Abstract
A measurement method of Müller matrix is proposed based on the trajectories on the Poincare sphere. A matrix obtained by the polar decomposition is related to its trajectories on the Poincare sphere. Therefore, the birefringence and dichroism matrices of a nondepolarizing unit are directly obtained via the position relationship of the input and output light polarization states on the Poincare sphere. The results show that, based on the proposed method, only one rotatable polarizer is required in the polarization control system. With the measurement of two polarization states determined before and after tests, not only the birefringence and dichroism matrices can be directly obtained, but also the state of the fiber to detect does not need change during the measuring process.

1 引言

随着光纤技术的迅猛发展,偏振效应在光纤领域内的重要性日益凸显[1-4]。由于光纤体积小、质量轻,因此光纤系统比传统光学系统的集成度高,但光纤中的偏振态很容易受到外界环境的影响。在高精度的光纤传感系统和高速光纤通信系统中,偏振问题更不可忽视[5-8]。通过琼斯矩阵法和穆勒矩阵法可以获得光纤器件的偏振模色散和偏振相关损耗等偏振信息[9]。琼斯矩阵法适用于处理偏振光的相干问题,但由于琼斯矢量为复数形式,不方便测量[10]。穆勒矩阵法可直接测量光强,是实验中经常采用的测量方法[11-12]。目前穆勒矩阵的测量方法主要包括直接测量法、半自动测量法和全自动测量法[13],这三种测量方法的测量原理基本相同,都是通过调整两个1/4波片的方位角获得以样品穆勒矩阵的各元素为未知数的线性无关方程组,求解方程组可获得穆勒矩阵。

穆勒矩阵并没有明确的物理意义,很难从中直接获得光纤器件的偏振模色散和偏振相关损耗等,因此很多学者研究了穆勒矩阵的分解方法。Cloude[14]提出所有的穆勒矩阵都是4个非消偏穆勒矩阵之和;Gil和Bernabeu[15]将极分解应用到穆勒矩阵中;Lu和Chipman[16]则较为系统地阐述了穆勒矩阵的极分解法,将穆勒矩阵分解为3个矩阵:二向色性矩阵(对应偏振相关损耗)、双折射矩阵和退偏矩阵。对于非消偏元件,其穆勒矩阵由双折射矩阵和二向色性矩阵组成。Ossikovski[17]利用奇异值分解,将非消偏穆勒矩阵分解为两个线性双折射元件和一个线性双折射衰减器。极分解法是目前比较主流的穆勒矩阵分解法。

光纤器件中的偏振态容易受到影响,且很难用肉眼判断光纤的内在偏振主轴,因此在测量时应尽量避免改变待测光纤器件的状态。全自动的偏振控制精度高、速度快,但系统复杂性高、成本昂贵。本文将极分解得到的矩阵与邦加球轨迹联系起来,通过输入光和输出光偏振态在邦加球上位置的关系,可直接获取双折射矩阵和二向色性矩阵。只需在偏振控制系统中安装一个可旋转起偏器,并在测试前和测试后提供两个不同方向的线偏振光,即可获得光纤器件的偏振模色散和偏振相关损耗等信息。

2 穆勒矩阵的极分解

斯托克斯空间法是描述偏振光的重要方法,该方法可将所有偏振光与空间内半径为1的一个球(邦加球)内的坐标一一对应(图1)。设任意偏振光的表达式为

图 1. 偏振光在斯托克斯空间内的表示

Fig. 1. Polarized light in Stokes space

下载图片 查看所有图片

Ex,y(t)=Pcosαexp()sinα+2-2P2expjδ~xtexpjδ~yt,(1)

式中P为完全偏振光能量与总能量之比,即光的偏振度;δx方向与y方向偏振光的相位差;αx方向与y方向偏振光的振幅比; δ~x(t)为x方向偏振光的随机相位; δ~y(t)为y方向偏振光的随机相位。根据斯托克斯矢量的定义,可以推导出(1)式所表示的偏振光在斯托克斯空间内坐标 S1,S2,S3下各个分量的表达式为

S1=Pcos2αS2=Psin(2α)cosδS3=Psin(2α)sinδ(2)

从(2)式可以看出,偏振光在斯托克斯空间内的矢径由其偏振度P决定,在半径为P的球面上的位置由偏振光中完全偏振光的参数αδ决定。设线段OS与平面S2OS3的夹角为β,则α=π/4-β/2。图1中,S1轴的两端点分别为水平偏振光(H)和垂直偏振光(V),S2轴的两端点分别为45°偏振光和-45°偏振光,S3轴的两端点分别为左旋圆偏振光(LC)和右旋圆偏振光(RC)。

对于非消偏元件,入射光和透射光的斯托克斯矢量SinSout都在邦加球的球面上,它们之间的关系可由穆勒矩阵表示为

1Sout=M1Sin,(3)

式中M为穆勒矩阵,可以分解为一个双折射矩阵MR和一个二向色性矩阵MD,即

M=MRMD,(4)

其中

MR=100mR,(5)MD=K1DTDmD,(6)

式中mRmD均为三阶子矩阵;D为衰减矢量;K为归一化系数,具体参见文献[ 16]。

3 单一功能器件的穆勒矩阵和邦加球轨迹

3.1 双折射矩阵

假设双折射矩阵MR的本征偏振态V1V2分别为

V1=cosξexp(),sinξT,(7)V2=-sinξexp(),cosξT,(8)

式中ξ为偏振态中x分量与y分量的振幅比,ϑ为两偏振分量的相位差。本征偏振态在邦加球上的坐标为

M1=cos(2ξ),sin(2ξ)cosϑ,sin(2ξ)sinϑ,(9)M2=-cos(2ξ),-sin(2ξ)cosϑ,-sin(2ξ)sinϑ(10)

如果将邦加球的OS1旋转到M1M2上(先将OS2OS3OS1逆时针旋转ϑ,得到OS1tOS2tOS3t;再将OS1tOS2tOS3t顺时针旋转2ξ),得到新坐标系S1nS2nS3n。偏振光在以V1V2为基矢时,相干矩阵C1,2由其在坐标系S1nS2nS3n中的坐标决定[18]:

C1,2=121+S1nS2n+jS3nS2n-jS3n1-S1n(11)

对于本征偏振态为V1V2的双折射元件,其在V1V2基矢下的琼斯矩阵为

JR=exp/200exp-/2,(12)

因此输出光的相干矩阵为

C1,2(out)=12exp(/2)00exp(-/2)1+S1nS2n+jS3nS2n-jS3n1-S1nexp(-/2)00exp(/2)=121+S1nS2n+jS3nexp()S2n-jS3nexp(-)1-S1n(13)

输出光在坐标系S1nS2nS3n中坐标(r1n,r2n,r3n)下的各个分量可表示为

r1n=S1n,(14)r2n=S2ncosδ-S3nsinδ,(15)r3n=S2nsinδ+S3ncosδ(16)

因此,对于本征偏振态为V1V2的双折射元件,双折射相位差δ 的引入将使入射光绕着OS1n旋转δ,如图2(a)所示。S1S2S3坐标系下的mR可通过坐标转换获得:

mR=m-1(ξ,ϑ)·1000cosδ-sinδ0sinδcosδm(ξ,ϑ),(17)

式中m(ξ,ϑ)为坐标系S1S2S3与坐标系S1nS2nS3n的转换矩阵。

图 2. 邦加球轨迹。(a)双折射元件;(b)二色性元件

Fig. 2. Trajectories on Poincare sphere. (a) Birefringent unit; (b) dichromatic unit

下载图片 查看所有图片

3.2 二向色性矩阵

对于本征偏振态为V1V2的二色性元件,衰减系数为D,其琼斯矩阵为

JD=k1+D001-D,(18)

式中k为能量系数。输出光的相干矩阵为

C1,2(out)=k21+D001-D1+S1nS2n+jS3nS2n-jS3n1-S1n1+D001-D=k2(1+D)1+S1n1-D2S2n+jS3n1-D2S2n-jS3n(1-D)1-Sin(19)

将输出光的相干矩阵归一化,可得到输出光偏振态在坐标系S1nS2nS3n中坐标(r1n,r2n,r3n)下的各个分量为

r1n=D+S1n1+DS1n,(20)r2n=1-D2S2n1+DS1n,(21)r3n=1-D2S3n1+DS1n(22)

从(20)~(22)式中可以看出,对于本征偏振态为V1V2的二色性元件,D的引入将使入射光沿着圆S1nOSin运动,运动弧度由DS1n决定,如图2(b)所示。S1S2S3坐标系下的穆勒矩阵可通过坐标转换获得:

MD=K100m-1(ξ,ϑ)1D00D100001-D200001-D2100mξ,ϑ=K1DTDmD,(23)

图 3. 偏振态测量装置

Fig. 3. Measurement setup for polarization states

下载图片 查看所有图片

其中K可保证入射光和出射光的坐标都在邦加球的球面上,其大小与入射光的坐标值S1n有关。

4 基于极分解的邦加球轨迹测量法

若入射光的斯托克斯参量为 1,SinT,则出射光的斯托克斯参量为

1Sout=K1DTmRDmRmD1Sin=1mRD+mRmDSin1+DTSin(24)

当入射光偏振态相同,入射波长取有微小差别的两个波长λ1λ2时,有

Soutλ1=11+DTSinmRλ1D+mRλ1mDSin=mRλ11+DTSinD+mDSin,(25)Soutλ2=11+DTSinmRλ2D+mRλ2mDSin=mRλ21+DTSinD+mDSin(26)

可以看出,Soutλ1Soutλ2在邦加球上的位置是由同一个矢量D+mDSin分别与mRλ1/ 1+DTSinmRλ2/ 1+DTSin相乘所得,其中1+DTSin为一个数,mRλ1mRλ2为双折射旋转矩阵,二者分别使D+mDSin绕着双折射元件本征态所确定的轴旋转δ1δ2。根据双折射的定义,可得

δ=-λΔλδ1-δ2,(27)

式中Δλ为波长变化量。因此,通过旋转的轨迹可以获得mR所对应的本征态及双折射引起的角度旋转δ,也就得到了MRMD由衰减矢量D0决定,根据前文的分析,MD使入射光沿着圆S1nOSin运动,其中S1n即为归一化衰减矢量D0。为了获得D0的坐标,需要提供两组输入/输出偏振态 Sin1,Sout1Sin2,Sout2,另一组输入/输出偏振态可通过旋转偏振片获取。中间态S'out1S'out2可由极分解获得:

1S'out1S'out2=MR-11Sout1Sout2=MD1Sin1Sin2(28)

D0即为圆S'out1OSin1和圆S'out2OSin2的交线,D值由Sin1S'out1Sin2S'out2的坐标获得。

5 实验和讨论

实验装置如图3所示,可调激光器通过一个偏振方向可旋转的线偏振光起偏器后,其出射的激光进入待测光纤,输出光偏振态由偏振态分析仪测量。实验中待测光纤器件由一段圆二向色性光纤和一段线性双折射光纤组成。调整起偏器的角度,使输入待测光纤器件的光为水平偏振光,改变可调激光器的波长,分别设定λ1=1.55 μm,λ2=1.5504 μm,λ3=1.5508 μm。测得3个输出光在邦加球上的坐标分别为P1=(0.0045,-0.6214,-0.7834),P2=(0.0112,-0.3827,-0.9237)和P3=(-0.0009,-0.8124,-0.5832)。根据空间几何关系,可以得到本征态对应的邦加球上的坐标为M1=(-0.9996,0.0283,0)和M2=(0.9996,-0.0283,0)。双折射轨迹的圆心坐标O'=(0.0221,-0.0006,0),半径为0.9997,O'P1O'P2间的弧度角为

θ=acos2×0.99972-0.276922×0.9997×0.9997=0.2779rad,(29)

则双折射旋转角δ=1077 rad,因此有

m'R=1000cosδ-sinδ0sinδcosδ=1000-0.97010.24260-0.2426-0.9701(30)

双折射矩阵MR

MR=10000 0.9985-0.0516-0.01520-0.0516-0.8424-0.53630 0.0152 0.5363-0.8439(31)

双折射矩阵所表示的本征态和邦加球轨迹如图4所示,图中蓝色的点表示两个本征态的位置,红色的点为测量值,黑色的点为双折射矩阵的轨迹,显然此时的双折射为线性双折射。

图 4. 测得双折射矩阵的邦加球轨迹

Fig. 4. Trajectories on Poincare sphere of measured birefringence matrix

下载图片 查看所有图片

一组输入/输出偏振态即为Sin1=[1,0,0]T,Sout1=[0.0045,-0.6214,-0.7834]T。旋转起偏器,测得Sout2=[0.0067,-0.6255,-0.7802]T,此时取下被测光纤器件,测得Sin2=[0.9511,0.3090,0]T。根据双折射矩阵,可以得到S'out1=[0.0250,0.1031,0.9943],S'out2=[0.0274,0.1082,0.9938]。根据空间几何可得平面Sin1OS'out1Sin2OS'out2的交线与邦加球的交点为D1=(0.0116,0.1031,0.9946)和D2=(-0.0116,-0.1031,-0.9946)。分析可得衰减矢量方向为D0=[0.0116,0.1031,0.9946],衰减值D=0.99990,由于衰减矢量基本接近圆偏振光,因此二向色性光纤为圆二色性光纤,其二向色性矩阵MD

MD=K1.00000.01160.10310.99450.01160.01450.00120.01140.10310.00120.02480.10110.99450.01140.10110.9894(32)

因此,根据(4)式可得待测光纤器件的穆勒矩阵为

M=K 1.0000 0.0116 0.1031 0.9945-0.0089 0.0142-0.0016-0.0089-0.6208-0.0078-0.0752-0.6163-0.7838-0.0087-0.0720-0.7806(33)

为了验证实验的正确性,将待测器件中的圆二色性光纤与线性双折射光纤分离。首先测量线性双折射光纤的拍长,采用宽谱光源干涉法[19],测得双折射光纤长度LB=575 mm,拍长L'B=3.4 mm。采用轨迹法计算得到光纤拍长L'B=3.35 mm。可以看出,两种方法得到的光纤拍长基本吻合。接着测量了圆二色性光纤的偏振相关损耗。输入水平方向线偏振光,测得圆二色性光纤输出光在邦加球上的坐标为(0.0192,0.1165,0.9930),两本征态的消光比大于40 dB,与用轨迹法得到的S'out1基本一致。实验表明:利用所提出的邦加球轨迹法测得的光纤器件的特性是可信的,但由于线性双折射光纤的相位差较易受外界环境影响,因此在测量过程中会引入一定的误差。另外,轨迹法中空间曲线和曲面方程以及矩阵的计算也会引入一定的计算误差。

6 结论

将极分解得到的矩阵与其邦加球轨迹联系起来,通过输入光和输出光偏振态在邦加球上位置的关系,直接获取了非消偏元件的双折射矩阵和二向色性矩阵。研究表明:该实验方法中采用的偏振控制系统只需要一个可旋转的起偏器,在测试前和测试后提供两个确定的偏振态,即可直接获得非消偏元件的双折射矩阵和二向色性矩阵,从而获得光纤器件的偏振模色散和偏振相关损耗等参数。该测量过程操作简单,无需改变待测光纤的状态,为光纤器件的偏振特性分析提供了一种便捷直观的方法。

参考文献

[1] 董自凯, 徐润亲, 田金荣, 等. 全光纤偏振输出锁模光纤激光器[J]. 激光与光电子学进展, 2017, 54(8): 081404.

    Dong Z K, Xu R Q, Tian J R, et al. All-fiber based mode-locked fiber laser with polarization output[J]. Laser & Optoelectronics Progress, 2017, 54(8): 081404.

[2] 杨军, 苑勇贵, 喻张俊, 等. 光学相干域偏振测量技术及其在高精度光纤陀螺器件测量中的应用[J]. 光学学报, 2018, 38(3): 0328007.

    Yang J, Yuan Y G, Yu Z J, et al. Optical coherence domain polarimetry technology and its application in measurement for evaluating components of high precision fiber-optic gyroscopes[J]. Acta Optica Sinica, 2018, 38(3): 0328007.

[3] Goossens J W, Yousefi M I, Jaouën Y, et al. Polarization-division multiplexing based on the nonlinear Fourier transform[J]. Optics Express, 2017, 25(22): 26437-26452.

[4] 李珊珊, 郝霞, 白晋军, 等. 偏振可调的太赫兹单偏振单模光子晶体光纤[J]. 中国激光, 2016, 43(9): 0901005.

    Li S S, Hao X, Bai J J, et al. Terahertz single-polarization single-mode photonic crystal fiber with tunable polarization[J]. Chinese Journal of Lasers, 2016, 43(9): 0901005.

[5] 李鹏程, 刘琨, 江俊峰, 等. 基于FPGA的分布式光纤传感系统偏振控制研究[J]. 中国激光, 2017, 45(5): 051002.

    Li P C, Liu K, Jiang J F, et al. Research on polarization control of distributed optical fiber sensing system based on FPGA[J]. Chinese Journal of Lasers, 2017, 45(5): 051002.

[6] 谢泽, 杨远洪, 闫晗, 等. 光纤陀螺光纤环分布偏振耦合误差预估[J]. 中国激光, 2017, 44(10): 1006004.

    Xie Z, Yang Y H, Yan H, et al. Estimation of distributed polarization coupling error in optic fiber coil of fiber optic gyroscope[J]. Chinese Journal of Lasers, 2017, 44(10): 1006004.

[7] 李传生, 张朝阳, 孙海江, 等. 保偏延迟光纤环偏振串音对光纤电流互感器的影响[J]. 中国激光, 2014, 41(4): 1105008.

    Li C S, Zhang C Y, Sun H J, et al. Effect of polarization crosstalk of polarization-maintaining delay optical fiber coil on the fiber-optic current sensor[J]. Chinese Journal of Lasers, 2014, 41(11): 1105008.

[8] 高洋洋, 周卫宁, 雷莉莉, 等. 光纤陀螺用超辐射发光二极管启动偏振特性及其影响研究[J]. 激光与光电子学进展, 2015, 52(11): 112302.

    Gao Y Y, Zhou W N, Lei L L, et al. Research on polarization characteristics of SLD start-up used in fiber optic gyroscope and its effect[J]. Laser & Optoelectronics Progress, 2015, 52(11): 112302.

[9] Yao X S, Chen X J, Liu T G. High accuracy polarization measurements using binary polarization rotators[J]. Optics Express, 2010, 18(7): 6667-6685.

[10] 黄山, 赵华凤, 俞涛, 等. 对琼斯矩阵测量PMD的改进[J]. 光电子·激光, 2004, 15(1): 120-122.

    Huang S, Zhao H F, Yu T, et al. Improved Jones matrix eigenanalysis method in measuring PMD[J]. Journal of Optoelectronics·Laser, 2004, 15(1): 120-122.

[11] Xu Z X, Yao X S, Ding Z Y, et al. Accurate measurements of circular and residual linear birefringences of spun fibers using binary polarization rotators[J]. Optics Express, 2017, 25(24): 30780-30792.

[12] 黄爱显, 张昊, 邱怡申, 等. 斯托克斯矢量法在全光纤电流传感器中的应用[J]. 光子学报, 2013, 42(6): 699-704.

    Huang A X, Zhang H, Qiu Y S, et al. Application of Stokes vector method on the full-fiber current sensor[J]. Acta Photonica Sinica, 2013, 42(6): 699-704.

[13] 陈立刚. 偏光器件Muller矩阵的智能化测量[D]. 曲阜: 曲阜师范大学, 2004.

    Chen LG. Smart measurement of the Muller matrix of polarization devices[D]. Qufu: Qufu Normal University, 2004.

[14] Cloude S R. Group theory and polarization algebra[J]. Optik, 1986, 31(11): 26-36.

[15] Gil J J, Bernabeu E. Obtainment of the polarizing and retardation parameters of a non-depolarizing optical system from the polar decomposition of its Mueller matrix[J]. Optik, 1987, 76(2): 67-71.

[16] Lu S Y, Chipman R A. Interpretation of Mueller matrices based on polar decomposition[J]. Journal of the Optical Society of America A, 1996, 13(5): 1106-1113.

[17] Ossikovski R. Interpretation of nondepolarizing Mueller matrices based on singular-value decomposition[J]. Journal of the Optical Society of America A, 2008, 25(2): 473-482.

[18] 肖悦娱, 杨辉祥, 徐怀宝, 等. 基于斯托克斯空间的任意正交基矢下琼斯矩阵的测量[J]. 激光与光电子学进展, 2016, 53(9): 092604.

    Xiao Y Y, Yang H X, Xu H B, et al. Measuring Jones matrix of arbitrary orthogonal basic vectors based on Stokes space[J]. Laser & Optoelectronics Progress, 2016, 53(9): 092604.

[19] 米剑, 张春熹, 李铮, 等. 利用宽谱光源及偏光干涉测量保偏光纤拍长[J]. 光电子·激光, 2006, 17(9): 1074-1077.

    Mi J, Zhang C X, Li Z, et al. Measuring the beatlength of polarization maintaining fiber by broadband light source and conoscopic interference[J]. Journal of Optoelectronics·Laser, 2006, 17(9): 1074-1077.

肖悦娱, 蒋晓勇, 陈华. 基于邦加球轨迹的穆勒矩阵测量法[J]. 激光与光电子学进展, 2018, 55(12): 122603. Yueyu Xiao, Xiaoyong Jiang, Hua Chen. Measurement Method of Müller Matrix Based on Trajectories on Poincare Sphere[J]. Laser & Optoelectronics Progress, 2018, 55(12): 122603.

本文已被 4 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!