《最优化方法》(研究生)期末考试练习题答案
二.简答题
min -5y19y2, s.t. 4y13y2 3, -2y1 y22, 1.
3y14y28, y1,y20;
x3 x40, (以x1为源行生成的割平面方程) 2.
注意:在x1为整数的情况下,因为x3,x40,该方程自然满足,这是割平面的退化情形
1
656
111
x3 x4, (以x2为源行生成的割平面方程)
442
3.
a10,b13
1a10.382(b1a1)00.382*31.146
1a10.618(b1a1)00.618*31.854 (1)(1.146)32*1.14610.2131(1)(1.854)32*1.85413.6648
事实上,不经计算也可以看出
(1)(1),所以a20,b21.854。
即:初始的保留区间为[0,1.854]。近似的最优解:x*
01.854
0.927.2
f1(x)x1ex2*(1)2.7x1ex22.7
4.令
f2(x)x1ex2*(0)1x11f1(x)x1e
x2*(1)
0.4x1e
x2
0.4
f1(x)x1ex2*(2)0.1x1e2x20.1
拟合问题等价于求解下列最小二乘问题:min
(f(x))
i
i1
4
2
三.计算题
1.分别用最速下降方法和修正的牛顿法求解无约束问题 minf(x)x14x2。 取初始点x
2
2
1
T
2,2,0.1.
解:f2x1
在初始点x1
2,2T
8x2,
,
f4
16
从而最速下降法的搜索方向为:d14
16
.
最优步长*满足f(x(1)*d1)min
f(x(1)d1)
其中
f(x(1)d1)f(24,216)(24)24(216)2.
求解得到:17/130, 从而
x22,2T
17/130*(-4,-16)T(96/65,6/65)T
在x2
96/65,6/65T
,f192/65-48/65,
从而最速下降法的搜索方向为:d1192/65
48/65
.
继续迭代,直至满足精度。
在初始点x12,2T
,
G201
1/2008,G01/8
从而修正的牛顿法的搜索方向为:
d1G1f1/204-2
01/816 -2
最优步长*满足
f(x(1)*d1)minf(x(1)d1
)
f(x(1)d1)f(22,22)5(22)2.1, 从而x222,22T
(0,0)T
在x20,0T
,f00,显然满足精度。
2
f20x
(2)正定,08
x(2)即为所求的极小点。(1分)
易见:
min
2.讨论约束极值问题
2
f(x)x12x26x16x28
s.t.
x1x24xx0
21
x10x02
的Kuhn-Tucker点。
f(x)[2x16,2x26]T
g1(x)x1x24, g1(x)[1,1]T g2(x)x1x2, g2(x)[1,1]T
g3(x)x1, g3(x)[1,0]T。 g4(x)x2, g4(x)[0,1]T
所以该约束问题的Kuhn-Tucker条件为
2x161110121304102x61
2
(xx4)0
2
11
2(x1x2)0
3x10
4x20
x1x240xx0
2
1
x10x021,2,3,40
若x10,x20,则x1x2440, 从而10。
代入第一个式子得到62306240
从而2460,矛盾。
若x10,x20,则x1x20从而20。又由于x20,从而40
代入第一个式子得到61302x2610
从而136,2x26362x230,与x20,30矛盾。
若x10,x20,则x1x20则20,30
代入第一个式子得到2x16106140
从而146,2x16462x140,与x10,40矛盾。
若x10,x20,且x1x20则20。又由于x10,x20从而30, 40,
代入第一个式子得到2x16102x2610
从而x1x2,与x1x20矛盾。
若x10,x20,x1x2,但x1x24则10。又由于x10,x20从而30, 40,
代入第一个式子得到2x16202x2620从而12,
x1x23,这与x1x24矛盾。
x10,
x20,x1x24,
讨论了所有的情况后,我们可以得出结论,从而(2,2)为该问题的一个KuhnTucker点.
3.构造増广函数
x)(x1)2(x)22
123k(x1x22)若x1x220
k((x2(x2
11)23)
若x1x220若x1x220,则
由2(x11)0
k(x)2(x3)
20,可得
x11,x23,与x1x220矛盾.若x1x220,则
由2(x11)2k(x1x22)0
k(x)2(x3)2k(x1x22)
20,可得
x11
2x12,x112.从而x22
k12k这是对于固定的k,minxR
2
k(x)的解.lim x1k
lim
1
120
kk
x2 x122从而
x*(0,2)为原问题的最优解。
minf(x)x2
16x192x2
4.用内点法求解非线性规划
s.t. g1(x)x130 g2(x)x230
构造増广函数
2k(x)x16x192x2k(
1x1
3
)13x2
2xk16由(x213)k(x)0k,可得20(x
23)2xk
13k
2x23
2
这是对于固定的k,min
xR
2
k(x)的解.limk3
k0
x1lim0
3k2
k
lim1k0
xlim3
k0
2
3从而
x*(3,3)为原问题的最优解。
(4分)
(2分)
(4分)
(4分)
(2分)
(4分)
5. 构造増广函数
1212
x1x2vk(x1x21)k(x1x21)226
x1vk2k(x1x21)0
由k(x)10,可得xv2(xx1)2kk12
3
2vk3(2kvk)
x23x1,x1k.从而x2
18k18k
k(x)
(4分)
(2分)
这是对于固定的k,mink(x)的解.2
xR
k
lim x1lim
2kvk1
k184k
(注意vk是一个数)
(4分)
3
x23x1
4从而
6.解:首先化成标准形式
13
x*(,)为原问题的最优解。
44
minZ3x12x2x3Mx7Mx8
x1x2x3x46
x3x5x74x1
x2x3x6x83x,x,x,x,x,x,x,x0
以x1为换入变量,根据最小比值原则确定x7为换出变量。
以x2为换入变量,根据最小比值原则确定x4为换出变量。
检验数全部为正,但人工变量没有完全换出,说明此优化问题没有可行解(可以验证原问题中包含矛盾的条件),此最优单纯形表的最优基是
1101101
B(p2,p1,p8)010 ,B010
101111
四.应用题
解:设x1,x2 分别为该厂生产甲乙两种产品的数量。该问题的目标规划模型为: (2分)
minzP1d1P2d2P3(2d35d4)
(3分)
12x130x2d1d12400
x12x2d2d2200
(5分)
x1d3d360 x2d4d480
x1,x2,di,di0
(i1,2,3,4)
其中在P3级目标中,因甲产品的利润与乙产品利润的比值为2:5,故取权系数为2:5. 求解过程见图.
(5分)
满足P1,P2目标的解空间为三角形ABC区域,考虑P3的目标要求时,因d3的权系数小于d4 的权系数,故先取d40,这时解空间为ACD区域,在此区域中,只有D点使d3取
值最小,故取D点为满意解,其坐标为(40,80),即该厂每年应生产甲产品40个单位,乙产品80个单位.
(5分)
《最优化方法》(研究生)期末考试练习题答案
二.简答题
min -5y19y2, s.t. 4y13y2 3, -2y1 y22, 1.
3y14y28, y1,y20;
x3 x40, (以x1为源行生成的割平面方程) 2.
注意:在x1为整数的情况下,因为x3,x40,该方程自然满足,这是割平面的退化情形
1
656
111
x3 x4, (以x2为源行生成的割平面方程)
442
3.
a10,b13
1a10.382(b1a1)00.382*31.146
1a10.618(b1a1)00.618*31.854 (1)(1.146)32*1.14610.2131(1)(1.854)32*1.85413.6648
事实上,不经计算也可以看出
(1)(1),所以a20,b21.854。
即:初始的保留区间为[0,1.854]。近似的最优解:x*
01.854
0.927.2
f1(x)x1ex2*(1)2.7x1ex22.7
4.令
f2(x)x1ex2*(0)1x11f1(x)x1e
x2*(1)
0.4x1e
x2
0.4
f1(x)x1ex2*(2)0.1x1e2x20.1
拟合问题等价于求解下列最小二乘问题:min
(f(x))
i
i1
4
2
三.计算题
1.分别用最速下降方法和修正的牛顿法求解无约束问题 minf(x)x14x2。 取初始点x
2
2
1
T
2,2,0.1.
解:f2x1
在初始点x1
2,2T
8x2,
,
f4
16
从而最速下降法的搜索方向为:d14
16
.
最优步长*满足f(x(1)*d1)min
f(x(1)d1)
其中
f(x(1)d1)f(24,216)(24)24(216)2.
求解得到:17/130, 从而
x22,2T
17/130*(-4,-16)T(96/65,6/65)T
在x2
96/65,6/65T
,f192/65-48/65,
从而最速下降法的搜索方向为:d1192/65
48/65
.
继续迭代,直至满足精度。
在初始点x12,2T
,
G201
1/2008,G01/8
从而修正的牛顿法的搜索方向为:
d1G1f1/204-2
01/816 -2
最优步长*满足
f(x(1)*d1)minf(x(1)d1
)
f(x(1)d1)f(22,22)5(22)2.1, 从而x222,22T
(0,0)T
在x20,0T
,f00,显然满足精度。
2
f20x
(2)正定,08
x(2)即为所求的极小点。(1分)
易见:
min
2.讨论约束极值问题
2
f(x)x12x26x16x28
s.t.
x1x24xx0
21
x10x02
的Kuhn-Tucker点。
f(x)[2x16,2x26]T
g1(x)x1x24, g1(x)[1,1]T g2(x)x1x2, g2(x)[1,1]T
g3(x)x1, g3(x)[1,0]T。 g4(x)x2, g4(x)[0,1]T
所以该约束问题的Kuhn-Tucker条件为
2x161110121304102x61
2
(xx4)0
2
11
2(x1x2)0
3x10
4x20
x1x240xx0
2
1
x10x021,2,3,40
若x10,x20,则x1x2440, 从而10。
代入第一个式子得到62306240
从而2460,矛盾。
若x10,x20,则x1x20从而20。又由于x20,从而40
代入第一个式子得到61302x2610
从而136,2x26362x230,与x20,30矛盾。
若x10,x20,则x1x20则20,30
代入第一个式子得到2x16106140
从而146,2x16462x140,与x10,40矛盾。
若x10,x20,且x1x20则20。又由于x10,x20从而30, 40,
代入第一个式子得到2x16102x2610
从而x1x2,与x1x20矛盾。
若x10,x20,x1x2,但x1x24则10。又由于x10,x20从而30, 40,
代入第一个式子得到2x16202x2620从而12,
x1x23,这与x1x24矛盾。
x10,
x20,x1x24,
讨论了所有的情况后,我们可以得出结论,从而(2,2)为该问题的一个KuhnTucker点.
3.构造増广函数
x)(x1)2(x)22
123k(x1x22)若x1x220
k((x2(x2
11)23)
若x1x220若x1x220,则
由2(x11)0
k(x)2(x3)
20,可得
x11,x23,与x1x220矛盾.若x1x220,则
由2(x11)2k(x1x22)0
k(x)2(x3)2k(x1x22)
20,可得
x11
2x12,x112.从而x22
k12k这是对于固定的k,minxR
2
k(x)的解.lim x1k
lim
1
120
kk
x2 x122从而
x*(0,2)为原问题的最优解。
minf(x)x2
16x192x2
4.用内点法求解非线性规划
s.t. g1(x)x130 g2(x)x230
构造増广函数
2k(x)x16x192x2k(
1x1
3
)13x2
2xk16由(x213)k(x)0k,可得20(x
23)2xk
13k
2x23
2
这是对于固定的k,min
xR
2
k(x)的解.limk3
k0
x1lim0
3k2
k
lim1k0
xlim3
k0
2
3从而
x*(3,3)为原问题的最优解。
(4分)
(2分)
(4分)
(4分)
(2分)
(4分)
5. 构造増广函数
1212
x1x2vk(x1x21)k(x1x21)226
x1vk2k(x1x21)0
由k(x)10,可得xv2(xx1)2kk12
3
2vk3(2kvk)
x23x1,x1k.从而x2
18k18k
k(x)
(4分)
(2分)
这是对于固定的k,mink(x)的解.2
xR
k
lim x1lim
2kvk1
k184k
(注意vk是一个数)
(4分)
3
x23x1
4从而
6.解:首先化成标准形式
13
x*(,)为原问题的最优解。
44
minZ3x12x2x3Mx7Mx8
x1x2x3x46
x3x5x74x1
x2x3x6x83x,x,x,x,x,x,x,x0
以x1为换入变量,根据最小比值原则确定x7为换出变量。
以x2为换入变量,根据最小比值原则确定x4为换出变量。
检验数全部为正,但人工变量没有完全换出,说明此优化问题没有可行解(可以验证原问题中包含矛盾的条件),此最优单纯形表的最优基是
1101101
B(p2,p1,p8)010 ,B010
101111
四.应用题
解:设x1,x2 分别为该厂生产甲乙两种产品的数量。该问题的目标规划模型为: (2分)
minzP1d1P2d2P3(2d35d4)
(3分)
12x130x2d1d12400
x12x2d2d2200
(5分)
x1d3d360 x2d4d480
x1,x2,di,di0
(i1,2,3,4)
其中在P3级目标中,因甲产品的利润与乙产品利润的比值为2:5,故取权系数为2:5. 求解过程见图.
(5分)
满足P1,P2目标的解空间为三角形ABC区域,考虑P3的目标要求时,因d3的权系数小于d4 的权系数,故先取d40,这时解空间为ACD区域,在此区域中,只有D点使d3取
值最小,故取D点为满意解,其坐标为(40,80),即该厂每年应生产甲产品40个单位,乙产品80个单位.
(5分)