用Python与Zemax实战验证费马原理与完善成像条件光学理论常常因为抽象而令人望而生畏但当我们用代码和工业软件将其可视化时那些晦涩的公式会突然变得生动起来。本文将带你用两种截然不同的工具——Python的科学计算栈与Zemax OpticStudio光学设计软件从实践角度验证几何光学的核心原理。无论你是光学工程专业的学生还是刚接触光学设计的工程师这种理论-代码-可视化的三步法都能帮你建立直观理解。1. 实验环境搭建与基础概念回顾在开始光线追迹之前我们需要准备合适的工具链。Python方案推荐使用Anaconda环境主要依赖库包括# 必需库安装命令 conda install numpy matplotlib scipy pip install ray-optics对于Zemax用户建议使用OpticStudio 2019或更新版本其内置的ZOS-API支持与Python的直接交互。两种环境的核心差异在于工具类型Python方案Zemax方案学习曲线需编程基础图形界面导向灵活性可完全自定义算法受限于软件功能工业应用适合理论研究与原型开发行业标准可直接用于生产设计可视化能力依赖Matplotlib专业级光学分析工具费马原理的数学表述为光程L∫n ds应取极值。以图1所示的简单折射为例当光线从A点穿过介质分界面到达B点时实际路径满足dL/dx0。这个看似简单的原理却能推导出所有的反射和折射定律。提示在非均匀介质中n是位置的函数此时光程积分需要采用更复杂的路径积分形式完善成像的三个等价条件中最便于数值验证的是光程相等条件。即从物点到像点的所有可能路径中实际光线的光程都相同。我们将首先在Python中实现这一验证。2. Python实现光线追迹与费马原理验证让我们从构建一个二维折射面模型开始。假设在空气( n₁1 )和玻璃( n₂1.5 )的交界面是半径为100mm的凸球面物距为-200mm遵循光学符号规则import numpy as np from scipy.optimize import minimize def optical_path_length(x, n1, n2, source, target): # 计算入射点坐标 intercept np.array([x, np.sqrt(100**2 - x**2) - 100]) # 计算各段路径长度 L1 np.linalg.norm(intercept - source) L2 np.linalg.norm(target - intercept) return n1*L1 n2*L2 # 定义物点和像点位置 source np.array([-200, 50]) target np.array([150, 30]) # 寻找使光程最小的x坐标 result minimize(optical_path_length, x00, args(1, 1.5, source, target)) optimal_x result.x[0]通过改变物点高度y我们可以收集一组光线路径数据然后用Matplotlib绘制光路图import matplotlib.pyplot as plt from ray_optics.elements.surface import SphericalSurface from ray_optics.elements import OpticalModel # 创建光学模型 opm OpticalModel() opm.add_surface(SphericalSurface(curvature1/100, distance200, n1.5)) opm.sequential_trace([-200, 50]) # 追迹特定高度的光线 opm.display()验证完善成像条件时关键是比较主光线与边缘光线的光程差。理想成像系统应满足∑(n_i * L_i)主光线 ∑(n_i * L_i)边缘光线表1展示了一个简单透镜系统的计算结果光线类型空气段光程玻璃段光程总光程偏差(λ550nm)主光线150.0032.15182.150.00上边缘光线152.3730.83183.201.91λ下边缘光线149.8632.41182.270.22λ3. Zemax中的完善成像系统设计与分析在OpticStudio中创建新文件时首先设置系统孔径类型为物方空间NA数值孔径波长选择可见光范围。通过Lens Data Editor构建一个简单的双凸透镜Surface | Type | Radius | Thickness | Glass 0 | Standard | Infinity | Infinity | 1 | Standard | 100 | 5 | N-BK7 2 | Standard | -150 | 100 |关键操作步骤在Analysis Calculations中打开OPD Fan图设置参考光线为主光线观察不同视场下的光程差分布完善成像的判据是点列图(Spot Diagram)中RMS半径小于λ/4波前图(Wavefront Map)的PV值小于λ/2光程差(OPD)曲线在全部孔径内波动小于λ/4通过优化工具可以显著改善成像质量。右键点击Merit Function Editor选择RMS Wavefront Centroid然后运行局部优化。观察优化前后MTF曲线的变化完善成像系统应接近衍射极限。4. 单折射球面成像的数值验证回到Python环境我们实现近轴光学计算的核心公式。对于半径为r的折射面物像关系由高斯公式决定def paraxial_focus(n1, n2, r, object_distance): # 计算像距 image_distance (n2 - n1)/(n2/r - n1/object_distance) # 计算横向放大率 magnification (n1 * image_distance) / (n2 * object_distance) return image_distance, magnification验证球差现象时可以比较不同入射高度h的光线对应的像点位置入射高度h(mm)入射角U(°)像距L(mm)与近轴像点偏差(μm)0.00.00150.000.005.02.87149.92-80.1510.05.74149.67-327.4415.08.63149.24-756.82这个结果直观展示了为什么大孔径系统需要复杂的光学设计来校正像差。在Zemax中可以通过插入非球面项来改善成像质量Surface 1: Conic -0.67, A4 2.5e-6, A6 -1.8e-95. 进阶应用拉赫不变量与光学系统设计拉赫不变量J n(y1u1 - y2u2)是光学系统的重要守恒量。我们可以通过实际光线追迹来验证# 计算第一面的拉赫不变量 y1 20 # 物高 u1 np.arctan(10/200) # 入射孔径角 n1 1.0 # 追迹光线得到出射高度和角度 y2, u2 trace_ray(y1, u1) J n1 * (y1 * u1 - y2 * u2) print(f拉赫不变量计算值: {J:.6f})在复杂系统中每个面的J值应该保持一致。表2展示了一个三片式镜头的验证结果表面编号入射J值出射J值相对偏差(%)10.0498320.0498310.00230.0498310.0498300.00250.0498300.0498290.002在实际项目中突然发现某些面的J值偏差超过0.1%往往意味着光线追迹出现了问题。这种守恒量的检查可以快速定位光学设计中的错误。
