MIT 18.701 — Lecture 16

Bilinear and Hermitian Forms

We now shift our focus to generalizations of the dot product known as bilinear forms.

§ Bilinear Forms: Definitions and Properties

Definition (Bilinear Form). Let VV be an FF-vector space. Then a bilinear form is a function ⟨∙,∙⟩:V×V→F\langle \bullet, \bullet \rangle: V \times V \to F that is linear in both of its parameters. In other words, for all x,y,z∈Vx, y, z \in V and λ∈F\lambda \in F, we have:

⟨λx,y⟩=λ⟨x,y⟩=⟨x,λy⟩  and  ⟨x+y,z⟩=⟨x,z⟩+⟨y,z⟩  and  ⟨x,y+z⟩=⟨x,y⟩+⟨x,z⟩\langle \lambda x, y \rangle = \lambda \langle x, y \rangle = \langle x, \lambda y \rangle ~ \text{ and } ~ \langle x + y, z \rangle = \langle x, z \rangle + \langle y, z \rangle ~ \text{ and } ~ \langle x, y + z \rangle = \langle x, y \rangle + \langle x, z \rangle

Example. Some examples and non-examples:

We can also use adjectives to describe bilinear forms.

Definition (Symmetric, Semidefinite, Definite). For a bilinear form ⟨∙,∙⟩\langle \bullet, \bullet \rangle over an R\mathbb{R}-vector space VV,

Definition (Inner Product). A bilinear form that is symmetric and positive definite is an inner product.

Example. Consider the following candidates for bilinear forms over R2\mathbb{R}^2.

§ Bilinear Form   ⟹  \implies Matrix (Only if You Choose a Basis!)

Theorem. (Bilinear Form Classification) Suppose VV is an FF-vector space with basis (b1,…,bn)(b_1, \dots, b_n). Then every bilinear form ⟨∙,∙⟩\langle \bullet, \bullet \rangle over VV must look like ⟨x,y⟩=x⊤Ay\langle x, y \rangle = x^{\top}Ay for some A∈Fn×nA \in F^{n \times n}.

Proof: Consider x=∑ixibix = \sum_{i} x_i b_i and y=∑jyjbjy = \sum_j y_j b_j. Then by the linearity of bilinear forms,

⟨x,y⟩=⟨∑ixibi, ∑jyjbj⟩=∑i,jxiyj⟨bi,bj⟩.\langle x, y \rangle = \left \langle \sum_i x_i b_i, \ \sum_j y_j b_j \right \rangle = \sum_{i, j} x_iy_j \langle b_i, b_j \rangle.

Then if we take A∈Fn×nA \in F^{n \times n} via Ai,j=⟨bi,bj⟩A_{i, j} = \langle b_i, b_j \rangle, the above reads ⟨x,y⟩=x⊤Ay\langle x, y \rangle = x^{\top}Ay, as desired.   ■\blacksquare

Notice the analogy with “Linear Transformations   ⟹  \implies Matrix (Only if You Choose a Basis!)”. To be explicit,

Theorem. (Change of Bilinear Basis) Consider a bilinear form ⟨∙,∙⟩:V×V→F\langle \bullet, \bullet \rangle: V \times V \to F.

Say B∈Fn×nB \in F^{n \times n} is the matrix with columns {b1,…,bn}\{b_1, \dots, b_n\}, meaning Bei=biBe_i = b_i for all ii. Then M′=B⊤MBM' = B^{\top}MB.

Proof: The entries Mi,j′M'_{i, j} are determined by the relationship Mi,j′=⟨bi,bj⟩M_{i, j}' = \langle b_i, b_j \rangle, since {b1,…,bn}\{b_1, \dots, b_n\} is the basis for M′M'.

Using the identity bi=Beib_i = Be_i, this relationship rewrites itself in the basis {e1,…,en}\{e_1, \dots, e_n\} as:

Mi,j′=⟨Bei,Bej⟩= (Bei)⊤M(Bej)  in the basis {e1,…,en}= ei⊤(B⊤MB)ej= (B⊤MB)i,j.\begin{align*}M'_{i, j} = \langle Be_i, Be_j \rangle = \ & (Be_i)^{\top}M (Be_j) ~ \text{ in the basis } \{e_1, \dots, e_n\} \\ = \ & e_i^{\top}(B^{\top}MB)e_j \\ = \ & (B^{\top}MB)_{i, j}.\end{align*}

So the entries of M′M' and B⊤MBB^{\top}MB match everywhere, meaning M′=B⊤MBM' = B^{\top}MB.   ■\blacksquare

— Warning: Matrices are Grids of Numbers —

Given a choice of basis, every bilinear form corresponds to a matrix.

This matrix is just a matrix. It is nothing more than a grid of numbers in Fn×nF^{n \times n}.

Do not try to interpret this matrix as a linear transformation—it won't work.

The statement M′=B⊤MBM' = B^{\top}MB is a claim about matrices, NOT linear transformations!

In particular, the expression B⊤MBB^{\top}MB means nothing irrespective of a basis,
because there is no meaning to “the transpose of a linear map” without a basis.

More explicitly, in the proof above, the line “⟨Bei,Bej⟩=(Bei)⊤M(Bej)\langle Be_i, Be_j \rangle = (Be_i)^{\top}M(Be_j)”
is only true if we define the transpose (∙)⊤(\bullet)^{\top} using the basis {e1,…,en}\{e_1, \dots, e_n\}.
The truth value of this line depends on our choice of basis,
because (Bei)⊤M(Bej)(Be_i)^{\top}M(Be_j) has a basis-dependent value.

§ Foreshadowing: Quadratic Polynomials are Bilinear Forms

Claim: Every quadratic polynomial in nn variables (x1,…,xn)(x_1, \dots, x_n) can be viewed in terms of bilinear forms.

Proof: Just note the injection:

quadratic polynomials=(∑i,jAi,jxixj)⏟degree 2 terms+(∑ibixi)⏟degree 1 terms + (c)=x⊤Ax+b⊤x+c.\text{quadratic polynomials} = \underbrace{\left(\sum_{i, j} A_{i, j} x_i x_j\right)}_{\text{degree 2 terms}} + \underbrace{\left(\sum_{i} b_i x_i \right)}_{\text{degree 1 terms}} \ + \ \left(c\right) = x^{\top}Ax + b^{\top}x + c.

So quadratic polynomials inject with ordered triples (A,b,c)∈Rn×n×Rn×R(A, b, c) \in \mathbb{R}^{n \times n} \times \mathbb{R}^n \times \mathbb{R}. Specifically, both Ai,jA_{i, j} and Aj,iA_{j, i} correspond to the coefficient of xixjx_ix_j, so the matrix AA better be symmetric.

Remark. Only if we're working in a field FF where 2≠02 \neq 0. If F=F2F = \mathbb{F}_2, for example, then there is no symmetric AA that yields x⊤Ax=x1x2x^{\top}Ax = x_1x_2.

This injection comes naturally with Taylor Series; the 2nd2^{\text{nd}} degree Taylor polynomial of f:Rn→Rf: \mathbb{R}^n \to \mathbb{R} at x=0x = 0 is:

f(x)≈12x⊤ ⁣[∂2f∂xi∂xj(0)]⏟Ax+[∇f(0)]⏟b ⁣⊤x+f(0)⏟c.f(x) \approx \dfrac{1}{2} x^{\top} \!\underbrace{\left[\dfrac{\partial^2 f}{\partial x_i \partial x_j}(0)\right]}_{A} x + \underbrace{\left[\nabla f(0)\right]}_{b} \!^{\top}x + \underbrace{f(0)}_{c}.

The funny double-partial-derivative AA is called the Hessian. The Hessian is always symmetric, expectedly.

