2000字范文,分享全网优秀范文,学习好帮手!
2000字范文 > matlab求多元函数驻点 机械优化设计课后习题答案.doc

matlab求多元函数驻点 机械优化设计课后习题答案.doc

时间:2021-09-18 23:47:36

相关推荐

matlab求多元函数驻点 机械优化设计课后习题答案.doc

第一章习题答案

1-1 某厂每日(8h制)产量不低于1800件。计划聘请两种不同的检验员,一级检验员的标准为:速度为25件/h,正确率为98%,计时工资为4元/h;二级检验员标准为:速度为15件/h,正确率为95%,计时工资3元/h。检验员每错检一件,工厂损失2元。现有可供聘请检验人数为:一级8人和二级10人。为使总检验费用最省,该厂应聘请一级、二级检验员各多少人?

解:(1)确定设计变量;

根据该优化问题给定的条件与要求,取设计变量为X = ;

(2)建立数学模型的目标函数;

取检验费用为目标函数,即:

f(X) = 8*4*x1+ 8*3*x2 + 2(8*25*0.02x1 +8*15*0.05x2 )

=40x1+ 36x2

(3)本问题的最优化设计数学模型:

min f (X) = 40x1+ 36x2 X∈R3·

s.t. g1(X) =1800-8*25x1+8*15x2≤0

g2(X) =x1 -8≤0

g3(X) =x2-10≤0

g4(X) = -x1 ≤0

g5(X) = -x2 ≤0

1-2 已知一拉伸弹簧受拉力,剪切弹性模量,材料重度,许用剪切应力,许用最大变形量。欲选择一组设计变量使弹簧重量最轻,同时满足下列限制条件:弹簧圈数,簧丝直径,弹簧中径。试建立该优化问题的数学模型。

注:弹簧的应力与变形计算公式如下

解: (1)确定设计变量;

根据该优化问题给定的条件与要求,取设计变量为X = ;

(2)建立数学模型的目标函数;

取弹簧重量为目标函数,即:

f(X) =

(3)本问题的最优化设计数学模型:

min f (X) = X∈R3·

s.t. g1(X) =0.5-x1 ≤0

g2(X) =10-x2 ≤0

g3(X) =x2-50 ≤0

g4(X) =3-x3 ≤0

g5(X) =≤0

g6(X) =≤0

1-3 某厂生产一个容积为8000 cm3的平底、无盖的圆柱形容器,要求设计此容器消耗原材料最少,试写出这一优化问题的数学模型。

解:根据该优化问题给定的条件与要求,取设计变量为X = ,

表面积为目标函数,即:

minf(X) = x12 + 2 x1 x2

考虑题示的约束条件之后,该优化问题数学模型为:

minf(X) = x12 + 2 x1 x2

X=[x1,x2]T∈R2

s.t. g1(X) = -x1 ≤0

g2(X) = -x2 ≤0

h1(X) = 8000 - x12 x2 = 0

1-4 要建造一个容积为1500 m3的长方形仓库,已知每平方米墙壁、屋顶和地面的造价分别为4元、6元和12元。基于美学的考虑,其宽度应为高度的两倍。现欲使其造价最低,试导出相应优化问题的数学模型。

解:(1)确定设计变量;

根据该优化问题给定的条件与要求,取设计变量为X = ;

(2)建立数学模型的目标函数;

取总价格为目标函数,即:

f(X) = 8(x1 x3 + x2 x3) + 6 x1 x2 + 12 x1 x2

(3)建立数学模型的约束函数;

1)仓库的容积为1500 m3。即:

1500-x1 x2 x3 =0

2)仓库宽度为高度的两倍。即:

x2 -2 x3 = 0

3)各变量取值应大于0,即:

x1 > 0, x2 .> 0.,则 -x1 ≤0,-x2 ≤0

(4)本问题的最优化设计数学模型:

min f (X) = 8(x1 x3 + x2 x3) + 18 x1 x2 X∈R3·

s.t. g1(X) = -x1 ≤0

g2(X) = -x2 ≤0

g3(X) = -x3 ≤0

h1(X) = 1500-x1 x2 x3 =0

h2(X) = x2 -2 x3 = 0

1-5 绘出约束条件:

; ; 所确定的可行域

1-6 试在三维设计空间中,用向量分别表示设计变量:

; ; 。

第二章习题答案

2-1 请作示意图解释:的几何意义。

2-2 已知两向量,求该两向量之间的夹角。

2-3 求四维空间内两点和之间的距离。

2-4 计算二元函数在处,沿方向的方向导数和沿该点梯度方向的方向导数。

2-5 已知一约束优化设计问题的数学模型为

求:

(1) 以一定的比例尺画出当目标函数依次为时的四条等值线,并在图上画出可行区的范围。

(2) 找出图上的无约束最优解和对应的函数值,约束最优解和;

