工程优化方法第二章_第1页
工程优化方法第二章_第2页
工程优化方法第二章_第3页
工程优化方法第二章_第4页
工程优化方法第二章_第5页
已阅读5页,还剩44页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

工程优化方法第二章第一页,共四十九页,2022年,8月28日§1.1多元函数的定义n元函数:

n元线性函数:

n元二次函数:

n元向量值线性函数:其中第二页,共四十九页,2022年,8月28日§1.2多元函数的可导性和可微性在点存在,的偏导数,记为的某邻域内极限则称此极限为函数设函数在点对第i个分量注意:(1)式也可写为其中定义1.2.1(偏导数)第三页,共四十九页,2022年,8月28日§1.2多元函数的可导性和可微性可表示成称为函数在点(x1,x2)的全微分,记作则称函数f(x1,x2

)在点(x1,x2)可微,定义(二元函数的可微性)如果二元函数z=f(x1,x2)在定义域D

的内点(x1,x2)处全增量第四页,共四十九页,2022年,8月28日§1.2多元函数的可导性和可微性定义中增量的表达式等价于记第五页,共四十九页,2022年,8月28日§1.2多元函数的可导性和可微性

若函数z=f(x1,x2)在点(x1,x2)可微,则该函数在该点偏导数必存在,称向量是函数z=f(x1,x2)在点(x1,x2)的梯度。且有二元多元可微第六页,共四十九页,2022年,8月28日§1.2多元函数的可导性和可微性定义(多元函数的可微性)

设若

使

有:则称f(x)在处可微。给定区域D上的n

元实值函数与二元函数可微的等价形式类似引入第七页,共四十九页,2022年,8月28日§1.2多元函数的可导性和可微性定理1.2.1(可微必可导)

若在处可微,则在该点处关于各变量的一阶偏导数存在,且

证明:令,依次取两边除以并取的极限有:

在处可微,则(3)对成立,第八页,共四十九页,2022年,8月28日§1.3多元函数的梯度定义(多元函数梯度)以的n个偏导数为分量的向量称为

f(x)在x处的梯度,若f在处可微,令p=x-x0,

由得记为注:梯度也可称为函数f(x)关于向量x

的一阶导数。这与一元函数展开到两项的Taylor

公式是相对应的。第九页,共四十九页,2022年,8月28日§1.3多元函数的梯度性质1的证明:过点的等值面方程为:

设f(x)在定义域内有连续偏导数,即有连续梯度,则梯度有以下两个重要性质:设是过点同时又完全在等值面(6)上的任一条光滑曲线L的方程,为参数,点对应的参数就是把此曲线方程代入(6),得到性质1:

函数在某点的梯度不为零,则必与过该点的等值面垂直。性质2:

梯度方向是函数具有最大变化率的方向。第十页,共四十九页,2022年,8月28日§1.3多元函数的梯度即函数f(x)在处的梯度与过该点在等值面上的任一条曲线L在此点的切线垂直。从而与过该点的切平面垂直,性质1成立。两边同时在处关于求导数,根据求导的链式法则有:向量恰为曲线L在处的切向量,则第十一页,共四十九页,2022年,8月28日§1.3多元函数的梯度定义1.3.2(方向导数)设在点x处可微,p=te为固定向量,其中t是向量p的模,e

为向量p的单位向量,则称极限:注:若则f(x)从出发在附近沿p方向是下降的。为说明性质2:

梯度方向是函数具有最大变化率的方向为函数f(x)在点处沿方向p的方向导数,记为,

若则f(x)从出发在附近沿p方向是上升的。引进方向导数当t>0充分小时,有第十二页,共四十九页,2022年,8月28日§1.3多元函数的梯度

若则f(x)从出发在附近沿p方向是下降的。

若则f(x)从出发在附近沿p方向是上升的。

因此又将方向导数称为f(x)在处沿方向p的变化率。

方向导数正负决定了函数升降;升降速度的快慢由方向导数绝对值大小来决定,绝对值越大升降速度越大;第十三页,共四十九页,2022年,8月28日§1.3多元函数的梯度定理

若在点处可微,则其中e

为p方向上的单位向量。证明:f在可微,则根据可微定义,容易看到:当时,有由前面证明即知p为下降方向。利用方向导数定义并将上式中的p换成te有:第十四页,共四十九页,2022年,8月28日§1.3多元函数的梯度由于,β为方向p与的夹角。

从而梯度方向是函数具有最大变化率的方向,性质2成立。推论

若,则p是函数f(x)在处的下降方向;若,则p是函数f(x)在处的上升方向。

可见梯度方向即为函数的最速上升方向;负梯度方向即为函数的最速下降方向。

当夹角为0(β=0o)

,即沿梯度方向()时,方向导数取得最大值;当夹角为180o(β=180o)

