摘要:分振幅型全Stokes同时偏振成像仪具有实时性好、空间分辨率高、精度高等优点, 有很高的应用价值. 分振幅型全Stokes同时偏振成像系统利用偏振分束器、1/2波片和1/4波片将入射光Stokes矢量调制在4幅图像中, 可解析入射光Stokes矢量. 1/2波片和1/4波片的相位延迟误差对Stokes矢量测量精度有着不可忽略的影响. 建立了包含上述两种误差的Stokes矢量测量误差方程, 分析了1/2波片和1/4波片相位延迟耦合误差对自然光、0°/45°线偏光、左旋圆偏光等典型基态入射光的Stokes矢量测量误差的影响, 推导了任意偏振态的Stokes矢量测量误差的表征方法. 在邦加球球面和球内选取不同偏振度的Stokes矢量作为入射光进行仿真. 结果表明, Stokes矢量测量误差和偏振度测量误差均随着入射光偏振度的增大而增大. 选取入射光偏振度为1时的偏振测量精度评估系统. 为满足2%的偏振测量精度, 1/2波片相位延迟误差应在±1.6°内, 1/4波片相位延迟误差应在±0.5°内. 这对提高系统的偏振测量精度具有重要意义, 为系统设计和研制提供了重要的理论指导.
关键词: 偏振成像/
波片/
相位延迟误差/
Stokes矢量

English Abstract

全文HTML

--> --> -->

偏振是光的一种本质属性, 偏振信息不仅包含强度信息还包含偏振度、偏振角和椭圆度角等信息. 光与大气或目标相互作用时会改变光的偏振特性, 光的偏振特性的变化可以表征大气或目标的物理属性, 因此偏振探测广泛用于仿生导航、水下目标探测、天文探测、大气遥感、空间目标3D重建和生物医学等领域^[1?6]. 为了实现目标偏振信息的探测, 需要研制高精度的偏振成像系统, 偏振成像系统可分为两类: 分时偏振成像系统和同时偏振成像系统. 分时偏振成像系统中普遍存在转动部件或分时调制器, 导致分时偏振成像系统无法偏振探测快速变化的目标. 为实现动态目标的偏振探测, 同时偏振成像系统的研制已成为研究热点. 1982年, Azzam^[7]最先提出一种基于振幅分割的同时偏振成像系统, 该系统利用镀膜分光器和Wollaston棱镜将入射光分为4束光, 由4个探测器实现了对入射光Stokes参数线偏振分量的同时测量. 2008年, Pezzaniti等^[8]成功研制了分振幅全Stokes同时偏振成像系统, 利用偏振分束器(polarization beam splitter, PBS)、1/2波片(half wave plate, HWP)和1/4波片(quarter wave plate, QWP)将入射光分为4束光, 由4个探测器实现对入射光全Stokes参数的同时测量. 2003年, Oka和Kaneko^[9]首次提出一种基于双折射楔形棱镜的通道调制型偏振成像仪, 随后国内外****对基于双折射楔形棱镜的通道调制型偏振成像仪进行了改进^[10?12], 提高了基于双折射楔形棱镜的通道调制型偏振成像仪的偏振探测性能. 2016年, 权乃承等^[13]提出一种基于孔径分割与视场分割的通道型成像光谱偏振技术, 并进行仿真分析验证了该方案的可行性. 2018年, 冯斌等^[14]研制了分焦平面同时偏振成像仪, 并给出了性能评估模型. 为了提高偏振测量精度, 国内外研究人员针对不同调制原理的同时偏振成像系统进行了系统参数误差分析等方面的工作^[15?19].
分振幅型全Stokes同时偏振成像系统没有时间和空间的失调, 具有实时性好、空间分辨率高、精度高等优点, 同时还能实现全Stokes参数的测量, 具有很高的应用价值. 波片是分振幅型全Stokes同时偏振成像系统的核心光学元件, 波片相位延迟误差对系统偏振测量精度有着不可忽略的影响, 为实现自主研制高偏振测量精度的分振幅型全Stokes同时偏振成像系统, 分析系统中波片相位延迟误差对偏振测量精度的影响已成为急需解决的关键问题.
本文介绍了分振幅型全Stokes同时偏振成像系统的工作原理, 建立了包含HWP和QWP相位延迟误差的Stokes矢量测量误差方程; 仿真了HWP和QWP相位延迟误差对自然光、0°/45°线偏光、左旋圆偏光等7种典型基态入射光的Stokes矢量测量误差的影响, 推导了任意入射光Stokes矢量测量误差的表征方法; 利用邦加球(Poincaré sphere)采样法仿真了不同偏振度的Stokes矢量作为入射光时, HWP和QWP相位延迟误差对偏振测量精度的影响, 为保证系统的偏振测量精度在2%内, 给出了波片相位延迟的加工精度要求, 这对系统参数设计、原理样机的研制具有重要的指导意义.

2

-->

2.1.偏振成像系统工作原理

分振幅型全Stokes同时偏振成像系统原理如图1所示, 入射光通过光学镜头后, 被部分偏振分束器(partial polarization beam splitter, PPBS)分成两路, 一路是从PPBS透射的光束, 被HWP调制后由PBS1分离为P偏振光和S偏振光, 从PBS1透射的P偏振光由CCD1接收, 反射的S偏振光由CCD2接收; 另一路是从PPBS反射的光束, 被QWP调制后由PBS2分离为P偏振光和S偏振光, 从PBS2透射的P偏振光由CCD3接收, 反射的S偏振光由CCD4接收.
图 1 分振幅型全Stokes同时偏振成像系统原理图
Figure1. Scheme of the division-of-amplitude full Stokes simultaneous polarization imaging system.

文献[20]研究了光学镜头的起偏效应. 以通常使用的玻璃进行分析, 当入射角很小(小于5°)时, 角度每相差1°, 造成的偏振度测量偏差为0.01%—0.03%. 本文研究的偏振成像系统是小视场成像系统, 入射主光束的入射方向与光学镜头法线的夹角很小, 所以不考虑光学镜头的起偏效应. 分振幅型全Stokes同时偏振成像系统4个通道的Mueller矩阵${{{M}}_{{\rm{path}}1}} $?$ {{{M}}_{{\rm{path}}4}}$为

