机械优化设计考试复习
1
2 3外推法确定搜索区间,函数值形成 高-低-高 区间 4数学规划法的迭代公式是 X计算最佳步长
5若n维空间中有两个非零向量d0,d1,满足(d0)TGd1=0,则d0、d1之间存在_共轭关系 6,与负梯度成锐角的方向为函数值 方向,与梯度成直角的方向为函数值 不变 方向。 外点;内点的判别
7 8、那种方法不需要要求一阶或二阶导数: 坐标轮换法
9、拉格朗日乘子法是 升维法 P37
10 11,.函数fx1,x2x1x2
2
2
k1
X
k
kd和
k
212
4x1x25在X0点处的梯度为,海赛矩阵为
40
2
44 2
12.目标函数是一项设计所追求的指标的数学反映,因此对它最基本的要求是能用 来评价设计的优劣,,同时必须是设计变量的可计算函数 。
13.建立优化设计数学模型的基本原则是确切反映。
14.约束条件的尺度变换常称 规格化,这是为改善数学模型性态常用的一种方法。 15,.随机方向法所用的步长一般按此法是指依次迭代的步长按一定的比例 递增的方法。
16.最速下降法以 负梯度 方向作为搜索方向,因此最速下降法又称为 梯度法,其收
敛速度较。
17,.二元函数在某点处取得极值的充分条件是fX00正定
18.拉格朗日乘子法的基本思想是通过增加变量将等式约束 优化问题变成 无 约束优化问题,这种方法又被称为 升维 法。
19
20坐标轮换法的基本思想是把多变量 的优化问题转化为 单变量 的优化问题
21.在选择约束条件时应特别注意避免出现 相互矛盾的约束, ,另外应当尽量减少不必要的约束 。
22.目标函数是n维变量的函数,它的函数图像只能在空间中描述出来,为了在n维空间中反映目标函数的变化情况,常采用 目标函数等值面 的方法。 23协调曲线法是用来解决 设计目标互相矛盾 的多目标优化设计问题的。
24.是首要和关键的一步,它是取得正确结果的前提。 二、解答题
1、试述两种一维搜索方法的原理,它们之间有何区别 答:搜索的原理是:区间消去法原理
区别:(1)、试探法:给定的规定来确定插入点的位置,此点的位置确定仅仅按照区间的缩短如何加快,而不顾及函数值的分布关系,如黄金分割法
(2)、插值法:没有函数表达式,可以根据这些点处的函数值,利用插值方法建立函数的某种近似表达式,近而求出函数的极小点,并用它作为原来函数的近似值。这种方法称为插值法,又叫函数逼近法。
2、在变尺度法中,为使变尺度矩阵Hk与Gk近似,并具有容易计算的特点,Hk必须附加哪些条件?
1
答:(1)必须是对称正定的(2)要求有简单的迭代形式 (3)必须满足拟牛顿条件 3,总结:无约束优化方法 只算函数值方法
1,坐标轮换法:小规模,收敛慢(无耦合问题快);2,单形替换法:中小规模,收敛较快, 3,格点法:非凸问题;4,Monte Carlo 法:非凸问题。 计算一阶导数方法
1, 梯度法:中小规模,开始快;2,共轭梯度法:中大规模,收敛快,程序简单; 2, 变尺度法:中大规模,收敛快;4,Powell方法:中大规模,收敛快。 计算二阶导数方法
1, Newton 方法:收敛快,计算难度大;2,共轭方向法:收敛快,计算难度大。 4.共轭梯度法中,共轭方向和梯度之间的关系是怎样的?试画图说明。 . 对于二次函数,fX
12
XGXbXc,从
k
TT
k
沿G的某一共轭方向d作一维搜索,到达X点出发,
X
k1
点,则Xk1点处的搜索方向d
T
j
应满足
dg
j
k1
gk0,即终点X
k1
与始点Xk的梯度
kj
之差gk1gk与d的共轭方向d正交。
3.为什么说共轭梯度法实质上是对最速下降法进行的一种改进?.
答:共轭梯度法是共轭方向法中的一种,在该方法中每一个共轭向量都依赖于迭代点处的负梯度构造出来
的。共轭梯度法的第一个搜索方向取负梯度方向,这是最速下降法。其余各步的搜索方向是将负梯度偏转一个角度,也就是对负梯度进行修正。所以共轭梯度法的实质是对最速下降法的一种改进。
4.简述随机方向法的基本思路
答:随机方向法的基本思路是在可行域内选择一个初始点,利用随机数的概率特性,产生若干个随机方向,并从中选择一个能使目标函数值下降最快的随机方向作为可行搜索方向。从初始点出发,沿搜索方向以一定的步长进行搜索,得到新的X值,新点应该满足一定的条件,至此完成第一次迭代。然后将起始点移至X,重复以上过程,经过若干次迭代计算后,最终取得约束最优解。
5.凸规划:对于约束优化问题 minfX s.t. gjX 若f
0 (j1,2,3,,m)
X、gjX
(j1,2,3,,m)都为凸函数,则称此问题为凸规划。
6目标函数值下降,且不会越出可行域。 7.设计空间:n个设计变量为坐标所组成的实空间,它是所有设计方案的组合 8.收敛性:是指某种迭代程序产生的序列X
k
k
0,1,收敛于limX
k
k1
X
9. 黄金分割法:是指将一线段分成两段的方法,使整段长与较长段的长度比值等于较长段与较短段长度的比值。
10.可行域:满足所有约束条件的设计点,它在设计空间中的活动范围称作可行域。 三,计算
1、求目标函数f(X)=2x12+3x22-x1x2-2x2-9在点X1 =1,1处的函数变化率最大的方向及其数值。 f(x)x14x1-x241解:▽f(x1)=f(x)=
6x2-x1-261x2
p
f(x1)
f(x1)
3 23
3,3
99
数值f(x1)2
23^23^232
2、求函数f(X)=x13+x22-4x1-2x2+27在点X1(2,1)处的二阶泰勒展开式。
Tx-x1f(x)f(x1)f(x1)
f(x1)27-1
3x1^2-48
f(x1)
2x2-206x10120
H(x1)
0202x1-2
f(x)27-180
x2-1
22
=6x1+x2-16x1+6+27
x-xTH(x)x-x
1
1
1
解:
x
1
120x1-2-2x2-1
02x2-1
3、用共轭梯度法求函数f(X)=2x1^2-x1x2+3x2^2+5的最优解,初始点(1,2),迭代精度0.02
初始点x
(0)
1,2,s
T
(0)
4x1x22
f(x0)
x16x211
T
x
(1)
x
(1)
(0)
(0)
s
(0)
1,2
(0)
(0)
12(0)x1(1)2
(0)(1)
11211x2
(0)
将x代入f(x)得f(
(0)
)2(12)(-12
2(0)
)(211
(0)
)
(3211令
df(
(0)
)5
(0)
2
)
(1)
(0)
0得5
xs
(1)
12(0)11
(0)
112211
(1)
2
2
f(x
(1)
f(x
(1)
)
f(x
(0)
(0)
s
(0)
解:
(0)
)
)
88309
f(x
4x1x268)
x16x2661
s
(1)
682176550
88309
-66111972060x
(1)
x
(2)
s
(1)
s
(1)
x
(1)
f(x
(1)
)
T
(1)
s
(1)T
H(x
(1)
)s
(1)
0.0001
(2)
111765507.650.000111297206014.80
7.65
x
14.80
最优解为:
2
4,求函数 , x 2 ) x 1 x 22 4 x 1 2 x 2 5 的极值。 f ( x1
解 首先,根据极值的必要条件求驻点 f
x2x14 001
f(x) f02x22x
x2x
x1020得驻点为 x
x201
再根据极值的充分条件,判断此点是否为极值点。由于 2
f
20的一阶主子式和二阶主子式分别为 2
x1x
2 2ff
2
2020x1x1x200H(x)H(x)40 22
f02f02
2xxx2x
21
故 H 0 ) 为正定矩阵 x 0 , 2 T 为极小点,相应的极值为 f ( x 0) (x 0 1
5.试用牛顿法求fX8x15x2的最优解,设X
2
2
0
1010。
T
初始点为X
f
0
1010,则初始点处的函数值和梯度分别为
T
X1700
f
X
16x14x2200,沿梯度方向进行一维搜索,有 4x10x14012
XX0f
10
X
10
0
10200102000
101401400
0为一维搜索最佳步长,应满足极值必要条件
fX
minfX
1
fX
2
2
min810200041020001014005101400
min
010600000596000, 从而算出一维搜索最佳步长 0
[1**********]0
1
0.0562264
1020001.2452830
则第一次迭代设计点位置和函数值X
101402.12830190
f
X24.4528302,从而完成第一次迭代。按上面的过程依次进行下去,便可求得最
1
优解。
6、试用黄金分割法求函数f
20
的极小点和极小值,设搜索区间
a,b0.2,1(迭代一次即可)
解:显然此时,搜索区间a,b0.2,1,首先插入两点1和2,由式
1b(ba)10.61810.20.5056 2a(ba)0.20.61810.20.6944
计算相应插入点的函数值f140.0626,f229.4962。 因为f1f2。所以消去区间a,1,得到新的搜索区间1,b, 即1,ba,b0.5056,1。
第一次迭代:
插入点10.6944, 20.50560.618(10.5056)0.8111
相应插入点的函数值f129.4962,f225.4690,
由于f1f2,故消去所以消去区间a,1,得到新的搜索区间1,b,则形成新的搜索区间1,ba,b0.6944,1。至此完成第一次迭代,继续重复迭代过程,最终可得到极小点。
7.用牛顿法求目标函数fX16x2
2
125x2+5的极小点,设X
0
22T
。
f
解:由 X
0
22T
,则f
X0
x
1
32x1f
6450x100 x22
2f2
f
2
f
X
x2
1
x
1x20
32
2
f2
f0
50,其逆矩阵
x2x1x2
2
1
0
21
fX0
3201 50
1
0
因此可得:X1X02
fX
1
f
X0
232216401000
50
fX15,从而经过一次迭代即求得极小点X0
0T
,f
X
5
为
§3.6 无约束优化设计方法小结(续)
Poweel法
迭代次数搜索次数
共轭梯度法牛顿法变尺度(DEF)法
2 2 1 26 2 1 2
搜索方向收敛速度
共轭方向一阶算法较慢
f(x
(k1)
)
(k)
S
(k)
H(x
(k)
)
1
f(x
(k)
)H
(k)
f(x
(k)
)
超线性二次收敛
中
x,f(x),f(x)
二阶算法收敛最快最大
x,f(x),f(x)H(x),H(x)
1
超线性二次收敛较大
x,f(x),f(x)E
(x)
存储量小
x,f(x)
,H
(x)
适用维数稳定性
n
好中差中下
黄金分割法 程序图
共轭梯度法 程序图
变尺度法程序图