简介一份聚焦M型Czerny-Turner光谱仪结构优化的技术资料面向光学系统设计与研发的科研人员和技术工程师也适合对光谱学原理有基础的进阶学习者。压缩包为单个docx文档大小约29KB内容紧凑既有理论推导也有可直接运行的Python代码和Zemax仿真验证流程。文档以反二次衍射条件为切入点用Python逐步演示光栅与准直镜、聚焦镜之间的几何关系求解并对牛顿迭代、scipy优化等数值方法做了说明光学仿真部分则详解了Zemax OpticStudio中的光路搭建、复色光源设置、光线追迹、点列图分析以及以Yspot为评价函数的优化策略。此外还给出了分辨率、球差等关键参数的计算函数便于读者量化评估成像质量。目前已有308人学习阅读后能快速复现光谱仪的优化计算逻辑形成从理论到仿真的完整方法链并可直接迁移至其他光学系统的性能分析中。1. M型Czerny-Turner光谱仪的结构优化为什么算法和仿真缺一不可M型Czerny-Turner光谱仪是当前便携式光栅光谱仪里最常见的紧凑布局入射狭缝、准直球面镜、平面衍射光栅、聚焦球面镜和线阵探测器光路在平面内折成M字形。好处是元件少、加工和装调相对宽容坏处是两块球面镜都以离轴方式工作像散和彗差无法从结构上消除只能靠参数配平来压低。结构优化做的事就是把镜距、两个反射镜的离轴角、光栅入射角当作变量先用Python写出像差模型和评价函数让优化器自动找到一组较优几何参数再把参数送进Zemax做光线追迹仿真确认点列图、能量集中度和波长标定都符合预期。这篇笔记适合准备做样机、或者刚拿到光栅不知道光路怎么排的工程师和学生目标是一步步带你跑通整条优化链路。2. 从光路布局到像差瓶颈M型Czerny-Turner的优化从拆解开始2.1 准直镜、光栅和聚焦镜的“M”字形光线到底怎么走经典Czerny-Turner结构中狭缝发出的发散光先被准直镜变成平行光投向平面光栅光栅按波长把光分散成不同方向的平行光束聚焦镜再把这些平行光束分别聚焦到探测器上。M型布局与经典结构的区别在于准直镜和聚焦镜放在光栅的同侧光从狭缝出发后先折向准直镜反射到光栅再从光栅折向聚焦镜最后落到探测器整体像字母M。这个折返排布让整机长度比单臂结构短差不多三分之一也正是微型光谱仪普遍采用它的原因。从优化角度看需要定的几何量包括准直镜焦距f1、聚焦镜焦距f2、两个镜面的离轴角θ1和θ2、光栅刻线密度、光栅入射角α以及狭缝到准直镜、光栅到聚焦镜的光程。狭缝到准直镜距离基本锁定在f1附近光栅到聚焦镜沿主光线的距离锁定在f2附近真正需要留给优化器的自由度主要是三个θ1、θ2以及f2与f1的比值。不要一开始就放开全部变量变量一多评价函数就会在局部极值里打转后面代码里你会看到我特意把变量控制在三四个以内。2.2 像散、彗差和场曲哪一项先把分辨率拖死先说最重要的像散。球面镜离轴使用时子午方向和弧矢方向的焦距不再相等。设球面镜焦距为f离轴角为θ子午焦距为f·cosθ弧矢焦距为f/cosθ两者之差近似为f·θ²。举个例子f100mm、θ5°时两个方向的焦点相差约0.38mm。对F数5左右的光束这个差值换算到探测器上的弥散斑会到几十微米直接超过常见线阵探测器10到25μm的像元尺寸。因此像散是M型Czerny-Turner里第一个要压下去的项。常见的做法是让准直镜和聚焦镜工作在相反的离轴斜方向上使两片镜子的像散部分抵消。彗差是第二个要处理的项。彗差与离轴角、光束口径都有关系准直镜把狭缝成像到无穷远离轴角大时会引入明显彗差聚焦镜工作在平行光入射状态彗差主要正比于离轴角θ2和口径。工程上常用Shafer条件确定初始点让两个镜子的焦距比和离轴角满足特定关系使两镜彗差在中间光栅处互相抵消。需要注意这只是近轴近似真实系统还要叠加光栅的像散最终要回到Zemax追迹结果确认。场曲在这个系统里通常不是第一优先级因为探测器是线阵或者短面阵宽度只有一两毫米。除非你的方案要求宽谱段平场否则先保证像散和彗差收敛场曲可以后面通过探测器倾斜或者加平场镜再补。2.3 初始结构参数比盲猜更可靠的起点下面是一组可见光波段常用的起步参数中心波长550nm工作带宽约200nm光栅刻线密度1200线/mm。组件参数起步值入射狭缝宽度25μm准直镜焦距 / 口径100mm / 30mm聚焦镜焦距 / 口径100mm / 30mm准直镜离轴角θ14.5°聚焦镜离轴角θ23.8°平面光栅刻线密度1200线/mm光栅入射角α23°探测器像元14μmθ1和θ2为什么不等因为光栅工作在非利特罗角度入射臂和出射臂的几何并不对称两个离轴角等大不代表像差平衡。起步值让θ1略大于θ2是给优化器一个明确的方向。可能你在不同资料里看到完全不同的起始角度这不奇怪因为M型结构对狭缝位置、光栅指向、探测器方向的定义方式不一样等效离轴角会有差异。只要初始结构能让光线顺利到达探测器就可以从这里开始优化不存在唯一正确答案。3. Python实现第一步光栅方程、像差模型和最小可运行代码3.1 坐标约定和光栅方程先把正负号统一坐标约定是这类项目里最容易翻车的点。我习惯把光栅面法线设为0°入射光和衍射光在法线同侧时角度取同号异侧时取异号光栅方程为mλ d(sin α sin β)其中d是光栅常数单位nmα是入射角β是衍射角m是衍射级次。为避免符号混乱下面的代码统一用角度制并在注释里写清楚“入射角相对光栅法线”。import numpy as np def grating_diffraction(wavelength_nm, lines_per_mm, alpha_deg, order1): 光栅方程: order * wavelength_nm d * (sin(alpha) sin(beta)) wavelength_nm : 波长, nm lines_per_mm : 刻线密度, lp/mm alpha_deg : 入射角, 度, 相对光栅法线 order : 衍射级次, 常用 1 或 -1 d_nm 1e6 / lines_per_mm sin_alpha np.sin(np.deg2rad(alpha_deg)) sin_beta order * wavelength_nm / d_nm - sin_alpha if abs(sin_beta) 1.0: raise ValueError(f波长 {wavelength_nm} nm 超出该级次自由光谱范围) return np.rad2deg(np.arcsin(sin_beta)) # 1200线/mm, 550nm, 入射角23° beta grating_diffraction(550.0, 1200.0, 23.0, order1) print(f一级衍射角: {beta:.3f}°)这段代码的输入是波长、刻线密度和入射角输出是衍射角同时会把无解的波长直接拦下来。做参数扫描时这个保护非常关键否则优化器会在无解区域里算出一堆NaN整个结果就废了。3.2 用Python估算离轴球面镜的像散和彗差后续优化需要一个能快速评估像质的评价函数。这里我不直接写完整光线追迹而是用离轴球面镜的子午和弧矢焦点公式估算像散再叠加一阶彗差项。它比Zemax粗糙但作为初筛完全够用。def mirror_astigmatism(focal_mm, off_axis_deg, na): 单个离轴球面镜的像散估算。 返回子午焦点、弧矢焦点和焦线弥散半径。 focal_mm : 近轴焦距, mm off_axis_deg : 离轴角, 度 na : 光束数值孔径近似值 theta np.deg2rad(off_axis_deg) ft focal_mm * np.cos(theta) fs focal_mm / np.cos(theta) # 在子午和弧矢焦点的中点放置探测器, # 弥散半径近似为焦距差的一半与NA的乘积 blur abs(fs - ft) / 2 * na return ft, fs, blur这里的物理依据很简单子午焦距变短弧矢焦距变长焦线间距就是二者差值。把探测器放在两者中点弥散斑半径约等于半焦距差乘以NA。这个公式只反映一阶像散没有计入球差和高阶彗差但对趋势判断足够了。3.3 把像散和彗差合成评价函数用差分进化做结构优化把准直镜和聚焦镜的像散、彗差合并成一个标量作为优化目标。变量设置成三个θ1、θ2和f2/f1。准直镜焦距f1固定为100mm光栅用1200线/mm中心波长550nm入射角23°。from scipy.optimize import differential_evolution def merit_function(params): 返回综合像差光斑半径, 单位mm. params [theta1_deg, theta2_deg, f2_ratio] theta1_deg, theta2_deg, f2_ratio params f1 100.0 f2 f1 * f2_ratio na 0.1 # 典型F/5光束对应的NA约0.1 # 两片镜子的像散残差 _, _, blur1 mirror_astigmatism(f1, theta1_deg, na) _, _, blur2 mirror_astigmatism(f2, theta2_deg, na) astigmatism_total abs(blur1 - blur2) # 一阶彗差项: 正比于离轴角, 方向和镜片配置有关, 这里取绝对值差 coma1 f1 * np.tan(np.deg2rad(theta1_deg)) * na ** 3 coma2 f2 * np.tan(np.deg2rad(theta2_deg)) * na ** 3 coma_total abs(coma1 - coma2) return float(np.hypot(astigmatism_total, coma_total)) bounds [(2.0, 8.0), (2.0, 8.0), (0.85, 1.15)] result differential_evolution(merit_function, bounds, seed42) print(最优结构参数:, result.x) print(最小综合像斑半径:, result.fun, mm)几个参数说明bounds里的角度范围限制在2°到8°过大离轴角会让彗差爆发过小则机械上放不下镜架和光栅。f2_ratio限制在0.85到1.15也就是聚焦镜焦距在准直镜焦距附近微调超出这个范围通常意味着光路几何已经不合理。differential_evolution是全局优化算法比L-BFGS这类局部算法更适合这个目标函数。当前函数虽然简单但真实使用中你很快会加入波长扫描项多峰特征会很明显全局搜索能避免困在局部极值。seed42用来固定随机数保证每次运行结果一致。如果你去掉seed两次优化结果可能因为初始种群不同而略有差异复现实验时要注意。跑完这段代码你可能得到一个类似[4.2, 3.6, 0.97]的结果。这意味着准直镜离轴角约4.2°聚焦镜离轴角约3.6°焦距比约0.97。它不是最终答案只告诉你应该往哪个方向去试。4. 把Python优化结果送进Zemax仿真搭建和验证4.1 在Zemax中搭建M型Czerny-Turner光路Zemax里搭离轴反射镜系统主流做法是“坐标断点加标准面”组合面型本身是标准球面镜镜子的位置和指向靠坐标断点控制。这样做的好处是后续优化可以直接把坐标断点的偏心或倾斜设成变量坏处是坐标断点顺序不能乱一旦回零不对光线就到处乱飞。我在序列模式里一般这样搭OBJ面物体类型选“物方数值孔径”物高取狭缝半宽约12.5μm。波长设置三个中心波长550nm再加450nm和650nm。STO面光阑放在准直镜位置。先插入一个坐标断点把坐标轴平移并倾斜到准直镜的入射方向再插入标准面曲率半径设为-200mm对应f1100mm口径30mm并勾选Mirror。之后再加一个坐标断点把坐标轴恢复到光栅入射方向。光栅面面型选“衍射光栅”把刻线密度填1200衍射级次填1。光栅面的刻线方向决定色散平面务必和坐标断点的旋转方向对应。聚焦镜再是一组坐标断点、标准面、坐标断点把聚焦镜放在衍射角对应的方向上曲率半径-200mmMirror。建议先把准直镜到光栅这一段搭通在光栅后面放一个临时探测器确认三个波长确实分开了再接聚焦镜和IMA面。分段调通比一次性搭完再排错快得多尤其是第一次做离轴系统时。4.2 Python变量和Zemax位置的对应关系Python变量Zemax中的位置准直镜离轴角θ1准直镜前坐标断点的Y倾斜聚焦镜离轴角θ2聚焦镜前坐标断点的Y倾斜准直镜焦距f1标准面曲率半径 -2×f1聚焦镜焦距f2聚焦镜标准面曲率半径 -2×f2光栅入射角α光栅面之前坐标断点的Y旋转还要注意符号约定。坐标断点里的“旋转”顺序会影响光线方向我只在YZ平面内做倾斜保持X方向为零。Python模型里算出来的角度都是正数到Zemax里可能某个角度需要取负具体看光路是从上方入射还是从下方入射。4.3 评价函数设置哪些变量放开哪些锁死Zemax优化向导里选“RMS光斑半径”模式视场取三个波长对应的小视场权重可以按中心波长1、边缘波长0.5来分配。变量设置建议如下可设变量聚焦镜前坐标断点的Y旋转角、准直镜前坐标断点的Y旋转角、光栅面之前的厚度间距。尽量锁死光栅刻线密度、衍射级次、镜面曲率半径。曲率半径也可以放开但通常焦距已经由机械结构决定放开后会得到非标准焦距反而难加工。探测器到聚焦镜的距离可以设成变量让Zemax自动找最佳像面位置。优化过程中如果出现光线追迹失败先把坐标断点的半口径加大到50mm等光路形状确认对再改回来。这个临时放大口径的做法几乎每个做离轴系统的人都会用算是通用的排错手段。4.4 Python和Zemax结果对不上先查这三个映射两套工具结果对不上不要急着怀疑代码先检查三处离轴角的定义不同。Python里θ1是准直镜入射角的绝对值Zemax坐标断点的倾斜值对应的是镜面法线的旋转反射后光线方向会以两倍角度改变符号很容易弄反。光栅衍射级次的正负。Zemax中级次的正负由坐标断点的旋转方向决定和Python里α的符号约定不一致时衍射角会整体反号。曲率半径正负。标准面设为Mirror时凹面镜的曲率半径在光线从左入射时通常取负值。设成正的话系统直接发散点列图会大得离谱。这三条逐项排查能解决大部分“仿真和计算对不上”的问题。5. 避坑与常见问题排查新手最容易翻车的五个点5.1 衍射角报“出界”优化器跑出NaN现象Python运行时报出自由光谱范围错误优化循环直接中断或者目标函数返回NaN优化结果毫无意义。原因光栅方程里mλ/d超过正弦函数的可解范围。通常是刻线密度、波长或入射角组合不当也可能是优化器跑到了边界外。解决在光栅衍射函数里加显式检查超出范围时返回一个很大的惩罚值而不是抛异常。这样差分进化能学会避开这个区域。刻线密度和中心波长确定后先手算一遍衍射角确认余量足够再开始扫描。5.2 Zemax追迹中途断线系统窗口一片空白现象优化一两步后光线在某个表面处突然消失点列图上没有能量。原因坐标断点的平移和倾斜顺序写错光线打到反射镜背面或者镜面半口径比实际光束尺寸小光线直接从旁边漏过去。解决先把所有反射镜半口径临时加到50mm确认整体光路形状正常。然后逐个检查坐标断点“倾斜、平移、标准面、平移、倾斜”的对称组是安全的基准写法。不要直接跳到非序列模式重做序列模式里养成分段调通的习惯排错更快。5.3 优化不收敛每次跑的结果都不一样现象用差分进化优化时目标函数不下降或者不同seed跑出完全不同的参数。原因目标函数里存在NaN或无穷大优化器在无效区域里反复试探另一个常见原因是边界设太宽比如θ1从0°到30°大多数区域内彗差都大得离谱优化器找不到有效梯度。解决确认所有分支都有出口无效区域返回一个大常数比如1e6。把不等式约束写进函数第一行。θ1和θ2的边界收窄到你愿意接受的结构范围即可。另外固定seed复现实验时记录下来这不是玄学是工程可溯源性的基本要求。5.4 一级和二级衍射光谱混叠现象仿真中点列图没有问题但实际装出样机后某些波长在探测器上出现两个像斑标定不对。原因一级衍射和二级衍射光谱在宽谱段下交叉重叠。只要光源覆盖超过一个倍频程就一定存在这个串扰。解决优化时把带宽作为输入约束检查整个范围内的自由光谱范围如果带宽超过一个倍频程加入透射式长通滤光片或者换更高刻线密度的光栅并重新算入射角。这个坑在Python模型里往往不会被注意到因为单波长优化看不到重叠。5.5 Python最优参数和Zemax最优参数偏差太大现象Python优化出θ14.2°、θ23.6°送进Zemax重新优化后变成θ15.1°、θ23.2°偏差超过10%。原因Python的简化模型漏掉了光栅在非利特罗角下的固有像散也忽略了镜面球差的非线性部分。一阶公式只能反映趋势不能直接作为最终装调参数。解决把Python结果当作初始点在Zemax中为每个变量设定±15%的范围再做一轮局部优化。最终以Zemax点列图为准。两套工具的正确关系是Python负责快速全局搜索Zemax负责精确落地缺一不可。6. 把优化过程自动化Python驱动Zemax做闭环迭代6.1 用zospy把参数循环写进一个脚本手动跑完一轮“Python初筛、Zemax验证、人工调整参数”之后你会发现真正的项目需要重复这个流程几十次。波长要覆盖多个目标值温度公差要扫一遍装调容差也要看。手工做下去效率实在太低常见做法是用zospy操作Zemax的接口直接改坐标断点的值、跑评价、读取RMS半径。import numpy as np import zospy as zp # 连接Zemax并打开已搭好的M型Czerny-Turner文件 oss zp.OpticStudioExtension() oss.load_file(CT_M_initial.zmx) # 假设第7面是聚焦镜前坐标断点, 把Y倾斜设为3.5度 oss.system_editor.set_cell(7, tilt_y, 3.5) # 运行点列图分析读取RMS半径 spot oss.analysis.spotdiagram() rms spot.data.rms_radius print(f当前结构RMS光斑半径: {rms:.3f} um)这段代码的价值不在单次修改而在于能用for循环把theta1从3°扫到6°。每次修改后读取对应的RMS值最后用numpy的argmin找到最佳结构。zospy的调用方式随版本会有差异但思路就是这个思路。跑自动化之前先手动确认一次单次修改能正常写入Zemax再去跑批量循环避免一次改几个变量导致系统崩掉还找不到原因。6.2 我会怎么验证优化结果而不是只看点列图除了RMS光斑半径我还会加两个验证指标这两个指标能避免后期装调时被公差打脸。第一个是像面离焦容差。把探测器位置前后各偏移50μm看RMS半径是否明显恶化。偏移后如果RMS翻倍说明系统对焦深太敏感这种结构生产一致性会很差。第二个是波长标定线性度。对工作带宽里的三个以上波长分别用光栅方程计算探测器上的理论位置再在Zemax里读取实际主光线位置两者的偏差应小于一个探测器像元。这样能确认光栅安装角度的误差不会让波长标定直接报废。把这些验证脚本嵌进zospy循环你会得到一张“角度-离焦-标定偏差”的表格。我在最初几次做这类优化时总以为RMS够小就算完事结果装样机后被离焦容差坑过好几回。后来养成一个习惯每次跑完仿真先手动在纸上写一遍光栅方程验算入射角确认数值一致再继续下一步。看似多花五分钟却能让后面调试少走很多弯路。希望这个习惯对你也管用。本文还有配套的精品资源点击获取
