Next Article in Journal
Recognition of Stator Winding Inter-Turn Fault in Interior-Mount LSPMSM Using Acoustic Signals
Next Article in Special Issue
An Efficient Parallel Extragradient Method for Systems of Variational Inequalities Involving Fixed Points of Demicontractive Mappings
Previous Article in Journal
Characterisation and Manipulation of Polarisation Response in Plasmonic and Magneto-Plasmonic Nanostructures and Metamaterials
Previous Article in Special Issue
Implicit Extragradient-Like Method for Fixed Point Problems and Variational Inclusion Problems in a Banach Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Generalized Trust-Region Sub-Problem with Additional Linear Inequality Constraints—Two Convex Quadratic Relaxations and Strong Duality

1
Department of Mathematics, Statistics and Physics, Qatar University, Doha 2713, Qatar
2
Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Guilan, Rasht 4199613776, Iran
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(8), 1369; https://doi.org/10.3390/sym12081369
Submission received: 26 July 2020 / Revised: 11 August 2020 / Accepted: 12 August 2020 / Published: 17 August 2020
(This article belongs to the Special Issue Symmetry in Nonlinear Functional Analysis and Optimization Theory)

Abstract

:
In this paper, we study the problem of minimizing a general quadratic function subject to a quadratic inequality constraint with a fixed number of additional linear inequality constraints. Under a regularity condition, we first introduce two convex quadratic relaxations (CQRs), under two different conditions, that are minimizing a linear objective function over two convex quadratic constraints with additional linear inequality constraints. Then, we discuss cases where the CQRs return the optimal solution of the problem, revealing new conditions under which the underlying problem admits strong Lagrangian duality and enjoys exact semidefinite optimization relaxation. Finally, under the given sufficient conditions, we present necessary and sufficient conditions for global optimality of the problem and obtain a form of S-lemma for a system of two quadratic and a fixed number of linear inequalities.

1. Introduction

Consider the following generalized trust-region subproblem with additional linear inequality constraints:
p * : = inf q 1 ( x ) : = x T A x + 2 a T x q 2 ( x ) : = x T B x + 2 b T x + β 0 , c i T x d i , i = 1 , , m ,
where A , B R n × n are symmetric matrices but not necessarily positive semidefinite, a , b , c i R n and β , d i R , i = 1 , , m . Model problem (1) arises for instance in nonlinear optimization problems with linear inequality constraints when the trust-region methods are applied to solve [1] or in general nonlinear programming for computing search directions when sequential quadratic programming methods are employed [2] and polyhedral data uncertainty [3,4] or in robust optimization problems under matrix norm [5].
Although some special cases of problem (1) are polynomially solvable, it is a difficult problem in general. When B = I , b = 0 , β < 0 and ( c i , d i ) = ( 0 , 0 ) , i = 1 , , m , problem (1) is known as the trust-region subproblem (TRS) which is fundamental in the trust-region methods for unconstrained optimization [1]. Several efficient algorithms have been introduced for TRS in the literature [6,7,8,9,10]. Specifically, TRS has many nice properties such as exact semidefinite optimization (SDO) relaxation and strong duality [7,11]. In general, these important features do not hold for the following extended trust-region subproblem (eTRS) with a fixed number of linear inequality constraints, even for eTRS with one linear inequality constraint [12,13,14]:
min x T A x + 2 a T x | | x | | 2 + β 0 , c i T x d i , i = 1 , , m .
Jeyakumar and Li [15] proved that the SDO-relaxation of eTRS is exact whenever the dimension condition
dim Null ( A λ min ( A ) I ) s + 1 ,
is satisfied where s = dim span { c 1 , , c m } and λ min ( A ) denotes the smallest eigenvalue of A. They also derived necessary and sufficient optimality conditions for eTRS under the Slater condition and the dimension condition (2). Later, in Reference [16], the authors obtained the exactness of SDO-relaxation of eTRS under the following condition
Rank ( [ A λ min ( A ) I c 1 c m ] ) n 1 ,
which is more general than the dimension condition of Jeyakumar and Li. Most recently, in Reference [17], the authors have studied variants of TRS having additional conic constraints,
min q 1 ( x ) = x T A x + 2 a T x | | x | | 1 , H x h K ,
where A R n × n is a symmetric matrix, H R m × n , h R m and K R m is a closed convex cone. Assuming λ min ( A ) < 0 , they introduced the following convex relaxation for (4):
min q 1 ( x ) + λ min ( A ) ( 1 | | x | | 2 ) | | x | | 1 , H x h K ,
and showed that this convex relaxation is exact if there exists nonzero z Null ( A λ min ( A ) I ) such that H z K and a T z 0 .
When ( c i , d i ) = ( 0 , 0 ) , i = 1 , , m , problem (1) reduces to the generalized trust-region subproblem (GTRS) [18]. The GTRS has been well studied in the literature and several methods have been proposed to solve it under various assumptions [11,18,19,20,21,22,23,24,25]. It has strong duality and exact SDO-relaxation under the Slater condition [11,22]. Under the assumption that the Hessian of quadratic functions are simultaneously diagonalizable, Ben-Tal and den Hertog [26] proved that GTRS admits a second order cone programming (SOCP) reformulation. Under the same assumption, they generalized this result to GTRS with two quadratic inequality constraints. They showed that under certain additional conditions, the optimal solution of the original problem can be recovered from the optimal solution of the SOCP relaxation [26]. However, in Reference [26], it has been illustrated that even in the simplest case where B 0 and the second constraint is linear (eTRS with m = 1 ), the SOCP relaxation may not be exact. After that, Locatelli [27] extended the SOCP relaxation to eTRS and gave conditions under which the SOCP relaxation is tight. Moreover, in Reference [28], it has been shown that the SOCP relaxation of eTRS and its SDO-relaxation are equivalent. Also, through this equivalence, new conditions are introduced that ensure the exactness of the SDO-relaxation of eTRS and are more general than the condition (3) [28]. Most recently, in Reference [29], the authors have proposed a new convex quadratic reformulation for GTRS that minimizes a linear objective function subject to two convex quadratic constraints. They also have shown that the optimal solution of GTRS can be recovered from the optimal solution of the new reformulation.
The main contributions of this paper are as follows:
(i)
Under a regularity condition, we present two convex quadratic relaxations (CQRs) under two different conditions for problem (1) that minimize a linear objective function subject to two convex quadratic constraints with a fixed number of additional linear inequality constraints. Our CQRs are inspired by the one proposed for GTRS in Reference [29]. Then we derive sufficient conditions under which problem (1) is equivalent to exactly one of the CQRs and the optimal solution of (1) can be recovered from an optimal solution of the CQRs. These sufficient conditions are easy to verify and involve only one (any) optimal solution of CQRs. We also show that under these conditions the attractive features of GTRS such as strong Lagrangian duality and exact SDO-relaxation hold for (1). It should be noted that in the case of GTRS, the CQRs reduce to the ones proposed in Reference [29]. The CQRs are always exact for GTRS but in the presence of linear constraints, they are not exact in general.
(ii)
Exploiting the results in (i), we also derive sufficient conditions that are expressed in terms of the data of the model problem (1) for exactness of the CQRs, strong Lagrangian duality and consequently for tightness of the SDO-relaxation. In the case of eTRS, these sufficient conditions reduce to the one presented in Reference [17] that is the existing best results in the literature. As a consequence, we present necessary and sufficient conditions for global optimality of problem (P) under the new condition together with the Slater condition. We also obtain a form of S-lemma for the system of two quadratic and a fixed number of linear inequalities.
(iii)
The sufficient conditions in (i) and (ii) ensure the exactness of the CQRs and the SDO-relaxation of problem (1). It is worth noting that solving large-scale semidefinite problems is still an intractable task. In contrast, the CQRs are significantly more tractable than SDOs, and advanced commercial software is available to solve them [30].
The rest of the paper is organized as follows—in Section 2, we introduce the CQRs and discuss when and how one can obtain an optimal solution of problem (1) from an optimal solution of the CQRs, revealing new sufficient conditions for strong duality of problem (1). In Section 3, we use the results in Section 2 to derive sufficient conditions based on the data of the original problem for exactness of the CQRs, strong Lagrangian duality and exact SDO-relaxation. We also present necessary and sufficient conditions for global optimality of problem (1) and an application of strong duality to S-lemma.
Notation 1.
Throughout this paper, for a symmetric matrix A, A 0 ( A 0 ) denotes that A is positive definite (positive semidefinite). Moreover, Null ( A ) and Rank ( A ) denote its Null space and Rank. Finally, A B : = trac ( A B ) is the usual matrix inner product of two symmetric matrices A and B.

