1 引言
在消色差系统中,残余色差以二级光谱的形式存在,当对成像质量要求更高时,需要对二级光谱进行校正[1-2]。已校正二级光谱的系统称为复消色差光学系统。在需要更进一步提高光学系统的分辨率时,需要对复消色差光学系统中的三级光谱进行校正。
校正色差可利用基于几何图形和数学方法求解出的最佳玻璃组合[3-7]。这些方法有助于由两种或三种玻璃组成的简单系统的二级光谱的校正。但是,当光学系统变得更复杂时,波段加宽,所使用的光学玻璃种类增多,这些方法不再适用。此外,还有研究者提出使用基因优化算法[8]来优化玻璃选择。尽管这种方法可以得到非常好的结果,但计算量大,即使对于简单的光学系统也需要对大量的数据进行处理。
为了能够快速选择复消色差,得到可对三级光谱进行校正的最佳光学玻璃组合,本文提出了一种子波段缝合的方法,对设计波长进行拆分和缝合,并针对拆分和缝合后产生的子波段进行复消色差以选择最佳的光学玻璃组合。
2 光学玻璃选择方法
单透镜与消色差系统的初级轴向色差曲线与二级光谱曲线分别如图1(a)、(b)所示,其中F为焦点。可以看出,消色差光学系统的轴向色差曲线呈U形,并具有一个拐点,可以使消色差波长对应相同的焦面位置。图1(c)、(d)所示为复消色差透镜系统的典型轴向色差曲线,可以看出,复消色差光学系统的轴向色差曲线有两种形式:一种为S形,该曲线具有两个拐点,可以使得两个边缘波长与中心波长对应同样的焦平面位置,两边缘波长与中心波长之间的波长对应的焦面位置相对复消色差焦平面位置的距离即为三级光谱;另一种为W形,该曲线具有三个拐点,故其对应的三级光谱残余量要小于S形的[9]。
校正三级光谱就是在指定的光谱范围内对4个波长校正轴向色差,使得4个波长对应的焦面位置重合[10]。对于图1(d)所示的这种情况,将焦平面向左移动一点就可以产生4个焦点,因此能够产生W形轴向色差曲线的复消色差光学系统,其三级光谱已经得到了校正。
图 1. 光谱。(a)初级色差;(b)二级;(c)三级,S形;(d) 三级,W形
Fig. 1. Spectra. (a) Primary; (b) secondary; (c) tertiary, S-shaped; (d) tertiary, W-shaped
下载图片 查看所有图片
为了能够产生W形的轴向色差曲线,除对边缘波长和中心内波长进行复消色差外,还要考虑其余波长。为能引入其余波长,采用子波段缝合型复消色差校正方法。
2.1 子波段缝合光学玻璃组合的选择方法
由j种玻璃构成的薄透镜模型如图2所示,为了使该系统能够对某一波段实现复消色差[11],有
式中ϕi为第i个透镜的光焦度;Pi为第i个透镜所使用的材料的相对部分色散系数;νi为第i个透镜所使用材料的色散系数。为了能够实现对多色焦点偏移曲线的控制,需要对整个设计波段进行拆分再缝合,过程如图3所示。设这个波段中含有m个特征波长(λ1,λ2,…,λm且λ1<λ2<…<λm),相邻的三个特征波长构成一个子波段区间,并将其作为缝合的基本元素。一级缝合是在基础子波段λi-1、λi、λi+1上引入一个相邻的特征波长λi+2,将新引入的波长作为边缘波长,以原有的中心波长λi和边缘波长λi+1分别作为新生成的范围更大的子波段区间各自的中心波长,生成两个新的子波段区间作为一级缝合的结果,即λi-1、λi、λi+2与λi-1、λi+1、λi+2。例如,以λ1、λ2、λ3所构成的子波段作为基础子波段,一级缝合在其基础上加入λ4,一级缝合产生的结果是以λ1、λ4作为边缘波长,以λ2、λ3作为中心波长的两个子波段区间,即λ1、λ2、λ4与λ1、λ3、λ4。二级缝合则是在一级缝合的基础之上进一步引入一个相邻的特征波长,并将其作为新的子波段区间的边缘波长,再重复缝合的过程。以此类推直至包含m个特征波长的m级缝合,覆盖整个设计波段。
图 2. 薄透镜系统示意图
Fig. 2. Schematic of thin lens system
下载图片 查看所有图片
图 3. 波段的拆分与缝合
Fig. 3. Splitting and stitching of frequency bands
下载图片 查看所有图片
对全部的光学玻璃重新计算相对于每个基础子波段区间以及缝合产生的新子波段区的相对部分色散P和色散系数ν,同时计算
与
的值,分别记为V与ξ。则由(2)、(3)式可以得到
式中各参数
=,
=,
=,…,
=,
=,
=,
=,…,
=;V与ξ的下角标代表计算波段,上角标代表不同的玻璃。
将(4)、(5)式中由V与ξ构成的矩阵分别记为
与
,故(4)、(5)式可合写为
式中
=[
]T;T表示矩阵的转置,下同;0为零矩阵。
故(1)式可写为
式中
=[
]。
将(6)、(7)式联立可以得到
式中
=[
]T。(8)式即为j块光学玻璃对m(m>j)个波长的复消色差条件。
对所设计波段的拆分和缝合使得方程的数量多于未知数的数量,故(8)式无解。在这种情况下使用最小二乘解作为方程的近似解。令矩阵参量
则有
由(10)式可以求解出j种薄透镜的光焦度近似解,这种近似解并不能严格满足m个主要波长上的复消色差条件,即不能使(8)式两边严格相等。此时计算各个子波段的色差残余值,将
代入(6)式中,令矩阵
则方差值D=
。D表示各个子波长的焦距对于各自中心波长焦距的偏离程度,该值可以用来评估该组光学玻璃的色差校正能力。如果D小于所设计的光学系统的焦深,则可认为该组光学玻璃可消除该系统内产生的色差。
除了对色差的校正能力的考量,还希望玻璃所承担的光焦度在一个合理的范围内。单个透镜如果承担很大光焦度,其表面曲率会很大,会给单色像差的校正和加工带来诸多问题。同时希望该组光学玻璃可产生较小的高级球差和色球差。
2.2 边界条件的计算
对于已知光焦度的一个简单薄透镜,其球差[12]可以表示为
式中y为半口孔径;ϕ为光焦度;β为角弥散,其表达式为
式中n为折射率;K用来表示透镜的形状,其表达式为
式中r1和r2分别为透镜的前、后表面半径。
由(12)~(14)式可知,单薄透镜的球差是关于折射率和形状的函数。
单薄透镜的球差最小形状为
此时,球差为只与折射率有关的函数:
此形状K下的色球差为
式中长波时的球差LC=-2βC/(y);短波时的球差LF=-2βF/(y),其中βC、βF分别为长波和短波的角弥散,ϕC、ϕF分别为长波和短波的光焦度。
根据高级球差估计式可知,此时的高级球差[13]为
式中L0.707为0.707孔径处的球差;Lm为边缘孔径处的球差。
通过计算可以获得j个薄透镜组中每个透镜的光焦度
、多色焦点偏离程度D、高级球差ΔL'Sn以及色球差ΔLFC。合适组合中单个透镜的光焦度不应过大,残余的色差值应小于系统焦深,高级球差和色球差应尽量小。
2.3 玻璃组合的选择步骤
要明确光学系统的使用波段和有效焦距(EFL)等设计参数,再根据经济性等因素对玻璃库进行初步的筛选,剔除不常见玻璃或成本较高的玻璃。将筛选后的玻璃库中的玻璃按照每一种组合中含有j种光学玻璃进行组合,并将所有可能的组合进行编号。对设计波段进行拆分和缝合,对每一个子波段计算每一种光学玻璃的P、ν值,利用(11)式得到光焦度分配
,令
=,
表示各个元件的光焦度绝对值之和。若
值过大,则认为该组合的光焦度分配不合理。对光焦度分配合理的组合计算D值,将D值乘以有效焦距(EFL)判断该乘积是否小于光学系统焦深,若是则将结果保留。对所有可能的光学玻璃组合进行计算,再对满足
、D值要求的组合计算ΔL'Sn与ΔLFC。选择ΔL'Sn与ΔLFC大小合适的组合进行光学系统设计,该过程如图4所示,其中N为循环次数, F number为相对孔径倒数。
图 4. 玻璃组合选择方法的流程图
Fig. 4. Flow chart of glass combination selection methods
下载图片 查看所有图片
3 设计实例
根据设计流程完成一个宽波段长焦相机镜头的优化设计,在设计过程中只使用球面镜。选择的设计实例参数见表1。
表 1. 光学系统要求
Table 1. Specifications for optical system
Parameter | Requirement |
---|
EFL /mm | 400 | Field of view /(°) | ±5 | F number | 6 | Spectral waveband /μm | 0.45-1.014 | MTF | Close to diffraction limit |
|
查看所有表
从给定波段中选择5个特征波长:0.45, 0.5461, 0.6438,0.852,1.014 μm,将0.6438 μm作为EFL的中心波长。计算所采用的光学玻璃库为从成都光明(CDGM)玻璃库中选取的综合性能优且常年在库的光学玻璃,计算光学玻璃在所有的子波段区间上的V与ξ。将材料库中的材料每三个一组进行排列组合,利用(10)式计算每一组光学玻璃中每一种材料的归一化光焦度Φ1、Φ2、Φ3。计算归一化光焦度绝对值之和
,通过
值的大小可以判断该组透镜的光焦度是否合理,记录
值小于10的玻璃组合,将所使用的光学玻璃牌号写入表2中。当测试完所有的玻璃组合,对满足
值小于10条件的玻璃组合按照(11)式计算残差值,再由残差值求出方差值D。由于方差D值是按照归一化的光焦度进行计算的,在与光学系统的焦深进行比较时需要乘以EFL。D值乘以光学系统的EFL得出的值为d。再利用(17)、(18)式计算ΔLFC与ΔL'Sn的值,结果见表2。
表 2. 光学玻璃组合的参数
Table 2. Optical glass combination parameters
N | Mark 1 | Mark 2 | Mark 3 | Φ1 | Φ2 | Φ3 | | d /(10-2 mm) | ΔL'/10-4 | Δ /10-4 |
---|
70 | H-FK61 | H-BAK7 | H-ZF62 | 3.70 | -2.79 | 0.1 | 6.59 | 3.07 | 2.33 | 5.89 | 74 | H-FK61 | H-BAK7 | H-ZF88 | 3.70 | -2.77 | 0.08 | 6.55 | 4.08 | 2.34 | 6.81 | 34 | H-FK61 | H-ZF62 | H-BAK8 | 4.23 | 0.14 | -3.38 | 7.75 | 4.40 | 2.49 | -12.5 | 45 | H-FK61 | H-ZF52 | H-BAK8 | 4.27 | 0.18 | -3.46 | 7.91 | 4.44 | 2.49 | -17.8 | 119 | H-FK61 | H-ZF88 | H-BAK8 | 4.22 | 0.11 | -3.34 | 7.67 | 6.10 | 2.50 | 20.0 | ︙ | ︙ | ︙ | ︙ | ︙ | ︙ | ︙ | ︙ | ︙ | ︙ | ︙ |
|
查看所有表
将表2中的数据按照d值由小到大排列。根据设计要求,计算出该光学系统的焦深为±0.046 mm,前4组的d值都小于焦深,再比较ΔL'Sn与ΔLFC。第1组的ΔL'Sn与ΔLFC明显比其他几组的都小,故使用第1组光学玻璃组合来完成设计。
光学系统图如图5所示,采用7片球面镜片、3种玻璃,其中第1、2、7片使用了牌号为H-FK61的光学玻璃,第3、4、6片使用了牌号为H-BAK7的光学玻璃,第5片使用了牌号为H-ZF62的光学玻璃。
图 5. 光学系统结构示意图
Fig. 5. Structural diagram of optical system
下载图片 查看所有图片
光学系统的成像性能如图6所示,其中RMS代表均方根,T表示子午方向,R表示弧矢方向。由图6(a)知,各个视场中的弥散斑尺寸与艾里斑大小十分接近,这与图6(c)中MTF接近衍射极限相对应,说明该光学系统的成像质量良好。由图6(b)可知,该系统的最大畸变小于0.25%。
图 6. 光学系统的成像性能。(a)点列图;(b)畸变;(c) MTF
Fig. 6. Imaging performances of optical system. (a) Spot diagram; (b) distortion; (c) MTF
下载图片 查看所有图片
图7所示为最佳像点位置随着波长变化的曲线图,最佳像点的最大偏移量为-0.03 mm,小于光学系统本身的焦深±0.046 mm。该曲线图呈W形,与零点出现了四处交点,说明该系统使得四种波长的最佳像点位置重合在一起,实现了三级光谱的校正。
图 7. 多色焦点漂移曲线
Fig. 7. Chromatic focal shift
下载图片 查看所有图片
4 结论
提出了一种针对宽波段长焦距光学系统的光学玻璃选择的优化计算方法。该方法可以快速得出符合三极光谱校正的最优光学玻璃组合,并给出了设计实例。研究结果表明,通过所提计算方法得到的光学玻璃组合实现了0.45~1.014 μm波段的复消色差,并且能够实现W形的多色焦点偏移曲线,即对三级光谱进行了校正,且各个视场的成像质量良好,MTF接近衍射极限。
王靖锋, 崔庆丰. 宽波段光学系统三级光谱的校正方法[J]. 光学学报, 2018, 38(11): 1122002. Jingfeng Wang, Qingfeng Cui. Calibration Method of Tertiary Spectrum for Broadband Optical Systems[J]. Acta Optica Sinica, 2018, 38(11): 1122002.