Tutorials › Statistical Theory I - Probability and Distribution Theory › Signs, Bounds, and Support of Quadratic Forms

Gaussian Quadratic Forms · Tutorial 754 of 994

Signs, Bounds, and Support of Quadratic Forms

Learn how eigenvalues control the sign, bounds, and support of Gaussian quadratic forms, and why support can include values the random variable never attains.

Advanced 10 min read

What You'll Learn

  • Express a Gaussian quadratic form in orthogonal eigen-coordinates.
  • Identify when the form is nonnegative, nonpositive, or indefinite.
  • Derive sharp eigenvalue bounds for vectors of fixed length.
  • Determine support from the range of a continuous quadratic map and Gaussian full support.
  • Distinguish almost-sure values from points in the distribution’s support.

Why the Eigenvalues Matter

A quadratic form can look complicated in its original coordinates. For a Gaussian vector, however, an orthogonal change of coordinates turns it into a weighted sum of squares. The signs of those weights—the eigenvalues of the matrix—then reveal whether the form can be negative or positive, give exact bounds when the vector’s length is fixed, and identify the support of its distribution.

Throughout, let \(U\sim N_p(\mu,I_p)\), where \(I_p\) is the \(p\times p\) identity matrix, and let \(A\) be a real symmetric \(p\times p\) matrix. We study \(Q=U^\mathsf{T}AU\). The identity covariance ensures that an orthogonal change of coordinates preserves the standard-normal structure; the mean may be nonzero.

Definition: A symmetric matrix \(A\) is positive semidefinite if \(x^\mathsf{T}Ax\geq 0\) for every \(x\in\mathbb{R}^p\), and negative semidefinite if \(x^\mathsf{T}Ax\leq 0\) for every \(x\). It is indefinite if its quadratic form takes both positive and negative values. The support of a real random variable is the set of points such that every open neighborhood of each point has positive probability.

Diagonal Coordinates and Eigenvalue Bounds

By the spectral theorem, there is an orthogonal matrix \(H\) and a diagonal matrix \(\Lambda=\operatorname{diag}(\lambda_1,\ldots,\lambda_p)\) such that \(A=H\Lambda H^\mathsf{T}\). Set \(V=H^\mathsf{T}U\) and \(\nu=H^\mathsf{T}\mu\). Then \(V\sim N_p(\nu,I_p)\), and

$$ Q=V^\mathsf{T}\Lambda V=\sum_{i=1}^p\lambda_iV_i^2. $$

Thus the eigenvalues are the weights on the squared coordinates. Since the squared coordinates are nonnegative, nonnegative eigenvalues give a nonnegative form; nonpositive eigenvalues give a nonpositive form. Eigenvalues of both signs permit both signs of the form.

Theorem: Let \(\lambda_{\min}\) and \(\lambda_{\max}\) denote the smallest and largest eigenvalues of \(A\). For every \(x\in\mathbb{R}^p\),
$$ \lambda_{\min}\lVert x\rVert^2\leq x^\mathsf{T}Ax\leq\lambda_{\max}\lVert x\rVert^2. $$
Both bounds are attained when \(x\) is a corresponding unit eigenvector. In particular, on the sphere \(\lVert x\rVert=r\), the minimum and maximum of \(x^\mathsf{T}Ax\) are \(r^2\lambda_{\min}\) and \(r^2\lambda_{\max}\).

Proof. Write \(x=Hy\), which preserves length because \(H\) is orthogonal. Then \(\lVert x\rVert^2=\sum_i y_i^2\) and \(x^\mathsf{T}Ax=\sum_i\lambda_i y_i^2\). For every \(i\), \(\lambda_{\min}\leq\lambda_i\leq\lambda_{\max}\); multiplying by \(y_i^2\geq 0\) and summing gives the two inequalities. If \(y\) is a unit coordinate vector for an eigenvalue \(\lambda_{\min}\) or \(\lambda_{\max}\), the corresponding bound is an equality. Scaling that eigenvector to length \(r\) proves the sphere statement. \(\square\)

Worked Example: A Positive Semidefinite Form With a Zero Eigenvalue