2. Convex Quadratic Relaxation, Global Minimization and Strong Duality

In this section, following the idea of Reference [29], first we present two new convex quadratic relaxations for problem (1) under two different conditions. Then we discuss cases where problem (1) is equivalent to one of the CQRs and its global optimal solution can be obtained from an optimal solution of the CQRs. This equivalence reveals new conditions under which problem (1) enjoys strong Lagrangian duality. We start by considering the following assumptions
Assumption 1.
There exists λ ^ 0 such that A + λ ^ B 0 .
Assumption 2.
The Slater condition holds for problem (1), that is, there exists x ^ with q 2 ( x ^ ) < 0 , c i T x ^ d i for i = 1 , , m .
Assumptions 1 and 2 ensure that problem (1) is solvable as proved in the following lemma.
Lemma 1.
Suppose that Assumptions 1 and 2 hold. Then problem (1) has an optimal solution, that is, the infimum in (1) is always attainable.
Proof. 
See Appendix A.1. □
Assumption 1 implies that matrices A and B are simultaneously diagonalizable by congruence [31], that is, there exists an invertible matrix S and diagonal matrices D = diag ( α 1 , , α n ) and E = diag ( e 1 , , e n ) such that S T A S = D and S T B S = E . Define I P S D : = { λ 0 | A + λ B 0 } . By Assumption 1, I P S D . It is easy to see that I P S D = [ λ 1 , λ 2 ] where
λ 1 = max { α i e i | e i > 0 } , λ 2 = min { α i e i | e i < 0 } .
We have two cases for the set I P S D as follows, where λ ^ 1 = max { 0 , λ 1 } :
Condition 1.
I P S D = [ λ ^ 1 , λ 2 ] as long as B is not positive semidefinite.
Condition 2.
I P S D = [ λ ^ 1 , ) as long as B is positive semidefinite.
In the sequel, we introduce two new CQRs (6) and (7) corresponding to Condition 1 and Condition 2, respectively, by defining h 1 ( x ) = q 1 ( x ) + λ ^ 1 q 2 ( x ) and h 2 ( x ) = q 1 ( x ) + λ 2 q 2 ( x ) :
p 1 * : = inf x , t t h 1 ( x ) t , h 2 ( x ) t , c i T x d i , i = 1 , , m ,
and
p 2 * : = inf x h 1 ( x ) q 2 ( x ) 0 , c i T x d i , i = 1 , , m .
In the case where both A and B are positive semidefinite, problem (7) is equal to problem (1) that is already a convex quadratic problem. Hence, from now on, we suppose that at least one of q 1 ( x ) and q 2 ( x ) is nonconvex. As it is shown in the proof of Lemma 1, problem (1) is equivalent to its epigraph as follows:
p * = inf t , x t q 1 ( x ) t , q 2 ( x ) 0 , c i T x d i , i = 1 , , m .
Problems (6) and (7) are both convex. Under Condition 1, problem (6) is a convex relaxation of problem (8) and hence p 1 * p * . To see this, let x be a feasible solution of problem (8). Since q 2 ( x ) 0 and λ ^ 1 , λ 2 0 , then
h 1 ( x ) = q 1 ( x ) + λ ^ 1 q 2 ( x ) q 1 ( x ) t ,
h 2 ( x ) = q 1 ( x ) + λ 2 q 2 ( x ) q 1 ( x ) t .
Therefore, the feasible region of problem (6) contains that of problem (8) and since the two problems have the same objective function, then p 1 * p * . Next, suppose that Condition 2 holds. Problems (7) and (1) have the same feasible region and since h 1 ( x ) q 1 ( x ) , we have p 2 * p * . The following lemma states that problems (6) and (7) are bounded from below and their optimal values are attained.
Lemma 2.
Under Assumptions 1 and 2, problems (6) and (7) are bounded from below and their optimal values are attained.
Proof. 
See Appendix A.2. □
In the case of GTRS, the CQRs (6) and (7) reduce to the ones introduced in Reference [29]. Under Assumptions 1 and 2, the CQRs (6) and (7) are always exact for GTRS [29] while in the presence of linear constraints, they are not exact in general as illustrated in the following example.
Example 1.
Consider the following GTRS:
p * = min q 1 ( x ) : = 1 2 x 2 1 2 x q 2 ( x ) : = x 2 1 0 .
The CQR relaxation of (9) is
p 2 * = min h 1 ( x ) = 1 2 x 1 2 q 2 ( x ) = x 2 1 0 .
It is easy to see that the CQR (10) is exact, that is, p * = p 2 * = 1 . If we add the linear constraint x 1 2 to problem (9), then the resulting CQR is not exact.
In the following theorems, we state cases where the nonconvex problem (1) is equivalent to one of the CQR (6) or (7), that is, p * = p 1 * or p * = p 2 * and the optimal solution of (1) can be obtained from an optimal solution of the CQRs. In particular, this equivalence results in sufficient conditions for the strong Lagrangian duality for problem (1).
Theorem 1.
Suppose that Assumptions 1, 2 and Condition 1 hold, ( x * , t * ) is an optimal solution of problem (6) and one of the following holds:
(1) 
h 1 ( x * ) = h 2 ( x * ) = t * .
(2) 
h 1 ( x * ) < t * and there exists nonzero z Null ( A + λ 2 B ) such that ( a + λ 2 b ) T z = 0 and c i T z 0 for i = 1 , , m .
(3) 
h 2 ( x * ) < t * and λ ^ 1 = 0 .
(4) 
h 2 ( x * ) < t * , λ ^ 1 > 0 and there exists nonzero z Null ( A + λ 1 B ) such that ( a + λ 1 b ) T z = 0 and c i T z 0 for i = 1 , , m .
Then problem (1) is equivalent to (6), that is, p * = p 1 * , strong duality holds for (1) and its Lagrangian dual problem is solvable. Also, in case (1) or (3), x * is an optimal solution to (1) and in case (2) or (4), x ¯ * : = x * + α * z is an optimal solution for (1) where α * is the positive root of the quadratic equation q 2 ( x * + α z ) = 0 .
Proof. 
See Appendix A.3. □
Theorem 2.
Suppose that Assumptions 1, 2 and Condition 2 hold, x * is an optimal solution of problem (7) and one of the following holds:
(1) 
q 2 ( x * ) = 0 .
(2) 
q 2 ( x * ) < 0 and there exists nonzero z Null ( A + λ ^ 1 B ) such that ( a + λ ^ 1 b ) T z = 0 and c i T z 0 for i = 1 , , m .
Then problem (1) is equivalent to (7), that is, p * = p 2 * , strong duality holds for (1) and its Lagrangian dual problem is solvable. Also, in cases (1), x * is an optimal solution to (1) and in case (2), x ¯ * : = x * + α * z is optimal for (1) where α * is the positive root of the quadratic equation q 2 ( x * + α z ) = 0 .
Proof. 
See Appendix A.4. □

