基于MATLAB的鲍威尔法求极值问题.doc

基于MATLAB的鲍威尔法求极值问题.doc

《基于MATLAB的鲍威尔法求极值问题.doc》由会员分享，可在线阅读，更多相关《基于MATLAB的鲍威尔法求极值问题.doc（13页珍藏版）》请在冰豆网上搜索。

机械优化设计报告（4）

基于MATLAB的鲍威尔法求极值问题

姓名：

xxx 学号：

xxx

（北京理工大学机械与车辆学院车辆工程，北京100081）

摘要：

无约束优化方法主要有七种，按照求导与否把这些方法分为间接法和直接法。

牛顿法的成败与初始点选择有极大关系，其可靠性最差；坐标轮换法、单纯形法和最速下降法对于高维优化问题计算效率很低，有效性差；由于编制变尺度法程序复杂，其简便性不足。

综合考虑后，鲍威尔法、共轭梯度法具有较好的综合性能。

本文首先对鲍威尔法的原理进行阐述，根据其迭代过程给出流程图，并编写MATLAB程序。

最后用此MATLAB程序求解实际的极值问题，并对求解结果进行简要分析。

1.鲍威尔法的基本思想

1.1其他优化方法对鲍威尔法形成的影响

通过对鲍威尔法的学习，可以很明显看出来其迭代思想中汲取了其他几种优化方法的核心思想。

为了更全面、更深入的学习鲍威尔法，很有必要对其他有影响的优化思想进行学习和梳理。

由最基本的数学基础知识可知，梯度方向是函数增加最快的方向，负梯度方向是函数下降最快的方向，于是，利用这个下降最快方向产生了最速下降法。

每次迭代都沿着负梯度方向进行一维搜索，直到满足精度要求为止。

其特点是相邻两个搜索方向互相正交，所以很明显的一个现象就是刚开始搜索步长比较大，愈靠近极值点其步长愈小，收敛速度愈慢，特别当二维二次目标函数的等值线是较扁的椭圆时，迭代速度更慢。

这时，倘若目标函数是等值线长、短轴都平行于坐标轴的椭圆形，则通过坐标轮换法可以很高效的解决问题。

通过两次分别沿坐标轴进行一维搜索，便可达到极值点。

但对于目标函数的等值线椭圆的长、短轴倾斜于坐标轴时，坐标轮换法的搜索效率也显得极低。

抛开这两种特殊情况，对于一般形态的目标函数，如果在某些明显可以直达最优点的情况下（一般为靠近极值点区域），迭代过程完全可以不沿负梯度方向搜索，取而代之的是找到直达最优点的方向，一步到位。

但这样的直达方向应该如何去找呢？

共轭梯度法由此产生。

其基本原理是：

任意形式的目标函数在极值点附近的特性都近似一个二次函数，其等值线在极值点附近为近似的同心椭圆簇，而同心椭圆簇有一个特性便是任意两条平行线与椭圆簇切点的连线必通过椭圆的中心。

而这个连线方向便是所寻找的直达方向。

通过对其迭代过程的分析，很明显可以看出需大量的求目标函数的一阶和二阶偏导数。

对于一些实际的机械优化问题，目标函数可能复杂到难以求取其偏导数或者根本就不存在，求取极值问题就显得无从下手。

所以，一个效率高的、适应性强的优化方法急需出现，而鲍威尔法便是这么一种综合的方法。

1964年，鲍威尔提出了对共轭方向法的改进方法——鲍威尔共轭方向法。

一维搜索法、共轭方向和坐标轮换法的思想在鲍威尔法中体现的淋漓尽致。

下面就对鲍威尔法的基本原理进行阐述。

1.2鲍威尔法的数学原理

通过前文可知，鲍威尔法也算一种共轭方向法，但与共轭梯度法相比，不需要对函数求导，而是在迭代过程中逐次构造出用于搜索的共轭方法。

1）.对于二维无约束优化问题，采用鲍威尔法求解的迭代过程如图1-1所示。

任选一初始点,令，按照坐标轮换法，选择两个单位向量和，以此作为搜索方向进行第一轮搜索得到点。

2）.用和的连线方向构成新的搜索方向。

从点出发，沿方向一维搜索得到点，作为下一轮搜索的初始点。

3）.从出发，依次沿和方向进行一维搜索，得到点和点。

4）.用点和点的连线方向构成新的搜索方向。

和是从两个不同点出发沿相同方向搜索得到的，所以与和的连线方向互为共轭方向。

从点出发，沿方向一维搜索得到点。

因是从点出发依次沿两个互为共轭的方向和进行两次一维搜索得到的，所以就是该二维二次函数的极小点。

图1-1二维情况下鲍威尔法的迭代过程

将上述二维优化问题扩展到n维的情况，得到鲍威尔法的基本迭代过程：

从初始点出发依次沿n个线性无关的方向组进行一维搜索得到一个终点，沿初始点和终点的连线方向一维搜索得到下一轮迭代的初始点，并以这个方向作为下一轮迭代方向组中的最后一个方向，同时去掉第一个方向，组成的新方向组进行第一轮迭代。

若目标函数是个n维的正定二次函数，则经过这样的n轮迭代以后，就可以收敛到最优点。

但是这种方法有一个缺陷，通过这种方法产生的n个新方向，有可能是线性相关或近似线性相关的。