,即沿负梯度方向()时,方向导数取得最小值。第十五页,共四十九页,2022年,8月28日§1.3多元函数的梯度

上升方向变化率为0方向下降方向结论:函数在与其梯度正交的方向上变化率为0;成锐角的方向上是上升的;成钝角的方向上是下降的。

第十六页,共四十九页,2022年,8月28日§1.3多元函数的梯度解:由于则函数在处的最速下降方向例1.3.1试求目标函数在点处的最速下降方向,并求沿这个方向移动一个单位长度后新点的目标函数值。此方向上的单位向量新点是第十七页,共四十九页,2022年,8月28日§1.3多元函数的梯度几个常用的梯度公式:第十八页,共四十九页,2022年,8月28日§1.4多元函数的Hesse矩阵定义(Hesse矩阵)多元函数记f(x)的Hesse矩阵为第十九页,共四十九页,2022年,8月28日§1.4多元函数的Hesse矩阵常用的梯度和Hesse阵公式:第二十页,共四十九页,2022年,8月28日§1.4多元函数的Hesse矩阵多元函数的Taylor展开:设

二阶可导。在x*的邻域内Lagrange余项

对x,,记xx*+(x-x*)一阶Taylor展开式二阶Taylor展开式一阶中值公式对x,,使第二十一页,共四十九页,2022年,8月28日第二章基本概念和理论基础本章主要内容:§1多元函数的梯度及其Hesse矩阵§2多元函数的极值及其判别条件§3等高线§4多元函数分析(二次函数)§5凸集、凸函数、凸规划§6几个重要的不等式第二十二页,共四十九页,2022年,8月28日§2多元函数的极值及其判别条件

对于一个极小化问题,我们希望知道的是全局极小点,而到目前为止的一些最优化算法却基本上是求局部极小值点的。因此一般要先求出所有局部极小值点,再从中找出全局极小点。

为了求出函数的局部极小值点,考察函数f在局部极小点处满足什么条件?反过来,满足什么条件的点是局部极小点?

这就是接下来我们要考虑的多元函数的极值条件。首先回顾二元函数的极值条件。第二十三页,共四十九页,2022年,8月28日§2.1二元函数的极值判别条件定理2.1.1(必要条件)

设(1)为D的一个内点;可微;(2)在处,则在的极值点;(3)为且注:可微的极值点一定是驻点,反之不一定成立。第二十四页,共四十九页,2022年,8月28日§2.1二元函数的极值判别条件定理2.1.2(充分条件)

设(1)为D的一个内点;二次连续可微;(2)在的驻点,即(3)为且令则(1)当时,具有极值取严格极大值取严格极小值(2)当时,不是的极值点,称为函数的鞍点;(3)当时,不能确定,需另行讨论。第二十五页,共四十九页,2022年,8月28日§2.2多元函数的极值判别条件定理2.2.1(必要条件)

设(1)x*为D的一个内点;可微;(2)在则的极值点;(3)为且定义

设是D的内点,若则称为f的驻点。第二十六页,共四十九页,2022年,8月28日§2.2多元函数的极值判别条件

设为任意单位向量,是的局部极小点。由定义知:当,即时,总有:令则而是D的内点,从而与之对应的t=0是

的局部极小点。

由一元函数极小点必要性条件知

,而由前述性质知

,由单位向量任意性,即知。证明:(若,取,则矛盾。)第二十七页,共四十九页,2022年,8月28日§2.2多元函数的极值判别条件例

在处梯度为,但只是双曲抛物面的鞍点,而不是极小点。f第二十八页,共四十九页,2022年,8月28日§2.2多元函数的极值判别条件定理2.2.2(充分条件)

设(1)x*为D的一个内点;(2)二次连续可微;在(3)则的严格局部极小点。为(4)正定;第二十九页,共四十九页,2022年,8月28日§2.2多元函数的极值判别条件证明:因正定,则使对,均有:(x充分接近时)。将f在处按Taylor公式展开注意,有:当x充分接近时,上式左端的符号取决于右端的第一项(为正)。故第三十页,共四十九页,2022年,8月28日课后作业P382.12.32.9--2.14第三十一页,共四十九页,2022年,8月28日第二章基本概念和理论基础本章主要内容:§1多元函数的梯度及其Hesse矩阵§2多元函数的极值及其判别条件§3等高线§4多元函数分析(二次函数)§5凸集、凸函数、凸规划§6几个重要的不等式第三十二页,共四十九页,2022年,8月28日§3等高线例3.1

求解,这是定义在平面上的无约束极小化问题,其目标函数在三维空间中代表一个曲面。