3. New Conditions for Strong Duality and Exact SDO-Relaxation

In the previous section, we derived sufficient conditions for exactness of the CQRs and strong Lagrangian duality of problem (1) based on an optimal solution of CQRs (6) and (7) (see Theorems 1 and 2). Here we exploit the results in Theorems 1 and 2 to derive new sufficient conditions that are expressed in terms of the data of the original problem for the exactness of the CQRs, strong Lagrangian duality and tightness of the SDO-relaxation of problem (P). Recall that we have assumed at least one of A and B is not positive semidefinite. Otherwise, by Assumption 2, problem (1) is a convex optimization problem that satisfies the Slater condition and hence, it has strong duality and exact SDO-relaxation.
The so-called SDO-relaxation of (1) is
p r * : = min M X M 0 X 0 , M i X 0 , i = 1 , , m , I 0 X = 1 , X 0 ,
where
M = A a a T 0 , M 0 = B b b T β , I 0 = O n × n O n × 1 O 1 × n 1 , M i = O n × n c i 2 c i T 2 d i , i = 1 , , m .
The dual of (11) is
d * : = max s M + i = 0 m y i M i s I 0 0 , y i 0 , i = 0 , , m ,
which is also the Lagrangian dual problem of (1). Note that by Assumption 1, problem (12) is strictly feasible and hence, d * = p r * . This together with the fact that d * p r * p * imply that the strong duality holds for (1) if and only if the SDO-relaxation for (1) is exact.
Consider 3.
Consider problem (1). We say that problem (1) satisfies Condition 3 whenever one of the following holds:
1. 
Condition 1 holds, λ ^ 1 = 0 and there exists nonzero z Null ( A + λ 2 B ) such that ( a + λ 2 b ) T z 0 and c i T z 0 for i = 1 , , m .
2. 
Condition 1 holds, λ ^ 1 > 0 , there exist nonzero z 1 Null ( A + λ 1 B ) and z 2 Null ( A + λ 2 B ) , such that ( a + λ 1 b ) T z 1 0 , ( a + λ 2 b ) T z 2 0 , c i T z 1 0 and c i T z 2 0 for i = 1 , , m .
3. 
Condition 2 holds and there exists nonzero z Null ( A + λ ^ 1 B ) such that ( a + λ ^ 1 b ) T z 0 and c i T z 0 for i = 1 , , m .
Lemma 3.
Suppose that Assumptions 1, 2 and Condition 3 hold for problem (1). Then the CQRs (6) and (7) are exact and problem (1) enjoys strong duality and exact SDO-relaxation.
Proof. 
Suppose that Condition 1 holds and let ( x * , t * ) be an optimal solution of (6). If h 1 ( x * ) = h 2 ( x * ) = t * , then by Theorem 1, the CQR (6) is exact, strong Lagrangian duality holds for problem (1) and the SDO-relaxation is exact. Otherwise, either h 1 ( x * ) < t * or h 2 ( x * ) < t * . Let h 1 ( x * ) < t * . We show that, in this case, for all z Null ( A + λ 2 B ) satisfying Condition 3, we have ( a + λ 2 b ) T z = 0 , implying that Item (2) in Theorem 1 holds and thus the CQR (6) is exact, strong Lagrangian duality holds for (1) and the SDO-relaxation is exact. To this end, suppose by contradiction that there exists z Null ( A + λ 2 B ) such that ( a + λ 2 b ) T z < 0 and c i T z 0 for i = 1 , , m . Set x ¯ * = x * + α * z where α * is the positive root of Equation (A17). We have q 2 ( x ¯ * ) = 0 , c i T x ¯ * 0 for i = 1 , , m and since ( a + λ 2 b ) T z < 0 ,
h 2 ( x ¯ * ) = x ¯ * T ( A + λ 2 B ) x ¯ * + 2 ( a + λ 2 b ) T x ¯ * + λ 2 β < h 2 ( x * ) = t * .
On the other hand, since q 2 ( x ¯ * ) = 0 , we have t ¯ : = h 1 ( x ¯ * ) = h 2 ( x ¯ * ) < t * . These mean that ( x ¯ * , t ¯ ) is a feasible solution of (6) with t ¯ < t * that contradicts the fact that ( x * , t * ) is optimal for (6). Next, suppose that h 2 ( x * ) < t * . If λ ^ 1 = 0 , then by Theorem 1, the CQR (6) is exact, strong Lagrangian duality holds for problem (1) and the SDO-relaxation is exact. Otherwise, a similar discussion as above proves the existence of vector z in Item (2) of Theorem 1. Similarly, we can prove the existence of vector z in Item (2) of Theorem 2 when Condition 2 holds. □
Remark 1.
Condition 3 can be verified easily by solving a linear programming problem. To verify condition
0 z Null ( A + λ 2 B ) s . t ( a + λ 2 b ) T z 0 and c i T z 0 , i = 1 , , m ,
it is sufficient to solve the following linear programming problem:
p ^ : = min ( a + λ 2 b ) T z ( A + λ 2 B ) z = 0 , c i T z 0 , i = 1 , , m .
Condition (13) holds if and only if problem (14) is either unbounded or has multiple optimal solutions. If problem (14) is unbounded, then obviously condition (13) holds. If problem (14) is bounded, then p ^ = 0 . In this case, since z = 0 is an optimal solution of (14), condition (13) holds if and only if (14) has multiple optimal solutions. The same discussion holds when λ 2 is replaced by λ 1 or λ ^ 1 .
Remark 2.
It is worth noting that the sufficient conditions for strong duality of problem (1) established in Theorems 1 and 2 are more general than Condition 3. The following is an example where Condition 3 does not hold, while the condition in Theorem 2 does.
Example 2.
Consider the following one-dimensional problem:
min q 1 ( x ) : = 1 2 x 2 1 2 x q 2 ( x ) : = x 2 1 0 , x 0 .
The CQR relaxation of (15) is
min h 1 ( x ) = 1 2 x 1 2 q 2 ( x ) = x 2 1 0 , x 0 .
The optimal solution of (16) is x * = 1 and q 2 ( x * ) = 0 . Hence, the sufficient condition in Theorem 2 is fulfilled. However, Condition 3 is not fulfilled. Moreover, it is easy to verify that strong duality holds for problem (15) and the SDO-relaxation and the CQR (16) are exact.
Remark 3.
For the case of eTRS which is a special case of (4), the convex relaxation (5) is equal to problem (7). Also, it is easy to see that in the case of eTRS, Condition 3 reduces to the condition in Reference [17].
We now present necessary and sufficient conditions for global optimality of (1) whenever Condition 3 and Assumptions 1 and 2 are satisfied.
Corollary 1.
For problem (1), suppose that Condition 3, Assumptions 1 and 2 hold. Let x * be a feasible solution of (1). Then x * is a global minimizer of (1) if and only if there exist nonnegative multipliers λ i * , i = 0 , , m such that the following conditions hold
( A + λ 0 * B ) x * = ( a + λ 0 * b + i = 1 m λ i * 2 c i ) ,
λ 0 * q 2 ( x * ) = 0 ,
λ i * ( c i T x * d i ) = 0 , i = 1 , , m ,
( A + λ 0 * B ) 0 .
Proof. 
Let x * be a global minimizer of (1). Recall that by Lemma 3, strong duality holds for problem (1). Suppose that ( λ 0 * , λ 1 * , , λ m * ) is an optimal solution of Lagrangian dual of (1) and d * denotes the dual optimal value. We have
p * = d * = min x { q 1 ( x ) + λ 0 * q 2 ( x ) + i = 1 m λ i * ( c i T x d i ) } q 1 ( x * ) + λ 0 * q 2 ( x * ) + i = 1 m λ i * ( c i T x * d i ) q 1 ( x * ) = p * ,
where the last inequality follows from λ i * 0 , i = 0 , , m and feasibility of x * . We conclude that the two inequalities in this chain hold with equality. Since the inequality in the second line is an equality, we conclude that x * is a minimizer of the minimization problem in the first line. This gives relations (17) and (20). Moreover, it follows from the last line that λ 0 * q 2 ( x * ) + i = 1 m λ i * ( c i T x * d i ) 0 which with the fact that each term in this sum is nonpositive, we obtain (18) and (19). Conversely, suppose that x * satisfies (17) to (20). We have the following chain of inequalities:
p * d * : = max λ i 0 , i = 0 , , m min { q 1 ( x ) + λ 0 q 2 ( x ) + i = 1 m λ i ( c i T x d i ) } min { q 1 ( x ) + λ 0 * q 2 ( x ) + i = 1 m λ i * ( c i T x d i ) } = q 1 ( x * ) + λ 0 * q 2 ( x * ) + i = 1 m λ i * ( c i T x * d i ) = q 1 ( x * ) p * ,
where the first inequality comes from weak duality property, the second equality follows from (17) and (20), the third equality follows from (18) and (19) and the last inequality comes from the fact that x * is a feasible solution of (1). Therefore, q 1 ( x * ) = p * and so x * solves (1). □
As a consequence of our strong duality result, we obtain the following form of the celebrated S-lemma [32] for a system of two quadratic and a fixed number of linear inequalities.
Lemma 4.
Let A , B R n × n , a , b , c i R n , β , γ , d i R , i = 1 , , m . Suppose that Assumptions 1, 2 and Condition 3 are satisfied. Then the following statements are equivalent:
(1) 
x T B x + 2 b T x + β 0 , c i T x d i , i = 1 , , m x T A x + 2 a T x + γ 0 .
(2) 
There exist λ i 0 , i = 0 , , m such that
( x T A x + 2 a T x + γ ) + λ 0 ( x T B x + 2 b T x + β ) + i = 1 m λ i ( c i T x d i ) 0 , x R n .
Proof. 
It is easy to see that ( 2 ) ( 1 ) . So, in the sequel, we show that ( 1 ) ( 2 ) . To see this, let ( 1 ) holds. Then,
min x { x T A x + 2 a T x | x T B x + 2 b T x + β 0 , c i T x d i , i = 1 , , m } γ .
By Lemma 3, we have strong duality for the above problem. Suppose that ( λ 0 , λ 1 , , λ m ) is the dual optimal solution. This means that
min x { x T A x + 2 a T x | x T B x + 2 b T x + β 0 , c i T x d i , i = 1 , , m } = min x { x T A x + 2 a T x + λ 0 ( x T B x + 2 b T x + β ) + i = 1 m λ i ( c i T x d i ) } γ ,
which implies (2). □

