Báo cáo hóa học: " Research Article A Note on Convergence Analysis of an SQP-Type Method for Nonlinear Semidefinite Programming" ppt

10 243 0
Báo cáo hóa học: " Research Article A Note on Convergence Analysis of an SQP-Type Method for Nonlinear Semidefinite Programming" ppt

Đang tải... (xem toàn văn)

Thông tin tài liệu

Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 2008, Article ID 218345, 10 pages doi:10.1155/2008/218345 Research Article A Note on Convergence Analysis of an SQP-Type Method for Nonlinear Semidefinite Programming Yun Wang, 1 Shaowu Zhang, 2 and Liwei Zhang 1 1 Department of Applied Mathematics, Dalian University of Technology, Dalian 116024, China 2 Department of Computer Science, Dalian University of Technology, Dalian 116024, China Correspondence should be addressed to Yun Wang, wangyun 3412@163.com Received 29 August 2007; Accepted 23 November 2007 Recommended by Kok Lay Teo We reinvestigate the convergence properties of the SQP-type method for solving nonlinear semidef- inite programming problems studied by Correa and Ramirez 2004.Weprove,underthestrong second-order sufficient condition with the sigma term, that the local SQP-type method is quadrati- cally convergent and the line search SQP-type method is globally convergent. Copyright q 2008 Yun Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 1. Introduction We consider the following nonlinear semidefinite programming: SDP min fx s.t.hx0,gx ∈S p  , 1.1 where x ∈ R n , f : R n →R, h : R n →R l ,andg : R n →S p are twice continuously differentiable functions, S p is the linear space of all p × p real symmetric matrices, and S p  is the cone of all p × p symmetric positive semidefinite matrices. Fares et al. 20021 studied robust control problems via sequential semidefinite pro- gramming technique. They obtained the local quadratic convergence rate of the proposed SQP- type method and employed a partial augmented Lagrangian method to deal with the problems addressed there. Correa and Ramirez 20042 systematically studied an SQP-type method for solving nonlinear SDP problems and analyzed the convergence properties, they obtained the global convergence and local quadratic convergence rate. Both papers used the same sub- problems to generate search directions, but employed different merit functions for line search. The convergence analysis of both papers depends on a set of second-order conditions without sigma term, which is stronger than no gap second-order optimality condition with sigma term. 2 Journal of Inequalities and Applications Comparing with the work by Correa and Ramirez 20042, in this note, we make some modifications to the convergence analysis, and prove that all results in 2 still hold under the strong second-order sufficient condition with the sigma term. It should be pointed out that the importance of exploring numerical methods for solving nonlinear semidefinite programming problems has been recognized in the optimization com- munity. For instance, Ko ˇ cvara and Stingl 3, 4 have developed PENNLP and PENBMI codes for nonlinear semidefinite programming and semidefinite programming with bilinear matrix inequality constraints, respectively. Recently, Sun et al. 20075 considered the rate of con- vergence of the classical augmented Lagrangian method and Noll 20076 investigated the convergence properties of a class of nonlinear Lagrangian methods. In Section 2, we introduce preliminaries including differential properties of the metric projector onto S p  and optimality conditions for problem 1.1.InSection 3, we prove, under the strong second-order sufficient condition with the sigma term, that the local SQP-type method has the quadratic convergence rate and the global algorithm with line search is convergence. 2. Preliminaries We use R m×n to denote the set of all the matrices of m rows and n columns. For A and B in R m×n , we use the Frobenius inner product A, B  trA T B , and the Frobenius norm A F   trA T A, where “tr” denotes the trace operation of a square matrix. For a given matrix A ∈S p , its spectral decomposition is A  PΛP T  P ⎛ ⎜ ⎝ λ 1 00 0 . . . 0 00λ p ⎞ ⎟ ⎠ P T , 2.1 where Λ is the diagonal matrix of eigenvalues of A and P is a corresponding orthogonal matrix. We can express Λ and P as Λ ⎛ ⎝ Λ α 00 000 00Λ γ ⎞ ⎠ ,P  P α P β P γ  , 2.2 where α, β, γ are the index sets of positive, zero, negative eigenvalues of A, respectively. 2.1. Semismoothness of the metric projector In this subsection, let X, Y ,andZ be three arbitrary finite-dimensional real spaces with a scalar product ·, · and its norm ·. We introduce some properties of the metric projector, especially its strong semismoothness. The next lemma is about the generalized Jacobian for composite functions, proposed in 7. Lemma 2.1. Let Ψ : X→Y be a continuously differentiable function on an open neighborhood  N of x and let Ξ : O⊆Y→Z be a locally Lipschitz continuous function on the open set O containing Yun Wang et al. 3 y :Ψx. Suppose that Ξ is directionally differentiable at every point in O and that J x Ψx : X→Y is onto. Then it holds that ∂ B Φx∂ B ΞyJ x Ψx, 2.3 where Φ :  N→Z is defined by Φx :ΞΨx, x ∈  N. The following concept of semismoothness was first introduced by Mifflin 8 for func- tionals and was extended by Qi and Sun in 9 to vector valued functions. Definition 2.2. Let Φ : O⊆X→Y be a locally Lipschitz continuous function on the open set O. One says that Φ is semismooth at a point x ∈Oif iΦis directionally differentiable at x; ii for any Δx ∈ X and V ∈ ∂Φx Δx with Δx→0, Φx Δx − Φx − V Δxo  Δx  . 2.4 Furthermore, Φ is said to be strongly semismooth at x ∈Oif Φ is semismooth at x and for any Δx ∈ X and V ∈ ∂Φx Δx with Δx→0, Φx Δx − Φx − V ΔxO  Δx 2  . 2.5 Let D be a closed convex set in a Banach space Z,andletΠ D : Z→Z be the metric projector over D.Itiswellknownin10, 11 that Π D · is F-differentiable almost everywhere in Z and for any y ∈ Z, ∂Π D y is well defined. Suppose A ∈S p , then it has the spectral decomposition as 2.1, then the merit projector of A to S p  is denoted by Π S p  A and Π S p  AP ⎛ ⎜ ⎝  λ 1   00 0 . . . 0 00  λ p   ⎞ ⎟ ⎠ P T , 2.6 where λ i    max {0,λ i },i 1, ,p. For our discussion, we know from 12 that the projection operator Π S p  · is directionally differentiable everywhere in S p  and is a strongly semismooth matrix-valued function. In fact, for any A ∈S p , H ∈S p  , there exists V ∈ ∂Π S p  A  H, satisfying Π S p  A  HΠ S p  AV HO  H 2  . 2.7 2.2. Optimality conditions Let the Lagrangian function of 1.1 be Lx, λ, μfx  λ, hx    μ, gx  . 2.8 4 Journal of Inequalities and Applications Robinson’s constraint qualificationCQ is said to hold at a feasible point x if 0 ∈ int  hx gx    Jhx Jgx  R n −  {0} S p   . 2.9 If x is a locally optimal solution to 1.1 and Robinson’s CQ holds at x, then there exist Lagrangian multipliers λ, μ ∈ R l ×S p , such that the following KKT condition holds: 0  ∇ x Lx, λ, μ∇fxJhx ∗ λ  Jgx ∗ μ, 0  hx, g xΠ S p   g xμ  , 2.10 which is equivalent to F x, λ, μ0, where Fx, λ, μ : ⎛ ⎜ ⎝ ∇fxJhx ∗ λ  Jgx ∗ μ hx gx − Π S p   gxμ  ⎞ ⎟ ⎠ . 2.11 Let Λ x be the set of all the Lagrangian multipliers satisfying 2.10.ThenΛx is a nonempty, compact convex set of R l ×S p if and only if Robinson’s CQ holds at x, see 13. Moreover, it follows from 13 that the constraint nondegeneracy condition is a sufficient con- dition for Robinson constraint qualification. In the setting of the problem 1.1, the constraint nondegeneracy condition holding at a feasible point x can be expressed as  Jh x Jg x  R n   {0} lin  T S p   g x     R l S p  , 2.12 where linT S p  gx is the lineality space of the tangent cone of S p  at gx.Ifx, a locally optimal solution to 1.1, is nondegenerate, then Λ x is a singleton. For a KKT point  x, λ, μ of 1.1, without loss of generality, we assume that gx and μ have the spectral decomposition forms g xP ⎛ ⎝ Λ α 00 000 000 ⎞ ⎠ P T ,μ P ⎛ ⎝ 00 0 00 0 00Λ γ ⎞ ⎠ P T . 2.13 We state the strong second-order sufficient condition SSOSC coming from 7. Definition 2.3. Let x be a stationary point of 1.1 such that 2.12 holds at x. One says that the strong second-order sufficient condition holds at x if  d, ∇ 2 xx Lx, λ, μd  − Υ g x   μ, Jgxd  > 0, ∀d ∈ aff Cx \{0}, 2.14 where { λ, μ} Λx ⊂ R l ×S p ,aff Cx is the affine hull of the critical cone Cx: aff C x  d : Jhxd  0,P T β  Jg xd  P γ  0,P T γ  Jg xd  P γ  0  . 2.15 And the linear-quadratic function Υ B : S p ×S p →R is defined by Υ B D, A : 2  D, AB † A  , D, A ∈S p ×S p , 2.16 B † is the Moore-Penrose pseudoinverse of B. Yun Wang et al. 5 The next proposition relates the SSOSC and nondegeneracy condition to nonsingularity of Clarke’s Jacobian of the mapping F defined by 2.11. The details of this proof can be found in 7. Proposition 2.4. Let  x, λ, μ be a KKT point of 1.1. If nondegeneracy condition 2.12 and SSOSC 2.14 hold at x, then any element in ∂Fx, λ, μ is nonsingular, where F is defined by 2.11. 3. Convergence analysis of the SQP-type method In this section, we analyze the local quadratic convergence rate of an SQP-type method and then prove that the SQP-type method proposed in 2 is globally convergent. The analysis is based on the strong second-order sufficient condition, which is weaker than the conditions used in 1, 2. 3.1. Local convergence rate Linearizing 1.1 at the current point x k ,λ k ,μ k , we obtain the following tangent quadratic problem: min Δx ∇f  x k  T Δx  1 2 Δx T ∇ 2 xx L  x k ,λ k ,μ k  Δx, s.t.h  x k   Jh  x k  Δx  0,g  x k   Jg  x k  Δx ∈S p  , 3.1 where ∇ 2 xx Lx k ,λ k ,μ k J x ∇ x Lx k ,λ k ,μ k .LetΔx k ,λ k QP ,μ k QP  be a KKT point of 3.1, then we have  FΔx k ,λ k QP ,μ k QP ; x k ,λ k ,μ k 0, where  F  ζ, η, ξ; x k ,λ k ,μ k  : ⎛ ⎜ ⎝ ∇f  x k   ∇ 2 xx L  x k ,λ k ,μ k  ζ  Jh  x k  ∗ η  Jg  x k  ∗ ξ hx k Jhx k ζ g  x k   Jg  x k  ζ − Π S P   g  x k   Jg  x k  ζ  ξ  ⎞ ⎟ ⎠ . 3.2 The following algorithm is an SQP-type algorithm for solving 1.1, which is based on computing at each iteration a primal-dual stationary point Δx k ,λ QP k ,μ QP k  of 3.1. Algorithm 3.1 Step 1. Given an initial iterate point x 1 ,λ 1 ,μ 1 . Compute hx 1 , gx 1 , ∇fx 1 , Jhx 1  and Jgx 1 . Set k : 1. Step 2. If ∇ x Lx k ,λ k ,μ k 0, hx k 0, gx k  ∈S P  , stop. Step 3. Compute ∇ 2 xx Lx k ,λ k ,μ k , and find a solution Δx k ,λ k QP ,μ k QP  to 3.1. Step 4. Set x k1 : x k Δx k ,λ k1 : λ k QP ,μ k1 : μ k QP . Step 5. Compute hx k1 , gx k1 , ∇fx k1 , Jhx k1  and Jgx k1 . Set k : k  1 and go to step 2. From item f of 7, Theorem 4.1, we obtain the error between Δx k ,λ QP k ,μ QP k  and  x, λ, μ directly. Theorem 3.2. Suppose that f, h, g are twice continuously differentiable and their derivatives are lo- cally Lipschitz in a neighborhood of a local solution x to 1.1. Suppose nondegeneracy condition 2.12 6 Journal of Inequalities and Applications and SSOSC 2.14 hold at x. Then there exists a neighborhood U of x, λ, μ such that if x k ,λ k ,μ k  in U, 3.1 has a local solution Δx k together with corresponding Lagrangian multiplies λ k QP ,μ k QP  satisfying   Δx k      λ k QP − λ      μ k QP − μ    O     x k ,λ k ,μ k  −  x, λ, μ     . 3.3 Now we are in a position to state that the sequence of primal-dual points generated by Algorithm 3.1 has quadratic convergence rate. Theorem 3.3. Suppose that f, h, g are twice continuously differentiable and their derivatives are lo- cally Lipschitz in a neighborhood of a local solution x to 1.1. Suppose nondegeneracy condition 2.12 and SSOSC 2.14 hold at x. Consider Algorithm 3.1,inwhichΔx k is a minimum norm station- ary point of the tangential quadratic problem 3.1. Then there exists a neighborhood U of  x, λ, μ such that, if x 1 ,λ 1 ,μ 1  ∈ U, Algorithm 3.1 is well defined and the sequence {x k ,λ k ,μ k } converges quadratically to  x, λ, μ. Proof. By Theorem 3.2,weknowAlgorithm 3.1 is well defined. Let δ k :   x k ,λ k ,μ k  − x, λ, μ   , 3.4 then Δx k  O  δ k  ,λ k1 − λ  O  δ k  ,μ k1 − μ  O  δ k  , 3.5 where Δx k is the minimum norm solution to 3.1,andλ k1  λ k QP , μ k1  μ k QP are the associated multipliers. Using Taylor expansion of 3.2 at  x, λ, μ,notingthat∇ x Lx, λ, μ0, x k1  x k Δx k ,and3.5,weobtain ∇ 2 xx Lx, λ, μ  x k1 − x   Jhx ∗  λ k1 − λ   Jgx ∗  μ k1 − μ   O  δ 2 k  , Jh x  x k1 − x   O  δ 2 k  . 3.6 As the projection operator Π S p  · is strongly semismooth, we have that there exists V ∈ ∂Π S p  gxμ such that Π S p   g xμ  Π S p   gx k Jgx k Δx k  μ k QP   V  g xμ − gx k  −Jgx k Δx k − μ k QP   O    g xμ − gx k  −Jgx k Δx k − μ k QP   2  . 3.7 Since g xμ − g  x k  −Jg  x k  Δx k − μ k QP  Jg  x k  x − x k1    μ − μ k QP   O  δ 2 k  , 3.8 we have Π S p  gx k Jgx k Δx k  μ k QP  Π S p  gxμ − V Jgx k x − x k1 μ − μ k QP   Oδ 2 k . 3.9 Yun Wang et al. 7 Noting the fact that g xΠ S p  gxμ, by Taylor expansion of the third equation of 3.2 at  x, μ,weobtain V − IJg xx k1 − xV μ k1 − μOδ 2 k . 3.10 Therefore, we can conclude that ⎛ ⎜ ⎜ ⎝ ∇ 2 xx Lx, λ, μ Jhx ∗ Jgx ∗ Jhx 00 −Jg xV Jgx 0 V ⎞ ⎟ ⎟ ⎠ ⎛ ⎜ ⎜ ⎝ x k1 − x λ k1 − λ μ k1 − μ ⎞ ⎟ ⎟ ⎠  O  δ 2 k  . 3.11 Since the nondegeneracy condition 2.12 and SSOSC 2.14 hold, we have from Propo- sition 2.4 that 3.11 implies the quadratic convergence of the sequence {x k ,λ k , μ k }. 3.2. The global convergence The tangential quadratic problem constrained here is slightly more general than 3.1 in the sense that the Hessian of the Lagrangian ∇ 2 xx Lx k ,λ k ,μ k  is replaced by some positive definite matrix M k . Thus the tangential quadratic problem in Δx now becomes min Δx ∇f  x k  T Δx  1 2 Δx T M k Δx, s.t.h  x k   Jh  x k  Δx  0, g  x k   Jg  x k  Δx ∈S p  . 3.12 The KKT systemof 3.12 is ∇f  x k   M k Δx k  Jh  x k  ∗ λ k QP  Jg  x k  ∗ μ k QP  0,h  x k   Jh  x k  Δx k  0, g  x k   Jg  x k  Δx k − Π S p   g  x k   Jg  x k  Δx k  μ k QP   0. 3.13 To obtain theglobal convergence, we use the Han penalty function given by 14, as a merit function and Armijo line search. For problem 1.1, the Han penalty function is defined by Θ σ xfxσ   hx   − σλ min  gx  − , 3.14 where λ min gx is the smallest eigenvalue of gx, · − denote min {·, 0} and σ>0 is a positive constant. The following proposition comes from 2 directly. Proposition 3.4. i If f, h, g have a directional derivative at x in the direction d ∈ R n ,thenΘ σ has also a directional derivative at x in the direction d. If, in addition, x is feasible for 1.1,wehave Θ  σ x; df  x; dσ   h  x; d   − σλ min  N T JgxdN  , 3.15 where N ν 1 , ,ν r  is the matrix whose columns ν i form an orthonormal basis of Kergx. ii If x is a feasible point of 1.1 and Θ σ has a local minimum at x,thenx is the local solution to 1.1. Furthermore, if f, h, g are differentiable at x and nondegeneracy condition 2.12 holds at x, then σ≥ max { λ, tr−μ}. iii If μ<0 and σ≥ max {λ, tr−μ},thenL·,λ,μ ≤ Θ σ ·. 8 Journal of Inequalities and Applications To discuss the conditions ensuring the exactness of Θ σ , we need the following lemma from 3.10. Lemma 3.5. Suppose nondegeneracy condition 2.12 and SSOSC 2.14 hold at x. Then there exists c 0 > 0, such that for any c>c 0 there exist a neighborhood V of x and a neighborhood U of λ, μ,for any λ, μ ∈ U, the problem min L c x, λ, μ s.t.x∈ V 3.16 has a unique solution denote x c λ, μ. The function x c ·, · is locally Lipschitz continuous and semis- mooth on U. Furthermore, there exists ρ>0, for any λ, μ ∈ U,   x − x c λ, μ   ≤ ρ   λ, μ −  λ, μ   /c, 3.17 where L c x, λ, μ : fx  hx,λ   c 2   hx   2  1 2c    Π S p   − μ − cgx    2 F −μ 2 F  3.18 is the augmented Lagrangian function with the penalty parameter c for 1.1. Theorem 3.6. Suppose that f, h, g are twice differentiable around a local solution x to 1.1,atwhich nondegeneracy condition 2.12 and SSOSC 2.14 hold. If σ>max { λ, tr−μ},thenΘ σ has a strict local minimum at x. Proof. For the definition of the projection operator Π S p  ·,wehave Π S p   − μ − cgx   −μ − cgxΠ S p   cgxμ  , 3.19 and for any W ∈S p  , c>0,   Π S p   cgxμ  −  cgxμ    2 F ≤   W −  cgxμ    2 F . 3.20 Then   Π S p   cgxμ  − cgx   2 F − 2  μ, Π S p   cgxμ  − cgx  ≤−2  μ, W − cgx     W − cgx   2 F 3.21 holds for any W ∈S p  .Sotakingμ  μ and W  cΠ S p  gx,weobtainthat   Π S p   cgx μ − cgx   2 F − 2  μ, Π S p   cgx μ − cgx  ≤−2c  μ, Π S p   − gx   c 2   Π S p   − gx    2 F , 3.22 which implies L c x, λ, μ ≤ fx  λ, hx   c 2   hx   2 −  μ, Π S p   − gx   c 2   Π S p   − gx    2 F ≤ fx   hx     λ  c 2   hx     λ max  Π S p   − gx   tr− μ c 2 p  i1 λ i  Π S p   − gx   . 3.23 Yun Wang et al. 9 Since σ>max { λ, tr−μ}, for any fixed c>0, there exists a neighborhood V c of x such that L c x, λ, μ ≤ fxσ   hx    σλ max  Π S p   − gx  Θ σ x, ∀x ∈ V c . 3.24 From Lemma 3.5, we know that there exist an r>c 0 and a neighborhood V r of x where x is a strict minimum of L r ·, λ, μ. So we can conclude that x is a strict minimum of Θ σ on V c ∩ V r . Let us outline the line-search SQP-type algorithm that uses the merit function Θ σ · de- fined in 3.14 and the parameter updating scheme from 14, which is a generalized version to the algorithm in 2. Algorithm 3.7 Step 1. Given a positive number σ>0,  ∈ 0, 1/2, β ∈ 0, 1/2. Choose an initial iterate x 1 ,λ 1 ,μ 1  ∈ R n × R l ×S p . Compute fx 1 , hx 1 , gx 1 , ∇fx 1 , Jhx 1  andJgx 1 . Set k : 1,σ 1  σ. Step 2. If ∇ x Lx k ,λ k ,μ k 0,hx k 0,gx k  ∈S p  , stop. Step 3. Compute a symmetric matrix M k and find a solution Δx k ,λ k QP ,μ k QP  to 3.12. Step 4. Adapt σ k . if σ k−1 ≥ max {tr−μ k1 , λ k1 }  σ then σ k  σ k−1 else σ k  max {1.5σ k−1 , max {tr−μ k1 , λ k1 }  σ} Step 5. Compute w k : −  Δx k ,M k Δx k    μ k QP ,g  x k    λ k QP ,h  x k  − σ k   h  x k     σ k λ min  g  x k  − . 3.25 Using backtracking line search rule to compute the step length α k : Step 6. set i  0, α k,0  1; Step 7. if Θ σ k  x k  αΔx k  ≤ Θ σ k  x k   αw k 3.26 holds for α  α k,i ,thenα k  α and stop the line search. Step 8. else, choose α k,i1 ∈ βα k,i , 1 − ββα k,i ; Step 9. set i : i  1, go to step 7 Step 10. Set x k1 : x k  α k Δx k , λ k1 : λ k QP , μ k1 : μ k QP . Step 11. Compute fx k1 ,hx k1 ,gx k1 , ∇fx k1 , Jhx k1  and Jgx k1 . Set k : k  1 and go to step 2. Now we are in a position to state the global convergence of the line search SQP Algorithm 3.7, whose proof can be found in 2. 10 Journal of Inequalities and Applications Theorem 3.8. Suppose that f, h, g are continuously differentiable and their derivatives are Lipschitz continuous. Consider Algorithm 3.7, if positive definite matrices M k and M −1 k are bounded, then one of the following situations occurs: i the sequence {σ k } is unbounded, in which case {λ k1 ,μ k1 } is also unbounded; ii there exists an index k 2 such that σ k  σ for any k≥k 2 , and one of the following situations occurs: aΘ σ x k →∞, b ∇ x Lx k ,λ k ,μ k →0,hx k →0,λ min gx k  − →0, and μ k1 ,gx k →0. Acknowledgments The research is supported by the National Natural Science Foundation of China under Project no. 10771026 and by the Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministry of China. References 1 B. Fares, D. Noll, and P. Apkarian, “Robust control via sequential semidefinite programming,” SIAM Journal on Control and Optimization, vol. 40, no. 6, pp. 1791–1820, 2002. 2 R. Correa and H. C. Ramirez, “A global algorithm for nonlinear semidefinite programming,” SIAM Journal on Optimization, vol. 15, no. 1, pp. 303–318, 2004. 3 M. Ko ˇ cvara and M. Stingl, “Pennon: a code for convex nonlinear and semidefinite programming,” Optimization Methods & Software, vol. 18, no. 3, pp. 317–333, 2003. 4 M. Ko ˇ cvara and M. Stingl, “Solving nonconvex SDP problems of structural optimization with stability control,” Optimization Methods & Software, vol. 19, no. 5, pp. 595–609, 2004. 5 D. Sun, J. Sun, and L. Zhang, “The rate of convergence of the augmented Lagrangian method for nonlinear semidefinite programming,” Mathematical Programming, 2008. 6 D. Noll, “Local convergence of an augmented Lagrangian method for matrix inequality constrained programming,” Optimization Methods and Software, vol. 22, no. 5, pp. 777–802, 2007. 7 D. Sun, “The strong second-order sufficient condition and constraint nondegeneracy in nonlinear semidefinite programming and their implications,” Mathematics of Operations Research, vol. 31, no. 4, pp. 761–776, 2006. 8 R. Mifflin, “Semismooth and semiconvex functions in constrained optimization,” SIAM Journal on Control and Optimization, vol. 15, no. 6, pp. 959–972, 1977. 9 L. Q. Qi and J. Sun, “A nonsmooth version of Newton’s method,” Mathematical Programming, vol. 58, no. 1–3, pp. 353–367, 1993. 10 E. H. Zarantonello, “Projections on convex sets in Hilbert space and spectral theory—I: projections on convex sets,” in Contributions to Nonlinear Functional Analysis (Proceedings of a Symposium held at the Mathematics Research Center, University of Wisconsin, Madison, Wis., 1971), E. H. Zarantonello, Ed., pp. 237–341, Academic Press, New York, NY, USA, 1971. 11 E. H. Zarantonello, “Projections on convex sets in Hilbert space and spectral theory—II: spectral the- ory,” in Contributions to Nonlinear Functional Analysis (Proceedings of a Symposium held at the Mathematics Research Center, University of Wisconsin, Madison, Wis., 1971), E. H. Zarantonello, Ed., pp. 343–424, Aca- demic Press, New York, NY, USA, 1971. 12 D. Sun and J. Sun, “Semismooth matrix-valued functions,” Mathematics of Operations Research,vol.27, no. 1, pp. 150–169, 2002. 13 J. F. Bonnans and A. S hapiro, Perturbation Analysis of Optimization Problems , Springer Series in Opera- tions Research, Springer, New York, NY, USA, 2000. 14 J. F. Bonnans, J. C. Gilbert, C. Lemar ´ echal, and C. A. Sagastiz ´ abal, Numerical Optimization. Theoretical and Practical Aspects, Universitext, Springer, Berlin, Germany, 2003. . augmented Lagrangian method for nonlinear semidefinite programming,” Mathematical Programming, 2008. 6 D. Noll, “Local convergence of an augmented Lagrangian method for matrix inequality constrained programming,”. Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 2008, Article ID 218345, 10 pages doi:10.1155/2008/218345 Research Article A Note on Convergence Analysis of an SQP-Type Method. con- vergence of the classical augmented Lagrangian method and Noll 20076 investigated the convergence properties of a class of nonlinear Lagrangian methods. In Section 2, we introduce preliminaries

Ngày đăng: 22/06/2014, 18:20

Tài liệu cùng người dùng

  • Đang cập nhật ...

Tài liệu liên quan