二元函数最优化问题,具有明显的几何特征,从几何图形上,可以直观了解函数的变化,我们把这种几何解释推广到n维空间中,对后面优化方法的研究是有益处的。第三十三页,共四十九页,2022年,8月28日§3等高线在平面上任给一点,就对应有一个目标函数值是过点作平面的垂线与S曲面交点的纵坐标。反之,任给一个值f0,使目标函数f(z)取值为f0的点z的个数就不相同了。可能没有,可能只有一个,可能有多个。

这一事实的几何意义是:过f

轴上坐标为f0的点作坐标平面的平行平面L,可能与曲面S无交点(f0<0),可能与S有一个交点(f0=0),可能与S交成一条曲线(f0>0).第三十四页,共四十九页,2022年,8月28日§3等高线定义3.1(等值线)平面L截曲面S得到一个圆,将它投影到平面上,仍为同样大小的圆,在这个圆上每一点的目标函数值均为f0。若一条曲线上任何一点的目标函数值等于同一常数,则称此曲线为目标函数的等值线。我们感兴趣的是至少有一个交点(f0≥0)的情形。注意:变动

f的值,得到不同等值线,这是一组同心圆;对应f0=0的等值线缩为一点G;对应f0<0的等值线为空集.随着f值变小,等值线圆半径变小,最后缩为一点,即为问题的最小值点

G,即第三十五页,共四十九页,2022年,8月28日§3等高线例3.2

用图解法求解解:先画出目标函数等值线,再画出约束曲线。对应的最优值为由图易见约束直线与等值线的切点是最优点,利用解析几何的方法得该切点为本处约束曲线是一条直线,这条直线就是可行域;而最优点就是可行域上使等值线具有最小值的点.第三十六页,共四十九页,2022年,8月28日§3等高线定义3.2(等值面)在三维和三维以上的空间中,使目标函数取同一常数值的面{x|f(x)=r,r是常数}

称为目标函数的等值面。定理3.1

若多元函数在其极小点处的

Hesse阵正定,则它在这个极小点附近的等值面近似地呈现为同心椭球面族。第三十七页,共四十九页,2022年,8月28日§3等高线证明:设是多元函数f的极小点。并设是充分靠近极小点的一个等值面,即充分小。将在点展开因为极小值点,

这是等值面的一个近似曲面。由于假设正定,则

是以为中心的椭球面方程。又是高阶无穷小量,则第三十八页,共四十九页,2022年,8月28日§3等高线等值面的性质:不同值的等值面之间不相交,因为目标函数是单值函数;除了极值点所在的等值面外,不会在区域内部中断,因为目标函数是连续的;等值面稠的地方,目标函数值变化得较快,而稀疏的地方变化得比较慢;一般地,在极值点附近,等值面(线)近似地呈现为同心椭球面族(椭圆族);二次函数的等值面是同心椭球面族,极值点是这个椭圆球面的共同中心。

第三十九页,共四十九页,2022年,8月28日第二章基本概念和理论基础本章主要内容:§1多元函数的梯度及其Hesse矩阵§2多元函数的极值及其判别条件§3等高线§4多元函数分析§5凸集、凸函数、凸规划§6几个重要的不等式第四十页,共四十九页,2022年,8月28日§4.1二次函数在n元函数中,除了线性函数:

或最简单最重要的一类就是二次函数。第四十一页,共四十九页,2022年,8月28日§4.1二次函数定义(二次型)

代数学中将特殊的二次函数称为二次型。二次函数的一般形式为其中均为常数,向量矩阵表示形式其中Q为对称矩阵第四十二页,共四十九页,2022年,8月28日§4.1二次函数定义

设Q为n×n对称矩阵,若,均有,则称矩阵Q是正定的。若,均有

,则称矩阵Q是半正定的。若-Q是正定的,则称Q是负定的。若-Q是半正定的,则称Q是半负定的。对于二次函数,我们更关心的是Q为正定矩阵的情形。第四十三页,共四十九页,2022年,8月28日§4.1二次函数A是正定矩阵等价于以下4个条件成立:

(1)存在非奇异矩阵G,使得A=GTG;

(2)A的所有特征根大于零;

(3)有满秩矩阵G,使A=GTG;

(4)A的所有顺序主子式都大于零.怎么判定一个对称矩阵Q是不是正定的?Sylvester(西尔维斯特)定理:一个n×n对称矩阵Q是正定矩阵的充要条件是矩阵Q的各阶顺序主子式都是正的。第四十四页,共四十九页,2022年,8月28日§4.1二次函数解:对称矩阵Q的三个顺序主子式依次为例

判定矩阵是否正定:矩阵Q是正定的。第四十五页,共四十九页,2022年,8月28日§4.1二次函数证明:作变换,代入二次函数式中:根据解析几何知识,Q为正定矩阵的二次型的等值面是以坐标原点为中心的同心椭球面族。由于上式中的是常数,所以的等值面也是以为中心的同心椭球面族,回到原坐标系中去,原

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论