无约束优化计算方法.ppt

第62页,共103页，星期六，2024年，5月第63页,共103页，星期六，2024年，5月第64页,共103页，星期六，2024年，5月第65页,共103页，星期六，2024年，5月第66页,共103页，星期六，2024年，5月鲍威尔方法是鲍威尔于1964年提出的，以后又经过他本人的改进。该法是一种有效的共轭方向法，它可以在有限步内找到二次函数的极小点。对于非二次函数只要具有连续二阶导数，用这种方法也是有效的。第67页,共103页，星期六，2024年，5月单纯形方法一、基本思想单纯形替换法也是一种不使用导数的求解无约束极小化问题的直接搜索方法，与前面几种方法不同的是，单纯形替换法不是利用搜索方向从一个点迭代到另一个更优的点，而是从一个单纯形迭代到另一个更优的单纯形。第68页,共103页，星期六，2024年，5月定义：单纯形n维空间中的恰好有n+1个顶点（极点）的有界的凸多面体称之为一个单纯形。根据定义，可知，一维空间中的单纯形是线段，二维空间中的单纯形是三角形，而三维空间中的单纯形则是四面体。在单纯形替换算法中，从一个单纯形到另一个单纯形的迭代主要通过反射、扩张、收缩和缩边这4个操作来实现。下面以二维问题为例来对4种操作进行说明（参见下图）。第69页,共103页，星期六，2024年，5月第70页,共103页，星期六，2024年，5月第71页,共103页，星期六，2024年，5月第72页,共103页，星期六，2024年，5月（1）反射——设除了最劣点X1以外的基余顶点的中心为X4，作X1关于点X4的对称点X5，称X5为X1的反射点。求反射点的过程称之为反射。（2）扩张——在得到反射点X5之后，如果X5优于原单纯形的最劣点（即有），表明反射方向（X5—X1）是有利方向，反射成功。若进一步有，可沿反射方向前进适当的距离到点X6。X6称之为扩张点，求扩张点的过程称之为扩张。扩张之后，若扩张点X6优于反射点X5，则扩张成功，以X6取代X1，得新单纯形{X6,X2,X3}；否则，扩张失败，舍弃扩张点，以反射X5点取代X1，得新单纯形{X5,X2,X3}。设当前的单纯形的顶点为X1，X2，X3，且有第73页,共103页，星期六，2024年，5月例:用黄金分割法求函数f(x)=3x3-4x+2的极小点

给定x0=0,h=1,ε=0.2。解：1）确定初始区间x1=x0=0,f1=f(x1)=2x2=x0+h=0+1=1,f2=f(x2)=1由于f1f2,应在原方向继续向前探测。x3=x2+h=1+1=2,f3=f(x3)=18由于f2f3,可知初始区间已经找到，即[a,b]=[x1,x2]=[0,2]2）用黄金分割法缩小区间第一次缩小区间：x1=0+0.382X(2-0)=0.764,f1=0.282x2=0+0.618X(2-0)=1.236,f2=2.72f1f2,新区间[a,b]=[a,x2]=[0,1.236],b-a0.2第30页,共103页，星期六，2024年，5月第二次缩小区间：令x2=x1=0.764,f2=f1=0.282x1=0+0.382X(1.236-0)=0.472,f1=0.317由于f1f2,故新区间[a,b]=[x1,b]=[0.472,1.236]因为b-a=1.236-0.472=0.7640.2,应继续缩小区间。第31页,共103页，星期六，2024年，5月第三次缩小区间：令x1=x2=0.764,f1=f2=0.282x2=0.472+0.618X(1.236-0.472)=0.944,f2=0.747由于f1f2,故新区间[a,b]=[a,x2]=[0.472,0.944]因为b-a=0.944-0.472=0.4720.2,应继续缩小区间。第32页,共103页，星期六，2024年，5月第四次缩小区间：令x2=x1=0.764,f2=f1=0.282x1=0.472+0.382X(0.944-0.472)=