Advances in Applied Mathematics
Vol.05 No.04(2016), Article ID:18977,10 pages
10.12677/AAM.2016.54072

Finite Convergence of the Solution Set for SQP Subproblem

Yajing Gu, Wenling Zhao*

School of Science, Shandong University of Technology, Zibo Shandong

Received: Oct. 28th, 2016; accepted: Nov. 10th, 2016; published: Nov. 18th, 2016

Copyright © 2016 by authors and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

ABSTRACT

Sequential Quadratic Programming (SQP) method is one of the most effective methods for solving constrained optimization problems. There is a class of subproblems during the process via SQP method. The subproblem which is called SQP multi-parameter subproblem is a multi-parametric quadratic programming. In this work, we introduce weak sharp solution set into SQP multi-pa- rameter subproblem. The character of weak sharp solution set is discussed. Under weak sharp conditions of solution set, the sufficient and necessary condition for finite convergence of feasible solution sequence via any algorithm is obtained.

Keywords:Subproblems of Sequential Quadratic Programming, Multi-Paramatic Programming, Feasible Solution Sequence, Weak Sharp Set, Finite Convergence

SQP子问题解集的有限收敛性

顾亚静,赵文玲*

山东理工大学理学院,山东 淄博

收稿日期:2016年10月28日;录用日期:2016年11月10日;发布日期:2016年11月18日

摘 要

序列二次规划方法(SQP)是求解约束优化问题的最有效的方法之一。SQP方法求解过程中产生的子问题是一个带参数的二次规划问题(SQP多参数规划子问题)。本文在SQP多参数规划子问题中,引入了其解集弱强的概念,讨论了弱强集的性质,并在其解集是弱强的条件下,给出了由任意算法所产生的可行解序列有限收敛的必要与充分条件。

关键词 :序列二次规划方法的子问题,多参数规划,可行解序列,弱强集,有限收敛

1. 引言

序列二次规划方法(SQP)是求解约束优化问题的一类最有效的方法。

考虑带约束的最优化问题(NP)为

其中是连续可微函数,是约束集。

求约束优化问题(NP)的SQP方法是一种迭代法,它的基本思想就是:将问题中的目标函数用一个二次模型近似,从而得到问题(NP)的一个近似二次规划。求解这个二次规划,便得到问题(NP)的一个近似解迭代点,再从这个点出发,重复以上步骤。这样通过求解一系列二次规划,便得到原问题的解的一串近似解点列。

是问题(NP)的当前的迭代点,通常情况下,下一步就是通过求解下列子问题来实现的。

其中,表示的梯度,为一对称正定矩阵。此问题是一个非线性的多参数规划问题。

多参数规划问题有着广泛的应用。仅在预测控制领域,就被广泛应用于炼油,石化,化工等生产过程中。目前线性多参数规划的发展相对来说比较成熟,国内外学者主要研究非线性的多参数规划。关于非线性多参数规划的研究,早期Fiacco [1] 考虑了Hilbert空间上的非线性多参规划,主要研究了最优性条件,含参最优解向量的局部灵敏性分析及方向导数的计算。在Fiacco的基础上,Grancharova [2] 考虑了空间上非线性多参规划,给出了典型的非线性规划算法SQP方法的KKT条件。M. Diehl [3] ,E. S. Mostafa [4] [5] 和V. Kungurtsev [6] [7] 进一步研究了SQP非线性多参规划,根据问题的性质给出了参数解。Diehl [3] 考虑了SQP在每一步迭代的子问题。V. Kungurtsev [6] 通过对经典SQP增加两次修正,给出了参数解,得到算法的局部收敛。

以上文献虽然利用一些方法对多参数规划问题求解,并得到了参数解,但是它们的结果都依赖于这些具体的方法,不能保证对任意算法所产生的迭代点列是收敛的,更不能保证是有限收敛的。本文将研究任意可行解序列的有限收敛性。

