内点惩罚函数法的改进算法.pdfVIP

第 l7卷第 3期 工 科 数 学 V0】．17_No．3 200]年 6月 JOURNAl，OFMATHEMATICSFORTECHNOLOGY Jun．2001 内点惩罚函数法的改进算法 邵陆寿， 吴怀孔 (安徽农业太学 ．台肥 230036) [摘 要]本文针对 内点惩罚函数法现有算法在计算中设计变量常常越出可行域边界，导致计算失效而 提 出了一种改进 的计算方法 该方法 的主导思想就是保证对 内点惩罚函数的求极值过程一直限定在可行设 计区域 内一从而保证各 ‘(7“，均在可行域内 并进行 了实例计算验证 ． [关链词]内点惩罚函数；计算方法；可行设计区域 [中圈分类号]0224 ：文献标识码：A [文章编号 1007—4120(2001)03．0026—05 1 前 言 惩罚函数法(SUMT法)中的内点法，具有解决机械约束优化设计问题的一个诱人的特点，它能计 算出一系列的逐步得到改进的可行的设计方案，因而为设计人员提供了根据具体的实际情况选择最恰 当方案的途径．但是根据笔者多年的教学实践及科研实践发现，该方法的现有算法 在计算中经常出现 失效 ，其表现形式为设计变量常常越出可行设计域 ，那么此时既使得到了一个最优解。此解也是一个 不可行方案． 2 现有算法设计变量越出可行域的原因 内点惩罚函法适于处理的数学模型为 min ，( )， s．1． ( )≤0． 内点惩罚函数可表示为 )=，(_) ， 其中 y。 y y ，… ， lira7 =0． (2) 现有算法t直接将约束优化模型转化为无约束优化模 型，即求 min~(R，)．求得 ()，然后以 ()为初点， 使 y 。=d (1)，再无约束寻优 ，求 ( )．…，直至 满足收敛精度为止．其求解过程如图1所示． 由图】可知，构造的惩罚函数 ≠( ，)若是一个凸函 数t即使设计变量 ，在迭代 中有时会越出可行域边界 ，但 与 ’()∈ ，所以该方法是正确的．但是若 ≠(，)不是 一 个凸函数 ，如图2所示． 圈 1 内点涪迭代示意围 [收稿 日期]2000．O2—28 第3期 部陆寿等 ：内点惩罚函数法的改进算法 27 此时现有算法在调用无约束优化方法求惩罚函数极值时，设计变量就有可能越出可行域边界．由图 2分析，只有保证搜索 区间[n， ：∈ ，如图 2 中， 一A，口3一B所示 ，方能求得 。(y。)，并且 。(y)∈ ．依次类推 ，只有保证[ ， ]∈寥 ，则 。(y)∈ ．但仔细分析一下一维搜索方法，对 非凸的函数 ，这一点根本无法保证．根据下式： X 一xf一．+n ， 将其代人惩罚函数构造公式 (2)，得 (x ，y)—， () (一1，2，… ，N)． 各 (n)的性状并非一致 ，且并非一定是 凸函数． 因此并不能保证各 ()的极值点都落在 域 中．如图 3所示． 图 3中，Y (。)能找到区间[．4，B]，且能找到 n ∈ ．但对 ^()若找到 的区间为[J4，B ，且 只能找到 告毋 ，设计变量就越 出了可行域边 界． 对 Ⅳ维非凸函数 ，k轮搜索的最终点 就不 图2 非凸惩罚函数选代示意田 能保证 ∈够．如 图2，有可能我到 而不是 ’而机械优化设计中．≠( ，)属于凸函数是有 限．因 此在机械优化设计中，使用 内点法常常失效． 圈3 一维区间搜索示意图 3 内点惩罚函数的改进算法 针对模型(1)及内点惩罚函数构造式(2)可转化为如下的模型’ ：