数模与最优化.pptVIP

设函数f(x)是一个单峰函数，其初始区间为[a0,b0]，在此区间内任取两点x1与x2 ， 且 x1x2 ，计算函数值f(x1)与f(x2)，并比较，则处于以下三种情况 f(x1)f(x2)，构成新区间[a0, x2] f(x1)f(x2) ，构成新区间[x1,b0] f(x1)=f(x2) ，构成新区间[a0, x2] 或[x1,b0] 因此，只需在搜索区间内取两点，计算它们的函数值并加以比较后，总可以将搜索区间缩小，这就是消去法的原理 存在问题： * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.4 平分法 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.4 平分法 迭代步骤 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.4 平分法 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.5 黄金分割法 2.5.1　适用范围与基本思想 1）适用范围 2）基本思想 2.5.2　基本原理 如何插入新点才是科学合理的？ 插入新点原则 所插入两点在搜索区间内应处于对称位置 所选内点应使区间缩短率恒定 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.5 黄金分割法 【小结】只要第一个点取在原始区间的0.618处，第二个点在它的对称位置，这就保证无论经过多少次舍去，保留的点始终在新区间的0.168处。要进一步缩短区间，在其保留点的对称位置再作一次计算就可以了。而且每次消去时，区间的缩短率不变。 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.5 黄金分割法 2.5.3　具体计算方法及步骤 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.5 黄金分割法 2.5.3　具体计算方法及步骤 算法框图 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.6 牛顿法（函数逼近法） 基本思想 牛顿法（又称为切线法）是一种函数逼近法（还有抛物线法等） 在迭代点附近用二阶泰勒展式近似目标函数，进而求出极值点的估计值 基本原理 利用牛顿迭代公式计算目标函数一阶导数的驻点 牛顿迭代法解方程 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.6 牛顿法（函数逼近法） * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 2.6 牛顿法（函数逼近法） 牛顿法寻优就是利用牛顿迭代解方程方法求解目标函数的一阶偏导数等于零的解，即驻点，此时将目标函数的一阶偏导数方程看成一个一般方程，即可利用上述迭代法求解 牛顿法迭代方程 步骤简单，收敛快 需计算一、二阶导数，计算量大 初始点选择不合适可能导致发散 优点 缺点 函数本身比较复杂（导数有弯曲） 曲线 上凹选a作为初始点失效 曲线 下凹选b作为初始点失效 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 3 无约束条件下多变量寻优 3.1 最速下降法(间接法) 3.2 牛顿法（间接法） 共轭梯度法、变尺度法（DFP法，最为有效方法） 3.3 变量轮换法（直接法） 3.4 单纯形替换法（直接法） * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 3.1 最速下降法 下降迭代算法包含两个关键步骤 从一个迭代点如何选择搜索方向？（本节将来介绍） 在搜索方向上如何进一步进行一维搜索？（上面几节已经介绍） 1847年 柯西Cauchy提出，梯度法 基本思想 每次沿目标函数在迭代点处函数值下降最快的方向（梯度反方向）进行一维搜索，逐步达到极小点。 步骤 (见黑板) * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 3.1 最速下降法 收敛速度 锯齿现象 * 西安建筑科技大学 冶金工程学院 金属材料加工研究所 压力加工数模与优化 * 3.2 牛顿法 基本思想 在目标函数f(x)具有连续二阶偏导数的条件下，用一个二次函数Q(x)去近似目标函数，然后求出这个二次函数的极小点，做为f(x)极小点的近似值 迭代公式推导 设X(k)是f(X)极小点X*的第k次近似，将f(X)在X(K)点