---
title: Quadratic Fractional Programming
url: https://www.emergentmind.com/topics/quadratic-fractional-programming
type: topic
---

# Quadratic Fractional Programming

Quadratic fractional programming (QFP) studies optimization problems whose objective or constraint functions are ratios of quadratic forms. These arise extensively in signal processing, communications, control, and combinatorial optimization, due to natural occurrence of figures of merit like signal-to-noise/interference ratios, Cramér-Rao bounds, and eigenvalue or margin-based metrics in machine learning. The nonconvex, nonhomogeneous rational structure fundamentally distinguishes QFP from classical quadratic programming, necessitating a specialized toolkit that spans parametric analysis, convexification, conic duality, minorization–maximization, semidefinite relaxations, and combinatorial enumeration.

## 1. Canonical Problem Formulations

The prototypical scalar QFP takes the form
\[
\min_{x\in X} \frac{f(x)}{g(x)} \qquad \text{with } f(x) = x^\top A x + 2a^\top x + \alpha, \;\; g(x) = x^\top B x + 2b^\top x + \beta,
\]
subject to $g(x) > 0$ for all $x \in X \subseteq \mathbb{R}^n$ and possibly additional constraints (polyhedral, quadratic, combinatorial, or constant-modulus) [1402.4198]. Multi-objective or vector QFPs generalize this to simultaneous minimization of $m$ coordinate-wise ratios $f_i(x)/g_i(x)$ over a common feasible region [1309.5508].

Problem classes of central importance include:
- Scalar single-ratio minimization/maximization
- Max–min or sum-of-ratios objectives
- QFP with multi-dimensional matrix terms (e.g., $\operatorname{Tr}\left(A(x)B(x)^{-1}A(x)^H\right)$)
- QFP under quadratic, polyhedral, or discrete constraints ($x \in \{-1,1\}^n$, $|x_k|=1$)
- Fractionally constrained fractional programming (FCFP), where both objective and constraints are ratio structured [2503.09922]

## 2. Classical and Semidefinite Programming Methodologies

For convex feasible sets, the standard approach replaces the ratio objective by a parametric family:
\[
F(\lambda) = \inf_{x \in X} \left\{ f(x) - \lambda g(x) \right\},
\]
with the optimal value $\lambda^* = \min_{x \in X} f(x)/g(x)$ characterized as the unique solution to $F(\lambda^*) = 0$ (the Dinkelbach root) [1402.4198, 1211.4664]. For single-ratio with quadratic numerator and denominator, this links QFP to the generalized trust region subproblem and invites strong duality and convex relaxation techniques.

Key results include:
- For two-sided quadratic constraints, $\ell \leq h(x) \leq u$, the exact solution for $\lambda^*$ is computable by semidefinite programming (SDP), using a suitably generalized S-lemma:
  \[
  \lambda^* = \max_{\lambda, \mu \in \mathbb{R}} \left\{ \lambda \,:\, 
  \begin{pmatrix}
    A - \lambda B + \mu C & a - \lambda b + \mu c \\
    (a - \lambda b + \mu c)^\top & \alpha - \lambda \beta + \mu \bar{h}
  \end{pmatrix} 
  \succeq 0
  \right\}
  \]
  with $\bar h = \ell$ or $u$ [1402.4198].

- In the case of a single quadratic constraint $h(x) \leq u$, the standard S-lemma suffices and the SDP relaxation remains exact and polynomial-time [1402.4198].

- Absence of positive-definite pencils (i.e., “hard cases” where no Slater point exists for some $A-\lambda B+pC \succ 0$) can be overcome by exploiting this extended S-lemma and associated constraint qualification (“Assumption B”) [1402.4198].

- Nonconvex quadratic-fractional programs (e.g., with nonconvex sets or concave denominator) can be globally addressed by canonical duality. The original problem is parametrized, and the dual reduces to a low-dimensional (often 2D) concave maximization, enjoying no duality gap under mild conditions [1211.4664].

- For discrete QFP ($x\in\{-1,1\}^n$), e.g.,
  \[
  \min_{x\in\{-1,1\}^n} \frac{x^\top A x + \alpha}{x^\top B x + \beta},
  \]
  with fixed spectral ranks, the global optimum can be found in $O(n^{r_a+r_b+1}\log^2 n)$ via an accelerated Newton–Dinkelbach scheme, provided the diagonal patterns meet appropriate sign conditions [2411.09190].

## 3. Quadratic Transform and Fractional Programming Algorithms

The "quadratic transform" (QT) offers a general-purpose, highly modular approach for scalar and matrix QFP with convex (or nonconvex) constraint sets [2503.09977, 2312.05726, 2301.11269]. The basic device is to decouple each ratio into a jointly biconcave surrogate:
\[
\frac{A(x)}{B(x)} \;\Longleftrightarrow\; \max_{u} \big[ 2u\sqrt{A(x)} - u^2 B(x) \big].
\]
Alternating optimization over auxiliary $u$ (closed form) and $x$ (convex subproblem when $A$ is concave and $B$ convex) yields monotonic ascent and eventual convergence to a stationary point [2503.09977, 2312.05726]. This framework extends to maximization and minimization of sums of (functions of) ratios, matrix ratios (by trace–Schur surrogates), and to mixed max–min formulations [2305.02704].

For multi-term and matrix-valued ratios, the quadratic transform yields subproblems amenable to standard convex or discrete programming techniques, or even matching and assignment solvers (for graph clustering and SVM margin maximization) [2503.09977].

Accelerated variants leveraging Nesterov-type extrapolation or Steffensen's fixed-point schemes can achieve $O(1/k^2)$ local convergence rates [2312.05726, 2601.00564]. Theoretical guarantees and practical performance are established for diverse applications (massive MIMO, radar, ISAC, AoI minimization, SVMs, graph cuts).

