maths.free › Linear Algebra › 2. Vectors › Orthogonality
Orthogonality
In this section we define a couple more operations with vectors, and prove a few theorems.
Orthogonality
In this section we define a couple more operations with vectors, and prove a few theorems. At first blush these definitions and results will not appear central to what follows, but we will make use of them at key points in the remainder of the course (such as , ). Because we have chosen to use $\complexes$ as our set of scalars, this subsection is a bit more, uh, complex than it would be for the real numbers. We will explain as we go along how things get easier for the real numbers ${\mathbb R}$. If you have not already, now would be a good time to review some of the basic properties of arithmetic with complex numbers described in . With that done, we can extend the basics of complex number arithmetic to our study of vectors in $\complex{m}$.
- \set{\colvector{1\\-1\\2},\,\colvector{5\\3\\-1},\,\colvector{8\\4\\-2}}
- What is the distinction between an orthogonal set and an orthonormal set?
- What is nice about the output of the Gram-Schmidt process?
Complex Arithmetic and Vectors
We know how the addition and multiplication of complex numbers is employed in defining the operations for vectors in $\complex{m}$ ( and ). We can also extend the idea of the conjugate to vectors.
With this definition we can show that the conjugate of a column vector behaves as we would expect with regard to vector addition and scalar multiplication.
These two theorems together tell us how we can push complex conjugation through linear combinations.
Inner products
This operation is a bit different in that we begin with two vectors but produce a scalar. Computing one is straightforward.
Example: Computing some inner products
The inner product of \vect{u}=\colvector{2+3i\\5+2i\\-3+i}&&\text{and}&& \vect{v}=\colvector{1+2i\\-4+5i\\0+5i} is \innerproduct{\vect{u}}{\vect{v}} &=(\conjugate{2+3i})(1+2i)+(\conjugate{5+2i})(-4+5i)+(\conjugate{-3+i})(0+5i)\\ &=(2-3i)(1+2i)+(5-2i)(-4+5i)+(-3-i)(0+5i)\\ &=(8+i)+(-10+33i)+(5-15i)\\ &=3+19i
The inner product of \vect{w}=\colvector{2\\4\\-3\\2\\8}&&\text{and}&& \vect{x}=\colvector{3\\1\\0\\-1\\-2} is \innerproduct{\vect{w}}{\vect{x}}&= (\conjugate{2})3+(\conjugate{4})1+(\conjugate{-3})0+(\conjugate{2})(-1)+(\conjugate{8})(-2)\\ &=2(3)+4(1)+(-3)0+2(-1)+8(-2)=-8.
In the case where the entries of our vectors are all real numbers (as in the second part of ), the computation of the inner product may look familiar and be known to you as a dot product or scalar product. So you can view the inner product as a generalization of the scalar product to vectors from $\complex{m}$ (rather than ${\mathbb R}^m$).
Note that we have chosen to conjugate the entries of the first vector listed in the inner product, while it is almost equally feasible to conjugate entries from the second vector instead. In particular, prior to Version 2.90, we did use the latter definition, and this has now changed to the former, with resulting adjustments propogated up through (only). However, conjugating the first vector leads to much nicer formulas for certain matrix decompositions and also shortens some proofs.
There are several quick theorems we can now prove, and they will each be useful later.
Condensed — the full section is in Beezer, A First Course in Linear Algebra.
Norm
If treating linear algebra in a more geometric fashion, the length of a vector occurs naturally, and is what you would expect from its name. With complex numbers, we will define a similar function. Recall that if $c$ is a complex number, then $\modulus{c}$ denotes its modulus ().
Computing a norm is also easy to do.
Example: Computing the norm of some vectors
The norm of \vect{u}=\colvector{3+2i\\1-6i\\2+4i\\2+i} is \norm{\vect{u}}&= \sqrt{\modulus{3+2i}^2+\modulus{1-6i}^2+\modulus{2+4i}^2+\modulus{2+i}^2}\\ &=\sqrt{13+37+20+5}=\sqrt{75}=5\sqrt{3}
The norm of \vect{v}=\colvector{3\\-1\\2\\4\\-3} is \norm{\vect{v}}= \sqrt{\modulus{3}^2+\modulus{-1}^2+\modulus{2}^2+\modulus{4}^2+\modulus{-3}^2} =\sqrt{3^2+1^2+2^2+4^2+3^2}=\sqrt{39}.
Notice how the norm of a vector with real number entries is just the length of the vector. Inner products and norms are related by the following theorem.
When our vectors have entries only from the real numbers says that the dot product of a vector with itself is equal to the length of the vector squared.
Notice that contains three implications: \vect{u}\in\complex{m}&\Rightarrow\innerproduct{\vect{u}}{\vect{u}}\geq 0\\ \vect{u}=\zerovector&\Rightarrow\innerproduct{\vect{u}}{\vect{u}}=0\\ \innerproduct{\vect{u}}{\vect{u}}=0&\Rightarrow\vect{u}=\zerovector
The results contained in are summarized by saying the inner product is positive definite.
Condensed — the full section is in Beezer, A First Course in Linear Algebra.
Orthogonal Vectors
Orthogonal is a generalization of perpendicular. You may have used mutually perpendicular vectors in a physics class, or you may recall from a calculus class that perpendicular vectors have a zero dot product. We will now extend these ideas into the realm of higher dimensions and complex scalars.
Example: Two orthogonal vectors
The vectors \vect{u}&=\colvector{2 + 3i\\4 - 2i\\1 + i\\1 + i} & \vect{v}&=\colvector{1 - i\\2 + 3i\\4 - 6i\\1} are orthogonal since \innerproduct{\vect{u}}{\vect{v}} &=(2-3i)(1-i)+(4+2i)(2+3i)+(1-i)(4-6i)+(1-i)(1)\\ &=(-1-5i)+(2+16i)+(-2-10i)+(1-i)\\ &=0+0i.
We extend this definition to whole sets by requiring vectors to be pairwise orthogonal. Despite using the same word, careful thought about what objects you are using will eliminate any source of confusion.
We now define the prototypical orthogonal set, which we will reference repeatedly.
Notice that $\vect{e}_j$ is identical to column $j$ of the $m\times m$ identity matrix $I_m$ () and is a pivot column for $I_m$, since the identity matrix is in reduced row-echelon form. These observations will often be useful. We will reserve the notation $\vect{e}_i$ for these vectors. It is not hard to see that the set of standard unit vectors is an orthogonal set.
Example: Standard Unit Vectors are an Orthogonal Set
Compute the inner product of two distinct vectors from the set of standard unit vectors (), say $\vect{e}_i$, $\vect{e}_j$, where $i\neq j$, \innerproduct{\vect{e}_i}{\vect{e}_j}&= \conjugate{0}0+ \conjugate{0}0+\cdots+ \conjugate{1}0+\cdots+ \conjugate{0}0+\cdots+ \conjugate{0}1+\cdots+ \conjugate{0}0+ \conjugate{0}0\\ &=0(0)+0(0)+\cdots+1(0)+\cdots+0(1)+\cdots+0(0)+0(0)\\ &=0
So the set $\set{\vectorlist{e}{m}}$ is an orthogonal set.
Condensed — the full section is in Beezer, A First Course in Linear Algebra.
Gram-Schmidt Procedure
The Gram-Schmidt Procedure is really a theorem. It says that if we begin with a linearly independent set of $p$ vectors, $S$, then we can do a number of calculations with these vectors and produce an orthogonal set of $p$ vectors, $T$, so that $\spn{S}=\spn{T}$. Given the large number of computations involved, it is indeed a procedure to do all the necessary computations, and it is best employed on a computer. However, it also has value in proofs where we may on occasion wish to replace a linearly independent set by an orthogonal set.
This is our first occasion to use the technique of mathematical induction for a proof, a technique we will see again several times, especially in . So study the simple example described in first.
Example: Gram-Schmidt of three vectors
We will illustrate the Gram-Schmidt process with three vectors. Begin with the linearly independent (check this!) set S=\set{\vect{v}_1,\,\vect{v}_2,\,\vect{v}_3}=\set{ \colvector{1\\1+i\\1},\, \colvector{-i\\1\\1+i},\, \colvector{0\\i\\i} }
Then \vect{u}_1&=\vect{v_1}=\colvector{1\\1+i\\1}\\ \vect{u}_2&=\vect{v}_2 -\frac{\innerproduct{\vect{u}_1}{\vect{v}_2}}{\innerproduct{\vect{u}_1}{\vect{u}_1}}\vect{u}_1 =\frac{1}{4}\colvector{-2-3i\\1-i\\2+5i}\\ \vect{u}_3&=\vect{v}_3 -\frac{\innerproduct{\vect{u}_1}{\vect{v}_3}}{\innerproduct{\vect{u}_1}{\vect{u}_1}}\vect{u}_1 -\frac{\innerproduct{\vect{u}_2}{\vect{v}_3}}{\innerproduct{\vect{u}_2}{\vect{u}_2}}\vect{u}_2 =\frac{1}{11}\colvector{-3-i\\1+3i\\-1-i} and T=\set{\vect{u}_1,\,\vect{u}_2,\,\vect{u}_3} =\set{ \colvector{1\\1+i\\1},\, \frac{1}{4}\colvector{-2-3i\\1-i\\2+5i},\, \frac{1}{11}\colvector{-3-i\\1+3i\\-1-i} } is an orthogonal set (which you can check) of nonzero vectors and $\spn{T}=\spn{S}$ (all by ). Of course, as a by-product of orthogonality, the set $T$ is also linearly independent ().
One final definition related to orthogonal vectors.
Once you have an orthogonal set, it is easy to convert it to an orthonormal set multiply each vector by the reciprocal of its norm, and the resulting vector will have norm 1. This scaling of each vector will not affect the orthogonality properties (apply ).
We will see orthonormal sets again in . They are intimately related to unitary matrices () through . Some of the utility of orthonormal sets is captured by in . Orthonormal sets appear once again in where they are key in orthonormal diagonalization.
Condensed — the full section is in Beezer, A First Course in Linear Algebra.
Symbols used here
Add a_k for k = 1 up to n.
The non-negative number whose square (n-th power) is x.
Logical connectives.
i² = −1.
Inequalities that allow equality; < and > exclude it.
The two sides are different.
A quantity with magnitude and direction; a column of numbers.
A rectangular array of numbers; a linear map.
Scaling factor of area/volume under A; zero means singular.
The matrix that undoes A; A with rows and columns swapped.
The factor by which an eigenvector is stretched: Av = λv.
Σ u_i v_i; the length of v, √(v·v).
Questions people ask
What does a determinant mean geometrically?
It is the factor by which the matrix scales area (2×2) or volume (3×3), with a negative sign if orientation flips. Zero means the matrix flattens space and cannot be undone.
What is an eigenvector?
A direction the matrix does not turn — it only stretches it by the eigenvalue. Along eigenvectors a complicated matrix acts like multiplication by a number.
Why is matrix multiplication not commutative?
Because a matrix is a transformation and AB means "do B, then A". Rotating then reflecting is not the same as reflecting then rotating.
Sankhani wanu
Parts of this page are adapted from Beezer, A First Course in Linear Algebra (GFDL 1.2). Condensed and re-explained here; errors are ours.
Zambiri pa Linear Algebra
DeterminantsMatrix inverseEigenvalues and eigenvectorsMatrix multiplicationRow reductionVector spaces, span and linear independenceOrthogonality, projections and least squaresDiagonalisation and matrix powersLinear transformations and change of basis