CN114545621B - 一种基于遗传算法的变形反射镜结构优化方法 - Google Patents
一种基于遗传算法的变形反射镜结构优化方法 Download PDFInfo
- Publication number
- CN114545621B CN114545621B CN202210252231.6A CN202210252231A CN114545621B CN 114545621 B CN114545621 B CN 114545621B CN 202210252231 A CN202210252231 A CN 202210252231A CN 114545621 B CN114545621 B CN 114545621B
- Authority
- CN
- China
- Prior art keywords
- deformable mirror
- deformed
- driver
- reflector
- rms
- 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.)
- Active
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 47
- 230000002068 genetic effect Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000005457 optimization Methods 0.000 title claims description 7
- 230000004075 alteration Effects 0.000 claims abstract description 62
- 230000000694 effects Effects 0.000 claims abstract description 35
- 201000009310 astigmatism Diseases 0.000 claims description 14
- 230000002093 peripheral effect Effects 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 11
- 239000010410 layer Substances 0.000 description 24
- 238000010586 diagram Methods 0.000 description 8
- 230000035772 mutation Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 239000002356 single layer Substances 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000002301 combined effect Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000003702 image correction Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000004083 survival effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B26/00—Optical devices or arrangements for the control of light using movable or deformable optical elements
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/0025—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for optical correction, e.g. distorsion, aberration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Biophysics (AREA)
- Optics & Photonics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Medical Informatics (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Optical Elements Other Than Lenses (AREA)
- Mechanical Light Control Or Optical Switches (AREA)
Abstract
本发明公开一种基于遗传算法的变形反射镜结构优化方法,包括:1)根据弹性力学变形原理确定影响变形反射镜发生形变的结构参数;2)根据变形反射镜结构要求初步确定结构参数范围;3)编写子函数,以变形反射镜结构参数为自变量,根据弹性力学方程得到子驱动影响函数,以标准Zernike像差为拟合目标,计算不同结构参数组合下变形反射镜拟合效果;4)编写遗传算法,以变形反射镜对Zernike像差的拟合残差值作为适应度函数,将变形反射镜结构参数作为算法自变量;5)利用遗传算法进行迭代计算得到变形反射镜最优结构参数组合;6)利用迭代后参数值进一步确定变形反射镜结构参数范围,循环计算至变形反射镜结构满足最终校正拟合要求。本发明能够提高变形反射镜设计效率。
Description
技术领域
本发明适用于自适应光学领域,具体涉及一种基于遗传算法的变形反射镜设计优化方法。
背景技术
对于自适应光学而言,变形反射镜作为系统波前校正器的一种,是自适应光学系统中不可或缺的,而变形反射镜对于像差校正的效果主要取决于变形反射镜结构,故确定变形反射镜的结构及其参数变得至关重要。
目前现有变形反射镜的结构设计过程如下:根据经验提出变形反射镜大致结构,随后建立变形反射镜三维模型,利用有限元分析形变,得到像差校正效果,根据像差校正效果对变形反射镜结构参数逐个进行单一变化直至达到像差校正要求。针对该设计过程,存在一个比较严重的问题,即变形反射镜各结构参数变化对校正效果影响不确定,且各变量之间相互影响造成变形反射镜结构参数修改后校正效果趋势不定,需连续不断重复修改模型参数、仿真、计算使工作量变得巨大,工作效率极其低下。针对这个问题,本发明提出利用遗传算法结合弹性力学来快速寻找满足变形反射镜像差校正要求下的结构及其参数。
遗传算法是一种基于自然选择和群体遗传机理的搜索算法,模拟自然选择和自然遗传过程中的繁殖、杂交和突变现象。遗传算法首先随机产生一些个体(即初始解),根据预定的适应度函数对每一个个体进行评估,给出一个适应度值,基于此适应度值,选择一些个体用来产生下一代,选择操作体现了“适者生存”的原理,“好”的个体被用来产生下一代,“坏”的个体则被淘汰,然后选择出来的个体,经过交叉和变异算子进行再组合生成新的一代,这一代的个体由于继承了上一代的一些优良性状,因而在性能上要优于上一代,这样逐步朝着最优解的方向进化。
因此,具体遗传算法可以看成是一个由可行解组成的群体逐步进化的过程。利用遗传算法可以探讨驱动器排布以及变形反射镜结构参数对像差校正效果的影响,从而快速得到像差校正能力较好的变形反射镜结构参数值。
发明内容
本发明要解决的技术问题为:克服现有变形反射镜结构参数样本量巨大导致工作效率低下问题,提供一种基于遗传算法的变形反射镜设计方法,给出变形反射镜驱动器排布对变形反射镜像差校正能力影响,快速确定变形反射镜结构参数的同时满足其像差校正要求。
本发明解决上述技术问题采用的技术方案是:
一种基于遗传算法的变形反射镜结构优化方法,根据变形反射镜待校正的Zernike像差面形初步确定变形反射镜驱动器排布,以变形反射镜结构参数作为遗传算法自变量,拟合效果作为算法适应度函数进行结构参数优化,得到满足像差校正要求的结构参数组合。
具体的步骤如下:
步骤一:根据弹性力学变形原理确定影响变形反射镜发生形变的结构参数;
步骤二:根据变形反射镜结构要求初步确定结构参数范围;
步骤三:编写子函数,以变形反射镜结构参数为自变量,根据弹性力学方程得到子驱动影响函数,以标准Zernike像差为拟合目标,计算不同结构参数组合下变形反射镜拟合效果;
步骤四:编写遗传算法,以变形反射镜对Zernike像差的拟合残差值作为适应度函数,将变形反射镜结构参数作为算法自变量;
步骤五:利用遗传算法进行迭代计算得到变形反射镜最优结构参数组合;
步骤六:利用迭代后参数值进一步确定变形反射镜结构参数范围,循环计算至变形反射镜结构满足最终拟合要求。
进一步地,步骤一中影响变形反射镜拟合能力的结构参数有:变形反射镜整体面形尺寸、变形反射镜厚度、变形反射镜驱动器截面尺寸、变形反射镜固定点位置、驱动器位置。
进一步地,步骤二中变形反射镜结构要求为两点,分别是变形反射镜有效校正面积尺寸与变形反射镜结构需满足弹性薄板要求。其中弹性薄板要求为变形反射镜厚度与横向最小尺寸比值介于0.0125-0.2。
进一步地,步骤三中影响函数指变形反射镜在单个驱动器下以单位力或者单位位移驱动得到的变形反射镜面形。驱动器位置不同,影响函数不同,根据弹性力学影响函数分两种,一种为驱动器在中心位置的影响函数,另一种为驱动器处于偏心位置下的影响函数,分别为公式(1)、(2):
其中,k=r/a,r为所求点到中心的距离,a为圆板半径,K=a/b,b为载荷圆中心到圆板中心的距离;m=a/t,t为圆板厚度,q为施加在中心载荷圆上均布载荷,P为施加在偏心载荷圆上均布载荷,E为弹性模量,为两条线之间夹角,一条线为偏心载荷圆中心与圆板中心连线;另一条线为所求点到圆板中心的连线,A0,A1,A2,b0,b1,b2,c0,c1,c2为相关系数,A0,A1,A2可由K、k表示,具体公式如下;
A0=0.19(K3-2.65K2+1.65)
A1=0.822(K4-1.95K3+0.955K)
A2=1.30(K4-2.47K3+1.47K2)
b0=2.65,b1=1.95,b2=2.47
c0=1.65,c1=0.955,c2=1.47
进一步地,步骤三中影响变形反射镜拟合效果主要有两个大方向的影响因素,分别为驱动器排布以及变形反射镜结构参数。驱动器排布主要会改变影响函数的叠加状态,而变形反射镜结构参数主要影响变形反射镜单驱动器下的影响函数。其中,在变形反射镜相同结构参数范围下,对比不同驱动器排布和数量对变形反射镜拟合效果的影响;在变形反射镜同一驱动器排布下,对比变形反射镜结构参数值对变形反射镜拟合效果的影响。
所述变形反射镜驱动器排布主要为中心驱动器加周边环布驱动器;周边环布驱动器由一层或多层驱动环组成,每一层驱动器数为均布4点或均布8点。
所述驱动器数量最少为9个,最多至41个。
所述驱动器外部层数最少为1层,最多为9层。
所述影响函数的叠加状态为不同驱动系数下驱动器的影响函数面形的线性叠加。
进一步地,步骤三中所述变形反射镜所校正的Zernike像差为离焦和像散像差。
进一步地,步骤三中所述变形反射镜的像差校正效果由拟合后的面形残差RMS值与原始面形RMS值的比值判定。
进一步地,步骤四中所述适应度函数为变形反射镜镜像校正后面形残差RMS与Zernike像差面形RMS比值的和,具体公式为:
Final(i)=rms_d/rms_defocus+rms_a/rms_astigmatism (3)
其中,rms_d为变形反射镜校正离焦像差后的残差面形RMS值;rms_defocus为待拟合Zernike离焦像差面形RMS值;rms_a为变形反射镜像散像差后的残差面形RMS值;rms_astigmatism为待拟合Zernike像散像差面形RMS值。
进一步地,步骤四中遗传算法的过程主要为:初始参数设定,适应度计算,选择、遗传、变异,判断是否满足迭代设定,满足则输出变量值否则返回适应度计算进而选择、遗传、变异,判断循环直至结束。
本发明的有益效果在于:
(1)本发明利用遗传算法和弹性力学对变形反射镜结构参数进行设计优化的方法,可以快速地确定变形反射镜的驱动器数和驱动器排布以及某种排布方式下最合适的变形反射镜结构参数组,极大地提高变形反射镜设计阶段的工作效率,为变形反射镜的设计思路上提供一定的参考价值,在变形反射镜结构设计方面具有较好的应用场景。
(2)本发明该方法克服了传统变形反射镜设计下控制变量法及二分法中工作量大,效率低下的缺点,极大地提高了工作效率。
(3)本发明该方法的拟合像差不拘泥于Zernike某几种像差,可根据实际情况确定遗传算法适应度函数,快速得到满足拟合要求的变形反射镜结构及参数。
附图说明
图1为本发明基于遗传算法的变形反射镜结构优化方法流程示意图;
图2为实施例中遗传算法流程示意图;
图3为实施例中一种变形反射镜单4驱动层驱动分布图;
图4为实施例中一种变形反射镜单8驱动层驱动分布图;
图5为实施例中一种变形反射镜混合驱动层驱动分布图;
图6为实施例中相同参数取值、不同驱动器排布、相同驱动器层数下像差校正残差比值变化图;
图7为实施例中相同参数取值、相同驱动器排布、不同驱动器层数下像差校正残差比值变化图;
图8为实施例中不同参数取值、相同驱动器排布、相同驱动器层数下像差校正残差比值变化图。
具体实施方式
下面结合附图及具体实施方式详细介绍本发明。但以下的实施例仅限于解释本发明,本发明的保护范围应包括权利要求的全部内容,而且通过以下实施例,本领域技术人员即可以实现本发明权利要求的全部内容。
本发明一种基于遗传算法的变形反射镜结构优化方法,根据变形反射镜待校正的Zernike像差面形初步确定变形反射镜驱动器排布,以变形反射镜结构参数作为遗传算法自变量,拟合效果作为算法适应度函数进行结构参数优化,得到满足像差校正要求的结构参数组合。
具体的步骤如下:
步骤一:根据弹性力学变形原理确定影响变形反射镜发生形变的结构参数;
步骤一中影响变形反射镜拟合能力的结构参数有:变形反射镜整体面形尺寸、变形反射镜厚度、变形反射镜驱动器截面尺寸、变形反射镜固定点位置、驱动器位置。
步骤二:根据变形反射镜结构要求初步确定结构参数范围;
步骤二中变形反射镜结构要求为两点,分别是变形反射镜有效校正面积尺寸与变形反射镜结构需满足弹性薄板要求。其中弹性薄板要求为变形反射镜厚度与横向最小尺寸比值介于0.0125-0.2。
步骤三:编写子函数,以变形反射镜结构参数为自变量,根据弹性力学方程得到子驱动影响函数,以标准Zernike像差为拟合目标,计算不同结构参数组合下变形反射镜拟合效果;
步骤三中影响函数指变形反射镜在单个驱动器下以单位力或者单位位移驱动得到的变形反射镜面形。驱动器位置不同,影响函数不同,根据弹性力学影响函数分为两种,一种为驱动器在中心位置的影响函数,另一种为驱动器处于偏心位置下的影响函数,分别为公式(1)、(2):
其中,k=r/a,r为所求点到圆板中心的距离,a为圆板半径,K=a/b,b为载荷圆中心到圆板中心的距离;m=a/t,t为圆板厚度,q为施加在中心载荷圆上均布载荷,P为施加在偏心载荷圆上均布载荷,E为弹性模量,为两条线之间夹角,一条线为偏心载荷圆中心与圆板中心连线;另一条线为所求点到圆板中心的连线;A0,A1,A2,b0,b1,b2,c0,c1,c2为相关系数,A0,A1,A2可由K、k表示,具体公式如下:
A0=0.19(K3-2.65K2+1.65)
A1=0.822(K4-1.95K3+0.955K)
A2=1.30(K4-2.47K3+1.47K2)
b0=2.65,b1=1.95,b2=2.47
c0=1.65,c1=0.955,c2=1.47
步骤三中影响变形反射镜拟合效果主要有两个大方向的影响因素,分别为驱动器排布以及变形反射镜结构参数。驱动器排布主要会改变影响函数的叠加状态,而变形反射镜结构参数主要影响变形反射镜单驱动器下的影响函数。其中,在变形反射镜相同结构参数范围下,对比不同驱动器排布和数量对变形反射镜拟合效果的影响;在变形反射镜同一驱动器排布下,对比变形反射镜结构参数值对变形反射镜拟合效果的影响。
所述变形反射镜驱动器排布主要为中心驱动器加周边环布驱动器;周边环布驱动器由一层或多层驱动环组成,每一层驱动器数为均布4点或均布8点。
所述驱动器数量最少为9个,最多至41个。
所述驱动器外部层数最少为1层,最多为9层。
所述影响函数的叠加状态为不同驱动系数下驱动器的影响函数面形的线性叠加。
步骤三中所述变形反射镜所校正的Zernike像差为离焦和像散像差。
步骤三中所述变形反射镜的像差校正效果由拟合后的面形残差RMS值与原始面形RMS值的比值判定。
步骤四:编写遗传算法,以变形反射镜对Zernike像差的拟合残差值作为适应度函数,将变形反射镜结构参数作为算法自变量;
步骤四中所述适应度函数为变形反射镜镜像校正后面形残差RMS与Zernike像差面形RMS比值的和,具体公式为:
Final(i)=rms_d/rms_defocus+rms_a/rms_astigmatism (3)
其中,rms_d为变形反射镜校正离焦像差后的残差面形RMS值;rms_defocus为待拟合Zernike离焦像差面形RMS值;rms_a为变形反射镜像散像差后的残差面形RMS值;rms_astigmatism为待拟合Zernike像散像差面形RMS值。
步骤四中遗传算法的过程主要为:初始参数设定,适应度计算,选择、遗传、变异,判断是否满足迭代设定,满足则输出变量值否则返回适应度计算进而选择、遗传、变异,判断循环直至结束。
步骤五:利用遗传算法进行迭代计算得到变形反射镜最优结构参数组合;
步骤五中所述变形反射镜最优结构组合是根据算法收敛稳定后得到的结构参数组合,是在特定驱动器排布和特定变形反射镜参数范围下得到的最优结构参数组合。
步骤六:利用迭代后参数值进一步确定变形反射镜结构参数范围,循环计算至变形反射镜结构满足最终拟合要求。
步骤六中所述迭代后参数值是指步骤五中不同驱动器排布的所有最优结构参数组合,根据结构参数组合下适应度函数值的变化趋势,适当改变变形反射镜结构参数范围,直到适应度函数值小于变形反射镜所要求的值。
实施例
一种基于遗传算法的变形反射镜参数确定流程如图1所示。具体的流程为:首先确定变形反射镜需要校正的Zernike像差,根据Zernike像差的面形特征与数值确定变形反射镜驱动器的排布方式与遗传算法的适应度函数,其中驱动器排布分别为图3、图4、图5,适应度函数函数则以各个面形拟合残差比的关系确定;其次则是根据弹性力学要求和变形反射镜结构要求确定影响变形反射镜像差校正能力的参数及其范围;确定结构参数范围后,选择图3、图4、图5中的驱动器排布,利用遗传算法得到相关变形反射镜结构参数,随后将结构参数值进行存储并分析变形反射镜驱动器排布与驱动器数量对变形反射镜拟合效果的影响,从而进一步确定变形反射镜的具体驱动器排布方式;最后观察变形反射镜的校正残差比值是否达到像差校正要求,达到则输出变形反射镜驱动排布方式及其结构参数值,未达到则根据结构参数值趋势进一步修改变形反射镜参数取值范围观察变形反射镜像差校正效果,直至变形反射镜像差校正效果达到要求。
变形反射镜驱动器排布根据离焦和像散像差的面形特征确定,除中心驱动器外,外层驱动层有两种,一种为单层环形均布4个,一种为单层环形均布8个,根据每层的驱动器数可将驱动器排布分为三种:驱动层均为单层为4驱动器的驱动器排布如图3;驱动层均为单层为8驱动器的驱动器排布如图4;以及驱动器驱动层中既有单层4驱动既有单层8驱动的驱动器排布如图5。
周期性微透镜列阵4中的各个子透镜单元的周期性分布,其子透镜单元的口径、焦距以及阵列数与中心离轴型微透镜列阵中的子透镜单元相同,分别为500μm、47.48mm和20×20,驱动分布图如图4所示。
如图6所示,在变形反射镜结构取相同参数取值、不同驱动器排布、相同驱动器层数下像差校正残差比值变化不大,对于此处的离焦和像散像差而言,真正影响变形反射镜像差校正效果的驱动位置不多,只要特别重要的驱动器确定,其他位置的驱动器则无关紧要,可视为失效驱动器。
如图7所示,在变形反射镜结构取相同参数取值、相同驱动器排布、不同驱动器层数下像差校正残差比值变化不大,和图6结果相似,主要的驱动器确定后,即使在原有基础上增加驱动器数目,只会影响变形反射镜各个驱动器上乘积系数而对于变形反射镜整体像差校正效果影响不大。
如图8所示,在变形反射镜结构取不同参数取值、相同驱动器排布、相同驱动器层数下像差校正残差比值呈现一定的趋势变化。在变形反射镜有效校正面积一定的情况下,适当的增大变形反射镜固定点离中心位置以及中心驱动器作用面积会使得变形反射镜的拟合能力得到大幅度提升,但是随着结构上固定点离中心距离越大,变形反射镜的拟合优化效果逐步变小,而中心驱动器作用面在增长到一定程度后则不会再增大。
结合数据分析,变形反射镜对于离焦和像散像差残差占比值均在10%以下,根据线性叠加原理,则可以保证变形反射镜对于任何可由离焦、像散像差线性组合的面形的校正残差占比均在10%以下,视为达到完全校正。
本发明结合弹性力学和遗传算法得到影响变形反射镜的拟合能力因素,且明确展示了各因素对变形反射镜像差拟合能力的影响效果,为变形反射镜前期的结构设计提供了一个参考;且该方法可以在确定变形反射镜的校正要求下快速确定变形反射镜的排布方式及其结构的具体参数值,很大程度上的减小了变形反射镜结构设计的工作量,极大地提高了工作效率。
Claims (6)
1.一种基于遗传算法的变形反射镜结构优化方法,其特征在于:根据变形反射镜待校正的Zernike像差面形初步确定变形反射镜驱动器排布,以变形反射镜结构参数作为遗传算法自变量,拟合效果作为算法适应度函数进行结构参数优化,得到满足像差校正要求结构参数组合;具体的步骤如下:
步骤一:根据弹性力学变形原理确定影响变形反射镜发生形变的结构参数;
步骤二:根据变形反射镜结构要求初步确定结构参数范围;其中变形反射镜结构要求为两点,分别是变形反射镜有效校正面积尺寸与变形反射镜结构需满足弹性薄板要求,其中弹性薄板要求为变形反射镜厚度与横向最小尺寸比值介于0.0125-0.2;
步骤三:编写子函数,以变形反射镜结构参数为自变量,根据弹性力学方程得到子驱动影响函数,以标准Zernike像差为拟合目标,计算不同结构参数组合下变形反射镜拟合效果;其中影响函数指变形反射镜在单个驱动器下以单位力或者单位位移驱动得到的变形反射镜面形,驱动器位置不同,影响函数不同,根据弹性力学影响函数分为两种,一种为驱动器在中心位置的影响函数,另一种为驱动器处于偏心位置下的影响函数,分别为公式(1)、(2):
其中,k=r/a,r为所求点到圆板中心的距离,a为圆板半径,K=a/b,b为载荷圆中心到圆板中心的距离;A0,A1,A2,b0,b1,b2,c0,c1,c2为相关系数,A0,A1,A2可由K、k表示;m=a/t,t为圆板厚度;q为施加在中心载荷圆上均布载荷,P为施加在偏心载荷圆上均布载荷,E为弹性模量,为两条线之间夹角,一条线为偏心载荷圆中心与圆板中心连线;另一条线为所求点到圆板中心连线,A0,A1,A2,b0,b1,b2,c0,c1,c2为相关系数,A0,A1,A2可由K、k表示,具体公式如下;
A0=0.19(K3-2.65K2+1.65)
A1=0.822(K4-1.95K3+0.955K)
A2=1.30(K4-2.47K3+1.47K2)
b0=2.65,b1=1.95,b2=2.47
c0=1.65,c1=0.955,c2=1.47
步骤四:编写遗传算法,以变形反射镜对Zernike像差的拟合残差值作为适应度函数评判值,将变形反射镜结构参数作为算法自变量;
步骤五:利用遗传算法进行迭代计算得到变形反射镜最优结构参数组合;
步骤六:利用迭代后参数值进一步确定变形反射镜结构参数范围,循环计算至变形反射镜结构满足最终拟合要求。
2.根据权利要求1所述的一种基于遗传算法的变形反射镜结构优化方法,其特征在于:步骤一中影响变形反射镜拟合能力的结构参数有:变形反射镜整体面形尺寸、变形反射镜厚度、变形反射镜驱动器截面尺寸、变形反射镜固定点位置、驱动器位置。
3.根据权利要求1所述的一种基于遗传算法的变形反射镜结构优化方法,其特征在于:步骤三中影响变形反射镜拟合效果主要有两个大方向的影响因素,分别为驱动器排布以及变形反射镜结构参数,驱动器排布主要会改变影响函数的叠加状态,而变形反射镜结构参数主要影响变形反射镜单驱动器下的影响函数,为此,在变形反射镜相同结构参数范围下,对比不同驱动器排布和数量对变形反射镜拟合效果的影响;在变形反射镜同一驱动器排布下,对比变形反射镜结构参数值对变形反射镜拟合能力的影响;
所述变形反射镜驱动器排布主要为中心驱动器加周边环布驱动器;周边环布驱动器由一层或多层驱动环组成,每一层驱动器数为均布4点或均布8点;
所述驱动器数量最少为9个,最多至41个;
所述驱动器外部层数最少为1层,最多为9层;
所述影响函数的叠加状态为不同驱动系数下驱动器的影响函数面形的线性叠加。
4.根据权利要求1所述的一种基于遗传算法的变形反射镜结构优化方法,其特征在于:步骤三中所述变形反射镜所校正的Zernike像差为离焦和像散像差;
步骤三中所述变形反射镜的像差校正效果由拟合后的面形残差RMS值与原始面形RMS值的比值判定。
5.根据权利要求1所述的一种基于遗传算法的变形反射镜结构优化方法,其特征在于:步骤四中所述适应度函数为变形反射镜像差校正后面形残差RMS与Zernike像差面形RMS比值的和,具体公式为:
Final(i)=rms_d/rms_defocus+rms_a/rms_astigmatism(3)
其中,rms_d为变形反射镜校正离焦像差后的残差面形RMS值;rms_defocus为待拟合Zernike离焦像差面形RMS值;rms_a为变形反射镜校正像散像差后的残差面形RMS值;rms_astigmatism为待拟合Zernike像散像差面形RMS值。
6.根据权利要求1所述的一种基于遗传算法的变形反射镜结构优化方法,其特征在于:步骤四中遗传算法的过程主要为:初始参数设定,适应度计算,选择、遗传、变异,判断是否满足迭代设定,满足则输出变量值否则返回适应度计算进而选择、遗传、变异,判断循环直至结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210252231.6A CN114545621B (zh) | 2022-03-15 | 2022-03-15 | 一种基于遗传算法的变形反射镜结构优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210252231.6A CN114545621B (zh) | 2022-03-15 | 2022-03-15 | 一种基于遗传算法的变形反射镜结构优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114545621A CN114545621A (zh) | 2022-05-27 |
CN114545621B true CN114545621B (zh) | 2023-09-22 |
Family
ID=81663839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210252231.6A Active CN114545621B (zh) | 2022-03-15 | 2022-03-15 | 一种基于遗传算法的变形反射镜结构优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114545621B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115079390B (zh) * | 2022-06-24 | 2023-05-02 | 中国科学院西安光学精密机械研究所 | 一种大形变高面形精度变曲率反射镜结构参数优化方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1904712A (zh) * | 2006-07-28 | 2007-01-31 | 中国科学院光电技术研究所 | 一种基于遗传算法的自适应校正激光器像差的装置 |
CN1987547A (zh) * | 2006-12-30 | 2007-06-27 | 中国科学院光电技术研究所 | 一种利用望远镜次镜自动校正望远镜像差的装置 |
CN1987546A (zh) * | 2006-12-30 | 2007-06-27 | 中国科学院光电技术研究所 | 一种利用变形镜实现光束自动整形的装置 |
JP2009086248A (ja) * | 2007-09-28 | 2009-04-23 | Nippon Hoso Kyokai <Nhk> | 光制御装置 |
JP2010040161A (ja) * | 2008-07-11 | 2010-02-18 | Nippon Hoso Kyokai <Nhk> | ホログラム記録装置および遺伝的アルゴリズムを用いたホログラムの参照光の制御方法 |
CN102722023A (zh) * | 2012-05-29 | 2012-10-10 | 中国科学院光电技术研究所 | 一种基于分离模式控制和优化算法的多变形镜控制方法 |
CN103901617A (zh) * | 2014-04-17 | 2014-07-02 | 淮海工学院 | 一种基于模型的无波前探测自适应光学系统 |
CN111580271A (zh) * | 2020-06-24 | 2020-08-25 | 哈工大机器人(中山)无人装备与人工智能研究院 | 自适应像差校正方法及基于该方法的光片显微成像装置 |
CN112363315A (zh) * | 2020-11-06 | 2021-02-12 | 长春理工大学 | 基于泽尼克模型系数优化的非球面窗口光学像差校正方法 |
WO2021035911A1 (zh) * | 2019-08-28 | 2021-03-04 | 青岛蓝海未来海洋科技有限责任公司 | 一种基于数据正反向驱动线性变参数遗传算法的无人艇路径规划方法及系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2003067306A1 (ja) * | 2002-02-06 | 2005-06-02 | 独立行政法人産業技術総合研究所 | 可変形ミラーの調整方法、調整装置および、その調整方法の実施のための処理プログラム |
CN104769481B (zh) * | 2012-10-12 | 2018-12-18 | 统雷有限公司 | 紧凑、低色散以及低像差自适应光学扫描系统 |
MX361245B (es) * | 2013-03-05 | 2018-11-30 | Halliburton Energy Services Inc | Sistema, metodo y producto de programa informatico para fortalecimiento con respecto al entorno y diseño de sistema fotometrico. |
JP6456085B2 (ja) * | 2014-09-25 | 2019-01-23 | キヤノン株式会社 | 形状可変鏡システム、その制御方法および眼科装置 |
DE112018002600B4 (de) * | 2017-05-19 | 2023-08-10 | Kawasaki Jukogyo Kabushiki Kaisha | Adaptiv-optische vorrichtung, optisches system und optisches wellenfrontkompensationsverfahren |
-
2022
- 2022-03-15 CN CN202210252231.6A patent/CN114545621B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1904712A (zh) * | 2006-07-28 | 2007-01-31 | 中国科学院光电技术研究所 | 一种基于遗传算法的自适应校正激光器像差的装置 |
CN1987547A (zh) * | 2006-12-30 | 2007-06-27 | 中国科学院光电技术研究所 | 一种利用望远镜次镜自动校正望远镜像差的装置 |
CN1987546A (zh) * | 2006-12-30 | 2007-06-27 | 中国科学院光电技术研究所 | 一种利用变形镜实现光束自动整形的装置 |
JP2009086248A (ja) * | 2007-09-28 | 2009-04-23 | Nippon Hoso Kyokai <Nhk> | 光制御装置 |
JP2010040161A (ja) * | 2008-07-11 | 2010-02-18 | Nippon Hoso Kyokai <Nhk> | ホログラム記録装置および遺伝的アルゴリズムを用いたホログラムの参照光の制御方法 |
CN102722023A (zh) * | 2012-05-29 | 2012-10-10 | 中国科学院光电技术研究所 | 一种基于分离模式控制和优化算法的多变形镜控制方法 |
CN103901617A (zh) * | 2014-04-17 | 2014-07-02 | 淮海工学院 | 一种基于模型的无波前探测自适应光学系统 |
WO2021035911A1 (zh) * | 2019-08-28 | 2021-03-04 | 青岛蓝海未来海洋科技有限责任公司 | 一种基于数据正反向驱动线性变参数遗传算法的无人艇路径规划方法及系统 |
CN111580271A (zh) * | 2020-06-24 | 2020-08-25 | 哈工大机器人(中山)无人装备与人工智能研究院 | 自适应像差校正方法及基于该方法的光片显微成像装置 |
CN112363315A (zh) * | 2020-11-06 | 2021-02-12 | 长春理工大学 | 基于泽尼克模型系数优化的非球面窗口光学像差校正方法 |
Non-Patent Citations (4)
Title |
---|
.Optimization esign and error distribution for secondary mirror adjusting mechanism of large optical payload.《Optics and precision engineering》.2016,第 24 卷(第 5 期),全文. * |
Han chunyang 等 * |
杨平 ; 许冰 ; 姜文汉 ; 陈善球 ; .遗传算法在自适应光学系统中的应用.《光学学报》.2007,(第09期),全文. * |
杨慧珍 ; 李新阳 ; 姜文汉 ; .自适应光学系统几种随机并行优化控制算法比较.《强激光与粒子束》.2018,(第01期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN114545621A (zh) | 2022-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105988213B (zh) | 自由曲面离轴光学系统的设计方法 | |
CN108152948B (zh) | 离轴非球面光学系统的设计方法 | |
CN114545621B (zh) | 一种基于遗传算法的变形反射镜结构优化方法 | |
CN110490326A (zh) | 一种自适应多变异策略的差分进化方法及系统 | |
CN113900249B (zh) | 基于粒子群算法的同轴反射式光学系统初始结构求解方法 | |
CN111507049B (zh) | 一种镜头像差仿真及优化方法 | |
TWI622800B (zh) | 離軸非球面三反光學系統 | |
CN102313989A (zh) | 一种单个非球面透镜像差校正的方法 | |
CN116009246A (zh) | 一种基于深度学习的偏振光学系统自动优化设计方法 | |
CN117332694A (zh) | 光学系统多目标优化方法 | |
CN113591621A (zh) | 基于改进的多尺度自适应特征融合的DC-YOLOv4算法 | |
TW202144849A (zh) | 自由曲面成像系統的設計方法 | |
CN116430583A (zh) | 基于预测神经网络的消色差超透镜设计方法及存储介质 | |
Al-Sakran et al. | Automated re-invention of a previously patented optical lens system using genetic programming | |
CN112363315B (zh) | 基于泽尼克模型系数优化的非球面窗口光学像差校正方法 | |
CN115657303B (zh) | 双级联超构透镜的设计方法及长波红外超构透镜成像系统 | |
CN116626885A (zh) | 一种超构光栅的设计方法 | |
WO2003067306A1 (fr) | Procede et dispositif de reglage d'un miroir de forme variable et programme de traitement servant a mettre en application ce procede de reglage | |
CN116362125A (zh) | 一种基于深度学习的涡扇发动机多级预测性维护方法 | |
Tsai et al. | Suppression of primary chromatic aberration by genetic algorithm in an advanced telephoto lens | |
CN118821603A (zh) | 基于深度学习的准零刚度低频隔振结构逆向设计方法 | |
Jones et al. | Automated synthesis of both the topology and numerical parameters for seven patented optical lens systems using genetic programming | |
CN118493405B (zh) | 基于关节独立决策的机械臂逆运动学优化求解方法及系统 | |
CN116027547A (zh) | 超表面透镜的逆向设计方法、系统、设备及介质 | |
CN115236848A (zh) | 一种全口径高精度变形镜控制方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |