基于Monte Carlo方法的散射介质对光束偏振态的影响分析.docVIP

基于Monte Carlo方法的散射介质对光束偏振态的影响分析 【摘 要】为了能准确分析光束的偏振特性在大气传输过程中所受到的影响，特别是在光束发生散射后含有决定大气光学特性的偏振信息这一目的。采用Stokes-Mueller方法处理光束偏振态，以及Monte Carlo方法对偏振光的传输特性进行建模分析。 【关键词】蒙特卡洛；偏振传输；散射 【Abstract】In order to accurate analysis the influence of polarized light beam through the atmosphere， especially the light beam which scattered with large amounts of information has decided the properties of atmospheric optical. The Stokes-Mueller method is adapted to fix the polarization state of light beam，and analysis the properties of polarized light transmit simulation with Monte Carlo method. Getting the polarization of incident light which is scattered. 【Key words】Monte Carlo；Polarization transmission；Scatter 0 引言 近些年，对偏振光束多次散射传输的模拟，主要采用的方法是Monte Carlo方法，Monte Carlo法则是通过大量随机试验来对散射传输的过程进行模拟。这样就能对任何环境、任意条件下的散射过程进行模拟，得到的结果也与通过实验方式所得的结果相同，这些正是其他方法所存在的瓶颈。当光束在大气的中传输时，许多因素影响着光束所携带的偏振信息，尤其是在与大气中的颗粒碰撞发生散射后，散射后的原本光束所携带偏振信息会发生极大的变化，又因为偏振光在散射介质中是参照某些规律传输的，因此分析偏振光大气传输的特性十分必要。 1 偏振光子散射的Monte Carlo模型 在用Monte Carlo法对偏振光束在均匀介质中传输的特性进行建模时，普遍是采用stokes矢量的三个宏观可以通过实验设备准确测得的独立参量来描述。它们是电矢量在坐标系中沿x，y轴的分量Ex，Ey和两振动分量在瞬间的相位差δ。所以当模拟偏振光子在大气中传输时，由以下的矩阵形式来表示stokes矢量： 其中，“”符号表示全部时间的平均值。stokes矢量的矩阵形式是四维列矩阵矢量，完全偏振光和部分偏振光都可以用它来表示，即可以表示包括偏振度在内的任意光的偏振态。在这个列矩阵中，I表示X轴和Y轴上的强度和，即光波的光强。Q表示X轴和Y轴的强度差，U表示+45°和-45°方向上的强度差，V表示右旋（或左旋）偏振光。完全偏振光则有如下的这种关系： 公式（2）说明了参量I、Q、U、V之间并不是完全独立的。 偏振光束在大气或不同介质中传输时，光束会与大气中的各种不均匀的颗粒、气溶胶等发生碰撞，产生多次散射，这些碰撞后的粒子会使入射光束的偏振态产生变化。为了完成对偏振光多次散射的模拟，我们将发生碰撞的散射不规则粒子均看作均匀粒子。而同样有一部分光子经过多次散后在射介质中被完全的吸收或溢出介质的参考平面。 2 介质中光子散射过程分析 模拟光子在介质中的多次散射模型其实质是以光子的单次散射为基础模型，建立如图1所示的多次散射模型图，散射粒子在空间中随机出现，同时粒子在空间中服从均匀分布。让入射光束从Z轴正方向进入，把XOZ 平面定义为参考平面，在Z轴上建立一个接收平面，光子在散射介质中发生多次碰撞，每当光子发生散射时，都要重新对stokes矢量进行更新以确保参考平面与入射光子在同一平面。统计所有经过接收平面光子的偏振态并计算其期望值，这样就可以得到接收光束的偏振度了[1]。 我们得到的偏振分量是以stokes矢量的形式得到的，光束散射结束后的总偏振度为： 其中的I、Q、U、V表示探测平面接收到的所有光子偏振分量的求和值[2]。 3 仿真结果验证 从图2中很清晰的发现各向同性介质的穆勒矩阵元素的散射强度分布图像都在某一程度上具有对称性、反对称性和轴对称性。其中对角线附近的元素散射强度分布特别接近。将左上角与右上角连接后，不但可以看出该矩阵有着轴对称性，而且还可以看出该矩阵还有着轴反对称性，元素m11和m44的图形形状相同均为饼状型但散射强度恰恰颠倒过来，同时元素m44的散射强度分布范围明显的略小于元素m11的散射强度分布范围，元