基于区域分解和逐级光线分裂的杂散光分析.pdfVIP

第29卷第6期 计 算 物 理 V01．29．NO．6 CHINESE OFCOMPUTATIONALPHYSICS Nov．．2012 2012年11月 JOURNAL 文章编号：1001—246XI 2012)06-0876-05 基于区域分解和逐级光线分裂的杂散光分析 陈 学， 孙 创， 夏新林 (哈尔滨工业大学能源科学与工程学院，黑龙江哈尔滨 150001) 摘要：针对蒙特卡罗法运算量大的问题，通过引入区域分解和逐级光线分裂技术，减少杂散光线与表面交点求解 时所涉及的表面数及跟踪光线数，降低计算量．区域分解时，考虑子区域与表面数的匹配关系；并基于杂散光传播 过程的模拟信息，逐级确定光线分裂倍数．以大衰减比抑制结构和多传感器系统为例，进行杂散光分析．结果表明， 对复杂系统杂散光传播的蒙特卡罗法模拟，应用区域分解和逐级光线分裂技术，可以在保持计算精度的前提下，有 效地提高计算效率． 关键词：杂散光；蒙特卡罗法；区域分解；逐级光线分裂 中图分类号：TKl24 文献标识码：A 0 引言 杂散光，也称杂散辐射，是光学系统中非正常传输的能量．在成像光学系统中，杂散光的产生导致系统信 噪比下降，影响成像质量，对杂散光的分析与抑制已成为光学工程关键技术之一¨。31．杂散光分析属于以小 量级光辐射能量传播及分布特性预测为目标的定向辐射传输计算，对复杂光学系统，广泛采用蒙特卡罗法进 行杂散光分析H。1． 由于蒙特卡罗法的计算量大，结合杂散光分析的特点，进行算法和求解技术改进研究一直受到重视，如 在杂散光分析商业软件中已应用的光线分裂和重点采样技术H1，杜胜华等M1对多光谱系统杂散辐射分析提 出的反向蒙特卡罗法和双向蒙特卡罗法等．这些研究，主要针对单传感器系统从系统入口面至成像面的杂散 光传输分析盯。8]，对开放式大型多传感器系统缺乏研究．因总体尺寸大、仪器间存在遮挡效应，开放式多传感 器系统的杂散光环境复杂，计算量大p“1．缺乏高效的计算方法，已成为此类系统杂散光分析与抑制设计的 关键制约因素． 根据杂散光传播的蒙特卡罗法模拟过程中，光线与表面交点求解是影响计算量的主要因素，本文引入区域分 解和逐级光线分裂技术，通过合理减少求解所涉及的表面数及跟踪光线数来有效降低计算量．以开放式多传感器 系统和大衰减比遮光罩抑制结构为例，采用该方法进行杂散光传播过程模拟，分析其可靠性和计算效率． 1 区域分解与逐级光线分裂方法 采用蒙特卡罗法分析杂散光时，抽样光线传播过程的模拟跟踪通过求解光线与系统表面的交点实现，是 主要的计算环节㈠“。2I．其交点求解，可表述为求解如下方程组…1 『f2厶+∞￡， ，1、 《 L、。1， 【西(f)=0， 式中，玉、f分别为光线源点和交点的位置矢量；∞为光线传播方向的单位矢量；奎(f)=0为系统表面方程； t为待求参数且满足t0． 由此可知，减少求解光线交点时涉及的表面数和跟踪光线数可减少方程组数量，降低计算量．传统的蒙 收稿日期：2012—03—05；修回日期：2012—05—27 基金项目：国家自然科学基金资助项目 作者简介：陈学(1985一)，男，博士研究生，主要从事红外仿真与杂散光分析研究 通讯作者：夏新林，E—mail：xiaxl@hit．edu．cn 万方数据 第6期 陈学等：基于区域分解和逐级光线分裂的杂散光分析 特卡罗法分析，以增大初始抽样光线数来保证结果的可靠性，导致计算量增加．如果根据系统整体布局采用 分界面(真实表面或者虚拟透明表面)将空间计算域划分为多个计算子区域，减少求解交点时涉及的表面 数；再将光线沿传播路径分解为多个能量级依次降低的区段(计算子区域或根据传播