Essential Student Algebra: Volume Four, Linear Algebra
T. S. Blyth & E. F. Robertson · First-year university · Study material generated by AI from this edition. Not affiliated with the publisher.
Loading…
Chapter 6 · pp. 79–96
Chapter 6: Dual Spaces
This chapter adds a way to measure lengths and angles to a vector space: the inner product. With it we can build orthonormal bases, find closest vectors, and see that every linear functional is an inner product with a fixed vector. That link between a space and its dual leads to adjoint maps, orthogonal complements and unitary matrices.
Inner products and their rules
An inner product on a complex vector space V gives each pair x,y a complex number (x∣y). It must obey four rules: (x+x′∣y)=(x∣y)+(x′∣y), (αx∣y)=α(x∣y), (x∣y)=(y∣x), and (x∣x)>0 unless x=0. Conjugate symmetry makes (x∣x) real, so the positivity rule makes sense. A real inner product space is the same idea over R, without conjugation. (p. 79)
Conjugate symmetry also controls the second slot. The product is additive there, but scalars come out conjugated: (x∣αy)=(αy∣x)=αˉ(x∣y). Also (x∣0)=0=(0∣x). Over R the product is linear in both slots. (pp. 79–80)
A common mistake is to think a complex inner product is linear in both arguments. It is not: (x∣iy)=−i(x∣y). Forgetting this gives sign errors when expanding ∥x+y∥2 or computing Fourier coefficients. (p. 80)
Examples of inner product spaces
On Cn the standard product is (z∣w)=∑jzjwj. On Rn it is ∑jxjyj, which is the dot product when n=2,3. For z=(1,i) and w=(i,1): (z∣w)=1⋅iˉ+i⋅1=−i+i=0. (p. 80)
Two more examples: continuous real functions on [a,b] with (f∣g)=∫abfg, where the integral replaces the sum, and matrices with (A∣B)=tr(B∗A). For n×n complex matrices the trace equals ∑aijbij, so this is the standard product on Cn2. For real matrices B∗ is just Bt. (pp. 80–81)
Norm, distance and two key inequalities
The norm is ∥x∥=(x∣x) and the distance is d(x,y)=∥x−y∥. In R2, ∥x∥2=x12+x22, so these agree with Pythagoras. Also ∥x∥=0 only when x=0. (p. 81)
Scaling gives ∥λx∥=∣λ∣∥x∥, because (λx∣λx)=λλˉ(x∣x). For example ∥(3+4i)x∥=5∥x∥. (p. 81)
Cauchy-Schwarz:∣(x∣y)∣≤∥x∥∥y∥. For x=0 put z=y−∥x∥2(y∣x)x, which is orthogonal to x. Then 0≤∥z∥2=∥y∥2−∥x∥2∣(x∣y)∣2, and rearranging gives the result. Check with x=(1,2), y=(3,1): 5≤510≈7.07. (pp. 81–82)
Triangle inequality:∥x+y∥≤∥x∥+∥y∥. Expand ∥x+y∥2=∥x∥2+2Re(x∣y)+∥y∥2, bound Re(x∣y) by ∣(x∣y)∣, apply Cauchy-Schwarz to get (∥x∥+∥y∥)2, then take square roots. (p. 82)
Orthogonal and orthonormal sets
Vectors are orthogonal when (x∣y)=0. An orthogonal set has all distinct pairs orthogonal, and an orthonormal set is an orthogonal set of unit vectors. In R2, orthogonal means x1y1+x2y2=0, i.e. perpendicular lines. The standard bases of Rn and Cn are orthonormal, and so are the functions 1,sinkx,coskx as an orthogonal set under ∫−ππfg (they need scaling to be unit length). (pp. 82–83)
To make an orthogonal set of nonzero vectors orthonormal, replace each x by x∗=x/∥x∥. The length becomes 1 and the direction is unchanged, so orthogonality stays. Example: (3,4)↦(3/5,4/5). (p. 83)
Orthogonal sets of nonzero vectors are linearly independent. If ∑λkxk=0, take the product with xi. Every term vanishes except λi∥xi∥2=0, so λi=0. Hence n orthogonal nonzero vectors in an n-dimensional space form a basis. (p. 83)
Bessel's inequality and expansions
If {e1,…,en} is orthonormal, then ∑k∣(x∣ek)∣2≤∥x∥2 (Bessel). Proof: let z=x−∑(x∣ek)ek; then 0≤∥z∥2=∥x∥2−∑∣(x∣ek)∣2. Equality holds exactly when z=0. (pp. 83–84)
When the ek span W, these statements are equivalent: x∈W; Bessel is an equality; x=∑(x∣ek)ek; and (x∣y)=∑(x∣ek)(ek∣y) for all y. The key step is that if x=∑λkek, taking the product with ej gives λj=(x∣ej). This gives a computable test for membership in the span. (pp. 83–84)
An orthonormal basis is an orthonormal set that is also a basis. Examples: the standard bases, and the matrices Epq (a single 1 in position (p,q)) under (A∣B)=tr(B∗A). (p. 84)
Relative to an orthonormal basis, three formulas hold: x=∑(x∣ek)ek (Fourier expansion), ∥x∥2=∑∣(x∣ek)∣2, and (x∣y)=∑(x∣ek)(y∣ek) (Parseval). The numbers (x∣ek) are Fourier coefficients. For x=(1,1,1) and the basis found below, they are 2, 2/3, 1/3, and 2+32+31=3=∥x∥2. (pp. 86–87)
The Gram-Schmidt process
Given independent vectors x1,…,xk, Gram-Schmidt produces an orthonormal list y1,…,yk with the same span. Set y1=x1∗. For each later t, subtract the components along earlier vectors, xt−∑i<t(xt∣yi)yi, then divide by the norm. The subtracted vector is orthogonal to the earlier yj by construction, and each yt is a combination of x1,…,xt. (p. 85)
Example in R3 with (0,1,1),(1,0,1),(1,1,0): y1=(0,1,1)/2. Then (1,0,1)−(0,21,21)=21(2,−1,1), giving y2=(2,−1,1)/6. Finally the third vector minus its components is 32(1,1,−1), giving y3=(1,1,−1)/3. Dropping scalar factors before normalising saves arithmetic. (p. 86)
Consequently every finite-dimensional inner product space has an orthonormal basis. Any orthonormal subset extends to one: extend it to an ordinary basis first, then apply Gram-Schmidt. The original vectors come out unchanged. (pp. 86–87)
Inner product isomorphisms
An inner product isomorphism is a vector space isomorphism f with (f(x)∣f(y))=(x∣y). It preserves lengths, distances and orthogonality; a rotation of R2 is an example. (p. 87)
A linear map is such an isomorphism exactly when it sends an orthonormal basis of V to an orthonormal basis of W. One direction: (f(ei)∣f(ej))=(ei∣ej)=δij. The other uses Parseval's identity in both spaces. (pp. 87–88)
Duality and adjoints
A conjugate transformation satisfies f(x+y)=f(x)+f(y) and f(λx)=λˉf(x). If it is also bijective it is a conjugate isomorphism. Over R these are just linear maps. An example is x↦xˉ on Cn. (p. 88)
In a finite-dimensional inner product space V, the map y↦y♭, where y♭(x)=(x∣y), is a conjugate isomorphism onto the dual space. It is injective because y♭=0 gives (y∣y)=0. It is surjective because any functional f equals x↦(x∣z) with z=∑f(ek)ek. So every linear functional is an inner product with a unique vector. (pp. 89–90)
This gives the adjoint. For linear f:V→W there is a unique linear f∗:W→V with (f(x)∣y)=(x∣f∗(y)). For x↦Ax on Cn it is y↦A∗y. (p. 90)
Rules: (f+g)∗=f∗+g∗, (λf)∗=λˉf∗, (h∘f)∗=f∗∘h∗, (f∗)∗=f. Each follows by checking the defining identity and using uniqueness. The scalar is conjugated and the order reverses. (p. 91)
When dimV=dimW, these are equivalent: f is an inner product isomorphism; f is an isomorphism with f−1=f∗; f∘f∗=idW; f∗∘f=idV. If f∗f=id then (f(x)∣f(y))=(x∣f∗f(y))=(x∣y). So the inverse is found by taking the adjoint. (pp. 91–92)
Orthogonal complements and best approximation
The orthogonal complement of a non-empty set E is E⊥={y:(x∣y)=0 for all x∈E}, always a subspace. In R3, the complement of a plane through the origin is the perpendicular line through the origin. (p. 92)
If W is a finite-dimensional subspace, V=W⊕W⊥. For x∈V take an orthonormal basis ei of W, let x′=∑(x∣ei)ei and x′′=x−x′. Then x′′∈W⊥. The sum is direct because a vector z in both has (z∣z)=0. (p. 92)
The closest vector in W to x is its component a=∑(x∣ek)ek. Write x=a+b with b∈W⊥. For y∈W, ∥x−y∥2=∥a−y∥2+∥b∥2≥∥x−a∥2, with equality only at y=a. This is orthogonal projection. (p. 93)
Application: with (f∣g)=∫−ππfg, the functions 1,sinkx,coskx (k≤n) are orthogonal and span a (2n+1)-dimensional space. The nearest trigonometric polynomial to f is 21a0+∑k≤n(akcoskx+bksinkx), where ak=π1∫fcoskxdx and bk=π1∫fsinkxdx. As n→∞ this gives the Fourier series, which converges to f when f is infinitely differentiable. (pp. 93–94)
In finite dimensions, dimW⊥=dimV−dimW and W⊥⊥=W. Clearly W⊆W⊥⊥, and both have the same dimension. (p. 94)
For subspaces A,B: A⊆B implies B⊥⊆A⊥; (A∩B)⊥=A⊥+B⊥; (A+B)⊥=A⊥∩B⊥. Taking complements reverses inclusion and swaps sums with intersections. (p. 94)
For linear f: Imf∗=(Kerf)⊥ and Kerf∗=(Imf)⊥. If x=f∗(y) and z∈Kerf, then (z∣x)=(f(z)∣y)=0. If f∗(y)=0, then (f(x)∣y)=0 for all x. A dimension count turns the inclusions into equalities. (pp. 94–95)
Matrices of adjoints and unitary matrices
The adjoint (conjugate transpose) of A=[aij] is A∗ with [A∗]ij=aji. For A=(12i3), A∗=(1−i23). For real matrices it is the transpose. (p. 95)
If A represents f relative to ordered orthonormal bases, then A∗ represents f∗ relative to the same bases. This works because aij=(f(vj)∣wi)=(f∗(wi)∣vj). Orthonormal bases are essential; otherwise the matrix of f∗ is generally not A∗. (pp. 95–96)
A square matrix U is unitary if U−1=U∗ (orthogonal if real). B is unitarily similar to A if B=U∗AU for some unitary U. Unitary matrices represent inner product isomorphisms. Two matrices represent the same map in different orthonormal bases exactly when they are unitarily similar, which is an equivalence relation. (p. 96)
Check yourself
Show that (x∣αy)=αˉ(x∣y) follows from the axioms, and compute (x∣iy) if (x∣y)=2+i.
Apply Gram-Schmidt to (1,1,0) and (1,0,1) in R3 with the standard product.
Why must an orthogonal set of nonzero vectors be linearly independent?
If W is a plane through the origin in R3, what are W⊥ and W⊥⊥, and how are their dimensions related?
Find A∗ for A=(201+i−i), and explain why orthonormal bases matter for the matrix of f∗.