因为新方向，如果其中出现了，则就可以表示为，，…，的线性组合，这样用新方向替换后，新的坐标轮换方向组就成为线性相关的一组向量，以后的各轮迭代计算将在维数下降了的空间内进行，这将导致算法收敛不到真正的最优点。

针对此现象，通过适当改进，产生了新的鲍威尔法。

1.3改进的鲍威尔法

由上一小节讨论可知，如果对新方向组的形成加以适当的选择，防止其线性相关，则可以避免鲍威尔法的“退化”。

在这里直接给出是否用方向替代来组成新的搜索方向组的判别条件，供后面编程所用。

设，，。

其中最大者。

若且成立，则用方向替代方向，否则仍用原来的n个搜索方向。

这样，改进后的鲍威尔法保证了对于非线性函数计算时收敛的可靠性。

2.迭代过程和算法流程图

2.1迭代过程

1）.给定初值点，n个线性无关的初始向量组及精度，置。

2）.从点出发沿中的方向进行一维搜索：

，其中

3）.判断是否满足收敛准则，若

成立，则迭代停止，输出，即为最优点；否则跳至第4步。

4）.计算第k轮迭代中相邻两点目标函数值的下降量，并找出下降量最大者。

相应的方向：

5）.沿共轭方向计算反射点。

6）.检验判别条件。

令，，，若同时满足：

则转第6步；否则转第7步。

7）.由出发沿方向进行一维搜索，求出此方向极小点，并令下一轮迭代初始点。

同时令，同时置。

并转第2步。

8）.若，则；否则。

置，转第2步。

2.2算法流程图

鲍威尔法的算法流程图如图2-1。

图2-1鲍威尔法的算法流程图

3.鲍威尔法的MATLAB程序

function[x,minf,k,allX,allY]=Powell_Graphic（f,x0,D,var,eps）

formatlong;

ifnargin==4%默认精度为1.0e-2

eps=1.0e-2;

end

n=length（var）+1;

k=0;

symsl;

while1

tx=zeros（size（D））;

tx（:

1）=x0;

fori=1:

n-1%在每个方向上开始一维搜索

xv=tx（:

i）+l*D（:

i）;

fx=Funval（f,var,xv）;

[a,b]=minJT（fx,0,0.1）;%一维搜索进退法确定搜索区间

tl=minHJ（fx,a,b）;%黄金分割法求最优点

tx（:

i+1）=tx（:

i）+tl*D（:

i）;

end

D（:

n）=tx（:

n）-tx（:

1）;%判断收敛性是否满足精度要求

ifnorm（D（:

n））<=eps

x=tx（:

n）;

allX（:

k*n+1:

（k+1）*n）=tx（:

1:

n）;

forj=1:

n

allY（:

k*n+j）=Funval（f,var,tx（:

j））;

end

break;

else

forj=1:

n

FX（j）=Funval（f,var,tx（:

j））;

end

maxDF=-inf;

m=0;

forj=1:

n-1

df=FX（j）-FX（j+1）;

ifdf>maxDF

maxDF=df;

m=j+1;

end

end

tmpF=Funval（f,var,2*tx（:

n）-tx（:

1））;

fl=（FX

（1）-2*FX（n-1）+tmpF）*（FX

（1）-FX（n-1）-maxDF）^2;

iffl<0.5*maxDF*（FX

（1）-tmpF）^2&&tmpF

（1）%检验判断条件决定换方向

xv=tx（:

n）+l*D（:

n）;

fx=Funval（f,var,xv）;

[a,b]=minJT（fx,0,0.1）;

tl=minHJ（fx,a,b）;

x0=tx（:

n）+tl*D（:

n）;

D（:

m:

（n-1））=D（:

（m+1）:

n）;

else

iftmpF

x0=2*tx（:

n）-tx（:

1）;

else

x0=tx（:

n）;

end

end

k=k+1;

allX（:

（k-1）*n+1:

k*n）=tx（:

1:

n）;

forj=1:

n

allY（:

（k-1）*n+j）=Funval（f,var,tx（:

j））;

end

end

end

minf=Funval（f,var,x）;

formatshort;

程序中用到的函数minJT、minHJ和Funval分别为进退法确定搜索区间函数、黄金分割法求极值函数和计算函数值函数，直接调用即可，由于篇幅限制，故不在此列出其代码。

4.利用MATLAB求解实例

4.1实例

本文以本章例题4.9为例来验证鲍威尔法算法在MATLAB里的实现。

已知目标函数，求其极小值，计算精度。

在MATLAB命令窗口里输入如下命令：

>>symsst;

>>f=t^2+s^2-t*s-10*t-4*s+60;

>>D=[100;010];

>>[x,mf,k]=Powell（f,[00],D,[ts]）

计算结果为：

x=

8.0000

6.0000

mf=

8.0000

k=

2

与书中提供的精确解完全吻合，说明此方法有效可行。

计算变量如图3-1。

图3-1计算结果（）

4.2实例结果分析

此处主要从以下三个方面来分析此算法及计算结果。

1）可靠性。

通过对其他多元多次函数进行求极值计算，结果都比较精确，说明此算法可靠性、通用性好。

另外，对于给定不同的初始点，其求解结果同样精确。

分别如下：

图3-2[-10-10]

图3-3[00]

图3-4[1010]

可见，不论给予什么样的迭代初始点，其最终迭代结果同样正确且精确。

2）有效性。

对于