LA

Linear Algebra

6 theorems

01 — Linear Algebra

Comprehensive study guide for Bocconi Math Module 2 (cod. 30063), General Exam. Source materials: Alice Sicconi's Proofs - Linear Algebra.pdf, Linear Algebra and Differential Calculus.pdf, lect{1,2,3,5}_lin_alg.pdf, TA41_lin_alg.pdf, General_24524_ENG_SOL.pdf, Final exercises solution.pdf.


§1. Overview & Exam Relevance

Linear Algebra is the first block of the first partial and accounts for approximately 27% of the general exam (typically 1–2 MCQs worth 5 pts each plus one open-ended question worth up to 20 pts — about 30–40 pts of the 150-pt exam).

Topic scope. The exam tests:

  • matrix arithmetic: products, determinants, minors, transposition, the trace;
  • eigenvalues, eigenvectors, eigenspaces, characteristic polynomial;
  • diagonalization of symmetric matrices via the Spectral Theorem;
  • quadratic forms f(x)=xTAxf(\mathbf{x}) = \mathbf{x}^T A \mathbf{x} and the bijection with symmetric matrices;
  • classification (sign) of quadratic forms using (i) eigenvalue signs or (ii) Sylvester–Jacobi leading principal minors.

Typical MCQ patterns (from May 2024 General exam).

  • MCQ1 (Mode A) / MCQ6 (Mode B): "Which of the following is an eigenvector of AA?" — plug each candidate into AvA\mathbf{v} and test Av=λvA\mathbf{v} = \lambda \mathbf{v}.
  • Similar MCQs: "The eigenvalues of AA are …" with multiplicity traps — shortcut via λi=trA\sum \lambda_i = \operatorname{tr} A and λi=detA\prod \lambda_i = \det A.

Typical open-ended pattern.

  • Q7 (Mode A) / Q12 (Mode B), May 2024: state the correspondence theorem (Prop 1201), write the quadratic form attached to a symmetric matrix, prove it is positive semi-definite but not positive definite, find a basis of R3\mathbb{R}^3 made of eigenvectors.
  • "Final exercises" Q1: quadratic form f(x,y,z)=kx2+ky24z22xzf(x,y,z) = kx^2 + ky^2 - 4z^2 - 2xz — find kk so that the origin is a strong global maximum, determine eigenvalue signs, find the kernel for k=0k=0.

Why this topic is high-leverage for the final partial.

  • The Spectral Theorem, Sylvester–Jacobi, and quadratic form signing show up again in Differential Calculus (Hessian classification of stationary points — the Hessian is a symmetric matrix and its quadratic form is the second-order Taylor remainder). Mastering linear algebra pays double dividends.

§2. Definitions

2.1 Vector and matrix