弱强性是优化问题扰动分析的一个重要工具。在数学规划和变分不等式问题中,解集的弱强性对迭代算法产生点列的收敛性起着重要作用。Burke and Ferris [8] 给出了数学规划的解集的弱强性的概念,并证明解集的弱强性是逼近点算法(见 [9] [10] [11] [12] )与某些重要的迭代算法(见 [13] [14] [15] [16] [17] )具有有限收敛性的充分条件。随后,Marcotte and Zhu [18] ,Xiu与Zhang [19] ,Zhou与Wang [12] 相继做出了推广与改进,并弱化条件,得到算法收敛的结果。

本文研究序列二次规划SQP多参数规划子问题mpSQP(x):

其中是闭凸集,,对上连续可微。

mpSQP(x)问题的最优解集记为

在此问题中,引进解集的弱强性的概念,并在其解集是弱强的条件下,给出了由任意可行解序列有限收敛的必要与充分条件。本文结构如下。在第2节中给出了一些预备知识。在第3节对mpSQP(x)的解集给出了弱强概念,讨论解集弱强的一些性质。在第4节中,当mpSQP(x)的解集是弱强集时,给出了可行解点列具有有限收敛性的充要条件。第5节为本文的总结。

2. 预备知识

在这一节里,我们将介绍本文用到的一些基本的概念。

是一个无限子序列,定义

据上述定义即知

在点的切锥定义为

在点的正则法锥定义为

在点一般意义下的法锥定义为

的极锥定义为

据 [20] 中的命题6.5,我们有

是凸集时,据 [20] 中的定理6.9,我们有

在闭集上的投影定义为

的距离定义为

如果是闭集,则有

设函数在点的次微分,则在点的投影次微分定义为

如果在点是可微的,据 [20] 的定理25.1知,即此时投影次微分即是投影梯度

我们称序列有限收敛于,如果存在,当时,有

3. mpSQP问题解集的弱强性

本节我们将给出mpSQP问题的弱强性,及与弱强性有关的一系列定理和推论。

考虑如下的序列二次规划SQP多参数子问题mpSQP(x):

其中,

是闭凸集,,对上连续可微。mpSQP(x)问题的最优解集记为

在给出关于mpSQP的弱强性定义之前,我们先来介绍一下对于约束问题,光滑的凸规划及非光滑的凸规划中有关弱强性的定义。

对带约束的优化问题:

在光滑的凸规划中,是正常闭凸的可微函数,为非空闭凸集,其解集上是具有模的弱强集,当且仅当对

在非光滑的凸规划中,是正常的闭凸函数,为非空闭凸集,其解集上是具有模的弱强集,当且仅当对

其中的次微分。

值得注意的是,光滑时,两个定义等价,此时上是常值向量(见 [21] ,推论6)。对于mpSQP问题我们得到相似的结果。基于mpSQP问题,下面给出其解集弱强的定义。

定义3.1 在mpSQP问题中,解集称为是具有模的弱强集,如果存在参数,使得对

。(3.1)

利用上面的定义,可以得到弱强集具有下面的性质。

定理3.1在mpSQP问题中,对,解集为具有模的弱强集的充要条件为

(3.2)

证明 (必要性)对,因为为mpSQP问题的弱强集,由(3.1)得到,

因此

(充分性)对,由(3.2)知

,单位球,满足

的任意性,

即(3.1)式成立。证毕。

定理3.2 若mpSQP问题的解集是模为的弱强集,当且仅当

。 (3.3)

证明 (必要性)

为弱强集,为单位球,则在点,使得

也就是

因此

,上式即为

得到

即(3.3)式成立。

(充分性)

由(3.3),对是单位球,存在满足

其中任意,得到(3.1)式。证毕。

定理3.3若mpSQP问题中的解集是参数为的弱强集,当且仅当对

。 (3.4)

证明 (必要性) 对于,且。解集为模为的弱强极小集,连续可微,又由定理3.2必要性,

满足

即得到(3.4)式。

(充分性) 对,如果,则

显然有

否则,任取

,有

推出

由此可知解集是法向量为的超平面的子集。

令序列收敛到,对正数序列,满足。可以得到

模为,再由(3.4)式,得到

,有

因此,由定理3.2 充分性即得为弱强集,弱强参数即为。证毕。

4. mpSQP问题可行解序列的有限收敛性

