用户名: 密码: 验证码:
无约束优化问题的修正拟牛顿非单调信赖域算法研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
信赖域方法是非线性优化的一类重要的数值计算方法.它在近二十年来受到非线性优化领域许多研究者的关注,是非线性优化的研究热点.与线搜索相比,信赖域有两个突出的优点:一是它有很强的稳定性和强适性,二是它具有很强的收敛性.由于信赖域的有界性,它可以处理非凸的近似模型.目前,信赖域方法已经和传统的线搜索方法并列为求解非线性规划问题的两类主要数值方法[1],
     与线性搜索方法相比,信赖域算法不仅具有很强的收敛性[2],而且对于病态问题也能有效地解决,需要的迭代次数少,但由于求解子问题花费代价高,往往不易求解新的迭代点;而线性搜索方法易于求得新的迭代点.为充分发挥两种方法的优势,1991年,Jorge Nocedal和袁亚湘[3]提出将信赖域算法和线性搜索方法相结合来构造新计算方法的思想,在文[25]中,采用回溯(backtracking)线搜索,优点是不需重解子问题,大大减少了计算量,但为了保证序列{B_k}的正定性,却使得一些B_k未能满足拟牛顿方程,这样做往往使B_k逼近(?)~2f(x_k)的效果不佳,从而信赖域子问题不能很好地逼近原问题.在文[5]中E.MichaelGertz提出了一种新的带线搜索的信赖域方法,它不仅继承了文[25]中方法不需重解子问题的优点,而且由于在每步都采用Wolfe线搜索,使得序列{B_k}满足拟牛顿方程且保证其正定性,充分开发了拟牛顿校正公式的性质,克服了文[25]中方法的缺点.
     本文主要研究应用修正拟牛顿方程的非单调信赖域方法.因为,在实际计算中,对于某些问题单调算法并不能保证算法的有效性. 1986年, Grippo等人[6,7]提出了一种非单调线搜索,并将此技术分别运用到Newton法和截Newton法中1993年,邓乃扬等人[8]首次将非单调技术应用到信赖域方法中,在一定条件下证明了其全局收敛性和超线性收敛性,数值试验表明对某些问题,非单调信赖域方法比相应的单调算法有更好的数值结果.以上提到的非单调技术都是以(?)为参考函数值来实现的,其中m(0)=0, 0≤m(k)≤min[m(k-1)+1,M], M是给定的正整数.鉴于以上工作的基础上,本文提出了两类新的非单调信赖域方法.
     第一章,我们首先介绍了最优化问题的研究背景和现状,以及信赖域子问题的求解方法.其次回顾求解无约束最优化问题的主要非精确线性搜索方法.
     第二章,我们提出了一类新拟牛顿非单调信赖域方法.采用加权的r_k用以调整信赖域半径,在适当的条件下,证明了此算法的全局收敛性.数值结果表明该算法的有效性.
     第三章,我们提出了一类带线搜索的修正拟牛顿非单调信赖域算法.不同于传统的非单调信赖域算法,此算法在每步都采用非单调Wolfe线搜索得到下一个迭代点.这样得到的新算法不仅不需要重复求解子问题,而且由于应用基于修正拟牛顿方程的校正,可以得到逼近Hessian阵精度更高的B_k.在适当的条件下,证明了该算法的全局收敛性.数值试验证实该算法是有效的.
Trust region method is one of the important numerical computational methods in conlinear optimization. It has been paid so much attention by many researchers in the field of nonlinear optimization. Compared with line search method, TR method has two outstanding advantages: one is it has very strong stability and robustness; the other is its strong convergence. Because of the boundness of trust region. it can deal with nonconvex models. So far. TR method enlists in one of the two major numerical options for solving nonlinear optimization [1].
     Compared with line search method, TR method not only has strong convergence [2], but can solve bad-scaled problems efficiently with less iteration numbers. However, because it is difficult to solve the subproblem. the new iterative point is hard to obtain by using TR method, while line search method is easier to obtain the new point. As a result, in order to make good use of the advantages of both methods. Jorge Nocedal and Yuan Yaxiang [3] proposed an idea of combining TR method and line search method in 1991. In paper [25]. with backtracking line search and no need to resolve the subprobleni of TR method. the computational amount is greatly reduced. But in order to maitain the positive definiteness of B_k, causing some B_k do not perfectly match (?)~2f(x_k). and then the subproblem of TR method does not match the original problem accordingly. In paper [5], E.Micheal Gertz propesed a new TR method with line search method. It not only no need to resolve the subproblem, but since using Wolfe line search at every iteration, B_k satisfies the quasi-Newton equation and maitains the positive definiteness. It fully developped the propertises of quasi-Newton formular and avoided the disadvantages of paper [25].
     In this paper, we mainly study the nonmonotonic TR method based on modified quasi-Newton equations. Since in real computations, the monotonic algorithms can not guarantee the validity of some problems, In 1986. Grippo,etc [6,7] lay out a nonmonotone line search technique, and applied it to Newton method and cutted Newton method. In 1993, Deng Naiyang and others [8] firstly applied nonmonotonic technique to trust region method and proved its global and superlinear convergence under suitable conditions. Numerical experiments showed that the nonmonotonic trust region method was superior to the corresponding monotonic method for some problems. The nonmonotonic techniquementioned above is realized by f_(?) = (?) f(x_(k-j)), where m(0) = 0,0≤m(k)≤min[m(k-1)+1,M],M is the given positive integer. Based on the proceeding work, we propese two kinds of new nonmonotonic TR algorithms.
     In chapter 1, we firstly give a brief introduction of the background and current situation of Optimization, and methods to solve TR subproblem. THen, we review the major nonexact line search methods for unconstrained optimization.
     In chapter 2, we propose a new quasi-Newton nonmonotonictrust region algorithm. We adjust the trust radius using not only (?). but also the previous ratios {r_(k-m)…,r_k}, where m is some positive integer. Under proper assumptions. we prove the global convergence of the algorithm and numerical experiments show the algorithm is competitive.
     In chapter 3, we present a modified quasi-Newton nonmonotonic trust region algorithm with Wolfe line search. Unlike traditonal nonmonotonic trust region algorithms, our algorithm gets the next point by the nonmonotonic Wolfe line search at each iteration. This new algorithm not only does not resolve the sub-problem but also satisfies the modified quasi-Newton condition at each iteration and simultaneously maintains a positive-definite approximation to the Hessian of the objective function. Under mild conditions, we prove the global convergence and numerical results show its efficiency.
