江西旅游网站建设方案_惠州市seo网站设计_长春网站优化平台_手机优化大师为什么扣钱

时间:2025/7/10 8:21:17来源：https://blog.csdn.net/sonnet00/article/details/139841265 浏览次数:0次

斯托克斯向量（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 npdef WavePlate(delta):Mwp = np.zeros([4, 4])Mwp[0, 0] = 1Mwp[1, 1] = 1Mwp[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] = 1R[3, 3] = 1R[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 Rdef rotate(M,theta):R1 = RotateMatrix(theta)R2 = RotateMatrix(-theta)return R2@M@R1Mp = 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.]]

关键字：江西旅游网站建设方案_惠州市seo网站设计_长春网站优化平台_手机优化大师为什么扣钱

本网仅为发布的内容提供存储空间，不对发表、转载的内容提供任何形式的保证。凡本网注明“来源：XXX网络”的作品，均转载自其它媒体，著作权归作者所有，商业转载请联系作者获得授权，非商业转载请注明出处。

我们尊重并感谢每一位作者，均已注明文章来源和作者。如因作品内容、版权或其它问题，请及时与我们联系，联系邮箱：809451989@qq.com，投稿邮箱：809451989@qq.com

责任编辑：