$\left\{ \begin{aligned}& {{{M}}_{{\rm{path}}1}}{\rm{ = }}{{P}}\left( {1,0} \right) \cdot {{Ret}}\left( { - {\rm{22}}.{\rm{5}}^\circ ,180^\circ } \right) \cdot {{P}}\left( {{T_{\rm{p}}},{T_{\rm{s}}}} \right), \\[-1pt]& {{{M}}_{{\rm{path}}2}}{\rm{ = }}{{P}}\left( {0,1} \right) \cdot {{Ret}}\left( { - {\rm{22}}.{\rm{5}}^\circ ,180^\circ } \right) \cdot {{P}}\left( {{T_{\rm{p}}},{T_{\rm{s}}}} \right), \\[-1pt]& {{{M}}_{{\rm{path}}3}}{\rm{ = }}{{P}}\left( {1,0} \right) \cdot {{Ret}}\left( {{\rm{45}}^\circ ,90^\circ } \right) \cdot {{P}}\left( {1 - {T_{\rm{p}}},1 - {T_{\rm{s}}}} \right), \\[-1pt] & {{{M}}_{{\rm{path}}4}}{\rm{ = }}{{P}}\left( {0,1} \right) \cdot {{Ret}}\left( {{\rm{45}}^\circ ,90^\circ } \right) \cdot {{P}}\left( {1 - {T_{\rm{p}}},1 - {T_{\rm{s}}}} \right), \end{aligned} \right.$

${{Ret}}\left( {\theta ,\varphi } \right){\rm{ = }}\left[ {\begin{array}{*{20}{c}} 1&0&0&0 \\[-1pt] 0&\!\!\!\!{{{\cos }^2}\left(2\theta\right) + {{\sin }^2}\left(2\theta\right) \cos \varphi }&\!\!\!\!\!\!\!\!\!\!{\sin \left(2\theta\right) \cos \left(2\theta\right) - \sin \left(2\theta\right) \cos \left(2\theta\right) \cos \varphi }&\!\!\!\!{\sin \varphi \sin \left(2\theta\right) } \\[-1pt] 0&\!\!\!\!{\sin \left(2\theta\right) \cos \left(2\theta\right) - \sin \left(2\theta\right) \cos \left(2\theta\right) \cos \varphi }&\!\!\!\!\!\!\!\!\!\!{{{\sin }^2}\left(2\theta\right) + {{\cos }^2}\left(2\theta\right) \cos \varphi }&\!\!\!\!{ - \sin \varphi \cos \left(2\theta\right)} \\[-1pt] 0&\!\!\!\!{ - \sin \varphi \sin \left(2\theta\right)}&\!\!\!\!\!\!\!\!\!\!{\sin \varphi \cos \left(2\theta\right)}&\!\!\!\!{\cos \varphi } \end{array}} \right],$

$\begin{aligned}&{{P}}\left( {{T_{\rm{p}}},{T_{\rm{s}}}} \right) =\\&\frac{1}{2}\left[ {\begin{array}{*{20}{c}} {{T_{\rm{p}}} + {T_{\rm{s}}}}&{{T_{\rm{p}}} - {T_{\rm{s}}}}&0&0 \\ {{T_{\rm{p}}} - {T_{\rm{s}}}}&{{T_{\rm{p}}} + {T_{\rm{s}}}}&0&0 \\ 0&0&{2\sqrt {{T_{\rm{p}}}{T_{\rm{s}}}} }&0 \\ 0&0&0&{2\sqrt {{T_{\rm{p}}}{T_{\rm{s}}}} } \end{array}} \right].\end{aligned}$

(1)—(3)式中, ${T_{\rm{p}}}$和${T_{\rm{s}}}$分别为PPBS的p光透射系数和s光透射系数, ${{P}}\left( {{T_{\rm{p}}},{T_{\rm{s}}}} \right)$和${{P}}(1-{T_{\rm{p}}},$$1-{T_{\rm{s}}})$分别为PPBS透射和反射时的Mueller矩阵; ${{P}}\left( {1,0} \right)$为PBS1和PBS2透射时的Mueller矩阵, ${{P}}\left( {0,1} \right)$为PBS1和PBS2反射时的Mueller矩阵; ${{Ret}}\left( {\theta ,\varphi } \right)$为波片(相位延迟量$\varphi $)快轴与x轴正向夹角为$\theta $时的Mueller矩阵.
入射光Stokes矢量${{{S}}^\dagger }_{{\rm{in}}} = {[{{S}}_0^\dagger ,{{S}}_1^\dagger ,{{S}}_2^\dagger ,{{S}}_3^\dagger ]^{\rm{T}}}$, 系统的测量矩阵为M, 则CCD1-CCD4中对应像元的光强度矢量${{I}} = {[{{{I}}_{{\rm{path1}}}},{{{I}}_{{\rm{path2}}}},{{{I}}_{{\rm{path}}3}},{{{I}}_{{\rm{path4}}}}]^{\rm{T}}}$为

$\!\!\!\begin{array}{l}{{I}}=\left[\!\!\!{\begin{array}{*{20}{c}}{{{{M}}_{{\rm{path}}1}}(1,1)}&\!{{{{M}}_{{\rm{path}}1}}(1,2)}&\!{{{{M}}_{{\rm{path}}1}}(1,3)}\!\!\!\!\!\!\!\!\!\\{{{{M}}_{{\rm{path2}}}}(1,1)}&\!{{{{M}}_{{\rm{path2}}}}(1,2)}&\!{{{{M}}_{{\rm{path2}}}}(1,3)}\!\!\!\!\!\!\!\!\!\\{{{{M}}_{{\rm{path3}}}}(1,1)}&\!{{{{M}}_{{\rm{path3}}}}(1,2)}&\!{{{{M}}_{{\rm{path3}}}}(1,3)}\!\!\!\!\!\!\!\!\!\\{{{{M}}_{{\rm{path4}}}}(1,1)}&\!{{{{M}}_{{\rm{path4}}}}(1,2)}&\!{{{{M}}_{{\rm{path4}}}}(1,3)}\!\!\!\!\!\!\!\!\!\end{array}} \right.\\\left. {\begin{array}{*{20}{c}}&{{{{M}}_{{\rm{path}}1}}(1,4)}\\&{{{{M}}_{{\rm{path2}}}}(1,4)}\\&{{{{M}}_{{\rm{path3}}}}(1,4)}\\&{{{{M}}_{{\rm{path4}}}}(1,4)}\end{array}} \right] \cdot \left[ {\begin{array}{*{20}{c}}{{{S}}_0^\dagger }\\{{{S}}_1^\dagger }\\{{{S}}_2^\dagger }\\{{{S}}_3^\dagger }\end{array}} \right] = {{M}} \cdot {{{S}}^\dagger }_{{\rm{in}}}.\end{array}$

若已知测量矩阵M和光强度矢量${{I}}$, 则可获得入射光Stokes矢量:

${{{S}}^\dagger }_{{\rm{in}}} = {{{M}}^{{\rm{ - }}1}} \cdot {{I}} = {{R}} \cdot {{I}}.$

2 -->

2.2.波片相位延迟误差建模

PPBS的分束比为T_p/T_s, HWP和QWP相位延迟误差分别为$\sigma $和$\delta $时, 分振幅型全Stokes同时偏振成像系统4个通道的Mueller矩阵$ {{M}}_{{\rm{path}}1}^\dagger$?$ {{M}}_{{\rm{path}}4}^\dagger $和系统测量矩阵${{{M}}^\dagger }$分别如(6)式和(7)式所示:

$\left\{ \begin{aligned}& {{M}}_{{\rm{path}}1}^\dagger ={{P}}\left( {1,0} \right) \cdot {{Ret}}\left( { - {\rm{22}}.{\rm{5}}^\circ ,\left( {180^\circ + \sigma } \right)} \right) \cdot {{P}}\left( {{T_{\rm{p}}},{T_{\rm{s}}}} \right), \\[-1pt]& {{M}}_{{\rm{path2}}}^\dagger = {{P}}\left( {0,1} \right) \cdot {{Ret}}\left( { - {\rm{22}}.{\rm{5}}^\circ ,\left( {180^\circ + \sigma } \right)} \right) \cdot {{P}}\left( {{T_{\rm{p}}},{T_{\rm{s}}}} \right), \\[-1pt]& {{M}}_{{\rm{path3}}}^\dagger = {{P}}\left( {1,0} \right) \cdot {{Ret}}\left( {{\rm{45}}^\circ ,\left( {90^\circ + \delta } \right)} \right) \cdot {{P}}\left( {1{\rm{ - }}{T_{\rm{p}}},1-{T_{\rm{s}}}} \right), \\[-1pt]& {{M}}_{{\rm{path4}}}^\dagger = {{P}}\left( {0,1} \right) \cdot {{Ret}}\left( {{\rm{45}}^\circ ,\left( {90^\circ + \delta } \right)} \right) \cdot {{P}}\left( {1-{T_{\rm{p}}},1-{T_{\rm{s}}}} \right),\end{aligned} \right.$

$\begin{aligned} {{{M}}^\dagger } & = \left[ {\begin{array}{*{20}{c}} {{{M}}_{{\rm{path}}1}^\dagger (1,1)}&{{{M}}_{{\rm{path}}1}^\dagger (1,2)}&{{{M}}_{{\rm{path}}1}^\dagger (1,3)}&{{{M}}_{{\rm{path}}1}^\dagger (1,4)} \\ {{{M}}_{{\rm{path}}2}^\dagger (1,1)}&{{{M}}_{{\rm{path}}2}^\dagger (1,2)}&{{{M}}_{{\rm{path}}2}^\dagger (1,3)}&{{{M}}_{{\rm{path}}2}^\dagger (1,4)} \\ {{{M}}_{{\rm{path}}3}^\dagger (1,1)}&{{{M}}_{{\rm{path}}3}^\dagger (1,2)}&{{{M}}_{{\rm{path}}3}^\dagger (1,3)}&{{{M}}_{{\rm{path}}3}^\dagger (1,4)} \\ {{{M}}_{{\rm{path}}4}^\dagger (1,1)}&{{{M}}_{{\rm{path}}4}^\dagger (1,2)}&{{{M}}_{{\rm{path}}4}^\dagger (1,3)}&{{{M}}_{{\rm{path}}4}^\dagger (1,4)} \end{array}} \right] \\&= \tiny \left[ {\begin{array}{*{20}{c}}{\dfrac{{\left( {3{T_{\rm{p}}} + {T_{\rm{s}}}} \right) - \cos \left( \sigma \right) \cdot \left( {{T_{\rm{p}}} - {T_{\rm{s}}}} \right)}}{8}}&{\dfrac{{\left( {3{T_{\rm{p}}} - {T_{\rm{s}}}} \right) - \cos \left( \sigma \right) \cdot \left( {{T_{\rm{p}}} + {T_{\rm{s}}}} \right)}}{8}}&{\dfrac{{ - \sqrt {{T_{\rm{p}}} \cdot {T_{\rm{s}}}} \cdot \left( {\cos \left( \sigma \right) + 1} \right)}}{4}}&{\sqrt {\dfrac{{{T_{\rm{p}}} \cdot {T_{\rm{s}}}}}{8}} \sin \left( \sigma \right)} \\ {\dfrac{{\left( {{T_{\rm{p}}} + 3{T_{\rm{s}}}} \right) + \cos \left( \sigma \right) \cdot \left( {{T_{\rm{p}}} - {T_{\rm{s}}}} \right)}}{8}}&{\dfrac{{\left( {{T_{\rm{p}}} - 3{T_{\rm{s}}}} \right) + \cos \left( \sigma \right) \cdot \left( {{T_{\rm{p}}} + {T_{\rm{s}}}} \right)}}{8}}&{\dfrac{{\sqrt {{T_{\rm{p}}} \cdot {T_{\rm{s}}}} \cdot \left( {\cos \left( \sigma \right) + 1} \right)}}{4}}&{ - \sqrt {\dfrac{{{T_{\rm{p}}} \cdot {T_{\rm{s}}}}}{8}} \sin \left( \sigma \right)} \\ {\dfrac{{2 - \left( {{T_{\rm{p}}} + {T_{\rm{s}}}} \right) + \sin \left( \delta \right) \cdot \left( {{T_{\rm{p}}} - {T_{\rm{s}}}} \right)}}{4}}&{\dfrac{{\left( {{T_{\rm{s}}} - {T_{\rm{p}}}} \right) + \sin \left( \delta \right) \cdot \left( {{T_{\rm{p}}} + {T_{\rm{s}}} - 2} \right)}}{4}}&0&{\dfrac{{\sqrt {\left( {1 - {T_{\rm{p}}}} \right) \cdot \left( {1 - {T_{\rm{s}}}} \right)} \cdot \cos \left( \delta \right)}}{2}} \\ {\dfrac{{2 - \left( {{T_{\rm{p}}} + {T_{\rm{s}}}} \right) - \sin \left( \delta \right) \cdot \left( {{T_{\rm{p}}} - {T_{\rm{s}}}} \right)}}{4}}&{\dfrac{{\left( {{T_{\rm{s}}} - {T_{\rm{p}}}} \right) - \sin \left( \delta \right) \cdot \left( {{T_{\rm{p}}} + {T_{\rm{s}}} - 2} \right)}}{4}}&0&{\dfrac{{ - \sqrt {\left( {1 - {T_{\rm{p}}}} \right) \cdot \left( {1 - {T_{\rm{s}}}} \right)} \cdot \cos \left( \delta \right)}}{2}} \end{array}} \right] . \end{aligned} $

实际的测量矩阵是包含HWP和QWP相位延迟误差的测量矩阵${{{M}}^{}}^\dagger $, 则实际的入射光Stokes矢量${{{S}}_{{\rm{in}}}}$为

${{{S}}_{{\rm{in}}}} = {\left( {{{{M}}^{}}^\dagger } \right)^{ - 1}} \cdot {{I}}.$

由(5)式可得计算的入射光Stokes矢量${{{S}}^\dagger }_{{\rm{in}}}$, 则Stokes矢量的测量误差${{{e}}_{{S}}}$为

$\begin{aligned}{{{e}}_{{S}}} = & \left[ {\begin{array}{*{20}{c}} {{e_{{{S}}{\rm{0}}}}} \\ {{e_{{{S}}{\rm{1}}}}} \\ {{e_{{{S}}{\rm{2}}}}} \\ {{e_{{{S}}{\rm{3}}}}} \end{array}} \right] = {{S}}_{{\rm{in}}}^\dagger - {{{S}}_{{\rm{in}}}} = {{R}} \cdot {{{M}}^\dagger }\\ & \times {{{S}}_{{\rm{in}}}} - {{S}} _{{\rm{in}}} =\left( {{{R}} \cdot {{{M}}^\dagger } - {{E}}} \right) \cdot {{{S}}_{{\rm{in}}}},\end{aligned}$

(9)式中, ${e_{{{S}}i}}$表示Stokes矢量${{{S}}_i}$分量的测量误差(其中i = 0, 1, 2, 3), E为4 × 4单位矩阵, 矩阵R为理论测量矩阵M的逆矩阵.
对于入射光Stokes矢量, 一般关心的是光强的相对值^[21], 因此, 本文在接下来的波片相位延迟误差分析中, 将入射光Stokes矢量进行归一化处理.

本文将采用基于典型基态入射光和基于邦加球采样的两种方法对系统中波片相位延迟误差进行仿真和分析. 基于典型基态入射光的方法通过研究典型基态入射光的Stokes矢量测量误差的变化规律便可分析出任意偏振态的入射光Stokes矢量测量误差的变化规律. 邦加球采样的方法可用于更加全面、完备地评估波片相位延迟误差对系统偏振测量精度的影响.
2

-->

3.1.基于典型基态入射光的波片相位延迟误差仿真及分析

任意偏振态的入射光Stokes矢量都可分解成7种典型基态入射光(自然光、0°/90°/45°/135°线偏振光、右旋和左旋圆偏振光) Stokes矢量的线性组合. 通过分析7种典型基态入射光的Stokes矢量测量误差, 便可解析任意偏振态的入射光Stokes矢量测量误差.
归一化的任意偏振态的Stokes矢量${{{S}}_{{\rm{in}}}}$分解成自然光、0°/90°/45°/135°线偏振光、右旋和左旋圆偏振光Stokes矢量的线性组合表示为

$\begin{array}{c}{{{S}}_{{\rm{in}}}} = \left[\! {\begin{array}{*{20}{c}} 1\\{{{{{{S}}_1}} \mathord{\left/ {\vphantom {{{{{S}}_1}} {{{{S}}_0}}}} \right. \kern-\nulldelimiterspace} {{{{S}}_0}}}}\\{{{{{{S}}_2}} \mathord{\left/ {\vphantom {{{{{S}}_2}} {{{\bf{S}}_0}}}} \right. \kern-\nulldelimiterspace} {{{{S}}_0}}}}\\{{{{{{S}}_3}} \mathord{\left/ {\vphantom {{{{{S}}_3}} {{{{S}}_0}}}} \right. \kern-\nulldelimiterspace} {{{{S}}_0}}}}\end{array}}\! \right] = \left( {1 - {l_1} - {l_2} - {l_3} - {l_4} - {l_5} - {l_6}} \right)\left[ {\begin{array}{*{20}{c}}1\\0\\0\\0\end{array}} \right] + {l_1}\left[ {\begin{array}{*{20}{c}}1\\1\\0\\0\end{array}} \right] \\\;\;\;\;\;+ {l_2}\left[ \!\!{\begin{array}{*{20}{c}}\;\;\; 1\\{ - 1}\\\;\;\; 0\\\;\;\; 0\end{array}} \right] + {l_3}\left[ {\begin{array}{*{20}{c}}1\\0\\1\\0\end{array}} \right] + {l_4}\left[\!\! {\begin{array}{*{20}{c}}\;\;\; 1\\\;\;\; 0\\{ - 1}\\\;\;\; 0\end{array}} \right] + {l_5}\left[{\begin{array}{*{20}{c}}1\\0\\0\\1\end{array}} \right] + {l_6}\left[\!\! {\begin{array}{*{20}{c}}\;\;\; 1\\\;\;\; 0\\\;\;\; 0\\{ - 1}\end{array}} \right],\end{array}$

(10)式中, ${l_{2i - 1}} \!-\! {l_{2i}} = {{{{{S}}_i}} / {{{{S}}_0}}}$且${l_{2i - 1}} \cdot {l_{2i}} \;{\rm{ = 0}}\left( {i = 1,}\right.$$\left. {2,3} \right)$.
由(9)式可分别求出自然光、0°/90°/45°/135°线偏振光、右旋和左旋圆偏振光Stokes矢量的测量误差${{e}}_{{S}}^{(0)}$?$ {{e}}_{{S}}^{(6)}$, 结合(10)式可得任意偏振态的入射光Stokes矢量测量误差${{ e}_{{S}}}$:

${{{e}}_{{S}}} = \left( {1 - {l_1} - {l_2} - {l_3} - {l_4} - {l_5} - {l_6}} \right) \cdot {{e}}_{{S}}^{(0)} + {l_1} \cdot {{e}}_{{S}}^{(1)} + {l_2} \cdot {{e}}_{{S}}^{(2)} + {l_3} \cdot {{e}}_{{S}}^{(3)} + {l_4} \cdot {{e}}_{{S}}^{(4)} + {l_5} \cdot {{e}}_{{S}}^{(5)} + {l_6} \cdot {{e}}_{{S}}^{(6)}.$

表1为分振幅型全Stokes同时偏振成像系统的设计参数, 以其为仿真参数, 则入射光分别为自然光、0°/90°/45°/135°线偏振光、右旋和左旋圆偏振光时, HWP和QWP相位延迟误差对入射光Stokes参数测量误差的影响如图2所示. 其中 x轴和y轴分别为HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $(单位: (°)), z轴为入射光Stokes参数测量误差.

参数名称	参数值
1/2波片相位延迟量 (Retardance of HWP)	180°
1/4波片相位延迟量 (Retardance of QWP)	90°
1/2波片快轴方位角 (Fast axis orientation of HWP)	?22.5°
1/4波片快轴方位角 (Fast axis orientation of QWP)	45°
部分偏振分束器分束比 (Splitting ratio of PPBS)	Tp/Ts = 0.8/0.2

表1分振幅型全Stokes同时偏振成像系统设计参数
Table1.Parameters of division-of-amplitude full Stokes simultaneous polarization imaging system.

图 2 不同入射光情况下的Stokes参数测量误差　(a) 自然光; (b) 0° 线偏光; (c) 90° 线偏光; (d) 45° 线偏光; (e) 135° 线偏光; (f) 右旋圆偏光; (g) 左旋圆偏光
Figure2. Errors of Stokes parameters of different incident light: (a) Unpolarized light; (b) 0° liner polarized light; (c) 90° liner polarized light; (d) 45° liner polarized light; (e) 135° liner polarized light; (f) right circularly polarized light; (g) left circularly polarized light.

由图2可知, 系统中同时存在HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $时, 7种典型基态入射光${{S}_0}$分量和${{S}_1}$分量的测量误差均为0, ${{S}_2}$分量的测量误差仅受$\sigma $的影响, ${{S}_3}$分量的测量误差仅受$\delta $的影响.
已知7种基态入射光Stokes矢量测量误差时, 利用(11)式便可解析出任意偏振态的入射光Stokes矢量测量误差. 由图2和(11)式分析可知, 系统中同时存在HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $时, 任意偏振态的入射光的${{S}_0}$分量和${{S}_1}$分量的测量误差均为0, 任意偏振态的入射光${{S}_2}$分量的测量误差仅受$\sigma $的影响, 任意偏振态的入射光${{S}_3}$分量的测量误差仅受$\delta $影响.
2 -->

3.2.基于邦加球采样的波片相位延迟误差仿真与分析

由(9)式可知, Stokes矢量测量误差随着入射光Stokes矢量的变化而变化, 为了更加完备地评估波片相位延迟误差对分振幅型全Stokes同时偏振成像系统的Stokes矢量测量精度的影响, 需要分析不同偏振态的Stokes矢量作为入射光时波片相位延迟误差对系统的Stokes矢量测量精度的影响.
本文采用在邦加球的球面上或球内均匀地选取不同偏振态的Stokes矢量作为入射光, 进而评估HWP和QWP相位延迟误差对系统的偏振测量精度的影响. 在邦加球球面上或球内均匀采样时的任意偏振态的Stokes矢量为

$\left\{ \begin{aligned}& {{{S}}_n}={\left[ {\begin{array}{*{20}{c}} 1&{{P}\cos \left( {2{\chi _n}} \right)\cos \left( {2{\psi _n}} \right)} \quad {{P}\cos \left( {2{\chi _n}} \right)\sin \left( {2{\psi _n}} \right)}&{{P}\sin \left( {2{\chi _n}} \right)} \end{array}} \right]^{\rm{T}}}, \\ & {{\chi _n}=\dfrac{{n \cdot {\text{π}}}}{{2 \cdot {N_\chi } \cdot {N_\psi }}} - \dfrac{{\text{π}}}{4};} \quad {{\psi _n} = \dfrac{{n \cdot {\text{π}}}}{{{N_\psi }}};}\quad {n = 1,2,\cdots,{N_\chi } \cdot {N_\psi };}\quad {{\psi _n} \in \left[ {0,{\text{π}}} \right];} \quad{{\chi _n} \in \left[ { - \dfrac{{\text{π}}}{4},\dfrac{{\text{π }}}{4}} \right],} \end{aligned} \right.$

式中${{{S}}_n}$为邦加球球面上或球内的第n个入射光采样点的归一化Stokes矢量, P为偏振度, ${\chi _n}$和${\psi _n}$分别为第n个入射光采样点${{{S}}_n}$的椭圆度角和偏振角, ${N_\chi }$表示从邦加球的南极到北极围绕球体采样的圆周个数, ${N_\psi }$表示每个采样圆周上的采样点数.
令(12)式中的P = 1, ${N_\chi }$ = 20, ${N_\psi }$ = 50, 则图3为邦加球球面上1000个不同偏振态的Stokes矢量的三维分布和Stokes参数的数值分布.
图 3 邦加球球面上选取1000个不同偏振态的Stokes矢量的　(a) 三维分布和(b) Stokes参数的数值分布
Figure3. (a) 3D distribution and (b) stokes parameters values of 1000 Stokes vectors different degrees of polarization selected on the Poincaré sphere.

以表1中的系统设计参数作为仿真参数, 以图3中的1000个Stokes矢量作为入射光采样点, 则图4(a)和图4(b)分别为仅存在1°的HWP相位延迟误差时和仅存在1°的QWP相位延迟误差时1000个入射光采样点的Stokes参数${{{S}}_0}$?${{{S}}_3}$测量误差. 图4中横坐标为入射光采样点从邦加球南极到北极的采样序数, 纵坐标为Stokes参数测量误差.
图 4 1000个邦加球球面上的入射光采样点的Stokes参数测量误差　(a) 仅存在1°的 HWP相位延迟误差; (b) 仅存在1°的QWP相位延迟误差
Figure4. The measurement errors of Stokes parameters of 1000 incident light sampling points selected on the Poincaré sphere is simulated: (a) There is only 1° phase delay error of HWP; (b) there is only 1° phase delay error of QWP in the system.

由图4(a)可知, HWP相位延迟误差$\sigma $仅影响入射光${{{S}}_2}$分量的测量精度, 从邦加球的南极到赤道的采样过程中(即采样序数n从1到500的过程中), ${{{S}}_2}$分量测量误差的绝对值逐渐变小; 从邦加球的赤道到北极的采样过程中(即采样序数n从500到1000的过程中), ${{{S}}_2}$分量测量误差的绝对值逐渐变大. 由此可知, 入射光的椭圆度角的绝对值越小, 则${{{S}}_2}$分量测量误差的绝对值越小.
由图4(b)可知, QWP相位延迟误差$\delta $仅影响入射光${{{S}}_3}$分量的测量精度, 从邦加球的南极到赤道的采样过程中(即采样序数n从1到500的过程中), ${{{S}}_3}$分量测量误差的振幅逐渐变大; 从邦加球的赤道到北极的采样过程中(即采样序数n从500到1000的过程中), ${{{S}}_3}$分量测量误差的振幅逐渐变小. 由此可知, 入射光的椭圆度角的绝对值越小, ${{{S}}_3}$分量测量误差的振幅越大; 入射光的椭圆度角恒定时, ${{{S}}_3}$分量的测量误差随入射光的偏振角变化而变化(在邦加球赤道处, 即入射光为线偏振光时, ${{{S}}_3}$分量的测量误差受入射光偏振角的影响最大).
HWP和QWP相位延迟误差分别为$\sigma $和$\delta $时, 入射光采样点(采样总数为${N_\chi } \cdot {N_\psi }$)的偏振度均为P时, 在${N_\chi } \cdot {N_\psi }$个入射光采样点中${{S}}_i$分量测量误差绝对值的最大值$\max \_{{{e}}_{{{S}}i}}\left( {\sigma ,\delta ,P} \right)$(其中 i = 0, 1, 2, 3)为

$\left\{ \begin{aligned}& \max \_{{{e}}_{{{S}}i}}\left( {\sigma ,\delta ,P} \right) = \max \left( {\left| {{e_{Si}}\left( {\sigma ,\delta ,P,n} \right)} \right|} \right), \\& i = 0,1,2,3;\;n = 1,2,\cdots,{N_\chi } \cdot {N_\psi },\end{aligned} \right.$

其中${e_{Si}}\left( {\sigma ,\delta ,P,n} \right)$表示HWP相位延迟误差为$\sigma $、QWP相位延迟误差为$\delta $且入射光采样点(采样总数为${N_\chi } \cdot {N_\psi }$)的偏振度均为P时第n个入射光采样点的${{{S}}_i}$分量的测量误差. 与(13)式对应的偏振测量精度${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P} \right)$定义为

$\begin{array}{c} {\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P} \right) = \max \left( {\max \_{e_{Si}}\left( {\sigma ,\delta ,P} \right)} \right) \times 100\% ,\\ i = 0,1,2,3.\end{array}$

由图4可知, 仅存在1° 的HWP相位延迟误差且入射光采样点均为完全偏振光时, 偏振测量精度${\rm{acc}}\_{{S}} \left( {\sigma ={1^ \circ },\delta =0,P = 1} \right)$为1.23%; 仅存在1° 的QWP相位延迟误差且入射光采样点均为完全偏振光时, 偏振测量精度${\rm{acc}}\_{{S}} \, (\sigma =0,\delta ={1^ \circ }, $ $P = 1) $为3.48%.
以表1中的系统设计参数作为仿真参数, 分别令(12)式中的偏振度$P$为1.0, 0.8, 0.5, 0.2和0.1, 不同偏振度对应的采样总数均为1000 (即(12)式中恒有${N_\chi } $ = 20, ${N_\psi } $ = 50), 则当入射光采样点分别为邦加球球面上的完全偏振光($P$ = 1)和邦加球球内的部分偏振光($P$为0.8, 0.5, 0.2和0.1)时, HWP相位延迟误差$\sigma $或QWP相位延迟误差$\delta $对偏振测量精度${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P} \right)$影响的仿真结果如图5所示.
图 5 不同偏振度的采样点作为入射光时对偏振测量精度的影响　(a) 仅HWP相位延迟误差; (b) 仅QWP相位延迟误差; (c) HWP和QWP相位延迟耦合误差
Figure5. When the sampling points with different degrees of polarization are used as incident light, the effect of measurement accuracy: (a) The phase delay error of the HWP; (b) the phase delay error of the QWP; (c) the phase delay errors of the HWP and the QWP on polarization.

在图5(a)中, 横坐标为HWP相位延迟误差$\sigma $(单位为 (°)), 纵坐标为偏振测量精度${\rm{acc}}\_{{S}} ( \sigma ,\delta =$$ 0,P )$, 当偏振度P恒定时, ${\rm{acc}}\_{{S}}\left ( {\sigma ,\delta = 0,P} \right)$的值随${\rm{|}}\sigma {\rm{|}}$增大而增大; 当$\sigma $恒定时, ${\rm{acc}}\_{{S}}\left( {\sigma ,\delta = 0,P} \right)$的值随偏振度P增大而增大.
在图5(b)中, 横坐标为QWP相位延迟误差$\delta $(单位为 (°)), 纵坐标为偏振测量精度${\rm{acc}}\_{{S}}( \sigma = $$0,\delta ,P)$, 当偏振度P恒定时, ${\rm{acc}}\_{{S}}\left( {\sigma = 0,\delta ,P} \right)$的值随$\left| \delta \right|$增大而增大; 当$\delta $恒定时, ${\rm{acc}}\_{{S}}( \sigma = $$0,\delta ,P)$的值随偏振度P增大而增大.
图5(c)为HWP和QWP相位延迟耦合误差对偏振测量精度${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P} \right)$的影响, $\sigma $和$\delta $恒定时, ${\rm{acc}}\_{{S}}\left( {\sigma,\delta ,P} \right)$的值随偏振度P增大而增大.
由图5可知, HWP相位延迟误差$\sigma $或QWP相位延迟误差$\delta $恒定时, 入射光采样点的偏振度P = 1时对应的偏振测量精度${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P = 1} \right)$为最大值, 因此, 采用${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P = 1} \right)$来评估分振幅型全Stokes同时偏振成像系统的偏振测量精度. 系统中同时存在HWP和QWP相位延迟误差时, 系统的偏振测量精度${\rm{acc}}\_{{{S}}^{\left( {\sigma ,\delta } \right)}}$定义为