引文
[1] 席少霖.非线性最优化方法[M].北京:高等教育出版社, 1992
    [2] 袁亚湘,孙文瑜.最优化理论与方法[M].北京:科学出版社, 1997
    [3] Jorge Nocedal, Ya-xiang Yuan. Combining trust region and line search tech-niques[R]. Technical Reprot, NAM06, Dept. of Computer Science, Northwestern University, Illinois, USA. 1991.
    [4] Jorge Nocedal, Ya-xiang Yuan. Combining Trust Region and Line Search Tech-niques[J]. Advances in Nonlinear Programming, 1998, 153-175.
    [5] E.Michael Gertz. A quasi-Newton trust-region method[J]. Math.Program, 200-1. Ser.A.
    [6] L.Grippo, F.lampariello, S.Lucidi. A nonmonotone line search technique for New-ton's method[J]. SIAM J.Numer.Anal., 1986, 23(4): 707-716.
    [7] L.Grippo, F.Lampariello, S.Lucidi. A truncated Newton method with nonmono-tone line search for unconstrained optimization [J]. J Optim Theory Anal, 1989, 60: 401-419.
    [8] Deng N Y, Xiao Y, Zhou F J. A nonmonotonic trust region algorithm [J]. JOTA, 1993, 26: 259-285.
    [9] J.J.More. The Levenberg-Marquardt algorithm: implementation and theory. in: G.A.Watson, ed.. Lecture Notes in Mathematics 630: Munerical Analysis. Springer-Verlag, Berlin. 1978. pp. 105-116.
    [10]M.J.D.Powell.A new algorithm for unconstrained optimization.In J.B.Rosen,O.L.Mangasarian,and K.Ritter[M],editors,Nonlinear Programming.AcademicPress,1970a.
    [11]M.J.D.Powell.A hybrid method for nonlinear equations.In P.Rabinowitz,editor,Numerical Methods for Nonlinear Algebrac Equations.Gordon and Breach,1970b.
    [12]J.E.Dennis,H.H.W.Mei.Two new unconstrained optimization algorithms which use function and gradient values[J].JOTA,1979,28:452-482.
    [13]R.H.Byrd,R.B.Schnabel,G.A.Schultz.Approximate solution of the trust regions problem by minimization over two-dimensional subspaces[J].Mathematical Programming,1988,40:247-263.
    [14]G.A.Schultz,R.B.Schnabel.R.H.Byrd.A family of trust-region-based algorithms for unconstrained minimization with strong global convergence properties[J].SIAM Journal on Numerical Analysis.1985,22:47-67.
    [15]Zhang J Z,Xu C X.A class of indefinite dogleg path methods for unconstrained minimization[J].SIAM J.on Optimization,1994,V.9.P.646-667.
    [16]T.Steihaug.The conjugate gradient method and trust regions in large scale optimization[J].SIAM Journal on Numerical Analysis,1983,20:626-637.
    [17]Yuan Y X.On the truncated conjugate gradient method[R].Reprot,ICM99-003.ICMSEC,Chinese Academy of Sciences.Beijing,China,1999.
    [18]M.D.Hebden.An algorithm for minimization using exact second derivatives[R].Atomic Energy Research Establishment Report TP515,Harwell,England.
    [19]D.M.Gay.Computing optimal local constrained step[J].SIAM J.Sci.Stat.Comp.,1981,2:186-197.
    [20]D.C.Sorensen.Newton's method with a model trust region modification.SIAM J.Numer.Anal.,1982,20:409-426.
    [21] J.J.More, Burton S. Garbow, Kenneth E. Hillstrom. Testing unconstrained opti-mization software[J]. ACM Trans.Math.Software, 1981, 7: 17-41.
    [22] Jorge Nocedal, Stephen J. Wright. Numerical Optimization[M]. Science Press, China, 2006.
    [23] E.Michael.Gertz. Combination trust-region line-search methods for uncon-strained optimization [J]. University of California San Diego, 1999.
    [24] Ph.L.Toint. Towards an efficient sparsity exploiting Newton method for mini-mization. Sparse Martrices and Their Uses[M]. Academic Press, New Yord, 1882. I.S. Duff 57-87.
    [25] J.Nocedal, Y.Yuan. Combining trust region and line search techniques[J]. Opti-mization Technology Center mar OTC. 1998, 04.
    [26] P.E.Gill, M.H.Wright. Course Notes for numerical Nonlinear Optimization. De-partment of Mathematics, University of Carlifornia, San Diego, 2001.
    [27] Fan J Y, Yuan Y X. A new trust region algorithm with trust region radius converging to zero[J]. Chinese Academy of Science. 2004.
    [28] M.J.D.Powell. On trust region methods for unconstrained minimization without derivatives [J]. In math. programming, Ser.B Springer verlag, 2003, 97: 605-623.
    [29] 李正锋,邓乃扬.一类新的非单调信赖域算法及其收敛性[J].应用数学学报,1999,22(3):290-294
    [30] P.Toint. A nonmonotone trust region algorithm for nonlinear optimization sub-ject to convex constraints [J]. Math. Prog., 1997, 77(1): 69-94.
    [31] 柯小伍,韩继业.一类新的信赖域算法的全局收敛性[J].应用数学学报,1995,18(4):608-615
    [32] Wei Z X, Li G Y, Qi L Q. New quasi-Newton Methods for unconstrained optimiza-tion problems[J]. Applied Mathematics and Computation, 2006, 175: 1156-1188.
    [33] Dai Y H, Xu D C. A new family of trust region algorithms for unconstrained optimization[J]. Journal of Computational Mathematics, 2003, 21(2): 221-228.
    [34] Ke X W, Han J Y. A class of nonmonotone trust region algorithms for un-constrained optimization(in Chinese)[J]. Science in China(Series A), 1998, 28, 488-492.
    [35] Li Z F, Deng N Y. A new family of nonmonotone trust region algorithms and its properties(in Chinese)[J]. ACTA Mathematicae Applicatae Sinica, 1999, 22, 457-465.
    [36] Ke X W, Liu G H, Xu D C. A nonmonotone trust region algorithm for uncon-strained optimization [J]. Chinese Science Bulletin, 1996, 41, 197-201.
    [37] Zhang J Z, Deng N Y, Chen J H. New quasi-Newton equation and related meth-ods for unconstrained optimization[J]. Journal of Optimization Theory and Ap-plication, 1999, 102(1): 147-167.
    [38] Yuan Y X. Convergence properties of trust region method[J]. Mathematica Nu-merica Sinica, 1994, 16, 333-346.
    [39] Li D H, M. Fukushima. A modified BFGS method and its global convergence in nonconvex minimization [J]. Journal of Computational and Applied Mathematics. 2001, 129, 15-35.
    [40] D.Li, M. Fukushima. On the global convergence of the BFGS method for non-convex unconstrained optimization problems[J]. SIAM Journal on Optimization,2001,11, 1054-1064.
    [41] D.Li, L.Qi. BFGS-trust region method for minimizations[R]. Technical Reprot. Department of Applied Mathematics, Hongkong Polytechnic University. July,2002, 1-12.
    [42] G.Yu, Z.Wei, L.Guan. A modified BFGS-trust region method for nonconvex min-imization[J]. Pacific Journal of Optimization, 2006,2,119-133.
    [43] 朱德通.非单调投影梯度方法解不等式约束优化[J].数学年刊(A).1998.19(6):749-760.
    [44] W.Sun. Nonmonotone trust region method for solving optimization problems[J]. Applied Mathematics and Computation, 2004, 156, 159-174.
    [45] J.Fu, W.Sun. Nonmonotone adaptive trust-region method for unconstrained op-timization problems[J]. Applied Mathematics and Computation, 2005, 163, 489-504.
    [46] 刘光辉,彭积明.一类非单调算法的收敛性质[J].计算数学,1994,16(1):65-71.
    [47] J.Z.Zhang, C.X.Xu. Properties and numerical performance of quasi-Newton methods with modified quasi-Newton equations[J]. Journal of Computational and Applied Mathematics, 2001, 137, 269-278.
    [48] 王岩.应用新拟牛顿方程的信赖域方法.南京理工大学学报,2007.
    [49] 刘培培.无约束最优化问题的一类非单调信赖域算法研究.首都师范大学学报,2007.
    [50] 刘培培,陈兰平.一类拟牛顿非单调信赖域算法及其收敛性[J].数学进展,2008.37(1):92-100.
    [51] 姚升保,施保昌,彭积明.一类带线搜索的非单调信赖域方法[J].数学杂志,2003.23(3):290-294.
    [52] 莫降涛,刘春燕,颜世翠.带有固定步长的非单调信赖域算法[J].曲阜师大学报,2006.32(3):30-34.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700