1－1．简述优化设计问题数学模型的表达形式。

答：优化问题的数学模型是实际优化设计问题的数学抽象。

在明确设计变量、约束条件、目标函数之后,优化设计问题就可以表示成一般数学形式。

求设计变量向量[]12Tn x x x x =使()min f x →且满足约束条件()0(1,2,)k h x k l ==()0(1,2,)j g x j m ≤=利用可行域概念,可将数学模型的表达进一步简练。

设同时满足()0(1,2,)j g x j m ≤=和()0(1,2,)k h x k l ==的设计点集合为R ,即Ｒ为优化问题的可行域，则优化问题的数学模型可简练地写成求x 使 min()x Rf x ∈符号“∈”表示“从属于”。

在实际优化问题中,对目标函数一般有两种要求形式：目标函数极小化()minf x →或目标函数极大化()max f x →。

由于求()f x 的极大化与求()f x -的极小化等价，所以今后优化问题的数学表达一律采用目标函数极小化形式。

１-２.简述优化设计问题的基本解法。

（不要抄书,要归纳) 答：求解优化问题可以用解析解法，也可以用数值的近似解法。

解析解法就是把所研究的对象用数学方程（数学模型)描述出来，然后再用数学解析方法（如微分、变分方法等）求出有化解。

但是，在很多情况下,优化设计的数学描述比较复杂,因而不便于甚至不可能用解析方法求解；另外，有时对象本身的机理无法用数学方程描述，而只能通过大量试验数据用插值或拟合方法构造一个近似函数式,再来求其优化解,并通过试验来验证;或直接以数学原理为指导，从任取一点出发通过少量试验(探索性的计算),并根据试验计算结果的比较,逐步改进而求得优化解。

这种方法是属于近似的、迭代性质的数值解法。

数值解法不仅可用于求复杂函数的优化解，也可以用于处理没有数学解析表达式的优化问题。

因此,它是实际问题中常用的方法,很受重视。

其中具体方法较多,并且目前还在发展。

但是,应当指出,对于复杂问题,由于不能把所有参数都完全考虑并表达出来,只能是一个近似的最后的数学描述。

由于它本来就是一种近似,那么，采用近似性质的数值方法对它们进行解算，也就谈不到对问题的精确性有什么影响了。

不管是解析解法，还是数值解法，都分别具有针对无约束条件和有约束条件的具体方法。

可以按照对函数倒数计算的要求，把数值方法分为需要计算函数的二阶导数、一阶导数和零阶导数(即只要计算函数值而不需计算其导数)的方法。

2-1.何谓函数的梯度？梯度对优化设计有何意义?答:二元函数f(x １,x2)在x0点处的方向导数的表达式可以改写成下面的形式⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡∂∂∂∂=∂∂+∂∂=∂∂2cos 1cos 212cos 21cos 1θθθθxo x f x f xo x f xo x f xo d f 令xo T x f x f x f x fx f ⎥⎦⎤⎢⎣⎡∂∂∂∂=∂∂∂∂=∇21]21[)0(并称它为函数f （x1，x2）在x0点处的梯度。

假设⎥⎦⎤⎢⎣⎡=2cos 1cos θθd为D 方向上的单位向量,则有d T x f xo d f)0(∇=∂∂ 即函数f （x1,x2）在x ０点处沿某一方向d 的方向导数xo df∂∂等于函数在该点处的梯度)0(x f ∇与d 方向单位向量的内积。

梯度方向是函数值变化最快的方向,而梯度的模就是函数变化率的最大值。

梯度与切线方向d 垂直,从而推得梯度方向为等值面的法线方向。

梯度)0(x f ∇方向为函数变化率最大方向，也就是最速上升方向。

负梯度-)0(x f ∇方向为函数变化率最小方向,即最速下降方向。

2－2.求二元函数2122212122),(x x x x x x f +-+=在T x ]0,0[0=处函数变化率最大的方向和数值。

解;由于函数变化率最大的方向就是梯度的方向,这里用单位向量p 表示,函数变化率最大和数值时梯度的模)0(x f ∇。

求f （ｘ1，x ２）在ｘ0点处的梯度方向和数值,计算如下：()⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡+-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡∂∂∂∂=∇120122214210x x x x fx f x f2221)0(⎪⎭⎫ ⎝⎛∂∂+⎪⎭⎫ ⎝⎛∂∂=∇x f x f x f =5⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-=∇∇=5152512)0()0(x f x f p 2－3．试求目标函数()2221212143,x x x x x x f +-=在点X ０=[1，０]T处的最速下降方向,并求沿着该方向移动一个单位长度后新点的目标函数值。

解：求目标函数的偏导数21221124,46x x x f x x x f +-=∂∂-=∂∂则函数在X 0＝[1,0]Ｔ处的最速下降方向是⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡-+-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡∂∂∂∂-=-∇=====462446)(012121012102121x x x xx x x x x f x f X f P 这个方向上的单位向量是：13]2,3[4)6(]4,6[T 22T -=+--==P P e 新点是⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=+=132133101e X X 新点的目标函数值1321394)(1-=X f2－4.何谓凸集、凸函数、凸规划?（要求配图）一个点集(或区域),如果连接其中任意两点x １、x2的线段都全部包含在该集合内，就称该点集为凸集，否则为非凸集。

函数f （ｘ)为凸集定义域内的函数，若对任何的01α≤≤及凸集域内的任意两点x １、ｘ２，存在如下不等式:称f （x ）是定义在图集上的一个凸函数。

