APP下载

一种改进的BFGS算法及其收敛性分析

2011-07-06陈奎林

重庆理工大学学报(自然科学) 2011年11期

陈奎林

(重庆大学数学与统计学院,重庆 401331)

1 问题的提出

本文针对无约束最优化问题:minf(x),x∈Rn开展研究,其中f:Rn→R是一个连续可微的函数。拟牛顿算法中的BFGS算法是一类解决无约束最优化问题的有效算法,这类算法的关键是Bk的选取,不同的选取方式,对应不同的BFGS算法。

为了使算法更具优越性,文献[1-5]在传统的拟牛顿方程的基础上提出了一个新的拟牛顿方程,即,其中是一个对称正定矩阵,从而得到了一类改进的BFGS方法,其迭代公式为

受以上思想的启发,提出了Ak的一种选取方法,即Ak=‖gk+1-gk‖I,进而提出了一类新的BFGS算法,并证明了它的全局收敛性和超线性收敛性。

本文采用的搜索准则(Wolfe准则)为:

其中0<σ1<σ2<1。

改进的BFGS算法步骤:

步骤1 给出初始点x0∈Rn和初始对称正定矩阵B0∈Rn×n,令ε>0,k=0。

步骤2 若梯度函数在迭代点xk处满足‖gk‖≤ε,则停止;否则计算Bkdk+gk=0,求出搜索方向dk。

步骤3 利用Wolfe准则求得步长αk,令xk+1=xk+αkdk。

步骤4 计算Ak=‖gk+1-gk‖I,并代入 式(1),修正Bk得Bk+1。

步骤5 令k:=k+1,转步骤2。

2 收敛性证明

为证明算法的全局收敛性和超线性收敛性需要以下假设:

①f(x)是二阶连续可微的,x*是f(x)的极小点。

②水平集Ω={x|f(x)≤f(x0)}是有界凸集,其中x0是给定的。

③f(x)在凸集Ω上连续可微,且存在一个常数L>0,使得下式成立:

其中:g(x)=▽f(x);‖·‖是Euclidean范数。

④f(x)一致凸,即存在常数m和M,使得

这里∀x∈Ω,z∈Rn,其中 G(x)是 f(x)的海色矩阵。

⑤ G(x)在Ω上Lipschitz连续。即存在L'>0使得∀x,y∈Ω,有

由上面的假设容易得到

2.1 全局收敛性证明

引理1[4]对任意给定的 k 和由改进的 BFGS 算法产生的(α,x,g,d),若>0,那么

kk+1k+1k+1Bk+1正定。

引理2 若假设(a)、(b)、(d)成立,由改进的BFGS算法产生的点列为 {xk},则,其中sk=xk+1-xk。

证明由Taylor公式,

由假设④及式(2)有

引理3 若假设②~④成立,则式(4)(5)成立

登录APP查看全文