4. Conclusions

Under a regularity condition, we introduced two convex quadratic relaxations (CQRs) corresponding to two different conditions for model problem (1) that are the problems of minimizing a linear objective function over two convex quadratic constraints with additional linear inequality constraints. We presented sufficient conditions based on an optimal solution of the CQRs under which the problem (1) is equivalent to exactly one of the CQRs. We also showed that this equivalence reveals the strong duality holds for (1) and consequently problem (1) enjoys exact SDO-relaxation. We also derived new sufficient conditions based on the data of the model problem (1) for strong Lagrangian duality, exact SDO-relaxation and exact CQRs.
Possible subjects for future research direction would be finding more general conditions for exactness of CQRs (maybe even necessary and sufficient conditions), and other relaxations that are as simple as problems (6) and (7), but tight under more general conditions.

Author Contributions

Conceptualization, A.T. and M.S.; Formal analysis, M.S.; Funding acquisition, T.A.A.; Investigation, A.H.; Methodology, T.A.A. and A.H.; Project administration, T.A.A.; Writing—original draft, A.T.; Writing—review & editing, A.H. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to express their thanks to Qatar University for supporting their project under Grant NCBP-QUCP-CAS-2020-1.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Proof of Lemma 1

Consider epigraph reformulation of (1) as following:
p 0 * : = inf t , x t q 1 ( x ) t , q 2 ( x ) 0 , c i T x d i , i = 1 , , m .
First we show that the infimum in (A1) is attainable. To this end, note that under Assumptions 1 and 2, problem (1) without the linear constraints is bounded from below (see Theorem 5 of Reference [33]). This implies that (A1) is equivalent to the following problem:
p 0 * = inf t , x t q 1 ( x ) t , q 2 ( x ) 0 , c i T x d i , i = 1 , , m , t ^ t M ,
where t ^ is the optimal (infimum) value of q 1 ( x ) over the constraint q 2 ( x ) 0 and M is a sufficiently large constant. Let S denote the feasible region of problem (A2). The set S is closed and the objective function in problem (A2) is continuous. Therefore, to prove that the infimum in (A1) is attainable, it is sufficient to establish that S is bounded. Since t ^ t M , we only need to show that there exists M ^ > 0 such that | | x | | M ^ for all ( x , t ) S . To do so, let h ( x ) : = q 1 ( x ) + λ ^ q 2 ( x ) where λ ^ is the same as in Assumption 1. The function h is strictly convex and for any ( x , t ) S , we have h ( x ) M , implying that S is bounded. Next, we show that problem (A1) is equivalent to problem (1) and hence the infimum in (1) is always attainable. Let ( x * , t * ) be an optimal solution of (A1). Since q 2 ( x * ) 0 and c i T x * d i , i = 1 , , m , x * is feasible for (1) and since q 1 ( x * ) t * , we have p * t * = p 0 * . We show that p * = p 0 * . Suppose by contradiction that p * < p 0 * . Then by definition of infimum, there exists a feasible solution of (1), x, such that q 1 ( x ) < p 0 * . Set t = q 1 ( x ) . Then ( x , t ) is a feasible solution of (A1) with a smaller objective value that contradicts the fact that p 0 * is the optimal value of (A1). Therefore, p * = p 0 * and since q 1 ( x * ) = t * = p * , x * is the optimal solution of (1).