对于SQP子问题,文献 [21] 得到由任意算法所产生的可行解序列与最优解之间的距离有一个全局误差界。本文在其解集是弱强的条件下,给出了由任意算法所产生的可行解序列有限收敛到弱强集的充分必要条件。

定理4.1若mpSQP问题的解集是模为的弱强集,时,有限收敛到的充要条件为

(4.1)

证明 (必要性)对充分大的时,有,此时有

为凸集,由 [20] 中命题6.5和命题6.9知

因此有

(充分性) 反证法。假设存在序列,对,但是仍然有

是闭凸集,对任一给定的,存在满足

的弱强性定义,有

其中

,因此

那么

时,我们得到,矛盾。因此,(4.1)成立时,对收敛到。证毕。

关于变分不等式弱强性的文献 [12] [18] 中,迭代点有限收敛的充要条件为目标函数负梯度在切锥上的投影收敛到0。在目标函数可微的凸规划 [8] 中,迭代点有限收敛等价于目标函数负梯度在切锥上的投影收敛到0;在目标函数不可微的凸规划 [12] 中,迭代点有限收敛要求目标函数负梯度在切锥上的投影收敛到0。同样,本文中以参数形式的SQP子问题为目标函数,要求函数连续可微,去掉了 [8] [18] 中

这一条件,得到迭代点有限收敛的充要条件为目标函数负梯度在切锥上的投影收敛到0。

对于mpSQP中的,可写为的形式,又

考虑变分不等式问题

mpSQP问题就可以看做变分不等式问题的一个推广。就有如下推论:

推论4.1若VIP问题的解集是模为的弱强集,时,有限收敛到的充要条件为

推论4.1即为变分不等式 [18] 中定理5.2关于VIP问题有限收敛的结果,且去掉了不必要的条件

收敛到0。

5. 结束语

本文在多参数规划的环境中分析SQP子问题,通过定义解集的弱强性,将弱强的概念推广到多参数规划问题中,不拘泥于具体算法的约束,得到多参数规划最优解序列的有限收敛性。应该指出的是,在我们这个结果的推论中,包括了现有的相关文献中在解集是弱强条件下相应的结果。此外,我们建立的弱强的概念为许多最优化算法的有限收敛提供了更弱的充分条件。

最近,N. H. Xiu研究了可行集S为低秩集时三种不同的法锥对可行解序列的包含关系影响,下一步我们可以针对这三种不同定义的法锥做分别的讨论。另外,弱强集与误差界关系密切,可做进一步探讨。

基金项目

国家自然科学基金资助项目(No.11271233)和山东省自然科学基金资助项目(ZR2012AM016)资助。

文章引用

顾亚静,赵文玲. SQP子问题解集的有限收敛性
Finite Convergence of the Solution Set for SQP Subproblem[J]. 应用数学进展, 2016, 05(04): 620-629. http://dx.doi.org/10.12677/AAM.2016.54072

