偏振测量——典型光学元件的穆勒矩阵

sonnet00

于 2025-01-04 18:31:12 发布

阅读量1.5k

点赞数 4

分类专栏：光学量测文章标签： python

本文链接：https://blog.csdn.net/sonnet00/article/details/139841265

版权

光学量测专栏收录该内容

16 篇文章

订阅专栏

斯托克斯向量（Stokes Vector）为1*4向量，可用于描述偏振光的状态，其四个元素分别代表总光强S0，垂直或水平方向线偏振分量S1，正负45度线偏振分量S2，圆偏振分量S3。

当一束光经过某个器件时，其偏振态发生变化，对应到斯托克斯向量的变化。在数学上即可考虑到用一个4*4的矩阵描述该器件对光的偏振态的影响，这个矩阵就是穆勒矩阵。这里记录一些典型光学元件的穆勒矩阵以供查询。

线性偏振片（偏振方向为0度）

$M_P=\begin{bmatrix} 1 & 1 & 0& 0 \\ 1 & 1 & 0& 0 \\ 0 & 0 & 0& 0 \\ 0 & 0 & 0& 0 \end{bmatrix}$

玻片（waveplate 或补偿器compensator 或延迟器retarder），其中phi为玻片的延迟量，快轴角度为0

$M_{WP}=\begin{bmatrix} 1 & 0 & 0& 0 \\ 0 & 1 & 0& 0 \\ 0 & 0 & cos\delta& -sin\delta\\ 0 & 0 & sin\delta& cos\delta\end{bmatrix}$

另外，穆勒矩阵的旋转矩阵如下，其中theta为旋转角度

$R(\theta)=\begin{bmatrix} 1 & 0 & 0& 0 \\ 0 & cos2\theta & sin2\theta& 0 \\ 0 & -sin2\theta & cos2\theta& 0\\ 0 & 0 & 0& 1\end{bmatrix}$

一个光学元件旋转一定角度后，其穆勒矩阵可以如下计算：

$M(\theta) = R(-\theta)M_0R(\theta)$

进而可以求得任意玻片和偏振片旋转一定角度后的穆勒矩阵。

以下python代码

import numpy as np

def WavePlate(delta):
    Mwp = np.zeros([4, 4])
    Mwp[0, 0] = 1
    Mwp[1, 1] = 1
    Mwp[2, 2] = np.cos(delta)
    Mwp[3, 3] = np.cos(delta)
    Mwp[2, 3] = -np.sin(delta)
    Mwp[3, 2] = np.sin(delta)
    return Mwp
def RotateMatrix(theta):
    R = np.zeros([4,4])
    R[0, 0] = 1
    R[3, 3] = 1
    R[1, 1] = np.cos(2 * theta)
    R[2, 2] = np.cos(2 * theta)
    R[1, 2] = np.sin(2 * theta)
    R[2, 1] = -np.sin(2 * theta)
    return R

def rotate(M,theta):
    R1 = RotateMatrix(theta)
    R2 = RotateMatrix(-theta)
    return R2@M@R1

Mp = np.zeros([4,4])
Mp[0:2,0:2]=1
print('角度为0度的线性偏振片穆勒矩阵： \n',Mp)
print('角度为45度的线性偏振片穆勒矩阵： \n',np.round(rotate(Mp,np.pi/4)))
print('角度为-45度的线性偏振片穆勒矩阵： \n',np.round(rotate(Mp,-np.pi/4)))
print('角度为90度的线性偏振片穆勒矩阵： \n',np.round(rotate(Mp,np.pi/2)))

delta1 = np.pi/4 #四分之一玻片
delta2 = np.pi/2 #半玻片

print('快轴角度为0度的四分之一玻片穆勒矩阵： \n',np.round(WavePlate(delta1),2))
print('快轴角度为90度的四分之一玻片穆勒矩阵： \n',np.round(rotate(WavePlate(delta1),np.pi/2),2))
print('快轴角度为0度的半玻片穆勒矩阵： \n',np.round(WavePlate(delta2),2))
print('快轴角度为90度的半玻片穆勒矩阵： \n',np.round(rotate(WavePlate(delta2),np.pi/2),2))

运行结果：

角度为0度的线性偏振片穆勒矩阵： 
 [[1. 1. 0. 0.]
 [1. 1. 0. 0.]
 [0. 0. 0. 0.]
 [0. 0. 0. 0.]]
角度为45度的线性偏振片穆勒矩阵： 
 [[1. 0. 1. 0.]
 [0. 0. 0. 0.]
 [1. 0. 1. 0.]
 [0. 0. 0. 0.]]
角度为-45度的线性偏振片穆勒矩阵： 
 [[ 1.  0. -1.  0.]
 [ 0.  0. -0.  0.]
 [-1. -0.  1.  0.]
 [ 0.  0.  0.  0.]]
角度为90度的线性偏振片穆勒矩阵： 
 [[ 1. -1.  0.  0.]
 [-1.  1. -0.  0.]
 [ 0. -0.  0.  0.]
 [ 0.  0.  0.  0.]]
快轴角度为0度的四分之一玻片穆勒矩阵： 
 [[ 1.    0.    0.    0.  ]
 [ 0.    1.    0.    0.  ]
 [ 0.    0.    0.71 -0.71]
 [ 0.    0.    0.71  0.71]]
快轴角度为90度的四分之一玻片穆勒矩阵： 
 [[ 1.    0.    0.    0.  ]
 [ 0.    1.   -0.    0.  ]
 [ 0.   -0.    0.71  0.71]
 [ 0.   -0.   -0.71  0.71]]
快轴角度为0度的半玻片穆勒矩阵： 
 [[ 1.  0.  0.  0.]
 [ 0.  1.  0.  0.]
 [ 0.  0.  0. -1.]
 [ 0.  0.  1.  0.]]
快轴角度为90度的半玻片穆勒矩阵： 
 [[ 1.  0.  0.  0.]
 [ 0.  1. -0.  0.]
 [ 0. -0.  0.  1.]
 [ 0. -0. -1.  0.]]