Let \(A=\begin{pmatrix}2&0\\0&0\end{pmatrix}\), and let \(U\sim N_2((1,-2),I_2)\). The eigenvalues are \(2\) and \(0\), so

$$ Q=U^\mathsf{T}AU=2U_1^2. $$

For every vector \(u=(u_1,u_2)^\mathsf{T}\), \(u^\mathsf{T}Au=2u_1^2\geq 0\). The eigenvalue bounds give \(0\leq u^\mathsf{T}Au\leq 2\lVert u\rVert^2\); on the circle \(\lVert u\rVert=3\), the exact range is \([0,18]\). The lower endpoint occurs at \(u=(0,3)^\mathsf{T}\), and the upper endpoint at \(u=(3,0)^\mathsf{T}\).

For the random vector \(U\), \(Q\geq 0\) always, and \(Q=0\) only if \(U_1=0\), an event of probability zero. Nevertheless, values of \(U_1\) arbitrarily close to zero have positive probability, so values of \(Q\) arbitrarily close to zero are possible with positive probability.

Worked Example: An Indefinite Form in Non-diagonal Coordinates

Take \(A=\begin{pmatrix}1&2\\2&1\end{pmatrix}\). Its eigenvectors are proportional to \((1,1)^\mathsf{T}\) and \((1,-1)^\mathsf{T}\), with eigenvalues \(3\) and \(-1\). For \(u=(u_1,u_2)^\mathsf{T}\), the form is

$$ u^\mathsf{T}Au=u_1^2+4u_1u_2+u_2^2 =3\left(\frac{u_1+u_2}{\sqrt{2}}\right)^2 -\left(\frac{u_1-u_2}{\sqrt{2}}\right)^2. $$

At \(u=(1,1)^\mathsf{T}\), the value is \(6>0\); at \(u=(1,-1)^\mathsf{T}\), it is \(-2<0\). Hence \(A\) is indefinite. On the unit circle, the theorem gives the sharp bounds \(-1\leq u^\mathsf{T}Au\leq 3\). On a circle of radius \(2\), they become \(-4\leq u^\mathsf{T}Au\leq 12\).

Sign and Support for a Gaussian Vector

The matrix determines the sign for every possible input vector. For a Gaussian vector with identity covariance, it also determines which signs occur with positive probability. A Gaussian density with positive-definite covariance is positive everywhere, so every nonempty open ball has positive probability. In particular, if the form is negative at even one vector, continuity makes it negative in a neighborhood of that vector, and the Gaussian vector enters that neighborhood with positive probability.

Theorem: For \(U\sim N_p(\mu,I_p)\) and symmetric \(A\):
  • \(Q=U^\mathsf{T}AU\geq 0\) almost surely if and only if \(A\) is positive semidefinite.
  • \(Q\leq 0\) almost surely if and only if \(A\) is negative semidefinite.
  • If \(A\) is indefinite, then \(P(Q>0)>0\) and \(P(Q<0)>0\), and the support of \(Q\) is all of \(\mathbb{R}\).
  • If \(A\) is nonzero positive semidefinite, the support of \(Q\) is \([0,\infty)\), although \(Q=0\) has probability zero. For nonzero negative semidefinite \(A\), the support is \((-\infty,0]\), and \(Q=0\) has probability zero.
  • If \(A=0\), then \(Q=0\) almost surely and its support is \(\{0\}\).

Proof. First, a positive semidefinite \(A\) makes \(u^\mathsf{T}Au\geq 0\) for every \(u\), so \(Q\geq 0\) almost surely. Conversely, if \(A\) is not positive semidefinite, there is a vector \(x\) with \(x^\mathsf{T}Ax<0\). Continuity gives an open ball around \(x\) on which the form remains negative. Since \(U\) has a density positive everywhere, that ball has positive probability. Thus \(Q\geq 0\) almost surely is impossible. Applying the same argument to \(-A\) proves the negative-semidefinite assertion and shows that an indefinite form has both signs with positive probability.

