[已解决]CV的STP优化导致出错？

隔壁的薛定谔

如图所示，一个很简单的模型(EPD= 20mm,单波长，0视场等均为默认值)。
运行以下优化指令

AUT
EFL=145
MNT 3
MXT 6

STP Y

DRA Y
GO


之后，则会出现错误，如下，透镜厚度、焦距等都不受控制。这是什么原因？STP我理解为 对某个变量 每次改变的步长更大一些 所以收敛更快一些，但是从目前的情况来看，STP完全跳出了我的约束范围，导致了错误。






另外，运行完上述优化侯，查看命令行，提示"NO.4" 变量，这是哪个变量？  这里的变量指的是曲率半径和厚度吗？ 这些变量是如何编号的？

最佳答案

月排行榜 / 总排行榜

dpstt

2024-3-18 09:06

初始的优化不建议用STP，个人感觉是优化差不多要Merit Function收敛时再用，STP有一定局限性。如果你在前期想快速收敛可以试试ROU。
变量4估计是像面厚度，可以用der lis指令列出变量及其定义

跳转到最佳答案楼层

dpstt

初始的优化不建议用STP，个人感觉是优化差不多要Merit Function收敛时再用，STP有一定局限性。如果你在前期想快速收敛可以试试ROU。
变量4估计是像面厚度，可以用der lis指令列出变量及其定义

隔壁的薛定谔

dpstt 发表于 2024-3-18 09:06
初始的优化不建议用STP，个人感觉是优化差不多要Merit Function收敛时再用，STP有一定局限性。如果你在前期 ...

1. stp导致程序报错中断优化，等等都可以理解，但是现在STP导致的是我的约束不受控制，那么这个情况是否正常？另外初期优化不建议使用STP那意思说是当系统的评价函数比较小的时候用STP才比较有效果吗？我看说明书上建议stp在早期优化中使用。说明书如下
Step optimization works best in the early stages of optimization when the starting design is relatively far from an error function minimum.


2. 用DER list可以展示出目前的变量，那么这些变量在优化的所有cycle中的排列序号都是固定的吗？

dpstt

隔壁的薛定谔 发表于 2024-3-18 09:53
1. stp导致程序报错中断优化，等等都可以理解，但是现在STP导致的是我的约束不受控制，那么这个情况是否 ...

1、约束不受控制是会遇到的，需要根据情况调整目标及变量。只是个人体会和建议，怎么用还是看你自己的用法。

2、固定的

隔壁的薛定谔

dpstt 发表于 2024-3-18 10:02
1、约束不受控制是会遇到的，需要根据情况调整目标及变量。只是个人体会和建议，怎么用还是看你自己的用法 ...

约束分为拉格朗日乘子约束和权重约束。

拉格朗日乘子约束我理解为是一定会满足的，也就是说在满足拉格朗日乘子的约束前提下 再去寻找MF极小值点。
然而我那一段优化代码中，本质上是拉格朗日乘子约束，例如EFL=145;MNT 3;MXT 6;都是拉格朗日乘子约束。


您说的约束不受控制，指的是权重约束吧？例如WTC给权重。这种约束 可能会无效，需要加大权重等。

大佬，指导下上述理解。。。

dpstt

隔壁的薛定谔 发表于 2024-3-18 10:59
约束分为拉格朗日乘子约束和权重约束。

拉格朗日乘子约束我理解为是一定会满足的，也就是说在满足拉格朗 ...

理论上是这样的，实际在用的时候不是所有条件一定会满足，会有拉格朗日乘子约束被突破的情况。

隔壁的薛定谔

dpstt 发表于 2024-3-18 11:07
理论上是这样的，实际在用的时候不是所有条件一定会满足，会有拉格朗日乘子约束被突破的情况。 ...

确认一下，您要表达的是“拉格朗日乘子被突破”  
而不是  “无法满足拉格朗日乘子约束，因此优化报错”

dpstt

隔壁的薛定谔 发表于 2024-3-18 11:17
确认一下，您要表达的是“拉格朗日乘子被突破”  
而不是  “无法满足拉格朗日乘子约束，因此优化报错”
...

不是报错，是边厚出现负数类似情况

隔壁的薛定谔

dpstt 发表于 2024-3-18 11:31
不是报错，是边厚出现负数类似情况

看了下MNE，MNE是根据近轴光线追迹 从而确定的，因此可能与实际情况不一样，所以不准确？
所以“边厚出现负数”  可能并没有突破拉格朗提乘子约束，只是MNE计算不准确，系统误认为“边厚出现负数”是满足MNE拉格朗日乘子约束的？

dpstt

隔壁的薛定谔 发表于 2024-3-18 11:40
看了下MNE，MNE是根据近轴光线追迹 从而确定的，因此可能与实际情况不一样，所以不准确？
所以“边厚出 ...

照这么说，MXT也是同样的情况，MXT 6可能约束不了