${\rm acc}\_{{{S}}^{\left( {\sigma ,\delta } \right)}}={\rm acc}\_{{S}}\left( {\sigma ,\delta ,P=1} \right).$

2 -->

3.3.波片相位延迟误差对偏振度测量精度的影响

偏振度的测量精度是评价分振幅型全Stokes同时偏振成像系统性能的重要指标, 入射光的Stokes矢量${{ S}_{{\rm{in}}}}{\rm{ = [}}{S_0},{S_1},{S_2},{S_3}{]^{\rm{T}}}$时, 偏振度$P = $$ {{\sqrt {{S_1}^2 + {S_2}^2 + {S_3}^2} } / {{S_0}}}$. 系统中同时存在HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $时, 入射光的偏振度为P时的偏振度测量误差$\Delta {P^{(\sigma ,\delta ,P)}}$为

$\Delta {P^{(\sigma ,\delta ,P)}} = \sum\limits_{i = 0}^3 {\frac{{\partial P}}{{\partial {S_i}}} \cdot \Delta {S_i}} = \sum\limits_{i = 0}^3 {\frac{{\partial P}}{{\partial {S_i}}} \cdot {e_{Si}}} .$

由图2和(11)式分析可知, 系统中同时存在HWP和QWP相位延迟误差时, 任意偏振态的入射光${S_0}$分量的测量误差${e_{S0}}$和${S_1}$分量的测量误差${e_{S1}}$均为0. 因此(16)式简化为

$\begin{aligned}[b]\Delta {P^{(\sigma ,\delta ,P)}} & = \frac{{\partial P}}{{\partial {S_2}}} \cdot {e_{S2}}{\rm{ + }}\frac{{\partial P}}{{\partial {S_3}}} \cdot {e_{S3}} \\& = \frac{{{S_2} \cdot {e_{S2}}{\rm{ + }}{S_3} \cdot {e_{S3}}}}{{{S_0} \cdot \sqrt {{S_1}^2 + {S_2}^2 + {S_3}^2} }}.\end{aligned}$

以表1中的系统设计参数作为仿真参数, 以图3中的1000个邦加球球面上(偏振度P均为1)的Stokes矢量作为入射光采样点, 则图6(a)和图6(b)分别为仅存在1° 的HWP相位延迟误差时和仅存在1° 的QWP相位延迟误差时1000个邦加球球面上的入射光采样点的偏振度测量误差.
图 6 1000个邦加球球面上的入射光采样点的偏振度测量误差　(a) 仅存在1°的HWP相位延迟误差; (b) 仅存在1°的QWP相位延迟误差
Figure6. The measurement errors of DOP of 1000 incident light sampling points selected on the Poincaré sphere is simulated: (a) There is only 1° phase delay error of HWP; (b) there is only 1° phase delay error of QWP in the system.

