汽车优化设计 第二章:优化方法的数学基础 王琥 湖南大学 机械与运载工程学院

Slides:



Advertisements
Similar presentations
一、 一阶线性微分方程及其解法 二、 一阶线性微分方程的简单应用 三、 小结及作业 §6.2 一阶线性微分方程.
Advertisements

目录 上页 下页 返回 结束 习题课 一、导数和微分的概念及应用 二、导数和微分的求法 导数与微分 第二章.
Yunnan University Chapt 5. 微分学基本定理及其应用 导 数导 数 函数性质 中值定理 §1. 中值定理 §2. 泰勒公式 §3. 函数的升降、凸性与极值 §4. 平面曲线的曲率 §5. 待定型.
第八章 习题课 多元函数微分学. 一 基本要求 1 理解二元函数的概念,会求定义域。 2 了解二元函数的极限和连续的概念。 3 理解偏导数的概念,掌握偏导数及高阶偏导 数的求法。 4 掌握多元复合函数的微分法。 5 了解全微分形式的不变性。 6 掌握隐函数的求导法。
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第八章 第四节 机动 目录 上页 下页 返回 结束 一个方程所确定的隐函数 及其导数 隐函数的微分法.
高等数学一 主讲 杨俊 演示文稿制作 杨俊. 高等数学一 第 3 章 一元函数微分学的应用 第 4 章 一元函数 积分学及应用 第 1 章 函数、极限与连续 第 2 章 导数与微分.
一、会求多元复合函数一阶偏导数 多元复合函数的求导公式 学习要求: 二、了解全微分形式的不变性.
2.6 隐函数微分法 第二章 第二章 二、高阶导数 一、隐式定义的函数 三、可微函数的有理幂. 一、隐函数的导数 若由方程 可确定 y 是 x 的函数, 由 表示的函数, 称为显函数. 例如, 可确定显函数 可确定 y 是 x 的函数, 但此隐函数不能显化. 函数为隐函数. 则称此 隐函数求导方法.
第十二章 第二节 一元函数 y = f (x) 的微分 机动 目录 上页 下页 返回 结束 对二元函数的全增量是否也有类似这样的性质? 全微分.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二章 导数与微分 一. 内 容 要 点 二. 重 点 难 点 三. 主 要 内 容 四. 例 题与习题.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
第四章 多元函数微分学 一元函数微分学 推广 多元函数微分学 注意: 善于类比, 区别异同 一元函数、极限与连续 一元函数的导数
《解析几何》 -Chapter 3 §7 空间两直线的相关位置.
第八章 向量代数 空间解析几何 第五节 空间直线及其方程 一、空间直线的点向式方程 和参数方程 二、空间直线的一般方程 三、空间两直线的夹角.
3.4 空间直线的方程.
= =9 =1 s s L 在 平面上任给一点 ,就对应有一个目标函数值 = 这个值就是过 点作 平面的垂线与S曲面交点的纵坐标。
圆的一般方程 (x-a)2 +(y-b)2=r2 x2+y2+Dx+Ey+F=0 Ax2+Bxy+Cy2+Dx+Ey+ F=0.
《解析几何》 乐山师范学院 0 引言 §1 二次曲线与直线的相关位置.
第五章 二次型. 第五章 二次型 知识点1---二次型及其矩阵表示 二次型的基本概念 1. 线性变换与合同矩阵 2.
§ 18.4 条件极值 一、极值 二、 条件极值拉格朗日乘数法.
第3节 二次型与二次型的化简 一、二次型的定义 二、二次型的化简(矩阵的合同) 下页.
一、二阶行列式的引入 用消元法解二元线性方程组. 一、二阶行列式的引入 用消元法解二元线性方程组.
§ 18.4 条件极值 一、极值 二、 条件极值拉格朗日乘数法.
第三章 函数逼近 — 最佳平方逼近.
数学建模方法及其应用 韩中庚 编著.
第四节 对数留数与辐角原理 一、对数留数 二、辐角原理 三、路西定理 四、小结与思考.
恰当方程(全微分方程) 一、概念 二、全微分方程的解法.
第五节 微积分基本公式 、变速直线运动中位置函数与速度 函数的联系 二、积分上限函数及其导数 三、牛顿—莱布尼茨公式.
一、原函数与不定积分 二、不定积分的几何意义 三、基本积分公式及积分法则 四、牛顿—莱布尼兹公式 五、小结
第四章 函数的积分学 第六节 微积分的基本公式 一、变上限定积分 二、微积分的基本公式.
第四章 一元函数的积分 §4.1 不定积分的概念与性质 §4.2 换元积分法 §4.3 分部积分法 §4.4 有理函数的积分
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
第六章 微分中值定理及其应用.
第三节 格林公式及其应用(2) 一、曲线积分与路径无关的定义 二、曲线积分与路径无关的条件 三、二元函数的全微分的求积 四、小结.
1.5 场函数的高阶微分运算 1、场函数的三种基本微分运算 标量场的梯度f ,矢量场的散度F 和F 旋度简称 “三度” 运算。
多元函数微分学学习辅导 一、内容提要 二、典型例题 首页 上页 返回 下页 结束.
全 微 分 欧阳顺湘 北京师范大学珠海分校
2-7、函数的微分 教学要求 教学要点.
第5章 §5.3 定积分的积分法 换元积分法 不定积分 分部积分法 换元积分法 定积分 分部积分法.
全国高校数学微课程教学设计竞赛 知识点名称: 导数的定义.
计算机数学基础 主讲老师: 邓辉文.
§2 求导法则 2.1 求导数的四则运算法则 下面分三部分加以证明, 并同时给出相应的推论和例题 .
第一章 函数 函数 — 研究对象—第一章 分析基础 极限 — 研究方法—第二章 连续 — 研究桥梁—第二章.
第八模块 复变函数 第二节 复变函数的极限与连续性 一、复变函数的概念 二、复变函数的极限 二、复变函数的连续性.
第4章 非线性规划 4.5 约束最优化方法 2019/4/6 山东大学 软件学院.
正切函数的图象和性质 周期函数定义: 一般地,对于函数 (x),如果存在一个非零常数T,使得当x取定义域内的每一个值时,都有
第16讲 相似矩阵与方阵的对角化 主要内容: 1.相似矩阵 2. 方阵的对角化.
§6.7 子空间的直和 一、直和的定义 二、直和的判定 三、多个子空间的直和.
3.1 变化率与导数   3.1.1 变化率问题 3.1.2 导数的概念.
函 数 连 续 的 概 念 淮南职业技术学院.
第一节 不定积分的概念与性质 一、原函数与不定积分的概念 二、不定积分的几何意义 三、基本积分表 四、不定积分的性质 五、小结 思考题.
第三章 函数的微分学 第二节 导数的四则运算法则 一、导数的四则运算 二、偏导数的求法.
4) 若A可逆,则 也可逆, 证明: 所以.
学习任务三 偏导数 结合一元函数的导数学习二元函数的偏导数是非常有用的. 要求了解二元函数的偏导数的定义, 掌握二元函数偏导数的计算.
1.非线性规划模型 2.非线性规划的Matlab形式
建模常见问题MATLAB求解  .
一元二次不等式解法(1).
第 六 章 约束最优化方法.
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
§2 方阵的特征值与特征向量.
滤波减速器的体积优化 仵凡 Advanced Design Group.
第三节 函数的微分 3.1 微分的概念 3.2 微分的计算 3.3 微分的应用.
选修1—1 导数的运算与几何意义 高碑店三中 张志华.
第四节 向量的乘积 一、两向量的数量积 二、两向量的向量积.
第三节 数量积 向量积 混合积 一、向量的数量积 二、向量的向量积 三、向量的混合积 四、小结 思考题.
第八节 第九章 多元函数的极值及其求法 一、多元函数的极值 二、最值应用问题 三、条件极值.
Presentation transcript:

