<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AM</journal-id><journal-title-group><journal-title>Applied Mathematics</journal-title></journal-title-group><issn pub-type="epub">2152-7385</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/am.2012.330198</article-id><article-id pub-id-type="publisher-id">AM-24125</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  A New Lagrangian Multiplier Method on Constrained Optimization
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>ou-Lin</surname><given-names>Shang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Sheng-Li</surname><given-names>Guo</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Xiang-Yi</surname><given-names>Jiang</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mathematics, Henan University of Science and Technology, Luoyang, China</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>mathshang@sina.com(OS)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>31</day><month>10</month><year>2012</year></pub-date><volume>03</volume><issue>10</issue><fpage>1409</fpage><lpage>1414</lpage><history><date date-type="received"><day>July</day>	<month>8,</month>	<year>2012</year></date><date date-type="rev-recd"><day>August</day>	<month>8,</month>	<year>2012</year>	</date><date date-type="accepted"><day>August</day>	<month>15,</month>	<year>2012</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  In this paper, a new augmented Lagrangian function with 4-piecewise linear NCP function is introduced for solving nonlinear programming problems with equality constrained and inequality constrained. It is proved that a solution of the original constrained problem and corresponding values of Lagrange multipliers can be found by solving an unconstrained minimization of the augmented Lagrange function. Meanwhile, a new Lagrangian multiplier method corresponding with new augmented Lagrangian function is proposed. And this method is implementable and convergent.
 