Appendix A.2. Proof of Lemma 2

Consider problem (6). By Theorem 2.9 of Reference [29], problem
inf x , t t h 1 ( x ) t , h 2 ( x ) t ,
is bounded from below and its optimal value is attained. This implies that problem (6) is equivalent to the following problem:
p 1 * = inf x , t t h 1 ( x ) t , h 2 ( x ) t , c i T x d i , i = 1 , , m , t * t M ,
where t * is the optimal value of (A3) and M is a sufficiently large constant. Let S denote the feasible region of (A4). The set S is closed and the objective function in (A4) is continuous. Therefore, to prove the assertion, it is sufficient to establish that S is bounded. Since t * t M , we only need to show there exists M ^ > 0 such that | | x | | M ^ for all ( x , t ) S . To this end, let h ( x ) : = α 1 h 1 ( x ) + α 2 h 2 ( x ) where α 1 , α 2 > 0 and α 1 + α 2 = 1 . The function h is strictly convex and for any ( x , t ) S , we have h ( x ) M . The proof for problem (7) is similar.

Appendix A.3. Proof of Theorem 1

Since problem (6) is convex and satisfies the Slater condition, there exist nonnegative multipliers μ 1 * , μ 2 * , s i * , i = 1 , , m , such that
( A + ( μ 1 * λ ^ 1 + μ 2 * λ 2 ) B ) x * = ( a + ( μ 1 * λ ^ 1 + μ 2 * λ 2 ) b + i = 1 m s i * 2 c i ) ,
μ 1 * ( h 1 ( x * ) t * ) = 0 ,
μ 2 * ( h 2 ( x * ) t * ) = 0 ,
s i * ( c i T x * d i ) = 0 , i = 1 , , m ,
μ 1 * + μ 2 * = 1 ,
h 1 ( x * ) t * ,
h 2 ( x * ) t * ,
c i T x * d i , i = 1 , , m .
It follows from (A10) and (A11) that there are three possible cases: (i) h 1 ( x * ) = h 2 ( x * ) = t * , (ii) h 1 ( x * ) < t * or (iii) h 2 ( x * ) < t * . In what follows, we discuss these cases and we consider the two possible cases h 2 ( x * ) < t * , λ 1 ^ = 0 and h 2 ( x * ) < t * , λ ^ 1 > 0 , separately.
(1)
h 1 ( x * ) = h 2 ( x * ) = t * implies that ( λ ^ 1 λ 2 ) ( x * T B x * + 2 b T x * + β ) = 0 . Furthermore, since λ ^ 1 λ 2 , we obtain
x * T B x * + 2 b T x * + β = 0 .
It follows from (A12) and (A13) that ( x * , t * ) is also feasible for (8) and since (6) is a relaxation of (8), then ( x * , t * ) solves (8), p * = p 1 * and thus x * solves (1). To prove strong duality, set λ * = μ 1 * λ ^ 1 + μ 2 * λ 2 . Since μ 1 * 0 , μ 2 * 0 and μ 1 * + μ 2 * = 1 , then λ * [ λ ^ 1 , λ 2 ] and thus
A + λ * B 0 .
Also, we have
p * d * : = max γ i 0 , i = 0 , , m min x q 1 ( x ) + γ 0 q 2 ( x ) + i = 1 m γ i ( c i T x d i ) min x q 1 ( x ) + λ * q 2 ( x ) + i = 1 m s i * ( c i T x d i ) = q 1 ( x * ) + λ * q 2 ( x * ) + i = 1 m s i * ( c i T x * d i ) = q 1 ( x * ) p * ,
where the first inequality follows from the weak duality property, the second equality follows from (A5) and (A14), the third equality follows from (A8) and (A13) and the last inequality follows from (A12) and (A13). Therefore, we have p * = d * , that is, the strong duality holds for problem (1) and the maximum in (A15) is attained.
(2)
In this case, μ 2 * = 1 and hence, h 2 ( x * ) = t * . Then h 1 ( x * ) < t * and h 2 ( x * ) = t * imply that ( λ ^ 1 λ 2 ) ( x * T B x * + 2 b T x * + β ) < 0 . Since λ ^ 1 < λ 2 , we obtain
x * T B x * + 2 b T x * + β > 0 .
By the assumption, there exists nonzero z Null ( A + λ 2 B ) such that ( a + λ 2 b ) T z = 0 and c i T z 0 for i = 1 , , m . Consider the following quadratic equation of variable α :
q 2 ( x * + α z ) = α 2 z T B z + 2 α ( z T B x * + b T z ) + x * T B x * + 2 b T x * + β = 0 .
The fact that z T B z < 0 (see Lemma 3.4 of Reference [22]) with (A16) imply that the above equation has a positive root α * . Set x ¯ * = x * + α * z . We have q 2 ( x ¯ * ) = 0 and since c i T z 0 , i = 1 , , m , we also have c i T x ¯ * d i , i = 1 , , m . Furthermore, since z Null ( A + λ 2 B ) and ( a + λ 2 b ) T z = 0 , we have
h 2 ( x ¯ * ) = x ¯ * T ( A + λ 2 B ) x ¯ * + 2 ( a + λ 2 b ) T x ¯ * + λ 2 β = h 2 ( x * ) = t * .
It follows from (A18) and q 2 ( x ¯ * ) = 0 that q 1 ( x ¯ * ) = t * and consequently h 1 ( x ¯ * ) = t * . These indicate that ( x ¯ * , t * ) is an optimal solution of (6) which is also feasible for (8). Since (6) is a relaxation of (8), ( x ¯ * , t * ) solves (8), p * = p 1 * and thus x ¯ * solves (1). The same approach as in part (1) can be applied to show that strong duality holds for (1) and the Lagrangian dual problem is solvable.
(3)
In this case, μ 1 * = 1 and hence, h 1 ( x * ) = t * . Also, h 2 ( x * ) < t * and h 1 ( x * ) = t * imply that ( λ 2 λ ^ 1 ) ( x * T B x * + 2 b T x * + β ) < 0 . Then λ 2 > λ ^ 1 results in
x * T B x * + 2 b T x * + β < 0 .
It follows from (A19) and (A12) that ( x * , t * ) is also feasible for (8) and since (6) is a relaxation of (8), x * solves (1) and p * = p 1 * . Then by setting λ * = 0 , the same approach as in part (1) can be applied to show that strong duality holds for (1) and the Lagrangian dual problem is solvable.
(4)
By the assumption, there exists nonzero z Null ( A + λ 1 B ) such that ( a + λ 1 b ) T z = 0 and c i T z 0 for i = 1 , , m . Consider the following quadratic equation of variable α :
q 2 ( x * + α z ) = α 2 z T B z + 2 α ( z T B x * + b T z ) + x * T B x * + 2 b T x * + β = 0 .
The fact that z T B z > 0 (see Lemma 3.4 of Reference [22]) with (A19) imply that the above equation has a positive root α * . Then following the same discussion as in part (2) where λ 2 is replaced by λ 1 and h 2 ( x ¯ * ) in (A18) is replaced by h 1 ( x ¯ * ) , it can be shown that x ¯ * : = x * + α z solves problem (1), p * = p 1 * , strong duality holds for problem (1) and the Lagrangian dual problem is solvable.