汽车优化设计 第二章:优化方法的数学基础 王琥 wanghu@hnu.edu.cn 湖南大学 机械与运载工程学院 汽车车身先进设计制造国家重点实验室

主要内容 优化设计=极值问题 多变量、多约束非线性优化 高等数学极值理论是求解基础,但是不能直接求出最优解。 对多变量约束优化问题的求解方法所涉及的数学概念及有关理论进行补充和扩展。 方向导数与梯度;二次函数及正定矩阵;无约束优化问题的极值条件 ;等式约束优化问题的极值条件 ;不等式约束优化问题的极值条件;优化设计问题的基本解法

方向导数与梯度 方向导数 二元函数在点x0处沿某一方向s的方向导数 方向导数是偏导数概念的推广 方向导数与偏导数之间的数量关系是 O x2 1 2

方向导数与梯度 一个三元函数f(x1, x2,x3)在x0(x10, x20, x30) 点处沿s方向的方向导数为 n元函数在点x0处沿s方向的方向导数

梯度 二元函数的梯度 为函数F(x1,x2)在x0点处的梯度。

梯度 设 梯度的模: 设: 为单位向量。 则有

梯度 若上式为0,则说明方向导数是沿着等值线的切线向,而梯度是沿着与等值线切线相垂直的方向,且这时方向导数达到最大值,这说明梯度是函数值变化最快的方向,而梯度的模就是函数变化率的最大值 。 当方向s与梯度方向的夹角为锐角时,上式大于0,当方向s与梯度的夹角为钝角时,上式小于0,这说明与梯度成锐角的方向是函数值增加(上升)的方向,而与梯度成钝角的方向是函数值减小(下降)的方向。 梯度方向与等值线的关系