参考文献 (References)

  1. 1. Fiacco, A.V. (1983) Introduction to Sensitivity and Stability Analysis in Nonlinear Programming. Academic Press, New York.

  2. 2. Grancharova, A. and Johansen, T.A. (2012) Multi-Parametric Programming. Lecture Notes in Control and Information Sciences, 429, 1-37. http://dx.doi.org/10.1007/978-3-642-28780-0_1

  3. 3. Diehl, M., Uslu, I., et al. (2001) Real-Time Optimi-zation for Large Scale Processes: Nonlinear Model Predictive Control of a High Purity Distillation Column. Online Optimization of Large Scale Systems. Springer Berlin, Heidelberg, 363-383.

  4. 4. https://www.researchgate.net/publication/27279305_Real-Time_Optimization_for_Large_Scale_Processes_Nonlinear_Model_Predictive_Control_of_a_High_Purity_Distillation_Column

  5. 5. Mostafa, E.S. (2012) An SQP Trust Region Method for Solving the Discrete-Time Linear Quadratic Control Problem. International Journal of Applied Mathematics and Computer Science, 22, 353-363. http://dx.doi.org/10.2478/v10006-012-0026-5

  6. 6. Mostafa, E.S., Ismail, H.G. and Al-Afandi, N.F. (2013) An SQP Augmented Lagrangian Method for Two Classes of Nonlinear Semidefinite Programming Problems Arising in Discrete-Time Feedback Control. Pacific Journal of Optimization, 9, 511-534.

  7. 7. Kungurtsev, V. (2013) Second-Derivative Sequential Quadratic Programming Methods for Nonlinear Optimization. Dissertations and Theses, Gradworks.

  8. 8. Kungurtsev, V. and Diehl, M. (2014) SQP Methods for Parametric Nonlinear Optimization. http://www.optimization-online.org/DB_FILE/2014/02/4255.pdf

  9. 9. Burke, J.V. and Ferris, M.C. (1993) Weak Sharp Minima in Mathematical Programming. SIAM Journal on Control and Optimization, 31, 1340-1359. http://dx.doi.org/10.1137/0331063

  10. 10. Ferris, M.C. (1991) Finite Termination of the Proximal Point Algorithm. Mathematical Programming, 50, 359-366. http://dx.doi.org/10.1007/BF01594944

  11. 11. Matsushita, S.Y. and Xu, L. (2012) Finite Termination of the Proximal Point Al-gorithm in Banach Spaces. Journal of Mathematical Analysis and Applications, 387, 765-769. http://dx.doi.org/10.1016/j.jmaa.2011.09.032

  12. 12. Rockafellar, R.T. (1976) Monotone Operators and the Proximal Point Algo-rithm. SIAM Journal on Control and Optimization, 14, 877-898. http://dx.doi.org/10.1137/0314056

  13. 13. Zhou, J.C. and Wang, C.Y. (2012) New Characterizations of Weak Sharp Minima. Optimization Letters, 6, 1773-1785. http://dx.doi.org/10.1007/s11590-011-0369-0

  14. 14. Gupta, A., Bhartiya, S. and Nataraj, P.S.V. (2011) A Novel Approach to Multi-Parametric Quadratic Programming. Automatica, 47, 2112-2117. https://www.researchgate.net/publication/220158033_A_novel_approach_to_multiparametric_quadratic_programming

  15. 15. Burke, J.V. and More, J.J. (1988) On the Identification of Active Constraints. SIAM Journal on Numerical Analysis, 25, 1197-1211. http://dx.doi.org/10.1137/0725068

  16. 16. Wang, C.Y., Liu, Q. and Yang, X.M. (2005) Convergence Properties of Nonmonotone Spectral Projected Gradient Methods. Journal of Computational and Applied Mathematics, 182, 51-66. http://dx.doi.org/10.1016/j.cam.2004.10.018

  17. 17. Wang, C.Y., Zhao, W.L., Zhou, J.H. and Lian, S.J. (2013) Global Convergence and Finite Termination of a Class of Smooth Penalty Function Algorithms. Optimization Methods and Software, 28, 1-25. http://dx.doi.org/10.1080/10556788.2011.579965

  18. 18. Xiu, N.H. and Zhang, J.Z. (2003) Some Recent Advances in Projec-tion-Type Methods for Variational Inequalities. Journal of Computational and Applied Mathematics, 152, 559-585. http://dx.doi.org/10.1016/S0377-0427(02)00730-6

  19. 19. Marcotte, P. and Zhu, D.L. (1998) Weak Sharp Solutions of Variational Inequalities. SIAM Journal on Optimization, 9, 179-189. http://dx.doi.org/10.1137/S1052623496309867

  20. 20. Xiu, N.H. and Zhang, J.Z. (2005) On Finite Convergence of Proximal Point Algorithms for Variational Inequalities. Journal of Mathematical Analysis and Applications, 312, 148-158. http://dx.doi.org/10.1016/j.jmaa.2005.03.026

  21. 21. Rockafellar, R.T. and Wets, R.J. (1998) Variational Analysis. Springer, New York. http://dx.doi.org/10.1007/978-3-642-02431-3

  22. 22. Zhao, W.L. and Song, D.J. (2007) A Global Error Bound via the SQP Me-thod for Constrained Optimization Problem. Journal of Industrial and Management Optimization, 3, 775-781. http://www.aimsciences.org/journals/displayArticles.jsp?paperID=2746

*通讯作者。

期刊菜单