Appendix A.4. Proof of Theorem 2

Since we have assumed that q 1 ( x ) and q 2 ( x ) are not both convex, then λ ^ 1 > 0 . Moreover, problem (7) is convex, satisfies the Slater condition and by Lemma 2 is solvable. Let x * be an optimal solution of (7). Therefore, there exist nonnegative multipliers μ 1 * , s i * , i = 1 , , m , such that
( A + ( λ ^ 1 + μ 1 * ) B ) x * = ( a + ( μ 1 * + λ ^ 1 ) b + i = 1 m s i * 2 c i ) ,
μ 1 * q 2 ( x * ) = 0 ,
s i * ( c i T x * d i ) = 0 , i = 1 , , m ,
q 2 ( x * ) 0 ,
c i T x * d i , i = 1 , , m .
(1)
In this case, since q 2 ( x * ) = 0 and (7) is a relaxation of (1), then q 1 ( x * ) = p 2 * p * q 1 ( x * ) and consequently p * = p 2 * = q 1 ( x * ) . Next by setting λ * = μ 1 * + λ ^ 1 , the same approach as in part (1) of Theorem 1 can be applied to show that strong duality holds for problem (1) and the Lagrangian dual problem is solvable.
(2)
By the assumption, there exists nonzero z Null ( A + λ ^ 1 B ) such that ( a + λ ^ 1 b ) T z = 0 and c i T z 0 for i = 1 , , m . Consider the following quadratic equation of variable α :
q 2 ( x * + α z ) = α 2 z T B z + 2 α ( z T B x * + b T z ) + x * T B x * + 2 b T x * + β = 0 .
The fact that z T B z > 0 with q 2 ( x * ) < 0 implies that the above equation has a positive root α * . Set x ¯ * = x * + α * z . We have q 2 ( x ¯ * ) = 0 , c i T x ¯ * 0 for i = 1 , , m and
h 1 ( x ¯ * ) = x ¯ * T ( A + λ ^ 1 B ) x ¯ * + 2 ( a + λ ^ 1 b ) T x ¯ * + λ ^ 1 β = h 1 ( x * ) ,
since ( a + λ ^ 1 b ) T z = 0 and z Null ( A + λ ^ 1 B ) . These imply that x ¯ * is an optimal solution of (7). Moreover, since q 2 ( x ¯ * ) = 0 , we have q 1 ( x ¯ * ) = h 1 ( x ¯ * ) , implying that x ¯ * solves problem (1), p * = p 2 * . Now the same approach as in part (1) shows that strong duality holds for (1) and the Lagrangian dual problem is solvable.