由图6(a)可知, 系统中仅存在HWP相位延迟误差$\sigma $时, 入射光从邦加球的南极到赤道的采样过程中(或者入射光从邦加球的赤道到北极的采样过程中), 偏振度测量误差$\Delta {P^{(\sigma ={1^ \circ },\delta =0,P = 1)}}$的振幅先增大后减小, 由此可知, 入射光为椭圆偏振光时, $\Delta {P^{(\sigma ={1^ \circ },\delta =0,P = 1)}}$随入射光偏振角的变化而变化(入射光的椭圆度角为 ±45° 时, $\Delta {P^{(\sigma ={1^ \circ },\delta =0,P = 1)}}$受入射光偏振角的影响最大). 由图6(b)可知, 系统中仅存在QWP相位延迟误差$\delta $时, 入射光从邦加球的南极到赤道的采样过程中(或者入射光从邦加球的赤道到北极的采样过程中), 偏振度测量误差$\Delta {P^{(\sigma = 0,\delta ={1^ \circ },P = 1)}}$的振幅先增大后减小, 入射光为椭圆偏振光时, $\Delta {P^{(\sigma = 0,\delta ={1^ \circ },P = 1)}}$随入射光偏振角的变化而变化.
系统中存在HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $时, 令(12)式中的${N_\chi }$ = 20, ${N_\psi }$ = 50, 对于偏振度均为P的1000个入射光采样点, 选择偏振度测量误差绝对值的最大值作为入射光偏振度为P时的偏振度测量精度 ${\rm{acc}}\_P (\sigma ,\delta ,P)$. 由图6可知, 系统中仅存在1°的HWP相位延迟误差时, 入射光偏振度为1时的偏振度测量精度${\rm{acc}}\_P (\sigma ={1^ \circ },\delta =0,P=1)$为0.61%; 系统中仅存在1°的QWP相位延迟误差时, 入射光偏振度为1时的偏振度测量精度${\rm{acc}}\_P (\sigma {\rm{ = }}0, \delta {\rm{ = }}{1^ \circ }, P{\rm{ = }}1)$为2.08%.
以表1中的系统设计参数作为仿真参数, 分别令(12)式中偏振度$P$为1, 0.8, 0.5, 0.2和0.1, 不同偏振度对应的入射光采样总数均为1000 (即(12)式中恒有${N_\chi }$ = 20, ${N_\psi }$ = 50), HWP相位延迟误差$\sigma \in \left( { - {2^ \circ },{2^ \circ }} \right)$且QWP相位延迟误差$\delta \in \left( { - {2^ \circ },{2^ \circ }} \right)$时, 则入射光的偏振度$P$分别为1, 0.8, 0.5, 0.2和0.1时的偏振度测量精度${\rm{acc}}\_P (\sigma ,\delta ,P)$的仿真结果如图7所示. 由图7可知, $\sigma $和$\delta $恒定时, 入射光的偏振度P越大, 偏振度测量精度${\rm{acc}}\_P (\sigma ,\delta ,P)$的值越大, 因此选择入射光的偏振度P = 1时对应的偏振度测量精度${\rm{acc}}\_P (\sigma ,\delta ,P = 1)$来评估系统的偏振度测量精度.
图 7 入射光的偏振度$P$分别为1.0, 0.8, 0.5, 0.2和0.1时, 偏振度测量精度${\rm{acc}}\_P (\sigma ,\delta ,P)$随HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $的变化关系
Figure7. Variation relation of measure accuracy ${\rm{acc}}\_P (\sigma ,\delta ,P)$ of DOP with the phase delay error of HWP and the phase delay error of QWP under the condition of $P$ = 1.0, 0.8, 0.5, 0.2 and 0.1.