(3) 若加入一个等式约束条件:

求此时的最优解,。

解:下图为目标函数与约束函数(条件)设计平面X1OX2 。其中的同心圆是目标函数依次为f(X)=1、2、3、4时的四条等值线;阴影的所围的部分为可行域。

由于目标函数的等值线为一同心圆,所以无约束最优解为该圆圆心即:

X1*=[3,4]T

函数值 f(X1*)= 0 。

而约束最优解应在由约束线g1(X)=0,g2(X)=0,g3(X)=0,g4(X)=0,组成的可行域(阴影线内侧)内寻找,即约束曲线g1(X)=0与某一等值线的一个切点X2*,可以联立方程: ,解得X2*=[2,3] 。

函数值 f(X2*)= (2-3)2 + (3-4)2 = 2 。

加入等式约束条件,则X3*为可行域上为h1(X)=0上与某一条等值线的交点,可以联立方程: , 解得X3*=[5/2,5/2] 。

函数值 f(X3*)= (5/2-3)2 + (5/2-4)2 = 2.5 。

2-6 试证明在点处函数具有极小值。

证明:求驻点:,

H(X)是正定的, 所以驻点必定是极小点。故在点处函数具有极小值。

2-7 求函数的极值点,并判断其极值的性质。

解:,

H(X)是正定的,所以,为凸函数。

2-8 试判断函数的凸性。

解:,

H(X)是正定的,

所以,为凸函数。

2-9 试用向量及矩阵形式表示并证明它在上是一个凸函数。

解:,

H(X)是正定的,

所以,为凸函数。

2-10 现已获得优化问题

的一个数值解,试判定该解是否上述问题的最优解。

第三章习题答案

3-1 函数,当初始点分别为及时,用进退法确定其一维优化的搜索区间,取初始步长。

解:当时

(1)取

=0.1

比较,因 ,所以应作前进搜索。

⑵步长加倍:

=0.3

再比较,因,所以还应再向前搜索,为此应舍去上一次的点。所以:。

(3) 步长加倍:

=0.7

.

比较,因,所以还应再向前搜索, 。

(4) 步长加倍:

=1.5

.

比较,因。已找到具有“高-低-高”特征的区间

即:时,

时,

时,。

所以,,单峰区间为:

当时

同理可得:

3-2 用黄金分割法求函数在区间中的极小点,要求计算到最大未确定区间长度小于0.05。

解:(1)在初始区间[a,b]=[-3,5]中取计算点并计算函数值

(2)比较函数值,缩短搜索区间

因有f1≤f2,则

(3)判断迭代终止条件

b-a>ε

不满足迭代终止条件,比较函数值f1、f2继续缩短区间。

将各次缩短区间的有关计算数据列于下表。

表 黄金分割法的搜索过程

区间缩短次数

a

b

α (1)

α (2)

f1

f2

(原区间)

-3

5

0.056

1.944

0.115

7.667

1

-3

1.944

-1.111

0.056

-0.987

0.115

2

-3

0.056

-1.832

-1.111

-0.306

-0.987

3

-1.832

0.056

-1.111

-0.665

-0.987

-0.888

4

-1.832

-0.665

-1.386

-1.111

-0.851

-0.987

(5-8)略

9

-1.11122

-0.94097

-1.046

-1.006

-0.997867

-0.999964

3-3 用二次插值法求函数的最优解。已知搜区间为,选代精度。

解:采用Matlab编程计算得:

3-4 函数,取初始点为,规定沿点的负梯度方向进行一次一维优化搜索,选代精度:。

(1)用进退法确定一维优化搜索区间;

(2)用黄金分割法求最优化步长及一维优化最优值;

(3)用二次插值法求最优化步长及一维优化最优值;

(4)上述两种一维优化方法在求解本题时,哪一个种方法收取更快,原因是什么?

解:最优点,最优值

二次插值法更快

3-5 求的极小点,选代精度。要求:

(1)从出发,为步长确定搜索区间;

(2)用黄金分割法求极值点;

(3)用二次插值法求极值点。

解:

(1) ①由已知条件可得,

因为,应作前进搜索。

②步长加倍,,

因为,所以还应再向前搜索,为此应舍去上一次的点。所以:

③步长加倍,,

因为,所以还应再向前搜索,为此应舍去上一次的点。所以:

④步长加倍,,

因为,所以还应再向前搜索,为此应舍去上一次的点。所以:

⑤步长加倍,,

因为,所以已找到具有“高—低—高”特征的区间

即时,;

时,;

时,。

(2)由(1)确定的搜索区间[0.7,3.1],利用Matlab进行黄金分割法一维优化搜索得:

(3)由(1)确定的搜索区间[0.7,3.1],利用Matlab进行二次插值法一维优化搜索得:

9

展开阅读全文

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。