References

  1. Conn, A.R.; Gould, N.I.M.; Toint, P.L. Trust Region Methods; SIAM: Philadelphia, PA, USA, 2000.
  2. Boggs, P.T.; Tolle, J.W. Sequential quadratic programming. Math. Program. Numer. 1995, 4, 1–51. [Google Scholar]
  3. Jeyakumar, V.; Li, G. Strong duality in robust convex programming: Complete characterizations. SIAM J. Optim. 2010, 20, 3384–3407. [Google Scholar] [CrossRef]
  4. Jeyakumar, V.; Li, G. A robust von-Neumann minimax theorem for zero-sum games under bounded payoff uncertainty. Oper. Res. Lett. 2011, 39, 109–114. [Google Scholar] [CrossRef]
  5. Bertsimas, D.; Pachamanova, D.; Sim, M. Robust linear optimization under general norms. Oper. Res. Lett. 2004, 32, 510–516. [Google Scholar] [CrossRef]
  6. Adachi, S.; Iwata, S.; Nakatsukasa, Y.; Takeda, A. Solving the trust region subproblem by a generalized eigenvalue problem. SIAM J. Optim. 2017, 27, 269–291. [Google Scholar] [CrossRef]
  7. Fortin, C.; Wolkowicz, H. The trust region subproblem and semidefinite programming. Optim. Methods Softw. 2004, 19, 41–67. [Google Scholar] [CrossRef]
  8. Gould, N.I.; Lucidi, S.; Roma, M.; Toint, P.L. Solving the trust-region subproblem using the Lanczos method. SIAM J. Optim. 1999, 9, 504–525. [Google Scholar] [CrossRef]
  9. Moré, J.J.; Sorensen, D.C. Computing a trust region step. SIAM J. Sci. Stat. Comput. 1983, 4, 553–572. [Google Scholar] [CrossRef] [Green Version]
  10. Rojas, M.; Santos, S.A.; Sorensen, D.C. A new matrix-free algorithm for the large-scale trust-region subproblem. SIAM J. Optim. 2001, 11, 611–646. [Google Scholar] [CrossRef]
  11. Ye, Y.; Zhang, S. New results on quadratic minimization. SIAM J. Optim. 2003, 14, 245–267. [Google Scholar] [CrossRef]
  12. Beck, A.; Eldar, Y. Strong duality in nonconvex quadratic optimization with two quadratic constraints. SIAM J. Optim. 2006, 17, 844–860. [Google Scholar] [CrossRef] [Green Version]
  13. Sturm, J.F.; Zhang, S. On cones of nonnegative quadratic functions. Math. Oper. Res. 2003, 28, 246–267. [Google Scholar] [CrossRef] [Green Version]
  14. Salahi, M.; Taati, A.; Wolkowicz, H. Local nonglobal minima for solving large-scale extended trust-region subproblems. Comput. Optim. Appl. 2017, 66, 223–244. [Google Scholar] [CrossRef] [Green Version]
  15. Jeyakumar, V.; Li, G. Trust-region problems with linear inequality constraints: Exact SDP relaxation, global optimality and robust optimization. Math. Program. 2014, 147, 171–206. [Google Scholar] [CrossRef]
  16. Hsia, Y.; Sheu, R. Trust region subproblem with a fixed number of additional linear inequality constraints has polynomial complexity. arXiv 2013, arXiv:1312.1398. [Google Scholar]
  17. Ho-Nguyen, N.; Kilinç-Karzan, F. A Second-order cone based approach for solving the trust-region subproblem and its variants. SIAM J. Optim. 2017, 27, 1485–1512. [Google Scholar] [CrossRef] [Green Version]
  18. Jiang, R.; Li, D.; Wu, B. SOCP reformulation for the generalized trust region subproblem via a canonical form of two symmetric matrices. Math. Program. 2018, 169, 531–563. [Google Scholar] [CrossRef]
  19. Ben-Tal, A.; Teboulle, M. Hidden convexity in some nonconvex quadratically constrained quadratic programming. Math. Program. 1996, 72, 51–63. [Google Scholar] [CrossRef]
  20. Feng, J.M.; Lin, G.X.; Sheu, R.L.; Xia, Y. Duality and solutions for quadratic programming over single non-homogeneous quadratic constraint. J. Glob. Optim. 2012, 54, 275–293. [Google Scholar] [CrossRef]
  21. Moré, J.J. Generalizations of the trust region problem. Optim. Methods Softw. 1993, 2, 189–209. [Google Scholar] [CrossRef]
  22. Pong, T.K.; Wolkowicz, H. The generalized trust region subproblem. Comput. Optim. Appl. 2014, 58, 273–322. [Google Scholar] [CrossRef] [Green Version]
  23. Adachi, S.; Nakatsukasa, Y. Eigenvalue-based algorithm and analysis for nonconvex QCQP with one constraint. Math. Program. 2019, 173, 79–116. [Google Scholar] [CrossRef] [Green Version]
  24. Taati, A.; Salahi, M. A conjugate gradient-based algorithm for large-scale quadratic programming problem with one quadratic constraint. Comput. Optim. Appl. 2019, 74, 195–223. [Google Scholar] [CrossRef] [Green Version]
  25. Salahi, M.; Taati, A. An efficient algorithm for solving the generalized trust region subproblem. Comput. Appl. Math. 2018, 37, 395–413. [Google Scholar]
  26. Ben-Tal, A.; Hertog, D.D. Hidden conic quadratic representation of some nonconvex quadratic optimization problems. Math. Program. 2014, 143, 1–29. [Google Scholar] [CrossRef]
  27. Locatelli, M. Some results for quadratic problems with one or two quadratic constraints. Oper. Res. Lett. 2015, 43, 126–131. [Google Scholar] [CrossRef]
  28. Locatelli, M. Exactness conditions for an SDP relaxation of the extended trust region problem. Optim. Lett. 2016, 10, 1141–1151. [Google Scholar] [CrossRef]
  29. Jiang, R.; Li, D. Novel reformulations and efficient algorithms for the generalized trust region subproblem. SIAM J. Optim. 2019, 29, 1603–1633. [Google Scholar] [CrossRef] [Green Version]
  30. Grant, M.; Boyd, S. CVX: Matlab Software for Disciplined Convex Programming, Version 2.1. Available online: http://cvxr.com/cvx (accessed on 15 June 2020).
  31. Lancaster, P.; Rodman, L. Canonical forms for hermitian matrix pairs under strict equivalence and congruence. SIAM Rev. 2005, 47, 407–443. [Google Scholar] [CrossRef]
  32. Ben-Tal, A.; Nemirovski, A. Lectures on Modern Convex Optimization: Analysis, Algorithms, and Engineering Applications; SIAM: Philadelphia, PA, USA, 2001.
  33. Hsia, Y.; Lin, G.X.; Sheu, R.L. A revisit to quadratic programming with one inequality quadratic constraint via matrix pencil. Pac. J. Optim. 2014, 10, 461–481. [Google Scholar]

