非球面光学系统的光路计算方法与编程实现.pdfVIP

维普资讯 第21卷 第4期 青 岛大学学报 (工程技术版) Vo1．21No．4 2006年 12月 JOURNAL0FlQINGDAOUNIVERSITY (ET) DeC．2006 ● 文章编号：1006—9798(2006)04—0044—05 非球面光学系统的光路计算方法与编程实现 张风生 (青岛大学机电工程学院，山东青岛266071) 摘要：针对非球面光学系统的复杂性，给出了一种有效的光路计算方法，并导出了光路计 算的通用公式和符号规则。该方法首先以光线与过折射面顶点的切平面的交点为起点， 采取步进逼近的方法确定入射点在光线上的位置区间；然后利用二分法在该区间内高精 度求解入射点的坐标。可高精度求得入射点法线和折射光线方位，并计算出累积光程和 系统透过率。基于该方法编制的光路计算程序已应用于OCT(光学层析成像)系统的仿 真研究，其计算精度和效率得到了验证。 关键词：非球面镜 ；光路计算；OCT系统仿真 中图分类号 ：0435 文献标识码 ：A 光路计算亦称光线追迹，是基于几何光学的基本定律计算光线在光学系统中的传播路径。对于球面光 学系统，可利用球面方程直接求解空间光线在其上的入射点，然后求得折射光线。但由于非球面方程远比球 面复杂 ¨，无法通过直观的数学关系式求解入射点，只有采用逼近法进行求解。此时，如何既保证光路计算 的精度又加快算法的速度就非常重要。因为光路计算过程中的误差是非线性累计的，光线在前一折射面上 折射的微小偏差，将被后续折射面放大，导致光路计算结果误差很大，甚至完全错误的结果。而光路计算的 速度在某些应用中也很关键，例如光学层析成像 (opticalcoherencetomography，OCT)系统的仿真口]，需要 对数以百万计的光子通过光学系统的轨迹进行计算。本文采用步进逼近与二分法相结合，实现了光线在非 球面上入射点坐标的高精度快速求解，保证了光路计算的精度与效率。 1 光线在非球面上折射的计算方法 i．1 非球面 非球面通常以十阶或十二阶多项式表示 。设z为曲面的弯曲量，^=== ；为光线入射高度，则 z—f(x,y)一 7iCh2 十A^十Ae^十As^。十A。·^。 (1) 式中，C为曲面与光轴交点处的曲率；Az为圆锥系数，A ，Ae， 、、 h A。和A。为不同阶次项的系数 该非球面示意图见图1所示， 当CO时，曲面弯向正光轴方 向，即 0(图中实线)；当C 光线 ／ 0，曲面弯向负光轴方向，即zO(图中虚线)。当Az一0时，曲 面为球面；当A。一1时，曲面为双曲面；当Az=一1时，曲面 ／ 为抛物面；当一IA。O时，焦点在光轴上的椭 圆面上；当A O时，焦点在垂直于光轴的直线上的椭圆面上。 圈1 非球面示意圈 收稿 日期 ：2006-06—28，修 回日期 ：2006—10—3O 作者筒介 ：张风生 (1962一)，男，河北廊坊人，工学博士，教授 ，主要研究方 向为测试技术与智能仪器及光学层析成像技术 维普资讯 第4期 张风生：非球面光学系统的光路计算方法与编程实现 45 1．2 光线入射点的求解 图2所示为光线入射到第 i个折射面的示意图。入射光 线R与过折射面顶点的切平面交于P 点，在折射面上的入 射点为Pz。因为入射光线的方位和折射面在光路中的位置 是已知的，P 点的坐标 ( ， ，z)可直接计算，因