</p></abstract><kwd-group><kwd>Nonlinear Programming; NCP Function; Lagrange Function; Multiplier; Convergence</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Considering the following nonlinear inequality constrained optimization Problem (NLP):</p><disp-formula id="scirp.24125-formula70078"><label>(1)</label><graphic position="anchor" xlink:href="22-7400953\2a2fdc2e-275c-44ba-bf95-39186dab2aef.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="22-7400953\d910101d-0efc-4304-9358-f73d28d4f7d4.jpg" /> and</p><p><img src="22-7400953\810341c1-5d0e-45c6-ab61-eb0438952c7e.jpg" /></p><p>are continuously differentiable functions.</p><p>We denote by</p><p><img src="22-7400953\5660ffc1-7b5d-487c-845e-13f2425478a4.jpg" /></p><p>the feasible set of the problem (NLP).</p><p>The Lagrangian function associated with the problem (NLP) is the function</p><p><img src="22-7400953\0804cd83-538c-4dcf-a70d-c3796a776d45.jpg" /></p><p>where</p><p><img src="22-7400953\5a4067b8-3c18-413f-a72f-b4188230240c.jpg" /></p><p>are the multiplier vectors, For simplicity, we use <img src="22-7400953\8b431e67-8051-4c8c-a670-39aee401f09f.jpg" /> to denote the column vector <img src="22-7400953\6803708a-d03e-4e92-b5df-70f4302bbba7.jpg" /></p><p>Defintion 1.1. A point <img src="22-7400953\8aa7e62a-4480-4336-9c74-63b465be2fba.jpg" /> is called a Karush-Kuhn-Tucker (KKT) point or a KKT pair of Problem (NLP), if it satisfies the following conditions:</p><disp-formula id="scirp.24125-formula70079"><label>(2)</label><graphic position="anchor" xlink:href="22-7400953\31382f63-9764-4d88-862c-dc0fdf175514.jpg"  xlink:type="simple"/></disp-formula><p>where<img src="22-7400953\d9b81d85-a311-425d-b3c5-1b84e69dcf0b.jpg" />, we also say <img src="22-7400953\d1682062-5945-4718-941f-3c535d94f7fb.jpg" /> is a KKT point if there exists a <img src="22-7400953\f603c13b-5c1f-47fe-9b9d-c79a152a4074.jpg" /> such that <img src="22-7400953\7c79f09c-1b19-4d9a-ad53-169d3e5de1a1.jpg" /> satisfies (2).</p><p>For the nonlinear inequality constrained optimization problem (NLP), there are many practical methods to solve it, such as augmented Lagrangian function method [1-6], Trust-region filter method [7,8], QP-free feasible method [9,10], Newton iterative method [11,12], etc. As we know, Lagrange multiplier method is one of the efficient methods to solve problem (NLP). Pillo and Grippo in [1-3] proposed a class of augmented Lagrange function methods which have nice equivalence between the unconstrained optimization and the primal constrained problem and get good convergence properties of the related algorithm. However, a max function is used for these methods which may be not differentiable at infinite numbers of points. To overcome this shortcoming, Pu in [<xref ref-type="bibr" rid="scirp.24125-ref4">4</xref>] proposed a augmented Lagrange function with FischerBurmeister nonlinear NCP function and Lagrange multiplier methods. Pu and Ding in [<xref ref-type="bibr" rid="scirp.24125-ref6">6</xref>] proposed a Lagrange multiplier methods with 3-piecewise linear NCP function. In this paper, a new class augmented Lagrange function with 4-piecewise linear NCP function and some Lagrange multiplier methods are proposed for the minimization of a smooth function subject to smooth inequality constraints and equality constrains.</p><p>The paper is organized as follows: In the next section we give some definitions and properties about NCP function, and then define a new augmented Lagrange function with 4-piecewise NCP function. In Section three, we give the algorithm. In Section four, we prove convergence of the algorithm. Some conclusions are given in Section five.</p></sec><sec id="s2"><title>2. Preliminaries</title><p>In this section, we recall some definitions and define a new Lagrange multiplier function with 4-piecewise NCP function.</p><p>Definition 2.1 (NCP pair and SNCP pair). We call a pair (a, b) to be an NCP pair if <img src="22-7400953\1bea4d3e-1e56-44ce-8c94-b5e25f71c9d7.jpg" /> and ab = 0; and call (a, b) to be an SNCP pair if (a, b) is a pair and<img src="22-7400953\3a602d72-16c1-4141-830e-c31a8cb03489.jpg" />.</p><p>Definition 2.2 (NCP function). A function <img src="22-7400953\668071b3-89b2-439c-93bd-925875850b9d.jpg" /> is called an NCP function if <img src="22-7400953\f0ead643-3266-4d39-b86c-c7d1e65b9b2a.jpg" /> if and only <img src="22-7400953\9f800d50-23fc-420a-9f1c-cef0e8e13f16.jpg" /> is an NCP pair.</p><p>In this paper, we propose a new 4-piecewise linear NCP function <img src="22-7400953\2e19af7e-0e94-4d83-9b3b-114ae36e8aa4.jpg" /> is as follows:</p><disp-formula id="scirp.24125-formula70080"><label>(3)</label><graphic position="anchor" xlink:href="22-7400953\8753c9b3-9126-4e28-8674-68be625463de.jpg"  xlink:type="simple"/></disp-formula><p>If<img src="22-7400953\dd74f42d-a039-432c-9883-99ab9c4e38d8.jpg" />, then</p><disp-formula id="scirp.24125-formula70081"><label>(4)</label><graphic position="anchor" xlink:href="22-7400953\92d5239f-8d73-41da-921b-958583be80b6.jpg"  xlink:type="simple"/></disp-formula><p>and</p><disp-formula id="scirp.24125-formula70082"><label>(5)</label><graphic position="anchor" xlink:href="22-7400953\46031ef0-594d-4859-a409-755187c4ff84.jpg"  xlink:type="simple"/></disp-formula><p>It is easy to check the following propositions:</p><p>1)<img src="22-7400953\2e1ba6a5-9fbb-44d4-834d-826f1ee08b68.jpg" />;</p><p>2) The square of <img src="22-7400953\143a8ab3-3eb1-4cc6-898e-88b089bfbd83.jpg" /> is continuously differentiable;</p><p>3) <img src="22-7400953\e71dc412-bcc6-44f9-9aa5-0765eff90d65.jpg" />is twice continuously differentiable everywhere except at the origin but it is strongly semi-smooth at the origin.</p><p>Let</p><p><img src="22-7400953\04626070-1898-43c3-a1b1-c39d5f81afbb.jpg" /></p><p>where<img src="22-7400953\73bb698d-1021-45e5-8c3e-2e2a308f1219.jpg" /> is a parameter. <img src="22-7400953\9e1ddd23-214a-4c4f-a7da-29e9a01d8338.jpg" />if and only if<img src="22-7400953\773ffaf8-ded3-44d5-bff0-9517cfb2465a.jpg" />, <img src="22-7400953\5447b0c3-308e-4516-9981-62a468713769.jpg" />and <img src="22-7400953\774fae47-57c9-4e44-bd55-b0d101c7fd34.jpg" /> for any<img src="22-7400953\2db08c8b-a960-4d1f-a87b-64b9c90415c1.jpg" />.</p><p>We construct function:</p><p><img src="22-7400953\1d3150df-2144-4e6e-80a1-9e3f7eb836ab.jpg" /></p><p>Clearly, the KKT point condition (2) is equivalently reformulate as the condition:</p><p><img src="22-7400953\81fa4b1a-0b5f-437a-a2a3-bfe2530523c3.jpg" /></p><p>If<img src="22-7400953\80ddbd8f-b44f-4e35-a8fc-e23a187bb2e3.jpg" />, then <img src="22-7400953\88373512-e2c4-48b9-87f3-baee8d5af213.jpg" /> is continuously differentiable at<img src="22-7400953\de059fbf-a891-4f07-985b-4e5588400fbf.jpg" />. We have</p><disp-formula id="scirp.24125-formula70083"><label>(6)</label><graphic position="anchor" xlink:href="22-7400953\37bb35a3-a8fb-4913-999b-a08675e5e584.jpg"  xlink:type="simple"/></disp-formula><p>where <img src="22-7400953\b846cc71-1b9d-4a00-ab13-5095efea9a6b.jpg" /> is the ith column of the unit matrix, its jth element is 1, and other elements are 0, in this paper take k = 1.</p><p>If<img src="22-7400953\0592b0ac-1edb-4a33-a34c-22178a5297e6.jpg" />, and then <img src="22-7400953\42acf896-ef83-4823-ad6f-c73a86dbcdce.jpg" /> is strongly semi-smooth and direction differentiable at<img src="22-7400953\3a4158c9-9ed1-4aab-9f74-2e3e674bd7be.jpg" />. We have</p><disp-formula id="scirp.24125-formula70084"><label>(7)</label><graphic position="anchor" xlink:href="22-7400953\85e05ef0-97d1-41bd-958b-110980ea080b.jpg"  xlink:type="simple"/></disp-formula><p>For Problem (NLP), we define a Di Pillo and Grippo type Lagrange multiplier function with 4-piecewise linear NCP function is as following:</p><disp-formula id="scirp.24125-formula70085"><label>(8)</label><graphic position="anchor" xlink:href="22-7400953\9824bcaf-e733-4808-a701-558adc31e0e7.jpg"  xlink:type="simple"/></disp-formula><p>where</p><p><img src="22-7400953\5cb9d408-812e-409d-ac8e-6b25264edd62.jpg" />are the Lagrange multiplier, C and D are positive parameters.</p><p>In this section, we gave some assumptions as follows:</p><p>Assumpion 1 f, <img src="22-7400953\9b22ef14-f33e-47e6-91b1-a1f46bbc696e.jpg" />, <img src="22-7400953\6e759caa-6f17-4ca6-8e66-95cb19447ad1.jpg" />, <img src="22-7400953\e14bbd16-cd65-4348-bae2-67dedab0a21a.jpg" /></p><p><img src="22-7400953\a3e7f851-d06b-4e22-9d6c-3fe85e40fc4d.jpg" />are twice Lipschitz continuously differentiable.</p><p>Define index set <img src="22-7400953\fe8fbdd8-4427-4b53-842f-1ce558d28d01.jpg" /> and <img src="22-7400953\d44d6ea3-32a0-432e-af7a-2ded37705959.jpg" /> as follows:</p><p><img src="22-7400953\356d1f31-b3f7-4909-80d1-fb94572df314.jpg" /></p><p>for any<img src="22-7400953\bd13dedd-c11e-4f75-af70-097669a43627.jpg" />, according to definition of<img src="22-7400953\a42470d6-4103-47f5-a5f3-994257d3ab6d.jpg" />, have</p><p><img src="22-7400953\6e556053-fd3d-43a1-bc87-f36fbbace0e6.jpg" /></p><p>for any<img src="22-7400953\3d34c09a-189b-4cc5-b9c3-8fc3a8242e16.jpg" />, we have 1) if<img src="22-7400953\a9a3b8f5-5b27-4762-8bb1-a8a20a7252a3.jpg" />, we have</p><disp-formula id="scirp.24125-formula70086"><label>(9)</label><graphic position="anchor" xlink:href="22-7400953\517a4a3d-3192-4112-8d6b-75c67b71d038.jpg"  xlink:type="simple"/></disp-formula><p>The gradient of <img src="22-7400953\53b66b7f-f230-4da5-b296-3009ec952081.jpg" /> is</p><disp-formula id="scirp.24125-formula70087"><label>(10)</label><graphic position="anchor" xlink:href="22-7400953\6943e8ab-1720-41e6-977f-800746df0cf5.jpg"  xlink:type="simple"/></disp-formula><p>The Henssian matrix of <img src="22-7400953\283ffac9-d216-4254-9204-ad40af5edaea.jpg" /> at KKT point <img src="22-7400953\6407eb86-b3f0-4569-9158-25d8fdd1287d.jpg" /> is</p><disp-formula id="scirp.24125-formula70088"><label>(11)</label><graphic position="anchor" xlink:href="22-7400953\72ba980c-cfcb-477d-af70-0ebe32cf9042.jpg"  xlink:type="simple"/></disp-formula><p>2) if<img src="22-7400953\fe93e0b0-3e94-400e-bdbb-72acacae3efd.jpg" />, then</p><disp-formula id="scirp.24125-formula70089"><label>(12)</label><graphic position="anchor" xlink:href="22-7400953\5113e67a-3a1a-4919-956a-167519342133.jpg"  xlink:type="simple"/></disp-formula><p>The gradient of <img src="22-7400953\d6a6a6eb-7e42-40c7-be32-8a432c0ef272.jpg" /> is</p><disp-formula id="scirp.24125-formula70090"><label>(13)</label><graphic position="anchor" xlink:href="22-7400953\4ea261a3-58bf-4e57-b10e-9ebb9149d74b.jpg"  xlink:type="simple"/></disp-formula><p>The Henssian matrix of <img src="22-7400953\391d3a54-25a1-4356-8d1f-e416efc0f220.jpg" /> at KKT point <img src="22-7400953\2e57132d-bec7-480e-a045-080413e67fbb.jpg" /> is</p><disp-formula id="scirp.24125-formula70091"><label>(14)</label><graphic position="anchor" xlink:href="22-7400953\fd4fc450-1267-46bc-8a21-3610ab78ac15.jpg"  xlink:type="simple"/></disp-formula><p>Definition 2.3 A point <img src="22-7400953\f34348f1-c410-4160-b9af-f8d01cb2e948.jpg" /> is said to satisfy the strong second-order sufficiency condition for problem (NLP) if it satisfies the first-order KKT condition and if <img src="22-7400953\f585da7b-9c03-42d6-ba7d-38f238edfde2.jpg" /> for all</p><p><img src="22-7400953\a97b69c8-ccd0-429f-91b0-c6c81eaf09e6.jpg" /></p><p>and<img src="22-7400953\68dc54dc-cb98-4b9a-b2df-ba792baea111.jpg" />.</p><p>Assumption 2 At any KKT point <img src="22-7400953\3e2dd9d0-a44f-4876-9b96-27e0292671b4.jpg" /> satisfied strong second-order sufficiency condition.</p><p>Lemma If <img src="22-7400953\88af61d0-eeb8-43aa-bed8-111b84e6c274.jpg" /> is a positive semi-definite matrix, for any<img src="22-7400953\a9a66921-23ac-4a99-b5d4-acfd3a9ed238.jpg" />, <img src="22-7400953\30cd7b78-8f09-437f-bc69-972eb140404d.jpg" />, matrix <img src="22-7400953\26af186c-79e0-4f8a-a0bd-ddc87a2af69e.jpg" /> satisfied<img src="22-7400953\530f3f40-c445-42ca-95db-cb80fb46e839.jpg" />, then exist<img src="22-7400953\3be2371b-2212-4e63-882e-480aa20331c8.jpg" />, for any<img src="22-7400953\c5b07fa0-3ec2-4d15-bd97-0fa7d087c823.jpg" />, <img src="22-7400953\685585d1-8cc4-4896-a3dc-e08dc89b0209.jpg" />is positive definite matrix (see [<xref ref-type="bibr" rid="scirp.24125-ref4">4</xref>]).</p><p>Theorem 2.1 If <img src="22-7400953\7fdfac68-a45b-4d54-b678-0704084e9eb7.jpg" /> is KKT point of problem (1), then for sufficiently large C and D, <img src="22-7400953\6eb6bc33-1241-4c8f-9451-ff6a35a39224.jpg" />is strong convex function at point<img src="22-7400953\53c8cf43-3bbd-4273-a266-12c3af4b896d.jpg" />.</p><p>Proof: Let <img src="22-7400953\9e824b0e-088a-4a8d-9fbc-113425985472.jpg" /></p><p><img src="22-7400953\21854b60-dd65-4cf6-a15f-e89cd251809b.jpg" /></p><p>for<img src="22-7400953\8abede21-5e66-42b4-832e-52bfd65547b6.jpg" />, we have</p><p><img src="22-7400953\337bd4e3-65bd-4e6c-beff-afcc3998da05.jpg" /></p><p>from A2, we have<img src="22-7400953\b49bb36d-962b-44f0-87f3-d3d17d17afee.jpg" />. Furthermore there is <img src="22-7400953\b5ca3367-ef42-4df5-a538-d5a691e30660.jpg" /> if<img src="22-7400953\3b2b066c-8adb-4eb2-b739-fafd3856388a.jpg" />, for any<img src="22-7400953\e5513397-5c50-42b9-a7a4-57d54fec7632.jpg" />, <img src="22-7400953\21988c8e-9379-4cbe-ae69-edbf06d6f291.jpg" />is positive definite matrix. And then for any <img src="22-7400953\d3af3053-d6fe-4325-b6b0-6514b2f2b2c8.jpg" /> and sufficiently large C and D have</p><p><img src="22-7400953\6ac18f8f-2408-417c-bd09-dc2b139a0c05.jpg" /></p><p>by its continuously, we may obtained that there is<img src="22-7400953\56a824cb-e84e-477a-a988-c74ca8779bb9.jpg" />, for all</p><p><img src="22-7400953\55befd38-65a9-42a1-adc6-d9701f37db09.jpg" /></p><p>we have <img src="22-7400953\7a6528b1-8d77-4a81-8de3-cc47ed78e0b1.jpg" /> the theorem hold.</p></sec><sec id="s3"><title>3. Lagrange Multiplier Algorithm</title><p>Step 0 Choose parameters<img src="22-7400953\75dc9882-a147-4432-9973-0f74d0db290c.jpg" />, <img src="22-7400953\c74e3122-e5d1-415d-af01-f51da3114e46.jpg" />, <img src="22-7400953\a0cd84b8-f020-445b-b1e7-82f20ca5e45f.jpg" />, <img src="22-7400953\950656ff-ea35-43ab-9069-10b3d88271e8.jpg" />, given point<img src="22-7400953\6ce71426-f4ab-4af1-8670-6a98b0d4a96b.jpg" />, and</p><p><img src="22-7400953\d6edc8cc-98e5-4bc2-88ff-311f2359959c.jpg" /></p><p>Let<img src="22-7400953\f676973e-7646-4a65-95f5-786839f66636.jpg" />.</p><p>Step 1 Solve following, we will obtain<img src="22-7400953\1d7e17e9-29db-4082-9ac8-77480b983cb8.jpg" />.</p><p><img src="22-7400953\63b51261-81ec-490e-a2a9-783478dd5a7e.jpg" /></p><p>if <img src="22-7400953\c78c145b-b5c2-41d3-9a53-dffb7d5e556c.jpg" /> and <img src="22-7400953\baa1f967-c5eb-456a-9737-8c64cda02d97.jpg" /> then stop.</p><p>Step 2 For<img src="22-7400953\c8e1bac9-2383-4deb-928c-42c8c230dc66.jpg" />, <img src="22-7400953\c48dcdcf-3680-47f3-b397-d8ab1748acc8.jpg" />, then <img src="22-7400953\1cf84e53-e909-4671-99e6-48f1b9845c4e.jpg" /> or<img src="22-7400953\be15512a-0ed5-4c88-83b0-ecb66d68aa16.jpg" />, for<img src="22-7400953\11545161-6a6b-4874-859a-73ed842673cf.jpg" />, if</p><p><img src="22-7400953\c8e5e723-bc60-49b3-924b-015923d10945.jpg" />then<img src="22-7400953\872a0e71-374b-442f-b2ab-c90baf49f5d8.jpg" />, or <img src="22-7400953\ff5861f3-e886-452d-af9d-7ea29e45f310.jpg" /></p><p>Step 3 Compute <img src="22-7400953\0cb401d7-2cff-4a7a-8f3f-f4e3d509ce1c.jpg" /> and <img src="22-7400953\a7220ec7-5170-4e01-96e0-5be1848694d8.jpg" /></p><p><img src="22-7400953\aae70818-ddc1-4910-8ec3-926f69e24a92.jpg" /></p><p>Step 4 Let k = k + 1, go to Step 1.</p></sec><sec id="s4"><title>4. Convergence of the Algorithm</title><p>In this section, we make a assumption follow as:</p><p>Assumption 3 For any<img src="22-7400953\98cc07bf-4be4-4041-8efe-23ead50b0463.jpg" />, <img src="22-7400953\bc7cdad0-71dd-445f-a333-da0a88ec0993.jpg" />, <img src="22-7400953\e1a99c1c-fc2a-4552-b9e0-cb8cc8c06864.jpg" />, <img src="22-7400953\13a9bef1-8eff-41f1-a783-8d81679bf702.jpg" />,</p><p><img src="22-7400953\81df30dd-a212-45bd-8297-45bb5640cbd3.jpg" />exists a minimizer point<img src="22-7400953\097c9c72-4f29-4146-a24e-868f69714878.jpg" />.</p><p>Theorem 4.1 Assume feasible set of problem (NLP) is non-empty set and <img src="22-7400953\f2efb087-3efe-45ce-a093-55adede6720d.jpg" /> is bounded, then algorithm is bound to stop after finite steps iteration.</p><p>Proof: Assume that the algorithm can not stop after finite steps iteration, by the sack of convenience, we define index set as following</p><p><img src="22-7400953\79312639-9692-4d1b-9f0d-3ee3bcd0dd23.jpg" /></p><p>according to assumption A3, it is clearly that <img src="22-7400953\c4fd6f2c-fc34-4c58-8073-54a40d7c87bc.jpg" /> or <img src="22-7400953\b2ef4237-44cd-4418-a3a5-779b76534864.jpg" /> are non-empty set. for any k, obtain</p><p><img src="22-7400953\27c31377-5c7a-48b6-a013-ffda79d5bbf4.jpg" /></p><p>from above assumption, we obtain that for any a <img src="22-7400953\1a2aa64f-b4b5-452a-a901-142207a3a965.jpg" /> there is<img src="22-7400953\09cc1ac6-d070-43ee-a9ff-d7989e17fdc0.jpg" />, for any <img src="22-7400953\694ccc47-32af-412b-adf7-dee88681efdf.jpg" /> and<img src="22-7400953\b94be4cd-a28c-4034-ac27-6958bf614ef5.jpg" />, <img src="22-7400953\d667049e-ea7b-4329-8228-9659d23c6b8d.jpg" />and<img src="22-7400953\a67386ea-0628-4783-8ca3-b4160c69b789.jpg" />, for sufficiently large k, it is not difficult to see that</p><p><img src="22-7400953\c18b5a44-2c74-4663-8266-c82b86510eb4.jpg" /></p><p>Or for any <img src="22-7400953\9e289384-4a97-4eea-9abf-fb82f5a8c957.jpg" /> and<img src="22-7400953\a7364add-d8a5-49fb-a3ea-85088074aa8f.jpg" />, <img src="22-7400953\fca72f6d-0a46-4b2b-adbc-31de2b746fa8.jpg" />and<img src="22-7400953\50c96d37-da72-4bcc-ab2c-580f4dbb5a2f.jpg" />, for sufficiently large k, have</p><p><img src="22-7400953\55597359-1cf8-4d9b-bf11-af502a4d5025.jpg" /></p><p>When <img src="22-7400953\0abeae01-83e5-480a-b6de-b73576927165.jpg" /> we can hold</p><p><img src="22-7400953\1e6f2622-e9c6-4da3-8668-e911bede43a3.jpg" /></p><p>Which contradicts A3, the theorem holds.</p><p>Theorem 4.2 Let <img src="22-7400953\02e2cbb3-375f-4c19-991f-b9245b853b4f.jpg" /> is a compact set, sequence <img src="22-7400953\e06d089c-995a-4d30-a8db-d92200207826.jpg" /> are generated by the algorithm, and<img src="22-7400953\6d9e2060-bbf7-4d64-b785-d8d1d8fc1f24.jpg" />, in algorithm, 0 take the place of<img src="22-7400953\18fe59b5-9125-4470-aa2f-fe6d9c657c34.jpg" />, either algorithm stops at its<img src="22-7400953\834995c1-568d-4b49-8b3a-fc3be2db9450.jpg" /> and <img src="22-7400953\91d27f1d-23fb-4e09-bcde-b4cf10c4211f.jpg" /> is solution of problem(NLP), or for any an accumulation <img src="22-7400953\40facc6e-5c6d-4b56-b6dc-6f75b5ea6f05.jpg" /> of sequence<img src="22-7400953\454c9334-1ceb-4564-9d94-765ff1a27cf7.jpg" />,<img src="22-7400953\7ae6afac-8a03-4983-b511-3a63644daa7e.jpg" /> is solution of problem (NLP).</p><p>Proof: Because the algorithm stops at its<img src="22-7400953\c3c6347f-21ed-4416-935c-f3acba3663df.jpg" />, then we have</p><disp-formula id="scirp.24125-formula70092"><label>(15)</label><graphic position="anchor" xlink:href="22-7400953\2c9b7dcf-781d-4e8e-b27d-a42fabb4812e.jpg"  xlink:type="simple"/></disp-formula><p>for<img src="22-7400953\d68f1bd9-6560-450c-9a0c-d0169bc48f27.jpg" />, <img src="22-7400953\88cc15c7-e4a9-4872-8873-87840018ce0a.jpg" />, it is easy to see that for any <img src="22-7400953\4741617d-ba01-42e1-84a5-47ebc727d622.jpg" /> have</p><p><img src="22-7400953\ce8b0ed6-83b0-4a1e-ab8e-a6498d5e6fc8.jpg" /></p><p>It is from Step 2 of the algorithm that we have</p><disp-formula id="scirp.24125-formula70093"><label>(16)</label><graphic position="anchor" xlink:href="22-7400953\ce9d265b-4b07-4216-bbe8-f9e019ac04d5.jpg"  xlink:type="simple"/></disp-formula><p>putting (15) (16) into (10) or (13), we can obtain</p><p><img src="22-7400953\9222fd56-b089-4e8d-9096-68a3a8f00b8a.jpg" /></p><p>for<img src="22-7400953\ecb7ac5b-d215-4c07-bb2f-7e6701082b46.jpg" />, <img src="22-7400953\c66cc703-6369-4955-bbee-30d0488f4453.jpg" />, according to definition of</p><p><img src="22-7400953\7f9a3ec3-ee16-42c0-b099-1007a08b2ab5.jpg" />, we can obtain, that</p><p><img src="22-7400953\6cef0ef8-20a8-4781-8620-25b1d7fbbbc0.jpg" /></p><p>First part of the theorem holds，<img src="22-7400953\9e2601c4-613f-4110-8438-9e6f87e2e7de.jpg" /> is solution of problem (NLP).</p><p>On the other hand, if the algorithm is not stop at<img src="22-7400953\2ecd7e06-d875-44e2-8a1a-195c374f967a.jpg" />, for any accumulation point <img src="22-7400953\5a3789df-2d9e-4ed4-88ee-eec0395fe8ff.jpg" /> of sequence<img src="22-7400953\e385c402-f4e7-4502-a9ba-403b1bb179bf.jpg" />, from theorem 4.1, we can obtain, for any positive number C, that</p><p><img src="22-7400953\9782b4a5-d9a9-414d-af4d-01fea810f4d4.jpg" /></p><p>for any<img src="22-7400953\3417fe50-4d34-4718-bb6e-925507010ae7.jpg" />, have</p><p><img src="22-7400953\6e2efefd-acf6-485d-932d-8ba742d7d164.jpg" /></p><p>Let<img src="22-7400953\73032011-1304-4163-9684-2ae6c592adab.jpg" />,have</p><p><img src="22-7400953\bbf5970d-fa5e-482b-b546-a01dad9d2a92.jpg" /></p><p>Clearly, second part of the theorem holds. <img src="22-7400953\87a9fd6d-f857-4bfc-9a0a-d6946a130eb1.jpg" />is solution of problem (NLP).</p></sec><sec id="s5"><title>5. Conclusion</title><p>A new Lagrange multiplier function with 4-piecewise linear NCP function is proposed in this paper which has a nice equivalence between its solution and solution of original problem. We can solve it to obtain solution of original constrained problem, the algorithm corresponding with it be endowed with convergence.</p></sec><sec id="s6"><title>6. Acknowledgements</title><p>This paper was partially supported by the NNSF of China under Grant No. 10971053, and NNSF of Henan under Grant No. 094300510050.</p></sec><sec id="s7"><title>REFERENCES</title></sec><sec id="s8"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.24125-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">G. Pillo and L. Grippo, “A New Class of Augmented Lagrangian in Nonlinear Programming,” SIAM Journal on Control and Optimization, Vol. 17, No. 5, 1979, pp. 617-628.</mixed-citation></ref><ref id="scirp.24125-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">G. Pillo and L. Grippo, “An Augmented Lageangian function for Inequlity Constrains in Nonlinear Programming Problems,” Optimization Theorem and Application, Vol. 36, No. 4, 1982, pp. 495-520. 
doi:10.1007/BF00940544</mixed-citation></ref><ref id="scirp.24125-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">G. Pillo and S. Lucidi, “An augmented Lagrange Function with Improved Exactness Properties,” SIAM Journal on Control and Optimization, Vol. 12, No. 2, 2001, pp. 376-406.</mixed-citation></ref><ref id="scirp.24125-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">D. G. Pu and Z. Jin, “A New Lagrange Multiplier Method,” Tongji University Journal, Vol. 9, No. 38, 2010, pp. 1384-1391.</mixed-citation></ref><ref id="scirp.24125-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">K. D. Li and D. G. Pu, “A class of New Lagrange Multiplier Methods,” OR Transaction, Vol. 10, No. 4, 2006, pp. 9-24.</mixed-citation></ref><ref id="scirp.24125-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">X.-W. Du and Y.-L. Shang, “Exact Augmented Lagrange Function for Nonlinear Programming Problems with Inequality and Equality Constraints,” Applied Mathematics and Mechanics, Vol. 26, No. 12, 2005, pp. 1651-1656.</mixed-citation></ref><ref id="scirp.24125-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">K. Su, “Trust-Region Filter Method with NCP Function,” Journal of Systems Science and Mathematical Science, Vol. 28, No. 12, 2008, pp. 1525-1534</mixed-citation></ref><ref id="scirp.24125-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">D. G. Pu, “A New Filter Method,” OR Transaction, Vol. 1, No. 15, 2011, pp. 46-58.</mixed-citation></ref><ref id="scirp.24125-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">W. Hua and D. G. Pu, “A Modified QP-Free Feasible Method,” Applied Mathmatical Modelling, Vol. 35, No. 4, 2011, pp. 1696-1708. doi:10.1016/j.apm.2010.10.002</mixed-citation></ref><ref id="scirp.24125-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">K. Li and D. G. Pu, “Filter QP-Free Method with 3-Piece Wise Linear NCP Function,” OR Transaction, Vol. 12, No. 2, 2008, pp. 49-57.</mixed-citation></ref><ref id="scirp.24125-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Y. X. Yuan, “Numerical Method for Nonlinear Programming,” Shanghai Science and Technology Press, Shanghai, 1993.</mixed-citation></ref><ref id="scirp.24125-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">B. L. Chen, “The Theory of Optimization and Calculate,” 2nd Edition, Tsinghua University Press, Beijing, 2005.</mixed-citation></ref></ref-list></back></article>