A Regularized Newton Method with Correction for Unconstrained Convex Optimization
ABSTRACT
In this paper, we present a regularized Newton method (M-RNM) with correction for minimizing a convex function whose Hessian matrices may be singular. At every iteration, not only a RNM step is computed but also two correction steps are computed. We show that if the objective function is LC2, then the method posses globally convergent. Numerical results show that the new algorithm performs very well.

Received 17 January 2016; accepted 12 March 2016; published 15 March 2016

1. Introduction

We consider the unconstrained optimization problem [1] -[3]

(1.1)

where is twice continuously differentiable, whose gradient and Hessian are denoted by and respectively. Throughout this paper, we assume that the solution set of (1.1) is nonempty and denoted by X, and in all cases refers to the 2-norm.

It is well known that is convex if and only if is symmetric positive semidefinite for all. Moreover, if is convex, then if and only if x is a solution of the system of nonlinear equations

(1.2)

Hence, we could get the minimizer of by solving (1.2) [4] -[8] . The Newton method is one of efficient solution methods. At every iteration, it computes the trial step

(1.3)

where and. As we know, if is Lipschitz continuous and nonsingular at the solution, then the Newton method has quadratic convergence. However, this method has an obvious disadvantage when the is singular or near singular. In this case, we may compute the Moore-Penrose step [7] . But the computation of the singular value decomposition to obtain is sometimes prohibitive. Hence, computing a direction that is close to may be a good idea.

To overcome the difficulty caused by the possible singularity of, [9] proposed a regularized Newton method, where the trial step is the solution of the linear equations

(1.4)

where I is the identity matrix. is a positive parameter which is updated from iteration to iteration.

Now we need to consider another question, “how to choose the regularized parameter?” Yamashita and

Fukushima [10] chose and showed that the regularized Newton method has quadratic convergence under the local error bound condition which is weaker than nonsingularity. Fan and Yuan [11] took

with and showed that the Levenberg-Marqulardt method preserves the quadratic convergence. Numerial results ([12] [13] ) show that the choice of performs more stable and preferable.

Inspired by the regularized Newton method [13] with correction for nonlinear equations, we propose a modified regularized Newton method in this paper. At every iteration, the modified regularized Newton method first solves the linear equations

(1.5)

to obtain the Newton step, where is updated from iteration to iteration, and solves the linear equations

(1.6)

to obtain the approximate Newton step.

It is easy to see

(1.7)

Then it solves the linear equations

(1.8)

to obtain the approximate Newton step.

The aim of this paper is to study the convergence properties of the above modified regularized Newton method and do a numerical experiment to test its efficiency.

The paper is organized as follows. In Section 2, we present a new regularized Newton algorithm with correction by trust region technique, and then prove the global convergence of the new algorithm under some suitable conditions. In Section 3, we test the regularized Newton algorithm with correction and compared it with a regularized Newton method. Finally, we conclude the paper in Section 4.

2. The Algorithm and Its Global Convergence

In this section, we first present the new modified regularized Newton algorithm by using trust region technique, then prove the global convergence. First, we give the modified regularized Newton algorithm.

Let and be given by (1.6) and (1.8), respectively. Since the matrix is symmetric and positive definite, is a descent direction of at, however may not be. Hence we prefer to use a trust region technique.

Define the actual reduction of at the kth iteration as

(2.1)

Note that the regularization step is the minimizer of the convex minimization problem

If we let

then it can be proved [4] that is also a solution of the trust region problem

By the famous result given by Powell in [14] , we know that

(2.2)

By some simple calculations, we deduce from (1.7) that

so, we have

(2.3)

Similar to, is not only the minimizer of the problem

but also a solution to the trust region problem

where. Therefore we also have

(2.4)

Based on the inequalities (2.2), (2.3) and (2.4), it is reasonable for us to define the new predicted reduction as

(2.5)

which satisfies

(2.6)

The ratio of the actual reduction to the predicted reduction

(2.7)

plays a key role in deciding whether to accept the trial step and how to adjust the regularized parameter.

The regularized Newton algorithm with correction for unconstrained convex optimization problems is stated as follows.

Algorithm 2.1

Step 1. Given, , ,. Set.

Step 2. If, then stop.

Step 3. Compute.

Solve

(2.8)

to obtain.

Solve

(2.9)

to obtain and set

