数学规划法

大家好
29
大家好
30
3.4 可行方向法
I. 概述
结构优化的一般数学规划表达式：
寻找一组设计变量 X = ( X1 , X2, ……, Xn )T min f( X ) X E n
s. t. g i (X) 0 i =1, m 设计变量的迭代公式 ----- X ( +1) = X ( ) + P ( ) 从 X ( ) 调参至 X ( +1) , 要求设计点可行, 并且目标函数还要下降, 即满足可用可行性条件:
= (1/2) (+1) ’0
大家好
33
2. 从约束边界点开始, 沿目标函数等值面(线)内向可行域 中侧移, 使设计点离开约束边界进入可行域.
大家好
34
两种调参方向交替进行，直到最优点为止。
这是一种较早期的方法，迭代次数多，代价大，步长 试凑，有太多的随意性。
大家好
35
III. 梯度投影法 Gradient projection method
Rosen 的梯度投影法适用于目标函数为非线性、约束 条件为线性的结构优化设计。
一般结构优化设计问题，多取重量为目标函数，它是 设计变量的线性函数，约束条件是应力、位移、稳定， 均为设计变量的非线性函数。
如果引入倒数设计变量，并把约束条件经过显式线性 近似处理，则问题正好符合梯度投影法要求。
大家好
P ( ) = - f ( X ( ) ) or P ( ) = f ( X ( ) )
X ( +1) = X ( ) - f ( X ( ) ) --- 最速下降 X ( +1) = X ( ) + f ( X ( ) ) --- 最速上升 其中步长按 采用试凑法. = 2 0 ( 在可行域内)
线性规划 非线性规划 如果目标函数与约束函数都是凸函数，则称为凸规划 如果目标函数是二次函数而约束函数是一次函数，则称
为二次规划 如果设计变量只允许取整数，则称为整数规划 如果在目标函数和约束函数中包含具有随机性质的参数
则称为随机规划
大家好
5
III. 数学规划问题的求解
对于线性规划问题，单纯形法十分有效
原问题：
min f ( x ) x = { x1, x2, x3,……, xn}T (1) 求解其最优化的必要条件：
▽f(x*)=0
(2)
但是式(2)是一个非线性方程组，与求解原问题同样困难。
在数学规划法中，是用迭代下降的算法找到极小值点。 即先假定一个初始设计x(0),然后在第k次迭代 (k=0,1,2,…),用x(k+1)代替x(k)，要求x(k+1)比x(k)更接近最 优解。对于无约束优化问题，也就是要求目标函数有 所下降，即
最速下降法是一阶线性收敛，收敛速度较慢。
最速下降法与变量长度有关，即与变量尺度关系很 大。
1) 最速下降法迭代过程简单，不受初始点位置限制。 因此虽然该方法有收敛慢，依赖于变量的尺度等缺 点，但这些缺点往往在最优点附近表现得才明显。
大家好
24
大家好
25
ii. Newton法
大家好
26
从Newton法迭代公式的推导过程可知，对任何二次 函数，用Newton法求极小点，只需迭代一次（如果 该二次函数有极小点存在的话）
量ε 前后两次迭代目标函数值下降的相对值已经足够小 工程结构优化时常采用固定的迭代次数作为停止迭代
准数学规划方法的一种基本 方法，也是其它优化方法的一种中间手 段
大家好
14
i. 0.618法
0.618法的基本思想 是通过取试探点和进 行函数值比较，使包 含极小点的搜索区间 不断缩短，从而各点 可以看作为极小点的 近似。这类方法仅需 计算函数值，用途广 泛，尤其适用于非光 滑函数形式。
结构优化设计
大家好
第三章 数学规划法
大家好
2
3.1 数学规划问题的分类及解法
I. 数学规划问题的一般提法是：
寻找一组设计变量变 X = { X1, X2, X3,……, Xn}T 使得 f ( x ) min
s. t. g i ( X ) 0 ge(X)= 0
i = 1,2, ……, m e = 1, 2, ……, n
36
mifn(X)(14)2(13)2
X1
X2
s.t. 1 6 12 X1 X2
34 12
X1 X2
X
1
1 3
m f( i z )n (z 1 4 )2 (z 2 3 )2
二点二次插值法I
大家好
19
大家好
20
二点二次插值法II（割线法）
大家好
21
不同搜索方法比较
连续性 收敛性
0.618法 0阶
牛顿法
2阶
二点二次插值法 1阶
割线法
1阶
0.618 2 1.618 1.618
大家好
适用范围
最优点在所选 区间中 靠近最优点
靠近最优点
靠近最优点
22
3.3 确定搜索方向的算法
由此构成三种不同形式的可行方向 法。
大家好
32
II. 最佳矢量法 Optimum vector method
交替使用两种行进方向:
1. 当设计点不在约束边界上时, 采用最速下降法 ( Steepest descent method )或最速上升法 ( Steepest ascent method ) .
f(x (k)+ P(k) )- f(x(k))≈ ▽Tf(x(k)) P(k) <0 由于是正标量，所以：
▽Tf(x(k)) P(k) <0 或者 -▽Tf(x(k)) P(k) >0 这说明搜索方向应该和目标函数负梯度方向夹角小于90， 这样的方向称之为下山方向。
大家好
8
基本的下降算法：
令k=0,给定初始解x(0)； *求搜索方向P(k),使▽Tf(x(k)) P(k) <0； *求搜索步长(k)，要求f(x (k)+ (k)P(k) )=min f(x (k)+
f(x(k+1))< f(x(k))
大家好
7
在数学规划中，一般迭代格式可以写成：
x (k+1) = x (k)+ P(k) , --- 步长( Step length )，正标量 P(k) --- 方向向量( Directional vector )
因此目标函数的下降要求可以改写成：
f(x (k)+ P(k) )< f(x(k)) 如果函数f(x)在x(k)是一次可微的，则对足够下小的有：
因此，如果定义一个一元函数φ():
φ()=f(x (k)+ P(k) ) 决定(k)方法就是估计φ的极小值，即φ()至少近似满足：
φ’()=0 这个方程是一个非线性方程，要用到一维搜索方法才能
确定步长，因此一维搜索对于提高下降算法非常重 要。
大家好
10
一维搜索算法
0.618法 Newton法 割线法 抛物线法 三次插值法
P(k) ) 修改x (k+1) = x(k)+ (k)P(k) 1) 检查收敛原则，不满足时令k=k+1,返回2)；满足则
停机。
大家好
9
一维搜索确定步长
步长(k)的决定方式常常是使目标函数在x (k+1) 点上达到 沿P(k)方向最小，即要求：
f(x (k )(k )P (k )) m 0f(x i(k ) n P (k ))
对Newton法H(k)=[▽2f(x(k))]-1 ，P(k)取负梯度， (k)=1
最速下降法收敛太慢， Newton法收敛最快，但要计算二 阶导数并要求求逆，工作量大，有时还会碰到不可克 服的困难。
因此人们想若H(k)的选取不需要计算二阶导数矩阵和求逆， 而又能逼近[▽2f(x(k))]-1 ，那么所确定的算法可能会类 似于Newton法收敛很快。基于这种想法，发展了一种 拟Newton法。 H(k)构造方法不同，有不同的拟Newton 法，其中DFP算法就是这类算法中最为著名的。
i. 最速下降法
前面已经知道，目标函数沿负梯度方向下降最快，因此 取负梯度方向为搜索方向，即：
P(k) = -▽f(x(k))
x (k+1) ＝x (k)- ▽f(x(k))
基本算法：
给定初始点x(0)，令k=0；
计算▽f(x(k)) 若▽f(x(k)) ≤ ε成立，则x*=x(k)，停机， 否则转下一步；
初始点要求较高，一定要在很靠近最优点附近选取
最优点。同时Newton法要求二阶导数矩阵，工作量 较大，且要求该矩阵的逆，这就要求▽2f(x(k))是非奇 异的
大家好
27
iii. 变尺度法
无约束优化的迭代公式为：x (k+1) = x(k)+ (k)H(k)P(k)
对最速下降法H(k)=I，P(k)取负梯度， (k)用一维搜索
其中, X ----设计变量 f ( x ) ----目标函数 g i ( X ) 和 g e ( X ) ---- 约束条件
大家好
3
II. 数学规划问题的分类
(1) 按约束的有无，可分为： 无约束最优化问题 有约束最优化问题 准无约束最优化问题
大家好
4
(2) 按目标函数和约束函数是否为线性，可分为：
1. 满足可用性条件 ( Usability condition )
f ( X ( ) )T p ( ) 0
2. 满足可行性条件 ( Feasibility condition )
g i (X ( ) + P ( ) ) 0
或者 g i (X ( ) )T P ( ) 0