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×VF\langle \bullet, \bullet \rangle: V \times V \to F that is linear in both of its parameters. In other words, for all x,y,zVx, 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=xAy\langle x, y \rangle = x^{\top}Ay for some AFn×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,jxiyjbi,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 AFn×nA \in F^{n \times n} via Ai,j=bi,bjA_{i, j} = \langle b_i, b_j \rangle, the above reads x,y=xAy\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×VF\langle \bullet, \bullet \rangle: V \times V \to F.

Say BFn×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=BMBM' = B^{\top}MB.

Proof: The entries Mi,jM'_{i, j} are determined by the relationship Mi,j=bi,bjM_{i, j}' = \langle b_i, b_j \rangle, since {b1,,bn}\{b_1, \dots, b_n\} is the basis for MM'.

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(BMB)ej= (BMB)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 MM' and BMBB^{\top}MB match everywhere, meaning M=BMBM' = 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=BMBM' = B^{\top}MB is a claim about matrices, NOT linear transformations!

In particular, the expression BMBB^{\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)=xAx+bx+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 202 \neq 0. If F=F2F = \mathbb{F}_2, for example, then there is no symmetric AA that yields xAx=x1x2x^{\top}Ax = x_1x_2.

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

f(x)12x ⁣[2fxixj(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×VC\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,zVx, 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 xx^{\dagger} and xx^* 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 BCn×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=BMBM' = 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:VVA: V \to V does x,y=Ax,Ay\langle x, y \rangle = \langle Ax, Ay \rangle for all x,yVx, y \in V always hold?

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

Remark. This implies A:VVA: V \to V preserves the dot product x,y=xy\langle x, y \rangle = x^{\top}y if and only if AA=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=xMy\langle x, y \rangle = x^{\dagger}My in this basis. The answer is all AA for which AMA=MA^{\dagger}MA = M.

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

Definition (Unitary). Take the Hermitian form to be x,y=xy\langle x, y \rangle = x^{\dagger}y. A matrix ACn×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,yCnx, y \in \mathbb{C}^n. Equivalently, AA is unitary if AA=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 ACn×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,iA_{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 ACn×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=xy\langle x, y \rangle = x^{\dagger}y is the basic inner product.

Proof: Just compute: Av,v=(Av)v=vAv=v,Av=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=xAy\langle x, y \rangle = x^{\dagger}Ay. Then v,v=v,v\langle v, v \rangle = \overline{\langle v, v \rangle}, so v,vR\langle v, v \rangle \in \mathbb{R}.

This reads vAvRv^{\dagger}Av \in \mathbb{R}. But vAv=λ(vv)v^{\dagger}Av = \lambda(v^{\dagger}v), and vvRv^{\dagger}v \in \mathbb{R} by inspection (i.e. it's ivivi=ivi2R\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 RC\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}”.