()()()121211fx x f x x αααα+-≤+-⎡⎤⎣⎦对于约束优化问题min ()s.t. ()0 (1,2,,)j f x g x j m ≤=若()j=j f x g x 、（） 1,2,...,m都是凸函数，则称此问题为凸规划。

３-1.简述一维搜索区间消去法原理。

（要配图)答:搜索区间（a,ｂ)确定之后,采用区间逐步缩短搜索区间,从而找到极小点的数值近似解。

假设搜索区间（a ，ｂ)内任取两点ａ1,ｂ1 ，a １《ｂ1,并计算函数值f （a1),f(ｂ１)。

将有下列三种可能情形； 1）f （a1）《f(ｂ１）由于函数为单谷，所以极小点必在区间(ａ,b １）内 2）f （a1)》f(b1）,同理,极小点应在区间(a1，b)内 3)f （a1)=ｆ（b １),这是极小点应在（ａ1，ｂ1)内3－2.简述黄金分割法中0．6１8的来由,搜索过程及程序框图。

黄金分割法适用于[],a b 区间上的任何单谷函数求极小值问题。

对函数除要求“单谷”外不作其他要求，甚至可以不连续。

因此,这种方法的适应面相当广。

黄金分割法也是建立在区间消去法原理基础上的试探方法，即在搜索区间[],a b 内适当插入两点1α、2α,并计算其函数值。

1α、2α将区间分成三段。

应用函数的单谷性质,通过函数值大小的比较，删去其中一段，使搜索区间得以缩短。

然后再在保留下来的区间上作同样的处置，如此迭代下去,使搜索区间无限缩小,从而得到极小点的数值近似解。

黄金分割法要求插入点1α、2α的位置相对于区间[],a b 两端点具有对称性，即1()b b a αλ=--2()a b a αλ=+-其中，λ为待定常数。

图a除对称要求外,黄金分割法还要求在保留下来的区间内再插入一点所形成的区间新三段，与原来区间的三段具有相同的比例分布。

设原区间[],a b 长度为1,如图a 所示，保留下来的区间[]2,a α长度为，区间缩短率为λ。

为了保持相同的比例分布,新插入点3α应在(1)λλ-位置上，1α在原区间的1λ-位置应相当于在保留区间的2λ位置。

故有21λλ-= 210λλ+-=取方程正数解，得510.618λ-=≈若保留下来的区间为[]1,b α，根据插入点的对称性，也能推得同样的λ值。

所谓“黄金分割”是指将一线段分成两段的方法,使整段长与较长段的长度比值等于较长段与较短段长度的比值，即 1::(1)λλλ=-同样算得0.618λ≈。

可见黄金分割法能使相邻两次搜索区间都具有相同的缩短率0.６18,所以黄金分割法又被称作0.618法。

图 ｂ黄金分割法的搜索过程是： （1） 给出初始搜索区间[],a b 及收敛精度ε，将λ赋以０.618。

（2）按坐标点计算公式1()b b a αλ=--、2()a b a αλ=+-计算1α和2α，并计算其对应的函数值1()f α，2()f α。

（3）根据区间消去法原理缩短搜索区间。

为了能用原来的坐标点计算公式，需进行区间名称的代换,并在保留区间中计算一个新的试验点及其函数值。

（4） 检查区间是否缩短到足够小和函数值收敛到足够近,如果条件不满足则返回到步骤（２)。

（5） 如果条件满足，则取最后两试验点的平均值作为极小点的数值近似解。

（6）黄金分割法的程序框图如图ｂ所示。

3－３．对函数ααα2)(2+=f ，当给定搜索区间55≤≤-α时，写出用黄金分割法求极小点*α的前三次搜索过程。

(要列表)解;此时的a=-5，b=5。

首先插入两点a １和a2。

可得ａ1=b －)(a b -λ=-1.18 , ａ2=a+)(a b -λ＝１.１８再计算相应插入点的函数值，得y １=f(a1)=-0.９６76,y2＝ｆ(ａ2)＝3．7524因为y ２>y １，所以消去区间[a2,b]，则新的搜索区间[a ，b]的端点ａ＝-５不变，而端点b=a ２=１．18第一次迭代;此时插入点a1＝ｂ-)(a b -λ=-2.６3９,a ２=-1．1８1。

相应插入点的函数值y1=f （ａ1)＝1．6８6，y2=f(a2)=-０.９67，由于ｙ1>y2，故消去区间[ａ,a １],新的搜索区间为[-２.639,1．18]，如此继续迭代下去 列出前三次迭代结果迭代序号 ａ a1 a ２ b Y1 比较 Y20 -5 -1.１8 １.１8 5 -０.９67６ < ３．7524 1 －５ -2.63９ -１．181 1.18 1．686 > -0．９67 2 －2.63９ -1.1８ -０.２79 1.1８ -０．9６7６ ＜ -０.48 3 -2.63９ －１.７37 -1.1８1 -0.２7９ -0.45７ > -０.４8２ 3-4.使用二次插值法求f (x ）=sin(x )在区间［2,6]的极小点，写出计算步骤和迭代公式,给定初始点x １＝2,x 2=4，x 3=6， ε=１0-４。

１ 2 ３ 4 x 1 2 4 4.55４57 4．55457迭代次数Ｋ= 4 ,极小点为 4.７１236 ,最小值为 -113131x x y y c --=，12122x x y y c --=,32123x x c c c --=)(213131c cx x x p -+=收敛的条件：ε<-22y y y p４－１.简述无约束优化方法中梯度法、共轭梯度法、鲍威尔法的主要区别。