To identify support, let \(g(x)=x^\mathsf{T}Ax\). Continuity and full support of \(U\) imply that the support of \(g(U)\) is the closure of \(g(\mathbb{R}^p)\): if \(t=g(x)\), every neighborhood of \(t\) contains values of \(g\) on an open neighborhood of \(x\), which has positive Gaussian probability. Conversely, if \(t\) lies outside the closure of the image, some neighborhood of \(t\) cannot be reached by \(g(U)\).

If \(A\) is nonzero positive semidefinite, its image under \(g\) is exactly \([0,\infty)\): zero is attained at \(x=0\), and for a unit eigenvector \(v\) with eigenvalue \(\lambda>0\), \(g(cv)=c^2\lambda\) ranges over all nonnegative values as \(c\) varies. The negative-semidefinite case gives \((-\infty,0]\). If \(A\) is indefinite, choose vectors with positive and negative form values. Scaling those vectors produces arbitrarily large values of both signs, and the continuous path between them takes every intermediate value. Its image is therefore \(\mathbb{R}\). Finally, for nonzero positive semidefinite \(A\), \(g(U)=0\) requires \(U\) to lie in the proper subspace \(\ker(A)\), which has probability zero under a density; the negative-semidefinite case follows by considering \(-A\). If \(A=0\), the form is identically zero. \(\square\)

Worked Example: Positive Definite Does Not Mean Support Excludes Zero

Let \(A=\begin{pmatrix}1&0\\0&4\end{pmatrix}\) and \(U\sim N_2((2,-1),I_2)\). Then \(Q=U_1^2+4U_2^2\). Both eigenvalues are positive, so \(Q>0\) unless \(U=(0,0)^\mathsf{T}\), an event of probability zero. Yet the support is \([0,\infty)\): for every \(t\geq 0\), the vector \(u=(\sqrt{t},0)^\mathsf{T}\) gives \(u^\mathsf{T}Au=t\), and the Gaussian density is positive near that vector. In particular, every interval \([0,\varepsilon)\) with \(\varepsilon>0\) has positive probability.

Worked Example: A Negative Semidefinite Form

Let \(A=\begin{pmatrix}-4&0\\0&-1\end{pmatrix}\) and \(U\sim N_2((0,1),I_2)\). Then \(Q=-4U_1^2-U_2^2\), with eigenvalues \(-4\) and \(-1\). For every \(u\),

$$ -4\lVert u\rVert^2\leq u^\mathsf{T}Au\leq-\lVert u\rVert^2. $$

On the circle \(\lVert u\rVert=3\), the sharp range is \([-36,-9]\). Without a restriction on length, the support of \(Q\) is \((-\infty,0]\). The value zero is in the support because vectors near the origin have positive probability, but it is not attained with positive probability: \(Q=0\) requires \(U=(0,0)^\mathsf{T}\).

Common Pitfall and Practical Use

A common mistake is to say that a positive-definite form has support \((0,\infty)\) because its value is strictly positive almost surely. That reasoning confuses values attained with positive probability and points in the support. A neighborhood of zero still has positive probability, however small it may be; therefore zero belongs to the support. The same distinction applies at the upper or lower endpoints of supports more generally.

The eigenvalue check is useful before attempting a more detailed distributional calculation. It establishes whether negative values are possible, supplies sharp bounds when the vector length is fixed, and reveals whether a proposed support is plausible. It also shows why the mean vector does not change these sign and support conclusions in this setting: a nondegenerate Gaussian vector still has positive density throughout \(\mathbb{R}^p\). The mean changes probabilities, but not the range available to the quadratic map.

Check Your Understanding

Use eigenvalues and the full-support property of the Gaussian vector to answer these questions.

  1. A symmetric matrix has eigenvalues \(0,2,5\). What can be said about the sign of its quadratic form, and what is the support of the corresponding noncentral Gaussian form?
  2. A symmetric matrix has eigenvalues \(-3,1\). What are the sharp bounds for its quadratic form on the circle of radius \(2\)?
  3. For a nonzero positive-definite matrix, can the associated Gaussian quadratic form equal zero with positive probability? Is zero in its support?
  4. Why does an indefinite matrix imply positive probability of both a positive and a negative quadratic-form value?
  5. If all eigenvalues are zero, what are the value and support of the quadratic form?