Solve

(2.10)

to obtain and set

Step 4. Compute. Set

(2.11)

Step 5. Choose as

(2.12)

Set and go step 2.

Before discussing the global convergence of the algorithm above, we make the following assumption.

Assumption 2.1. and are both Lipschitz continuous, that is, there exists a constant, such that

(2.13)

and

(2.14)

It follows from (2.14) that

(2.15)

The following lemma given below shows the relationship between the positive semidefinite matrix and symmetric positive semidefinite matrix.

Lemma 2.1. A real-valued matrix A is positive semidefinite if and only if is positive semidefinite.

Proof. See [4] . ♢

Next, we give the bounds of a positive definite matrix and its inverse.

Lemma 2.2. Suppose A is positive semidefinite. Then,

and

hold for any.

Proof. See [13] . ♢

Theorem 2.1. Under the conditions of Assumption 2.1, if f is bounded below, then Algorithm 2.1 terminates in finite iterations or satisfies

(2.16)

Proof. We prove by contradiction. If the theorem is not true, then there exists a positive and an integer such that

(2.17)

Without loss of generality, we can suppose. Set. Then

Now we will analysis in two cases whether T is finite or not.

Case (1): T is finite. Then there exists an integer such that

By (2.11), we have

Therefore by (2.12) and (2.17), we deduce

(2.18)

Since, , we get from (2.8) and (2.18) that

(2.19)

Duo to (1.7), we get

From (2.10), we obtain

(2.20)

where is a positive constant.

It follows from (2.1) and (2.5) that

(2.21)

Moreover, from (2.6), (2.17), (2.13) and (2.19), we have

(2.22)

for sufficiently large k.

Duo to (2.21) and (2.22), we get

(2.23)

which implies that. Hence, there exists positive constant such that, holds for all large k, which contradicts to (2.18).

Case (2): T is infinite. Then we have from (2.6) and (2.17) that

(2.24)

which implies that

(2.25)

The above equality together with the updating rule of (2.12) means

(2.26)

Similar to (2.20), it follows from (2.25) and (2.26) that

for some positive constant. Then we have

This equality together with (2.24) yields

which implies that

(2.27)

It follows from (2.8), (2.27), (2.26) and (2.20) that

(2.28)

Since from (2.8), we have from (2.13), (2.17) and (2.28) that

which means

(2.29)

By the same analysis as (2.23) we know that

(2.30)

Hence, there exists a positive constant such that holds for all sufficiently large k, which gives a contradiction to (2.29). The proof is completed. ♢

3. Numerical Experiments

In this section, we test the performance of Algorithm 2.1 on the unconstrained nonlinear optimization problem, and compared it with a regularized Newton method without correction. The function to be minimized is

(3.1)

where are constants. It is clear that function is convex and the minimizer set of is

The Hessian is given by

where,. Matrix is positive semidefinite for all x, but

singular as the sum of every column is zero. Since the Hessian is always singular, the Newton method cannot be used to solve nonlinear Equations (1.2). But by using the regularization technique, both regularized Newton method and Algorithm 2.1 work quite well.

The aims of the experiments are as follows: to check whether Algorithm 2.1 converges quadratically as stated in Section 3 and also to see how well the technique of correction works. We set, , , , , and for Algorithm 2.1.

Table 1 reports the norms of at every iteration when, ,

and. Algorithm 2.1 only take four iterations to obtain the minimizer of; decreases very quickly. The results show the sequence quadratic convergence. The iteration is as follows

We may observe that the whole sequence converges to.

We also ran the regularized Newton algorithm (RNA) without correction, that is, we do not solve the linear equations (2.9)-(2.10) and just set the solution of (2.8) to be the trial step. Then, we tested the regularized Newton algorithm without correction and Algorithm 2.1 for various of n, and different choices of the starting point. The results are listed in Table 2.: the selected value of; Dim: the dimension n of the problem;: the ith element; niter: the number of iterations required;: the final value of;: the final value of. We use as the stopping criterion.

Table 1. Results of Algorithm 2.1 to test quadratic convergence.

Table 2. Results of RNA and Algorithm 2.1.

Moreover, we can see for the same, n and, the number of iterations of Algorithm 2.1 is always less than that of RNA. And the correction term does help to improve RNA when the initial point is far away from the minimizer. These facts indicate that the introduction of correction is really useful and could accelerate the convergence of the regularized Newton method.