## 4. Convexification, Copositive and Combinatorial Approaches

For QFP over polytopes or discrete sets, convexification via homogenization and copositive programming forms a central paradigm [2310.08424]. The convex hull of the feasible set can be characterized in lifted (moment or completely positive) spaces:
\[
\min_{\rho, y, Y} \, b_0 \rho + b^\top y + \operatorname{Tr}(B Y) \quad \text{s.t.} \quad \langle A, Y\rangle + a^\top y + a_0 \rho = 1, \; Z = \begin{pmatrix} \rho & y^\top \\ y & Y \end{pmatrix} \in \mathcal{CP}_{n+1}
\]
where $\mathcal{CP}_{n+1}$ is the cone of completely positive $(n+1)\times(n+1)$ matrices. This reformulation is exact under mild conditions and admits strong relaxations via the RLT or moment hierarchies.

For 0-1 QFP over cardinality-constrained polytopes, hierarchy-based relaxations close optimality gaps several times faster than McCormick or CEF relaxations and are computationally efficient up to $n \sim 800$ [2310.08424]. In low-rank convex–convex cases, global optimality can be ensured by explicit region checking on sign patterns of rank vectors, yielding $O(2^r)$ worst-case complexity, but far better practical scaling [2301.11269].

## 5. Pareto Optimality, Vector QFP, and Multi-objective Extensions

For vector quadratic fractional programming (VQFP), necessary and sufficient Pareto optimality conditions generalize KKT by exploiting the structure of the ratio gradients. For $f_i(x)/g_i(x)$, the stationarity reads
\[
\sum_{i=1}^m t_i \left[\nabla f_i(x^*) - \frac{f_i(x^*)}{g_i(x^*)}\nabla g_i(x^*)\right] + \sum_{j=1}^{\ell} \mu_j \nabla h_j(x^*) = 0,\quad t_i>0,\; \mu_j \geq 0
\]
with appropriate feasibility and complementarity [1309.5508, 1309.5515].

Strong duality holds under Guignard–generalized constraint qualification (GGCQ). Second-order sufficient optimality reduces to positive definiteness of a weighted sum of Hessians on the critical cone. Additional geometric criteria (“radius of efficiency”) certify global vs. local Pareto solutions: efficient computational tests are available for algorithmic termination and local–global optimality distinction [1309.5515].

## 6. Specialized Domains: Fractionally Constrained Problems and Applications

In advanced scenarios such as FCFP—fractionally constrained fractional programming—when both the objectives and constraints are ratios (e.g., RIS-aided ISAC), the quadratic transform can be recursively applied to both constraint and objective ratios [2503.09922]. For constant-modulus or combinatorial constraints, penalty-based or dual variable-based linearization delivers tractable convex or linear subproblems, preserving polynomial solvability in problem dimension. This framework is essential for modern joint sensing–communications optimization, robust MIMO, and statistical waveform design.

The quadratic-fractional paradigm underlies a wide spectrum of engineering and data science problems:
- Power/rate allocation in wireless networks (SINR-based QFP) [2503.09977]
- Robust beamforming under imperfect CSI [2001.00655]
- Radar waveform optimization (CRLB/KLD minimization) [2503.09977, 2601.00564]
- Graph partitioning/object clustering (discrete matrix QFP) [2503.09977]
- SVM margin maximization and multi-class classifier design [2503.09977]
- Joint sensing-communication (RIS-based, ISAC) [2503.09922]
- Age-of-information minimization and queueing systems [2305.02704]

## 7. Algorithmic and Theoretical Summary

Quadratic fractional programming unifies a range of nonconvex, nonhomogeneous optimization tasks for which global stationarity can be certified via majorization-minimization, parametric analysis, or semidefinite duality. The critical techniques include:
- Parametric Dinkelbach-type root-finding for single-ratio cases
- Semidefinite relaxation and S-lemma-based SDP for quadratically constrained QFP, even in “hard cases” without Slater points [1402.4198]
- Canonical duality transforming nonconvex QFP into low-dimensional concave dual maximization [1211.4664]
- Quadratic transform for sum-of-ratios or matrix-ratio FP, yielding efficient AO/MM style alternations with convergence guarantees [2503.09977, 2312.05726, 2305.02704]
- Hierarchical convexification and combinatorial enumeration in discrete/low-rank settings [2310.08424]
- Pareto and multi-objective optimality criteria enabling certified local and global efficiency [1309.5508, 1309.5515]

A directed selection of approaches is provided in the following table:

| Formulation Type             | Representative Method          | Key Reference            |
|------------------------------|-------------------------------|--------------------------|
| Single-ratio, quadratic      | Parametric, S-lemma + SDP     | [1402.4198]              |
| Multi-term (sum-of-ratios)   | Quadratic Transform           | [2503.09977, 2312.05726] |
| Matrix-ratio QFP             | Matrix Quadratic Transform    | [2503.09977, 2601.00564] |
| Mixed max/min (multi-task)   | Unified QT + MM               | [2305.02704]             |
| Discrete/low-rank/fractional | Copositive/convexification    | [2310.08424, 2301.11269] |
| Pareto multi-objective       | KKT/critical cone analysis    | [1309.5508, 1309.5515]   |
| FCFP (objective + constr. FP)| Recursive QT + penalty/LP     | [2503.09922]             |

Despite the diversity of formulations, a recurring backbone is the reduction of ratio structures to biconcave or convexified surrogates amenable to efficient numerical solution, and the availability of theoretically sound optimality and duality guarantees for verifying stationarity or globality.

Source: https://www.emergentmind.com/topics/quadratic-fractional-programming