§ Hermitian Forms: Doing it all again…

What if bilinear forms were not actually bilinear and were forced to live in C\mathbb{C}?

Definition (Hermitian Form). Let VV be a C\mathbb{C}-vector space. Then a Hermitian form is a function ⟨∙,∙⟩:V×V→C\langle \bullet, \bullet \rangle : V \times V \to \mathbb{C} that is “linear” and “symmetric” in both of its parameters. In other words, for all x,y,z∈Vx, y, z \in V and μ,λ∈C\mu, \lambda \in \mathbb{C}, we have:

⟨x+y,z⟩=⟨x,z⟩+⟨y,z⟩  and  ⟨λx,μz⟩=λ‾μ⟨x,z⟩  and  ⟨y,x⟩=⟨x,y⟩‾.\langle x + y, z \rangle = \langle x, z \rangle + \langle y, z \rangle ~ \text{ and } ~ \langle \lambda x, \mu z \rangle = \overline{\lambda} \mu \langle x, z \rangle ~ \text{ and } ~ \langle y, x \rangle = \overline{\langle x, y \rangle} .

Remark. Hermitian forms are not bilinear forms. They're a completely separate thing.

Remark. The properties “positive semidefinite” and “positive definite” carry over to Hermitian forms. However, the “symmetric” condition is already forcefully embedded within the definition.

Definition (Dagger†\textbf{Dagger}^\dagger and Asterisk∗\textbf{Asterisk}^*). Let x†x^{\dagger} and x∗x^* both be shorthands for the conjugate transpose x‾⊤\overline{x}^{\top}.

Our discussion of “Bilinear form   ⟹  \implies Matrix (Only if You Choose a Basis!)” carries over. In particular,

Theorem. (Hermitian Forms w/ Bases) Suppose VV is a C\mathbb{C}-vector space with Hermitian form ⟨∙,∙⟩\langle \bullet, \bullet \rangle.

Say B∈Cn×nB \in \mathbb{C}^{n \times n} is the matrix with columns {b1,…,bn}\{b_1, \dots, b_n\}, meaning Bei=biBe_i = b_i for all ii. Then M′=B†MBM' = B^{\dagger}MB.

§ Transformations that Preserve a Bilinear / Hermitian Form

Which linear transformations preserve a bilinear form? Phrased more precisely,

Question: For which linear transformations A:V→VA: V \to V does ⟨x,y⟩=⟨Ax,Ay⟩\langle x, y \rangle = \langle Ax, Ay \rangle for all x,y∈Vx, y \in V always hold?

Answer: Work in a fixed basis, and say ⟨x,y⟩=x⊤My\langle x, y \rangle = x^{\top}My in this basis. The answer is all AA for which A⊤MA=MA^{\top}MA = M. The proof follows immediately from the change-of-basis formula.

Remark. This implies A:V→VA: V \to V preserves the dot product ⟨x,y⟩=x⊤y\langle x, y \rangle = x^{\top}y if and only if A⊤A=InA^{\top}A = I_n, as expected.
A subgroup of matrices that preserve a more general symmetric bilinear form is called an indefinite orthogonal group.

We can repeat the same thing for Hermitian forms, too.

Answer: Work in a fixed basis, and say ⟨x,y⟩=x†My\langle x, y \rangle = x^{\dagger}My in this basis. The answer is all AA for which A†MA=MA^{\dagger}MA = M.

Orthogonal matrices preserve bilinear forms ⟨x,y⟩=x⊤y\langle x, y \rangle = x^{\top}y. Analogously…

Definition (Unitary). Take the Hermitian form to be ⟨x,y⟩=x†y\langle x, y \rangle = x^{\dagger}y. A matrix A∈Cn×nA \in \mathbb{C}^{n \times n} is unitary if ⟨x,y⟩=⟨Ax,Ay⟩\langle x, y \rangle = \langle Ax, Ay \rangle for all x,y∈Cnx, y \in \mathbb{C}^n. Equivalently, AA is unitary if A†A=IA^{\dagger}A = I.