梯度 多元函数梯度 多元函数的方向导数表示为

梯度 梯度 模: 对于二元函数来说,函数的梯度方向与函数等值线面相垂直,也就是和等值面上过x0的一切曲线相垂直。 梯度 模: 对于二元函数来说,函数的梯度方向与函数等值线面相垂直,也就是和等值面上过x0的一切曲线相垂直。 由于梯度的模因点而异,即函数在不同点处的最大变化率是不同的。因此,梯度是函数的一种局部性质。

梯度 函数的梯度具有以下特征: 函数一点的梯度是由函数在该点上的所有一阶偏导数组成的向量。梯度的方向是该点函数值上升最快的方向,梯度的大小就是它的模。 函数在一点的梯度方向与函数过该点的等值线(面)的切线(平面)相垂直,或者说是该点等值线(面)的外法线方向。 梯度是函数在一点邻域内局部性态的描述。在邻域内上升得快的方向,离开邻域后就不一定上升得快,甚至可能下降。

梯度的计算实例 求函数 在点[3,2]T 的 梯度。 在点x(1)=[3,2]T处的梯度为: 解:

二次函数及正定矩阵 二次函数的一般形式为: 其中 均为常数。 其向量矩阵表示形式是: 其中 Q为对称矩阵 其中 均为常数。 其向量矩阵表示形式是: 其中 Q为对称矩阵 在代数学中将特殊的二次函数 称为二次型。对于二次函数,我们更关心的是Q为正定矩阵的情形。

二次函数及正定矩阵 定义:设Q为n×n对称矩阵 若 ,X≠0,均有 ,则称矩阵Q是正定的。 若 ,且 X≠0,均有 ,则称Q是负定的。 一个n×n对称矩阵Q是正定矩阵的充要条件是矩阵Q的各阶主子式都是正的。 一个n×n对称矩阵Q是负定矩阵的充要条件是矩阵Q的各阶主子式的值负、正相间。

算例 判定矩阵 是否正定 解:对称矩阵Q的三个主子式依次为: 因此知矩阵Q是正定的。

无约束优化问题的极值条件 多元函数的泰勒展开 二元函数: 二元函数:在点Xk处

无约束优化问题的极值条件 多元函数泰勒展开 海赛矩阵 (Hessian)

无约束优化问题的极值条件 对二次函数: Q 为二次函数的海赛(Hessian)矩阵,常量矩阵。 二次函数的梯度为:

算例 求目标函数 的梯度和Hessian矩阵。 解:因为 则 又因为: 故Hessian阵为:

算例 用泰勒展开将函数 在点 简化成线性函数与二次函数。 解:函数在点 的函数值、梯度和二阶导数矩阵:

算例 简化的线性函数 简化的二次函数

无约束优化问题的极值条件 1.F(x)在 处取得极值,其必要条件是: 即在极值点处函数的梯度为n维零向量。 函数的梯度为零的条件仅为必要的,而不是充分的。 例: 在 处梯度为 但 只是双曲抛物面的鞍点,而不是极小点。 则称 为 f 的驻点。 定义:设 是D的内点,若

Rastrigin function 为了判断从上述必要条件求得的 是否是极值点,需建立极值的充分条件。 为了判断从上述必要条件求得的 是否是极值点,需建立极值的充分条件。 根据函数在 点处的泰勒展开式,考虑上述极值必要条件,可得相应的充分条件。

Rastrigin function 处取得极值充分条件 Hessian矩阵 正定,即各阶主子式均大于零,则X*为极小点。