以表1中的系统设计参数作为仿真参数, 以图3中1000个偏振度P均为1的Stokes矢量作为入射光采样点, 则表2和表3分别为系统的偏振度测量精度${\rm{acc}}\_P(\sigma ,\delta ,P{\rm{ = }}1)$和系统偏振测量精度${\rm{acc}}\_{{{S}}^{\left( {\sigma ,\delta } \right)}}$随HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $的变化关系.

$\sigma $	$\delta$
	?1.0°	?0.9°	?0.6°	?0.5°	0°	0.5°	0.6°	0.9°	1.0°
?3.2°	3.22%	3.07%	2.65%	2.53%	2.03%	2.54%	2.68%	3.11%	3.26%
?3.1°	3.17%	3.02%	2.59%	2.47%	1.97%	2.48%	2.62%	3.05%	3.20%
?1.7°	2.50%	2.32%	1.82%	1.66%	1.06%	1.67%	1.83%	2.33%	2.50%
?1.6°	2.46%	2.28%	1.77%	1.61%	1.00%	1.62%	1.78%	2.28%	2.46%
?0.5°	2.13%	1.92%	1.31%	1.12%	0.31%	1.12%	1.31%	1.92%	2.13%
0°	2.09%	1.88%	1.25%	1.04%	0	1.04%	1.25%	1.88%	2.09%
0.5°	2.13%	1.92%	1.32%	1.12%	0.31%	1.12%	1.32%	1.92%	2.13%
1.6°	2.45%	2.27%	1.77%	1.62%	1.00%	1.61%	1.77%	2.27%	2.45%
1.7°	2.49%	2.31%	1.82%	1.67%	1.06%	1.66%	1.82%	2.31%	2.49%
3.1°	3.20%	3.05%	2.64%	2.50%	1.95%	2.46%	2.59%	3.02%	3.17%
3.2°	3.25%	3.11%	2.70%	2.57%	2.02%	2.52%	2.65%	3.07%	3.22%

