ES 21 First Exam Notes
ES 21 First Exam Notes
ES 21 Notes
MATRICES
c1
c
2
Definition
A matrix is a rectangular array of numbers or functions
arranged in rows and columns usually designated by a
capital letter and enclosed by brackets, parentheses or
double bars. A matrix may be denoted by:
A
:
:
:
1 i m
C.
a1 j
a
2j
j n
D.
We refer to aij as the entry or the element in the ith row and
jth column of the matrix.
.
.
.
a3 j
A.
c
m
ai 2 ... ain
ci
ai1
.
.
.
E.
bn]
l11 0 0
L l 21 l 22 0
l 31 l32 l 33
Example:
F.
0
d11 0
D 0 d22 0
0
0 d33
ES 21 Notes
Page 1 of 13
A.
Example:
s 0 0
S 0 s 0
0 0 s
H.
Example:
1 2 4
A
3 2 1
1 0 0
I3 0 1 0
0 0 1
B.
SCALAR MULTIPLICATION
2 0 1 3
6 0 3 9
3
4 3 2 5
12 9 6 15
Example:
4
1 2
S 2 2 5
4 5 3
K.
2
4 5
2
Example:
0 2 4
T 2 0 5
4 5 0
1 3
8 4 12
2 1 20 8
4
8 16 28
4 7
MATRIX SUBTRACTION
EQUALITY OF MATRICES
Example:
3 4 2
A
5 7 4
Two matrices A = (aij) and B = (bij) are equal if and only if the
following conditions are satisfied:
a)
b)
c)
1
c 2
y 1
2 andB
b
z
2
1
x 1 b 2
2 1 4
3 8 3
3 4 2
2 1 4
A B
5 7 4
3 8 3
4 (1) 2 4
3 2
7 8
4 3
5 (3)
x
A a
3
3 2 2
4 1 1
4 (2)
1 1
4 0 2
1 1 2
0 0 ... 0
0 0 ... 0
: :
:
0 0 ... 0
J.
3 2 2
B
4 1 1
1 2 4
A B
3 2 1
2 2
1 3
3 4 2 (1)
Example:
I.
1
2
1 5 6
8 1 7
c 1
MATRIX MULTIPLICATION
cij
aikbkj
MATRIX EXPONENTIATION
1(2) 2(3)
AB
3(2) 4(3)
6
8
AB
18 16
then A
1 1 2
2
3 1
2 4
3
SCALAR MULTIPLICATION
0 xA= O
1 x A= O
k l (A) = k (l A) = l (k A)
(k + l) A = k A + l A
k(A + B) = k A + k B
A(BC) = (AB)C
A(B + C) = AB + AC
(A + B)C = AC + BC
AI = IA = A
2 4
B
3 1
Associative Property
Left Distributive Property
Right Distributive Property
Existence of Multiplicative
Identity
1(4) 2(1)
3(4) 4(1)
2(1) (4)(3)
BA
3(1) (1)(3)
10 12
BA
2
0
2
4
1 3
3 2
2
3
MATRIX MULTIPLICATION
1 2
A
3 4
k1
Example:
for i = 1 to m and j = 1 to p
Note:
If
1
1
Note:
is
not
MATRIX TRANSPOSITION
(AT)T = A
(A + B)T = AT + BT
(k A)T = k AT
(AB)T = BTAT
2(2) (4)(4)
3(2) (1)(4)
Example:
Permutation:
2 1
A 5
3
4 6
1.
2.
3.
E.
6 2 3
B
1 4 5
A x B is a 3 x 3 matrix while B x A is
a 2 x 2 matrix.
A x C is a 3 x 2 matrix but C x A is not
defined.
B x C is not defined but C x B is
defined (2 x 3).
MATRIX TRANSPOSITION
Let
3 2, 3, , n} be the set of integers from 1 to n,
1 S={1,
C arranged
2 1 in increasing order. A rearrangement a 1a2a3an of
the
elements
in S is called a permutation of S.
By the Fundamental Principle of Counting we can put any
one of the n elements of S in the first position, any one of
the remaining (n-1) elements in the second position, any
one of the remaining (n-2) elements in the third position,
and so on until the nth position. Thus there are n(n-1)(n-2)
3*2*1 = n! permutations of S. We refer to the set of all
permutations of S by Sn.
Examples:
If S = {1, 2, 3} then S3 = {123, 132, 213, 231, 312, 321}
If S = {1, 2, 3, 4} then there are 4! = 24 elements of S4.
Odd and Even Permutations
A permutation a1a2a3an is said to have an inversion if a
larger number precedes a smaller one. If the total number of
ES 21 Notes
Page 3 of 13
a11
2.
DETERMINANT
det A A
a1j a2 j ...anj
1
B.
a11a22 a21a12s
3 x 3 Matrices
a11 a12
DEFINITION:
a12
a21 a22
2 x 2 Matrices
a13
a11 a12
a21 a22
a31 a32
METHOD OF COFACTORS
If A =
a11 a12
a
21 a22
If A =
21 a22 a23
a31 a32 a33
To evaluate
1
2
0
1
3
4
0
1
numbers
of (zeros.
course
(2)(2)(2) (3)(1)(3) (0)(1)(1) most
(0)(
2)(3)
2)(1The
)(1)optimal
(3)(
1)(2)of action is to
8 9 0 0 2 6
9
DIAGONAL METHOD
ES 21 Notes
Page 4 of 13
2.
1
2
2
3
0
00
3.
1
4
2
2
3
3
2
1
2
4
2
1
4
Example:
4
4.
a
x
b
1
y 5 then4
a
2x
2a
2b
2a
2xy
2
4x y
6xy
5.
3x
2x
6y
18
8 y 2xy2x
2xy
3
6y
18
6.
2x y
6x2
4z
3x y
5z
xyz
3x
2z
3y
3
a
2a
4a
6b
a b
3
6
1
2
4 2 and
2xy
5
5 3
2b2
4a
6b
a b
5
5
2b2
x2
0
xy
x2
x y z y
0
0 x2 4
1
3x ofya square
5z matrix
xzyA=[aij] of order n is that
Adjoint
square matrix with the same order n denoted by adj(A)=[A ji]
2x
y Aij is
3x
2z of the xelement a ij of matrix A. The
where
the cofactor
0
adjoint 3
ofx
a matrix
2x
y is the
5ztranspose
xyz of the matrix of cofactors
of the elements of A.
2
2x
The
3x
xy2
4z
2
4
6y
8.
0
0
1
2x
7.
x2
0
xy
x2
x y z y
a11 a12
Example:
0
0
1
a13
3
a
2x
8 y 2xy6
(2y
) x 13 y 2xy
4 y 2a
2
3
= 6
2xy
3. 18 2xy
xy
x2
x2
0
x y z y
a12
then 4 y 6z3x 4
2x
3x
a23
a33
2
4 0
2xy
3y
xz
xy
a13
a11
b21 c21 b22 c22 b23 c23 b21 b22 b23 c21
a31
a32
a33
a31 a32 a33
a31
2x
Example:
a22
a32
b
2x
3
2 y 2(5) 10
6z
If 4
2b
6y
1
a21
a31
9.
a11
1
If 2
A ij (1)i j Mij
Step 2: Set up the adjoint matrix by taking the transpose
of the matrix of cofactors.
ES 21 Notes
Page 5 of 13
adjA A ij
Example:
a
c
If A =
b
d
then adj(A) =
d b
c a
Inverse of a Matrix
A11 (1)2
-1
A12 (1)3
-1
A(A ) = (A )A = In.
Not all matrices has its inverse. However, if the inverse of a
matrix exists, it is unique.
If the inverse of a matrix exists, we say that the matrix is
invertible or non-singular. Otherwise, we say that the matrix
is non-invertible or singular.
Matrix Inversion Using the Adjoint
and the Inverse
adjA
A
1
3
1
3
A31 (1)4
A 33 (1)6
1
2
4
A 22 (1)4
A 32 (1)5
A 21 (1)3
A 23 (1)5
A 1
A13 (1)4
1 1
1 3
2
5
1 1
1
10
1 1
2
12
12 9 10
Aij 1 1 2
,thus
adj
A
5 2
3
adj A
1
Consequently, A
A
1 1 1
A 2 5 2
1 7 7
12
9
10
1
2
1 5
12
9
1 2
10 2 3
2
3
A 35 2 14 5 14 14 0
Since matrix A is singular, as evidenced by its zero
determinant, it can thus be concluded that the
Inverse of A (or A-1) does not exist.
Example 2:
matrix
1
A 1
4
1
2
2
1
3
3
A 6 12 2 8 6 3 7
ES 21 Notes
Page 6 of 13
a11
a
21
a31
a m1
a12
a13
a 22
a32
a 23
a33
am2
a m3
a1n
a 2 n
a 3n
a mn
x1
x 2
x3
x n
b1
b
2
b3
bm
12
9
1
1
10
5
2
3
x1
12 1 5
1
X x 2
9 1 2 *
7
x3
10 2 3
1
6
5
Where:
A
X
B
Coefficient Matrix
Column Matrix of Unknowns/Variables
Column Matrix of Constants
matrix:
Take
1
1
1
AX B
x1
X x2
x3
note
that the
assumes
that A-1 exists. If A-1
( A 1A ) * X A 1 * B
does not exist, we can
not find the solution to
I * X A 1B
the system AX = B.
A 1 AX A 1 * B derivation
x1 2 x 2 3 x3 6
an1 an2 an3 ann
4 x1 2 x 2 3 x3 5
Let:
X A 1B
x1 x 2 x3 1
A
xi
B
AI
Solution:
The above system of equations can be written in matrix
form:
1 1 1 x1
1 2 3 x
2
4 2 3 x3
1
6
5
1
A 1
4
1
2
2
1
3
3
x1
x
2
b1
b
2
x 3 b3
xn
bn
coefficient matrix
ith variable
right hand side constants
matrix resulting from replacing the
ith column of A by the column
vector of constants B
xi
Ai
A
Getting A-1
ES 21 Notes
Page 7 of 13
1 2 1
3 3 1
1 2 1
x1
x
2
x3
1
2
3
Solution:
Compute for the determinant of A first:
or
1 2 1
A 3 3 1
1 2 1
A 3 2 6 3 2 6 6
A1
A
1
2
3
or
or
l 42 a 42 l 41 (u12 )
A
0
x1 1 0
A
6
Applying the same process to solve x2 and x3:
1 1 1
A2
1
x2
* 3 2 1 1
A
6
1 3 1
a14 l11 0
a 24 l 21 l 22
a34 l 31 l32
a 44 l 41 l 42
or
or
1 u2 1[a l (u )] / l
24
24
21
14
22
A3
1
x3
* 3 3 2 1
A
6
1 a2
3l 31 (u13 ) l 32 (u 23 ) l 33 (1) 0(0)
33
a13
a 23
a33
a 43
A1 3 6 4 9 2 4 0
a12
a 22
a 32
a 42
l 32 a 32 l 31 (u12 )
1 2 1
A1 2 3 1
3 2 1
a11
a
21
a31
a 41
l 41 a 41
x1
or
0
0
l 33
l 43
0
0
*
0
l 44
or
l 33 a33 l 31 (u13 ) l 32 (u 23 )
a 43 l 41 (u13 ) l 42 (u 23 ) l 43 (1) l 44 (0)
l 43 a 43 l 41 (u13 ) l 42 (u 23 )
or
a 34 l 31 (u14 ) l 32 (u 24 ) l 33 (u 34 ) 0(1)
u 34 [ a34 l31 (u14 ) l32 (u 24 )] / l33
or
4 2 1
3 5 1
1 2 1
Knowing that
x1
x
2
x3
3
4
8
4 2 1
A 3 5 1
1 2 1
0 0 y1 34 y1 y31 3 4
3 y 7 y y4 1
7
3 2 0 y 2 y 4 51 y 2 52 y2 8y23 3
1 2 2 2 3
5
5
y
8
1
2 2 3
Note that the computed values of yi's here are not yet the
solution since the original system of equations is in terms of
xi's.
Stage 2: Back substitution using UX = Y
0
0
x1
x
2
x3
x3 3
1
1
1
as simply
1
2
4
1
1
2
0
1
1 3 x
2 6 y
2 4 1 z
2
1 3 | 1
1
2 6 | 1 .
2 4 1 | 1
2
1
, therefore
4 0 0
L 3 72 0
5 5
1 2 2
1 1 1
2
4
1
and
U 0 1 2
0 0
1
1 1
1 2 5
2 1
1 6 and B 0
2
If A
7
3 3 8
7 4
1 2 5 | 2 1
7
1 6 | 0
2
then A : B is
3 3 8 | 7 4
x 3
3
3
x 2 1 4x 3 1 x 2 2
x 11 2 x 1 2 x x13 1
1
4
22 2 4 3
3
This time (x1, x2, x3) = (1, 2, 3) is the solution to the original
system of equations.
AUGMENTED MATRIX OF A AND B
If A is an m x n matrix and B is a p x n matrix, then the
augmented matrix of A and B denoted by [A : B] is the
matrix formed by the elements of A and B separated by
pipes.
Example:
The following matrices are not in row echelon form. (Why
not?)
1
0
A 0
0
0
2
1
0
0
0
0 1 3
0 1 3
0 1
9
0 0
0
0 0
0
2
5
0
1
1
0
B 0
0
0
0 0
1 2
0 1
0 0
0 0
0
0
1
0
1
0
D
0
2
1
0
0
3
2
1
0
4
3
1 0 5 6 2
E 0 1 0 0 8
0 0 1 0 0
Example:
ES 21 Notes
Page 9 of 13
1
0
F 0
0
0
1
0
G
0
0
1
0
0
0
0
1
0
0
0
0 1 0 0 2
H 0 0 1 0 3
0 0 0 1 0
1 0 3
ON MATRIX EQUIVALENCE
THEOREMS
0 1 0
1. Every nonzero
0 0 0
2.
Every nonzero m x n matrix A = [aij] is row (column)
0 equivalent
0 0 to a matrix in reduced row (column) echelon
form.
3.
4.
0 2 0
1
A 3 1
6 4
2 8 2
2
a11
a
21
a31
2
2 8 2
B 3 1 6 4
1 0 2 0
a m1
0
1
C 3 1
1 4
2 0
6 4
1
1
am2
a m3
x1
b1
x
b
2
2
x3 b3
bm
a mn x n
a12
a22
a13 a1n| b1
a23 a2n| b2
a31
a32
a33 a3n| b3
am3 amn| bm
a 23
a33
am1 am2
TY0PE III
a 22
a32
a1n
a 2 n
a 3n
a11
a
21
1 0 2 0
A 0 1
0 4
2 8 2
2
a13
TYPE I
TYPE II
a12
ES 21 Notes
Page 10 of 13
2y
3z
2x
3x
z
z
8
3
1
A : B 2
3
9
8
Example: Let
1 |
9
2
3
3
4
y 2 3 1
x 2 y 3z 9
x 9 2(1) 3(3) 2
3.
4.
2y
3z
2x
3x
z
z
8
3
9
8
3
1 0 0 |
: B* 0 1 0 |
0 0 1 |
2
1
3
9
5
1 8
0 3 8 2
11 5 7 2
3 0
0 3
6
9
11 5
(remove
thirdrow)
8
2 (remove
first
, third
2
andfifthcolumn
)
1
8 (fourthco
5
remain
RANK OF A MATRIX
The rank of a matrix A = [aij] is the order of the largest
square submatrix of A with a non-zero determinant. We
denote the rank of A by rank(A) or simply r(A).
Example: What is the rank of A?
4
1 2 3
A 2 1 4
3
3 0 5 10
Solution:
Checking out first the determinants of 3x3 submatrices:
1
2
2
1
3
4 0
1
2
2
1
4
3
10
1
2
3
4
4
3
10
3 |
1 2
A : B 2 1 1 |
3 0 1 |
6
11
z3
y z 2
2.
8
2
5 2
7
2
1
8
1.
3 0
0 3
1 2 3 |
*
*
A : B 0 1 1 |
0 0 1 |
2
3
Solution:
2
3
4
1
2 4
6
8
B
3
6
9
12
4 8 12 16
ES 21 Notes
Page 11 of 13
14
e.g.
2
3
4
6
8
12 0 (Rows are proportional)
9
1
16
e.g.
2.
24 0
3.
3.
2
3
4
5
6
3
4
5
6
7
4
5
6
7
8
5
6
8
9
1 2 3 4 5
1 1 1 1 1
C' 1 1 1 1 1
1 1 1 1 1
1 1 1 1 1
:
:
R3' R3 R2
R4' R4 R3
R5' R5 R4
1 :
1 :
0
1
3 0
0
6
0
1
0
0
:
:
R2' R2 R1
1
1
1 0 0 : 5
2 1 1 : 3
4
4 2 2 : 2 0 1 0 : 3
6 1 1 : 6
0
0
0
:
2x
3y
3z
2k
2y
y
4z
5z
5k
8k2
We could easily see that all 5x5, 4x4 and 3x3 submatrices
of C have determinants equal to zero (THEOREM: Identical
rows). But for at least one 2x2 submatrix of C has a nonzero determinant.
1 0 )
2
1
3 3 | 2k
2 4 | 5k
1 5 | 8k 2
0
1
Solution:
Operating on rows of matrix C, we obtain the equivalent
matrix C
(e.g.
1
0
Example:
1
2
C 3
4
5
1
3
THEOREMS ON RANKS
2.
1
1
27 0
0
1
1 :
Example:
Therefore r(B) = 2.
1.
3 1 :
1 1 :
1
0
2 3 3 | 2k R1 R 2 1 2 4 | 5k
1 2 4 | 5k
2 3 3 |Performing
Gaussian
Elimination Method:
2
k
0 1 5 | 8k 2
0 1 5 | 81k 2 1 1 | 2
'
R 2 R 2 2R1 1 2 4 | 5k
0 1 5 | 8k
0 1 5 | 8k 2
'
R3 R3 R2 1 2 4 |
5k
0
1
5
|
8
k
0 0 0 | 8k 2 8k
m2 0
1 m2
0
| m 4
|
'
R 2 R 2 ( m 2)R1 1
1
1
|
2
0 m 2 m 1 | 2m 4
'
0
R 3 R 3 R1
0
m 2 1 | m 2
Therefore we have the following conclusions:
Therefore we have the following conclusions:
a) For a unique solution, r(A) = r[A:B] = n
This will be satisfied if m2-10 and m-20. Thus
we have a unique solution if m 1,2.
m2 - 1 = 0 and m+2 = 0.
There is no value of m that will satisfy both
8k2 - 8k = 0 k = 0, 1.
equation.
0
0
8k - 8k 0 k 0, 1.
2
Example 2:
For what values of m will the system of equations have
d) a unique solution
e) a non-unique solution
f) no solution
0
0
1
3
|
|
0
0
c
m 2c
m4
(m 2)a
a
1
1 1
m2 0 1
1
1 m2
|
2
|
0
| m 4
ES 21 Notes
Page 13 of 13