4. Concluding Remarks

In this paper, we propose a regularized Newton method with correction for unconstrained convex optimization. At every iteration, not only a RNM step is computed but also two correction steps are computed which make use of the previous available Jacobian instead of computing the new Jacobian. Numerical experiments suggest that the introduction of correction is really useful.

Acknowledgements

This research is supported by the National Natural Science Foundation of China (11426155) and the Hujiang Foundation of China (B14005).

NOTES

*Corresponding author.

Cite this paper
Li, L. , Qin, M. and Wang, H. (2016) A Regularized Newton Method with Correction for Unconstrained Convex Optimization. Open Journal of Optimization, 5, 44-52. doi: 10.4236/ojop.2016.51006.
References
[1]   Yang, W.W., Yang, Y.T., Zhang, C.H. and Cao, M.Y. (2013) A Newton-Like Trust Region Method for Large-Scale Unconstrained Nonconvex Minimization. Abstract and Applied Analysis, 2013.

[2]   Polyak, R.A. (2009) Regularized Newton Method for Unconstrained Convex Optimization. Mathematical Programming Series B, 120, 125-145.
http://dx.doi.org/10.1007/s10107-007-0143-3

[3]   Shen, C.G., Chen, X.D. and Liang, Y.M. (2012) A Regularized Newton Method for Degenerate Unconstrained Optimization Problems. Optimization Letters, 6, 1913-1933.
http://dx.doi.org/10.1007/s11590-011-0386-z

[4]   Sun, W. and Yuan, Y. (2006) Optimization Theory and Methods. Springer Science and Business Media, LLC, New York.

[5]   Zhou, W. and Li, D. (2008) A Globally Convergent BFGS Method for Nonlinear Monotone Equations without Any merit Functions. Mathematics of Computation, 77, 2231-2240.
http://dx.doi.org/10.1090/S0025-5718-08-02121-2

[6]   Argyros, I.K. and Hilout, S. (2012) On the Convergence of Damped Newton Method. Applied Mathematics and Computation, 219, 2808-2824.
http://dx.doi.org/10.1016/j.amc.2012.09.011

[7]   Nashed, M.Z. and Chen, X. (1993) Convergence of Newton-Like Method for Singular Operator Equations Using Outer Inverses. Numerische Mathematik, 66, 235-257.
http://dx.doi.org/10.1007/BF01385696

[8]   Kelley, C.T. (1999) Iterative Methods for Optimization. In: Frontiers in Applied Mathematics, Vol. 18, SIAM, Philadelphia, 2.
http://dx.doi.org/10.1137/1.9781611970920

[9]   Sun, D. (1999) A Regularization Newton Method for Solving Nonlinear Complementarity Problems. Applied Mathematics and Optimization, 40, 315-339.
http://dx.doi.org/10.1007/s002459900128

[10]   Li, D.H., Fukushima, M., Qi, L. and Yamashita, N. (2004) Regularized Newton Methods for Convex Minimization Problems with Singular Solutions. Computational Optimization and Applications, 28, 131-147.
http://dx.doi.org/10.1023/B:COAP.0000026881.96694.32

[11]   Fan, J.Y. and Yuan, Y.X. (2005) On the Quadratic Convergence of the Levenberg-Marquardt Method without Nonsingularity Assumption. Computing, 74, 23-39.
http://dx.doi.org/10.1007/s00607-004-0083-1

[12]   Fan, J.Y. and Pan, J.Y. (2009) A Note on the Levenberge-Marquardt Parameter. Applied Mathematics and Computation, 207, 351-165.
http://dx.doi.org/10.1016/j.amc.2008.10.056

[13]   Fan, J.Y. and Yuan, Y.X. (2014) A Regularized Newton Method for Monotone Nonlinear Equations and Its Application. Optimization Methods and Software, 29, 102-119.
http://dx.doi.org/10.1080/10556788.2012.746344

[14]   Powell, M.J.D. (1975) Convergence Properties of a Class of Minimization Algorithms. In: Mangasarian, O.L., Meyer, R.R. and Robinson, S.M., Eds., Nonlinear Programming, Vol. 2, Academic Press, New York, 1-27.
http://dx.doi.org/10.1016/b978-0-12-468650-2.50005-5

Top