A vector xRn\mathbf{x} \in \mathbb{R}^n is a column of nn real numbers: x=(x1x2xn).\mathbf{x} = \begin{pmatrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{pmatrix}.

A matrix AA of order m×nm \times n is a rectangular array of real entries aija_{ij} with 1im1 \le i \le m, 1jn1 \le j \le n. A square matrix of order nn has m=nm = n.

Dot product of two vectors: xy=i=1nxiyi=xTy\mathbf{x} \cdot \mathbf{y} = \sum_{i=1}^n x_i y_i = \mathbf{x}^T \mathbf{y}. Orthogonal vectors: xy=0\mathbf{x} \cdot \mathbf{y} = 0. Orthonormal vectors: orthogonal AND x=y=1\|\mathbf{x}\| = \|\mathbf{y}\| = 1. Norm: x=xx\|\mathbf{x}\| = \sqrt{\mathbf{x} \cdot \mathbf{x}}.

2.2 Transpose, symmetric, orthogonal

Transpose ATA^T: swap rows and columns — (AT)ij=aji(A^T)_{ij} = a_{ji}.

Symmetric matrix: AT=AA^T = A, equivalently aij=ajia_{ij} = a_{ji} for all i,ji,j. Symmetric matrices are always square.

Orthogonal matrix BB (of order nn): BBT=BTB=IB \cdot B^T = B^T \cdot B = I. Equivalently B1=BTB^{-1} = B^T. Key properties. An orthogonal matrix is invertible, its rows and columns each form an orthonormal set, and detB=±1\det B = \pm 1. An orthogonal matrix is not necessarily symmetric.

2.3 Determinant and minors

The determinant det(A)\det(A) is a scalar function on square matrices.

  • 2×22 \times 2: det(abcd)=adbc\det \begin{pmatrix} a & b \\ c & d \end{pmatrix} = ad - bc.
  • 3×33 \times 3: Laplace expansion along any row/column using the sign pattern (+++++)\begin{pmatrix} + & - & + \\ - & + & - \\ + & - & + \end{pmatrix}.
  • Rules: det(AT)=det(A)\det(A^T) = \det(A), det(A1)=1/det(A)\det(A^{-1}) = 1/\det(A), det(AB)=det(A)det(B)\det(AB) = \det(A) \det(B).
  • Singular / non-singular: AA is singular iff detA=0\det A = 0 (rows/columns are linearly dependent); AA is non-singular iff detA0\det A \neq 0.

Minor: determinant of a submatrix obtained by removing any number of rows and columns. Principal minor: minor obtained by removing the same rows as columns (e.g., remove row 2 and column 2). Leading (North-West) principal minor detAk\det A_k: the k×kk \times k minor obtained by keeping the first kk rows AND first kk columns (i.e., removing the last nkn-k rows and columns). For a 3×33 \times 3 matrix there are three leading principal minors: detA1=a11\det A_1 = a_{11}, detA2=det(a11a12a21a22)\det A_2 = \det \begin{pmatrix} a_{11} & a_{12} \\ a_{21} & a_{22} \end{pmatrix}, detA3=detA\det A_3 = \det A.

2.4 Rank and linear independence

A set of vectors {v1,,vk}\{\mathbf{v}_1, \ldots, \mathbf{v}_k\} is linearly independent iff the only solution to α1v1++αkvk=0\alpha_1 \mathbf{v}_1 + \cdots + \alpha_k \mathbf{v}_k = \mathbf{0} is α1==αk=0\alpha_1 = \cdots = \alpha_k = 0. Otherwise linearly dependent.

Rank of a matrix = number of linearly independent rows = number of linearly independent columns. For a square matrix of order nn, rank(A)=n    detA0\operatorname{rank}(A) = n \iff \det A \neq 0.

2.5 Basis, dimension, subspace

A vector subspace VRnV \subseteq \mathbb{R}^n is a non-empty subset closed under linear combinations: x,yV\mathbf{x}, \mathbf{y} \in V and α,βR\alpha, \beta \in \mathbb{R} imply αx+βyV\alpha \mathbf{x} + \beta \mathbf{y} \in V (in particular 0V\mathbf{0} \in V).

A spanning set of VV is a set of vectors whose linear combinations generate all of VV. A basis is a spanning set that is also linearly independent. Dimension dimV\dim V is the number of vectors in any basis.

In R2\mathbb{R}^2 there are 3 kinds of subspaces: {0}\{\mathbf{0}\} (dim 0), lines through 0\mathbf{0} (dim 1), and all of R2\mathbb{R}^2 (dim 2). In R3\mathbb{R}^3: {0}\{\mathbf{0}\}, lines through 0\mathbf{0}, planes through 0\mathbf{0}, and all of R3\mathbb{R}^3.

2.6 Eigenvalue, eigenvector, eigenpair

Let AA be a symmetric matrix of order nn. A scalar λR\lambda \in \mathbb{R} is an eigenvalue of AA and a non-zero vector vRn\mathbf{v} \in \mathbb{R}^n is an associated eigenvector if Av=λv,v0.\boxed{A \mathbf{v} = \lambda \mathbf{v}, \qquad \mathbf{v} \neq \mathbf{0}.}

The pair (λ,v)(\lambda, \mathbf{v}) is an eigenpair. The spectrum of AA is σ(A)={λR:det(AλI)=0}\sigma(A) = \{\lambda \in \mathbb{R} : \det(A - \lambda I) = 0\}.

Remark. The restriction to symmetric matrices is Sicconi's convention and matches the Bocconi syllabus — all eigenvalue theorems in §3 assume AA symmetric, which guarantees every root of the characteristic polynomial is real.

2.7 Eigenspace

The eigenspace associated with eigenvalue λ\lambda is Wλ={xRn:(AλI)x=0}.W_\lambda = \{ \mathbf{x} \in \mathbb{R}^n : (A - \lambda I)\mathbf{x} = \mathbf{0} \}. It contains every eigenvector associated with λ\lambda, plus the zero vector. By Theorem 1181 (§3.2), WλW_\lambda is a vector subspace of Rn\mathbb{R}^n.

Dimension / multiplicity: dimWλ\dim W_\lambda = algebraic multiplicity of λ\lambda as a root of pA(λ)p_A(\lambda) = nrank(AλI)n - \operatorname{rank}(A - \lambda I). For a symmetric matrix the algebraic and geometric multiplicities coincide.

2.8 Characteristic polynomial

The characteristic polynomial is pA(λ)=det(AλI)p_A(\lambda) = \det(A - \lambda I), a polynomial of degree nn. Its roots are exactly the eigenvalues of AA (Theorem 1179).

Shortcut for 2×22\times 2: pA(λ)=λ2(trA)λ+detAp_A(\lambda) = \lambda^2 - (\operatorname{tr} A)\lambda + \det A, so λ1+λ2=trA\lambda_1 + \lambda_2 = \operatorname{tr} A and λ1λ2=detA\lambda_1 \lambda_2 = \det A. For general nn: iλi=trA\sum_i \lambda_i = \operatorname{tr} A, iλi=detA\prod_i \lambda_i = \det A (each counted with multiplicity).

2.9 Eigenbasis

An eigenbasis of Rn\mathbb{R}^n is a basis made of orthonormal eigenvectors of a symmetric matrix AA. Any pair of eigenvectors from different eigenspaces is automatically orthogonal (Theorem 1184); eigenvectors from the same eigenspace may need Gram–Schmidt to be orthogonalised and then normalised.

2.10 Diagonalization

A matrix AA is orthogonally diagonalizable if there exists an orthogonal matrix BB (i.e., BT=B1B^T = B^{-1}) such that BTAB=Λ,B^T A B = \Lambda, where Λ=diag(λ1,,λn)\Lambda = \operatorname{diag}(\lambda_1, \ldots, \lambda_n) is the diagonal matrix of eigenvalues. Equivalently, A=BΛBTA = B \Lambda B^T (the spectral decomposition). The Spectral Theorem (§3.7) says every symmetric real matrix is orthogonally diagonalizable.

2.11 Quadratic form

A form f:RnRf : \mathbb{R}^n \to \mathbb{R} is a sum of monomials of the same degree. A quadratic form is a form of degree 2: a sum of monomials each of degree 2, e.g., f(x,y,z)=3x25xy+7yz+z2f(x,y,z) = 3x^2 - 5xy + 7yz + z^2.

Key identity (Prop 1201, §3.5): every quadratic form can be written as f(x)=xTAxf(\mathbf{x}) = \mathbf{x}^T A \mathbf{x} for a unique symmetric AA. The entries are read off as:

  • aiia_{ii} = coefficient of xi2x_i^2;
  • aij=ajia_{ij} = a_{ji} = 12\tfrac12 \cdot (coefficient of xixjx_i x_j).

Note. f(0)=0f(\mathbf{0}) = 0 always.

2.12 Positive/negative (semi)definite — indefinite

Let AA be a symmetric matrix of order nn and f(x)=xTAxf(\mathbf{x}) = \mathbf{x}^T A \mathbf{x}. Then AA (equivalently, ff) is:

| Class | Condition on ff | Geometric picture at 0\mathbf{0} | |-------|--------------------|--------------------------| | positive definite | f(x)>0  x0f(\mathbf{x}) > 0\;\forall \mathbf{x}\neq\mathbf{0} | strict (unique) global minimum | | positive semi-definite (not def.) | f(x)0  xf(\mathbf{x}) \ge 0\;\forall\mathbf{x}, x0\exists \mathbf{x}\neq\mathbf{0} with f(x)=0f(\mathbf{x})=0 | weak global minimum (infinitely many minimizers along a line) | | negative definite | f(x)<0  x0f(\mathbf{x}) < 0\;\forall \mathbf{x}\neq\mathbf{0} | strict global maximum | | negative semi-definite (not def.) | f(x)0  xf(\mathbf{x}) \le 0\;\forall\mathbf{x}, x0\exists \mathbf{x}\neq\mathbf{0} with f(x)=0f(\mathbf{x})=0 | weak global maximum | | indefinite | takes both signs | saddle at 0\mathbf{0} |

Pictures (Sicconi's lecture notes). Bowl f=x2+y2f=x^2+y^2 = positive definite. Half-pipe f=x2f=x^2 = positive semi-definite (not definite). Pringle f=x2y2f=x^2-y^2 = indefinite.


§3. Theorems, Propositions & Proofs

3.1 Theorem 1179 — Characterization of eigenvalues as roots of the characteristic equation

Statement. Let AA be a symmetric matrix of order nn. A real scalar λR\lambda \in \mathbb{R} is an eigenvalue of AA if and only if det(AλI)=0\det(A - \lambda I) = 0.

Source: Proofs - Linear Algebra.pdf p.3 (handwritten theorem 1179).

Proof.

(\Rightarrow) λ\lambda is an eigenvalue det(AλI)=0\Longrightarrow \det(A - \lambda I) = 0.

By definition there exists x0Rn\mathbf{x} \neq \mathbf{0} \in \mathbb{R}^n such that Ax=λx.A \mathbf{x} = \lambda \mathbf{x}. This is the definition of eigenvalue/eigenvector. Subtracting λx\lambda \mathbf{x} from both sides: Axλx=0,AxλIx=0,(AλI)x=0.A \mathbf{x} - \lambda \mathbf{x} = \mathbf{0}, \qquad A \mathbf{x} - \lambda I \mathbf{x} = \mathbf{0}, \qquad (A - \lambda I)\mathbf{x} = \mathbf{0}. Set B:=AλIB := A - \lambda I. Then Bx=0B \mathbf{x} = \mathbf{0} is a homogeneous linear system. Since x=0\mathbf{x} = \mathbf{0} is a trivial solution AND the eigenvector x0\mathbf{x}\neq\mathbf{0} is also a solution, the homogeneous system has infinitely many solutions. This forces detB=0\det B = 0, i.e., det(AλI)=0.\det(A - \lambda I) = 0.

(\Leftarrow) det(AλI)=0λ\det(A - \lambda I) = 0 \Longrightarrow \lambda is an eigenvalue.

The homogeneous system (AλI)x=0(A - \lambda I)\mathbf{x} = \mathbf{0} with det(AλI)=0\det(A - \lambda I) = 0 has infinitely many solutions. Pick one non-trivial solution x0Rn\mathbf{x} \neq \mathbf{0} \in \mathbb{R}^n with (AλI)x=0(A - \lambda I)\mathbf{x} = \mathbf{0}. Then AxλIx=0    Axλx=0    Ax=λx.A \mathbf{x} - \lambda I \mathbf{x} = \mathbf{0} \;\Longrightarrow\; A \mathbf{x} - \lambda \mathbf{x} = \mathbf{0} \;\Longrightarrow\; A \mathbf{x} = \lambda \mathbf{x}. This is the definition of an eigenpair. Hence λ\lambda is an eigenvalue of AA. \blacksquare


3.2 Theorem 1181 — Eigenspace is a subspace of Rn\mathbb{R}^n

Statement. Let AA be a symmetric matrix of order nn and λR\lambda \in \mathbb{R} an eigenvalue of AA. Then the eigenspace Wλ={xRn:(AλI)x=0}W_\lambda = \{\mathbf{x} \in \mathbb{R}^n : (A - \lambda I)\mathbf{x} = \mathbf{0}\} is a vector subspace of Rn\mathbb{R}^n.

Source: Proofs - Linear Algebra.pdf p.4 (handwritten theorem 1181).

Proof.

To show WλW_\lambda is a subspace, we verify it is closed under linear combinations. Take x1,x2Wλ\mathbf{x}_1, \mathbf{x}_2 \in W_\lambda and scalars α,βR\alpha, \beta \in \mathbb{R}. We compute: (AλI)(αx1+βx2)  =  (AλI)αx1+(AλI)βx2  =  α(AλI)x1=0+β(AλI)x2=0  =  α0+β0  =  0.(A - \lambda I)(\alpha \mathbf{x}_1 + \beta \mathbf{x}_2) \;=\; (A - \lambda I)\alpha \mathbf{x}_1 + (A - \lambda I)\beta \mathbf{x}_2 \;=\; \alpha\underbrace{(A - \lambda I)\mathbf{x}_1}_{=\mathbf{0}} + \beta\underbrace{(A - \lambda I)\mathbf{x}_2}_{=\mathbf{0}} \;=\; \alpha \mathbf{0} + \beta \mathbf{0} \;=\; \mathbf{0}. Linear combination; then apply the definition of eigenpair (AλI)xi=0(A - \lambda I)\mathbf{x}_i = \mathbf{0} to each term. Therefore αx1+βx2Wλ\alpha \mathbf{x}_1 + \beta \mathbf{x}_2 \in W_\lambda. In particular choosing α=β=0\alpha = \beta = 0 shows 0Wλ\mathbf{0} \in W_\lambda. Hence WλW_\lambda is a vector subspace of Rn\mathbb{R}^n. \blacksquare


3.3 Theorem 1184 — Eigenvectors of distinct eigenvalues are orthogonal

Statement. Let AA be a symmetric matrix of order nn and (λ1,x1)(\lambda_1, \mathbf{x}_1), (λ2,x2)(\lambda_2, \mathbf{x}_2) be eigenpairs with λ1λ2\lambda_1 \neq \lambda_2. Then x1x2\mathbf{x}_1 \perp \mathbf{x}_2, i.e., x1x2=0\mathbf{x}_1 \cdot \mathbf{x}_2 = 0.

Source: Proofs - Linear Algebra.pdf p.5 (handwritten theorem 1184).

Proof.

Start from the two eigen-equations Ax1=λ1x1A \mathbf{x}_1 = \lambda_1 \mathbf{x}_1 and Ax2=λ2x2A \mathbf{x}_2 = \lambda_2 \mathbf{x}_2. Taking the scalar x1T(Ax2)\mathbf{x}_1^T (A \mathbf{x}_2): \mathbf{x}_1^T (A \mathbf{x}_2) = \mathbf{x}_1^T (\lambda_2 \mathbf{x}_2) = \lambda_2 (\mathbf{x}_1 \cdot \mathbf{x}_2). \tag{I} And similarly: \mathbf{x}_2^T (A \mathbf{x}_1) = \mathbf{x}_2^T (\lambda_1 \mathbf{x}_1) = \lambda_1 (\mathbf{x}_1 \cdot \mathbf{x}_2). \tag{II}

Now the key manipulation. Since x1T(Ax2)\mathbf{x}_1^T (A \mathbf{x}_2) is a 1×11 \times 1 scalar α\alpha, we have α=αT\alpha = \alpha^T. Transposing: x1T(Ax2)=(x1T(Ax2))T=(Ax2)Tx1=x2TATx1.\mathbf{x}_1^T (A \mathbf{x}_2) = \bigl(\mathbf{x}_1^T (A \mathbf{x}_2)\bigr)^T = (A \mathbf{x}_2)^T \mathbf{x}_1 = \mathbf{x}_2^T A^T \mathbf{x}_1. When transposing a product, the order switches: (PQ)T=QTPT(PQ)^T = Q^T P^T. Because AA is symmetric, AT=AA^T = A, and so \mathbf{x}_1^T (A \mathbf{x}_2) = \mathbf{x}_2^T A \mathbf{x}_1. \tag{III}

Combining (I), (II), (III): λ2(x1x2)  =  x1TAx2  =  x2TAx1  =  λ1(x1x2).\lambda_2 (\mathbf{x}_1 \cdot \mathbf{x}_2) \;=\; \mathbf{x}_1^T A \mathbf{x}_2 \;=\; \mathbf{x}_2^T A \mathbf{x}_1 \;=\; \lambda_1 (\mathbf{x}_1 \cdot \mathbf{x}_2). Hence (λ2λ1)(x1x2)=0(\lambda_2 - \lambda_1)(\mathbf{x}_1 \cdot \mathbf{x}_2) = 0. Since λ1λ2\lambda_1 \neq \lambda_2, we conclude x1x2=0\mathbf{x}_1 \cdot \mathbf{x}_2 = 0, i.e., x1x2\mathbf{x}_1 \perp \mathbf{x}_2. \blacksquare

Why this matters. Theorem 1184 is the engine that lets us build orthonormal eigenbases: eigenvectors from distinct eigenspaces are free orthogonality. Only inside a single eigenspace do we sometimes have to run Gram–Schmidt.


3.4 Theorem 1191 — Invertible symmetric matrices and reciprocal eigenvalues

Statement. Let AA be a symmetric matrix of order nn. (1) AA is invertible if and only if all of its eigenvalues are non-zero. (2) If AA is invertible, then λσ(A)    1λσ(A1)\lambda \in \sigma(A) \iff \tfrac{1}{\lambda} \in \sigma(A^{-1}).

Source: Proofs - Linear Algebra.pdf p.6 (handwritten theorem 1191).

Proof of (1). The product rule for eigenvalues gives λ1λ2λn=detA.\lambda_1 \cdot \lambda_2 \cdots \lambda_n = \det A. Therefore λ1,,λn0    detA0    A is invertible.\lambda_1,\ldots,\lambda_n \neq 0 \iff \det A \neq 0 \iff A \text{ is invertible}. \quad\checkmark

Proof of (2), direction λσ(A)1/λσ(A1)\lambda \in \sigma(A) \Rightarrow 1/\lambda \in \sigma(A^{-1}).

AA is invertible so detA0\det A \neq 0, hence λ0\lambda \neq 0 by (1). Let x0\mathbf{x} \neq \mathbf{0} be an eigenvector for λ\lambda: Ax=λxA \mathbf{x} = \lambda \mathbf{x}. Divide by λ\lambda: x=1λAx.\mathbf{x} = \tfrac{1}{\lambda} A \mathbf{x}. Multiply both sides on the left by A1A^{-1}: A1x=A1(1λAx)=1λ(A1A)x=1λIx=1λx.A^{-1} \mathbf{x} = A^{-1}\bigl(\tfrac{1}{\lambda} A \mathbf{x}\bigr) = \tfrac{1}{\lambda}(A^{-1} A)\mathbf{x} = \tfrac{1}{\lambda} I \mathbf{x} = \tfrac{1}{\lambda} \mathbf{x}. Thus A1x=1λxA^{-1}\mathbf{x} = \tfrac{1}{\lambda}\mathbf{x}, which says 1λ\tfrac{1}{\lambda} is an eigenvalue of A1A^{-1} with the same eigenvector x\mathbf{x}.

Proof of (2), converse. If 1λ\tfrac{1}{\lambda} is an eigenvalue of A1A^{-1}, apply the result just proved to the matrix A1A^{-1} (which is symmetric when AA is symmetric): 11/λ=λ\tfrac{1}{1/\lambda} = \lambda is an eigenvalue of (A1)1=A(A^{-1})^{-1} = A. \blacksquare

Remark (lecture note). The eigenvectors stay the same when we pass from AA to A1A^{-1}; only the eigenvalues are reciprocated.


3.5 Proposition 1201 — Bijective correspondence between quadratic forms and symmetric matrices

Statement. There is a bijective correspondence between quadratic forms f:RnRf : \mathbb{R}^n \to \mathbb{R} and symmetric matrices of order nn, given by f(x)=xTAx=i=1nj=1naijxixjxRn.f(\mathbf{x}) = \mathbf{x}^T A \mathbf{x} = \sum_{i=1}^n \sum_{j=1}^n a_{ij} x_i x_j \qquad \forall\,\mathbf{x}\in\mathbb{R}^n.

Source: Proofs - Linear Algebra.pdf p.7 (handwritten prop 1201). Also: lect3 p.6 proposition 1295 in the textbook numbering.

Proof.

(\Rightarrow) Symmetric AA yields a unique quadratic form ff.

Expand xTAx\mathbf{x}^T A \mathbf{x}: xTAx=i=1naiixi2+21i<jnaijxixj.\mathbf{x}^T A \mathbf{x} = \sum_{i=1}^n a_{ii} x_i^2 + 2\sum_{1 \le i < j \le n} a_{ij} x_i x_j. The doubling 2aijxixj2 \cdot a_{ij} x_i x_j comes from combining the (i,j)(i,j) and (j,i)(j,i) entries, which are equal because AA is symmetric. This is a sum of monomials of degree 2, i.e., a quadratic form. Setting f(x):=xTAxf(\mathbf{x}) := \mathbf{x}^T A \mathbf{x} gives a well-defined quadratic form.

(\Leftarrow) Quadratic form ff yields a unique symmetric AA.

Given a quadratic form ff, we want to write f(x)=xTAxf(\mathbf{x}) = \mathbf{x}^T A \mathbf{x} for some symmetric AA. Reading coefficients:

  • The coefficient of xi2x_i^2 must equal aiia_{ii}, fixing the diagonal.
  • The coefficient of the monomial xixjx_i x_j (for iji \neq j) must equal aij+ajia_{ij} + a_{ji}.

Requiring AA to be symmetric (aij=ajia_{ij} = a_{ji}) forces aij=aji=12(coefficient of xixj)a_{ij} = a_{ji} = \tfrac{1}{2}\bigl(\text{coefficient of } x_i x_j\bigr). So the symmetric matrix AA is uniquely determined by ff. \blacksquare

Practical rule. To build AA from ff: put coefficients of squared terms on the diagonal; split each mixed term evenly across the two symmetric off-diagonal positions.

Example (Sicconi lect3 p.7). A=(3112)f(x1,x2)=3x122x22+2x1x2A = \begin{pmatrix} 3 & 1 \\ 1 & -2 \end{pmatrix} \leadsto f(x_1, x_2) = 3x_1^2 - 2x_2^2 + 2 x_1 x_2 (the off-diagonal 1+11 + 1 becomes a +2x1x2+2 x_1 x_2 term).


3.6 Theorem 1207 — Positive definite matrices are invertible (and inverse is PD)

Statement. If AA is positive definite, then (1) AA is invertible; (2) A1A^{-1} is also positive definite.

Source: Proofs - Linear Algebra.pdf p.8 (handwritten theorem 1207).

Proof of (1). Suppose by contradiction AA is not invertible, so detA=0\det A = 0. Then the homogeneous system Ax=0A \mathbf{x} = \mathbf{0} admits a non-trivial solution x0\mathbf{x} \neq \mathbf{0}. But then xTAx=xT0=0\mathbf{x}^T A \mathbf{x} = \mathbf{x}^T \mathbf{0} = 0 with x0\mathbf{x} \neq \mathbf{0}, contradicting the definition of positive definite (which demands xTAx>0\mathbf{x}^T A \mathbf{x} > 0 for all non-zero x\mathbf{x}). Absurd. Hence detA0\det A \neq 0 and AA is invertible.

Proof of (2). Take any x0\mathbf{x} \neq \mathbf{0} and set y=A1x\mathbf{y} = A^{-1}\mathbf{x} (equivalently x=Ay\mathbf{x} = A \mathbf{y}). Since detA0\det A \neq 0, A1A^{-1} is also invertible, so y0\mathbf{y} \neq \mathbf{0}. We now compute xTA1x\mathbf{x}^T A^{-1} \mathbf{x}: xTA1x=xy=i=1nxiyi=i=1n(j=1naijyj)yi=i=1nj=1naijyiyj=yTAy.\mathbf{x}^T A^{-1} \mathbf{x} = \mathbf{x} \cdot \mathbf{y} = \sum_{i=1}^n x_i y_i = \sum_{i=1}^n \Bigl(\sum_{j=1}^n a_{ij} y_j\Bigr) y_i = \sum_{i=1}^n \sum_{j=1}^n a_{ij} y_i y_j = \mathbf{y}^T A \mathbf{y}. The substitution xi=jaijyjx_i = \sum_j a_{ij} y_j uses x=Ay\mathbf{x} = A\mathbf{y}. Since AA is positive definite and y0\mathbf{y} \neq \mathbf{0}: xTA1x=yTAy>0.\mathbf{x}^T A^{-1} \mathbf{x} = \mathbf{y}^T A \mathbf{y} > 0. This holds for every x0\mathbf{x} \neq \mathbf{0}, so A1A^{-1} is positive definite. \blacksquare

Semi-definite companion (Proof p.8, bottom). If AA is positive semi-definite, then AA is invertible iff AA is positive definite. (Proof via the equivalence xTAx=0    Ax=0\mathbf{x}^T A \mathbf{x} = 0 \iff A \mathbf{x} = \mathbf{0} which holds for semi-definite matrices.)


3.7 Spectral Theorem (Theorem 1291) — Orthogonal diagonalization of symmetric matrices

Statement. A symmetric matrix AA of order nn is orthogonally diagonalizable. That is, there exists an orthogonal matrix BB (so BT=B1B^T = B^{-1}) such that BTAB=Λ,B^T A B = \Lambda, where Λ=diag(λ1,,λn)\Lambda = \operatorname{diag}(\lambda_1, \ldots, \lambda_n) is the diagonal matrix of eigenvalues, each repeated according to its multiplicity. Equivalently, A=BΛBTA = B \Lambda B^T (spectral decomposition).

Source: lect3_linalg.pdf p.1 — Sicconi reproduces textbook Theorem 1291. (Proof is course-level, constructed by assembling an orthonormal eigenbasis as the columns of BB — see algorithm M-LA-6 in §5.)

Idea of construction.

  1. Find all eigenvalues λ1,,λn\lambda_1, \ldots, \lambda_n of AA (Theorem 1179).
  2. For each distinct eigenvalue, find a basis of its eigenspace (Theorem 1181); eigenvectors from different eigenspaces are automatically orthogonal (Theorem 1184).
  3. Inside each eigenspace, apply Gram–Schmidt if necessary to turn the basis orthogonal, then normalise every vector to length 1.
  4. Place the resulting nn orthonormal eigenvectors as the columns of BB. Then BB is orthogonal and BTAB=ΛB^T A B = \Lambda with the λi\lambda_i on the diagonal (ordered to match the columns).

Useful consequence. detA=λ1λn\det A = \lambda_1 \cdots \lambda_n (since A=B2Λ=Λ|A| = |B|^2 |\Lambda| = |\Lambda|, noting B=±1|B| = \pm 1 for orthogonal BB).


3.8 Sign of eigenvalues ↔ definiteness (Theorem 1210)

Statement. Let AA be a symmetric matrix of order nn. (i) AA is positive definite     \iff all eigenvalues are >0> 0. (ii) AA is positive semi-definite     \iff all eigenvalues are 0\ge 0. (iii) AA is negative definite     \iff all eigenvalues are <0< 0. (iv) AA is negative semi-definite     \iff all eigenvalues are 0\le 0. (v) AA is indefinite     \iff AA has at least one positive and one negative eigenvalue.

Source: Proofs - Linear Algebra.pdf pp.9–10 (handwritten theorem 1210).

Proof of (i), \Rightarrow: positive definite implies eigenvalues >0> 0.

Let λ\lambda be an eigenvalue with a normalised eigenvector x\mathbf{x} (so xx=1\mathbf{x} \cdot \mathbf{x} = 1 and x0\mathbf{x}\neq\mathbf{0}). Then λ=λ(xx)=xT(λx)=xTAx.\lambda = \lambda \cdot (\mathbf{x} \cdot \mathbf{x}) = \mathbf{x}^T (\lambda \mathbf{x}) = \mathbf{x}^T A \mathbf{x}. Because AA is positive definite and x0\mathbf{x} \neq \mathbf{0}, xTAx>0\mathbf{x}^T A \mathbf{x} > 0, hence λ>0\lambda > 0.

Proof of (i), \Leftarrow: eigenvalues >0> 0 implies AA positive definite.

Assume λ1,,λn>0\lambda_1, \ldots, \lambda_n > 0. By the Spectral Theorem there exists an orthogonal BB with BTAB=ΛB^T A B = \Lambda. Take any x0\mathbf{x} \neq \mathbf{0} and set y=BTx\mathbf{y} = B^T \mathbf{x}. Since BB (and hence BTB^T) is orthogonal, BTB^T is invertible and detBT0\det B^T \neq 0, so y0\mathbf{y} \neq \mathbf{0}. Note that By=BBTx=Ix=xB \mathbf{y} = B B^T \mathbf{x} = I \mathbf{x} = \mathbf{x}. Substituting: xTAx=(By)TA(By)=yT(BTAB)y=yTΛy=i=1nλiyi2.\mathbf{x}^T A \mathbf{x} = (B\mathbf{y})^T A (B\mathbf{y}) = \mathbf{y}^T (B^T A B) \mathbf{y} = \mathbf{y}^T \Lambda \mathbf{y} = \sum_{i=1}^n \lambda_i y_i^2. At least one yi0y_i \neq 0 (since y0\mathbf{y}\neq\mathbf{0}), and every λi>0\lambda_i > 0, so λiyi2>0\sum \lambda_i y_i^2 > 0. Hence xTAx>0\mathbf{x}^T A \mathbf{x} > 0 for every x0\mathbf{x}\neq\mathbf{0}, which says AA is positive definite.

For semi-definiteness, replace >> by \ge throughout (the proof only uses the sign of the λi\lambda_i). Negative (semi)definite is obtained by applying the result to A-A: AA negative definite     \iff A-A positive definite     \iff all eigenvalues of A-A are >0>0     \iff all eigenvalues of AA are <0<0. \blacksquare


3.9 Sylvester–Jacobi's criterion (Theorem 1217)

Statement. Let AA be a symmetric matrix of order nn and denote its leading (North-West) principal minors detA1,detA2,,detAn\det A_1, \det A_2, \ldots, \det A_n. Then: (1) AA is positive definite     \iff all leading principal minors are strictly positive: detA1>0, detA2>0, , detAn>0.\det A_1 > 0,\ \det A_2 > 0,\ \ldots,\ \det A_n > 0. (2) AA is negative definite     \iff the leading principal minors are all non-zero and alternate in sign starting with a strictly negative sign: detA1<0, detA2>0, detA3<0,detAn{>0n even<0n odd\det A_1 < 0,\ \det A_2 > 0,\ \det A_3 < 0, \ldots \qquad \det A_n \begin{cases} > 0 & n\text{ even}\\ < 0 & n\text{ odd}\end{cases} (3) If the leading principal minors are all non-zero but the sign sequence is neither (1) nor (2), then AA is indefinite.

Source: Proofs - Linear Algebra.pdf pp.13–14 (handwritten theorem 1217, with its preliminary result 1216 on general principal minors).

Proof of (1), \Rightarrow: positive definite implies all leading principal minors >0> 0.

By Theorem 1216 (preliminary), if AA is positive definite then every principal minor of AA is strictly positive: for any principal submatrix AiiA_{ii} and any xRn1\mathbf{x} \in \mathbb{R}^{n-1} extended to x~Rn\tilde{\mathbf{x}} \in \mathbb{R}^n by inserting a 00 in the ii-th position, xTAiix=x~TAx~>0x0,\mathbf{x}^T A_{ii} \mathbf{x} = \tilde{\mathbf{x}}^T A \tilde{\mathbf{x}} > 0 \quad \forall \mathbf{x}\neq\mathbf{0}, so AiiA_{ii} is positive definite and (by Theorem 1207) detAii>0\det A_{ii} > 0. In particular every leading principal minor is strictly positive.

Proof of (1), \Leftarrow: all leading principal minors >0> 0 implies positive definite.

Take any x0Rn\mathbf{x} \neq \mathbf{0} \in \mathbb{R}^n. Let CC be the upper-triangular matrix produced by Briosch's theorem (a course result we import), which is invertible since detC0\det C \neq 0. Set z=Cx0\mathbf{z} = C \mathbf{x} \neq \mathbf{0}. Then xTAx=k=1ndetAkdetAk1zk2,\mathbf{x}^T A \mathbf{x} = \sum_{k=1}^n \frac{\det A_k}{\det A_{k-1}}\, z_k^2, with the convention detA0=1\det A_0 = 1. By hypothesis every ratio detAk/detAk1\det A_k / \det A_{k-1} is strictly positive, and at least one zk2>0z_k^2 > 0, so xTAx>0\mathbf{x}^T A \mathbf{x} > 0. Hence AA is positive definite.

Proof of (2), \Rightarrow. Suppose AA is negative definite and set B:=AB := -A. Then BB is positive definite, so detBk>0\det B_k > 0 for every kk. But detAk=det(Bk)=(1)kdetBk,\det A_k = \det(-B_k) = (-1)^k \det B_k, which gives detAk<0\det A_k < 0 for kk odd and detAk>0\det A_k > 0 for kk even — exactly the alternating pattern (starting negative).

Proof of (2), \Leftarrow. Reverse the argument: given AA with the alternating sign pattern, set B:=AB := -A; check that every detBk=(1)kdetAk>0\det B_k = (-1)^k \det A_k > 0; conclude BB is positive definite by part (1); hence A=BA = -B is negative definite.

(3) Residual case. If the signs match neither (1) nor (2), then AA cannot be positive/negative (semi)definite (those imply the specific patterns), so AA must be indefinite. \blacksquare

Caveat (boxed by Sicconi). If some detAk=0\det A_k = 0, the leading-minors test is inconclusiveAA may be positive semi-definite, negative semi-definite, or indefinite. Then you must check all principal minors (not just the leading ones) or compute eigenvalues directly.


3.10 Auxiliary results worth remembering

Result on signs of diagonal entries (Theorem 1213). If AA is positive (semi)definite then aii>0a_{ii} > 0 (resp. 0\ge 0) for every ii. Proof: evaluate ff on the coordinate vector ei\mathbf{e}_i: f(ei)=eiTAei=aiif(\mathbf{e}_i) = \mathbf{e}_i^T A \mathbf{e}_i = a_{ii}.

Product rule detA=λi\det A = \prod \lambda_i. Follows from the Spectral Theorem: detA=det(BΛBT)=B2detΛ=1λ1λn\det A = \det(B\Lambda B^T) = |B|^2 \cdot \det \Lambda = 1 \cdot \lambda_1 \cdots \lambda_n.

Trace rule trA=λi\operatorname{tr} A = \sum \lambda_i. Trace is invariant under ABTABA \mapsto B^T A B for orthogonal BB, and trΛ=λi\operatorname{tr} \Lambda = \sum \lambda_i.

Invertibility of symmetric AA. If AA is symmetric and invertible, then A1A^{-1} is symmetric too: (A1)T=(AT)1=A1(A^{-1})^T = (A^T)^{-1} = A^{-1}.


§4. Worked Examples

Example 4.1 — Eigenvalues of a 2×22\times2 matrix via the shortcut

Let A=(1224)A = \begin{pmatrix} 1 & 2 \\ 2 & 4 \end{pmatrix}. (Source: Linear Algebra and Differential Calculus.pdf p.10.)

Using trace/determinant.

  • trA=1+4=5λ1+λ2=5\operatorname{tr} A = 1 + 4 = 5 \Rightarrow \lambda_1 + \lambda_2 = 5.
  • detA=1422=0λ1λ2=0\det A = 1 \cdot 4 - 2\cdot 2 = 0 \Rightarrow \lambda_1 \lambda_2 = 0.

So one eigenvalue is 00 and the other is 55: σ(A)={0,5}\sigma(A) = \{0, 5\}.

Cross-check with characteristic polynomial. pA(λ)=(1λ)(4λ)4=λ25λ=λ(λ5)p_A(\lambda) = (1-\lambda)(4-\lambda) - 4 = \lambda^2 - 5\lambda = \lambda(\lambda - 5), confirming λ{0,5}\lambda \in \{0, 5\}.

Example 4.2 — Eigenspaces of the same matrix

λ=0\lambda = 0. Solve (A0I)x=0(A - 0\cdot I)\mathbf{x} = \mathbf{0}: {x+2y=02x+4y=0x=2y.\begin{cases} x + 2y = 0 \\ 2x + 4y = 0 \end{cases} \Longrightarrow x = -2y. The second equation is dependent. Eigenspace W0={xR2:x=α(21), αR}.W_0 = \Bigl\{\mathbf{x} \in \mathbb{R}^2 : \mathbf{x} = \alpha \begin{pmatrix} -2 \\ 1 \end{pmatrix},\ \alpha \in \mathbb{R}\Bigr\}.

λ=5\lambda = 5. Solve (A5I)x=0(A - 5I)\mathbf{x} = \mathbf{0}: 4x+2y=0y=2x-4x + 2y = 0 \Rightarrow y = 2x. Eigenspace W5={xR2:x=α(12), αR}.W_5 = \Bigl\{\mathbf{x} \in \mathbb{R}^2 : \mathbf{x} = \alpha \begin{pmatrix} 1 \\ 2 \end{pmatrix},\ \alpha \in \mathbb{R}\Bigr\}.

Eigenbasis. Check orthogonality: (2,1)(1,2)=2+2=0(-2,1)\cdot(1,2) = -2+2 = 0 — orthogonal (as guaranteed by Theorem 1184 since 050\neq 5). Normalise (each vector has norm 5\sqrt{5}): u1=(2/51/5),u2=(1/52/5).\mathbf{u}_1 = \begin{pmatrix} -2/\sqrt{5} \\ 1/\sqrt{5} \end{pmatrix}, \qquad \mathbf{u}_2 = \begin{pmatrix} 1/\sqrt{5} \\ 2/\sqrt{5} \end{pmatrix}. These form an eigenbasis of R2\mathbb{R}^2. The orthogonal matrix B=[u1 u2]B = [\mathbf{u}_1\ \mathbf{u}_2] gives the spectral decomposition BTAB=Λ=(0005).B^T A B = \Lambda = \begin{pmatrix} 0 & 0 \\ 0 & 5 \end{pmatrix}.

Example 4.3 — Eigenvalues of a 3×3 (textbook TA question)

A=(304050403)A = \begin{pmatrix} -3 & 0 & -4 \\ 0 & -5 & 0 \\ -4 & 0 & 3 \end{pmatrix}. (Source: Linear Algebra and Differential Calculus.pdf p.12 — Question 4.)

Shortcut via trace + answer-choice elimination. trA=35+3=5\operatorname{tr} A = -3 - 5 + 3 = -5, so λ1+λ2+λ3=5\lambda_1 + \lambda_2 + \lambda_3 = -5. Among the MCQ options, only "5-5 with multiplicity 2 and 55 with multiplicity 1" satisfies the sum: 5+(5)+5=5-5 + (-5) + 5 = -5.

Cross-check with determinant. Expanding along row 2 (only the central term survives), detA=5det(3443)=5(916)=5(25)=125\det A = -5 \cdot \det\begin{pmatrix}-3&-4\\-4&3\end{pmatrix} = -5(-9 - 16) = -5(-25) = 125. And indeed (5)(5)5=125(-5)\cdot(-5)\cdot 5 = 125. \checkmark

Example 4.4 — Writing a quadratic form from a matrix (and vice versa)

From matrix. Given A=(3112)A = \begin{pmatrix} 3 & 1 \\ 1 & -2 \end{pmatrix}. (Source: lect3 p.7.) Reading off: a11=3a_{11} = 3 (coefficient of x12x_1^2), a22=2a_{22} = -2 (coefficient of x22x_2^2), a12+a21=2a_{12}+a_{21} = 2 (coefficient of x1x2x_1 x_2). Hence f(x1,x2)=3x122x22+2x1x2.f(x_1, x_2) = 3x_1^2 - 2x_2^2 + 2 x_1 x_2.

From form. Given f(x1,x2,x3)=x12+2x1x3+x223x2x3f(x_1, x_2, x_3) = -x_1^2 + 2 x_1 x_3 + x_2^2 - 3 x_2 x_3. (Source: lect3 p.8.)

  • Diagonal: a11=1a_{11} = -1, a22=1a_{22} = 1, a33=0a_{33} = 0.
  • Off-diagonal: coefficient of x1x2x_1 x_2 is 0a12=00 \Rightarrow a_{12} = 0. Coefficient of x1x3x_1 x_3 is 2a13=12 \Rightarrow a_{13} = 1. Coefficient of x2x3x_2 x_3 is 3a23=3/2-3 \Rightarrow a_{23} = -3/2. A=(101013/213/20).A = \begin{pmatrix} -1 & 0 & 1 \\ 0 & 1 & -3/2 \\ 1 & -3/2 & 0 \end{pmatrix}.

Example 4.5 — Classifying a quadratic form two ways

Given f(x,y)=3x2+2xy+5y2f(x,y) = 3x^2 + 2xy + 5y^2, i.e., A=(3115)A = \begin{pmatrix} 3 & 1 \\ 1 & 5 \end{pmatrix}. (Source: lect2 p.8.)

Via Sylvester–Jacobi.

  • detA1=3>0\det A_1 = 3 > 0
  • detA2=detA=151=14>0\det A_2 = \det A = 15 - 1 = 14 > 0

Both leading principal minors positive \Rightarrow positive definite \Rightarrow 0\mathbf{0} is a strict (unique) global minimum of ff.

Via eigenvalue signs. trA=8\operatorname{tr} A = 8, detA=14\det A = 14. Both >0>0, so both eigenvalues are positive, giving the same conclusion.

Example 4.6 — The three 3×33\times 3 classification cases (lecture demo)

(Source: Linear Algebra and Differential Calculus.pdf pp.18–19.)

(a) Indefinite. f(x,y)=2x26xy+y2f(x,y) = 2x^2 - 6xy + y^2, A=(2331)A = \begin{pmatrix}2 & -3 \\ -3 & 1\end{pmatrix}. detA1=2>0\det A_1 = 2 > 0, detA2=29=7<0\det A_2 = 2 - 9 = -7 < 0. Sign pattern (+,)(+, -) matches neither positive nor negative definite, so AA is indefinite and 0\mathbf{0} is a saddle point.

(b) Indefinite with a zero leading minor. f(x,y,z)=x2+3xzyzf(x,y,z) = x^2 + 3xz - yz, A=(103/2001/23/21/20)A = \begin{pmatrix}1 & 0 & 3/2 \\ 0 & 0 & -1/2 \\ 3/2 & -1/2 & 0\end{pmatrix}. detA1=1\det A_1 = 1, detA2=0\det A_2 = 0, detA3=1/4\det A_3 = -1/4. Since detA30\det A_3 \neq 0 but the Sylvester pattern fails (a zero appears in the leading sequence yet the final determinant is non-zero and negative), AA is indefinite — not definite or semi-definite.

(c) Positive semi-definite (not definite). f(x,y,z)=x2+2y2+3z2+2xz+4yzf(x,y,z) = x^2 + 2y^2 + 3z^2 + 2xz + 4yz, A=(101022123)A = \begin{pmatrix}1 & 0 & 1 \\ 0 & 2 & 2 \\ 1 & 2 & 3\end{pmatrix}. Leading: detA1=1\det A_1 = 1, detA2=2\det A_2 = 2, detA3=0\det A_3 = 0. The zero means Sylvester–Jacobi is inconclusive. Check all principal minors: order-1 minors are 1,2,31, 2, 3 (all 0\ge 0); order-2 principal minors are det(2223)=2\det\begin{pmatrix}2&2\\2&3\end{pmatrix}=2, det(1113)=2\det\begin{pmatrix}1&1\\1&3\end{pmatrix}=2, det(1002)=2\det\begin{pmatrix}1&0\\0&2\end{pmatrix}=2 (all 0\ge 0); order-3 is detA=0\det A = 0. All principal minors 0\ge 0 \Rightarrow positive semi-definite \Rightarrow 0\mathbf{0} is a non-strict (weak) minimizer.

Example 4.7 — Full diagonalization of a 3×33\times 3 (lecture example)

A=(112112224)A = \begin{pmatrix}1 & -1 & 2 \\ -1 & 1 & -2 \\ 2 & -2 & 4\end{pmatrix}. Find an eigenbasis of R3\mathbb{R}^3. (Source: Linear Algebra and Differential Calculus.pdf pp.16–17.)

Step 1. Characteristic polynomial. Expanding and simplifying (full calculation in the source): pA(λ)=λ3+6λ2=λ2(λ+6)p_A(\lambda) = -\lambda^3 + 6\lambda^2 = \lambda^2 (-\lambda + 6). Roots: λ=0\lambda = 0 (multiplicity 2), λ=6\lambda = 6 (multiplicity 1).

Step 2. Eigenspaces.

  • λ=0\lambda = 0: solve Ax=0A \mathbf{x} = \mathbf{0}. Reduces to xy+2z=0x - y + 2z = 0, so x=y2zx = y - 2z. Two-parameter family: (xyz)=y(110)+z(201).\begin{pmatrix}x\\y\\z\end{pmatrix} = y \begin{pmatrix}1\\1\\0\end{pmatrix} + z\begin{pmatrix}-2\\0\\1\end{pmatrix}. Basis of W0W_0: x1=(1,1,0)T\mathbf{x}_1 = (1,1,0)^T, x2=(2,0,1)T\mathbf{x}_2 = (-2,0,1)^T.
  • λ=6\lambda = 6: solve (A6I)x=0(A - 6I)\mathbf{x} = \mathbf{0}. Gaussian elimination gives x=yx = -y, z=2yz = -2y. Basis of W6W_6: x3=(1,1,2)T\mathbf{x}_3 = (-1, 1, -2)^T.

Step 3. Eigenbasis. x1\mathbf{x}_1 and x3\mathbf{x}_3 sit in different eigenspaces so are orthogonal (and so are x2,x3\mathbf{x}_2, \mathbf{x}_3). But x1,x2\mathbf{x}_1, \mathbf{x}_2 are in the same eigenspace W0W_0 and happen NOT to be orthogonal: x1x2=2\mathbf{x}_1 \cdot \mathbf{x}_2 = -2. Apply Gram–Schmidt:

  • Normalise x1\mathbf{x}_1: x1=2\|\mathbf{x}_1\| = \sqrt{2}, so u1=(1/2,1/2,0)T\mathbf{u}_1 = (1/\sqrt{2}, 1/\sqrt{2}, 0)^T.
  • Orthogonal auxiliary: z2=x2(x2u1)u1=(2,0,1)T(2/2)(1/2,1/2,0)T=(1,1,1)T\mathbf{z}_2 = \mathbf{x}_2 - (\mathbf{x}_2 \cdot \mathbf{u}_1)\mathbf{u}_1 = (-2, 0, 1)^T - (-2/\sqrt{2})(1/\sqrt{2}, 1/\sqrt{2}, 0)^T = (-1, 1, 1)^T. Then z2=3\|\mathbf{z}_2\| = \sqrt{3}, giving u2=(1/3,1/3,1/3)T\mathbf{u}_2 = (-1/\sqrt{3}, 1/\sqrt{3}, 1/\sqrt{3})^T.
  • Normalise x3\mathbf{x}_3: x3=6\|\mathbf{x}_3\| = \sqrt{6}, u3=(1/6,1/6,2/6)T\mathbf{u}_3 = (-1/\sqrt{6}, 1/\sqrt{6}, -2/\sqrt{6})^T.

{u1,u2,u3}\{\mathbf{u}_1, \mathbf{u}_2, \mathbf{u}_3\} is the eigenbasis of R3\mathbb{R}^3 for AA, and B=[u1 u2 u3]B = [\mathbf{u}_1\ \mathbf{u}_2\ \mathbf{u}_3] satisfies BTAB=diag(0,0,6)B^T A B = \operatorname{diag}(0, 0, 6).


§5. Solution Methods

Each method is a named algorithm with input → steps → output → common pitfalls. Cross-references in the right column point to the May 2024 exam problems where the method is applied.


M-LA-1 — Verify a candidate vector is an eigenvector

Used on: May 2024 MCQ1 (Mode A), MCQ6 (Mode B).

Input. A symmetric matrix AA of order nn and a candidate vector vRn\mathbf{v} \in \mathbb{R}^n.

Steps.

  1. Compute the product AvA \mathbf{v}.
  2. Check whether the resulting vector is a scalar multiple of v\mathbf{v}: look for a single λR\lambda \in \mathbb{R} so that Av=λvA \mathbf{v} = \lambda \mathbf{v}.
  3. If yes, v\mathbf{v} is an eigenvector with eigenvalue λ\lambda. If no, v\mathbf{v} is not an eigenvector.

Output. "Yes, eigenvector with eigenvalue λ\lambda" or "No".

Pitfalls.

  • Forgetting that the zero vector is never an eigenvector (even though A0=0A \mathbf{0} = \mathbf{0} trivially), by definition.
  • On MCQs, always test the shortest vector first (fewer arithmetic operations). The May 2024 MCQ1 correct answer is (0,1,0)T(0,1,0)^T because A(0,1,0)T=(0,2,0)T=2(0,1,0)TA(0,1,0)^T = (0,-2,0)^T = -2 \cdot (0,1,0)^T.

M-LA-2 — Compute the eigenvalues of a 3×33\times 3 (or larger) matrix

Used on: Final exercises Q1, May 2024 Q7/Q12 (eigenvalue determination for A=11TA = \mathbf{1}\mathbf{1}^T).

Input. Symmetric matrix AA of order nn.

Steps.

  1. Write down AλIA - \lambda I (subtract λ\lambda from every diagonal entry).
  2. Compute pA(λ)=det(AλI)p_A(\lambda) = \det(A - \lambda I) by Laplace expansion along the simplest row/column (prefer rows/columns with zeros).
  3. Factor pA(λ)p_A(\lambda) over R\mathbb{R}. Useful tricks:
    • If a row/column of AλIA - \lambda I becomes identically zero at λ=λ\lambda = \lambda^\star, then λ\lambda^\star is a root and pAp_A is divisible by (λλ)(\lambda - \lambda^\star).
    • Rational-root test: integer roots of pAp_A divide the constant term pA(0)=detAp_A(0) = \det A (up to sign).
    • Sum/product shortcut: λi=trA\sum \lambda_i = \operatorname{tr} A and λi=detA\prod \lambda_i = \det A; use as a sanity check after solving.
  4. Solve pA(λ)=0p_A(\lambda) = 0 and list the eigenvalues with multiplicities. The total of multiplicities must equal nn.

Output. Spectrum σ(A)={λ1,,λk}\sigma(A) = \{\lambda_1, \ldots, \lambda_k\} with algebraic multiplicities.

Pitfalls.

  • On May 2024 Q7, A=(111111111)A = \begin{pmatrix}1&1&1\\1&1&1\\1&1&1\end{pmatrix} has trA=3\operatorname{tr} A = 3 and detA=0\det A = 0, and since AA is rank-1, λ=0\lambda = 0 has multiplicity n1=2n-1 = 2 and the remaining eigenvalue is trA=3\operatorname{tr} A = 3: σ(A)={0 (mult. 2),3 (mult. 1)}\sigma(A) = \{0\text{ (mult.\ 2)}, 3\text{ (mult.\ 1)}\}. Always check your calculation against these totals.
  • When the characteristic polynomial is a cubic with an obvious factor of λ\lambda (happens whenever detA=0\det A = 0), factor it out first and solve the remaining quadratic.

M-LA-3 — Find a basis of an eigenspace

Used on: May 2024 Q7 (Mode A) / Q12 (Mode B).

Input. Symmetric matrix AA, a specific eigenvalue λ\lambda.

Steps.

  1. Write the system (AλI)x=0(A - \lambda I)\mathbf{x} = \mathbf{0}.
  2. Gaussian eliminate. Expect rank(AλI)=nμ\operatorname{rank}(A - \lambda I) = n - \mu where μ\mu is the algebraic multiplicity of λ\lambda; you'll get μ\mu free parameters.
  3. Express the solution set as a span of μ\mu independent vectors — those are the basis of WλW_\lambda.

Output. A basis {v1,,vμ}\{\mathbf{v}_1, \ldots, \mathbf{v}_\mu\} of WλW_\lambda.

Pitfalls.

  • dimWλ\dim W_\lambda must equal the algebraic multiplicity of λ\lambda (this is the symmetric-matrix-always-diagonalizable fact). If you find fewer independent vectors, you've made an arithmetic mistake.
  • Don't forget that eigenspaces contain 0\mathbf{0} even though 0\mathbf{0} is not an eigenvector.

M-LA-4 — Build a basis of eigenvectors for Rn\mathbb{R}^n

Used on: May 2024 Q7(c) / Q12(c) — explicit request "determine a basis of R3\mathbb{R}^3 constituted by eigenvectors of AA".

Input. Symmetric matrix AA of order nn.

Steps.

  1. Find all eigenvalues (M-LA-2).
  2. For each distinct eigenvalue, find a basis of its eigenspace (M-LA-3).
  3. Union of all eigenspace bases is a basis of Rn\mathbb{R}^n (since bases from distinct eigenspaces are automatically linearly independent by Theorem 1184, and the total count equals nn by the symmetric spectral structure).
  4. If the question asks "not necessarily orthogonal or orthonormal" (as in May 2024 Q7(c)), you can stop here.
  5. If the question asks for an eigenbasis (orthonormal), continue with M-LA-6.

Output. A list of nn linearly independent eigenvectors of AA spanning Rn\mathbb{R}^n.

Pitfalls.

  • Inside a single eigenspace of multiplicity 2\ge 2, the basis vectors are not automatically orthogonal. Theorem 1184 only applies to eigenvectors from distinct eigenspaces.

M-LA-5 — Classify a quadratic form (eigenvalue method vs Sylvester–Jacobi)

Used on: May 2024 Q7/Q12, Final exercises Q1.

Input. A symmetric matrix AA (or a quadratic form from which you extract AA via M-LA-rules in §2.11).

Decision tree — which test to use.

  • Prefer Sylvester–Jacobi when the matrix is small (2×22\times2 or 3×33\times3) and the leading principal minors are easy to compute.
  • Prefer the eigenvalue-sign test when (i) the matrix has an obvious rank-deficiency (so some eigenvalue is automatically 00), (ii) you're told the eigenvalues or can get them cheaply, (iii) a leading principal minor is 00 (Sylvester is inconclusive in that case), or (iv) you're also asked the sign of eigenvalues.

Steps (Sylvester–Jacobi branch).

  1. Compute detA1,detA2,,detAn\det A_1, \det A_2, \ldots, \det A_n.
  2. Match against the pattern:
    • All >0>0 \Rightarrow positive definite.
    • Alternating starting with <0<0 (,+,,+,-, +, -, +, \ldots) \Rightarrow negative definite.
    • All non-zero but neither pattern \Rightarrow indefinite.
    • Some leading minor is 00 and detA0\det A \neq 0 \Rightarrow indefinite (since neither positive nor negative definite, and not semi-definite since detA0\det A\neq 0 for semi-definite-not-definite is impossible).
    • Some leading minor is 00 and detA=0\det A = 0 \Rightarrow inconclusive — check ALL principal minors.
  3. If inconclusive, evaluate every principal minor (not just leading): all 0\ge 0 \Rightarrow positive semi-definite; all of the signed form (0(\le 0 for odd order, 0\ge 0 for even order) \Rightarrow negative semi-definite; otherwise indefinite.

Steps (Eigenvalue-sign branch).

  1. Compute the eigenvalues (M-LA-2).
  2. Apply Theorem 1210: all >0>0 → positive definite, all 0\ge 0 → positive semi-definite, mixed → indefinite, etc.

Output. One of: positive definite, positive semi-definite (non-def), negative definite, negative semi-definite (non-def), indefinite.

Pitfalls (critical — these are exam traps!).

  • Negative definite needs ALTERNATING signs, not all negative. (,+,,)(-, +, -, \ldots) — starting negative, alternating.
  • Sylvester fails with leading zeros. If detAk=0\det A_k = 0 for any kk, do NOT conclude indefinite from that alone; it may be semi-definite. Switch to the full principal-minor test (all minors, not just leading).
  • Semi-definite ≠ definite. May 2024 Q7 is the canonical trap: A=11TA = \mathbf{1}\mathbf{1}^T has eigenvalues 0,0,30, 0, 3, all 0\ge 0, so AA is positive semi-definite — but not positive definite because λ=0\lambda = 0 occurs. The quadratic form f(x)=(x1+x2+x3)2f(\mathbf{x}) = (x_1+x_2+x_3)^2 vanishes on the hyperplane x1+x2+x3=0x_1+x_2+x_3 = 0, confirming it is not strictly positive.

M-LA-6 — Diagonalize a symmetric matrix (full spectral decomposition)

Used on: conceptual backbone of quadratic-form sign analysis; also direct exam topics when a question asks for BTAB=ΛB^T A B = \Lambda.

Input. Symmetric matrix AA of order nn.

Steps.

  1. Compute eigenvalues (M-LA-2).
  2. For each eigenvalue, find a basis of the eigenspace (M-LA-3).
  3. Orthogonalise within each eigenspace (Gram–Schmidt): keep the first basis vector, and for each subsequent one vk\mathbf{v}_k subtract its projections onto the already-processed vectors.
  4. Normalise every vector to unit length. The resulting nn vectors form an eigenbasis.
  5. Assemble BB by stacking the eigenbasis as columns; BB is automatically orthogonal (BT=B1B^T = B^{-1}).
  6. The diagonal matrix Λ=BTAB\Lambda = B^T A B has the eigenvalues on the diagonal, ordered to match the column order in BB.

Output. Orthogonal BB and diagonal Λ\Lambda with BTAB=ΛB^T A B = \Lambda (equivalently A=BΛBTA = B \Lambda B^T).

Pitfalls.

  • Don't forget to normalise after Gram–Schmidt. Orthogonalising preserves orthogonality but doesn't scale to unit length.
  • If an eigenspace already has an orthonormal basis (e.g., the spectral theorem delivered pre-orthogonal basis vectors), skip Gram–Schmidt — it's idempotent but wasted work.
  • Order matters: the Λ\Lambda you produce has the eigenvalues in the order the columns of BB are arranged. Swap two columns of BB ↔ swap two diagonal entries of Λ\Lambda.

§6. Practice Problems with Solutions

Problem 6.1 — TA session 1, exercise 1: positive-definite quadratic form parametric

Consider f:R2Rf : \mathbb{R}^2 \to \mathbb{R} with f(x)=(α21)x12+(2α)x22f(\mathbf{x}) = (\alpha^2 - 1)x_1^2 + (2 - \alpha)x_2^2, αR\alpha \in \mathbb{R}. The quadratic form is positive definite in R2\mathbb{R}^2 if and only if … ?

(A) α>1|\alpha| > 1, (B) 1<α<2-1 < \alpha < 2, (C) α<2\alpha < 2, (D) none of the preceding.

Solution. Associated matrix A=(α21002α)A = \begin{pmatrix} \alpha^2 - 1 & 0 \\ 0 & 2 - \alpha \end{pmatrix}. Apply Sylvester–Jacobi:

  • detA1=α21>0    α<1 or α>1\det A_1 = \alpha^2 - 1 > 0 \iff \alpha < -1 \text{ or } \alpha > 1.
  • detA2=detA=(α21)(2α)>0\det A_2 = \det A = (\alpha^2 - 1)(2 - \alpha) > 0.

Combining: α21>0\alpha^2 - 1 > 0 AND 2α>02 - \alpha > 0     \iff (α<1\alpha<-1 or α>1\alpha>1) AND α<2\alpha < 2     \iff α<1\alpha<-1 OR 1<α<21 < \alpha < 2.

So AA is positive definite iff α<1\alpha < -1 or 1<α<21 < \alpha < 2. Answer: (D) none of the preceding.

Source: TA41_lin_alg.pdf pp.1–2.


Problem 6.2 — May 2024 MCQ1 (Mode A) / MCQ6 (Mode B): which vector is an eigenvector?

Which of the following vectors is an eigenvector of A=(201021103)A = \begin{pmatrix}2&0&1\\0&-2&-1\\1&0&3\end{pmatrix}?

(A) (2,0,1)T(2, 0, -1)^T, (B) (3,1,2)T(3, 1, 2)^T, (C) (0,1,0)T(0, 1, 0)^T, (D) none of the others.

Solution. Apply M-LA-1. Testing v=(0,1,0)T\mathbf{v} = (0,1,0)^T: Av=(20+01+1000+(2)1+(1)010+01+30)=(020)=2(010)=2v.A \mathbf{v} = \begin{pmatrix} 2\cdot 0 + 0\cdot 1 + 1\cdot 0 \\ 0\cdot 0 + (-2)\cdot 1 + (-1)\cdot 0 \\ 1\cdot 0 + 0\cdot 1 + 3\cdot 0\end{pmatrix} = \begin{pmatrix}0\\-2\\0\end{pmatrix} = -2 \begin{pmatrix}0\\1\\0\end{pmatrix} = -2 \mathbf{v}. So (0,1,0)T(0,1,0)^T is an eigenvector with eigenvalue 2-2. Answer: (C).

Source: General_24524_ENG_SOL.pdf p.1 MCQ1.


Problem 6.3 — May 2024 Q7 (Mode A) / Q12 (Mode B): Sicconi's flagship exam problem

(a) State the theorem on the correspondence between symmetric matrices AA of order nn and quadratic forms q:RnRq : \mathbb{R}^n \to \mathbb{R}. (b) Consider the matrix A=(111111111).A = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{pmatrix}. Write the quadratic form associated to AA and prove that this form is positive semi-definite but not positive definite. (c) Determine a basis of R3\mathbb{R}^3 constituted by eigenvectors of AA (not necessarily orthogonal or orthonormal).

Solution.

(a) Proposition 1201: there is a bijective correspondence between quadratic forms q:RnRq : \mathbb{R}^n \to \mathbb{R} and symmetric matrices AA of order nn, given by q(x)=xTAx=i,jaijxixjq(\mathbf{x}) = \mathbf{x}^T A \mathbf{x} = \sum_{i,j} a_{ij} x_i x_j.

(b) Read off the quadratic form: q(x1,x2,x3)=x12+x22+x32+2x1x2+2x1x3+2x2x3=(x1+x2+x3)2.q(x_1, x_2, x_3) = x_1^2 + x_2^2 + x_3^2 + 2 x_1 x_2 + 2 x_1 x_3 + 2 x_2 x_3 = (x_1 + x_2 + x_3)^2. Eigenvalues via characteristic polynomial. Note AA is rank 1 (every row equals (1,1,1)(1,1,1)), so λ=0\lambda = 0 has multiplicity 2. The remaining eigenvalue equals trA=3\operatorname{tr} A = 3. Formally: det(AλI)=(1λ)3+23(1λ)=3λ+3λ2λ3+3λ=λ2(λ3),\det(A - \lambda I) = (1-\lambda)^3 + 2 - 3(1-\lambda) = -3\lambda + 3\lambda^2 - \lambda^3 + 3\lambda = -\lambda^2(\lambda - 3), giving λ1=0\lambda_1 = 0 (multiplicity 2), λ2=3\lambda_2 = 3 (multiplicity 1). Since all eigenvalues are 0\ge 0, by Theorem 1210(ii) AA is positive semi-definite. Since λ=0\lambda = 0 occurs (not all strictly positive), AA is not positive definite. (Concretely, q(1,1,0)=0q(1, -1, 0) = 0 with (1,1,0)T0(1,-1,0)^T \neq \mathbf{0}.)

(c) Bases of the eigenspaces.

  • W0W_0: solve Ax=0A\mathbf{x} = \mathbf{0}, which reduces to x1+x2+x3=0x_1 + x_2 + x_3 = 0, equivalently x3=x1x2x_3 = -x_1 - x_2. A basis is {(1,0,1)T,(0,1,1)T}\{(1, 0, -1)^T, (0, 1, -1)^T\}.
  • W3W_3: solve (A3I)x=0(A - 3I)\mathbf{x} = \mathbf{0}: {2x1+x2+x3=0x12x2+x3=0x1+x22x3=0\begin{cases} -2x_1 + x_2 + x_3 = 0 \\ x_1 - 2x_2 + x_3 = 0 \\ x_1 + x_2 - 2x_3 = 0 \end{cases} From the first: x2=2x1x3x_2 = 2x_1 - x_3. Substitute into the second: x12(2x1x3)+x3=0x3=x1x_1 - 2(2x_1 - x_3) + x_3 = 0 \Rightarrow x_3 = x_1. Then x2=x1x_2 = x_1. Basis: (1,1,1)T(1, 1, 1)^T.

Eigenbasis of R3\mathbb{R}^3 (union by M-LA-4): {(1,0,1)T, (0,1,1)T, (1,1,1)T}\{(1, 0, -1)^T,\ (0, 1, -1)^T,\ (1, 1, 1)^T\}.

Source: General_24524_ENG_SOL.pdf p.5 Q7 (Mode A) / p.11 Q12 (Mode B).


Problem 6.4 — Final exercises Q1: quadratic form with parameter kk

Given f(x,y,z)=kx2+ky24z22xzf(x, y, z) = k x^2 + k y^2 - 4 z^2 - 2 x z, kRk \in \mathbb{R}. (a) Find conditions on kk such that the origin is a strong global maximum point. (b) Determine the sign of the eigenvalues for k=1k = -1. (c) Set k=0k = 0, verify that λ=0\lambda = 0 is an eigenvalue and find the kernel of the linear function associated to the matrix AA.

Solution.

(a) The origin is always a stationary point of a quadratic form (f(0)=0f(\mathbf{0}) = 0; f(0)=0\nabla f(\mathbf{0}) = \mathbf{0}). For it to be a strong global maximum, we need f(x)=xTAx<0f(\mathbf{x}) = \mathbf{x}^T A \mathbf{x} < 0 for every x0\mathbf{x} \neq \mathbf{0}, i.e., AA must be negative definite. The associated matrix (half the mixed coefficient on the off-diagonal): A=(k010k0104).A = \begin{pmatrix} k & 0 & -1 \\ 0 & k & 0 \\ -1 & 0 & -4 \end{pmatrix}. Apply Sylvester–Jacobi (negative-definite pattern: detA1<0\det A_1 < 0, detA2>0\det A_2 > 0, detA3<0\det A_3 < 0):

  • detA1=k<0\det A_1 = k < 0.
  • detA2=kk0=k2>0\det A_2 = k \cdot k - 0 = k^2 > 0 — requires k0k \neq 0.
  • detA3=detA\det A_3 = \det A. Expand along row 2: detA=kdet(k114)=k(4k1)=4k2k\det A = k \cdot \det\begin{pmatrix}k & -1 \\ -1 & -4\end{pmatrix} = k(-4k - 1) = -4k^2 - k. Need 4k2k<0-4k^2 - k < 0, i.e., k(4k+1)>0k(4k+1) > 0     \iff k>0k > 0 or k<1/4k < -1/4.

Combining all three: k<0k < 0 AND k0k \neq 0 AND (k>0k > 0 or k<1/4k < -1/4) k<14\Longrightarrow k < -\tfrac{1}{4}.

Answer. AA is negative definite     k<14\iff k < -\tfrac{1}{4}, in which case the origin is a strong global max.

(b) For k=1k = -1, the three conditions above become 1<0 -1 < 0\ \checkmark, 1>0 1 > 0\ \checkmark, 41(1)=3<0 -4 \cdot 1 - (-1) = -3 < 0\ \checkmark. AA is negative definite. By Theorem 1210(iii), all three eigenvalues are <0< 0: λ1,λ2,λ3<0\lambda_1, \lambda_2, \lambda_3 < 0.

(c) For k=0k = 0, A=(001000104)A = \begin{pmatrix} 0 & 0 & -1 \\ 0 & 0 & 0 \\ -1 & 0 & -4 \end{pmatrix}. Characteristic polynomial: expanding along row 2, pA(λ)=λdet(λ114λ)=λ(λ(4+λ)1)=λ(λ2+4λ1)p_A(\lambda) = -\lambda \cdot \det\begin{pmatrix}-\lambda & -1 \\ -1 & -4-\lambda\end{pmatrix} = -\lambda\bigl(\lambda(4+\lambda) - 1\bigr) = -\lambda(\lambda^2 + 4\lambda - 1). Clearly λ=0\lambda = 0 is a root (factor of λ\lambda), so λ=0\lambda = 0 is an eigenvalue of AA.

Kernel of AA = W0W_0: solve Ax=0A\mathbf{x} = \mathbf{0}: {x3=00=0x14x3=0x3=0, x1=0, x2 free.\begin{cases} -x_3 = 0 \\ 0 = 0 \\ -x_1 - 4 x_3 = 0 \end{cases} \Longrightarrow x_3 = 0,\ x_1 = 0,\ x_2 \text{ free}. Hence W0={xR3:x=x2(010), x2R}.W_0 = \Bigl\{\mathbf{x} \in \mathbb{R}^3 : \mathbf{x} = x_2 \begin{pmatrix}0\\1\\0\end{pmatrix},\ x_2 \in \mathbb{R}\Bigr\}.

Source: Final exercises solution.pdf pp.1–3 Problem 1.


Problem 6.5 — TA session 1 exercise 2: characteristic polynomial & eigenbasis of a 2×22\times 2

Given A=(1223)A = \begin{pmatrix} 1 & -2 \\ -2 & 3 \end{pmatrix}, compute the characteristic polynomial, find the eigenvalues and eigenvectors, determine a basis of R2\mathbb{R}^2, and decide whether it is an eigenbasis.

Solution.

Characteristic polynomial. pA(λ)=λ2(trA)λ+detA=λ24λ+(34)=λ24λ1p_A(\lambda) = \lambda^2 - (\operatorname{tr} A)\lambda + \det A = \lambda^2 - 4\lambda + (3 - 4) = \lambda^2 - 4\lambda - 1.

Eigenvalues. Solving λ24λ1=0\lambda^2 - 4\lambda - 1 = 0 by the quadratic formula: λ=4±16+42=2±5\lambda = \tfrac{4 \pm \sqrt{16+4}}{2} = 2 \pm \sqrt{5}. So σ(A)={25, 2+5}\sigma(A) = \{2 - \sqrt{5},\ 2 + \sqrt{5}\}.

Eigenvectors. For λ1=25\lambda_1 = 2 - \sqrt{5}: (Aλ1I)x=0(A - \lambda_1 I)\mathbf{x} = \mathbf{0} gives (1+5)x12x2=0(-1 + \sqrt{5})x_1 - 2 x_2 = 0, i.e., x1=251x2x_1 = \tfrac{2}{\sqrt{5}-1} x_2. Eigenvector x=(2/(51), 1)T\mathbf{x} = (2/(\sqrt{5}-1),\ 1)^T.

For λ2=2+5\lambda_2 = 2 + \sqrt{5}: similarly (15)x12x2=0x1=21+5x2(-1 - \sqrt{5})x_1 - 2 x_2 = 0 \Rightarrow x_1 = -\tfrac{2}{1+\sqrt{5}} x_2. Eigenvector x=(2/(1+5), 1)T\mathbf{x} = (-2/(1+\sqrt{5}),\ 1)^T.

Basis of R2\mathbb{R}^2. Since λ1λ2\lambda_1 \neq \lambda_2, the two eigenvectors are linearly independent (and in fact orthogonal by Theorem 1184). They therefore form a basis of R2\mathbb{R}^2.

Is it an eigenbasis? An eigenbasis requires the vectors to be orthonormal. Check norms — clearly neither has norm 1. So this basis is not an eigenbasis, but an orthonormal eigenbasis is obtained by dividing each vector by its norm.

Source: TA41_lin_alg.pdf pp.3–7 exercise 2.


Problem 6.6 — TA session 1 exercise 7: write the quadratic form from a matrix

Given A=(111120103)A = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 2 & 0 \\ 1 & 0 & 3 \end{pmatrix}, provide the analytical expression of the quadratic form associated to AA.

Solution. AA is symmetric. Read coefficients: diagonal gives x12+2x22+3x32x_1^2 + 2 x_2^2 + 3 x_3^2; each off-diagonal contributes 21xixj2 \cdot 1 \cdot x_i x_j because aij+aji=2a_{ij} + a_{ji} = 2 for every mixed pair (i,j)=(1,2),(1,3)(i,j) = (1,2), (1,3), and 202 \cdot 0 for (2,3)(2,3). Therefore f(x1,x2,x3)=x12+2x22+3x32+2x1x2+2x1x3+0x2x3=x12+2x22+3x32+2x1x2+2x1x3.f(x_1, x_2, x_3) = x_1^2 + 2 x_2^2 + 3 x_3^2 + 2 x_1 x_2 + 2 x_1 x_3 + 0\cdot x_2 x_3 = x_1^2 + 2 x_2^2 + 3 x_3^2 + 2 x_1 x_2 + 2 x_1 x_3.

Source: TA41_lin_alg.pdf p.8 exercise 7.


Problem 6.7 — TA session 1 exercise 8: classify four quadratic forms

Establish the sign of: (a) f(x)=x12+x22+αx32f(\mathbf{x}) = x_1^2 + x_2^2 + \alpha x_3^2, (b) f(x)=2x12+2x22+2x322x1x22x1x32x2x3f(\mathbf{x}) = 2x_1^2 + 2x_2^2 + 2x_3^2 - 2x_1 x_2 - 2 x_1 x_3 - 2 x_2 x_3, (c) f(x)=2x128x1x2+2x22+2x32f(\mathbf{x}) = 2x_1^2 - 8 x_1 x_2 + 2 x_2^2 + 2 x_3^2, (d) f(x)=2x12+5x22+5x32+4x1x24x1x38x2x3f(\mathbf{x}) = 2x_1^2 + 5 x_2^2 + 5 x_3^2 + 4 x_1 x_2 - 4 x_1 x_3 - 8 x_2 x_3.

Solution (part (a) — worked fully; others use the same method).

A=diag(1,1,α)A = \operatorname{diag}(1, 1, \alpha). Eigenvalues (trivially, diagonal matrix): 1,1,α1, 1, \alpha, so σ(A)={1,1,α}\sigma(A) = \{1, 1, \alpha\} with algebraic multiplicity of 11 equal to 22.

  • α>0\alpha > 0: all eigenvalues >0>0 \Rightarrow positive definite.
  • α=0\alpha = 0: eigenvalues {1,1,0}\{1, 1, 0\}, all 0\ge 0, not all >0>0 \Rightarrow positive semi-definite (but not positive definite). Confirm by exhibiting a non-zero vector on which ff vanishes: f(0,0,x3)=0f(0, 0, x_3) = 0.
  • α<0\alpha < 0: eigenvalues {1,1,α}\{1, 1, \alpha\}, mixed signs \Rightarrow indefinite.

Parts (b), (c), (d) follow the same pattern: compute AA from ff (diagonal entries = xi2x_i^2 coefficients; off-diagonal = half the xixjx_i x_j coefficient), then either compute eigenvalues or apply Sylvester–Jacobi.

Source: TA41_lin_alg.pdf pp.9–10 exercise 8.


Problem 6.8 — Invertibility inferred from signs (Theorem 1207 applied)

Let AA be a symmetric 4×44 \times 4 matrix that is positive definite. What can you say about detA1\det A^{-1} and the eigenvalues of A1A^{-1}?

Solution. By Theorem 1207 (§3.6), A1A^{-1} exists and is also positive definite. By Theorem 1210, its four eigenvalues are all strictly positive. By Theorem 1191, they are exactly the reciprocals of the eigenvalues of AA. Therefore:

  • detA1=1detA>0\det A^{-1} = \tfrac{1}{\det A} > 0.
  • All four eigenvalues of A1A^{-1} are strictly positive — in fact, they are 1/λi1/\lambda_i where λi\lambda_i are the eigenvalues of AA, each >0>0.

(This is the kind of short theoretical question that Sicconi's practice materials refer to as "theorem-chaining" problems; useful as warm-up.)


§7. Common Pitfalls

Below are the errors Sicconi flags most often and the traps the Bocconi exam-writers exploit. Read this list the day before the exam.

  1. Eigenvector vs eigenspace dimension.

    • An eigenvector is a single non-zero vector. An eigenspace is a subspace (possibly multi-dimensional) containing all eigenvectors for a given λ\lambda, plus the zero vector. Saying "the eigenspace is (1,0,1)T(1, 0, -1)^T" is wrong; it's the span of (1,0,1)T(1,0,-1)^T.
    • Eigenvalue multiplicity = dimension of the corresponding eigenspace (for symmetric matrices).
  2. Sylvester–Jacobi for negative definite — alternating signs, not all negative.

    • Pattern is (,+,,+,)(-, +, -, +, \ldots), not (,,,,)(-, -, -, -, \ldots). A matrix with all leading principal minors negative is NOT negative definite in general (only in the trivial 1×11 \times 1 case).
    • Mnemonic: negative definite     \iff A-A is positive definite; so det(A)k=(1)kdetAk>0\det(-A)_k = (-1)^k \det A_k > 0, forcing detAk\det A_k to alternate.
  3. "Positive semi-definite" ≠ "positive definite". (May 2024 Q7 trap.)

    • Positive definite: xTAx>0\mathbf{x}^T A \mathbf{x} > 0 for every x0\mathbf{x} \neq \mathbf{0}. Requires all eigenvalues strictly >0>0, or all leading principal minors strictly >0>0.
    • Positive semi-definite: xTAx0\mathbf{x}^T A \mathbf{x} \ge 0 for every x\mathbf{x}. Allows some eigenvalues =0=0.
    • If detA=0\det A = 0, AA is NOT positive definite — even if every other leading minor is positive. At best it's semi-definite.
  4. Sylvester is inconclusive when a leading minor is zero.

    • If detAk=0\det A_k = 0 for some kk, you cannot conclude definiteness or indefiniteness from leading minors alone.
    • Fallback: compute all principal minors (not just leading), or find the eigenvalues directly.
  5. Gram–Schmidt is only needed inside a single eigenspace.

    • Eigenvectors from different eigenspaces of a symmetric matrix are automatically orthogonal (Theorem 1184). Trying to re-orthogonalize them wastes time and introduces arithmetic errors.
    • Inside one eigenspace of multiplicity 2\ge 2, arbitrary basis vectors may NOT be orthogonal — apply Gram–Schmidt there.
  6. Normalise after, not before, Gram–Schmidt.

    • Order: orthogonalise \to normalise. Reversing can leave vectors off the unit sphere.
  7. Quadratic form coefficients: half, not whole.

    • When building AA from ff, the off-diagonal entries are half the coefficient of xixjx_i x_j, not the whole coefficient. aij+aji=a_{ij} + a_{ji} = coefficient of xixjx_i x_j; symmetry (aij=ajia_{ij} = a_{ji}) then gives aij=12a_{ij} = \tfrac{1}{2}\cdot(coef).
  8. Eigenvalues of AA are eigenvalues of the shifted matrix, with a shift.

    • If B=A+cIB = A + cI, then σ(B)={λ+c:λσ(A)}\sigma(B) = \{\lambda + c : \lambda \in \sigma(A)\}, with the SAME eigenvectors. Useful for quick checks.
  9. Spectrum shortcut traps.

    • λi=trA\sum \lambda_i = \operatorname{tr} A and λi=detA\prod \lambda_i = \det A always hold (with multiplicity!). On MCQs with multiplicities, use the sum and product to eliminate wrong options fast (Sicconi's favourite shortcut on Linear Algebra and Differential Calculus.pdf p.12).
    • But: knowing only the sum does NOT determine the eigenvalues; you might need to compute det\det or the full characteristic polynomial.
  10. Eigenbasis vs basis of eigenvectors.

    • May 2024 Q7(c) explicitly says "not necessarily orthonormal" — meaning a basis of eigenvectors suffices (no need for Gram–Schmidt + normalisation). Read the problem statement carefully.
    • "Eigenbasis" in Sicconi's notes means orthonormal basis of eigenvectors; "basis made of eigenvectors" is the weaker requirement.
  11. Forgetting that "symmetric" matters.

    • All the eigenvalue theorems here (1179, 1181, 1184, 1191, 1210, 1217, Spectral, 1207) assume AA symmetric. Non-symmetric matrices may have complex eigenvalues, and their eigenvectors from different eigenspaces need not be orthogonal.
  12. Confusing the kernel of a linear map with the image.

    • Kernel of AA = null space = {x:Ax=0}\{\mathbf{x} : A\mathbf{x} = \mathbf{0}\} = eigenspace W0W_0 (when 00 is an eigenvalue).
    • For a symmetric matrix, dimker(A)+rank(A)=n\dim \ker(A) + \operatorname{rank}(A) = n (rank–nullity).
  13. Misreading the exam MCQ format.

    • Wrong answer: 1-1 pt. Blank: 00 pt. Correct: +5+5 pt. On MCQs where you can eliminate 2\ge 2 options, guess between the remaining. On a pure 50-50, skip.

Final checklist before the exam (linear algebra slice)

  • [ ] I can state and prove Theorem 1179 (eigenvalue     \iff root of char. poly.).
  • [ ] I can state and prove Theorem 1181 (eigenspace is a subspace).
  • [ ] I can state and prove Theorem 1184 (orthogonality of eigenvectors from distinct eigenvalues).
  • [ ] I can state and prove Theorem 1191 parts (1) and (2).
  • [ ] I can state (proof not required for some questions) Proposition 1201 and derive the matrix-form bijection.
  • [ ] I can state and prove Theorem 1207 (positive definite \Rightarrow invertible and A1A^{-1} PD).
  • [ ] I can state the Spectral Theorem and sketch the construction of BB.
  • [ ] I know the Sylvester–Jacobi sign patterns: (+,+,+,)(+, +, +, \ldots) for PD, (,+,,)(-, +, -, \ldots) for ND, and when the test is inconclusive.
  • [ ] I can classify any quadratic form both by eigenvalue signs and by Sylvester–Jacobi, and I know when to pick each method.
  • [ ] I can run M-LA-1 through M-LA-6 on a 3×33 \times 3 symmetric matrix in under 15 minutes.

End of Linear Algebra study guide.