《Photonics》：Error Analysis of Optical System in Laser Warning Receivers Based on Monte Carlo Method
编辑推荐：
激光告警接收机（Laser Warning Receiver, LWR）通过测量入射激光辐射的方向来定位激光威胁源。其光学系统固有的参数误差是影响角度测量精度的关键因素。现有误差分析方法多局限于单因素研究或经验估计，缺乏对多源耦合误差的系统处理。为解决这一局限
激光告警接收机（Laser Warning Receiver, LWR）通过测量入射激光辐射的方向来定位激光威胁源。其光学系统固有的参数误差是影响角度测量精度的关键因素。现有误差分析方法多局限于单因素研究或经验估计，缺乏对多源耦合误差的系统处理。为解决这一局限，本文提出一种基于蒙特卡洛（Monte Carlo）模拟的LWR光学系统误差分析方法。对包括主点偏移、焦距误差、探测器焦平面倾斜、畸变以及光斑中心提取误差在内的主要误差源进行了系统分析，并建立了一个综合误差传播模型，该模型避免了繁琐的理论推导，直观地揭示了各误差因素对角度测量精度的影响规律与耦合效应。仿真结果表明，在120° × 120°视场（Field of View, FOV）内，综合误差的3σ值为4716.46角秒（arcsec）。实验验证表明，实测最大误差为4901.12 arcsec，与仿真预测偏差184.66 arcsec（约0.05°）。所提出的方法可为LWR的误差标定、公差优化及高精度制造提供理论基础和工程指导。
激光告警接收机（Laser Warning Receiver, LWR）是探测来袭激光辐射入射方向、波长和功率密度等关键参数的核心设备，在激光威胁告警与定位中具有重要作用。随着光电对抗技术的不断发展，对LWR角度测量精度和系统稳定性的要求日益提高。作为LWR的核心组成，光学系统的参数误差直接影响入射方向估计的准确性。然而，现有的误差研究大多局限于单一误差因素或经验性粗略估计，缺乏系统全面的误差传播模型。实际中，光学测量系统普遍面临多源耦合误差和强非线性问题，传统基于一阶泰勒展开的误差分析方法难以准确表征系统真实的不确定性。虽然自适应光学（Adaptive Optics, AO）波前传感和夏克-哈特曼波前传感器（Shack–Hartmann Wavefront Sensor, SHWFS）在光斑质心提取、动态范围扩展和非线性误差建模方面积累了经验，但LWR在大视场条件下的误差源在空间分布和耦合机制上与上述研究存在本质差异。因此，针对LWR大视场下多源耦合误差传播与量化问题开展系统研究具有重要意义。
研究人员提出了一种基于蒙特卡洛（Monte Carlo）模拟的LWR光学系统误差分析方法。该方法建立了从真实入射角到测量像点的完整前向映射模型，严格整合了主点偏移、焦距误差、探测器焦平面倾斜、畸变和光斑中心提取误差五大误差源，并采用逆模型由测量坐标估计入射角，定义了方位角和仰角测量误差。通过单因素与综合因素蒙特卡洛仿真，获得了各误差源单独及组合作用下的角度测量误差概率分布。仿真结果表明，主点偏移误差、畸变误差和焦距误差是影响角度测量精度的主导因素，探测器焦平面倾斜误差次之，光斑中心提取误差贡献可忽略。综合误差的3σ值为4716.46角秒（arcsec）。实验验证中，实测最大误差为4901.12 arcsec，与仿真预测偏差仅约0.05°，验证了所提模型的有效性。该研究可为LWR的误差标定、公差优化和高精度制造提供理论依据和工程指导，并具有推广至其他成像式光学测量系统的潜力。
研究采用的主要关键技术方法如下：蒙特卡洛随机模拟中，根据制造和装配公差设置各误差源的标准差（主点偏移σ=1.5 pixel，焦距误差σ=0.1/3 mm，倾斜σ=0.5/3°），畸变误差σ=2/3 pixel来自Optiland v0.6.0光学设计软件公差分析，光斑中心提取误差σ=0.1 pixel来自实验统计。仿真采用直接采样策略，输入误差向量各分量服从零均值高斯分布，以角度误差均值μ和标准差σ的相对变化作为收敛判据（δμ<0.5%且δσ<0.5%并连续5000次满足），在120°×120°视场（Field of View, FOV）内均匀随机生成真实入射角，共进行100,000次独立抽样。实验验证搭建了由LWR、1064 nm脉冲Nd:YAG激光器（CNI，长春，中国）、高精度两轴转台（华中测控，武汉，中国，定位精度优于0.001°）和信息处理计算机组成的测试系统，在暗室环境（22±2°C）中测量169个测试点，每点重复10次取平均。
研究结果如下：
4.1 基于蒙特卡洛的方法
该方法基于概率统计理论，通过大量独立重复随机试验的样本统计来近似目标量的数学期望或概率分布。本文将光学系统参数误差视为随机变量，通过前向成像模型耦合传播，最终得到角度测量误差的统计特性。
4.1.1 随机抽样策略
根据中心极限定理，当误差源相互独立且效应较弱时，组合误差趋于正态分布。但本系统中探测器焦平面倾斜和畸变引入强非线性，因此采用直接蒙特卡洛采样，对输入误差向量θ=(α,β,Δf,Δsx,Δsy,Δxe,Δye)的各分量进行零均值高斯分布随机抽样。
4.1.2 收敛判据
采用角度测量误差的均值μ和标准差σ作为双收敛监测量，其相对变化均小于0.5%且该状态持续5000次连续采样时视为收敛。结果表明，当采样次数N≥5×10
4时，均值μ和标准差σ的相对变化均低于0.5%，故选择N=10
5次以保证统计稳定性。
4.1.3 单因素与综合因素仿真流程
设计了两种仿真模式：单因素模式仅随机采样目标误差项，综合因素模式同时采样所有误差项。每次试验在视场内均匀随机生成真实入射角，共进行100,000次。综合仿真得到角度测量误差均值为1542.33 arcsec，标准差为1058.04 arcsec，3σ范围为[0, 4716.46] arcsec。单因素仿真表明，综合误差均值远低于各单因素误差均值线性叠加（2876.16 arcsec），说明误差源之间存在非线性耦合和部分相互抵消。按平方均值误差比例计算贡献率，主点偏移、畸变和焦距误差占比较高，探测器焦平面倾斜次之，光斑中心提取误差贡献可忽略。
4.2 最大误差法仿真
采用最大误差法验证蒙特卡洛分析。当主点偏移在?4.5 pixel至+4.5 pixel范围内变化时，LWR角度测量误差随仰角和偏移量绝对值的增大而单调增加；焦距误差在?0.1 mm至+0.1 mm范围内变化时，其影响随仰角增大而增大。探测器焦平面倾斜固定为0.5°时，其影响随入射角变化；畸变+2 pixel和?2 pixel时的模拟结果也显示了其对角度测量精度的影响。光斑中心提取误差±0.1 pixel造成的影响可忽略。最大误差法的预测结果与蒙特卡洛仿真结果一致，验证了误差传播模型的有效性。
4.3 实验验证
实验装置包括LWR、激光源、两轴转台和信息处理计算机。激光源为输出1064 nm、脉宽10 ns、重复频率200 Hz的脉冲Nd:YAG激光器，扩束后光斑直径约100 mm，覆盖LWR入瞳。实验在暗室中进行，环境温度22±2°C。在视场中心（φ=0°，ε=0°）间隔30分钟进行10次重复测量，角度测量误差标准差为3.2 arcsec，表明系统稳定性良好。测试覆盖φ范围?60°至60°，ε范围0°至60°，共169个测试点，每点10次重复取平均。实测最大误差为4901.12 arcsec，与蒙特卡洛仿真预测值4716.46 arcsec相差约0.05°（184.66 arcsec），进一步验证了模型的有效性和可靠性。
总结讨论部分：研究人员讨论了各误差源的贡献和耦合效应，指出传统“最坏情况叠加”估计过于保守，而蒙特卡洛方法能够揭示多因素非线性耦合和部分抵消效应。实验系统表现出良好的重复性，随机测量噪声远小于系统误差，不影响验证结论。
研究结论可翻译为：本文针对LWR光学系统参数误差对角度测量精度的影响，建立了包含探测器焦平面倾斜和畸变等强非线性效应的完整前向映射链的综合误差传播模型，无需一阶近似，能够准确预测系统和随机误差，为后续标定和补偿提供了基础。系统地将蒙特卡洛方法应用于LWR光学系统多因素耦合误差分析，通过单因素和综合因素仿真获得了各误差源单独及组合作用下的角度测量误差概率分布。结果表明，主点偏移、畸变和焦距误差是总体不确定性的主要贡献项，而光斑中心提取误差影响可忽略。综合误差均值远低于各因素误差均值的线性叠加，表明误差源之间存在复杂的非线性耦合和部分抵消，凸显了蒙特卡洛方法相对传统最坏情况估计的优越性。通过120°×120°大视场下的最大误差法仿真和实验测量，验证了模型的有效性和准确性。实验最大误差（4901.12 arcsec）与蒙特卡洛预测（4716.46 arcsec）吻合良好，偏差仅约0.05°，证实了模型的可靠性。实验系统在视场中心的标准差为3.2 arcsec，表明随机测量噪声远小于系统误差。此外，本文建立的前向误差传播建模、蒙特卡洛随机采样和实验统计验证框架，不仅为LWR的误差标定、公差优化和高精度制造提供了工程指导，而且具有普适性，可推广至星敏感器、四象限探测器、视觉导航相机等其他成像式光学测量系统。