表2系统偏振度测量精度${\rm{acc}}\_P(\sigma ,\delta ,P {\rm{ = }} 1)$随HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $的变化关系
Table2.Variation relation of measure accuracy ${\rm{acc}}\_P(\sigma ,\delta ,P {\rm{ = }} 1)$ of DOP with the phase delay error $\sigma $ of HWP and the phase delay error $\delta $ of QWP.

$\sigma$	$\delta $
	?1.0°	?0.9°	?0.6°	?0.5°	0°	0.5°	0.6°	0.9°	1.0°
?3.2°	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%
?3.1°	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%
?1.7°	3.48%	3.14%	2.11%	2.11%	2.11%	2.11%	2.11%	3.14%	3.48%
?1.6°	3.48%	3.14%	2.09%	1.99%	1.99%	1.99%	2.09%	3.14%	3.48%
?0.5°	3.48%	3.14%	2.09%	1.74%	0.62%	1.74%	2.09%	3.14%	3.48%
0°	3.48%	3.14%	2.09%	1.74%	0	1.74%	2.09%	3.14%	3.48%
0.5°	3.48%	3.14%	2.09%	1.74%	0.62%	1.74%	2.09%	3.14%	3.48%
1.6°	3.48%	3.14%	2.09%	1.99%	1.99%	1.99%	2.09%	3.14%	3.48%
1.7°	3.48%	3.14%	2.11%	2.11%	2.11%	2.11%	2.11%	3.14%	3.48%
3.1°	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%	3.88%
3.2°	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%	4.01%

表3系统偏振测量精度${\rm{acc}}\_{{{S}}^{\left( {\sigma ,\delta } \right)}}$随HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $的变化关系
Table3.Variation relation of system polarization measurement accuracy ${\rm{acc}}\_{{{S}}^{\left( {\sigma ,\delta } \right)}}$ with the phase delay error $\sigma $ of HWP and the phase delay error $\delta $ of QWP.

为了满足气溶胶光学参数探测、目标和背景偏振探测等应用需求, 要求分振幅型同时偏振成像系统的偏振测量精度优于2%. 由表2和表3可知, 系统中同时存在HWP相位延迟误差$\sigma $和QWP相位延迟误差$\delta $时, 恒有${\rm{acc}}\_{{S}}^{\left( {\sigma ,\delta } \right)} \succ {\rm{acc}}\_$$ P (\sigma ,\delta ,P=1)$. 为保证系统偏振测量精度${\rm{acc}}\_{{{S}}^{\left( {\sigma ,\delta } \right)}}$和系统的偏振度测量精度${\rm{acc}}\_P (\sigma ,\delta ,P{\rm{ = }}1)$均在2%以内, 则HWP相位延迟误差应在 ±1.6°内, QWP相位延迟误差应在 ±0.5°内.

为了验证本文对分振幅型全Stokes偏振成像系统波片相位延迟误差分析结论的正确性, 搭建了如图8所示的实验光路, 实验光路由偏振态发生器(polarization state generator, PSG)、四分束偏振分析器(four-paths polarization analyzer)和陷阱探测器(trap detector)组成. 其中, PSG由He-Ne激光器、激光功率控制器(laser power controller, LPC)、退偏器、安装在电动转台上的偏振片和
图 8 波片相位延迟误差分析实验光路
Figure8. Experimental optical path of wave plate phase delay error analysis.

QWP组成; 四分束偏振分析器由PPBS, HWP, QWP, PBS1和PBS2组成.
光路中PSG和四分束偏振分析器的主要参数分别如表4和表5所列. 实验中四分束偏振分析器4个通道的Mueller矩阵${ M}_{\tau ,{\rm{path}}1}^{}$—${ M}_{\tau ,{\rm{path}}4}^{}$为

参数名称	参数值
He-Ne激光器输出波长	632.99 nm
He-Ne激光器光强稳定性	± 0.1%
电动转台旋转精度	0.005°
偏振片消光系数	≥ 10000∶1
零级QWP相位延迟量	89.87°@632.99 nm

表4实验光路中PSG的主要参数
Table4.Parameters of PSG in the experimental optical path.

参数名称	参数值
PPBS分束比${{T_{\rm p}^{({\rm PPBS})}} / {T_{\rm s}^{({\rm PPBS})}}}$	0.788/0.191
零级HWP相位延迟量	179.74°@632.99 nm
零级QWP相位延迟量	89.87°@632.99 nm
零级HWP快轴方位角	?22.5°
零级QWP快轴方位角	45°
PBS1分束比${{T_{\rm{p}}^{({\rm{PBS1}})}} / {T_{\rm{s}}^{({\rm{PBS1}})}}}$	0.981/0.0007
PBS2分束比${{T_{\rm{p}}^{({\rm{PBS2}})}} / {T_{\rm{s}}^{({\rm{PBS2}})}}}$	0.988/0.0008

表5实验光路中四分束偏振分析器的主要参数
Table5.Parameters of the four-paths polarization analyzer in the experimental optical path.