Share and Cite

MDPI and ACS Style

Almaadeed, T.A.; Taati, A.; Salahi, M.; Hamdi, A. The Generalized Trust-Region Sub-Problem with Additional Linear Inequality Constraints—Two Convex Quadratic Relaxations and Strong Duality. Symmetry 2020, 12, 1369. https://doi.org/10.3390/sym12081369

AMA Style

Almaadeed TA, Taati A, Salahi M, Hamdi A. The Generalized Trust-Region Sub-Problem with Additional Linear Inequality Constraints—Two Convex Quadratic Relaxations and Strong Duality. Symmetry. 2020; 12(8):1369. https://doi.org/10.3390/sym12081369

Chicago/Turabian Style

Almaadeed, Temadher A., Akram Taati, Maziar Salahi, and Abdelouahed Hamdi. 2020. "The Generalized Trust-Region Sub-Problem with Additional Linear Inequality Constraints—Two Convex Quadratic Relaxations and Strong Duality" Symmetry 12, no. 8: 1369. https://doi.org/10.3390/sym12081369

APA Style

Almaadeed, T. A., Taati, A., Salahi, M., & Hamdi, A. (2020). The Generalized Trust-Region Sub-Problem with Additional Linear Inequality Constraints—Two Convex Quadratic Relaxations and Strong Duality. Symmetry, 12(8), 1369. https://doi.org/10.3390/sym12081369

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop