CN104036071B - 大动态范围光学系统计算机辅助装调方法 - Google Patents
大动态范围光学系统计算机辅助装调方法 Download PDFInfo
- Publication number
- CN104036071B CN104036071B CN201410216269.3A CN201410216269A CN104036071B CN 104036071 B CN104036071 B CN 104036071B CN 201410216269 A CN201410216269 A CN 201410216269A CN 104036071 B CN104036071 B CN 104036071B
- Authority
- CN
- China
- Prior art keywords
- optical system
- adjustable optical
- free degree
- aberration
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
本发明涉及一种大动态范围光学系统计算机辅助装调方法,该方法包括以下步骤:1)计算并建立待装调光学系统随机失调点的灵敏度矩阵数据库;2)根据步骤1)中所建立的待装调光学系统随机失调点的灵敏度矩阵数据库以及待装调光学系统不同视场的像差测量值获取待装调光学系统各自由度调整量的计算结果;3)根据步骤2)中的装调光学系统各自由度调整量的计算结果对装调光学系统的各自由度进行调整。本发明提出了一种可实现大动态范围光学系统计算机辅助装调、能降低对失调光学系统粗装调的要求的大动态范围光学系统计算机辅助装调方法。
Description
技术领域
本发明属于光学领域,涉及一种大动态范围光学系统计算机辅助装调方法,尤其涉及一种基于蒙特-卡罗法和动态数据交互的大动态范围光学系统计算机辅助装调方法。
背景技术
计算机辅助装调就是通过对待装调光学系统的波像差进行实时检测,并通过计算方法对得到的波像差进行处理和优化,给出系统各个光学元件明确的调整方向和具体的量值,从而通过各种定位工具和调整工具使光学元件回到其实际的最佳位置,使光学系统成像质量以及各元件的相互位置可以接近于理论设计指标。
目前计算机辅助装调过程主要采用灵敏度矩阵法,其近似认为待装调光学系统像差与调节量之间为线性关系,即它们之间可以表示为下式:
A×ΔX+H=ΔF (1)
式中,A为灵敏度矩阵,ΔX为待装调光学系统主镜和次镜的自由度调整量,ΔF为系统不同视场像差偏差,H为噪声,通过对(1)式求解可得到系统调整量。该方法具有以下缺点:
1)由于灵敏度矩阵A为奇异矩阵,(1)式没有精确解,利用传统的最小二乘法,求解的结果对灵敏度矩阵性质的要求比较高,且失调光学系统各光学元件自由度调整量之间的近似相关性,都会导致解的收敛性无法满足。
2)常规的灵敏度矩阵法是在线性近似的前提下进行计算,其在以系统理论基准点(各视场波像差最小)为中心的很小范围内,(1)式求解的结果才能保持很好的收敛性,故其动态范围很小。
发明内容
为了解决背景技术中存在的上述技术问题,本发明提出了一种可实现大动态范围光学系统计算机辅助装调、能降低对失调光学系统粗装调的要求的大动态范围光学系统计算机辅助装调方法。
本发明采用的技术解决方案是:本发明提供了一种大动态范围光学系统计算机辅助装调方法,其特殊之处在于:所述大动态范围光学系统计算机辅助装调方法包括以下步骤:
1)计算并建立待装调光学系统随机失调点的灵敏度矩阵数据库;
2)根据步骤1)中所建立的待装调光学系统随机失调点的灵敏度矩阵数据库以及待装调光学系统不同视场的像差测量值获取待装调光学系统各自由度调整量的计算结果;
3)根据步骤2)中的装调光学系统各自由度调整量的计算结果对装调光学系统的各自由度进行调整。
上述步骤1)的具体实现方式是:
1.1)通过蒙特-卡罗法生成待装调光学系统各自由度调整变量的正态分布的随机数;
1.2)通过DDE建立计算机辅助装调方法与ZEMAX软件的动态数据交互,并通过DDE将各自由度调整变量的正态分布的随机数反馈到ZEMAX软件中;
1.3)以待装调光学系统各自由度调整变量随机数为基准点,计算并建立待装调光学系统随机失调点的灵敏度矩阵数据库。
上述步骤2)的具体实现方式是:
2.1)输入待装调光学系统不同视场的像差测量值;所述待装调光学系统n个视场、m个调节自由度变量、阻尼因子P以及迭代次数j;所述阻尼因子0<P<i,所述i是大于0的整数;
2.2)计算待装调光学系统像差测量值与随机失调点的理论像差值偏差的平方和S;
2.3)寻找满足步骤2.2)中待装调光学系统像差测量值与随机失调点的理论像差值偏差的平方和S的最小的随机失调点,在步骤1)中所建立的待装调光学系统随机失调点的灵敏度矩阵数据库中确定最佳灵敏度矩阵B;
2.4)对步骤2.3)所确定得到的最佳灵敏度矩阵B进行Gram-schmidt正交化处理得到矩阵B′,矩阵B′中列向量彼此正交,保证列之间的非相关性;
2.5)计算待装调光学系统各自由度调整量ΔX;待装调光学系统各自由度调整量ΔX的计算方式是:
B′×ΔX=ΔF
其中:
ΔF为待装调光学系统像差测量值与理论值的偏差;
2.6)通过DDE将步骤2.5)中所计算得到的待装调光学系统各自由度调整量反馈到ZEMAX软件中,获取待装调光学系统在调整后的像差;
2.7)计算待装调光学系统在调整后的像差与理论像差的偏差平方和MFi;
2.8)判断MFi是否最小;若否,阻尼因子P的值加1,重新选择灵敏度因子,并执行步骤2.4)至步骤2.8);若是,确定待装调光学系统最佳调整量和像差残余量;同时将迭代次数j的值加1,同时执行步骤2.4)至步骤2.8);
2.9)完成迭代次数的计算,输出待装调光学系统各自由度调整量的计算结果。
本发明的优点在于:
本发明提供了一种大动态范围光学系统计算机辅助装调方法,该方法通过动态数据交互(DDE)将其与ZEMAX软件进行动态数据交互,同时通过蒙特-卡罗法实现了大动态范围计算机辅助装调;本发明实现了试验与仿真的相结合,从而极大的提高了计算结果的准确度;降低了对待装调光学系统粗装调的要求;重复性好,计算结果置信度高;改善了非线性问题处理时局部解的收敛性,并提高了收敛速度。
附图说明
图1是本发明所提供的大动态范围光学系统计算机辅助装调方法的流程图。
具体实施方式
参见图1,本发明提供了一种大动态范围光学系统计算机辅助装调方法,该方法结合常规的灵敏度矩阵法小范围快速收敛和采用蒙特卡-洛法处理非线性问题动态范围大的优点,具体而言,该方法包括以下步骤:
1)预处理模块:
Step1:通过蒙特-卡罗法生成待装调光学系统各自由度调整变量的正态分布的随机数;
Step2:通过DDE建立计算机辅助装调方法与ZEMAX软件的动态数据交互,并通过DDE将各自由度调整变量的正态分布的随机数反馈到ZEMAX软件中;
Step3:以各自由度调整变量随机数为基准点,计算并建立待装调光学系统随机失调点的灵敏度矩阵数据库;
2)计算主模块:
Step4:输入待装调光学系统不同视场的像差测量值;
Step5:设置视场个数n、调节自由度变量个数m、阻尼因子P的范围为(0~i,i为大于0的整数)和迭代次数j;
Step6:计算待装调光学系统像差测量值与随机失调点的理论像差值偏差的平方和S;
Step7:寻找满足S最小的随机失调点,从而在灵敏度矩阵数据库中确定最佳灵敏度矩阵B;
Step8:对灵敏度矩阵B进行Gram-schmidt正交化处理得到矩阵B′,矩阵B′中列向量彼此正交,保证了列之间的非相关性;
Step9:计算待装调光学系统各自由度调整量ΔX;待装调光学系统各自由度调整量ΔX的计算方式是:
B′×ΔX=ΔF
其中:
ΔF为待装调光学系统像差测量值与理论值的偏差;
Step10:通过DDE将待装调光学系统各自由度调整量反馈到ZEMAX软件中,获取系统调整后的像差;
Step11:计算系统调整后的像差与理论像差的偏差平方和MFi;
Step12:判断MFi是否最小;
Step13:如果MFi不是最小,i加1,重新选择灵敏度因子,并从step8执行到Step12;如果MFi最小,确定待装调光学系统最佳调整量和像差残余量。迭代次数j加1,再从step8执行到Step13;
Step14:完成迭代次数的计算,输出待装调光学系统各自由度调整量的计算结果。
3)根据步骤2)中的装调光学系统各自由度调整量的计算结果对装调光学系统的各自由度进行调整。
通过理论仿真分析,此方法极大的提高了计算机辅助装调的动态范围,并保证了求解结果的收敛性。
Claims (1)
1.一种大动态范围光学系统计算机辅助装调方法,其特征在于:所述大动态范围光学系统计算机辅助装调方法包括以下步骤:
1)计算并建立待装调光学系统随机失调点的灵敏度矩阵数据库;
2)根据步骤1)中所建立的待装调光学系统随机失调点的灵敏度矩阵数据库以及待装调光学系统不同视场的像差测量值获取待装调光学系统各自由度调整量的计算结果;
3)根据步骤2)中的装调光学系统各自由度调整量的计算结果对装调光学系统的各自由度进行调整;
步骤1)的具体实现方式是:
1.1)通过蒙特-卡罗法生成待装调光学系统各自由度调整变量的正态分布的随机数;
1.2)通过DDE建立计算机辅助装调方法与ZEMAX软件的动态数据交互,并通过DDE将各自由度调整变量的正态分布的随机数反馈到ZEMAX软件中;
1.3)以待装调光学系统各自由度调整变量随机数为基准点,计算并建立待装调光学系统随机失调点的灵敏度矩阵数据库;
步骤2)的具体实现方式是:
2.1)输入待装调光学系统不同视场的像差测量值;所述待装调光学系统n个视场、m个调节自由度变量、阻尼因子P以及迭代次数j;所述阻尼因子0<P<i,所述i是大于0的整数;
2.2)计算待装调光学系统像差测量值与随机失调点的理论像差值偏差的平方和S;
2.3)寻找满足步骤2.2)中待装调光学系统像差测量值与随机失调点的理论像差值偏差的平方和S的最小的随机失调点,在步骤1)中所建立的待装调光学系统随机失调点的灵敏度矩阵数据库中确定最佳灵敏度矩阵B;
2.4)对步骤2.3)所确定得到的最佳灵敏度矩阵B进行Gram-schmidt正交化处理得到矩阵B′,矩阵B′中列向量彼此正交,保证列之间的非相关性;
2.5)计算待装调光学系统各自由度调整量ΔX;待装调光学系统各自由度调整量ΔX的计算方式是:
B′×ΔX=ΔF
其中:
ΔF为待装调光学系统像差测量值与理论值的偏差;
2.6)通过DDE将步骤2.5)中所计算得到的待装调光学系统各自由度调整量反馈到ZEMAX软件中,获取待装调光学系统在调整后的像差;
2.7)计算待装调光学系统在调整后的像差与理论像差的偏差平方和MFi;
2.8)判断MFi是否最小;若否,阻尼因子P的值加1,重新选择灵敏度因子,并执行步骤2.4)至步骤2.8);若是,确定待装调光学系统最佳调整量和像差残余量;同时将迭代次数j的值加1,同时执行步骤2.4)至步骤2.8);
2.9)完成迭代次数的计算,输出待装调光学系统各自由度调整量的计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410216269.3A CN104036071B (zh) | 2014-05-21 | 2014-05-21 | 大动态范围光学系统计算机辅助装调方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410216269.3A CN104036071B (zh) | 2014-05-21 | 2014-05-21 | 大动态范围光学系统计算机辅助装调方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104036071A CN104036071A (zh) | 2014-09-10 |
CN104036071B true CN104036071B (zh) | 2017-02-15 |
Family
ID=51466842
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410216269.3A Expired - Fee Related CN104036071B (zh) | 2014-05-21 | 2014-05-21 | 大动态范围光学系统计算机辅助装调方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104036071B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109002567B (zh) * | 2017-06-07 | 2023-02-03 | 中国航空工业集团公司洛阳电光设备研究所 | 基于实际光学系统装调检测光路的计算机辅助装调方法 |
CN107271150B (zh) * | 2017-06-27 | 2019-05-07 | 南京理工大学 | 一种基于干涉测量及灵敏度矩阵的计算机辅助装调方法 |
CN109164591B (zh) * | 2018-08-10 | 2021-01-05 | 宁波舜宇仪器有限公司 | 显微镜物镜计算机辅助装调方法 |
CN112394508A (zh) * | 2020-11-30 | 2021-02-23 | 中国科学院光电技术研究所 | 一种基于二阶灵敏度矩阵法的装调方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902857A (zh) * | 2012-10-08 | 2013-01-30 | 华中科技大学 | 一种光学元件支撑参数的蒙特卡洛分析方法及系统 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003114402A (ja) * | 2001-10-02 | 2003-04-18 | Sumitomo Electric Ind Ltd | 光合分波器およびその調整方法 |
-
2014
- 2014-05-21 CN CN201410216269.3A patent/CN104036071B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902857A (zh) * | 2012-10-08 | 2013-01-30 | 华中科技大学 | 一种光学元件支撑参数的蒙特卡洛分析方法及系统 |
Non-Patent Citations (4)
Title |
---|
三反射镜光学系统的计算机辅助装调技术研究;杨晓飞;《中国优秀博硕士学位论文全文数据库(博士) 工程科技II辑》;20060515(第05期);第C030-2页 * |
基于DDE接口技术的计算机辅助装调方法;刘志祥等;《应用光学》;20090330;第30卷(第3期);第486-490页 * |
计算机辅助装调的代理模型方法;张东阁等;《红外与激光工程》;20130330;第42卷(第3期);第680-685页 * |
面向装配过程的精密光学系统计算机辅助装调技术;薛晓光等;《计算机集成制造系统》;20111030;第17卷(第10期);第2163-2170页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104036071A (zh) | 2014-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104036071B (zh) | 大动态范围光学系统计算机辅助装调方法 | |
US11714174B2 (en) | Parameter calibration method and apparatus of multi-line laser radar, device and readable medium | |
CN106767933A (zh) | 深度相机误差的测量系统、测量方法、评价方法及补偿方法 | |
CN103323846B (zh) | 一种基于极化干涉合成孔径雷达的反演方法及装置 | |
CN104102836B (zh) | 一种电力系统快速抗差状态估计方法 | |
CN103837160B (zh) | 基于恒星匹配的光电经纬仪参量优化方法 | |
CN104990501B (zh) | 一种三维激光扫描装置的系统参数校准方法 | |
Gong et al. | OFLID: Simple method of overlap factor calculation with laser intensity distribution for biaxial lidar | |
KR101651893B1 (ko) | 신경망 회로를 이용한 원자로 노심의 축방향 출력분포 합성방법 및 그 방법이 적용된 노심감시계통 | |
CN109035345A (zh) | 基于高斯过程回归的tof相机距离校正方法 | |
CN107632964A (zh) | 一种平面地磁异常场向下延拓递归余弦变换法 | |
CN104048605B (zh) | 激光扫描测量标靶平差方程式构建方法 | |
CN105823441A (zh) | 一种基于双光敏传感器的光束偏差测量方法 | |
CN103438872A (zh) | 一种基于大坝三维前方交会测量的内外业一体化系统 | |
CN104316082A (zh) | 一种经纬仪外场无穷远距离校正方法 | |
CN103148779B (zh) | 位置测量设备中光源的调整装置 | |
CN103092815B (zh) | 对监测装置中的传递函数进行校准的方法 | |
CN108844647B (zh) | 一种确定温度传感器最佳安装位置的方法 | |
CN112013964A (zh) | 一种红外测温距离自动补偿方法、系统、装置和储存介质 | |
CN103674000B (zh) | 电子罗盘实时校准算法 | |
CN113820720B (zh) | 基于多参考基点的三维激光中心测距方法、系统及终端 | |
CN105387805A (zh) | 基于远场光斑能量曲线的拼接误差检测方法 | |
CN105136071B (zh) | 一种修正有效面积法测量数据的方法 | |
CN110474315A (zh) | 一种用于柔性直流小干扰稳定性判定的方法及系统 | |
CN103926561B (zh) | 一种用于超短基线安装误差校准的奇异值消除的参数估计权值设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 Termination date: 20180521 |
|
CF01 | Termination of patent right due to non-payment of annual fee |