外点罚函数优化实例.pdf

实用标准文案 外点罚函数优化实例 一、优化问题 如图 1 所示，为某一桁架的一部分，杆 2 距 O 点 30cm 处有一支点 C。为了 固定桁架，现欲在杆 1 和 2 上设置支点 A 和 B，用来连接杆 3 （可拆卸）。已知当 桁架固定时，杆 1 和 2 成直角；而且，杆 1 右边有一段长为 20cm 的重要部位， 不能设置支点。卸去杆 3 、收起桁架时，支点 A 的位置不能高于 BC 段中点 D 。求 取支点 A 、B 的位置，使得杆 3 的长度尽量小，以节省材料。 图 1 桁架结构示意图 二、数学模型 设 A 、B 两点距离 O 点的长度分别为 x 和 x ，而桁架固定时杆 1 和 2 成直角。 1 2 所以杆 3 的长度为 l x 2 x2 。 3 1 2 1 由图 1 可知， x1 20 且 x1 (x 2 30) ，即 x1 20 0 且 2x1 x2 30 0 。 2 设 X [x , x ]T ，取 f (X ) l 2 x 2 x 2 。因此，数学模型为： 1 2 3 1 2 极小化目标函数 2 2 min f (X ) x1 x2 约束条件 精彩文档 实用标准文案 g (X ) x 20 0 1 1 g2 (X ) 2x1 x2 30 0 三、求解数学模型 (1) 外点罚函数法求解 构造外点法罚函数，如下： ( k ) 2 2 (k ) 2 2 (X , M ) x x M {[(m ax( 20 x ,0)] [max( 2 x x 30,0)] } 1 2 1 1 2 程序流程图如图 2 所示： 精彩文档 实用标准文案 给定 X (0 ) 、 M (0 ) 、c 、 k=0 i=0 k=k+1 求 (X ( i) ) 与 i=i+1 Hessian 矩阵