7优化设计-4多维优化之约束优化方法.pptVIP

约束优化方法;一、约束优化问题的一般形式 ;二、约束优化方法分类　　;惩罚函数法　　; 　基本思想： 依据原约束优化问题的约束条件构建可限制其目标函数值脱离可行域之外的约束函数，并将其与原目标函数共同组成一个新目标函数，进而通过对新目标函数的求解实现约束优化问题向无约束优化问题的转化和求解．;2、惩罚函数法的内涵和本质;３、惩罚函数法原理　;　　　　惩罚函数构成　;　　惩罚函数应具备的特性　;　4、惩罚函数法的分类　　　　;内点法（障碍函数法）;或：;　3）内点法惩罚因子rk rk是一随迭带次数依序递减的正数序列;4）内点法惩罚函数性质 １）惩罚函数可行域与目标函数相同 ２）点越靠近目标函数约束边界所受惩罚越严重 ３）惩罚函数因惩罚因子不同而不同，惩罚因子越小，惩罚函数求解精度越高 ４）伴随迭代次数增加，惩罚因子不断减小，当惩罚因子趋于０时，惩罚函数的极小点就是目标函数的最优点;5）内点法收敛条件 　内点法收敛性的判断相对于无约束优化方法稍微复杂一些，具体就是要从两迭代点之间的距离以及两迭代点所对应的目标函数值的差距入手同时进行判断。;6）内点法计算步骤 １）选取初始点 x0 及惩罚因子初值r0、降低系数c、计算精度 ε1、ε2 ，令 k=0； ２） 构造惩罚函数； ３）选用一种无约束优化方法，求惩罚函数的最优解 xk*和对应目标函数值 ４）进行收敛判断，满足收敛条件，计算结束，xk* 为最优点，否则令 　　　X0=Xk*；rk+1=C*rk；k=k+1 　返回步骤3 进行下一轮迭代计算。;7）内点法使用的注意事项 1、初始点在原目标函数可行域内选取且不要离约束边界太近。可依据设计变量取值区间的上下限,按如下方法进行选取 x0=a+r(b-a) ；r为[0，1]区间的随机数 2、惩罚因子初值不要太小，否则影响精度 3、递减系数C应该在0.02-0.1之间取值;特点：　 　１、可用于目标函数比较复杂，或在可行域外无定义的场合下 　２、优化始终在可行域内进行，至少可以得到满足工程要求的近似解 不足： 　１、初始点和中间迭代的极值点均需严格满足所有约束条件 　２、不能解决等式约束问???　;外　点　法;２） 外点法惩罚项及惩罚函数的形式;　可直接作为惩罚项;　　注意！外点法的惩罚因子rk是一个伴随迭代依次递增的正数序列．;4）外点法收敛条件;5) 外点法惩罚函数性质;性质３：当迭代次数足够大，惩罚函数中各项违反约束的函数取值趋于０,惩罚函数的极小点就是目标函数的最优点;6）外点法计算步骤 １：给定初始点x0 以及初始惩罚因子 r0、递增系数a、收敛精度 ε1 ε2 ，令 k=0； ２： 构造惩罚函数； ３：用无约束优化方法求惩罚函数的最优解 xk* 和对应函数值 ４：运用终止准则进行收敛判断，满足收敛条件，计算结束，xk* 为最优点，否则 令X0=xk*；rk+1=a*rk*；k=k+1,返回步骤3 继续计算;7）外点法注意事项;　　２、约束裕量 　　由于外点法迭代在目标可行域之外进行，所以迭代最优点依然是不可行解，特别是对于仅仅由不等约束构成的惩罚函数，所以应该加强原约束条件所对应约束函数的约束强度，即留出约束裕量　　　　;8）外点法的优势和不足;混合法;三. 特点和注意事项 　１、可求解包含等式和不等约束的优化问题 　２、初始点便于选择，初始点可以不在可行域内，只要满足所有不等约束即可 　３、罚因子递减系数按照内点法选取： 　４、计算步骤和流程与外点法相同;惩罚函数优化设计实例; 已知条件： 受力P=294300N 跨度B=1520mm 壁厚T=2.5mm； 弹性模量E=2.119*105Mpa 比重ρ=8.142*10-5N/mm3； 许用压应力[σ]=680Mpa； ;设计变量：钢管直径D、长L或人字架高h; 该人字架设计问题可归结为如下的优化问题数学模型：; 内点法求解： １、建立约束函数：X=[D,L]　;２０ ≤ D≤ 140