We propose in this paper a general D.C. decomposition scheme for constructing SDP relaxation formulations for a class of nonconvexquadratic programs with a nonconvexquadratic objective function and convex quadratic ...
详细信息
We propose in this paper a general D.C. decomposition scheme for constructing SDP relaxation formulations for a class of nonconvexquadratic programs with a nonconvexquadratic objective function and convex quadratic constraints. More specifically, we use rank-one matrices and constraint matrices to decompose the indefinite quadratic objective into a D.C. form and underestimate the concave terms in the D.C. decomposition formulation in order to get a convex relaxation of the original problem. We show that the best D.C. decomposition can be identified by solving an SDP problem. By suitably choosing the rank-one matrices and the linear underestimation, we are able to construct convex relaxations that dominate Shor's SDP relaxation and the strengthened SDP relaxation. We then propose an extension of the D.C. decomposition to generate an SDP bound that is tighter than the SDP+RLT bound when additional box constraints are present. We demonstrate via computational results that the optimal D.C. decomposition schemes can generate both tight SDP bounds and feasible solutions with good approximation ratio for nonconvexquadraticallyconstrainedquadratic problems.
In this paper, we present new convex relaxations for nonconvex quadratically constrained quadratic programming (QCQP) problems. While recent research has focused on strengthening convex relaxations of QCQP using the r...
详细信息
In this paper, we present new convex relaxations for nonconvex quadratically constrained quadratic programming (QCQP) problems. While recent research has focused on strengthening convex relaxations of QCQP using the reformulation-linearization technique (RLT), the state-of-the-art methods lose their effectiveness when dealing with (multiple) nonconvexquadratic constraints in QCQP, except for direct lifting and linearization. In this research, we decompose and relax each nonconvex constraint to two second order cone (SOC) constraints and then linearize the products of the SOC constraints and linear constraints to construct some new effective valid constraints. Moreover, we extend the reach of the RLT-like techniques for almost all different types of constraint-pairs (including valid inequalities by linearizing the product of a pair of SOC constraints, and the Hadamard product or the Kronecker product of two respective valid linear matrix inequalities), examine dominance relationships among different valid inequalities, and explore almost all possibilities of gaining benefits from generating valid constraints. We also successfully demonstrate that applying RLT-like techniques to additional redundant linear constraints could reduce the relaxation gap significantly. We demonstrate the efficiency of our results with numerical experiments.
The objective in cancer radiotherapy is to maximize tumor-kill while limiting toxic effects of radiation dose on nearby organs-at-risk (OAR). Given a fixed number of treatment sessions, planners thus face the problem ...
详细信息
The objective in cancer radiotherapy is to maximize tumor-kill while limiting toxic effects of radiation dose on nearby organs-at-risk (OAR). Given a fixed number of treatment sessions, planners thus face the problem of finding a dosing sequence that achieves this goal. This is called the fractionation problem, and has received steady attention over a long history in the clinical literature. Mathematical formulations of the resulting optimization problem utilize the linear-quadratic (LQ) framework to characterize radiation dose-response of tumors and OAR. This yields a nonconvexquadraticallyconstrainedquadratic program. The optimal dosing plan in this forward problem crucially depends on the parameters of the LQ model. Unfortunately, these parameters are difficult to estimate via in vitro or in vivo studies, and as such, their values are unknown to treatment planners. The clinical literature is thus replete with debates about what parameter values will make specific dosing plans effective. This paper formulates this as an inverse optimization problem. The LQ dose-response parameters appear in the objective function, the left hand side, and the right hand side of the forward problem, and none of the existing generic methods can provide an exact solution of the inverse problem. This paper exploits the structure of the problem and identifies all possible parameter values that render the given dosing plan optimal, in closed-form. This closed-form formula is applied to dosing-plans from three clinical studies published within the last two years.
In this paper, a full-duplex (FD) amplify-and-forward (AF) relay is designed to compensate for the duplexing loss of the half-duplex (HD) AF relay. In particular, when there is no direct link between a source and a de...
详细信息
In this paper, a full-duplex (FD) amplify-and-forward (AF) relay is designed to compensate for the duplexing loss of the half-duplex (HD) AF relay. In particular, when there is no direct link between a source and a destination, joint analog domain self-interference suppression and digital domain residual self-interference cancellation is considered with an FD-AF relay having single receive antenna but multiple transmit antennas. Unlike previous approaches, a nonconvex quadratically constrained quadratic programming problem is formulated to find the optimal solution. The end-to-end spectral efficiency or, equivalently, the end-to-end signal-to-interference-plus-noise ratio from the source to the destination is chosen as the objective function to be maximized subject to the average transmit power constraint at the relay. In addition, an average power constraint is imposed on the output of the relay's receive antenna to avoid the nonlinear distortion in the low noise amplifier and the excessive quantization noise in the analog-to-digital converter. Through the systematic reduction and the partitioning of the constraint set, the optimal solution is derived in a closed algorithmic expression and shows how it allocates the transmission power not only in the direction of maximal performance improvement but also in the orthogonal direction in order to balance the system performance and the amount of self interference. It is shown that the optimal FD-AF relay significantly outperforms the optimal HD-AF relay even with the hardware limitations in the RF chain of the relay's receiver being well taken into account.
作者:
Runge, VincentUniv Evry
Univ Paris Saclay CNRS Lab Math & Modelisat Evry F-91037 Evry France
Considering a finite intersection of balls and a finite union of other balls in an Euclidean space, we propose an exact method to test whether the intersection is covered by the union. We reformulate this problem into...
详细信息
Considering a finite intersection of balls and a finite union of other balls in an Euclidean space, we propose an exact method to test whether the intersection is covered by the union. We reformulate this problem into quadraticprogramming problems. For each problem, we study the intersection between a sphere and a Voronoi-like polyhedron. That way, we get information about a possible overlap between the frontier of the union and the intersection of balls. If the polyhedra are non-degenerate, the initial nonconvex geometric problem, which is NP-hard in general, is tractable in polynomial time by convex optimization tools and vertex enumeration. Under some mild conditions, the vertex enumeration can be skipped. Simulations highlight the accuracy and efficiency of our approach compared with competing algorithms in Python for nonconvex quadratically constrained quadratic programming. This work is motivated by an application in statistics to the problem of multidimensional changepoint detection using pruned dynamic programming algorithms.
In this paper, we propose a branch-and-cut algorithm for solving a nonconvex quadratically constrained quadratic programming (QCQP) problem with a nonempty bounded feasible domain. The problem is first transformed int...
详细信息
In this paper, we propose a branch-and-cut algorithm for solving a nonconvex quadratically constrained quadratic programming (QCQP) problem with a nonempty bounded feasible domain. The problem is first transformed into a linear conic programming problem, and then approximated by semidefinite programming problems over different intervals. In order to improve the lower bounds, polar cuts, generated from corresponding cut-generation problems, are embedded in a branch-and-cut algorithm to yield a globally epsilon(r)-epsilon(z)-optimal solution (with respect to feasibility and optimality, respectively) in a finite number of iterations. To enhance the computational speed, an adaptive branch-and-cut rule is adopted. Numerical experiments indicate that the number of explored nodes required for solving QCQP problems can be significantly reduced by employing the proposed polar cuts.
This paper develops new semidefinite programming(SDP)relaxation techniques for two classes of mixed binary quadraticallyconstrainedquadratic programs and analyzes their approximation *** first class of problems fin...
详细信息
This paper develops new semidefinite programming(SDP)relaxation techniques for two classes of mixed binary quadraticallyconstrainedquadratic programs and analyzes their approximation *** first class of problems finds two minimum norm vectors in N-dimensional real or complex Euclidean space,such that M out of 2M concave quadratic constraints are *** employing a special randomized rounding procedure,we show that the ratio between the norm of the optimal solution of this model and its SDP relaxation is upper bounded by 54πM2 in the real case and by 24√Mπin the complex *** second class of problems finds a series of minimum norm vectors subject to a set of quadratic constraints and cardinality constraints with both binary and continuous *** show that in this case the approximation ratio is also bounded and independent of problem dimension for both the real and the complex cases.
Accurately monitoring the system's operating point is central to the reliable and economic operation of an autonomous energy grid. Power system state estimation (PSSE) aims to obtain complete voltage magnitude and...
详细信息
Accurately monitoring the system's operating point is central to the reliable and economic operation of an autonomous energy grid. Power system state estimation (PSSE) aims to obtain complete voltage magnitude and angle information at each bus given a number of system variables at selected buses and lines. Power flow analysis amounts to solving a set of noise-free power flow equations, and is cast as a special case of PSSE. Physical laws dictate quadratic relationships between available quantities and unknown voltages, rendering general instances of power flow and PSSE nonconvex and NP-hard. Past approaches are largely based on gradient-type iterative procedures or semidefinite relaxation (SDR). Due to nonconvexity, the solution obtained via gradienttype schemes depends on initialization, while SDR methods do not perform as desired in challenging scenarios. This paper puts forth novel feasible point pursuit (FPP)-based solvers for power flow analysis and PSSE, which iteratively seek feasible solutions for a nonconvex quadratically constrained quadratic programming reformulation of the weighted least-squares (WLS). Relative to the prior art, the developed solvers offer superior numerical performance at the cost of higher computational complexity. Furthermore, they converge to a stationary point of the WLS problem. As a baseline for comparing different estimators, the Cramer-Rao lower bound is derived for the fundamental PSSE problem in this paper. Judicious numerical tests on several IEEE benchmark systems showcase markedly improved performance of our FPP-based solvers for both power flow and PSSE tasks over popular WLS-based Gauss-Newton iterations and SDR approaches.
暂无评论