Matrix Algebra
Matrix Algebra
Matrix Algebra
INTRODUCTION TO MATRICES
Reference : Croft, A., & Davison, R. (2008). Mathematics for
Engineers - A Modern Interactive Approach, Pearson
Education.
where aij represents the number or element in the ith row and
jth column.
find A + B, B + C and B - C.
3 5 7 6 10 11
B + C = 2 1 3 9 5 8
4 2 1 5 5 7
3 5 7 6 4 1
2 1 3 9 1 10
B–C=
4 2 1 5 3 3
Introduction to Matrix Algebra Page 5
Multiplication of a Matrix by a Number
Any matrix can be multiplied by a number. To do this, each
element of the matrix is multiplied by that number.
7 13 5
e.g.2 If A 9 8 2 , find 2A, -A.
4 11 8
2 * 7 2 * 13 2 * 5 14 26 10
2A = 2 * 9 2 * 8 2 * 2 18 16 4
2 * 4 2 * 11 2 * 8 8 22 16
1* 7 1*13 1* 5 7 13 5
-A = 1* 9 1* 8 1* 2 9 8 2
1* 4 1*11 1* 8 4 11 8
cij a b
k 1
ik kj
3
i.e. c21 a2 k bk 1 8 * 1 2 * 3 6 * 2 26
k 1
Determinant
All square matrices, A, possess a determinant denoted by :
det(A), |A|.
Determinant of a 2 x 2 matrix
a b a b
If A , then det(A) = |A| = = ad - bc
c d c d
The minor of aij is the determinant formed by crossing out the ith
row and jth column of det(A).
Note that the term (-1)i+j is called the place sign of the element
on the ith row and jth column. The following may help you to
memorize this.
Introduction to Matrix Algebra Page 11
Determinant of a 3 x 3 Matrix
a11 a12 a13
Consider a general 3 x 3 matrix, A = a21 a22 a23
a a33
31 a32
1 4 17 24 5 4 17 4 17
11 24 5 1* 11 * 6*
31 15 31 15 24 5
1* 205 11 * 467 6 * 388 3014
6 31 15
Repeat the 1st and 2nd column to right hand side of 3rd
column to form a 3 x 5 matrix.
14
5 7 3
Introduction to Matrix Algebra Page 19
3 1 0
e.g. 3 Find the inverse of B 5 2 1 .
1 6 3
det( B ) 14
T
2 1 5 1 5 2
6 3 1 3 1 6
0 3 1
1 0 3 0 3 1
adj ( B) 14 9 3
6 3 1 3 1 6
1 28 17 1
0 3 3 3 1
2 1 5 1 5 2
0 3 1
1
B 14 9 3
1
14
28 17 1
Introduction to Matrix Algebra Page 20
TE CHNI Q UE S O F S O LV I NG
AL G E BRAI C E Q UATI O NS
2 x 3 y 3 2 3 3 x 3
i. (unique solution)
0 x y 1 0 1 1 y 1
a2 b2 a2 b2
Techniques of Solving Algebraic Equations Page 24
ii. 3 equations in 3 unknowns
a1 x b1 y c1 z k1 a1 b1 c1
If a2 x b2 y c2 z k 2 where a2 b2 c2 0
a x b y c z k
3 3 3 3 a3 b3 c3
then
k1 b1 c1 a1 k1 c1 a1 b1 k1
k2 b2 c2 a2 k 2 c2 a2 b2 k2
k3 b3 c3 a3 k3 c3 a3 b3 k3
x , y , z
a1 b1 c1 a1 b1 c1 a1 b1 c1
a2 b2 c2 a2 b2 c2 a2 b2 c2
a3 b3 c3 a3 b3 c3 a3 b3 c3
a1 x b1 y k1
Note that
a2 x b2 y k 2
a1 b1 x k1
can be written as
a2 b 2 y k2
a1 b1 c1 x k1
can also be written as a2 b2 c2 y k 2
a c3 z k3
3 b3
AX = B
A-1AX = A-1 B
I X = A-1 B
X = A-1 B
Introduction
Gaussian Elimination is a systematic way of simplifying a
system of equations.
A matrix, called an augmented matrix, which captures all the
properties of the equations, is used.
A sequence of elementary row operations on this matrix
eventually brings it into a form known as echelon form (to be
discussed later).
From this, the solution to the original equations is easily found.
a1 x b1 y k1
a2 x b2 y k 2
a1 b1 k1
a b k2
2
2
coefficien ts constants
this vertical line can be
omitted as in your textbook
a1 x b1 y c1 z k1
a2 x b2 y c2 z k 2
a x b y c z k
3 3 3 3
a1 b1 c1 k1
a2 b2 c2 k 2
a c3 k3
3 b3
9 x 7 y 65 z 15 9 7 65 15
b. 5 x 12 y 2 z 64 5 12 2 64
4 x 22 y 3 z 7 4 22 3 7
17 x 3 y 6 17 3 0 6
c. 4 x 8 y 3z 12 4 8 3 12
0 3 5 7
3 y 5z 7
a. 1 7 15 x 7 y 15 x 1
0 1 2
y 2 y 2
1 2 1 5 x 2 y z 5 x 2
b. 0 1 3 1 y 3z 1 y 2
0 0 z 1 z 1
1 1
1 2 5 78 1 4 1 5
a. 0 1 8 14 Yes b. 0 0 1 2 No
0 0 1 22 0 1 1 27
2 1 2 8
The augmented matrix is 1 3 3 4
4 2 1 1
1 3 3 4
Interchange row 1 and row 2
2 1 2 8
4 2 1 1
1 3 3 4 row 2 – 2*row 1
0 7 4 16 row 3 – 4*row 1
0 14 13 17
1 3 3 4
row 2 / 7
0 1 7
4 16
7
0 0 row 3 / -5
1 3
Hence z 3
4 16
y z or y 4
7 7
x 3 y 3z 4 or x 1
The solution is x z 1 4 3
T T
y
Techniques of Solving Algebraic Equations Page 45
4 x 3 y 5
e.g.3 Use Gaussian elimination to solve
3 x 2 y 3
4 3 5
3 2 3
4 3 5
0 1 3 3R 4 R
1 2
1 3 4 5 4 R1 4
0 1 3
x 3 y 4 5 4
y 3
x 1
y 3
Techniques of Solving Algebraic Equations Page 46
Class Practice
47
1 2 2 1 0 0 1 2 0 11 6 2 R1 2 R3
0 1 2 2 1 0 0 1 0 8 5 2 R2 2 R3
0 0 1 5 3 1 3R R 0 0 1 5 3 1
2 3
1 0 0 27 16 6 R1 2 R2
0 1 0 8 5 2
0 0 1 5 3 1
27 16 6
1
Hence A 8 5 2
5 3 1
Techniques of Solving Algebraic Equations Page 50
,
Example
51
A=
1 2 1 1 0 0
0 3 0 1 1 0
0 0 1 1 2 1 R 2 R
3 2
1 2 0 0 2 1 R1 R3
0 1 0 1 3 1 3 0 R2 3
0 0 1 1 1
2
1 0 0 2 3 4 3 1 R1 2 R2
0 1 0 1 3 1 3 0
0 0 1
1 2 1
2 3 4 3 1 2 4 3
1
1
Hence A 1 3 1 3 0 1 1 0
1 3
2 1 3 6 3
Techniques of Solving Algebraic Equations Page 52
EIGENVALUES AND EIGENVECTORS
a b x 0 or simply, AX 0
c d y 0
AX
AX
a b x x AX X
or
c d y y
solutions.
The system can be written as A I X 0
The system has non-trivial solutions if A I 0
4 5
det A I 2 7 2 0
2 3
7 7 2 412 7 41
Hence
21 2
AX X .
x 4
The eigenvector is t
y 1
Eigenvalues and Eigenvectors Page 63
1 , 5 1 4 4 4
When ( A I )
1 2 1 1 1
i.e. x y 0 or x y . In general, x t , y t
x 1
The eigenvector is t
y 1
3 1
A
1 5
1 2 0
A 1 1 1
3 2 2
3 4 1 1 1
When 4 , ( A I )
1 5 4 1 1
i.e. x y 0 or x y . In general, x t , y t
x 1
The eigenvector is t , for any number t
y 1
Eigenvalues and Eigenvectors Page 66
e.g.3 Find the eigenvalues and eigenvectors of A
1 2 0
A 1 1 1
3 2 2
1 2 0
det (A-I) = 1 1 1 0
3 2 2
1 1 1 1
1 2 0
2 2 3 2
1 2 3 2 1 0
1 2 3 2 0
1 1 2 0
Hence 1 or - 2
Eigenvalues and Eigenvectors Page 67
e.g.3 Find the eigenvalues and eigenvectors of A
1 1 2 0 0 2 0
When 1, A I 1 1 1 1 1 2 1
3 3
2 2 1 2 3
i.e. y 0 and x z 0 or x z
In general, x t , y 0, z t
x 1
The eigenvector is y 0 t , for any number t
z 1
Eigenvalues and Eigenvectors Page 68
e.g.3 Find the eigenvalues and eigenvectors of A
1 1 2 0 2 2 0
When 1, A I 1 1 1 1 1 0 1
3 3 2 1
2 2 1
i.e. x y 0 or x y and x z 0 or x z
In general, x t , y -t , z t
x 1
The eigenvector is y 1t , for any number t
z 1
Eigenvalues and Eigenvectors Page 69
e.g.3 Find the eigenvalues and eigenvectors of A
1 2 2 0 3 2 0
When 2, A I 1 1 2 1 1 1 1
3 3 2 0
2 2 2
2
i.e. 3 x 2 y 0 or x y and x y z 0 or z x y
3
3 5
In general, x t , y - t , z t
2 2
x 1
The eigenvector is y 3 2 t , for any number t
z 5 2
Eigenvalues and Eigenvectors Page 70
APPLICATION
Stiffness method is commonly used to analyse a structure as
matrix analysis can be used to solve the problem.
i j
4 EI 2 EI
qi i j 4 EI 2 EI
L L
qi L L i 2 EI 2 1 i k
q 2 EI 4 EI j 1 2
2 EI 4 EI j L j
qj i j L L
L L
qi
[q] - force matrix
q j
5 3
det A I 2 9 14 0
2 4
Hence 2 or 7
3 1 4 3 4 4 3
When 1, A I 5 6 1 3 5 5 3
3 3
3 2 1 3 1
i.e. z 0 and x y 0 or x y
In general, x t , y -t , z 0
x 1
The eigenvector is y 1t , for any number t
z 0
Eigenvalues and Eigenvectors Page 75
Hence 1 or 2 (repeated )
3 2 4 3 5 4 3
When 2, A I 5 62 3 5 4 3
3 3
3 2 2 3 0
5x 4 y
i.e. x y 0 and 5 x 4 y 3 z 0 or z
3
t
In general, x t , y -t , z
3
x 1
The eigenvector is y 1 t , for any number t
z 1 3
Eigenvalues and Eigenvectors Page 76