§ Symmetric Bilinear / Hermitian Forms

We'll conclude lecture by discussing what happens when bilinear forms are symmetric.

Definition (Hermitian Matrix). A matrix A∈Cn×nA \in \mathbb{C}^{n \times n} is Hermitian if A†=AA^{\dagger} = A.

Theorem. (Symmetry Begets Symmetry) Consider a symmetric bilinear / Hermitian form ⟨∙,∙⟩\langle \bullet, \bullet \rangle. Pick a basis.

Proof: For the bilinear form, the symmetry implies ⟨ei,ej⟩=⟨ej,ei⟩\langle e_i, e_j \rangle = \langle e_j, e_i\rangle, which reads Ai,j=Aj,iA_{i, j} = A_{j, i}. So A⊤=AA^{\top} = A.

For the Hermitian form, the symmetry implies ⟨ei,ej⟩=⟨ej,ei⟩‾\langle e_i, e_j \rangle = \overline{\langle e_j, e_i \rangle}, which reads Ai,j=Aj,i‾A_{i, j} = \overline{A_{j, i}}. So A†=AA^{\dagger} = A.   ■\blacksquare

We now conclude with the following fun fact:

Theorem. (Eigenvalues are Real) If A∈Cn×nA \in \mathbb{C}^{n \times n} is Hermitian, then all eigenvalues of AA are real.

Proof: Suppose vv has eigenvalue λ\lambda.

Lemma. Since AA is Hermitian, we have ⟨Av,v⟩=⟨v,Av⟩\langle Av, v \rangle = \langle v, Av \rangle, where ⟨x,y⟩=x†y\langle x, y \rangle = x^{\dagger}y is the basic inner product.

Proof: Just compute: ⟨Av,v⟩=(Av)†v=v†A†v=⟨v,A†v⟩=⟨v,Av⟩\langle Av, v \rangle = (Av)^{\dagger}v = v^{\dagger}A^{\dagger}v = \langle v, A^{\dagger}v \rangle = \langle v, Av \rangle.   □\square

The lemma says ⟨v,Av⟩=⟨Av,v⟩\langle v, Av \rangle = \langle Av, v \rangle, which reads ⟨v,λv⟩=⟨λv,v⟩\langle v, \lambda v \rangle = \langle \lambda v, v \rangle. Thus λ⟨v,v⟩=λ‾⟨v,v⟩\lambda \langle v, v \rangle = \overline{\lambda} \langle v, v \rangle, so λ=λ‾\lambda = \overline{\lambda}, done.   ■\blacksquare

For fun, here's another proof of the same theorem.

Proof: Suppose vv has eigenvalue λ\lambda, and consider the Hermitian form ⟨x,y⟩=x†Ay\langle x, y \rangle = x^{\dagger}Ay. Then ⟨v,v⟩=⟨v,v⟩‾\langle v, v \rangle = \overline{\langle v, v \rangle}, so ⟨v,v⟩∈R\langle v, v \rangle \in \mathbb{R}.

This reads v†Av∈Rv^{\dagger}Av \in \mathbb{R}. But v†Av=λ(v†v)v^{\dagger}Av = \lambda(v^{\dagger}v), and v†v∈Rv^{\dagger}v \in \mathbb{R} by inspection (i.e. it's ∑ivi‾vi=∑i∣vi∣2∈R\sum_i \overline{v_i}v_i = \sum_i |v_i|^2 \in \mathbb{R}). So λ∈R\lambda \in \mathbb{R}.   ■\blacksquare

By corollary, every symmetric real matrix has real eigenvalues! This is because R⊆C\mathbb{R} \subseteq \mathbb{C}. The intuition is that the complex-valued Hermitian setting will yield results in the real setting.

Remark. The analogue of “Hermitian over C\mathbb{C}” is “symmetric over R\mathbb{R}”.