Lecture Notes On Matrices
Lecture Notes On Matrices
Students
Dr. K. M. Udayanandan
Associate Professor
Nehru Arts and Science College
Kanhangad.
Syllabus
1
Orthogonal Matrices
AT = A−1
which requires
AT A = AAT = I
where I is the identity matrix. An orthogonal matrix A is necessarily in-
vertible (with inverse A−1 = AT ), unitary (A−1 = A† ) and therefore normal
(A† A = A A† ). The determinant of any orthogonal matrix is either +1 or 1.
Questions
Then
cos θ − sin θ 0
S T = sin θ cos θ 0
0 0 1
So
cos θ sin θ 0 cos θ − sin θ 0
SS T = − sin θ cos θ 0 sin θ cos θ 0
0 0 1 0 0 1
cos2 θ + sin2 θ
− sin θ cos θ + sin θ cos θ 0
= − sin θ cos θ + sin θ cos θ sin2 θ + cos2 θ 0
0 0 1
1 0 0
= 0 1 0 =I
0 0 1
2
SS T = I
So direction cosines of a three dimensional coordinates constitutes an
orthogonal matrix.
ST S = I
B = AT
Then
B T B = (AT )T AT
= AAT = I
BT B = I
Therefore transpose of an orthogonal matrix is also orthogonal.
Hermitian matrices
Properties
3
2. The difference of a square matrix and its conjugate transpose (C − C † )
is skew-Hermitian (also called anti hermitian, A = −A† )
Problems
note:
A column matrix is represented by |Xi called ket.
A row matrix is represented by hX| called bra.This method is called
Bracket method.
The conjugate transpose (also called Hermitian conjugate) of a bra is
the corresponding ket and vice versa: hA|† = |Ai, |Ai† = hA|
Proof:
We have
AX = λX
for any matrix with eigenvalue λ. Representing by bracket method
(A |Xi)† = (λ |Xi)†
|Xi† A† = λ∗ |Xi†
hX| A† = λ∗ hX| (2)
4
Multiplying eqn (1) by hX| from left side, We get
Here A is Hermitian, A† = A
Then (4) ⇒
hX| A |i = λ∗ hX|Xi (5)
Comparing eqn (3) and (5),LHS are equal, hence RHS must be equal.
λ hX|Xi = λ∗ hX|xi
λ = λ∗
Thus for any Hermitian matrices, eigenvalues are real.
2. Show that for any Hermitian matrix eigenvectors corresponding to dis-
tinct eigenvalues are orthogonal to each other.
A X 1 = λ1 X 1
5
Multiplying eqn (7) by hX| from left side,
λ1 − λ2 6= 0
Then
hX1 |X2 i = 0
That is eigenvectors are orthogonal.
A = (A∗ )T
B = (B ∗ )T
Substituting for AB+BA we get
AB + BA = (A∗ )T (B ∗ )T + (B ∗ )T (A∗ )T
= (B ∗ A∗ )T + (A∗ B ∗ )T
= (B ∗ A∗ + A∗ B ∗ )T
= [(BA)∗ + (AB)∗ ]T
= [(AB)∗ + (BA)∗ ]T
AB + BA = [(AB + BA)∗ ]T
which is definition of Hermitian matrix. Hence AB+BA is a Hermitian
matrix.
6
Consider two Hermitian matrices A and B.then by definition
A = (A∗ )T
B = (B ∗ )T
if AB is Hermitian,
AB = [(AB)∗ ]T
But (AB)∗ = A∗ B ∗ ⇒
AB = (A∗ B ∗ )T
also (AB)T = B T AT ⇒
AB = (B ∗ )T (A∗ )T
ie,
AB = BA
hence A and B commute.
Unitary matrices
U †U = U U † = I
or
U −1 = U †
where I is the identity matrix and U † is the conjugate transpose of U.
7
problems
√ √
√ 2/2 −i√ 2/2 0
1. Show that A = i 2/2 − 2/2 0 is a unitary matrix.
0 0 1
we have for a matrix A to be unitary AA† = I = A† A
√ √
√2/2 −i √ 2/2 0
A = i 2/2 − 2/2 0
0 0 1
Taking conjugate
√ √
2/2
√ i √2/2 0
A∗ = −i 2/2 − 2/2 0
0 0 1
Diagonalisation of matrices
a1 b1
Consider a square matrix say of order 2, A = . Let λ1 and λ2 be its
a2 b2
x1 x2
eigen value and , be corresponding eigen vectors.Constructing
y1 y2
8
a matrix P by writing the eigen vector as columns, We get
x1 x2
P =
y1 y2
Then
a1 x1 + b1 y1 a1 x2 + b1 y2
AP =
a2 x1 + b2 y1 a2 x2 + b2 y2
Now consider the eigen value equations,
AX = λX
a1 b1 x1 x1
= λ1
a2 b2 y1 y1
and
a1 b1 x2 x2
= λ2
a2 b2 y2 y2
Equating both sides, we get the equations
a1 x1 + b1 y1 = λ1 x1
a2 x1 + b2 y1 = λ1 y1
a1 x2 + b1 y2 = λ2 x2
a2 x2 + b2 y2 = λ2 y2
Substituting these values
λ1 x 1 λ2 x 2
AP =
λ1 y 1 λ2 y 2
x1 x2 λ1 0
AP =
y1 y2 0 λ2
λ1 0
Let = D, Then
0 λ2
AP = P D
Then
−1 λ1 0
P AP = D =
0 λ2
, which is nothing but diagonalised matrix.
9
So, If you want to diagonalise a diagonalisable matrix, find its eigen values
and write it as a diagonal elements of corresponding dimension.
Problems
6 −2 2
1. Diagonalise the matrix A = −2 3 −1
2 −1 3
We have
AX = λX
(A − λI)X = 0
|A − λI| = 0
6 − λ −2 2
−2 3 − λ −1 =0
2 −1 3 − λ
Simplifying, we get λ1 = 2, λ2 = 2 and λ3 = 8
Then diagonal matrix is
λ1 0 0 2 0 0
D= 0 λ2 0 = 0 2 0
0 0 λ3 0 0 8
0 0 1
2. Find the eigenvalues and eigenvectors of the matrix A given by A= 0 1 0
1 0 0
−1
and obtain the matrix S such that S AS is diagonal.
We have
0 0 1
A= 0 1 0
1 0 0
AX = λX
(A − λI)X = 0
|A − λI| = 0
−λ 0 1
0 1−λ 0 =0
1 0 −λ
10
We get
λ1 = 1 λ2 = 1 λ3 = −1
where λ1 ,λ2 and λ3 are the eigen values. Then diagonal matrix is
λ1 0 0 1 0 0
D= 0 λ2 0 = 0 1 0
0 0 λ3 0 0 −1
U †U = U U †U = I (10)
U X = λX (11)
(U X)† = (λX)†
X † U † = λ∗ X † (12)
Multiplying (10) and (11) ⇒
X † U † (U X) = λ∗ X † (λX)
X † (U † U )X = λ∗ λX † X
Substituting equation (10) ⇒
X † (I)X = λ∗ λX † X
X † X = λ∗ λX † X
(λ∗ λ − I) X † X = 0
X † X 6= 0 ⇒
λ∗ λ − 1 = 0
λ∗ λ = 1
|λ|2 = 1
4. Show that the eigenvalues of a real symmetric matrix are real and
11
eigenvectors with distinct eigenvalues are orthogonal.
We have AX = λX for any matrix with eigenvalue λ. Representing by
bracket method
A |Xi = λ |Xi (13)
(A |Xi)† = (λ |Xi)†
|Xi† A† = λ∗ |Xi†
since A is real matrix,A† = AT
λ hX|Xi = λ∗ hX|xi
λ = λ∗
Thus for any real symmetric matrices, eigenvalues are real.
Consider a matrix A. Let λ1 and λ2 be two distinct eigenvalues of A
and X1 and X2 be corresponding eigenvectors. Then we can write
A X 1 = λ1 X 1
12
Multiplying eqn (18) by hX| from left side,
λ1 − λ2 6= 0
Then
hX1 |X2 i = 0
That is ,eigenvectors are orthogonal.
x + 2y + 3z = λx
3x + y + 2z = λy
2x + 3y + z = λz
The equations can be represented using matrices as follows,
AX = λX
13
1 2 3 x x
3 1 2 y = λ y
2 3 1 z z
AX − λX = 0
(A − λI)X = 0
since λ 6= 0 we can write
|A − λI| = 0
ie,
1−λ 2 3
3 1−λ 2 =0
2 3 1−λ
√
−3 ± 9 − 12
λ=6 and λ=
2
√
−3 ± i 3
λ=
2
Hence real value of λ is 6.
Answer:
AX = λX
(A − λI)X = 0
|A − λI| = 0
7−λ −2 1
−2 10 − λ −2 =0
1 −2 7−λ
λ1 = 6 λ2 = 6 λ3 = 12
14
Then diagonal matrix is
λ1 0 0 6 0 0
D = 0 λ2 0 = 0 6 0
0 0 λ3 0 0 12
3. A 3×3 matrix has elements such that its trace is 11 and its determinant
is 36. The eigenvalues of the matrix are all known to positive integers.
What is the largest eigen value of the matrix?
2 3 0
4. Find the eigenvalues of 3 2 0
0 0 1
0 i
5. Find the eigenvalues of
i 0
15
aeiα b
6. is a unitary matrix. If a, b, b, d, α, β are real , what is the
ceiβ d
inverse of the matrix?
16