$\left\{ \begin{aligned}& {{M}}_{\tau ,{\rm{path}}1}^{}\left( {\sigma ,\delta } \right)={{P}}\left( {T_{\rm{p}}^{({\rm{PBS1}})},T_{\rm{s}}^{({\rm{PBS1}})}} \right) \cdot {{Ret}}\left( { - {\rm{22}}.{\rm{5}}^\circ ,\left( {180^\circ + \sigma } \right)} \right) \cdot {{P}}\left( {T_{\rm{p}}^{({\rm{PPBS}})},T_{\rm{s}}^{({\rm{PPBS}})}} \right), \\ & {{M}}_{\tau ,{\rm{path2}}}^{}\left( {\sigma ,\delta } \right)={{P}}\left( {1 - T_{\rm{p}}^{({\rm{PBS1}})},1 - T_{\rm{s}}^{({\rm{PBS1}})}} \right) \cdot {{Ret}}\left( { - {\rm{22}}.{\rm{5}}^\circ ,\left( {180^\circ + \sigma } \right)} \right) \cdot {{P}}\left( {T_{\rm{p}}^{({\rm{PPBS}})},T_{\rm{s}}^{({\rm{PPBS}})}} \right), \\ & {{M}}_{\tau ,{\rm{path3}}}^{}\left( {\sigma ,\delta } \right)={{P}}\left( {T_{\rm{p}}^{({\rm{PBS2}})},T_{\rm{s}}^{({\rm{PBS2}})}} \right) \cdot {{Ret}}\left( {{\rm{45}}^\circ ,\left( {90^\circ + \delta } \right)} \right) \cdot {{P}}\left( {1-T_{\rm{p}}^{({\rm{PPBS}})},1{\rm{ - }}T_{\rm{s}}^{({\rm{PPBS}})}} \right), \\ & {{M}}_{\tau ,{\rm{path4}}}^{}\left( {\sigma ,\delta } \right)={{P}}\left( {1 - T_{\rm{p}}^{({\rm{PBS2}})},1 - T_{\rm{s}}^{({\rm{PBS2}})}} \right) \cdot {{Ret}}\left( {{\rm{45}}^\circ ,\left( {90^\circ + \delta } \right)} \right) \cdot {{P}}\left( {1{\rm{ - }}T_{\rm{p}}^{({\rm{PPBS}})},1{\rm{ - }}T_{\rm{s}}^{({\rm{PPBS}})}} \right).\end{aligned} \right.$

将表5中PPBS, PBS1和PBS2的分束比代入(18)式, 并令$\sigma $ = ?0.26°和$\delta $ = ?0.13°, 得到四分束偏振分析器的实际测量矩阵${{{M}}_1}^\dagger $. 将表5中PPBS, PBS1和PBS2的分束比代入(18)式, 并令$\sigma $ = 0和$\delta $ = 0, 得到估算的测量矩阵${{{M}}_0}$.
实验中定义入射光的传播方向为z轴的正方向, 垂直于实验平台向上为y轴的正方向(坐标系满足右手法则). 利用PSG产生偏振度为1的标准偏振光, PSG起始状态为: QWP快轴与x轴的夹角为0°, 偏振片的透光轴与x轴的夹角为?45°. 将PSG中的QWP以步长$\omega $ = 6°进行逆时针旋转, PSG中的偏振片以步长$\varpi\; {\rm{ = }}\;\omega {\rm{ \;+\; 0}}{\rm{.}}{{\rm{3}}^ \circ }\;{\rm{ = }}\;{6.3^ \circ }$进行逆时针旋转, 同时旋转n次共产生300个不同偏振态的完全偏振光${\mathop {{S}} }_{{\rm{in}}}^{{\rm{PSG}}}(n)$(n = 1, 2, ···, 300), 相当于(12)式中的P = 1, ${N_\chi }$ = 10, ${N_\psi }$ = 30, 均匀分布在邦加球球面上的300个采样点.
利用陷阱探测器对4个通道的出射光强进行测量, 可得与${{S}}_{{\rm{in}}}^{{\rm{PSG}}}(n)$对应的光强度矢量${{ I}_1}(n){\rm{ = }}$${[{I_{1,{\rm{path1}}}}(n),{I_{1,{\rm{path}}2}}(n),{I_{1,{\rm{path}}3}}(n),{I_{1,{\rm{path4}}}}(n)]^{\rm{T}}}$(n = 1, 2, ···, 300), 与${{S}}_{{\rm{in}}}^{{\rm{PSG}}}(n)$对应的真实值为${{S}}_{{\rm{in}}}^{'{\rm{(1,real)}}}(n) =$$ {\left( {{{{M}}_1}^\dagger } \right)^{ - 1}} \cdot {{{I}}_1}(n)$(n = 1, 2, ···, 300), 归一化后的真实值为${{S}}_{{\rm{in}}}^{(1,{\rm{real}})}(n)$(n = 1, 2, ···, 300).
由(9)式可得入射光Stokes矢量的估算值为${{S}}_{{\rm{in}}}^{(1,{\rm{estimated}})}(n){\rm{ = }}{\left( {{{{M}}_0}} \right)^{{\rm{ - }}1}} \cdot {{{M}}_1}^\dagger \cdot {{S}}_{{\rm{in}}}^{{\rm{(1,real)}}}(n)$(n = 1, 2, ···, 300), 入射光Stokes矢量的测量误差为${{e}}_{{S}}^{(1)}(n) =$$ {{S}}_{{\rm{in}}}^{{\rm{(1,estimated)}}}(n) - {{S}}_{{\rm{in}}}^{{\rm{(1,real)}}}(n)$(n = 1, 2, ···, 300). 偏振度$P = {{\sqrt {{{{S}}_1}^2 + {{{S}}_2}^2 + {{{S}}_3}^2} } / {{{{S}}_0}}}$, 由${{S}}_{{\rm{in}}}^{(1,{\rm{real}})}(n)$和${{S}}_{{\rm{in}}}^{{\rm{(1,estimated)}}}(n)$可分别获得入射光的真实偏振度$P^{(1,{\rm{real}})}(n)$和估算偏振度$P^{(1,{\rm{estimated}})}(n)$, 偏振度的测量误差为$P^{(1,{\rm{error}})}(n) = P^{{\rm{(1,estimated)}}}(n)-$$ P^{{\rm{(1,real)}}}(n)$(n = 1, 2, ···, 300).
处理数据得, 光路中HWP相位延迟误差$\sigma $ = ?0.26°和QWP相位延迟误差$\delta $ = ?0.13°同时存在时, 入射光Stokes矢量和偏振度的测量结果如图9所示, 入射光从邦加球的南极到北极的采样过程中, ${{{S}}_0}$分量的测量误差和${{{S}}_1}$分量的测量误差为0 , ${{{S}}_2}$分量测量误差的绝对值先减小后增大, ${{{S}}_3}$分量测量误差的振幅先增大后减小, 与图4仿真结果的变化规律一致. 入射光从邦加球的南极到赤道的采样过程中(或从邦加球的赤道到北极的采样过程中), 偏振度测量误差的振幅先增大后减小, 与图6仿真结果的变化规律相符.
图 9 实验中HWP相位延迟误差$\sigma$ = ?0.26°和QWP相位延迟误差$\delta $ = ?0.13°时的测量结果 入射光　(a) ${{{S}}_0}$分量; (b) ${{{S}}_1}$分量; (c) ${{{S}}_2}$分量; (d) ${{{S}}_3}$分量; (e) 偏振度
Figure9. Measurement results: (a) ${{{S}}_0}$ component; (b) ${{{S}}_1}$ component; (c) ${{{S}}_2}$ component; (d) ${{{S}}_3}$ component; (e) DOP of the incident light under the condition of $\sigma $ = ?0.26° and $\delta$ = ?0.13°.

图9结果显示, 偏振测量精度${\rm{acc}}\_{{S}}(\sigma =$$ - {{0.26}^ \circ },\delta = - 0.13^ \circ,P=1)$为0.44%, 偏振度测量精度${\rm{acc}}\_P(\sigma = -{0.26^ \circ },\delta = - 0.13^ \circ, P=1)$为0.33%, 与表2和表3的仿真结果相符合, 即HWP相位延迟误差在 ±1.6°内, 且QWP相位延迟误差在 ±0.5°内时, 偏振测量精度和偏振度测量精度均小于2%.

本文针对分振幅型全Stokes同时偏振成像仪中HWP和QWP相位延迟误差对入射光Stokes矢量测量精度影响较大的问题, 建立了包含HWP和QWP相位延迟误差的Stokes矢量测量误差模型, 仿真了自然光、0°/90°/45°/135°线偏振光、右旋和左旋圆偏振光分别作为入射光时, HWP和QWP相位延迟耦合误差对Stokes矢量测量误差的影响, 给出了求解任意偏振态入射光的Stokes矢量测量误差的方法. 仿真结果表明, 同时存在HWP和QWP相位延迟误差时, 任意偏振态的入射光的${{{S}}_0}$分量和${{{S}}_1}$分量的测量误差均为0, 任意偏振态的入射光${{{S}}_2}$分量的测量误差仅受HWP相位延迟误差的影响, 任意偏振态入射光${{{S}}_3}$分量的测量误差仅受QWP相位延迟误差的影响. 为了更加完备地分析HWP和QWP相位延迟误差对系统偏振测量精度的影响, 本文提出分别在邦加球的球面和球内选取不同偏振度的Stokes矢量作为入射光, 在此基础上仿真了HWP和QWP相位延迟误差对偏振测量精度${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P} \right)$和偏振度测量精度${\rm{acc}}\_P(\sigma ,\delta ,P)$的影响. 仿真结果表明, HWP和QWP相位延迟误差恒定时, 偏振测量精度${\rm{acc}}\_{{S}}\left( {\sigma ,\delta ,P} \right)$的值和偏振度测量精度${\rm{acc}}\_P(\sigma ,\delta ,P)$的值都随着偏振度P的增大而增大, 因此, 选取入射光偏振度为1时对应的偏振测量精度${\rm{acc}}\_{{S}} $$\left( {\sigma ,\delta ,P = 1} \right)$和偏振度测量精度${\rm{acc}}\_P(\sigma ,\delta ,$$P=1)$来评估系统. 为保证系统偏振测量精度和偏振度测量精度均在2%内, HWP相位延迟误差应在 ±1.6°内, QWP相位延迟误差应在 ±0.5°内. 这对分振幅型全Stokes同时偏振成像系统的光学参数设计、波片相位延迟的加工精度要求和系统研制具有重要的意义.