算例 求函数 的极值。 解:由极值的必要条件 由此知X1是函数的极小值点, X4是函数的极大值点,X2和X3均为非极值点。

等式约束优化问题的极值条件 求解等式约束优化问题为: min f(x) s.t. hk(x)=0 (k=1,2,…,m) 对这一问题,在数学上有两种处理方法: 消元法(降维法) 拉格朗日乘子法(升维法)

等式约束优化问题的极值条件 消元法 对n维函数的优化问题: min f(x1,x2,…,xn) s.t. hk(x1,x2,…,xn)=0 (k=1,2,…,l) 由l个约束方程将n个变量中的前l个变量用其余n-l个变量表示,即有 将这些函数关系代入到目标函数中,从而得到只xl+1,xl+2,…,xl+n共有n-l个变量的函数(n-l维),这样就可以利用无约束优化问题的极值条件求解,此法称为消元法或降维法。

等式约束优化问题的极值条件 拉格朗日乘子法 对于具有l个等式约束的n维优化问题 minf(x1,x2,…,xn) s.t. hk(x1,x2,…,xn)=0 (k=1,2,…,l) 把原目标函数f(x)改造成为如下形式的新的目标函数 式中的hk(x)就是原目标函数f(x)的等式约束条件,而待定系数k称为拉格朗日乘子。这种方法称为拉格朗日乘子法。 在极值点处,有 共有n+l个方程,足以解出这n+l个变量,此法也称为升维法。

算例 用拉格朗日乘子法计算在约束条件 h(x1,x2)=2x1+3x2-6=0的情况下,目标函数f(x1,x2)=4x12+5x22的极值点坐标。 解 改造的目标函数是 则 联立得极值点x*(1.071,1.286)

不等式约束优化问题的极值条件

不等式约束优化问题的极值条件

不等式约束优化问题的极值条件

不等式约束优化问题的极值条件 对于一元函数在给定区间[a, b]上的极值条件,上式中的第一式为 分析极值点x*在区间[a, b]中所处的位置,将会出现三种可能的情况 当a<x*<b时,μ1= μ2 =0,则极值条件为 2)当x*=a时, μ1≥0 , μ2=0 ,则极值条件为 3)当x*=b时, μ1 =0 , μ2 ≥ 0 ,则极值条件为

不等式约束优化问题的极值条件

库恩-塔克条件 要么目标函数的负梯度等于起作用约束梯度的非负线性组合(此时 )。 不等式约束的多元函数极值的必要条件是著名的 库恩—塔克(Kuhn-Tucker)条件,它是非线性优化问 题的重要理论 (1)库恩—塔克条件 (K-T条件) 对于多元函数不等式的约束优化问题: 库恩—塔克条件表明:如点 是函数 的极值点,要么 (此时 ) 要么目标函数的负梯度等于起作用约束梯度的非负线性组合(此时 )。

库恩-塔克条件 起作用约束: 极值点处于等值线的中心 x2 x1 x2 x1 极值点处于两个约束曲线的交点上 g1 (x)=0 x﹡ O x1 x2 x﹡ g1 (x)=0 g2 (x)=0 g3 (x)=0 起作用约束: 极值点处于等值线的中心 O x1 x2 x﹡ g1(x)=0 g2(x)=0 极值点处于两个约束曲线的交点上

库恩-塔克条件 库恩—塔克条件的几何意义是: 在约束极小值点 处,函数 的负梯度一定能表示成所有起作用约束在该点梯度(法向量)的非负线性组合。

库恩-塔克条件 同时具有等式和不等式约束的优化问题 : K-T条件:

库恩-塔克条件 K-T条件是多元函数取得约束极值的必要条件,以用来作为约束极值的判断条件,又可以来直接求解较简单的约束优化问题。 对于目标函数和约束函数都是凸函数的情况, 符合K-T条件的点一定是全局最优点。这种情况K-T条件即为多元函数取得约束极值的充分必要条件。

库恩-塔克条件 库恩—塔克(K-T)条件应用举例 s.t 判断[1 0]T是否为约束最优点。 (1)当前点 为可行点,因满足约束条件

库恩-塔克条件 (2)在 起作用约束为g1和g2 , 因 (3) 各函数的梯度:

库恩-塔克条件 (4)求拉格朗日乘子 由于拉格朗日乘子均为非负,说明 是一个局部最优点,因为它满足K-T条件。

库恩-塔克条件 s.t.