A、其缩短率为0.618
B、a=b-入(b-a)
C、a=a+入(b-a)
D、在该方法中缩短搜索区间采用的是外推法。
11、与梯度成锐角的方向为函数值A方向,与负梯度成锐角的方向为函数值_B方向,与梯度成直角的方向为函数值C方向。
A、上升
B、下降
C、不变
D、为零
12、二维目标函数的无约束极小点就是B。
A、等值线族的一个共同中心
B、梯度为0的点
C、全局最优解
D、海塞矩阵正定的点
13、最速下降法相邻两搜索方向dk和dk+1必为B向量。
A相切
B正交
C成锐角
D共轭
14、下列关于内点惩罚函数法的叙述,错误的是A。
A可用来求解含不等式约束和等式约束的最优化问题。
B惩罚因子是不断递减的正值
C初始点应选择一个离约束边界较远的点。
D初始点必须在可行域内
三、问答题(看讲义)
1、试述两种一维搜索方法的原理,它们之间有何区别?
2、惩罚函数法求解约束优化问题的基本原理是什么?
3、试述数值解法求最佳步长因子的基本思路。
4、试述求解无约束优化问题的最速下降法与牛顿型方法的优缺点。
5、写出用数学规划法求解优化设计问题的数值迭代公式,并说明公式中各变量的意义,并说明迭代公式的意义。
6、什么是共轭方向?
满足什么关系?
共轭与正交是什么关系?
四、解答题
1、试用梯度法求目标函数f(X)=1.5xi2+0.5x22-xix2-2xi的最优解,设初始点x(0)=[-2,4]T,选代精度&=0.02(迭代一步)。
解:
首先计算目标函数的梯度函数(7],
Lx2—xl
计算当前迭代点的梯度向量值
梯度法的搜索方向为;疋=&」n,因此在迭代点x(0)的搜索方向为[12,-6]t
在此方向上新的迭代点为:
當沬皿仝3+於M=兀十昭匚喚
=
把新的迭代点带入目标函数,目标函数将成为一个关于单变量的函数)瑚躅
)=f=1.5C-2-12ct):
_0,5(4—6⑴:
-(-2_1'2^)(4-
—2(—2+=F(r)
-'■■",可以求出当前搜索方向上的最优步长
新的迭代点为
[1.5Z921
L2.23541
当前梯度向量的长度一-.一丨一.-'i44■■,因此继续进行迭代。
第一迭代步完成
2、试用牛顿法求f(X)=(xi-2)2+(xi-2x2)2的最优解,设初始点x(0)=[2,1]T。
解1:
(注:
题目出题不当,初始点已经是最优点,解2是修改题目后解法。
)牛顿法的搜索方向为祖淙「炉價厂也毬I,因此首先求出当前迭代点x(0)
的梯度向量、海色矩阵及其逆矩阵
0=沪矿如(°=盼
不用搜索,当前点就是最优点。
解2:
上述解法不是典型的牛顿方法,原因在于题目的初始点选择不当。
以下修改求解题目的初始点,以体现牛顿方法的典型步骤。
以非最优点x(0)=[1,2]T作为初始点,重新采用牛顿法计算
牛顿法的搜索方向为护—班时F,因此首先求出当前迭代点x(0)
的梯度向量、以及海色矩阵及其逆矩阵
梯度函数:
初始点梯度向量:
海色矩阵:
海色矩阵逆矩阵:
.Ja
当前步的搜索方向为:
骷眼j.护闵沁罚仁川詐[;]
新的迭代点位于当前的搜索方向上:
塞很HO二家3十沛:
二峯5十冷简
把新的迭代点带入目标函数,目标函数将成为一个关于单变量的函数阴倒
—-■=■,可以求出当前搜索方向上的最优步长
新的迭代点为噩如=工-—「_=0-171=13
当前梯度向量的长度—:
;—二-^二:
因此继续进行迭代。
第二迭代步:
H»zl——41
I-4*xlJ
F(f(xw)〕=:
IVffl=0因此不用继续计算,第一步迭代已经到达最优点。
这正是牛顿法的二次收敛性。
对正定二次函数,牛顿法一步即可求出最优点。
3、设有函数f(X)=xi2+2x22-2xix2-4xi,试利用极值条件求其极值点和极值。
解:
首先利用极值必要条件
D找出可能的极值点:
[2—2«—I4*x2—2*xl
求得辭卜IS,是可能的极值点。
再利用充分条件爭气疗正定(或负定)确认极值点
|2|=2>0
t71=8-4=4>°
因此正定,是极小点,极值为f(X*)=-8
4、求目标函数f(X)=xi2+xix2+2x22+4xi+6x2+10的极值和极值点。
解法同上
5、试证明函数f(X)=2xi2+5x22+x32+2x3X2+2x3xi-6x2+3在点[1,1,-2]t处具有极小值。
解:
必要条件:
+2*x3
VCO=H-2x3-*
2+2*x2+2*x3
将点[1,1,-2]T带入上式,可得
(J
vCO=o
0-
充分条件
縊Q21
沪00=0102
1.222)
4
0
0
10
=40
4D2
0102=80—40-16=24->0
'萨:
询正定。
因此函数在点[1,1,-2]t处具有极小值6、给定约束优化问题
minf(X)=(x1-3)2+(x2-2)2
s.t.g1(X)=—X12—X22+5>0
g2(X)=—X1—2X2+4>0
g3(X)=x1>0
g4(X)=x2>0
验证在点X[2,i]TKuhn-Tucker条件成立。
解:
首先,找出在点X[2,i]T起作用约束:
gi(x)二0
g2(x)二0
g3(x)二2
g4(x)二1
因此起作用约束为gi(X)、g2(X)。
然后,计算目标函数、起作用约束函数的梯度,检查目标函数梯度是否可以表示为起作用约束函数梯度的非负线性组合。
T[二洽匚1‘二1
求解线性组合系数临3:
:
沁慮⑪沁%碗
得到.—-=-均大于0
因此在点X[2,i]TKuhn-Tucker条件成立
7、设非线性规划问题
min
f(X)
(X1
2)2
2
X2
s.t.
g(x)
X
0
g2(X)
X2
0
ga(X)
2
X1
xf1
0
*T
用K-T条件验证X1,0为其约束最优点
解法同上
&已知目标函数为f(X)=x1+X2,受约束于:
gi(X)=-xi2+X2>0
g2(X)=xi为
写出内点罚函数。
解:
内点罚函数的一般公式为
罰心宀8F5+N习说
r(k)=cr(k-1),0vCv1
4
因此罚函数为:
9、已知目标函数为f(X)=(xi-1)2+(x2+2)2受约束于:
gi(X)=-x2-xi-1>0g2(X)=2-xi-x2>0g3(X)=xi>0g4(X)=x2>0试写出内点罚函数。
解法同上
i0、如图,有一块边长为6m的正方形铝板,四角截去相等的边长为x的方块并折转,造一个无盖的箱子,问如何截法(x取何值)才能获得最大容器的箱子。
试写出这一优化问题的数学模型以及用MATLAB软件求解的程序。
料最少,试写出这一优化问题的数学模型以及用MATLAB软件求解的程序
12、一根长I的铅丝截成两段,一段弯成圆圈,另一段弯折成方形,问应以怎样的比例截断铅丝,才能使圆和方形的面积之和为最大,试写出这一优化设计问题的数学模型以及用MATLAB软件求解的程序。
13、求表面积为300m2的体积最大的圆柱体体积。
试写出这一优化设计问题的数学模型以及用MATLAB软件求解的程序。
14、薄铁板宽20cm,折成梯形槽,求梯形侧边多长及底角多大,才会使槽的断面积最大。
写出这一优化设计问题的数学模型,并用matlab软件的优化工具箱求解(写出M文件和求解命令)。
15、已知梯形截面管道的参数是:
底边长度为c,高度为h,面积A=64516mm2,斜边与底边的夹角为见图1。
管道内液体的流速与管道截面的周长s的倒数成比例关系
(s只包括底边和两侧边,不计顶边)。
试按照使液体流速最大确定该管道的参数。
写出这一优化设计问题的数学模型。
并用matlab软件的优化工具箱求解(写出M文件和求解命令)。
16、某电线电缆车间生产力缆和话缆两种产品。
力缆每米需用材料9kg,3个工时,消耗电能4kW・h,可得利润60元;话缆每米需用材料4kg,10个工时,消耗电能5kW・h,可得利润120元。
若每天材料可供应360kg,有300个工时消耗电能200kWh可利用。
如要获得最大利润,每天应生产力缆、话缆各多少米?
写出该优化问题的数学模型以及
用MATLAB软件求解的程序。