Optical Axis Static Sensitivity Analysis for Infrared Thermal Imaging Folding Optical System
-
摘要: 红外热成像折转光学系统在复杂环境条件下,光轴容易因为光学元件的偏心或倾斜而发生漂移,影响系统对目标的指示精度。在红外热成像系统设计之初对光学系统开展光轴静态敏感度分析,能够识别出光学系统的敏感点,为满足光轴稳定性的结构优化设计提供约束条件。通过基于旋转矩阵的坐标变换,建立了光学元件旋转过程量和倾斜状态量的转换关系,从而实现了光学元件在任意方向倾斜的空间姿态模拟,确保了光轴敏感度公差分析中的蒙特卡罗采样与结构设计的约束条件相对应,并在此基础上搭建了对红外折转光学系统光轴静态敏感度分析的流程,编制了程序。用所编程序对某典型红外热成像折转光学系统进行实例分析,根据光轴稳定性的指标要求,依次对光学系统中各光学件的偏心量和倾斜量进行了光轴的灵敏度和反灵敏度分析,得出了初始公差限,再针对初始公差限数据进行了任意方向采样的蒙特卡罗分析,最终得出了各光学元件能够满足光轴稳定性指标的偏心和倾斜公差限数据,通过建立多重坐标系的方法验证了所得数据的准确性,为指导光机热优化设计奠定了基础。Abstract: The optical axis of an infrared thermal imaging folding optical system is prone to shift owing to decentering of the tilt of optical components under complex environmental conditions, which affects the indication accuracy of the system for the target. Static sensitivity analysis of the optical axis for the optical system at the beginning of the design of the infrared thermal imaging system is useful for identifying the sensitive points of the optical system and provides constraints for the structural optimization design to meet the stability of the optical axis. The conversion relationship between the rotation process and spatial state quantities of the optical components was established by coordinate transformation based on the rotation matrix to simulate the spatial attitude of the optical component tilted in any direction and ensure that the Monte Carlo sampling in the optical axis sensitivity analysis corresponds to the constraint conditions of the structural design. On this basis, the flow of the static sensitivity analysis of the optical axis of the infrared folding optical system was established, and a program was compiled. A typical infrared thermal imaging folding optical system was analyzed using this program. According to the index requirements of the optical axis stability, the optical axis sensitivity and inverse sensitivity of the decenter and tilt of each optical component in the optical system were analyzed, and the initial tolerance limit was obtained. Then, Monte Carlo analysis sampling could be performed in any direction according to the initial tolerance limit data; thus, the decenter and tilt tolerance limit data that meet the optical axis stability index could be obtained, and the accuracy of the obtained data was verified by establishing a multi-coordinate system. Static sensitivity analysis provides a foundation for guiding the design of optical–mechanical thermal optimization.
-
0. 引言
变电站电力设备长期处于工况运行状态,易发生各类故障而威胁电网安全稳定运行。研究显示[1],在电气设备故障中超过半数均伴有设备异常发热现象。
20世纪中期国外的一些研究机构与电力企业尝试使用热像仪记录设备温度进行故障检测[2]。发展至今,国内外的技术人员通过对运行设备定期热拍照,由检测者观察设备热图像的温度分布并判定设备中存在的异常温升,确定热故障类型[3],该方法可有效检测出电力设备的异常温升。加上其无接触检测的优点,适应电力系统要求,因而在电力系统故障检测中得到推广运用。然而当前电力设备热图谱普遍采用人工观察,判断较主观[4-5],且其准确程度对检测人员经验依赖大,准确度波动大且不够高效,并且此类传统的人工图像特征分析方法较难提炼异常发热区域特征与分布。近年来计算机技术和人工智能技术飞速发展,研究者们提出了一些具有自动化和智能化功能的设备热像图故障识别方法,以提高判定与故障处理效率与准确性。对于自动化算法来说,如何利用电力设备的各类特点敏锐地发现热像图中的异常温升区域是热故障识别的关键。胡世征等人通过检测设备对应测量点间温差与其中较热点温升的比值,实现相对温差判断[6]。郭英军等提出使用红外摄像仪在线监测与计算变压器铁心的温升变化情况[7]。王如意等根据设备热图谱在异常温升处的色差进行图像分割而提取故障点[8]。孙怡等人提出利用不同电气设备出现异常温升的经验数值进行分类阈值判定的方法[9]。王淼等人通过分析导线热像图诊断输电线路[10]。冯振新等采用局部区域聚类的方法对热图谱温升区进行提取[11]。而现有研究较多关注于故障点程度的分割、计算和分析工作,而运用三相设备自有特征进行故障区域自动查找定位的研究却较少涉及。
在目前现状下,本文提出一种三相自搜寻比较的电力设备热像图自动诊断算法,变电站三相设备或并列运行的同类设备通常具有完全相同的外观特点,其电气运行状态相同或相近(严重三相不平衡和短路故障情况除外),因而正常运行时三相设备一般具有相同的热状态,对应部位的发热状态一致,这是热检测的重要判据。当其中某相设备发生局部热故障时,正常相设备可成为其余两相设备发热异常判定的良好参照物。通过计算机自动搜寻算法提取并利用三相对应的热状态信息对比,可自动高效地发现热故障并准确判定其发生部位,进行诊断分析。
1. 三相自搜寻比较算法
三相自搜寻比较法的思路是以热点区域为追踪目标,对三相设备热像图进行计算机处理,自动化地得到可直接相比较分析的三相图片。其方法实施流程如图 1所示。
为实现自动热区域搜寻和比较,本路线中先利用数学形态学的方法对图像中视图相连的三相设备进行分割处理,进一步将每相提取为单独图片,另外由于取景角度和杂物的影响,需对单相图像进行图像调整使其保持相同图像大小、角度,通过这样的自动化处理就得到了可进行对应区域映射对比的单相图片。接着为了最优程度地进行对应处比较,采用数据模块化方法,将处理完成的单相图像按照同等规格进行区块划分,由于先前的单相图同等化处理和同规格的分区,每相图片相同坐标的区块温度很大程度上反应了三相设备相同部位的温度水平。对每区温度取平均值,并对三相设备全区块进行温度分布统计,可很好地反应出各相的温度水平和温升情况,基于此相同的数据模型,可以便捷地实现设备温升判定和热数据信息分析,完成自动温升搜寻和故障比较识别。
2. 三相设备热像图的提取分割
2.1 热像图形态学处理
实现三相设备或同类设备之间的热像图对比,需要将图片中三相设备按相进行分离。由于图像中存在连接体、环境异物影响,以及拍摄角度不规范导致设备局部重叠等原因,致使照片中的设备间常有连接物存在,增加了三相设备主体形域的提取分割难度,如图 2所示。
图 2为电力断路器三相设备热像图及其提取的二值化图像。热像图中最亮处温度达68.1℃,最暗处为-0.2℃。从图中可看出,对图像二值化后,这些干扰一般体现为设备主体形状间的纤细连接,为去除这些连接干扰,可以采用对三相设备热像图的二值图像进行形态学开操作运算。
形态学是由动植物形态研究发展而来。数学形态学的处理可实现从图像中提取边界、骨架等表达描绘图形的图像分量,常用操作有腐蚀膨胀,开闭运算等。
形态学上称类似的图形间的纤细连接为不可靠连接,可以使用形态腐蚀膨胀等操作将其切断去除。其中,在二维图像Z2(即一个二维整数空间)以一个一定像素大小的结构元B对图形对象A进行空间卷积,$\hat B$为B的补集,B对A腐蚀定义为:
$$ A⊖B\{ z|{\left( B \right)_z} \subseteq A\} $$ (1) 而膨胀则是:
$$A \oplus B = \left\{ {z\left| {{{\left( {\hat B} \right)}_z} \cap A \ne \emptyset } \right.} \right\}$$ (2) 可得出,腐蚀操作可消除运算对象接触的所有边界点,使其边界向内收缩,该过程能够消除目标对象附近小而无意义的物体(可干扰体)[2],而膨胀操作则相反,其使目标区域范围“变大”。将目标区域接触的背景点合并到该目标物中,使目标边界向外部扩张[2],可以用来填补目标区域中某些空洞以及消除包含在目标区域中的小颗粒噪声。
在此基础上,利用腐蚀和膨胀操作可以衍生出开启运算:
$$ A \circ B = (A⊖B) \oplus B $$ (3) 即对图形对象进行先腐蚀后膨胀的操作。
以及闭合运算:
$$ A\cdot B=(A \oplus B)⊖B $$ (4) 即对图形对象进行先膨胀后腐蚀的操作。它们对图形操作具有不同的效果,都能部分平滑对象轮廓,开操作断开较窄的狭颈并消除细的突出物,闭操作则一般弥合较窄的间断和细长的沟壑,消除小的孔洞,填补轮廓线中的断裂等[12]。
结合三相设备热像图及其得到的二值图像的特点,可对其分别进行合适的开闭类操作,实现对不可靠连接干扰的处理。
如图 3所示,可看出开运算对短小纤细处等不可靠连接进行了“切断”处理,并去除了周边的小干扰体,而闭运算弥合了狭窄的目标间缝隙,使得目标扩大连接在一起。
针对去除不可靠连接,分离三相设备主体形态的目的,可使用合适程度的开运算操作,现采用合适大小的圆状结构元,对图 4所示三相断路器设备的热像图进行开运算操作。
可以看出,开运算操作对三相设备热像图二值图像中不可靠连接具有较好的“切断”作用。且开运算运用像素半径不同的结构元执行,将得到有区别的执行效果,结构元尺寸减小,图像将由从连接处“猛烈截断”到保留了更多区域细节,应根据图像情况调整合适的结构元大小。
2.2 三相设备主体的标定与提取
经过形态处理后还需进一步将三相设备主体提取出来,观察形态处理后的二值图像,可以看出白色目标区域主要是设备主体区域以及无用的杂散对象,在图像中,把这些具有相同像素值且相互邻接的一片像素区域作为连通区域或连通分量,在二值图像中即前景白色的数个区域。为进一步提取设备区域,可将图像中的连通区域分配不同数字标签并标定。根据其分布特性,将每个白色(像素值1)的连通分量像素标为同一值,不同连通分量分别标定1, 2, 3…等标签,形成新的标签矩阵,并将其用不同颜色表示,如图 5所示。
根据分析与标签结果,图像中标定的区域主要是设备主体区域以及附带设备、背景或未完全去除的线管等。依据开运算操作时处理尺度不同,标定区域的大小与数量也不同,而在所有的标签区域中,需要提取三相设备主体所在的连通分量。
根据设备热像图的拍摄特点分析,设备本体是拍摄的主要目标,对合格的设备热像图来说,电气设备主体应占据照片主要区域。利用这个特点,三相设备热像图中,标定连通区域中像素数量最大的3个区域即是三相设备的主体区域,即设备目标区域。
如图 6所示,搜集所有标定区域的像素数与标签号,找到具有前三多像素数的区域标签,并将其标记到新的矩阵图像中并进行颜色标定。可以看出,图像中标定了三相设备主体所在的最大的3个连通区域,去除了其他干扰区域的影响。
因二值图像目标区域像素值为1,而背景为0,进一步可将新的三相区域二值矩阵与原设备热像图执行对应元素相乘,如图 7所示,能够得到去除所有干扰与背景元素的三相设备热像对比图。
2.3 设备图像分相及调整
得到三相设备主体的热图像后,为实现每相设备之间的热状态对比,应将其每相分离至单张图片上,并裁切调整至一致。由于每相所在的连通区域已标定不同标签,能够确定标签及所在位置,因此可为每相区域创建相同大小图像,将三相区域分别标定,即能得到分离的三相设备热图像。
将每相设备图像如图 8单独分离后,由于其大小有差异以及背景位置不同,为进行对应位置的比对,应将其裁切并调整为相同图片大小与方向。利用每相二值图像记录的设备区域位置裁切,并可通过图像分辨率调整将其调整为等大图像,如图 9所示。
3. 三相设备区块对比与温升区域判定
采用同类比较法,判定设备某区域是否出现异常温升,在完成了设备提取和图像调整后,三相热像图应是大小一致的图像。热像图中由颜色亮暗代表温度高低,由比色条色彩标定,由此将热像图矩阵转变为对应像素的温度分布矩阵,即可进行温度对比。可将彩色热像图的RGB像素值转换为反应亮暗程度的灰度图像:
$$ {\rm{gray}} = 0.299R + 0.587G + 0.114B $$ (5) 由热像图比色条中给出的图像温度区间[Tmin, Tmax],将灰度图像素值矩阵映射至该温度区间,即:
$$ T = ({T_{{\rm{max}}}} - {T_{{\rm{min}}}}){\rm{gray}} + {T_{{\rm{min}}}} $$ (6) 由此可得相同大小设备温度分布矩阵。而根据电气设备的故障发热特点,某处的故障可能会造成该处一片区域发热而产生异常温升,因此,可将每相的温度分布矩阵划分为数个区块,每个区块温度取其均值,以反映每处的温度分布水平。如图 10是实验设备三相的分区块均温分布表。
根据各个区块的平均温度值分布情况,还可对各相不同温度范围内的区块数量进行统计,更直观地分析其温度分布情况,如图 11将示例三相断路器的每相的温度分布的区块数量以直方图统计。由图可看出,排除背景的低温区块,该设备温度区块主要分布在30℃~40℃温度范围内。而在更高的温度层级,B相、C相有少量区块的温度跨过了一般范围直接达到60℃~70℃温度范围,相对大多数区块所在温度范围属于较明显的温升,而A相则相对正常,无少量区块相较有明显的温度跃升。
除讨论不同温度范围的区块数量外,各区块相对位置的温度分布情况是判定温升区域的重点。为了能够清晰地看出每相设备的区块温度分布情况,使用常规的红外图像比色条对温度分布表中每区块温度进行映射,如图 12所示。
将所有区块按温度均值映射回红外图像标记色。如图 13可见,可清晰获取感知每区块的均值温度,而按红外图像温度颜色标记的区块图可清晰地体现区块的温度水平及每相对应区块的温差。为更清晰体现温度比对水平,可进一步将红外比色区块图制为柱状图,如图 14。
至此,每相设备分区块的温度对比即已清晰明了,只要根据设定的阈值标准,即可判定异常温升的存在,在热成像电力设备故障检测领域中,即可根据同类比较,如同组三相设备、同相设备之间及同类设备之间对应部位的温差进行比较分析,这类判定方法中的相对温差法是一种良好的根据周边温度差值进行温升判定的检测方法[6],利用相对温差δt即同类设备对应测量点间的温差与其中较高温点温升的比值:
$${\delta _{\rm{t}}} = \frac{{{\tau _1} - {\tau _2}}}{{{\tau _1}}} \times 100\% = \frac{{{T_1} - {T_2}}}{{{T_1} - {T_0}}} \times 100\% $$ (7) 式中:τ1和T1为发热点的温升和温度;τ2和T2为正常相对应点的温升和温度;T0为环境温度参照体的温度。以示例三相断路器图谱分析,设定异常温升参考值为20 K,进行区块温差比较计算发现,B相区块中部和C相顶部连接处出现异常温度高值,计算相对四周的δt过大、超限,结果符合图 15的温度区块数分布图的统计分析结果。由此能判定B相和C相均存在区域发生了超限温升热故障。
如图 15所示,准确地分析出设备存在的发热区域,是快速识别诊断热故障的关键。运用该方法可有效实现三相设备异常温升区域的自动查找。图 16和图 17分别为用该方法实现两组不同角度拍摄的三相电流互感器设备发热故障搜寻效果图。
基于此,三相自搜寻比较法按照图 1本文方法流程对电力设备热像图进行内部比较判断,自动分析识别,从而能更有效地识别设备温升与热故障。
该方法适用于具有成套相同设备、多相并列运行等易于进行相间温度比较的电力设备或者相内温升明显易确定相对温差的设备。若针对拍摄角度随意的多相设备热图谱,其记录的各相热信息缺失,较难准确地对各部位的温升情况做出故障判定,另外对于形态较为独特或不连续的设备,其判定能力可能受限。因此本方法对热像图拍摄质量有一定要求,后续研究可通过积累大量热故障判定实例来提高多相判定的准确性和对各种热图谱的兼容性,以提升方法智能化水平。
4. 结语
以电力设备红外热成像检测中同类比较法与相对温差法为基础,本文提出了一种基于电气设备热信息的三相自搜寻比较检测算法。对电网常见的三相设备、同类设备的检测热像图,本算法针对三相设备热像图的故障显示特征,基于图像处理和温差分析实施对热故障的自动搜寻与处理,获得三相温升比较对象,并进行分区块温度分布统计与自动温升比较分析,实现以同类比较法和相对温差法为分析基础的电力设备热故障高效自动化检测。与传统的设备热像图人工检测分析方法相比,本方法能够实现设备热故障的计算机自动搜寻检测,相对人工检测具有更高的效率和准确度,且能够实现三相设备热故障的全自动化搜寻检测。不足之处是对热像图拍摄角度、质量有一定要求,工程实际中需要对红外检测作业人员的拍摄角度进行必要规范,形成符合要求的三相设备热像图。
-
表 1 各光学件倾斜与偏心最终公差限
Table 1 Final tolerances of tilt and decenter of each optical component
LENS Tilt tolerance/° Decenter tolerance/mm Lens 1 0.17811389 0.01035527 Lens 2 0.19549745 0.00539823 Lens 3 0.15172749 0.00718520 Lens 4 1.50604475 0.01545493 Lens 5 1.55067910 0.02021238 Mirror 1 0.02377837 - Mirror 2 0.01895890 - Lens 6 0.23078320 0.00937822 Lens 7 0.11583835 0.00953447 Lens 8 1.05775137 0.01556741 -
[1] 郭夏锐, 廖志波, 王春雨, 等. 光轴一致性误差对空间透射式系统像差和质量的影响[J]. 红外与激光工程, 2012, 41(2): 437-441. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201202031.htm GUO Xiarui, LIAO Zhibo, WANG Chunyu, et al. Aberration and alignment quality investigation of space transmission optical system with centering error of optical axes[J]. Infrared and Laser Engineering, 2012, 41(2): 437-441. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201202031.htm
[2] 马宏川, 范宏波, 林宇, 等. 热像仪光机热集成分析综述[J]. 红外技术, 2019, 41(2): 134-141. http://hwjs.nvir.cn/article/id/hwjs201902005 MA Hongchuan, FAN Hongbo, LIN Yu, et al. Review of thermal-structural-optical integration analysis of thermal imager[J]. Infrared Technology, 2019, 41(2): 134-141. http://hwjs.nvir.cn/article/id/hwjs201902005
[3] 王阳, 黄煜, 李占峰, 等. 利用恒星对天文观测系统光轴平行性检校[J]. 红外与激光工程, 2017, 46(5): 0517003. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201705018.htm WANG Yang, HUANG Yu, LI Zhanfeng, et al. Calibration of optical axis parallelism by using star for astronomical observation system[J]. Infrared and Laser Engineering, 2017, 46(5): 0517003. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201705018.htm
[4] 陈志斌, 肖文健, 马东玺, 等. 大间距多光轴一致性野外在线检测方法[J]. 光学学报, 2017, 37(1): 0112006. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201701024.htm CHEN Zhibin, XIAO Wenjian, MA Dongxi, et al. A method for large distance multi-optical axis parallelism online detection[J]. Acta Optica Sinica, 2017, 37(1): 0112006. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB201701024.htm
[5] 杨耀山, 赵泉, 王建龙, 等. 红外光学系统光轴稳定性建模与分析[J]. 电光与控制, 2017, 24(12): 82-84. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ201712018.htm YANG Yaoshan, ZHAO Quan, WANG Jianlong, et al. Optical axis stability modeling and analysis of infrared optical system[J]. Electronics Optics & Control, 2017, 24(12): 82-84. https://www.cnki.com.cn/Article/CJFDTOTAL-DGKQ201712018.htm
[6] 隋杰, 程会艳, 余成武, 等. 星敏感器光轴热稳定性仿真分析方法[J]. 空间控制技术与应用, 2017, 43(4): 38-41. https://www.cnki.com.cn/Article/CJFDTOTAL-KJKZ201704006.htm SUI Jie, CHENG Huiyan, YU Chengwu, et al. A thermal stability analysis and simulation method for boresight axis of star sensor[J]. Aerospace Control and Application, 2017, 43(4): 38-41. https://www.cnki.com.cn/Article/CJFDTOTAL-KJKZ201704006.htm
[7] 董雪岩, 许方宇, 陈骥, 等. 1 m红外太阳望远镜镜面位置误差致光轴偏移分析[J]. 红外技术, 2016, 38(10): 870-876. http://hwjs.nvir.cn/article/id/hwjs201610011 DONG Xueyan, XU Fangyu, CHEN Ji, et al. Analysis on optic axis deviation induced by lens position disorder of 1 m new vacuum solar telescope[J]. Infrared Technology, 2016, 38(10): 870-876. http://hwjs.nvir.cn/article/id/hwjs201610011
[8] Kalikivayi V, Kumar V C P, Kannan K, et al. Tolerance analysis of misalignment in an optical system using Shack–Hartmann wavefront sensor: experimental study[J]. Optical Engineering, 2015, 54(7): 075104.
[9] 高磊. 基于蒙特卡洛法的公差分析及优化设计方法研究[D]. 哈尔滨: 哈尔滨理工大学, 2015. GAO Lei. Research on the Tolerance Analysis Based on Monte Carlo Method and Optimal Design Method[D]. Harbin: Harbin University of Science and Technology, 2015.
[10] 刘琳, 张兴德, 贺谊亮. 基于蒙特卡洛模拟法的红外光学系统公差分析[J]. 激光与红外, 2020, 40(5): 497-499. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201005012.htm LIU Lin, ZHANG Xingde, HE Yiliang. Monte Carlo simulation and its application in the IR optical system[J]. Laser & Infrared, 2020, 40(5): 497-499. https://www.cnki.com.cn/Article/CJFDTOTAL-JGHW201005012.htm
[11] 袁贵星, 王平. 蒙特卡洛模拟及其在公差设计中的应用[J]. 天津科技大学学报, 2008, 23(2): 61-64. https://www.cnki.com.cn/Article/CJFDTOTAL-TQYX200802019.htm YUAN Guixing, WANG Ping. Monte Carlo simulation and its application in tolerance design[J]. Journal of Tianjin University of Science & Technology, 2008, 23(2): 61-64. https://www.cnki.com.cn/Article/CJFDTOTAL-TQYX200802019.htm
[12] YANG K Y, JIN N, YANG D, et al. Monte Carlo calculation on stray radiation scattered by the baffle with a heterotypic surface for uncooled infrared system[C]// Proceedings of SPIE, 2018, 10815: 108150R.
[13] 苏朋. 365 nm光刻照明系统中变焦系统的设计及公差分析[J]. 红外与激光工程, 2022, 51(7): 20210524. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202207027.htm SU Peng. Design and tolerance analysis of the zoom system in 365 nm UV lithography illumination system[J]. Infrared and Laser Engineering, 2022, 51(7): 20210524. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202207027.htm
[14] 王川川, 岑兆丰, 李晓彤, 等. 基于线性模型的聚焦光学系统静态公差优化算法[J]. 光子学报, 2017, 46(11): 1122001. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201711025.htm WANG Chuanchuan, CEN Zhaofeng, LI Xiaotong, et al. Optimization algorithm of static tolerance in focusing system based on linear model[J]. Acta Photonica Sinica, 2017, 46(11): 1122001. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201711025.htm
[15] 蔡东浩, 李雅灿, 魏立冬, 等. 高光谱成像仪中曲面棱镜的装调公差研究[J]. 光学学报, 2021, 41(6): 60-68. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202106007.htm CAI Donghao, LI Yacan, WEI Lidong, et al. Fabricating tolerance of curved prisms in hyperspectral spectrometer[J]. Acta Optica Sinica, 2021, 41(6): 60-68. https://www.cnki.com.cn/Article/CJFDTOTAL-GXXB202106007.htm
[16] 李美萱, 董连和. 浸没式光刻照明系统中会聚镜的设计及公差分析[J]. 应用光学, 2019, 40(5): 876-881. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX201905027.htm LI Meixuan, DONG Lianhe. Design and tolerance analysis of converging lens in immersion Lithography illumination system[J]. Journal of Applied Optics, 2019, 40(5): 876-881. https://www.cnki.com.cn/Article/CJFDTOTAL-YYGX201905027.htm
[17] 曹宇泽, 马文礼. 两步式灵敏度矩阵法在卡塞格林望远镜装调中的应用[J]. 光电工程, 2020, 47(2): 5-12. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC202002001.htm CAO Yuze, MA Wenli. Application of two step sensitivity matrix method in Cassegrain telescope alignment[J]. Opto-Electronic Engineering, 2020, 47(2): 5-12. https://www.cnki.com.cn/Article/CJFDTOTAL-GDGC202002001.htm
[18] Synopsys Inc. Core Module User's Guide, Version 9.1[M]. Mountain View: Synopsys, 2020.
[19] ZEMAX LLC. Zemax Optic Studio 21.3 User Manual, Version September 2021[M]// Kirkland: ZEMAX, 2021.
[20] 王增伟, 赵知诚, 杨溢, 等. 基于刚体运动完备方程的光机热集成分析方法[J]. 红外与激光工程, 2022, 51(6): 20210617. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202206021.htm WANG Zengwei, ZHAO Zhicheng, YANG Yi, et al. Thermal -structural-optical integrated analysis method based on the complete equations of rigid body motion[J]. Infrared and Laser Engineering, 2022, 51(6): 20210617. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202206021.htm
[21] 张坤, 钟兴, 孟遥, 等. 大视场纳型星敏光学系统公差灵敏度研究[J]. 光子学报, 2019, 48(5): 0522001. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201905010.htm ZHANG Kun, ZHONG Xing, MENG Yao, et al. Tolerance sensitivity research of nano-star sensor optical system with large field[J]. Acta Photonica Sinica, 2019, 48(5): 0522001. https://www.cnki.com.cn/Article/CJFDTOTAL-GZXB201905010.htm