LectureNotes (1)
LectureNotes (1)
Lecture Notes
c 2009-2017
Abdol-Reza Mansouri
These notes are meant to be detailed and expanded versions of the class-
room lectures you will be attending. Each section in these lecture notes will be
covered in one or (often) more classroom lectures.
These lecture notes will be slightly and gradually updated and completed
as the classroom lectures progress. Some parts of the lecture notes may even be
updated and slightly modified after the corresponding classroom lectures have
been delivered. It is therefore recommended that you always refer to the latest
version of these lecture notes, which will be posted on the course website.
It is strongly recommended that you read these lecture notes carefully
as a complement to the classroom lectures. Almost every section in these notes
contains examples which have been worked out in detail. It is strongly recom-
mended that you examine those in detail, and it is equally strongly recommended
that you work out in detail the examples that have been left to the reader. Last
but not least, it is strongly recommended that you attempt to solve the problems
at the end of each section.
I hope you will find these notes useful, and I would appreciate your
feedback. So please feel free to e-mail me at mansouri@[Link] . I look
forward to hearing from you.
Abdol-Reza Mansouri
January 2017
Contents
Section 0 1
Sets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Quantifiers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
Mappings and Functions . . . . . . . . . . . . . . . . . . . . . . . . . . 9
Complex Numbers . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
Section 1 19
Systems of Linear Equations . . . . . . . . . . . . . . . . . . . . . 20
Number of Solutions of Systems of Linear Equations . . . . . . . 21
Applications of Systems of Linear Equations . . . . . . . . . . . . 24
Section 2 33
Section 3 45
Section 4 51
Section 5 59
Section 6 69
Section 7 79
Section 8 95
Section 9 107
Section 10 117
Section 11 131
Section 12 151
Section 13 157
Applications of Matrix Multiplication . . . . . . . . . . . . . . . 165
3
4 CONTENTS
Section 14 175
Section 15 193
Applications of eigenvalues and eigenvectors . . . . . . . . . . . . 203
Section 0
Study Topics
• A Review of Basic Notions
1
2 SECTION 0
In this section, we shall review some of the basic mathematical notions that we
will be using throughout the course.
Sets
One of the basic notions that we shall deal with throughout the course
is that of set. By a set we mean nothing other than a collection of objects;
these objects are then called the elements or points of the set. Some simple
examples of sets are:
• the set of all integers greater than or equal to 0 and less than or equal to
5,
• the set of all real numbers greater than or equal to 0 and less than or equal
to 1,
• ...
Note that in the list of examples given above, the first 5 sets in the list have a
finite number of elements; the last 5, on the other hand, have infinitely many
elements.
Let us now denote the first set in the list above by the capital letter A. By
definition of the set A, the elements (or points) of A are all the integers greater
than or equal to 0, and less than or equal to 5. We can therefore write the set
A explicitely as:
A = {0, 1, 2, 3, 4, 5}.
The curly brace “{” to the left indicates the beginning of the set, and the curly
brace “}” to the right indicates its end. The symbols enclosed by the curly
braces and separated by commas “,” are the elements of the set (namely the
integers 0, 1, 2, 3, 4 and 5).
SETS 3
As far as a set goes, the order in which its elements are written is of no conse-
quence; hence, we can equivalently write the set A as:
A = {2, 5, 0, 4, 3, 1}
or as
A = {0, 5, 2, 1, 3, 4}
or any other permutation of the elements. On the other hand, when writing
down the elements of the set, we shall not allow any repetition of any element;
for example, we shall never write A as {0, 0, 1, 2, 3, 4, 5} or as {0, 1, 2, 3, 4, 5, 4},
or in any other such manner: Each element of the set must appear only once.
Let now S be any set; if some entity, say x, is an element of the set S, then we
write
x∈S
x∈
/S
to denote this. For example, for our set A above, we can write 0 ∈ A, since the
integer 0 is an element of A. On the other hand, the integer 8 is not an element
of A, and therefore we write 8 ∈
/ A. Similarly, the integer 5 is an element of the
set N, and we can therefore write 5 ∈ N.√On the other hand, the real number
√
2 is not an element of N, and we write 2 ∈ / N.
Now that we have defined what we mean by a set and what we mean by an
element of a set, we can define what we mean by two sets being equal:
Definition 1. Let A and B be two sets. A is said to be equal to B if A and B
have exactly the same elements. We write this as A = B.
In other words, we have A = B if and only if the following two conditions are
met:
(i) Every element of A is also an element of B, and
(ii) every element of B is also an element of A.
For example, defining the set E as E = {0, 1, 2, 3} and the set F as F =
{2, 1, 3, 0}, we can write
E=F
4 SECTION 0
since E and F have exactly the same elements, and hence, by our definition, are
equal (recall again that, for a set, the order in which its elements are listed is
of no importance). On the other hand, defining the set G as G = {3, 0, 2, 7, 1},
we can write
E 6= G,
since E and G do not have exactly the same elements; in particular, we can see
that 7 ∈ G whereas 7 ∈/ E, and this shows that the elements of E and G are
not exactly the same.
B = {x ∈ N : x ≥ 5}.
The expression to the right of the equality sign should be read as “the set of
all elements x in N such that x ≥ 5”; in that expression, the colon “:” is to be
read as “such that”. It is worth noting that in some texts the vertical bar “|”
is used in place of the colon “:”, and hence the set B could also be written as
B = {x ∈ N | x ≥ 5},
with the vertical bar “|” to be read as “such that” just like we did with the
semicolon. These are some of the standard ways of defining sets. Note that
this way of defining sets is not restricted only to sets having infinitely many
elements; for example, defining the set A as:
A = {5, 2, 3, 4, 6},
A = {x ∈ N : x ≥ 2 and x ≤ 6};
since every element of the set to the left is also an element of the set to the
right, and vice-versa. Similarly, we also have the equalities
and
and so on. These are all equivalent ways of defining the set A.
As another example, consider the set [0, 1] consisting of all real numbers greater
than or equal to 0 and less than or equal to 1. We can write:
N ⊂ R,
On the other hand, it is not true that every element of R is also an element of
N, and hence R is not a subset of N.
Using the “subset” relation, we can therefore state that two sets A and B are
equal if and only if we have both A ⊂ B and B ⊂ A; we can write this more
formally as:
A = B ⇔ A ⊂ B and B ⊂ A.
In the expression above, the symbol ⇔ should be read as “if and only if” or “is
equivalent to”.
Given two sets S and T , we can construct a new set from them as follows:
Definition 3. Let S and T be sets. We denote by S ∩ T the set of all elements
which are both in S and in T . We call S ∩ T the intersection of the sets S and
T.
For example, if A = {0, 1, 2, 3}, B = {2, 3, 4, 5}, and C = {5, 6, 7}, then we can
write:
The following properties of set intersection are easy to verify (and strongly
recommended to the reader):
Lemma 1. Let S, T, U be sets. We have:
S ∩T = T ∩ S,
S ∩ (T ∩ U ) = (S ∩ T ) ∩ U.
6 SECTION 0
Given two sets S and T , we can construct a new set from them in yet another
way:
Definition 4. Let S and T be sets. We denote by S ∪ T the set of all elements
which are in either S or T or both. We call S ∪ T the union of the sets S and T .
For example, for the sets A, B, C defined just above, we have:
The following properties of set union are easy to verify (and their verification is
again strongly recommended to the reader):
Lemma 2. Let S, T, U be sets. We have:
S∪T = T ∪ S,
S ∪ (T ∪ U ) = (S ∪ T ) ∪ U.
S \ (T ∩ U ) = (S \ T ) ∪ (S \ U ),
S \ (T ∪ U ) = (S \ T ) ∩ (S \ U ).
Definition 6. Let S and T be sets. We denote by S × T the set of all pairs of the
form (s, t) where s ∈ S and t ∈ T . We call S × T the Cartesian product of the
sets S and T .
E×E = {(0, 0), (0, 1), (0, 2), (1, 0), (1, 1), (1, 2), (2, 0), (2, 1), (2, 2)},
E×F = {(0, 4), (0, 5), (1, 4), (1, 5), (2, 4), (2, 5)},
F ×E = {(4, 0), (5, 0), (4, 1), (5, 1), (4, 2), (5, 2)},
F ×F = {(4, 4), (4, 5), (5, 4), (5, 5)}.
Note that if S and T are two sets, the set S × T is in general distinct from the
set T × S (as you can see on the previous example). Note also that if either S
or T is the empty set then S × T will also be the empty set.
We can use the Cartesian product operation to construct yet more sets from
a given set S as follows: We can construct not just S × S (as we did above),
but also (S × S) × S, and ((S × S) × S) × S, and ... We shall write these last
ones simply as S × S × S, and S × S × S × S, and ... (i.e. by removing the
parentheses). So, just as we had S × S be the set of all pairs (s1 , s2 ) with s1
and s2 elements of S, we will have:
• ···
• R2 , the set of all pairs of real numbers, i.e. the set of all pairs of the form
(a, b) where a, b ∈ R (shorthand for “a and b are both elements of R”),
• R3 , the set of all triples of real numbers, i.e. the set of all triples of the
form (a, b, c) where a, b, c ∈ R,
• and more generally, Rn , the set of all n-tuples of real numbers, i.e. the
set of all n−tuples of the form (a1 , a2 , · · · , an ) with a1 , a2 , · · · , an ∈ R.
Quantifiers
Let S be a set, and suppose all elements of S satisfy some property. We want
to express this fact formally. To make this more concrete, consider the set N of
all integers greater than or equal to 0; it is clear that every element x of N is
such that x + 1 is also an element of N. We can write this statement formally
(and succintly) as:
∀x ∈ N : x + 1 ∈ N.
The symbol ∀ in the expression above is to be read as “for any”, or, equivalently,
as “for all”. The symbol ∀ is called the universal quantifier. The colon symbol
“:” in the above expression should be read as “we have” or “the following holds”.
Using the universal quantifier, many otherwise length statements can be written
formally in a very succint way. For example, if S, T are two sets, we know that
S ⊂ T if and only if every element of S is also an element of T . We can write
this formally as:
S ⊂ T ⇔ ∀x ∈ S : x ∈ T
where again the symbol of equivalence ⇔ is meant to be read as “if and only
if” or “is equivalent to”.
Let now S be a set, and suppose there is at least one element in S which
satisfies some property. We want to express this fact formally. Again, to make
this more concrete, consider the set N of all integers greater than or equal to 0;
it is clear that there are elements in N (namely all multiples of 2) which, after
division by 2, yield again an integer in N. Clearly, this is not true of all elements
in N (consider any odd integer). We can write the statement “there exists at
MAPPINGS AND FUNCTIONS 9
Using the existential quantifier also allows writing expressions formally and
succintly. For example, if S, T are two sets, then we know their intersection
S ∩ T is non-empty if and only if (by definition of intersection), there is an
element in S which is also an element in T . We can write this formally as:
S ∩ T 6= ∅ ⇔ ∃x ∈ S : x ∈ T
Definition 7. Let S and T be two sets. A mapping (also called function) from
S to T is a rule which assigns to each element of S one and only one
element of T .
f :S → T
x 7→ f (x)
The above expression should be read as: “f is a mapping from the set S to the
set T which assigns to every element x of S the element f (x) of T ”.
In the expression above, S is called the domain set of the function f , T is
called the target set of f . Each function is uniquely defined by specifying its
domain set, its target set, and the rule by which it associates to each element
of the domain set a unique element of the target set.
1. The mapping
f :Z → N
x 7 → x2
is the mapping with domain set Z and target set N, defined by assigning
to each integer x its squared value x2 .
10 SECTION 0
2. The mapping
f :R → R
x 7 → sin(x)
is the mapping with domain set and target set both equal to R, defined
by assigning to each real number x its sine sin(x).
3. The mapping
f :N×N → R
x
(x, y) 7→
y2 +1
is the mapping with domain set the Cartesian product N × N (i.e. N2 ) and
target set R which assigns to every pair (x, y) of integers the real number
x
y 2 +1 .
4. The mapping
f : R3 → R2
(x, y, z) 7→ (x + y, x − z)
is the mapping with domain set R3 and target set R2 which assigns to every
triple (x, y, z) of real numbers the pair (x + y, x − z) of real numbers.
It is important to note that given two sets S and T , not every rule of assignement
from S to T defines a function; the rule of assignment will define a function from
S to T if to every element of S it associates a unique element of T . To make
this clear, consider the set N for both the domain set and the target set. Suppose
to each element x of N we assign an element y of N if x is a multiple of y (i.e.
if y divides x). For example, we would assign to the element 2 of N both the
element 1 of N (since 2 is a multiple of 1) and the element 2 of N (since 2 is a
multiple of 2), and so on. Does this rule of assignment define a function from
N to N ? NO! For the simple reason that some elements of the domain set (for
example 2) are assigned to more than one element of the target set.
Consider now the rule of assignment defined by assigning to every element x of
N the element y of N if the relation y = x − 5 is satisfied. Does this rule of
assignment define a function from N to N ? NO! For the simple reason that some
elements of the domain set (namely the integers 0, 1, 2, 3, 4 are not assigned to
any element at all of the target set.
and
g:N → R
x 7→ x3
are considered to be distinct mappings (even though they both map each
integer to its cube) since they do not have the same target set in their definition.
In what follows, we shall assume that all the domain sets and target sets involved
in the discussion are non-empty (since there would be not much to say if they
were!).
g:N → N
x 7→ x + 1
is not.
g◦f :S → U
x 7→ g(f (x))
idS : S → S
x 7→ x
i.e., the mapping which assigns to each element of S itself. idS is called the
identity mapping of S for the obvious reason that it maps each element of S
to itself. It is easy to verify that for any set S, the identity mapping idS of S is
both injective and surjective.
and hence f ◦idS = f , i.e. the two mappings f ◦idS and f are one and the same,
since they have identical domain and target sets, and since they have the same
COMPLEX NUMBERS 13
exact rule of assignement. It is possible to verify in the same exact manner that
the composition idT ◦ f is defined, and that idT ◦ f = f .
Complex Numbers
We will be dealing with complex numbers only at the very end of this course
(when eigenvalues and eigenvectors begin to appear), and since it is assumed
here that you have learned or soon will learn complex numbers in different
settings, we keep the treatment of complex numbers to the bare minimum.
Formally, a complex number can be considered as a pair of real numbers, i.e.
an element of R2 , with addition of pairs and multiplication of pairs defined as
follows:
It can be verified that two operations satisfy the usual properties of addition
and multiplication that are satisfied for real numbers. The set of all complex
numbers is denoted by C.
Informally, you can think of a complex number as a number of the form a + ib,
where a and b are real numbers, and where the symbol i satisfies i2 = −1 (i.e.
i is the square root of −1). Adding and multiplying complex numbers then
proceeds just as with real numbers, with the proviso that i2 = −1.
14 SECTION 0
For example, adding the complex numbers a+ib and c+id (where again a, b, c, d
are real numbers) yields (using the same grouping rules as for real numbers)
and multiplying them yields (using again the same grouping rules as for real
numbers)
For the complex number z = a + ib, where a, b ∈ R, the real number a is called
the real part of z, whereas the real number b is called the imaginary part of
z. Note that the complex number i can be written as
i = 0 + i1,
which shows that i has real part 0 and imaginary part 1. On the other hand, if
a ∈ R is any real number, we can write a as
a = a + i0,
which shows that the real number a has real part a (i.e. is equal to its own real
part) and has imaginary part 0.
As a very simple example of manipulation of real numbers, we all know that for
a, b ∈ R, we have the factorization:
a2 − b2 = (a − b)(a + b).
a2 + b2 = (a − ib)(a + ib).
there are only two complex numbers which satisfy the equation z 2 + 1 = 0,
namely the complex numbers i and −i. Similarly, if we wish to find a complex
number z which satisfies z 2 + 2 = 0, we can write
√ √
z 2 + 2 = (z + i 2)(z − i 2),
COMPLEX NUMBERS 15
√ √
and hence the roots of the equation z 2 + 2 = 0 are given by i 2 and −i 2.
Let now a, b ∈ R, and consider the complex number z = a + ib. The complex
conjugate of z is the complex number denoted by z̄ and defined by z̄ = a − ib.
Note that we have:
z z̄ = (a + ib)(a − ib) = a2 + b2 .
1 1
= = −i,
i 0 + i1
1
whereas 1+i yields
1 1 1 1
= = −i .
1+i 1 + i1 2 2
PROBLEMS:
6. For each of the cases below, determine whether the given rule of assign-
ment defines a valid function with domain set equal to N2 and target set
R; in case it does, determine whether the function is injective, surjective,
bijective.
(a)
f :Z → Z
x 7→ x + 1
and
g:Z → R
x2 − 1
x 7→
x2 + 1
(b)
f :Z → Z
x 7→ x2
and
g:N → R
1
x 7→
x+1
18 SECTION 0
(c)
f :R → R
x 7 → x2
and
g:R → R
1
x 7→
x2 +1
(d)
f :N → Z
x 7 → x2
and
g:Z → N
x 7→ x3
(e)
f :R → R
1
x 7→
x2 +4
and
g:Z → R
1
x 7→
x2 + 2
Section 1
Study Topics
• Systems of Linear Equations
• Number of Solutions to Systems of Linear Equations
19
20 SECTION 1
where we now wish to solve for the real numbers x, y and z, or, equivalently,
for the triple of real numbers (x, y, z); this is a system of linear equations
with two equations (i.e. m = 2) and three unknowns (i.e. n = 3).
where we wish to solve for the unknown real number x. Multiplying both sides
of the equation 2x = 3 by 21 , we obtain the equation:
1 1
2x = 3,
2 2
which, upon simplification, yields:
3
x= .
2
We conclude that the system of linear equations (A) has a unique solution, and
it is given by the real number x = 32 .
Let us still deal with the simplest case of systems of linear equations with
one equation and one unknown (i.e. m = n = 1), and consider now the system:
(B) 0x = 0,
where we again wish to solve for the real number x. It is clear that any real
number x satisfies the equation 0x = 0; we conclude therefore that the system
(B) has infinitely many solutions.
Still remaining with the simplest case of systems of linear equations with one
equation and one unknown (i.e. m = n = 1), consider now instead the system:
(C) 0x = 2,
where we again wish to solve for the real number x. It is clear that there is no
real number x that satisfies the equation 0x = 2; we conclude therefore that
the system (C) has no solution.
We can recapitulate our findings for the three examples that we considered as
follows:
22 SECTION 1
Let us now increase the complexity by one notch and examine systems of
linear equations with 2 equations and 2 unknowns (i.e. m = n = 2). To begin
with, let us consider the following system:
x + y = 2,
(A′ )
x − y = 1,
where we wish to solve for the pair (x, y) of real numbers. If the pair (x, y)
is a solution of the system (A′ ), then we obtain from the first equation that
x = 2 − y and from the second equation that x = 1 + y. Hence, if the pair (x, y)
is a solution of the system (A′ ), then we must have x = 2 − y and x = 1 + y;
this in turn implies that we must have 2 − y = 1 + y, which then implies that
2y = 1, from which we obtain that y = 12 . From the relation x = 2 − y, we
then obtain that x = 2 − 12 = 23 . Note that we could have also used the relation
x = 1 + y instead, and we would still have obtained x = 32 . What we have
therefore established is that if the pair (x, y) is a solution of the system (A′ ),
then we must have x = 32 and y = 12 , i.e., the pair (x, y) must be equal to the
pair ( 23 , 21 ); conversely, it is easy to verify, by simply substituting in values, that
the pair ( 32 , 21 ) is actually a solution to the system (A′ ); indeed, we have:
3 1 4
+ = = 2,
2 2 2
3 1 2
− = = 1,
2 2 2
as expected. We conclude therefore that the system (A′ ) of linear equations has
exactly one solution, and that this solution is given by the pair ( 32 , 12 ).
Consider now the following system of linear equations, again having 2 equa-
tions (i.e. m = 2) and 2 unknowns (i.e. n = 2):
x + y = 1,
(B ′ )
2x + 2y = 2,
where again we wish to solve for the pair (x, y) of real numbers. Applying the
same “method” as previously, we obtain the relation x = 1 − y from the first
equation and the relation 2x = 2 − 2y from the second equation; multiplying
both sides of this last relation by 12 yields, upon simplification x = 1 − y.
Proceeding as before, we use the two relations obtained, namely x = 1 − y and
x = 1 − y, obtaining as a result the equation 1 − y = 1 − y, which then yields
23
0y = 0. Clearly, any real number y satisfies the relation 0y = 0, and for each
such y, the corresponding x is given by x = 1 − y. In other words, if the pair
(x, y) is a solution of the system (B ′ ), then we must have x = 1−y, i.e., we must
have (x, y) = (1 − y, y). Conversely, for any real number y, the pair (1 − y, y)
is a solution of the system (B ′ ), as can be directly verified. Indeed,
(1 − y) + y = 1,
2(1 − y) + 2y = 2 − 2y + 2y = 2,
where again we wish to solve for the pair (x, y) of real numbers. Applying the
same “method” as previously, we obtain the relation x = 1 − y from the first
equation, and the relation 2x = −2y from the second equation; multiplying both
sides of this last relation by 12 yields, upon simplification, the relation x = −y.
We have therefore that if the pair (x, y) is a solution of system (C ′ ), then we
must have x = 1 − y and x = −y, and hence, we must also have as a result
that 1 − y = −y, which yields 0y = 1. Since there exists no real number y such
that 0y = 1, we conclude that there exists no pair (x, y) of real numbers which
would be a solution to system (C ′ ). We conclude therefore that system (C ′ ) has
no solution.
We recapitulate our findings concerning systems (A′ ), (B ′ ) and (C ′ ) as follows:
Note that these are exactly the cases that we encountered for the systems
(A), (B) and (C), which consisted of only one equation in one unknown.
Some very natural questions that come to mind at this point are:
• Can we come up with a system of linear equations that instead has exactly
2 solutions ? or exactly 17 solutions ? or 154 solutions ? ...
24 SECTION 1
We shall see in the next few lectures that these are the only cases that
we can encounter with systems of linear equations; that is, a system of linear
equations, with whatever number of equations and whatever number of un-
knowns, can have either exactly one solution, infinitely many solutions, or
no solution at all; it is instructive to contrast this with the case of polyno-
mial equations of degree 2 in one real variable, i.e. equations of the form
ax2 + bx + c = 0, with a 6= 0, which we know can never have infinitely many
solutions (recall that they can have exactly 0, 1, or 2 solutions).
In order to get to a point where we can prove this non-trivial result, we
shall first develop the necessary tools. This will be done in the next few
lectures.
Later on, we will also see how to quickly and systematically solve systems
of linear equations (as opposed to the ad hoc scheme we have been using) ...
x1 HCl + x2 N a → x3 N aCl + x4 H2 ,
x1 − 2x4 = 0,
x1 − x3 = 0,
x2 − x3 = 0.
A
x1 x5
1 5
3
B D
x3 x4
x2 2
consisting of four vertices (or nodes), labelled A,B,C,D, and edges connect-
ing them. The numeric weight attached to each edge denotes the number
of units flowing through that edge, and we write the so-called “node bal-
ance” equations by equating at each node the total weights of incoming
and outgoing edges. If at a given node, there are no incoming (resp. out-
going) edges, then the total weight of incoming (resp. outgoing) edges
for that node is 0. For the graph depicted above, x1 , x2 , x3 , x4 , x5 are the
unknown weights that have to be determined from the node balance equa-
tions. Here, we have a total of 1+x5 units flowing into node A, and a total
of 5 + x1 units flowing out of node A; equating these yields the equation
1+x5 = 5+x1 , or, equivalently, the equation x1 −x5 = −4. At node B, we
26 SECTION 1
x1 − x5 = −4,
x1 − x2 + x3 = 4,
x2 + x3 − x4 = 2,
x4 + x5 = 10,
c = 1,
a+b+c = 2,
4a + 2b + c = 4,
10 20 20 10
x1 x2
15 18
20 30 30 20
The temperature at each mesh point is indicated next to that mesh point.
There are two points on this mesh which are in the interior of the plate
(i.e. not on its boundary), with respective temperatures x1 and x2 ; we
wish to determine these two temperatures. The fundamental laws of heat
propagation imply that at steady-state, the temperature at each interior
mesh point is the average of the temperatures of its neighboring mesh
points. Applying this law to the two interior mesh points (with respective
temperatures x1 and x2 ) yields:
1
x1 = (15 + 20 + x2 + 30),
4
1
x2 = (x1 + 20 + 18 + 30),
4
which is equivalent to the system of linear equations:
4x1 − x2 = 65,
−x1 + 4x2 = 68,
from which x1 , x2 can be obtained. Note that the temperatures obtained
at the two interior points are merely approximations to the exact temper-
atures; however, the finer the mesh (i.e. the more grid points), the better
the approximation.
PROBLEMS:
For each of the following systems of linear equations, identify the number of
solutions using the same “procedure” as in the examples treated in this lecture
(this should convince you that there can be only the three cases mentioned in
this Lecture):
1.
2x = 3,
4x = 5,
where we wish to solve for the real number x.
2.
2x = 3,
4x = 6,
where we wish to solve for the real number x.
28 SECTION 1
3.
2x + y = 1,
4x + 2y = 3,
10.
x + y − z − w = 0,
−x − y = 2,
z + w = 3,
11. Methane (CH4 ) combines with oxygen gas (O2 ) to yield carbon dioxide
(CO2 ) and water (H2 O); write the chemical balance equations for this
reaction as a system of linear equations.
12. Tin oxide (SnO2 ) combines with hydrogen gas (H2 ) to yield tin (Sn) and
water (H2 O); write the chemical balance equations for this reaction as a
system of linear equations.
13. Iron (F e) combines with sulfuric acid (H2 SO4 ) to yield ferric sulfate
(F e2 (SO4 )3 ) and hydrogen gas (H2 ); write the chemical balance equa-
tions for this reaction as a system of linear equations.
14. Propane (C3 H8 ) combines with oxygen gas (O2 ) to produce water (H2 O)
and carbon dioxide (CO2 ); write the chemical balance equations for this
reaction as a system of linear equations.
15. Write the node balance equations (as a system of linear equations in the
unknowns x1 , x2 , x3 ) for the following graph:
A x1
4
x2
B x3 D E
x1
x2
x3
C
16. Write the node balance equations (as a system of linear equations in the
unknowns x1 , x2 , x3 , x4 , x5 ) for the following graph:
30 SECTION 1
A x1
2 3
3 x5
B D E
x4
x2 7
x3
C
17. Write the node balance equations (as a system of linear equations in the
unknowns x1 , x2 , x3 , x4 , x5 , x6 ) for the following graph:
x6
A F
x1 x5
x4 5x4
3x2
x1 2x3
B D E
7
2x6
x5 2x1
x3
C
18. The voltage V at the output of an electric device is related to the input
current I to the device by the polynomial function V (I) = aI 3 +bI 2 +cI+d,
where the real coefficients a, b, c, d are unknown; we wish to determine
a, b, c, d from experimental data. We experimentally measure the output
voltage at the input current values I = 0, 1, 2, 3, and we determine from
measurement that V (0) = 0.2, V (1) = 2.1, V (2) = 5.7, and V (3) =
10.4. Write the system of linear equations (in the unknowns a, b, c, d)
corresponding to these measurements.
19. The elongation l of a metal bar is modelled as a function of temperature
T by the polynomial function l(T ) = aT 4 + bT 3 + cT 2 + dT + e, where the
real coefficients a, b, c, d, e are unknown; we wish to determine a, b, c, d, e
from experimental data. We experimentally measure the elongation of the
bar at temperatures T = 0, 1, 2, 3, 4, and we determine from measurement
that l(0) = 0.11, l(1) = 2.21, l(2) = 4.35, l(3) = 7.26, l(4) = 9.77. Write
the system of linear equations (in the unknowns a, b, c, d, e) corresponding
to these measurements.
20. Consider a square thin metal plate with temperature at steady state and
31
x1 x2
15 45
20 x3 x4 50
20 30 40 50
The steady-state temperature at each mesh point is indicated next to
that mesh point. The four mesh points interior to the plate have respec-
tive steady-state temperatures x1 , x2 , x3 , x4 ; write the system of linear
equations that governs the relation between these temperatures.
21. Consider a square thin metal plate with temperature at steady state and
with known boundary temperature; we represent this thin metal plate by
the following square mesh:
10 20 30 40 50
x1 x2 x3
15 50
x4 x5 x6
20 55
x7 x8 x9
25 60
30 35 40 50 60
The steady-state temperature at each mesh point is indicated next to that
mesh point. The nine mesh points interior to the plate have respective
steady-state temperatures x1 , x2 , x3 , x4 , x5 , x6 , x7 , x8 , x9 ; write the system
of linear equations that governs the relation between these temperatures.
32 SECTION 1
Section 2
Study Topics
• Real vector spaces
• Examples of real vector spaces
33
34 SECTION 2
The basic notion in linear algebra is that of vector space. In this lecture, we
give the basic definition and review a number of examples.
Definition 12. Let V be a set, with two operations defined on it:
(i) An operation denoted by “+” and called addition, defined formally as a
mapping + : V × V → V which maps a pair (v, w) in V × V to the element
v + w of V;
(ii) An operation denoted by “·” and called multiplication by a scalar, defined
formally as a mapping · : R × V → V which maps a pair (α, v) in R × V
(i.e. α ∈ R and v ∈ V) to the element α · v of V.
With these two operations in place, V is said to be a real vector space if the
following properties are verified:
1. The operation + is associative, i.e. for any x, y, z in V, we have:
x + (y + z) = (x + y) + z
x+0=0+x=x
x + (−x) = (−x) + x = 0
x+y=y+x
α · (β · x) = (αβ) · x
α · (x + y) = α · x + α · y
(α + β) · x = α · x + β · x
1·x=x
35
We shall usually denote the vector space by (V, +, ·) and sometimes only by V,
when there is no risk of confusion about what the addition and multiplication
by scalar operation are. Each element of V is called a vector.
IMPORTANT NOTE: If given some set V and operations “+” and “·” de-
fined in some way, one wants to verify whether or not (V, +, ·) is a real vector
space, before going through the axioms one by one it is important to first
verify that the operations “+” and “·” are both well-defined. To illus-
trate this point with an example, if we let V = N2 and we define as the addition
operation on V the operation “+̃” on V by stipulating that ∀(x1 , y1 ) ∈ N2 and
√
∀(x2 , y2 ) ∈ N2 we define (x1 , y1 )+̃(x2 , y2 ) to be the pair (x1 x2 , y1 y2 ), then
this operation will not be well-defined since for some pairs (x1 , y1 ), (x2 , y2 ) in
N2 the pair (x1 , y1 )+̃(x2 , y2 ) will not be in N2 ; for
√ example, if we take the pairs
(1, 1) and (2, 2), we will have (1, 1)+̃(2, 2) = (2, 2) which is not an element of
N2 . Hence, we can stop right here and declare that since the operation +̃ is not
well-defined, (V, +̃, ·) will not be a real vector space (no matter how we define
the operation “·”).
It is easy to verify that endowed with these two operations, F (R; R) sat-
isfies all the axioms of a real vector space. We can write therefore that
(F (R; R), +, ·) is a real vector space.
(d) Let now F0 (R; R) denote the set of all functions f : R → R which satisfy
f (0) = 0, i.e. the set of all real valued functions of a real variable which
vanish at 0. Formally, we write:
Note that F0 (R; R) is a subset of F (R; R). Note also that we can define
the addition operation “+” on F0 (R; R) the same way we defined it on
F (R; R); indeed, if f, g are two elements of F0 (R; R), then the function
f + g satisfies (f + g)(0) = f (0) + g(0) = 0 and hence f + g is an element
of F0 (R; R) as well. Similarly, if we define the multiplication by scalar
operation “·” on F0 (R; R) the same way we defined it on F (R; R), we
have that if α is in R and f is in F0 (R; R), then the function α · f satisfies
(α · f )(0) = αf (0) = 0, and as a result, the function α · f is an element of
F0 (R; R) as well. Just as easily as with F (R; R), it is immediate to verify
that endowed with these two operations, F0 (R; R) satisfies all the axioms
of a real vector space. We can write therefore that (F0 (R; R), +, ·) is a
real vector space.
(e) Let Rn denote the set of all n−tuples (x1 , x2 , · · · , xn ) of real numbers
(where n is any integer ≥ 1); Rn is defined as the nth Cartesian product
of R with itself, i.e. Rn = R × R × · · · × R (n times). Note that we saw
the special case of this construction for n = 1 in Example (a) above, and
for n = 2 in Example (b) above. We define on Rn the addition operation
“+” as follows: For any n−tuples (x1 , x2 , · · · , xn ) and (y1 , y2 , · · · , yn ) of
real numbers, we define:
It is easy to verify that endowed with these two operations, Rn satisfies all
the axioms of a real vector space. We can write therefore that (Rn , +, ·)
is a real vector space.
37
(f) Let (Rn )0 denote the set of all n−tuples (x1 , x2 , · · · , xn ) of real numbers
which satisfy x1 + x2 + · · · + xn = 0. Note that (Rn )0 is a subset of Rn .
Furthermore, we can define on (Rn )0 the addition operation “+” and the
multiplication by scalar operation “·” in exactly the same way that we
defined them on Rn . Indeed, if (x1 , x2 , · · · , xn ) and (y1 , y2 , · · · , yn ) are
two elements of (Rn )0 , then the n−tuple (x1 + y1 , x2 + y2 , · · · , xn + yn ) is
also in (Rn )0 since
It is easy to verify that endowed with these two operations, (Rn )0 sat-
isfies all the axioms of a real vector space. We can write therefore that
((Rn )0 , +, ·) is a real vector space.
Let us now consider examples of sets with operations on them which do not
make them real vector spaces:
(g) On the set Rn , which we have defined in Example (e) above, define an
addition operation, which we denote by +̃ to distinguish it from the one
defined in Example (e), as follows: For any n−tuples (x1 , x2 , · · · , xn ) and
(y1 , y2 , · · · , yn ) of real numbers, define:
It is easy to verify that endowed with these two operations F (R; R) does
not satisfy all the axioms of a real vector space. We can write therefore
that (F (R; R), +,˜·) is not a real vector space.
38 SECTION 2
Now that the concept of real vector space is hopefully getting more concrete,
let us prove the following simple (and very intuitive) results for general vector
spaces:
Proof. Let us first prove the first statement. Since 0 = 0 + 0, we have, of course:
0 · v = (0 + 0) · v = 0 · v + 0 · v,
where we have used Property (7) of a vector space to get this last equality; now
by property (3) of a vector space, there exists an element in V, which we denote
by −0 · v, such that −0 · v + 0 · v = 0; adding −0 · v to both sides of the above
equality yields:
−0 · v + 0 · v = −0 · v + (0 · v + 0 · v),
−0 · v + 0 · v = (−0 · v + 0 · v) + 0 · v,
0 = 0 + 0 · v,
0 = 0 · v,
i.e. 0 · v = 0.
Let us now prove the second statement. Since 0 = 0 + 0 (Property (2)), we
have:
α · 0 = α · (0 + 0) = α · 0 + α · 0,
where we have used Property (6) to get this last equality; By property (3), there
exists an element in V, which we denote by −α · 0, such that −α · 0 + α · 0 = 0;
adding −α · 0 to both sides of the above equality yields:
−α · 0 + α · 0 = −α · 0 + (α · 0 + α · 0),
−α · 0 + α · 0 = (−α · 0 + α · 0) + α · 0,
39
0 = 0 + α · 0,
0 = α · 0,
i.e., α · 0 = 0.
v + (−v) = −v + v = 0
(where 0 denotes the zero vector of V). All we know is that such an element
does exist in V; we don’t know yet whether for each choice of v in V such an
element −v is unique or not, nor do we know how to compute it. The following
theorem answers these two questions:
Proof. Let us prove (a) first. Let then v ∈ V be any element of V, and assume
the element w ∈ V satisfies the defining conditions for −v, i.e. it satisfies:
w + v = v + w = 0.
We shall show that we necessarily have w = −v; this will prove that for any
v ∈ V, the element −v is uniquely defined. Since
w + v = 0,
(w + v) + (−v) = 0 + (−v),
and using axioms (1) and (2) of a real vector space, we obtain:
w + (v + (−v)) = −v,
w + 0 = −v,
40 SECTION 2
To recap the previous theorem: If you are given an element v in a real vector
space (V, +, ·), you can compute its inverse −v simply by computing (−1) · v
(i.e. by “multiplying the vector v by the real number −1”); as simple as that!
and
x1 αx1
x2 αx2
α· .. = ..
. .
xn αxn
It is easily verified (just as for Rn in Example (e) above) that endowed with
these two operations, R cn is a real vector space.
A WORD ON NOTATION:
PROBLEMS:
2. Show that the set of all integers Z = {· · · , −3, −2, −1, 0, 1, 2, 3, · · ·}, with
addition and multiplication defined as usual is not a real vector space,
and explain precisely why.
a1 x1 + a2 x2 + · · · + an xn = 0.
a1 x1 + a2 x2 + · · · + an xn = 1.
a1 x1 + a2 x2 + · · · + an xn ≥ 0.
Study Topics
• Vector subspaces
• Examples and properties of vector subspaces
45
46 SECTION 3
it follows that the n−tuple (αx1 , αx2 , · · · , αxn ) is in (Rn )0 , i.e. α(x1 , x2 , · · · , xn )
is in (Rn )0 . Hence, Property (iii) is verified as well. We conclude that
(Rn )0 is a vector subspace of the real vector space (Rn , +, ·).
(b) Let (V, +, ·) be a real vector space. Then V is itself a vector subspace of
V.
47
(c) Let (V, +, ·) be a real vector space. Then the set {0} consisting of the
zero element 0 of V alone is a vector subspace of V.
(d) Recall the real vector space (F (R; R), +, ·) of all real-valued functions of
a real variable defined in Lecture 2, as well as the subset F0 (R; R) of
F (R; R) defined in the same lecture (recall F0 (R; R) consists of all real-
valued functions f of a real variable for which f (0) = 0). It is easy to
verify that F0 (R; R) is a vector subspace of F (R; R).
For completeness, let us also consider examples of subsets which are not vector
subspaces:
(e) Recall again the real vector space (Rn , +, ·) of all real n−tuples (x1 , x2 , · · · , xn )
which was defined in Lecture 2; let (Rn )1 denote the subset of Rn consist-
ing of all n−tuples (x1 , x2 , · · · , xn ) of Rn for which x1 + x2 + · · · + xn = 1
(i.e. the entries add up to 1). It is easy to verify that (Rn )1 is not a
vector subspace of (Rn , +, ·).
(f) Recall now again the real vector space (F (R; R), +, ·) defined in Lecture
2, and let F1 (R; R) denote the subset of (F (R; R), +, ·) consisting of all
functions f ∈ (F (R; R), +, ·) for which f (0) = 1. It is easy to verify that
F1 (R; R) is not a vector subspace of (F (R; R), +, ·).
Vector subspaces of a given vector space have the following important property:
Theorem 3. Let (V, +, ·) be a real vector space, and let W1 ⊂ V and W2 ⊂ V
be two vector subspaces of V; then their intersection W1 ∩ W2 is also a vector
subspace of V.
Proof: To prove that W1 ∩ W2 is a vector subspace of V , we have to verify
that W1 ∩ W2 satisfies the three properties of a vector subspace.
(i) We begin by showing the first property, namely that the zero element 0
of V is in W1 ∩ W2 . Since W1 is by assumption a vector subspace of V,
the zero element 0 is in W1 ; similarly, since W2 is by assumption a vector
subspace of V, the zero element 0 is in W2 . Hence, the zero element 0
is in W1 and in W2 , that is, 0 is in the intersection W1 ∩ W2 of W1
and W2 .
(ii) Let us now prove the second property of a vector subspace, namely that
for any x and y in W1 ∩ W2 , x + y is also in W1 ∩ W2 . For this, let us
take x and y in W1 ∩ W2 ; we have to show that x + y is also in W1 ∩ W2 .
Since x is in W1 ∩ W2 , it is in W1 ; similarly, since y is in W1 ∩ W2 , it is
also in W1 . Since by assumption W1 is a vector subspace of V, x + y is
also in W1 . We now repeat this procedure for W2 instead of W1 ; since
x is in W1 ∩ W2 , it is in W2 ; similarly, since y is in W1 ∩ W2 , it is also
in W2 . Since by assumption W2 is a vector subspace of V, x + y is also
in W2 . Hence, we have obtained that x + y is in W1 and in W2 ; hence,
x + y is in W1 ∩ W2 .
48 SECTION 3
(ii) Let us now prove the third and last property of a vector subspace, namely
that for any α in R and any x in W1 ∩ W2 , then αx is in W1 ∩ W2 .
Since x is in W1 ∩ W2 , it is in W1 ; since by assumption W1 is a vector
subspace of V, αx is also in W1 . Similarly, since x is in W1 ∩ W2 , it is
also in W2 ; since by assumption W2 is a vector subspace of V, αx is also
in W2 . We have therefore obtained that αx is in W1 and in W2 ; hence,
αx is in W1 ∩ W2 .
PROBLEMS:
1. Consider the real vector space (R2 , +, ·) defined in Example (b) of Lecture
2. For each of the following subsets of R2 , determine whether or not they
are a vector subspace of R2 :
2. Consider the real vector space (F (R; R), +, ·) defined in Example (c) of
Lecture 2. For each of the following subsets of F (R; R), determine whether
or not they are a vector subspace of F (R; R):
Study Topics
• Linear combinations of vectors
• Linear span of a finite set of vectors
51
52 SECTION 4
We know that in a real vector space, we can add two vectors and multiply
a vector by a real number (aka scalar). A natural question that comes to
mind is: if we have a real vector space V and we take, say, two vectors v1 and
v2 of V, what can we get by doing all these possible operations on these
two vectors, i.e. by taking the set of all vectors of the form αv1 + βv2 (with
α and β being real numbers) ?
We shall examine this question in this lecture, and we shall see at the end of this
lecture how it relates to our original problem, that of understanding systems
of linear equations.
Definition 14. Let (V, +, ·) be a real vector space, and let v1 , · · · , vp be a finite
number of elements of V (with p ≥ 1). The expression
α1 v1 + α2 v2 + · · · + αp vp ,
v = α1 v1 + α2 v2 + · · · + αp vp ,
It is clear that S(v1 ,v2 ,··· ,vp ) is a subset of V; indeed, since v1 , v2 , · · · , vp are in
V , multiplying them by real numbers and adding them up gives us something
still in V, since V is a vector space. So we can write:
But there is more! Indeed, the following result shows that S(v1 ,v2 ,··· ,vp ) is not
just any old subset of V; rather, it is a vector subspace of V:
Proposition 1. Let (V, +, ·) be a real vector space, and let v1 , · · · , vp be a finite
number of elements of V (where p ≥ 1). The subset S(v1 ,v2 ,··· ,vp ) of V consisting
of all linear combinations of the vectors v1 , · · · , vp is a vector subspace of
V.
Proof. To prove this result, we have to show that S(v1 ,v2 ,··· ,vp ) satisfies the
three properties that a vector subspace of V should satisfy. We verify these
properties one by one:
(i) We have to show that the zero vector 0 of V is also an element of
S(v1 ,v2 ,··· ,vp ) ; to do this, we have to show that the zero vector 0 can be
53
and this shows that the zero vector 0 can indeed be expressed as a linear
combination of the vectors v1 , v2 , · · · , vp ; hence the zero vector 0 is an
element of S(v1 ,v2 ,··· ,vp ) , i.e. 0 ∈ S(v1 ,v2 ,··· ,vp ) .
(ii) Let x and y be two elements of S(v1 ,v2 ,··· ,vp ) ; we have to show that the
vector x + y is also in S(v1 ,v2 ,··· ,vp ) . But since by assumption x is in
S(v1 ,v2 ,··· ,vp ) , it must be that
x = α1 v1 + α2 v2 + · · · αp vp
for some real numbers α1 , α2 , · · · , αp , by definition of the set S(v1 ,v2 ,··· ,vp )
itself. Similarly, since by assumption y is in S(v1 ,v2 ,··· ,vp ) , it must be that
y = β1 v1 + β2 v2 + · · · βp vp
for some real numbers β1 , β2 , · · · , βp , by definition of S(v1 ,v2 ,··· ,vp ) . Hence,
the sum x + y of x and y can be written as:
x = α1 v1 + α2 v2 + · · · αp vp
for some real numbers α1 , α2 , · · · , αp , by definition of the set S(v1 ,v2 ,··· ,vp )
itself. Hence, we can write:
γx = γ(α1 v1 + α2 v2 + · · · αp vp )
= (γα1 )v1 + · · · (γαp )vp ,
Before going any further, we give a name to S(v1 ,v2 ,··· ,vp ) :
Definition 15. The vector subspace S(v1 ,v2 ,··· ,vp ) of V is called the subspace of
V generated by the vectors v1 , v2 , · · · , vp ; it is also called the linear span
of the vectors v1 , v2 , · · · , vp .
Now that S(v1 ,v2 ,··· ,vp ) has an honest name (actually two!), we examine a num-
ber of examples:
(a) Let (V, +, ·) be a real vector space, and let v1 = 0 be the zero vector
of V. The vector subspace of V generated by v1 is clearly seen to be the
zero subspace {0} of V.
(b) Consider the real vector space (R2 , +, ·) defined previously, and let v1 =
(1, 0) and v2 = (0, 1); it is easy to verify that the subspace of (R2 , +, ·)
generated by v1 , v2 is V itself.
(c) Consider the real vector space (R3 , +, ·) defined previously, and consider
the two elements v1 = (1, −1, 0) and v2 = (0, 1, −1) of R3 . Let us show
that the vector subspace S(v1 ,v2 ) of R3 spanned by v1 and v2 is the set of
all (a, b, c) in R3 such that a + b + c = 0, that is, we wish to show that
To show this equality, we will show first that we have the inclusion
these two inclusions will, together, show the desired equality. Let us then
begin by showing the first inclusion. For this, let us take an arbitrary
element x in S(v1 ,v2 ) , and show that it is also in the set {(a, b, c) ∈ R3 |a +
b + c = 0}. By definition of S(v1 ,v2 ) , x can be written as:
x = αv1 + βv2 ,
x = αv1 + βv2
= α(1, −1, 0) + β(0, 1, −1)
= (α, −α, 0) + (0, β, −β)
= (α, −α + β, −β);
where we wish to solve for the pair (x1 , x2 ) of real numbers. Let us now
recall c3 , +, ·) of all column vectors of the form
the real vector space (R
a
b where a, b, c are real numbers. Consider the elements v1 , v2 , w of
c
c 3
R defined as follows:
1 −1 1
v1 = 2 , v2 = 1 , w = 0 ,
1 −2 2
Hence, the pair (x1 , x2 ) of real numbers satisfies the system (A) of linear
equations above if and only if we have:
x1 · v1 + x2 · v2 = w;
where we wish to solve for the pair (x1 , x2 ) of real numbers. In order
to understand whether this system has a solution or not, consider the
c2 :
following vectors in R
1 −1 1
v1 = , v2 = ,w = .
−1 1 1
It follows that the pair (x1 , x2 ) is a solution of system (B) if and only if
we have:
x1 · v1 + x2 · v2 = w;
So we now have some basic “picture” of what it means for a system of linear
equations to have, or not to have, a solution: If some given vector happens to lie
in some given subspace, then the system has a solution (at least one); otherwise,
it has no solution.
In order to get to a point where we can explain why that number of solutions
is always 0,1, or ∞, we have to further sharpen our tools. This is what we will
do in the following sections.
57
PROBLEMS:
1. Consider the (by now familiar) real vector space R2 , consisting of all pairs
(x, y) of real numbers, with the usual addition and scalar multiplication
operations. Consider the vectors v1 , v2 , v3 , v4 , v5 in R2 defined as follows:
v1 = (0, 0), v2 = (1, 0), v3 = (−1, 0), v4 = (0, 3), v5 = (2, 1).
For each of these four vectors, determine whether they are in the linear
span of the other three.
3. Consider the (familiar) real vector space R3 consisting of all triples of the
form (x, y, z) with x, y, z real numbers, with the usual addition and scalar
multiplication operations. For each of the following list of vectors in R3 ,
determine whether the first vector is in the linear span of the last two:
(a) (1, 1, 1), (1, 2, 1), (1, 3, 1)
(b) (0, 0, 0), (1, 2, 1), (1, 3, 1)
58 SECTION 4
For each of these four elements (i.e. vectors) in F (R; R), verify whether
or not they are in the linear span of the other three.
Section 5
Study Topics
• Linear dependence of a set of vectors
• Linear independence of a set of vectors
59
60 SECTION 5
Consider the vector space R2 of all pairs (x, y) of real numbers, which we have
defined previously and used numerous times. Consider the following subsets of
R2 :
S1 = {(3, 2), (12, 8)}, S2 = {(3, 2), (0, 1)}.
What are the similarities and differences between S1 and S2 ? Well, they both
contain two vectors each; in that, they are similar; another similarity is that
they both contain the vector (3, 2). What about differences ? Notice that the
second vector in S1 , namely (12, 8), is a scalar multiple of (3, 2), i.e. (12, 8)
can be obtained by multiplying (3, 2) by some real number (in this case, the real
number 4); indeed, we have:
4(3, 2) = (4 × 3, 4 × 2) = (12, 8).
We can also obviously claim that (3, 2) is a scalar multiple of (12, 8) (this time
with a factor of 41 ), since we can write:
1 1 1
(12, 8) = ( × 12, × 8) = (3, 2).
4 4 4
On the other hand, no such thing is happening with S2 ; indeed, (0, 1) is not
a scalar multiple of (3, 2), since mutiplying (3, 2) by a real number can never
yield (0, 1) (prove it!). Similary, multiplying (0, 1) by a real number can never
yield (3, 2) (prove it also!).
So we have identified a fundamental difference between S1 and S2 : Whereas the
two elements of S1 are related, in the sense that they are scalar multiples of
one another, those of S2 are completely unrelated (in the sense that they are
not scalar multiples of each other).
Before going further, we generalize this observation into a definition:
Definition 16. Let (V, +, ·) be a real vector space, and let S = {v1 , v2 , · · · , vp }
be a finite subset of V.
(i) The subset S is said to be linearly independent if for any α1 , · · · , αp ∈ R,
the relation
α1 v1 + α2 v2 + · · · + αp vp = 0
implies that α1 = α2 = · · · = αp = 0.
(ii) The subset S is said to be linearly dependent if it is not linearly in-
dependent.
Equivalent ways to re-state these definitions are as follows:
(i) The subset S is linearly independent if the only real numbers α1 , α2 , · · · , αp
which yield
α1 v1 + α2 v2 + · · · + αp vp = 0
are given by α1 = α2 = · · · = αp = 0.
61
Before going any further, let us show that the set S1 above is a linearly de-
pendent subset of R2 , and that the set S2 above is a linearly independent
subset of R2 . Let us begin with S1 : Let α1 = −4 and α2 = 1; we have:
(recall that the zero element 0 of R2 is the pair (0, 0)). Hence we have found real
numbers α1 , α2 which are not both zero but such that the linear combination
α1 (3, 2) + α2 (12, 8) is the zero vector 0 = (0, 0) of R2 . This proves that S1 is a
linearly dependent subset of R2 .
Let us now show that the set S2 above is a linearly independent subset of
R2 . To do this, we have to show that if for some real numbers α1 , α2 we have
c3 : Let then
Let us show that S is a linearly independent subset of R
α1 , α2 , α3 be three real numbers such that
1 0 0 0
α1 0 + α2 1 + α3 0 = 0 ;
0 0 1 0
α1 f1 + α2 f2 = 0;
α1 f1 (t) + α2 f2 (t) = 0,
64 SECTION 5
for all t in R; in particular, this relation should hold for t = 1 and t = −1.
But for t = 1, we have:
α1 f1 (1) + α2 f2 (1) = α2 ,
α1 f1 (−1) + α2 f2 (−1) = α1 .
Consider now a real vector space (V, +, ·), and let S and T be two finite subsets
of V such that S ⊂ T (i.e. S is itself a subset of T ); suppose we know that S
is a linearly dependent subset of V; what can we then say about T ? The
following lemma does answer this question.
Lemma 4. Let (V, +, ·) be a real vector space, and S, T two finite subsets of V
such that S ⊂ T . If S is linearly dependent, then T is also linearly dependent.
Proof. Assume S has p elements and T has q elements (necessarily, we have
q ≥ p since S is a subset of T ). Let us then write:
S = {v1 , v2 , · · · , vp },
T = {v1 , v2 , · · · , vp , vp+1 , vp+2 , · · · , vq }.
α1 v1 + α2 v2 + · · · + αp vp
α1 v1 + α2 v2 + · · · + αp vp = 0.
α1 v1 + α2 v2 + · · · + αp vp + 0vp+1 + · · · + 0vq
is also equal to 0, and yet, not all of the coefficients in the above linear com-
bination are zero (since not all of the α1 , α2 , · · · , αp are zero); this proves that
T is a linearly dependent subset of V.
65
Remark 2. Let (V, +, ·) be a real vector space. We have seen in Example (c) above
that the subset {0} of V consisting of the zero vector alone is linearly dependent;
it follows from the previous lemma that if S is any finite subset of V such that
0 ∈ S, then S is linearly dependent.
Remark 3. The previous lemma also implies the following: If a finite subset S of a
vector space V is linearly independent, then any subset of S is linearly independent
as well.
Remark 4. Let (V, +, ·) be a real vector space, and let {v1 , · · · , vp } be a finite
subset of V. We will sometimes write “v1 , · · · , vp are linearly independent (resp.
dependent) elements of V” instead of “{v1 , · · · , vp } is a linearly independent (resp.
dependent) subset of V”; the two statements in quotes are meant to say exactly the
same thing. One reason we will occasionally prefer to express linear dependence or
independence of v1 , · · · , vp by writing a statement such as “v1 , · · · , vp are linearly
independent (resp. dependent)” is because this does not require v1 , · · · , vp to
be all distinct; if, however, we were to write “the subset {v1 , · · · , vp } is linearly
independent (resp. dependent)”, this would assume that v1 , · · · , vp are all distinct.
To give a concrete example, let v1 ∈ V be any element of V. Then we can say that
the vectors v1 and v1 are linearly dependent (since, for example, (1)·v1 +(−1)·v1 =
0V , i.e. there is a linear combination of the vectors v1 and v1 that gives the
zero vector of V but in which at least one coefficient is non-zero); on the other
hand, the subset of V consisting of the vectors v1 and v1 is the subset {v1 } (not
{v1 , v1 } Remember that repetitions are not allowed when writing sets by listing
their elements), so we would have no way of discussing the linear dependence or
independence of the vectors v1 and v1 by constructing a subset of V from them.
Remark 5. Continuing the previous remark, it should be pointed out that when
we write a statement such as “v1 , · · · , vp are linearly independent (resp. depen-
dent) elements of V” we mean that the “collection” or “family” consisting of the
vectors v1 , · · · , vp is linearly independent (resp. dependent), and not that each of
v1 , · · · , vp is separately linearly independent (resp. dependent). In other words, the
usage of the plural here is not meant in the usual sense it is commonly employed.
This is because linear dependence/independence is a collective property of all the
vectors involved. To give an example, a statement such as “the numbers 2, 4, 6 are
even” means that 2 is even, 4 is even, and 6 is even, whereas a statement such as
“v1 , v2 , v3 are linearly independent” means that the collection or family consisting
of the vectors v1 , v2 , v3 is linearly independent, i.e. for any α1 , α2 , α3 ∈ R such
that α1 · v1 + α2 · v2 + α3 · v3 = 0V it necessarily follows that α1 = α2 = α3 = 0.
PROBLEMS:
v1 = (1, 0), v2 = (3, 0), v3 = (0, 0), v4 = (1, 1), v5 = (2, 1).
66 SECTION 5
For each of the following finite subsets of R2 , specify whether they are
linearly dependent or linearly independent:
(a) {v1 }
(b) {v1 , v2 }
(c) {v2 }
(d) {v3 }
(e) {v1 , v4 }
(f) {v2 , v4 }
(g) {v3 , v4 }
(h) {v4 , v5 }
(j) {v1 , v4 , v5 }
(k) {v2 , v4 , v5 }
(l) {v1 , v2 , v4 , v5 }
(m) {v1 , v5 }
For each of the following finite subsets of R2 , specify whether they are
linearly dependent or linearly independent:
(a) {v1 }
(b) {v1 , v2 }
(c) {v2 }
(d) {v3 }
(e) {v1 , v4 }
(f) {v2 , v4 }
(g) {v3 , v4 }
(h) {v4 , v5 }
(j) {v1 , v4 , v5 }
(k) {v2 , v4 , v5 }
(l) {v1 , v2 , v4 , v5 }
(m) {v1 , v5 }
(n) {v1 , v2 , v3 , v4 , v5 }
(o) {v3 , v4 , v5 }
67
(p) {v2 , v3 , v4 , v5 }
3. Consider the real vector space (F (R; R), +, ·) defined in Lecture 2, and
consider the elements f1 , f2 , f3 , f4 in F (R; R), defined as follows:
f1 (t) = 1, ∀t ∈ R,
f2 (t) = t, ∀t ∈ R,
f3 (t) = (t + 1)2 , ∀t ∈ R,
f4 (t) = t2 + 5, ∀t ∈ R,
f5 (t) = t3 , ∀t ∈ R,
f6 (t) = t3 + 2t2 , ∀t ∈ R,
For each of the following finite subsets of F (R; R), specify whether they
are linearly dependent or linearly independent:
(a) {f1 }
(b) {f1 , f2 }
(c) {f1 , f2 , f3 }
(d) {f1 , f2 , f3 , f4 }
(e) {f2 , f3 , f4 }
(f) {f4 , f5 , f6 }
(g) {f3 f4 , f5 , f6 }
(h) {f2 , f3 , f4 , f5 , f6 }
(j) {f1 , f2 , f3 , f4 , f5 , f6 }
4. Let (V, +, ·) be a real vector space, and let S be a finite subset of V. Show
directly (without using Lemma 1 or Remark 2 of these lectures) that if 0
is in S then S is a linearly dependent subset of V .
68 SECTION 5
Section 6
Study Topics
• Linear dependence/independence and linear combinations
• Application to systems of linear equations (or, we now have
the complete answer to our initial question!)
69
70 SECTION 6
α1 v1 + α2 v2 + · · · αp vp = v;
A very natural question at this point is: In how many different ways can
this same v be expressed as a linear combination of v1 , v2 , · · · , vp ? Does there
exist another p−tuple of real numbers, say (γ1 , γ2 , · · · , γp ), distinct from the
p−tuple (α1 , α2 , · · · , αp ) such that we also have
γ1 v1 + γ2 v2 + · · · γp vp = v ?
α1 v1 + α2 v2 + · · · αp vp = v;
α1 v1 + α2 v2 + · · · αp vp = v
and
γ1 v1 + γ2 v2 + · · · γp vp = v.
α1 v1 + α2 v2 + · · · αp vp = γ1 v1 + γ2 v2 + · · · γp vp ,
which then implies (by putting everything on the left-hand side of the “=” sign):
α1 = γ1 ,
α2 = γ2 ,
···
αp = γp ,
α1 v1 + α2 v2 + · · · αp vp = v
and
γ1 v1 + γ2 v2 + · · · γp vp = v.
α1 v1 + α2 v2 + · · · αp vp = v;
α1 v1 + α2 v2 + · · · + αp vp = 0;
to show that S is linearly independent, we have to show that this last equality
implies that α1 , α2 , · · · , αp are all zero. Note first that 0 is in the linear span of
v1 , v2 , · · · , vp ; indeed, we can write:
Since we have assumed that any element in the linear span of v1 , v2 , · · · , vp can
be written only in a unique way as a linear combination of v1 , v2 , · · · , vp , and
since we have written 0 as the following linear combinations of v1 , v2 , · · · , vp :
0 = α1 v1 + α2 v2 + · · · + αp vp ,
0 = 0v1 + 0v2 + · · · + 0vp ,
72 SECTION 6
it follows that these two linear combinations must be one and the same, i.e.
their respective coefficients must be equal, i.e., we must have:
α1 = 0,
α2 = 0,
···
αp = 0.
α1 v1 + α2 v2 + · · · + αp vp = 0
α1 v1 + α2 v2 + · · · αp vp = v;
α1 v1 + α2 v2 + · · · αp vp = v.
λ1 v1 + λ2 v2 + · · · λp vp = v.
β1 v1 + β2 v2 + · · · + βp vp = 0.
Consider now for each real number µ the p−tuple of real numbers given by:
yields:
v = α1 v1 + α2 v2 + · · · αp vp ,
v = β1 v1 + β2 v2 + · · · βp vp ;
and since the two n−tuples (α1 , α2 , · · · , αp ) and (β1 , β2 , · · · , βp ) are assumed
distinct, there should be some integer i in {1, 2, · · · , p} such that αi 6= βi , i.e.
such that (αi − βi ) 6= 0; but this implies that the above linear combination of
v1 , v2 , · · · , vp is equal to the zero vector 0 but does not have all its coefficients
zero; as a result, S is linearly dependent.
74 SECTION 6
• why is it that a system of linear equations can have only 0,1, or infinitely
many solutions ?
where the aij and the bk are given real numbers. We would like to know how
many distinct n−tuples (x1 , x2 , · · · , xn ) of real numbers are a solution to that
system of equations (i.e. satisfy all the equalities above).
In order to answer this question using the tools we have developed, consider the
real vector space (R d m , +, ·) which we are by now familiar with, and consider the
a11 a12 a1n
a21 a22 a2n
v1 = 31 a , v 2 = a
32
, · · · , vn = a3n ,
··· ··· ···
am1 am2 amn
d
as well as the vector w in R m defined as:
b1
b2
w= b3 .
···
bm
x1 v1 + x2 v2 + · · · + xn vn = w.
PROBLEMS:
1. For each of the following systems of linear equations, study the number of
solutions by formulating and analyzing the corresponding linear algebra
problem, as done in the last section of this Lecture under the heading
“Application to Systems of Linear Equations”.
(a)
2x = 4,
4x = 5,
(b)
2x = 4,
4x = 8,
(c)
0x = 1,
(d)
0x = 0,
(e)
2x + y = 1,
10x + 5y = 5,
x − y = 1,
(f)
x + y − z = 0,
x + y = 1,
y + z = 2,
76 SECTION 6
(g)
x + y − z − w = 0,
x + y = 2,
(h)
x + y − z − w = 0,
−x − y = 2,
z + w = 3,
(i)
2x + y = 1,
4x + 2y = 2,
(j)
2x + y = 1,
4x + 2y = 3,
(k)
2x + y = 1,
x = 1,
(l)
2x + y = 1,
y = 0,
(m)
x + y + z = 1,
(n)
x + y + z = 1,
x=0
(o)
x + y + z = 1,
x=0
y=0
z=0
77
(p)
x + y + z = 1,
x=0
y=0
(q)
x + y + z = 1,
2x + 2y + 2z = 2
y=0
(r)
x + y + z = 1,
2x + 2y + 2z = 2
x=0
y=0
(s)
x + y + z = 1,
2x + 2y + 2z = 2
x=0
y=0
z=0
78 SECTION 6
Section 7
Study Topics
• Systems of linear equations in upper triangular form
• Solution by back-substitution
• Gaussian Elimination to put systems of linear equations in up-
per triangular form
• Augmented matrix of a system of linear equations
• Elementary Row operations
79
80 SECTION 7
We now take a little pause from theory and examine a quick and efficient
approach to solving systems of linear equations. Recall that in Section 1,
we solved systems of linear equations using the ad hoc (and quite primitive!)
procedure of isolating one variable from one of the equations, substituting its
expression in terms of the other variables in all other equations, and repeating
the procedure. Such a scheme could be suitable for small systems of linear
equations (i.e. having few equations and few unknowns), but quickly becomes
unwieldy for larger systems of linear equations. In this section, we will learn a
systematic approach to solving such systems. The approach we will learn is often
called “Gaussian Elimination and Substitution”. To motivate the approach, we
begin with a few examples.
Consider first the system of linear equations given by
2x4 = 6,
x1 + 3x2 − 4x3 + 6 = 5,
x2 + 3x3 − 6 = 1,
x3 + 12 = 2,
Note that this system of linear equations now involves only the variables x1 , x2 , x3
(since we already solved for x4 ), and is also in upper-triangular (aka row-echelon)
form. We begin again with the last row. The equation on the last row is given
by
x3 = −10,
which yields directly that x3 = −10. We now substitute this value of x3 that
we just obtained in all other equations, and we obtain:
x1 + 3x2 + 40 = −1,
x2 − 30 = 7,
x1 + 3x2 = −41,
x2 = 37.
Note that this system of linear equations now involves only the variables x1 , x2
(since we already solved for x3 and x4 ), and is also in upper-triangular form.
We begin again with the last row. The equation on the last row is given by
x2 = 37,
which yields directly that x2 = 37. We now substitute this value of x2 that
we just obtained into the first equation (since it is the only equation left!), and
we obtain:
x1 + 111 = −41,
x1 = −152.
We have therefore solved our original system, and we have found that the
unique solution to our system is given by (x1 , x2 , x3 , x4 ) = (−152, 37, −10, 3).
The important thing to note is how systematic and painless the whole procedure
was.
We now examine yet another example. Consider the system of linear equations
given by:
x3 = x4 + 6,
82 SECTION 7
Let us multiply the first row by −2 and add the result to the second row;
denoting row 1 by R1, row 2 by R2, and row 3 by R3, the operation we are
describing consists of replacing R2 by −2R1+R2. We shall denote this operation
as follows:
−2R1 + R2 → R2.
2x1 + x2 + x3 = 5,
0x1 − x2 + x3 = −8,
6x1 + x2 + 4x3 = 0.
To see that this system has exactly the same solution (or solutions) as our
original system, note that we can go back to our original system from this
system by replacing the second row with 2 times row 1 + row 2. Whenever two
systems of linear equations have exactly the same solutions, we say that they
are equivalent. By the row operation described by
−2R1 + R2 → R2,
we have therefore gone from our original system of linear equations to an equiv-
alent system. Consider now the (equivalent) system we have obtained:
2x1 + x2 + x3 = 5,
0x1 − x2 + x3 = −8,
6x1 + x2 + 4x3 = 0.
Let us multiply the first row by −3 and add the result to the third row; the
operation we are describing consists of replacing R3 by −3R1 + R3, and can be
denoted as follows:
−3R1 + R3 → R3.
2x1 + x2 + x3 = 5,
0x1 − x2 + x3 = −8,
0x1 − 2x2 + x3 = −15,
which we can also rewrite (by eliminating the terms having coefficient 0) as:
2x1 + x2 + x3 = 5,
−x2 + x3 = −8,
−2x2 + x3 = −15.
Note again that this system is equivalent to the previous one, and therefore,
to our original system. The solutions to this last system are exactly those of
84 SECTION 7
our original system, and vice-versa. Note also that this system is not upper-
triangular, but we’re almost there! To get there, let us now multiply the
second row by −2 and add the result to the third row; the operation we are
describing consists of replacing R3 by −2R2 + R3, and we shall denote it as
follows:
−2R2 + R3 → R3.
After this operation, our system becomes:
2x1 + x2 + x3 = 5,
−x2 + x3 = −8,
−x3 = 1.
Note again that this system is equivalent to the previous one, and therefore to
our original system, and that this last system is now indeed in upper-triangular
form; we can now easily solve this system by back-substitution, and we obtain
that the unique solution of this system, and hence of our original system, is
given by the triple (x1 , x2 , x3 ) = (−1/2, 7, −1).
Let us now describe a systematic approach that captures what we have done;
this process is known as Gaussian Elimination. Before we get there however,
we need some more terminology.
Definition 17. Consider the system of linear equations in m equations and n un-
knowns given by
a1,1 x1 + a1,2 x2 + · · · a1,n xn = b1 ,
a2,1 x1 + a2,2 x2 + · · · a2,n xn = b2 ,
···
am,1 x1 + am,2 x2 + · · · am,n xn = bm .
The augmented matrix of this system is the table of real numbers with m rows
and n + 1 columns given by
a1.1 a1,2 ··· a1,n | b1
a2.1 a2,2 ··· a2,n | b2
.
···
am.1 am,2 ··· am,n | bm
(Note: The vertical bars “|” in this table are meant to visually separate the coef-
ficients that multiply the unknowns from those that appear on the right-hand side
of the equation).
For example, the augmented matrix corresponding to the system of linear equa-
tions
2x1 + x2 + x3 = 5,
4x1 + x2 + 3x3 = 2,
6x1 + x2 + 4x3 = 0.
85
Hence, we can represent any system of linear equations by its augmented matrix,
and vice-versa.
We now define what it means for a system of linear equations – or, equivalently,
its augmented matrix, to be in upper-triangular, aka row-echelon form:
1. Each row with all entries equal to 0 is below every row having at least one
nonzero entry,
2. the leftmost non-zero entry on each row is to the right of the leftmost non-zero
entry of the preceding row.
x2 + x3 = 1,
2x1 + x2 − x3 = 0,
4x1 + x2 + 2x3 = 5.
Step 1: The augmented matrix of this system is the 3 × 4 table given by:
0 1 1 | 1
2 1 −1 | 0 .
4 1 2 | 5
Let us now multiply row 1 by −2 and add the result to row 3 (we denote
this by −2R1 + R3 → R3); we obtain the augmented matrix
2 1 −1 | 0
0 1 1 | 1 .
0 −1 4 | 5
2x1 + x2 − x3 = 0,
x2 + x3 = 1,
5x3 = 6.
This augmented matrix is not in row-echelon form yet; let us now multiply
row 1 by −1 and add the result to row 3 (we denote this by −R1 + R3 →
R3); we obtain the augmented matrix
1 0 1 −1 | 2
0 0 −2 4 | −4 .
0 0 −1 0 | −1
This augmented matrix is still not in row-echelon form yet; let us now
multiply row 2 by −1/2 and add the result to row 3 (we denote this by
− 12 R2 + R3 → R3); we obtain the augmented matrix
1 0 1 −1 | 2
0 0 −2 4 | −4 .
0 0 0 −2 | 1
x1 + x3 − x4 = 2,
−2x3 + 4x4 = −4,
−2x4 = 1.
Solving for x3 from the last equation yields x3 = 1, and substituting this
value of x3 in the first equation yields
3
x1 + 1 = ,
2
i.e.
1
x1 = .
2
PROBLEMS:
1. Solve each of the following systems of linear equations using the Gaussian
Elimination and Back-Substitution method described in this Section, i.e.
for each system:
(a)
x1 + x2 + x3 = 0,
2x1 + x3 = 1,
x1 − 3x3 = 2.
91
(b)
x1 + x2 + x3 = 0,
x1 + x3 = 1,
x1 − 3x3 = 2,
x1 + 2x3 = 1.
(c)
x1 + x2 + x3 = 0,
x1 + 2x3 = 1,
(d)
x1 + x2 + x3 = 0,
x1 + x2 = −1,
(e)
x1 + x2 + x4 = 0,
x1 + x3 = 1,
x1 − x4 = 2,
x1 + x2 = 1.
(f)
2x3 + x4 = 2,
x1 + x2 = 0,
x1 − x3 = 1,
x3 + x4 = 0.
(g)
2x3 + x4 = 2,
x1 + x2 = 0,
x1 − x3 = 1,
x3 + x4 = 0,
x1 + x2 + x3 + x4 = 1.
(h)
x2 + x3 = 1,
x1 − x4 = 1,
92 SECTION 7
x1 + x2 + x3 = 1 + a,
2x1 + 2x2 + 4x3 = 2,
3x1 + 3x2 + 2x3 = a,
x1 + x2 + x3 = a,
x1 + x3 = 2a,
4x1 + 4x3 = a,
x1 + x2 + 2x3 + 2x4 = a,
x1 + x2 + 2x4 = a,
2x1 + 4x4 = 1,
−x1 + x2 − 2x4 = 2.
5. For each of the following chemical reactions, solve the corresponding chem-
ical balance equations using Gaussian elimination and back-substitution:
(b) SnO2 + H2 Sn + H2 O
(d) C3 H8 + O2 H2 O + CO2
6. Solve the node balance equations (as a system of linear equations in the
unknowns x1 , x2 , x3 ) for the following graph:
93
A x1
−2
x2
B x3 D E
7
2x2
x3
C
7. Solve the node balance equations (as a system of linear equations in the
unknowns x1 , x2 , x3 , x4 ) for the following graph:
A
3x1
2 3
3x4 2x3
B D E
7
2x2 7x1
x3
C
8. Solve the node balance equations (as a system of linear equations in the
unknowns x1 , x2 , x3 , x4 , x5 , x6 ) for the following graph:
x6
A F
x1 x5
x4 5x4
3x2
x1 2x3
B D E
7
2x6
x5 2x1
x3
C
x1 x2
15 15
15 x3 x4 15
20 30 30 20
The steady-state temperature at each mesh point is indicated next to that
mesh point. The four mesh points interior to the plate have respective
steady-state temperatures x1 , x2 , x3 , x4 ; solve the system of linear equa-
tions (in the unknowns x1 , x2 , x3 , x4 ) that governs the relations between
these temperatures.
Section 8
Study Topics
• Generating set for a vector space
• Basis of a vector space
• Coordinates of a vector with respect to a basis
95
96 SECTION 8
α1 v1 + α2 v2 + · · · + αp vp
Definition 20. Let (V, +, ·) be a real vector space, and let S = {v1 , v2 , · · · , vp }
be a finite subset of V. S is said to be a generating set for the vector space V
if any vector v in V can be written as a linear combination of v1 , v2 , · · · , vp ; that
is, for any v in V, we can find real numbers α1 , α2 , · · · , αp such that:
v = α1 v1 + α2 v2 + · · · + αp vp .
It is easy to verify that none of the sets {v1 , v2 }, {v2 , v3 }, and {v1 , v3 }
is a generating set for R c3 , but that the set {v1 , v2 , v3 } is a generating
c 3
set for R . To be sure, let us do some of these verifications. Let us first
prove that {v1 , v2 , v3 } is a generating set for Rc3 : To do this we have to
c3
show that any element of R can be written as a linear combination of
v1 , v2 , v3 . Let then v be any element in R c3 ; by definition of Rc3 , v is a
real column vector, with three entries, i.e. is of the form
a
v = b ,
c
for some real numbers a, b and c; but then, it is clear that we have the
equality
equality would imply that 1 = 0). This shows that there is no linear
combination of v1 and v2 that equals v; hence, {v1 , v2 } is not a generating
c3 .
set for R
(b) Consider now the (equally familiar) real vector space Rc2 of all real column
vectors with two entries. Consider the subset S = {v1 , v2 , v3 } of R c2 ,
c2 are defined as:
where the vectors v1 , v2 , v3 in R
1 0 −1
v1 = , v2 = , v3 = .
0 1 0
c2 . To do this, we
Let us show that S = {v1 , v2 , v3 } is a generating set for R
have to show that any element c2
in R can be written as a linear combination
a
of v1 , v2 , v3 . Let then (with a, b real numbers) be any element of
b
c2 ; it is clear that we can write:
R
a
= av1 + bv2 + 0v3 ,
b
a c2 can be written as a linear combination
that is, any element of R
b
c2 .
of v1 , v2 , v3 . This shows that S is indeed a generating set for R
a
Let us note now that can also be expressed as the following linear
b
combination of v1 , v2 , v3 :
a
= 0v1 + bv2 − av3 ,
b
and so on ... (you get the idea). This shows that any element in Rc2 can
be written as many distinct linear combinations of v1 , v2 , v3 .
The case where a given generating set of a vector space is such that any vector
in that vector space can be expressed as a linear combination of vectors in the
generating set in a unique way is of particular interest, since then there is no
ambiguity as to how the vector should be written as a linear combination of the
99
elements of the generating set. Note that this was the case with Example (a)
above (prove it!), but not with Example (b).
A very natural question at this point is therefore the following: Suppose we
have a vector space V, and we are given a finite subset S of V which is also
a generating set for V. We then know of course that any element of V can
be written as some linear combination of the elements of S; how then, can we
determine whether that is the only possible linear combination that yields the
desired element ? We already have the answer to that question from Theorem
4 of Lecture 6; we can therefore state:
Theorem 11. Let (V, +, ·) be a real vector space, and let S = {v1 , · · · , vp } be a
finite subset of V such that:
1. Consider the real vector space (R2 , +, ·), and consider the elements (0, 1),(1, 1),(2, 2),(1, 0)
of R2 . Consider the following few tuples of vectors made from these 4 vec-
100 SECTION 8
tors:
B1 = ((0, 1)), B2 = ((0, 1), (1, 0)),
B3 = ((1, 0), (0, 1)), B4 = ((0, 1), (1, 1)),
B5 = ((0, 1), (2, 2)), B6 = ((1, 1), (2, 2)),
B7 = ((1, 1), (1, 0), (2, 2)), B8 = ((1, 1), (1, 0)),
B9 = ((1, 1), (0, 1)), B10 = ((1, 1), (1, 0), (2, 2)),
B11 = ((2, 2), (1, 0)), B12 = ((1, 0), (2, 2)),
B13 = ((0, 1), (2, 2)), B14 = ((0, 1), (2, 2), (1, 1)),
B15 = ((2, 2), (0, 1)), B16 = ((2, 2)),
B17 = ((2, 2), (1, 1)), B18 = ((2, 2), (0, 1), (1, 1), (0, 1)),
B19 = ((2, 2), (2, 2), (2, 2), (2, 2)), B20 = ((0, 1), (2, 2), (2, 2), (2, 2));
It is easy to verify that among these 20 tuples, only B2 , B3 , B4 , B5 , B8 ,
B9 , B11 , B12 , B13 , B15 , are a basis of the real vector space R2 . Again it
is important to remember that the two tuples B2 and B3 are not the
same since, although their constituent elements are the same (namely the
two vectors (1, 0) and (0, 1)), their order in B2 and B3 is not the same;
similarly, the two tuples B4 and B9 are not the same, the two tuples B5
and B15 are not the same, ...
c3 , +, ·), and consider the vectors
2. Consider the real vector space (R
1 0 0
v1 = 0 , v2 = 1 , v3 = 0 ,
0 0 1
1 1 0
v4 = 1 , v5 = 1 , v6 = 1 ,
0 1 1
c3 . Here again, it is easy to verify that each of the 3−tuples (v1 , v2 , v3 ),(v2 , v1 , v3 ),(v3 , v1 , v2 ),(v
of R
(and still many more!) are distinct bases of the real vector space R c3 ; on
the other hand, the tuples (v1 ), (v1 , v2 ), (v3 ), (v1 , v1 , v1 ), (v5 , v6 ) (to
give only a few examples) do not form a basis of R c3 .
The attentive reader may have noticed that in the case of the real vector space
R2 given above, all the tuples which did form a basis of R2 happened to have
exactly the same number of vectors (namely 2); similarly, in the case of the real
vector space Rc3 , all the tuples which did form a basis of R
c3 happened to have
exactly the same number of vectors (namely 3). We will see in the next lecture
that this is not accidental, and that in a vector space, all bases have the same
number of elements.
Now that we have defined the notion of basis for a real vector space, we define
another important notion, that of components of a vector with respect to
a basis:
101
α1 v1 + α2 v2 + · · · + αp vp = v.
The real number α1 is called the first component (or first coordinate) of v with
respect to the basis B, the real number α2 is called the second component (or
second coordinate) of v with respect to the basis B, ... and so on.
It is absolutely essential to note the following point:
• If B = (v1 , v2 , · · · , vp ) is a basis of V and v is an element of V, then v has
at least one component vector with respect to B (since {v1 , v2 , · · · , vp }
is a generating set for V); furthermore, since the vectors v1 , v2 , · · · , vp
are linearly independent, there is only a unique p−tuple (α1 , · · · , αp )
of real numbers for which α1 v1 + α2 v2 + · · · + αp vp is equal to v; i.e. to
each v in V there corresponds only a unique component vector. This is
why in the previous definition, we wrote “the component vector of v with
respect to B” and not “a component vector of v with respect to B”.
It is also important to make note of the following point:
• If B = (v1 , v2 , · · · , vp ) is a basis of V and v is an element of V, we have
defined the component vector of v with respect to B to be some p−tuple
of real numbers, i.e. some element of Rp ; we will find it occasionally
convenient to denote the component vector of v by an element of R cp
instead, i.e. by a column vector with p real entries.
We close this section with three examples:
(a) Consider the real vector space (R2 +, ·), and let W denote the subset of
R2 consisting of all pairs (x, y) of real numbers such that x = y. It is easy
to verify that W is a vector subspace of R2 , and hence, by Problem 4 of
Lecture 3, W (with the operations of addition and scalar multiplication
it receives from R2 ) is itself a real vector space. Consider the vector
v1 = (1, 1) of R2 ; clearly, v1 is an element of W, since W is the subset
of R2 consisting of all pairs (x, y) with x = y, and the pair (1, 1) satisfies
this condition. Consider the “1−tuple” B = (v1 ); let us show that B is a
basis of W. To do this, we have to show two things:
(i) That {v1 } is a generating set for W, and
(ii) that {v1 } is a linearly independent subset of W.
To show (i), namely that {v1 } is a generating set for W, we have to show
that any v ∈ W can be written as some linear combination of v1 ; let
then v = (x, y) be an element of W; since (x, y) is in W, we have (by
definition of W) x = y. Hence, v = (x, x) = x(1, 1) = xv1 , which shows
102 SECTION 8
which shows that the component vector of the vector (a, a) of W with
respect to the basis B is the 1−tuple (a). (Note: It may sound a bit
akward to talk about 1−tuples, but think of them as a list with only one
entry!).
(b) Consider the real vector space (R3 , +, ·), and let v1 = (1, 0, 0), v2 =
(0, 1, 0), v3 = (0, 0, 1), v4 = (1, 1, 0), v5 = (1, 1, 1). It is easy to verify
that each of the 3−tuples of vectors B1 = (v1 , v2 , v3 ), B2 = (v2 , v3 , v1 ),
and B3 = (v1 , v4 , v5 ) is a basis of R3 . Let now v = (a, b, c) be an element
of R3 . It is easy to verify that the component vector of v with respect to
B1 is (a, b, c) (i.e. v itself!), whereas its component vector with respect to
the basis B2 is (b, c, a), and its component vector with respect to the basis
B3 is (a − b, b − c, c).
(c) Consider again the real vector space (R3 , +, ·), and let W be the subset
of R3 consisting of all triples (x, y, z) with x + y + z = 0. It is easy to
verify that W is a vector subspace of R3 . By Problem 4 of Lecture 3, we
know that W (with the operations of addition and scalar multiplication
it receives from R3 ) is itself a real vector space. Consider the elements
v1 , v2 of W defined by:
v1 = (1, 0, −1),
v2 = (0, 1, −1).
PROBLEMS:
1. Consider the real vector space (R3 , +, ·), and consider the following vectors
in R3 :
v1 = (1, 0, 0), v2 = (1, 1, 0), v3 = (1, 2, 1), v4 = (0, 0, 3),
104 SECTION 8
2. Consider the real vector space c3 , +, ·), and let W be the subset of R
(R c3
x
consisting of all elements y of Rc3 for which x + y − z = 0. It is easy
z
to verify that W is a vector subspace of R c3 , and hence, is itself a real
vector space.
Consider now the following vectors in W:
1 0 1 0
v1 = 0 , v2 = 2 , v3 = −1 , v4 = 0 .
1 2 0 0
3. Consider the real vector space (R4 , +, ·), and let W be the subset of R4
consisting of all 4−tuples (x, y, z, w) of real numbers for which 3x+y−w =
0 and z − 2w = 0. It is easy to verify that W is a vector subspace of W,
and hence, is itself a real vector space. Consider the following vectors in
W:
v1 = (1, 0, 0, 1), v2 = (0, 1, 1, 2), v3 = (1, 1, 1, 3), v4 = (1, −1, −1, −1).
Study Topics
• Finite-dimensional vector spaces
• Dimension of a real vector space
107
108 SECTION 9
We know from elementary geometry that we can represent real numbers (i.e.
elements of the real vector space R) by points on a line, pairs of real numbers
(i.e. elements of the real vector space R2 ) by points in a plane, triples of real
numbers (i.e. elements of the real vector space R3 ) by points in space, ... The
question we may ask at this point is how the difference between, say, a line and
a plane, manifests itself in a difference between the vector spaces R and R2 .
Or how the difference between a plane and space manifests itself in a difference
between the vector spaces R2 and R3 . In other words, how is the “line-like”
feature of the real vector space R captured in the vector space structure of R,
how is the “plane-like” feature of the real vector space R2 captured in the vector
space structure of R2 , how is the “space-like” feature of the real vector space
R3 captured in the vector space structure of R3 , ... In one word: How do these
“geometric” concepts get expressed in terms of linear algebra ? We shall soon
see that one point of contact between linear algebra and geometry is given by
the notion of dimension.
Definition 23. Let (V, +, ·) be a real vector space.
• V is said to be finite-dimensional if there exists an integer N ≥ 0 such
that any subset of V containing N + 1 elements is linearly dependent.
The smallest integer N for which this holds is then called the dimension
of V (equivalently, V is said to have dimension N ).
• V is said to be infinite-dimensional if it is not finite-dimensional.
It is important to note the following points:
• By the above definition, if V has dimension N , then any subset of V
containing N + 1 elements is linearly dependent, and as a result, any
subset of V containing N + 2 or more elements is also linearly dependent.
• By the above definition, if V has dimension N , then any subset of V
containing N + 1 elements is linearly dependent, but there exists at least
one linearly independent subset of V containing exactly N elements,
since otherwise the dimension of V would be strictly less than N .
• V is infinite-dimensional if for any positive integer N , no matter how
large, there exists a linearly independent subset of V containing ex-
actly N elements.
Let’s recapitulate:
• V is said to have dimension N if and only if the following two conditions
are met:
(i) There exists a linearly independent subset of V containing ex-
actly N elements,
(ii) Any subset of V containing N +1 elements is linearly dependent.
Let us now examine some simple examples:
109
1. Recall the real vector space (F (R; R), +, ·) consisting of all functions f :
R → R (i.e. all real-valued functions of a real variable). Let us show that
F (R; R) is an infinite-dimensional vector space. To do this, we have
to show that for any integer N ≥ 0 (no matter how large), there exists a
linearly independent subset of F (R; R) containing exactly N elements.
Let then N be any integer ≥ 0. Consider the following elements f1 , f2 , · · · , fN
in F (R; R) defined as follows:
(You may want to draw the graphs of these functions just to see what
they look like). In other words, if k is an integer in the range 1 ≤ k ≤ N ,
fk is the function from R to R defined to be 0 everywhere except at
k (where it takes the value 1). Let us show that {f1 , f2 , · · · , fN } is a
linearly independent subset of F (R; R). Let then α1 , α2 , · · · , αN be
real numbers such that
α1 f1 + α2 f2 + · · · + αN fn = 0.
(Recall that the zero vector 0 of the vector space F (R; R) is the function
from R to R which maps every real number to zero). We have to show
that α1 ,α2 ,· · · ,αN must all be zero. Now, the equality
α1 f1 + α2 f2 + · · · + αN fn = 0.
is equivalent to
and in particular, this last equality must hold for particular choices of t,
such as t = 1, t = 2, · · · , t = N . But since f1 (1) = 1 and f2 (1) = f3 (1) =
· · · = fN (1) = 0, the above equality, for the particular choice of t = 1
becomes:
α1 = 0.
α2 = 0.
110 SECTION 9
α1 f1 + α2 f2 + · · · + αN fn = 0,
ξ+ξ = ξ
α·ξ = ξ, ∀α ∈ R.
It is easy to verify that (V, +, ·) is a real vector space and that ξ itself is
the zero vector of this vector space. Let us show that V has dimension
0. Let N = 0; there is only one possible subset of V containing N +
1 = 1 elements, namely {ξ}, which is V itself, and since that subset
contains the zero vector (since ξ itself is the zero vector), that subset is
linearly dependent. It follows from the definition of dimension that V is
0−dimensional.
3. Let now V = R, with the usual addition and multiplication operations (we
have already seen that with these operations, R is a real vector space). Let
us show that V has dimension 1. For this, we have to show that we can
find a linearly independent subset of R containing exactly one element,
but that any subset of R containing two elements is linearly dependent.
Consider then the subset of R given by {1} (i.e. the subset containing
only the vector 1 of R). Let α ∈ R, and assume α · 1 = 0; since α · 1 = α,
this implies α = 0, which shows that {1} is a linearly independent subset
of R. Let us now show that any subset of R containing two elements is
linearly dependent. Let then {v1 , v2 } be any subset of R; we have to show
that there exists a linear combination αv1 + βv2 of v1 , v2 that is equal to
0 but such that α, β are not both zero. Consider the following three cases:
(i) If v1 = 0, then with α = 1, β = 0, we obtain αv1 +βv2 = 0, and since
α, β are not both zero (since α = 1), this proves linear dependence
of {v1 , v2 } in this case.
(ii) If v2 = 0, then with α = 0, β = 1, we obtain αv1 +βv2 = 0, and since
α, β are not both zero (since β = 1), this proves linear dependence
of {v1 , v2 }.
111
The reader may feel at this point that the dimension of a real vector space may
not be such an easy thing to compute: After all, it took us quite a few lines
above to compute the dimension of a simple vector space such as R. Recall
again that to prove that a real vector space has dimension N , we have to do
two things:
The first item is usually easy to accomplish; it is enough to find one linearly
independent subset containing exactly N elements. The difficulty comes from
the second item: How does one prove that any subset containing N +1 elements
is necessarily linearly dependent ?
The following theorems will yield an extremely simple way to compute the
dimension of a finite-dimensional real vector space.
Proof. The proof follows directly from the definition of dimension. If we had p >
N , then that would mean that p ≥ N + 1; since V is assumed to have dimension
N , we know (by definition of dimension) that any subset of V containing N + 1
or more elements has to be linearly dependent; hence, if p > N , i.e. equivalently
p ≥ N + 1, the subset {v1 , · · · , vp } has N + 1 or more elements, and hence
is linearly dependent. Hence, if we assume that {v1 , · · · , vp } is linearly
independent, then we must have p ≤ N .
Proof. We will show that, under the assumption that {v1 , · · · , vp } is a gen-
erating set for V, we cannot have p < N . So let us begin by assuming that
{v1 , · · · , vp } is a generating set for V, and that p < N ; we will show that this
leads to a contradiction.
Since V is assumed to have dimension N , there does exist (by definition of
dimension) a linearly independent subset of V containing exactly N ele-
ments; let {e1 , · · · , eN } denote that subset. Since {v1 , · · · , vp } is a generating
set for V, each of e1 , e2 , · · · , eN can be expressed as a linear combination of
v1 , · · · , vp . Let us start with e1 ; we can write:
e1 = α1 v1 + α2 v2 + · · · + αp vp ,
e2 = λ1 e1 + λ2 v2 + λ3 v3 + · · · + λp vp ,
We can now combine the previous two theorems in the following theorem:
We can establish the following two very important corollaries of the previous
theorem:
(i) In a finite-dimensional real vector space, all bases have the same num-
ber of elements (and that number is the dimension of the vector space);
1. Consider the (by now extremely familiar) real vector space (R2 , +, ·). Let
us compute the dimension of R2 . As we have seen above, all we need to
do is find a basis for R2 ; the dimension of R2 will then be equal to the
number of elements of that basis. Consider then the following vectors in
R2 :
2. Consider the (equally familiar) real vector space (R3 , +, ·). Let us compute
the dimension of R3 . As we have seen above, all we need to do is find a
basis for R3 ; the dimension of R3 will then be equal to the number of
elements of that basis. Consider the following vectors in R3 :
3. More generally, let n be an integer ≥ 1 and consider the real vector space
(Rn , +, ·). Let us compute the dimension of Rn . Proceeding as before, let
us construct a basis for Rn ; the dimension of Rn will then be given by the
number of elements in that basis. For this, consider the following elements
of Rn :
v1 = (1, 0, 0, 0, · · · , 0),
v2 = (0, 1, 0, 0, · · · , 0),
v3 = (0, 0, 1, 0, · · · , 0),
···
vn = (0, 0, 0, 0, · · · , 1).
(i.e. vk is the n−tuple with the k th entry equal to 1 and all other entries
equal to 0). It is easy to verify that (v1 , v2 , v3 , · · · , vn ) is a basis of Rn .
Since it has exactly n elements, this shows that Rn has dimension n.
5. Consider now the vector subspace W of R3 defined as the set of all (x, y, z)
in R3 with x + y + z = 0. Consider the vectors v1 = (1, 0, −1) and
v2 = (0, 1, −1) of W. We have already seen (in the previous lecture) that
(v1 , v2 ) is a basis of W. Since it has 2 elements, we conclude that W is
a real vector space of dimension 2.
PROBLEMS:
1. Let C(R; R) be the real vector space of all continuous functions from R to
R (with addition and scalar multiplication defined as in F (R; R)). Show
that C(R; R) is an infinite-dimensional vector space.
10. Consider the vector subspace W of R3 defined as the set of all (x, y, z) ∈ R3
such that y = 0. Show that W has dimension 2.
11. Consider the vector subspace W of R3 defined as the set of all (x, y, z) ∈ R3
such that z = 0. Show that W has dimension 2.
12. Consider the vector subspace W of R3 defined as the set of all (x, y, z) ∈ R3
such that x + y = 0 and z = 0. Show that W has dimension 1.
13. Consider the vector subspace W of R3 defined as the set of all (x, y, z) ∈ R3
such that x − z = 0 and x + y + z = 0. Show that W has dimension 1.
14. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x + y + z + w = 0. Show that W has dimension 3.
15. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x + w = 0. Show that W has dimension 3.
16. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x − y = 0. Show that W has dimension 3.
17. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x = 0. Show that W has dimension 3.
18. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that w = 0. Show that W has dimension 3.
19. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x + y = 0 and z + w = 0. Show that W has dimension 2.
20. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x+y +z +w = 0 and z −w = 0. Show that W has dimension
2.
116 SECTION 9
21. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x + 2w = 0 and 2z + w = 0. Show that W has dimension 2.
22. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x − y = 0 and z + w = 0 and y + w = 0. Show that W has
dimension 1.
23. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x + 2y + z = 0 and z = 0 and 2y + w = 0. Show that W
has dimension 1.
24. Consider the vector subspace W of R4 defined as the set of all (x, y, z, w) ∈
R4 such that x + 2y + z + w = 0 and x + y = 0 and y − z = 0. Show that
W has dimension 1.
Section 10
Study Topics
• Linear Transformations
• Range and Kernel of a Linear Transformation
117
118 SECTION 10
In the previous lectures, we have studied vector spaces by themselves; for ex-
ample, we have started from a real vector space, and we have found subspaces
of that vector space, generating sets for that vector space, linearly dependent or
independent subsets for that vector space, bases for that vector space, ... and so
on. We now consider functions between vector spaces, i.e. functions from
one vector space to another. Among these functions, there are some which have
a desirable property (called linearity) which we will precisely define shortly;
these functions are called linear, and are the main object of study of this and
the next few lectures.
A WORD ON TERMINOLOGY: We will often use the words mapping
and transformation instead of function – They will all mean the same thing.
Definition 24. Let V and W be two real vector spaces, and let L : V → W
be a mapping from V to W. L is said to be a linear mapping (also called a
linear transformation or a linear function) if the following two properties are
verified:
1. For any v1 , v2 ∈ V, L(v1 + v2 ) = L(v1 ) + L(v2 );
2. For any α ∈ R and any v ∈ V, L(αv) = αL(v).
Before going further, let us give some examples of linear mappings:
1. Let V be a real vector space, and let L : V → V be the identity mapping
of V, defined by L(v) = v for all v ∈ V; it is easy to verify that L is a
linear mapping.
2. Let V be a real vector space, let α ∈ R be any real number, and let
L : V → V be the mapping defined by L(v) = αv for all v ∈ V; it is easy
to verify that L is a linear mapping.
3. Let V and W be real vector spaces, and let 0W denote the zero vector of
W. Let L : V → W denote the zero mapping from V to W, defined
by L(v) = 0W for all v ∈ V (i.e. everything in V is mapped to the zero
vector of W); here again, it is easy to verify that L is a linear mapping.
4. Consider the familiar vector spaces R2 and R, and let L : R2 → R be the
mapping defined by L((x, y)) = 2x + 3y for all (x, y) ∈ R2 . Let us show
that L is a linear mapping:
(i) Let v1 = (x1 , y1 ) ∈ R2 and v2 = (x2 , y2 ) ∈ R2 . We have to show
that L(v1 + v2 ) = L(v1 ) + L(v2 ). We have:
L(v1 + v2 ) = L((x1 , y1 ) + (x2 , y2 ))
= L((x1 + x2 , y1 + y2 ))
= 2(x1 + x2 ) + 3(y1 + y2 )
= (2x1 + 3y1 ) + (2x2 + 3y2 )
= L((x1 , y1 )) + L((x2 , y2 ))
= L(v1 ) + L(v2 ),
119
6. Consider again the familiar real vector spaces R2 and R3 , and let L : R2 →
R3 be the mapping defined by L((x, y)) = (x, y, 0) for all (x, y) ∈ R2 . It
is easy to verify that L is a linear mapping.
7. Consider once more the familiar real vector spaces R2 and R3 , and let
L : R3 → R2 be the mapping defined by L((x, y, z)) = (x, y) for all
(x, y, z) ∈ R3 . It is easy to verify that L is a linear mapping.
8. Consider now the familiar real vector spaces R3 and R, and let L : R3 → R
be the mapping defined by L((x, y, z)) = x + y + z for all (x, y, z) ∈ R3 .
It is easy to verify that L is a linear mapping.
9. Consider now the real vector spaces F (R; R) and R (recall F (R; R) is the
real vector space of all functions from R to R), and let L : F (R; R) → R
be the mapping defined by L(f ) = f (0) for all f ∈ F (R; R). Let us prove
that L is a linear mapping:
Before going further and giving examples of mappings between vector spaces
which are not linear, we prove the following useful and simple theorem:
Theorem 15. Let V and W be real vector spaces; let 0V denote the zero vector
of V, and let 0W denote the zero vector of W. let L : V → W be a linear
mapping. Then, we have:
L(0V ) = 0W .
Proof. Let α = 0 (i.e. α is the real number zero); we have α0V = 0V (we have
already shown that any real number times the zero vector is equal to the zero
vector, and in particular if that real number happens to be zero), and therefore,
we can write:
L(0V ) = L(α0V )
= αL(0V )
= 0L(0V )
= 0W
since, as we have already shown, the real number 0 times any vector of a vector
space is equal to the zero vector of that vector space.
Let us now examine mappings between vector spaces, which are not linear.
L(0V ) = w 6= 0W ,
i.e. L maps the zero vector of V to something other than the zero vector
of W; it follows therefore from the previous theorem that the mapping L
is not linear.
L(0V ) = λ0V + v1
= 0V + v1
= v1 6= 0V ,
and since L maps the zero vector of V to something other than the zero
vector of V, it follows from the previous theorem that the mapping L is
not linear.
3. Consider the real vector spaces R2 and R, and let L : R2 → R be the
mapping defined by L((x, y)) = x2 + y 2 for all (x, y) ∈ R2 . Let us show
that the mapping L is not linear. The previous theorem is not of much
help here, since L does indeed map the zero vector of R2 (namely the pair
(0, 0)) to the zero vector of R (namely the real number 0); indeed, we
have:
L((0, 0)) = 02 + 02 = 0 + 0 = 0;
so we cannot use the previous theorem to show that L is not linear, unlike
what we did in the previous examples.
So how do we proceed ? Well, we have to go back to the definition. From
there, it is clear that it is enough to find two vectors v1 , v2 in R2 for which
L(v1 + v2 ) is not equal to L(v1 ) + L(v2 ), or a real number α and a vector
v in R2 such that L(αv) is not equal to αL(v), in order to conclude that
L is not a linear mapping.
Let us choose for example α = 2 and v = (1, 1). For these choices of α
and v, we have:
Hence, we have found some real number α (namely α = 2) and some vector
v in R2 (namely v = (1, 1)) for which L(αv) is not equal to αL(v). This
shows that the second property in the definition of a linear mapping does
not hold for our mapping L, and we conclude from this that the mapping
L is not linear.
4. Consider the real vector space R2 , and let L : R2 → R2 be the mapping
defined by L((x, y)) = (x3 + y, x − y) for all (x, y) ∈ R2 . It is easy to verify
here as well that the mapping L defined here is not linear.
5. Consider the real vector space C([0, 1]; R) of all continuous functions from
[0, 1] to R, and define the mapping L : C([0, 1]; R) → R from the real
vector space C([0, 1]; R) to the real vector space R as follows: L(f ) =
R1
0
(f (t))2 dt for all f ∈ C([0, 1]; R). Let us show that L is not a linear
123
mapping. For this, let us show that the second property in the definition
of a linear mapping fails in this case. Let then α = 2, and define the vector
f ∈ C([0, 1]; R) as follows (recall that in the real vector space C([0, 1]; R),
a vector is nothing other than a continuous function from [0, 1] to R):
f (t) = 1, for all t ∈ R (in other words, f is defined to be the constant
function on [0, 1] which maps every t ∈ [0, 1] to the real number 1). Clearly,
f is continuous on [0, 1] with values in R, and hence is an honest element
of C([0, 1]; R). With these choices of α ∈ R and f ∈ C([0, 1]; R), let us
compute L(αf ) and αL(f ). We have:
Z 1
L(αf ) = ((αf )(t))2 dt
0
Z 1
= (αf (t))2 dt
0
Z 1
= α2 (f (t))2 dt
0
Z 1
2
= α (f (t))2 dt
0
Z 1
= 22 (1)2 dt
0
Z 1
= 4 dt
0
= 4,
whereas
Z 1
αL(f ) = α (f (t))2 dt
0
Z 1
= 2 (1)2 dt
0
Z 1
= 2 dt
0
= 2,
which shows that for these choices of α ∈ R and f ∈ C([0, 1]; R), L(αf )
is not equal to αL(f ), and this shows that the second property in the
definition of a linear mapping does not hold for our mapping L, and we
conclude from this that the mapping L is not linear.
Now that we have seen examples of linear (and non-linear) mappings between
vector spaces, let us examine an interesting property of linear mappings: Con-
sider three real vector spaces V, W, and Z, let L1 : V → W be a mapping
from the real vector space V to the real vector space W, and L2 : W → Z a
a mapping from the real vector space W to the real vector space Z. We know
124 SECTION 10
Theorem 16. Let V, W, and Z be real vector spaces, and let L1 : V → W and
L2 : W → Z be linear mappings. Then, the mapping L2 ◦ L1 : V → Z is also
linear.
Let us now consider two real vector spaces V and W and a linear mapping
L : V → W. Consider the set of all v in V which are mapped to the zero vector
0W of W by L, i.e. the set of all v ∈ V for which L(v) = 0W . Let us give a
name to the set of all such elements of V; we denote by ker(L) (and we call this
“Kernel of L or Null Space of L”) the set of all v ∈ V for which L(v) = 0W ;
we can write this statement more formally as:
ker(L) = {v ∈ V|L(v) = 0W }.
What do we know about ker(L) ? Well, the first thing we know about ker(L),
which follows straight from its definition, is that it is a subset of V. But as
the following theorem will show, ker(L) is not just any old subset of V; rather,
it is a vector subspace of V.
Theorem 17. Let V and W be real vector spaces, and let L : V → W be a linear
mapping. Then, the kernel ker(L) of L is a vector subspace of V.
Proof. To show that ker(L) is a vector subspace of the real vector space V, we
have to show that:
(i) The zero vector 0V of V is in ker(L),
(ii) for any v1 and v2 in ker(L), v1 + v2 is also in ker(L),
(iii) for any v in ker(L) and any α in R, αv is also in ker(L).
• We know from the first theorem we proved in this lecture that (by linearity
of L) L(0V ) = 0W , from which we obtain that 0V is indeed in ker(L)
(by definition of ker(L)); i.e. we have shown property (i), namely that
0V ∈ ker(L).
• Let us now prove property (ii). Let then then v1 , v2 ∈ ker(L); we have
to show that v1 + v2 is also in ker(L), i.e. that L(v1 + v2 ) = 0W . But,
by linearity of L and the fact that v1 and v2 are assumed to be in ker(L)
(and therefore L(v1 ) = 0W and L(v2 ) = 0W ):
L(v1 + v2 ) = L(v1 ) + L(v2 )
= 0W + 0W
= 0W ,
which proves that v1 + v2 is also in ker(L). This proves property (ii).
• Let us now prove property (iii). Let then then v ∈ ker(L) and α ∈ R; we
have to show that αv is also in ker(L), i.e. that L(αv) = 0W . But, by
linearity of L and the fact that v is assumed to be in ker(L) (and therefore
L(v) = 0W ):
L(αv) = αL(v)
= α0W
= 0W ,
which proves that αv is also in ker(L). This proves property (iii).
126 SECTION 10
We have shown that ker(L) satisfies all three properties of a vector subspace;
this shows that ker(L) is a vector subspace of V.
Recall that a mapping L : V → W is said to be injective or one-to-one if for
any two elements v1 , v2 ∈ V with v1 6= v2 it must follow that L(v1 ) 6= L(v2 );
equivalently, L is injective if for any v1 , v2 ∈ V, the equality L(v1 ) = L(v2 )
implies v1 = v2 (i.e two distinct elements in V cannot map to the same value
under L). The following theorem relates the property of L being injective to
the kernel ker(L) of L.
Theorem 18. Let L : V → W be a linear mapping from the real vector space V
to the real vector space W. We have:
• L is injective if and only if the kernel ker(L) of L is equal to {0V }, i.e.
ker(L) = {0V }.
Proof. (i) Assume first that L is injective; we have to show then that if v ∈
ker(L), then v = 0V (i.e. the only element of ker(L) is the zero vector of
V). Let then v ∈ ker(L). Then, L(v) = 0W (by definition of v being an
element of the kernel ker(L) of L). On the other hand, we also know that
L(0V ) = 0W . We therefore have:
L(v) = L(0V ).
Since L is assumed injective, it must follow that v = 0V . Hence, we have
shown that if v is any element in ker(L), then it must follow that v = 0V ;
this shows that ker(L) = {0V }.
(ii) Let us now assume that ker(L) = {0V }; we have to show then that for
any v1 , v2 ∈ V, the equality L(v1 ) = L(v2 ) implies that v1 = v2 . Let
then v1 , v2 ∈ V and assume that L(v1 ) = L(v2 ). Hence, by linearity of
L:
L(v1 − v2 ) = L(v1 + (−1)v2 )
= L(v1 ) + L((−1)v2 )
= L(v1 ) + (−1)L(v2 )
= L(v1 ) − L(v2 )
= 0W ,
which shows that v1 − v2 ∈ ker(L) (since L(v1 − v2 ) = 0W ). Since we
have assumed that ker(L) = {0V }, it must follow that v1 − v2 = 0V ,
i.e. that v1 = v2 . This is what we wanted to show, and this establishes
injectivity of L.
Theorem 19. Let V and W be real vector spaces, and let L : V → W be a linear
mapping. Then, the range Im(L) of L is a vector subspace of W.
Proof. We have to prove that Im(L) satisfies the three properties that a vector
subspace of W should satisfy, namely, we have to show that:
• We know from the first theorem we proved in this lecture that (by linearity
of L) L(0V ) = 0W , from which we obtain that there exists an element in
V (namely the zero vector of V) which is mapped by L to the zero vector
0W of W; this shows that 0W is in the range of L, i.e. 0W ∈ Im(L).
Hence, we have shown property (i).
• Let us now prove property (ii). Let then then w1 , w2 ∈ Im(L); we have
to show that w1 + w2 is also in Im(L), i.e. we have to show that w1 + w2
is L of something in V. Since we have w1 ∈ Im(L), by definition of
Im(L), there must exist a vector in V, call it v1 , such that L(v1 ) = w1 .
Similarly, since we have w2 ∈ Im(L), by definition of Im(L), there must
exist a vector in V, call it v2 , such that L(v2 ) = w2 . Consider now the
vector v in V defined by v = v1 + v2 (i.e. v is defined to be the sum of
the two vectors v1 and v2 ). We have:
L(v) = L(v1 + v2 )
= L(v1 ) + L(v2 )
= w1 + w2 ,
• Let us now prove property (iii). Let then then w ∈ Im(L) and α ∈ R; we
have to show that αw is also in Im(L), i.e. we have to show that αw is L
of something in V. We have:
L(αv) = αL(v)
= αw,
We have shown that Im(L) satisfies all three properties of a vector subspace;
this shows that Im(L) is a vector subspace of W.
We close this section with some important terminology and an important the-
orem which we shall merely state and not prove:
rank(L) + nullity(L) = N.
PROBLEMS:
Show which of the following mappings between real vector spaces are linear and
which are not linear:
1. L : R → R, defined by L(x) = 2x + 1.
2. L : R → R, defined by L(x) = x2 + x.
Study Topics
• Real Matrices
• Linear Transformations Defined by Matrices
• Range and Kernel of a Matrix
131
132 SECTION 11
In the previous lecture, we have studied a special class of functions from one
real vector space to another, namely those that were linear, and we called them
linear functions, or, equivalently, linear mappings, or equivalently, linear
transformations. In this lecture, we introduce the notion of a real “matrix”.
What is a real matrix ? Think of it as a table, with, say m rows and n columns,
of real numbers. That’s it! So what’s the relation between a matrix and a linear
transformation ? As we will see in this lecture, a real matrix with m rows and
n columns will allow us to define a linear transformation from R cn to Rdm.
Definition 26. Let m and n be integers ≥ 1. A real matrix with m rows and
n columns (also called a real m × n matrix) is a table (or array) of the form:
a1,1 a1,2 · · · a1,n
a2,1 a2,2 · · · a2,n
.. .. .. ,
. . ··· .
am,1 am,2 ··· am,n
where all the entries (i.e. a1,1 , a1,2 , a1,3 , · · · ) are real numbers.
ai,j is called the entry of the matrix on row i and column j.
Let us look immediately at some examples:
√
2
1. The matrix is a real 2 × 1 real matrix (2 rows, 1 column).
−π
2 0√
2. The matrix is a 2 × 2 real matrix (2 rows, 2 columns).
−1 1/ 3
1 0
3. The matrix −1 1 is a 3 × 2 real matrix (3 rows, 2 columns).
2 3
1 0 −2 1
4. The matrix is a 2×4 real matrix (2 rows, 4 columns).
−1 1 2 3
5. The matrix 1 3 2 1 5 is a 1 × 5 real matrix (1 row, 5 columns).
6. The matrix (−7) is a 1 × 1 real matrix (1 row, 1 column).
It is important to point out that the place of the elements in a matrix matters;
exchanging the place of two elements in the matrix will change the matrix
if
1 0
those two elements are not equal. For example, the matrices and
0 0
0 0
are not equal. More generally, if A is defined to be the m × n real
0 1
matrix
a1,1 a1,2 · · · a1,n
a2,1 a2,2 · · · a2,n
A= . .. .. ,
.. . ··· .
am,1 am,2 · · · am,n
133
then we have A = B if and only if ai,j = bi,j for all 1 ≤ i ≤ m and all 1 ≤ j ≤ n,
i.e. the corresponding entries of A and B must be equal.
Definition 27. Let
a1,1 a1,2 ··· a1,n
a2,1 a2,2 ··· a2,n
A= .. .. .. ,
. . ··· .
am,1 am,2 ··· am,n
a1,j
a2,j
A;j = a3,j .
..
.
am,j
Similarly, letting i ∈ {1, 2, · · · , m}, the ith row vector of A, denoted by Ai; is
defined to be the element of Rn given by:
We shall very often write the row vector Ai; without the separating commas,
i.e., as:
Ai; = ai,1 ai,2 ai,3 · · · ai,n ;
despite the lack of the separating commas, Ai; will still be understood as an n−tuple
of real numbers, i.e. as an element of Rn .
3 2
Let us look at an example: Consider the 3 × 2 real matrix −1 0 ; its first
0 1
3
column vector is the element of R c3 given by −1 , and its second column
0
2
vector is the element of R c3 given by 0 . Its first row vector is the element
1
134 SECTION 11
of R2 given by (3, 2), its second row vector is the element of R2 given by (−1, 0),
and finally, its third row vector is the element of R2 given by (0, 1). As was
mentionned above, we shall often omit the separating commas when writing
row vectors of matrices;
for example, we shall write the first row vector of this
matrix as 3 2 instead, and its second and third row vectors as −1 0
and 0 1 , respectively, and we shall still consider them as elements of R2
(i.e. pairs of real numbers), despite the absence of the separating commas.
Definition 28. Let m, n be integers ≥ 1. We denote the set of all m × n real
matrices by Mm,n (R).
In other words, an element of Mm,n (R) is nothing other than a real m × n
matrix, and any real m × n matrix is an element of Mm,n (R).
We now define two operations on Mm,n (R), namely addition of real matri-
ces and multiplication of a real matrix by a real number. Let us define
the addition operation first. Let then
a1,1 a1,2 · · · a1,n
a2,1 a2,2 · · · a2,n
A= . .. .. ,
.. . ··· .
am,1 am,2 ··· am,n
and
b1,1 b1,2 ··· b1,n
b2,1 b2,2 ··· b2,n
B= .. .. .. ,
. . ··· .
bm,1 bm,2 ··· bm,n
be two real m × n matrices. We define the sum of A and B, and we denote by
A + B, the real m × n matrix given by:
a1,1 + b1,1 a1,2 + b1,2 · · · a1,n + b1,n
a2,1 + b2,1 a2,2 + b2,2 · · · a2,n + b2,n
A+B = .. .. .. ;
. . ··· .
am,1 + bm,1 am,2 + bm,2 ··· am,n + bm,n
in other words, A + B is the real m × n matrix obtained by adding together the
corresponding entries of the matrices A and B.
Let now α be a real number; we define the product of the matrix A and the real
number α, and we denote by α · A (the “·”, denoting the scalar multiplication
operation – as with vectors, we shall drop the “·” in the notation very soon ...),
the real m × n matrix given by:
αa1,1 αa1,2 · · · αa1,n
αa2,1 αa2,2 · · · αa2,n
α·A= .. .. .. ;
. . ··· .
αam,1 αam,2 · · · αam,n
135
Theorem 21. (Mm,n (R), +, ·) is a real vector space. The zero vector 0 of
this vector space is the real m × n matrix with all entries equal to 0, i.e.,
0 0 ··· 0
0 0 ··· 0
0= . . . ;
.. .. · · · ..
0 0 ··· 0
The proof of this theorem is immediate and is left to the reader. The bottom
line is that we can think of real m × n matrices themselves as vectors (in the
real vector space Mm,n (R)) ...
Let now again m, n be integers ≥ 1 and let A ∈ Mm,n (R) be the matrix defined
by:
a1,1 a1,2 · · · a1,n
a2,1 a2,2 · · · a2,n
A= . .. .. ,
.. . ··· .
am,1 am,2 · · · am,n
cn
where the ai,j are some given real numbers. Consider the function LA from R
d
to R cn → R
m , i.e. the function L : R
A
d m defined as follows: Let
x1
x2
v = x3
..
.
xn
a1,1 x1 + a1,2 x2 + · · · + a1,n xn
a2,1 x1 + a2,2 x2 + · · · + a2,n xn
LA (v) = .. .
.
am,1 x1 + am,2 x2 + · · · + am,n xn
1. Consider the real 1 × 4 matrix A given by A = −1 1 0 3 ; The
mapping LA defined by A is the mapping LA : R c4 → R (we consider R
c1
to be just R itself!) given by the following rule: Let
x1
x2
v= x3
x4
Av = −x1 + x2 + 3x4 .
−1 2
2. Consider now the real 3 × 2 matrix A given by A = 3 0 ; The
4 5
c2 c3
mapping LA defined by A is the mapping LA : R → R given by the
following rule: Let
x1
v=
x2
and again, we shall often simply write Av instead of LA (v), i.e., we shall
write
−x1 + 2x2
Av = 3x1 .
4x1 + 5x2
i.e. the matrix A has 1’s on the diagonal, and 0’s everywhere else. As we
cn → R
have seen, A defines therefore a function LA : R cn as follows: Let
x1
x2
x3
v=
..
.
xn
i.e., for any v ∈ Rcn , we have: LA (v) = v. This shows that LA is nothing
other than the identity mapping of R cn (i.e. the function which maps
every vector in Rcn to itself); for this reason, the matrix A in this example
is usually called the n×n identity matrix. Hence, with A being the n×n
identity matrix, we can write (in keeping with our simplified notation):
Av = v for every v ∈ R cn .
Let now A be a real m × n matrix (i.e. an element of Mm,n (R); what can we
say about the function LA : Rcn → R dm ? Does it have any special property ?
cn → R
Theorem 22. Let A ∈ Mm,n (R); the mapping LA : R d m defined by A is
linear.
Proof. Let
a1,1 a1,2 ··· a1,n
a2,1 a2,2 ··· a2,n
A= .. .. ..
. . ··· .
am,1 am,2 ··· am,n
be any two elements in Rcn ; we have to show that LA (v1 +v2 ) = LA (v1 )+LA (v2 ).
Let us first calculate LA (v1 + v2 ), and before we do this, let us first calculate
v1 + v2 . We have:
x1 + y1
x2 + y2
v1 + v2 = x3 + y3 ,
..
.
xn + yn
and therefore, by definition of LA ,
a1,1 (x1 + y1 ) + a1,2 (x2 + y2 ) + · · · + a1,n (xn + yn )
a2,1 (x1 + y1 ) + a2,2 (x2 + y2 ) + · · · + a2,n (xn + yn )
LA (v1 + v2 ) = .. .
.
am,1 (x1 + y1 ) + am,2 (x2 + y2 ) + · · · + am,n (xn + yn )
Separating the xi ’s from the yj ’s in the expression obtained, we can write:
(a1,1 x1 + · · · + a1,n xn ) + (a1,1 y1 + · · · + a1,n yn )
(a2,1 x1 + · · · + a2,n xn ) + (a2,1 y1 + · · · + a2,n yn )
LA (v1 + v2 ) = ..
.
(am,1 x1 + · · · + am,n xn ) + (am,1 y1 + · · · + am,n yn )
a1,1 x1 + · · · + a1,n xn a1,1 y1 + · · · + a1,n yn
a2,1 x1 + · · · + a2,n xn a2,1 y1 + · · · + a2,n yn
= .. + ..
. .
am,1 x1 + · · · + am,n xn am,1 y1 + · · · + am,n yn
= LA (v1 ) + LA (v2 ).
cn , we have LA (v1 +
This shows that property (i) holds, namely, that ∀v1 , v2 ∈ R
v2 ) = LA (v1 ) + LA (v2 ).
Let us now prove property (ii); for this, let α ∈ R be any real number, and let
x1
x2
v = x3
..
.
xn
139
i.e. for i ∈ {1, · · · , n}, ei has entry 1 on row i and entry 0 everywhere else.
cn , and it is called the canonical basis
Note that (e1 , e2 , · · · , en ) is a basisof R
x1
x2
cn
of R . Note also that any element x3 ∈ R
cn can be written as
..
.
xn
x1
x2
x3
= x1 · e1 + x2 · e2 + · · · + xn · en .
..
.
xn
d
L(e1 ), L(e2 ), · · · , L(en ) are all elements of R m , i.e. are column vectors with m
Let now C ∈ Mm,n (R) be the real m × n matrix constructed from the entries
of L(e1 ), L(e2 ), · · · , L(en ) as follows:
c1,1 c1,2 ··· c1,n
c2,1 c2,2 ··· c2,n
C= .. .. .. .
. . ··· .
cm,1 cm,2 ··· cm,n
Note that the first column vector of C is nothing other than L(e1 ), the second
column vector of C is nothing other than L(e2 ), and so on.
141
x1
x2
x3 cn
We have, for any element ∈R :
..
.
xn
x1
c1,1 x1 + c1,2 x2 + · · · + c1,n xn
x2
c2,1 x1 + c2,2 x2 + · · · + c2,n xn
LC ( x3 ) = ..
.. .
.
cm,1 x1 + cm,2 x2 + · · · + cm,n xn
xn
c1,1 c1,2 c1,n
c2,1 c2,2 c2,n
c3,1 c3,2 c3,n
= x1 · + x2 · + · · · + xn ·
.. .. ..
. . .
cm,1 cm,2 cm,n
= x1 · L(e1 ) + x2 · L(e2 ) + · · · + xn · L(en )
x1
x2
= L( x3 ).
..
.
xn
Hence, LC = L.
Now that we have shown that for any real m × n matrix A the mapping LA :
cn → R
R d m is linear (and the converse result), we can ask some further questions;
in particular, we can try (in the spirit of the previous lecture) to identify the
kernel and range of LA ...
Before we go any further, a few words on notation and terminology:
• As has been stated a number of times already, we shall often write Av
instead of LA (v); in that notation, the linearity of LA that we have just
established means that A(v1 + v2 ) = Av1 + Av2 and A(αv) = αAv.
• We shall denote the kernel of LA by ker(A) instead of the more cumber-
some (but correct!) notation ker(LA ), and we shall often just say “kernel
of the matrix A” instead of the more cumbersome “kernel of the linear
mapping LA defined by the matrix A”.
• Similarly, we shall denote the range of LA by Im(A) instead of the more
cumbersome (but correct!) notation Im(LA ), and we shall often just say
“range (or image) of the matrix A” instead of the more cumbersome
“range (or image) of the linear mapping LA defined by the matrix A”.
142 SECTION 11
Let then
a1,1 a1,2 ··· a1,n
a2,1 a2,2 ··· a2,n
A= .. .. .. ,
. . ··· .
am,1 am,2 ··· am,n
be a real m × n matrix (i.e. an element of Mm,n (R)), and let us first try to
identify the range of A, i.e. the range of the linear mapping LA : R cn →
d m
R . Recall that the range of LA , which we shall simply write as Im(A), is
d
the set of all LA (v) in R cn ; i.e., writing this formally, we have:
m with v in R
cn }.
Im(A) = {LA (v)|v ∈ R
Let now
x1
x2
x3
v=
..
.
xn
cn ; we can write:
be any element in R
a1,1 x1 + a1,2 x2 + · · · + a1,n xn
a2,1 x1 + a2,2 x2 + · · · + a2,n xn
LA (v) = ..
.
am,1 x1 + am,2 x2 + · · · + am,n xn
a1,1 a1,2 a1,n
a2,1 a2,2 a2,n
= x1 . + x2 . + · · · + xn .
.. .. ..
am,1 am,2 am,n
= x1 A;1 + x2 A;2 + · · · xn A;n ,
where, as we have already seen, A;1 denotes the first column vector of A, A;2
denotes the first column vector of A, ..., and so on. Recall that the column
vectors of A, namely A;1 , A;2 , · · · , A;n are all elements of the real vector space
d
R m (since A is m × n, i.e. has m rows!).
x1
x2
cn , then LA (
in R x3
) is nothing other than the linear combination
..
.
xn
of A;1 , A;2 , · · · , A;n , i.e. assume there exist real numbers α1 , α2 , · · · , αn such
that:
We therefore have:
a1,1 a1,2 a1,n
a2,1 a2,2 a2,n
w = α1 .. + α2 .. + · · · + αn ..
. . .
am,1 am,2 am,n
a1,1 α1 + a1,2 α2 + · · · + a1,n αn
a2,1 α1 + a2,2 α2 + · · · + a2,n αn
= ..
.
am,1 α1 + am,2 α2 + · · · + am,n αn
α1
α2
= LA ( α3 ),
..
.
αn
d
vectors of A is nothing other than the vector subspace of R m generated by the
column vectors of A, i.e. the linear span of the column vectors of A. We have
therefore proved the following theorem:
Theorem 24. Let A be a real m × n matrix; the range Im(A) of the matrix A
(i.e. the range of the linear mapping LA : R cn → R
d m ) is the linear span of the
The observations made in these two examples can be packaged into a theorem:
145
Theorem 25. Let A be a real m × n matrix (i.e. with m rows and n columns). We
have:
d
• If Im(A) = R m , then we must have n ≥ m. Hence, if n < m, then Im(A)
d
cannot be equal to R m.
Proof. We know that Im(A) is equal to the linear span of the n column vectors
d
A;1 , A;2 , · · · , A;n of A. Recall that each of A;1 , A;2 , · · · , A;n is a vector in R m . If
Im(A) = R d, then the subset {A;1 , A;2 , · · · , A;n } is a generating set for R
m d m;
that generating subset has n elements and R d m has dimension m (we have
seen this in Lecture 8), and it follows therefore from Theorem 10 of Lecture 8
that if Im(A) = R d m then we must have n ≥ m.
Now that we have investigated the range Im(A) of a real matrix A, let us examine
its kernel ker(A) (again, by ker(A) we mean the kernel of the linear map LA
defined by A, i.e. ker(LA )). Let then
a1,1 a1,2 ··· a1,n
a2,1 a2,2 ··· a2,n
A= .. .. .. ,
. . ··· .
am,1 am,2 ··· am,n
be a real m × n matrix, and let us try to get some information about the kernel
ker(A) of A. By definition of ker(A), cn
a vector v of R is in ker(A) if and only if
0
.. dm ); let us be
LA (v) = 0Rdm (where 0 d
R m = . denotes the zero vector of R
0
more precise. Let then
x1
x2
x3
v=
..
.
xn
LA (v) = 0Rd
m
146 SECTION 11
0
(where 0Rd
m = ... denotes the zero vector of R
d m ) is equivalent to the
0
equality
cn be
we want to show that α1 , α2 , · · · , αn must then all be zero. Now let v ∈ R
given by:
α1
α2
v = α3 .
..
.
αn
The equality
can be written
LA (v) = 0Rd
m
which implies v ∈ ker(A) (by definition of ker(A)), and since we have assumed
that ker(A) = {0Rcn }, it follows that v = 0Rcn , i.e.
α1 0
α2 0
α3 0
= ,
.. ..
. .
αn 0
Theorem 26. Let A ∈ Mm,n (R) (i.e. A is a real m×n matrix). We have ker(A) =
{0Rcn } if and only if the column vectors of A are linearly independent.
We close this section with some important terminology and an important the-
orem which is a direct adaptation of the Rank-Nullity Theorem of Lecture
9:
rank(A) + nullity(A) = n,
Proof. The result follows directly from the Rank-Nullity Theorem of Lecture
9, taking into account the fact that:
PROBLEMS:
1. For each of the following choices for the matrix A, do the following:
• Specify
the linear
transformation LA that it defines (for example, if
1 0 c2 → R
c2 defined by
A= , then LA is the mapping LA : R
2 1
x1 x1 x1 c2 ,
LA ( )= , for every ∈R
x2 2x1 + x2 x2
• specify its range Im(A) as well as the dimension of Im(A),
148 SECTION 11
5 0 0 0
0 5 0 0
(p) A = .
0 0 5 0
0 0 0 5
3 0 0 0
(q) A = .
0 −2 0 0
0 3 0 0
(r) A = .
0 0 −1 0
1 0 2
2 −1 4
(s) A =
−2
.
0 −4
−1 1 1
1 0 2
2 −1 4
(t) A =
−2
.
0 −4
−1 1 −2
1 2 2
2 4 4
(u) A = .
−2 −4 −4
−1 −2 −2
1 2 2
2 4 4
(v) A =
−2
.
−4 −4
−1 2 −2
1 0 2 −1
(x) A = 3 0 6 −3 .
−2 0 −4 2
1 1 2 −1
(y) A = 3 0 6 −3 .
−2 0 −4 2
1 1 2 −1
(z) A = 3 0 6 3 .
−2 0 −4 2
150 SECTION 11
Section 12
Study Topics
• Systems of Linear Transformations Revisited
• Solving Systems of Linear Equations
151
152 SECTION 12
where the real numbers a1,1 , a1,2 , · · · , am,n and b1 , b2 , · · · , bm are given real
numbers, and we wish to solve for the real numbers x1 , x2 , · · · , xn .
Can we write this system in terms of matrices and vectors ? Before answering
this question, you may ask what the point of writing such a system in terms of
matrices and vectors would be. The answer to this last question is quite simple:
If we are able to write that system of equations in terms of matrices and vectors,
then with the linear algebra knowledge we have accumulated so far, we may be
able to say something interesting about that system of equations: For example,
whether or not it has a solution, how many solutions it does have, and so on,
and even beyond this, in case it does have a solution, how to go about finding
that solution. So let us not waste any time and let us try to express that system
of linear equations using matrices and vectors.
Let then A ∈ Mm,n (R) be the real m × n matrix defined by
a1,1 a1,2 · · · a1,n
a2,1 a2,2 · · · a2,n
A= . .. .. ,
.. . ··· .
am,1 am,2 ··· am,n
d
and let b ∈ R m be the real vector defined by
b1
b2
b = . .
. .
bm
It is easy to verify that the real numbers x1 , x2 , · · · , xn satisfy the above system
of linear equations if and only if the vector v ∈ R cn defined by
x1
x2
v = . .
..
xn
satisfies the vector equation (so called because the unknown is now a vector
cn ):
in R
LA (v) = b.
153
Note that in keeping with our simplified notation (i.e. writing Av instead of
LA (v)), we can write this last equation simply as
Av = b.
Before going further, let us now see how easily we can prove that our original
system of linear equations (E) can have only 0, 1 or infinitely many solutions.
We already proved this in Lecture 6, but let us prove it again here.
From the discussion above, we need only show therefore that the equation
LA (v) = b, where A ∈ Mm,n (R) and b ∈ R d m are given, and where we wish
to solve for the unknown vector v ∈ R cn , can have only 0, 1 or infinitely many
solutions. We have the following cases:
• Let A ∈ Mm,n (R); under what conditions on A does there exist, for
d
each choice of b ∈ R m , a unique solution to the equation L (v) = b ?
A
LA (v) = b.
Proof. (a) Assume that for each b ∈ R d m there exists a unique v ∈ R cn which
satisfies LA (v) = b; since for each b ∈ R d m the equation L (v) = b
A
d m
has a solution, it follows that Im(A) = R , which implies that m ≤ n;
furthermore, since that solution is unique, it follows that ker(A) = {0Rcn },
which implies that m ≥ n. Since we have both m ≤ n and m ≥ n, it
follows that m = n.
(b) Assume now that m = n and that the column vectors of A are linearly
independent. Since the n column vectors A;1 , A;2 , · · · , A;n of A are ele-
ments of R cn (since m = n), since they are assumed to be linearly inde-
pendent, and since R cn has dimension n, it follows (from Lecture 8) that
{A;1 , A;2 , · · · , A;n } is also a generating set for R cn , i.e. that Im(A) = R cn .
Hence, for each b ∈ R d m , there is at least one vector v ∈ R cn which satisfies
LA (v) = b. Furthermore, since the column vectors of A are assumed to
be linearly independent, it follows (from Lecture 10) that ker(A) = {0Rcn },
and therefore (again from lecture 10), for each b ∈ R d m , there is a unique
Now that, thanks to this last theorem, we have the answer to the question we
posed earlier, let us examine how we could go about solving an equation of the
form LA (v) = b, with A a square n × n real matrix.
NOTE: For simplicity, we will denote the set of all square n × n real matrices
simply by Mn (R).
cn . We wish to find v ∈ R
Let then A ∈ Mn (R), and let b ∈ R cn which satisfies
LA (v) = b.
Assume that there exist some matrix B ∈ Mn (R) such that LB (LA (v)) = v
cn , i.e., such that LB ◦ LA (v) = v for all v ∈ R
for all v ∈ R cn ; how can the
155
knowledge of this matrix B help us solve our equation ? Well, if some vector of
cn , say w, satisfies the equation
R
LA (w) = b,
then, applying LB to both sides, it must also satisfy the equation
LB (LA (w)) = LB (b),
cn , it follows
and since we assumed B was such that LB (LA (v)) = v for all v ∈ R
that w must also satisfy
w = LB (b),
i.e. w must be equal to LB (b), i.e. the solution to the equation LA (v) = b is
given by v = LB (b); in other words, we have solved the equation!
Let us examine this idea on an example before going further:
4 1
• Let A ∈ M2 (R) be the square 2 × 2 matrix defined by A = , and
3 1
c2 be the vector defined by b = b1 c2
let b ∈ R . We wish to find v ∈ R
b2
such that LA (v) = b. Note that since the column vectors of A are linearly
independent (as can be easily verified), it follows that for any b ∈ Rc2 , the
equation LA (v) = b has a unique solution.
Letnow B ∈ M2 (R)be the
1 −1 x1
square 2 × 2 matrix defined by B = . Let now v =
−3 4 x2
c2 . We have:
be any vector in R
LB ◦ LA (v) = LB (LA (v))
4x1 + x2
= LB ( )
3x1 + x2
(4x1 + x2 ) − (3x1 + x2 )
=
−3(4x1 + x2 ) + 4(3x1 + x2 )
x1
=
x2
= v,
Hence, thanks to the fact that we got our hands on that special matrix
B, we were able to solve our vector equation.
The natural questions one may ask at this point are:
• How did this matrix B show up ?
• How does one find such a matrix B in a systematic way ?
We will have the complete answer to these questions in the next few lectures.
Section 13
Study Topics
• Matrix Multiplication
157
158 SECTION 13
Before going further, let us pay attention to the fact that A has m rows and n
columns whereas B has n rows and p columns; in particular the number
of columns of A is equal to the number of rows of B.
We know that A defines the linear mapping
cn → R
LA : R d m,
cp → R
LB : R cn .
defined by:
cp .
LA ◦ LB (v) = LA (LB (v)), ∀v ∈ R
then, we obtain:
b1,1 x1 + · · · + b1,p xp
..
LB (v) = . ,
bn,1 x1 + · · · + bn,p xp
159
and hence
a1,1 (b1,1 x1 + · · · + b1,p xp ) + · · · + a1,n (bn,1 x1 + · · · + bn,p xp )
..
LA ◦ LB (v) = .
am,1 (b1,1 x1 + · · · + b1,p xp ) + · · · + am,n (bn,1 x1 + · · · + bn,p xp )
(a1,1 b1,1 + · · · + a1,n bn,1 )x1 + · · · + (a1,1 b1,p + · · · a1,n bn,p )xp
..
= .
(am,1 b1,1 + · · · + am,n bn,1 )x1 + · · · + (am,1 b1,p + · · · am,n bn,p )xp
= LC (v),
and let v be the real n × 1 matrix (i.e. the column vector) given by:
x1
x2
v = . ;
..
xn
which is nothing other than LA (v) !!!! This is why we shall often simply
write Av instead of LA (v) ...
(i.e. 1’s on the diagonal, and 0’s everywhere else), and let B ∈ Mn,p (R)
be a real n × p matrix. Since the number of columns of A (namely n) is
equal to the number of rows of B, it makes sense to consider the matrix
product AB. AB has n rows and p columns, just like B. Furthermore, a
simple computation shows that AB = B for any B ∈ Mn,p (R). For this
reason, (as we have seen before), this particular matrix A is called the
n × n identity matrix.
Note that since in this particular case the number of columns of B (namely
3) is equal to the number of rows of A, it also makes sense to compute
the matrix product BA, and we obtain after a simple computation that
BA is the 2 × 2 real matrix given by:
1 7
BA = .
1 5
Note that in this particular example where both AB and BA made sense,
AB and BA ended up being of different types: AB is 3 × 3 whereas BA is
2 × 2. The next example shows that even when AB and BA are of the
same type, they are not necessarily equal ...
5. Let now A and B both be real 2 × 2 matrices (so both matrix products
AB and BA make sense in this case) given by:
0 1 0 0
A= , B= .
0 0 1 0
whereas
0 0
BA = ,
0 1
which shows that AB is not equal to BA. So keep in mind that ma-
trix multiplication is NOT like ordinary multiplication of real
numbers! (Another way to put this is that matrix multiplication is not
commutative).
6. Let A ∈ Mm,n (R) be a real m × n matrix, and let B ∈ Mn,p (R) be
a real n × p matrix. Denote the p column vectors of B by (in order)
162 SECTION 13
B;1 , B;2 , · · · , B;p , i.e. we can write B = (B;1 B;2 · · · B;p ). It is easy to
verify that the matrix product AB has columns vectors given (in order)
by AB;1 , AB;2 , · · · , AB;p , i.e., we can write AB = (AB;1 AB;2 · · · AB;p ).
This is often useful in proofs ...
Now that we have seen the “mechanics” of matrix multiplication, let us examine
some of its properties:
(AB)C = A(BC),
A(B + C) = AB + AC
(A + B)C = AC + BC
NOTE: What statement (i) of the theorem says is that we can first multiply A
and B to obtain AB, and then multiply this with C to obtain (AB)C; we can
also first multiply B and C to obtain BC and then multiply A with BC to obtain
A(BC). Both give the same result. Hence, we can remove the parentheses and
simply write ABC instead of either (AB)C or A(BC).
We shall only prove (i); Properties (ii),(iii),(iv) can be directly verified by com-
putation and are left as an exercise to the reader.
Let us now see on some simple examples how we can apply matrix opera-
tions to the solution of systems
of linear equations:
1. Consider the system of linear equations given by
2x + 3y = 5
3x + y = 7,
where we wish to solve for the pair (x, y) of real numbers which satisfies
that system, if such a pair does exist. We can write this equation in matrix
form as:
2 3 x 5
= ,
3 1 y 7
2 3
and since the column vectors of the 2 × 2 matrix are linearly
3 1
independent we know that this system hasa unique solution. Consider
− 17 3
7
now the 2 × 2 matrix given by 3 . It is easy to verify that
7 − 72
1 3
−7 7 2 3 1 0
3 = ,
7 − 72 3 1 0 1
which is nothing other than the 2 × 2 identity matrix. Hence, multiplying
both sides of the equation
2 3 x 5
= ,
3 1 y 7
1 3
−7 7
by the matrix 3 , we obtain
7 − 72
1 3
1 3
−7 7 2 3 x −7 7 5
3 = ,
7 − 72 3 1 y 3
7 − 27 7
i.e., equivalently,
16
1 0 x 7
= 1 ,
0 1 y 7
and since
1 0 x x
= ,
0 1 y y
we finally obtain the solution to our system of equations as
16
x 7
= 1 .
y 7
164 SECTION 13
a1,1 x1 + · · · a1,n xn = b1
a2,1 x1 + · · · a2,n xn = b2
..
.
an,1 x1 + · · · an,n xn = bn
Av = b,
and, finally,
x1
x2
v= .. ,
.
xn
Av = b,
by B yields
BAv = Bb,
165
8 5 12 12 15 0 23 15 18 12 4 27 ,
8 5 12 12 15 0 23 15 18 12 4 27 ;
37 29 60 48 60 45 107 84 84 66 43 39 .
What we have done is “scramble” the original message; the code cannot
be easily recovered from the occurence statistics of the integers. But how
do we recover our original message ? Well, the party we are sending our
message to is aware that they need matrix B to recover the message. All
they need to do is break down the message they have received (i.e. the
1 × 12 matrix above) into blocks of length 2, and multiply each block on
the right by B. Since AB is the 2 × 2 identity matrix, the original message
can be recovered. For example, multiplying
the first block of length 2 of
the received message, i.e. 37 29 by matrix B on the right yields the
1 × 2 matrix 8 5 . Continuing this procedure for all subsequent 1 × 2
blocks yields the original message.
There are more sophisticated variations of this approach; we could for ex-
ample scramble the message even more by arranging the original message
as the following 2 × 6 matrix:
8 5 12 12 15 0
;
23 15 18 12 4 27
167
8 5
We would then take the first 2 × 2 block, i.e. the submatrix ,
23 15
37 29
multiply it on the right by A, obtaining . We would then do
107 84
the same with the next 2 × 2 block, and so on. We would then arrange
all the 2 × 2 blocks thus obtained as a 2 × 6 matrix which we would send
to the other party; the receiving party would then take the consecutive
2 × 2 blocks of this 2 × 6 matrix, multiply them on the right by B, thereby
recovering the original 2 × 6 matrix containing the message. Compared
to our original approach, this would scramble the original message even
more, making it even harder to recover the correspondence between the
letters and the integers.
Note that we could also scramble the original message more by taking A
and B to be 3 × 3 (or even larger) square matrices; the key properties of
A and B are that A should have integer entries (greater than or equal to
0) and that the product AB should be equal to the identity matrix. The
requirement that A have integer entries greater than or equal to 0 will
ensure that if the original unscrambled message is encoded using integers
greater than or equal to 0, so will the scrambled message. The requirement
that the product AB be equal to the identity matrix ensures that the
original message can be recovered from the scrambled message.
• Process Control
Consider two reservoirs containing a certain liquid, and assume the level
of the liquid in the
reservoirs
at time n (n = 0, 1, 2, . . .) is represented by
yn 2
the vector xn = ∈ R̂ , with yn representing the level in the first
zn
reservoir at time n, and zn representing the level in the second reservoir
at time n. Assume also that from time n to time n + 1, a certain fraction
of the liquid in both reservoirs is lost due to evaporation. We also assume
that a certain fraction of the liquid in the first reservoir is transferred to
the second reservoir through some channel in going from time n to time
n + 1. Finally, we assume that a certain amount of liquid is injected in
the first reservoir at times n = 0, 1, 2, . . .. This process can be modelled
mathematically as follows:
We have:
x1 = Ax0 + Bu0 ,
and hence
and finally
Since we have assumed that the reservoirs are initially empty, i.e. x0 is
the zero vector, we actually have:
The question
is now: Can we choose u0 , u1 , u2 such that x3 be equal to the
10
vector ? A simple calculation shows that the vectors B, AB, A2 B
20
span R̂2 ; hence not only can
we choose u0 , u1 , u2 to ensure that x3 is
10
equal to the vector – for example choosing u0 = 0, u1 = 100
20
and u2 = −65 yields the desired outcome – we can actually, by a proper
choice of u0 , u1 , u2 , make sure that x3 is any desired vector of R̂2 . In the
language of process control, the vector xn is called the state of the system
at time n, whereas un is called the control at time n.
PROBLEMS:
1. For each of the following choices for the matrices A and B, do the following:
−1 2 1 0
(i) A = , B= .
3 5 0 1
−1 2 0 1
(j) A = , B= .
3 5 1 0
−1 2
1 1
(k) A = 3 5 , B= .
0 1
7 4
−1 2
1 0
(l) A = 3 5 , B= .
0 1
7 4
−1 2
0 1
(m) A = 3 5 , B= .
1 0
7 4
−1 2 1 2 2
(n) A = 3 5 , B = 0 1 −3 .
7 4 −2 0 1
−1 2 1 0 0
(o) A = 3 5 , B= 0 2 0 .
7 4 0 0 −1
−1 2 0 1 0
(p) A = 3 5 , B = 1 5 7 .
7 4 0 0 1
−1 2 5 1 0 0
(q) A = 3 −2 5 , B = 0 1 0 .
7 0 4 0 0 1
−1 2 5 1 0 0
(r) A = 3 −2 5 , B = 0 0 1 .
7 0 4 0 1 0
−1 2 5 0 0 1
(s) A = 3 −2 5 , B = 0 1 0 .
7 0 4 1 0 0
−1 2 5 0 1 0
(t) A = 3 −2 5 , B = 1 0 0 .
7 0 4 0 0 1
−2 1
1 2 1
(u) A = 1 2 , B= .
0 3 5
3 0
−2 1 1 0 0
(v) A = 1 0 , B= 0 1 0 .
3 2 0 0 1
171
−1 3 0 1 0
(w) A = 0 0 , B = 1 0 0 .
1 5 0 0 1
−3 3
1
(x) A= 2 7 , B= .
1
−2 −5
−1 2 1
(y) A = 3 5 , B = 2 .
7 4 1
−5 −2 1 0 1 0
(z) A = −3 1 , B = 1 2 0 0 .
−2 0 0 3 0 1
16 12 5 1 19 5
M= .
0 3 1 12 12 27
M= 16 12 5 1 19 5 0 3 1 12 12 27 .
.5 .1 0
xn+1 = xn + un , n = 0, 1, 2, · · ·
.2 .7 1
(a) Assuming
that the reservoirs are empty at time n = 0, i.e. assuming
0
x0 = , can we suitably choose u0 to ensure that x1 will be
0
equal to
1
?
1
(b) Assuming
that the reservoirs are empty at time n = 0, i.e. assuming
0
x0 = , can we suitably choose u0 , u1 to ensure that x2 will be
0
equal to
1
?
1
1
(c) Assuming that the initial fluid levels are given by x0 = , can
0
2
we suitably choose u0 to ensure that x1 will be equal to ?
2
1
(d) Assuming that the initial fluid levels are given by x0 = , can
0
2
we suitably choose u0 , u1 to ensure that x2 will be equal to ?
2
1 1 0 1
xn+1 = 0 1 0 xn + 0 un , n = 0, 1, 2, · · ·
0 0 1 0
(a) Assuming
that
the reservoirs are empty at time n = 0, i.e. assuming
0
x0 = , can we suitably choose u0 to ensure that x1 will be
0
equal to
1
1 ?
1
(b) Assuming
that the reservoirs are empty at time n = 0, i.e. assuming
0
x0 = , can we suitably choose u0 , u1 to ensure that x2 will be
0
equal to
1
1 ?
1
(c) Assuming
that
the reservoirs are empty at time n = 0, i.e. assuming
0
x0 = , can we suitably choose u0 , u1 , u2 to ensure that x3 will
0
be equal to
1
1 ?
1
Section 14
Study Topics
• Invertible Square Matrices
• Determinant of a Square Matrix
175
176 SECTION 14
where the ai,j are the same real numbers as in the matrix A above and are
assumed given, the real numbers b1 , · · · , bn are assumed given, and we wish
to solve for the real numbers x1 , · · · , xn . We can write this system of linear
equations in matrix form as:
Av = b,
where A is the same matrix as defined above, b is the n × 1 matrix (i.e. column
vector) defined by
b1
b2
b = . ,
..
bn
1 0 ··· 0
0 1 ··· 0
where I = .. .. .. is the n × n identity matrix, we obtain, after
. . .
0 0 ··· 1
multiplying on the left both sides of the equation Av = b by the matrix B:
B(Av) = Bb,
(BA)v = Bb,
Iv = Bb,
v = Bb,
Remark 7. The reader may rightly wonder why we require in the above definition
that both AB and BA be equal to the n × n identity matrix I (since all we cared
about in solving our matrix equation was that BA be equal to I; we didn’t care
about AB). It is easy to show (but we won’t do it here) that AB = I implies
BA = I, and vice-versa. This is due to the fact that A and B are both square
matrices (and as an application of the rank-nullity theorem).
To recapitulate, we could have defined “A invertible” just by requiring the existence
of B ∈ Mn (R) such that BA = I; this would have been equivalent to the above
definition.
B = BI = B(AC) = (BA)C = IC = C,
Definition 31. Let A ∈ Mn (R) and assume A is invertible. The unique matrix
B ∈ Mn (R) which satisfies AB = BA = I is called the inverse of A and is
denoted A−1 .
(ii) it does not make any sense to talk about invertibility of an m × n matrix
with m 6= n.
2. Let A be the n × n zero matrix (i.e. the n × n matrix with all entries
equal to zero); since for any B ∈ Mn (R) the matrix products BA and
AB are equal to the zero matrix, it follows that there is no B ∈ Mn (R)
for which the matrix products BA and AB be equal to the n × n identity
matrix; hence the n × n zero matrix is not invertible.
5 1
3. Let A be the 2 × 2 matrix given by A = ; it is easy to verify that
4 1
1 −1
the 2 × 2 matrix given by B = satisfies AB = BA = I (I
−4 5
being the 2 × 2 identity matrix); hence A is invertible and B is its inverse.
Av = 0Rcn ,
and multiplying both sides of the above equation by the matrix B (the
inverse of A), we obtain:
(BA)v = 0Rcn ,
i.e.
v = 0Rcn ,
LB ◦ LA (v) = LA ◦ LB (v) = v, cn ,
∀v ∈ R
i.e.
BA = AB = I,
Before going further, let us summarize how we can apply this last ex-
tremely important theorem: Suppose you are given a real n × n matrix A;
then,
180 SECTION 14
which shows that the column vectors of A are not linearly independent;
hence, by the previous theorem, A is not invertible.
1 2 2
3. Consider the 3 × 3 real matrix A = 0 1 0 ; Is A invertible ?
2 3 5
Again, by the previous theorem, we have to see whether or not the column
vectors of A are linearly independent. Let A;1 , A;2 , A;3 denote the three
column vectors of A; it is easy to verify that {A;1 , A;2 , A;3 } is a linearly
independent subset of R c3 ; hence, it follows from the previous theorem
that A is invertible.
We now introduce an extremely important function of matrices, the deter-
minant. Before, however, we need to introduce some notation:
Definition 32. Let A ∈ Mn (R) be a square n × n matrix. Let i, j be integers ≥ 1
and ≤ n. We denote by [A]i,j the (n − 1) × (n − 1) matrix obtained from A by
deleting row i and column j.
Let us see consider a few examples:
2 1
1. Consider the real 2 × 2 matrix A = ; we have:
3 4
5 1 0
2. Consider the real 3 × 3 matrix A = −1 2 3 ; we have:
7 4 1
2 3 −1 3 −1 2
[A]1,1 = , [A]1,2 = , [A]1,3 = ,
4 1 7 1 7 4
1 0 5 0 5 1
[A]2,1 = , [A]2,2 = , [A]2,3 = ,
4 1 7 1 7 4
1 0 5 0 5 1
[A]3,1 = , [A]3,2 = , [A]3,3 = .
2 3 −1 3 −1 2
Now that we have seen how to use the determinant function to verify whether
or not a given real n × n matrix A is invertible, let us see how we can compute
the inverse A−1 of A. Before going further, we need an important definition:
In other words, the entry of M T on row i and column j is exactly the entry of
M on row j and column i. A very simple way to remember how to construct
the transpose AT of A is as follows: The first row of A becomes the first column
of AT , the second row of A becomes the second column of AT , ..., the last row
of A becomes the last column of AT .
Let us illustrate this on a few examples:
1 2 −3 7
1. Let A = be a real 2 × 4 matrix; its transpose AT is
−1 4 5 0
the real 4 × 2 matrix given by
1 −1
2 4
AT = −3
.
5
7 0
7 3 −3
2. Let A = −1 0 5 be a real 3 × 3 matrix; its transpose AT is the
2 2 1
real 3 × 3 matrix given by
7 −1 2
AT = 3 0 2
−3 5 1
1
3. Let A = be a real 2 × 1 matrix; its transpose AT is the real 1 × 2
2
matrix given by
AT = 1 2 .
Let us verify that we indeed have A−1 A = AA−1 = I (I being the 2 × 2 identity
matrix); we have:
1 d −b a b 1 ad − bc 0
A−1 A = =
ad − bc −c a c d ad − bc 0 ad − bc
1 0
= ,
0 1
a b 1 d −b 1 a b d −b
AA−1 = =
c d ad − bc −c a ad − bc c d −c a
1 ad − bc 0
=
ad − bc 0 ad − bc
1 0
= .
0 1
det(A) = det(AT ).
It follows from this last theorem that if A ∈ Mn (R) is invertible, then det(A−1 ) =
1
det(A) , i.e. the determinant of the inverse is the inverse of the determinant. In-
deed, since we have AA−1 = I (where I is the n × n identity matrix), we obtain
det(AA−1 ) = det(A) det(A−1 ) = det(I), and since we have already computed
that det(I) = 1, we obtain det(A) det(A−1 ) = 1, from which it follows that
1
det(A−1 ) = det(A) .
where ai,j denotes the entry of A on row i and column j, and (as before) [A]i,j
denotes the (n− 1)× (n− 1) matrix obtained from A by removing row i and column
j.
Remark 9. Note that the definition we gave for det(A) corresponds to the first
formula above for i = 1, i.e. “expansion along the first row of A”. What the first
formula in the above theorem says is that we can compute det(A) by “expanding
along any row of A” (not just the first one), and the second formula says that we
can compute det(A) by “expanding along any column of A”.
The previous result together with the previous remark lead us to the following
observation: Since we can compute det(A) by expanding along any row or col-
umn of A that we please, we may as well choose the row or column of A having
the most number of zero entries. We illustrate this idea on a few examples:
such a matrix is called lower triangular (since all the entries above the
diagonal are zero). Let us compute det(A). In line with the preceding
observation, it is to our advantage to expand along the row or column of
A having the most number of zero entries. Let us then compute det(A)
by expanding along the first row of A. Since the only possibly non-zero
entry on the first row of A is a1,1 , we obtain, upon applying the formula
for the determinant:
a2,2 0 ··· 0
a3,2 a3,3 · · · 0
det(A) = a1,1 det . . .. .. ,
.. .. . .
an,2 an,3 ··· an,n
and again computing this new determinant by row expansion along the
first row yields:
a3,3 0 ··· 0
a4,3 a4,4 ··· 0
det(A) = a1,1 a2,2 det . .. .. .. ,
.. . . .
an,3 an,4 ··· an,n
and iterating this procedure yields:
det(A) = a1,1 a2,2 a3,3 · · · an,n ,
i.e., in other words, the determinant of a lower triangular matrix A is
the product of the diagonal entries of A.
2. Let now A ∈ Mn (R) be given by
a1,1 a1,2 ··· a1,n
0 a2,2 ··· a2,n
A= . .. .. .. ;
. . . . .
0 0 ··· an,n
such a matrix is called upper triangular (since all the entries below the
diagonal are zero). Note that for this particular matrix, it would be most
advantageous to compute the determinant by expanding along the first
column. Note also that we can use our previous computation for lower
triangular matrices to compute det(A). Indeed, note that the transpose
AT of matrix A is lower-triangular, and hence det(AT ) is the product of
the diagonal entries of AT . We have also seen that transposing a matrix
does not change its determinant, i.e. det(AT ) = det(A). Putting all this
together, we obtain:
a1,1 0 ··· 0
a1,2 a2,2 · · · 0
det(A) = det(AT ) = det . . . .. = a1,1 a2,2 a3,3 · · · an,n ,
.. .. .. .
a1,n a2,n · · · an,n
187
such a matrix is called diagonal (since all the entries outside the diag-
onal are zero). Note that a diagonal matrix is both upper- and lower-
triangular. Hence, applying our previous results, we obtain:
where A;j denotes the j th column vector of A. If we now multiply the j th column
vector of A by some real number α (but only the j th column vector of A), we
obtain the matrix
A;1 · · · αA;j · · · A;n ,
which has exactly the same column vectors as A except for the j th one (which
has been multiplied by α. The following theorem indicates how the determinants
of those matrices are related:
Let now A ∈ Mn (R) be a real n × n matrix, and let α ∈ R; recall that αA is the
matrix obtained by multiplying each entry of A by α. How is det(αA) related
to det(A) ? The following theorem answers this question:
det(αA) = αn det(A).
we obtain:
αA = αA;1 ··· αA;j ··· αA;n ,
Remark 10. Note that it follows from the previous theorem that for any A ∈
Mn (R), we have:
where again A;i denotes the ith column vector of A and A;j its j th column
vector; what happens to the determinant if we interchange these two column
vectors, i.e. how is the determinant of the matrix
A;1 · · · A;j · · · A;i · · · A;n ,
and let à be the real n× n matrix obtained from A by interchanging column vectors
i and j (where i 6= j), i.e.,
à = A;1 · · · A;j · · · A;i · · · A;n .
Then:
det(Ã) = − det(A);
in other words, interchanging two distinct column vectors changes the sign of the
determinant.
189
Let now A ∈ Mn (R) be a real n × n matrix given (again in terms of its column
vectors) by
A = A;1 · · · A;j · · · A;n ,
which differs from A only in the j th column vector. How is the determinant of
this last matrix related to the determinant of A ?
Theorem 40. Let A ∈ Mn (R) be a real n × n matrix given in terms of its column
vectors by
A = A;1 · · · A;j · · · A;n ,
and let v ∈ Rcn be a column vector; consider the matrix obtained from A by adding
th
v to the j column vector of A, i.e. the matrix
A;1 · · · A;j + v · · · A;n ,
and let
A;1 ··· v ··· A;n ,
PROBLEMS:
1. For each of the following choices for the matrix A, establish whether or
not A is invertible in the following two distinct ways:
(i) Check whether or not the column vectors of A are linearly indepen-
dent.
(ii) Compute the determinant det(A) of A to see whether or not it is
non-zero.
−1 2
(a) A = .
2 1
190 SECTION 14
1 2
(b) A = .
2 0
1 2
(c) A = .
0 2
1 0
(d) A = .
0 1
0 1
(e) A = .
1 0
0 1
(f) A = .
0 1
1 1
(g) A = .
0 0
1 1
(h) A = .
0 1
1 1
(i) A = .
1 1
3 2
(j) A = .
6 4
−2 0
(k) A = .
0 −2
−3 0
(l) A = .
7 0
−1 2
(m) A = .
2 −4
1 0 3
(n) A = 2 1 2 .
3 3 3
3 2 1
(o) A = 2 0 1 .
0 0 1
1 2 1
(p) A = 2 4 2 .
3 5 2
1 1 0
(q) A = 2 7 0 .
3 −3 0
1 0 0
(r) A = 1 1 0 .
1 1 1
191
1 0 0
(s) A = 0 1 0 .
0 0 1
1 1 1
(t) A = 0 1 1 .
0 0 1
1 1 1
(u) A = 0 1 1 .
1 0 1
2 1 3
(v) A = 7 −1 −1 .
5 −2 2
0 1 0
(w) A = −1 2 3 .
−1 5 4
1 0 0 0
1 −1 0 0
(x) A= 1
.
2 −3 0
1 5 −7 −2
1 0 0 −1
1 −1 1 2
(y) A= 0 −2 2
.
3
1 0 0 −1
2 0 1 0
0 −2 5 0
(z) A= 0
.
1 0 2
1 1 −1 −1
Study Topics
• Eigenvalues and eigenvectors of an endomorphism
• Eigenvalues and eigenvectors of a square matrix
193
194 SECTION 15
d αt
(L(f ))(t) = (e ) = αeαt = αf (t) = (αf )(t),
dt
which shows that
L(f ) = αf.
Since f is not the zero vector of C ∞ (R; R), and since L(f ) = αf , we
conclude that f is an eigenvector of L associated to the eigenvalue α.
LA (v) = λv,
Av = λv.
λv − Av = 0Rcn ,
196 SECTION 15
i.e., that
(λI − A)v = 0Rcn ,
where I denotes the n × n identity matrix. But then, this last equality means
nothing other than that v is in the kernel ker(λI − A) of the matrix λI − A.
Since v is assumed to be distinct from the zero vector of R cn and since it lies
in the kernel ker(λI − A) of λI − A, it follows that ker(λI − A) must contain
more than just the zero vector of R cn ! But then this means that λI − A is not
invertible (otherwise ker(λI − A) would have contained only the zero vector of
cn ), i.e. that det(λI − A) = 0 !
R
Conversely, assume the real number λ ∈ R is such that det(λI −A) = 0; this then
implies that the kernel ker(λI − A) of the matrix λI − A is not equal to {0Rcn },
i.e. there does exist a vector v 6= 0Rcn in ker(λI − A); in other words, we have
both v 6= 0Rcn and (λI − A)v = 0Rcn , and this last equality implies Av = λv.
To recapitulate: We have shown that if λ ∈ R is such that det(λI − A) = 0,
then there exists a v ∈ Rcn such that v 6= 0 cn and Av = λv; this shows that λ
R
is an eigenvalue of A.
We have therefore proved the following important theorem:
Theorem 41. Let A ∈ Mn (R) be a real n × n matrix. Let λ ∈ R. We have:
• λ is an eigenvalue of A if and only if det(λI − A) = 0 (where again I
denotes the n × n identity matrix).
The importance of this theorem lies in that it will provide us a systematic way
to compute the eigenvalues of any n × n matrix. This theorem has the following
useful corollary:
Theorem 42. Let A ∈ Mn (R) be a real n × n matrix. Then 0 is an eigenvalue of
A if and only if A is not invertible.
Proof. 0 being an eigenvalue of A is equivalent (by the last theorem) to det(0I −
A) = 0, i.e. to det(−A) = 0, i.e. to (−1)n det(A) = 0, i.e to det(A) = 0,
which itself is equivalent to A being not invertible (by what we saw in the last
lecture).
Let now A ∈ Mn (R) be given by
a1,1 a1,2 ··· a1,n
a2,1 a2,2 ··· a2,n
A= . .. .. .
.. . ··· .
an,1 an,2 ··· an,n
We can then write, ∀λ ∈ R:
λ − a1,1 −a1,2 ··· −a1,n
−a2,1 λ − a2,2 ··· −a2,n
det(λI − A) = det .. .. .. .
. . ··· .
−an,1 −an,2 ··· λ − an,n
197
where the real numbers c0 , c1 , · · · , cn−2 , cn−1 are determined by the entries of
matrix A (i.e. by the ai,j ).
It is important to note the following key points:
2. the coefficients of the polynomial det(λI −A) are determined by the entries
of matrix A,
Due to the importance of the polynomial (in λ) defined by det(λI − A), we give
it a special name:
(ii) Find the roots of the characteristic polynomial of A, i.e. find all the
scalars λ which satisfy det(λI − A) = 0 ; these are the eigenvalues of A.
det(λI − A) = 0,
i.e. by solving
λ2 − 4λ + 3 = 0;
det(λI − I) = 0,
i.e. by solving
(λ − 1)n = 0
3. Let now A ∈ Mn (R) denote the n × n zero matrix (i.e. all entries of A
are zero); let us try to find the eigenvalues of A. Again, as detailed above,
all we have to do is:
(i) Compute the characteristic polynomial det(λI − A) of A,
(ii) find the roots of that characteristic polynomial, i.e. all λ which satisfy
det(λI − A) = 0.
The characteristic polynomial det(λI − A) of A is easily computed to be:
det(λI − A) = det(λI)
= λn det(I)
= λn .
det(λI − A) = 0,
i.e. by solving
λn = 0.
det(λI − A) = 0,
i.e. by solving
(λ − a)(λ − d)(λ − f ) = 0
The eigenvalues of A are then given by all the scalars λ which satisfy
det(λI − A) = 0,
λ2 + 1 = 0.
(λI − A)v = 0,
which implies
Av = λIv,
which implies
Av = λv,
c2 be defined by v1 = 1
Let v1 ∈ R . The above calculation shows
−1
that {v1 } is a generating set for the eigenspace of A associated with
eigenvalue 1; furthermore, it is easy to verify that {v1 } is a linearly inde-
pendent set. Hence, (v1 ) is a basis for the eigenspace of A associated
to the eigenvalue 1, and we conclude therefore that that eigenspace is one-
dimensional. In particular, v1 itself is an eigenvector of A associated
with eigenvalue 1. Let us verify this last claim: Clearly v1 6= 0, and we
have
2 1 1 2−1 1
Av1 = = = = v1 = 1v1 ,
1 2 −1 1−2 −1
as expected.
Let us now compute the eigenspace of A associated to the eigenvalue 3
of A. The eigenspace associated to the eigenvalue 3 of A is the vector
c2 given by:
subspace of R
3 − 2 −1 1 −1
ker(3I − A) = ker( ) = ker( ),
−1 3 − 2 −1 1
as expected.
203
Theorem 43. Let A ∈ Mn (R) and let λ be an eigenvalue of A. Then for any
integer k ≥ 1, λk is an eigenvalue of Ak .
cn , v 6= 0,
Proof. Since by assumption λ is an eigenvalue of A, there exists v ∈ R
such that Av = λv; hence:
(where we have used both the associativity and linearity of matrix multiplica-
tion) and upon iterating this procedure, we obtain
PROBLEMS:
3 2
(s) A = .
0 3
1 0 0
(t) A= 0 1 0 .
0 0 1
1 1 0
(u) A= 0 1 1 .
0 0 1
1 0 1
(v) A= 0 1 0 .
0 0 1
1 1 1
(w) A= 0 1 1 .
0 0 1
2 0 0
(x) A= 0 5 0 .
0 0 7
2 1 0
(y) A= 0 5 1 .
0 0 7
2 0 1
(z) A= 0 5 0 .
0 0 7
2 0 0 0
0 2 0 0
(z1) A=
0 0
.
7 1
0 0 0 7
2 1 0 0
0 2 0 0
(z2)
A= .
0 0 7 1
0 0 0 7
2 1 0 0
0 2 0 0
(z3) A=
0 0
.
7 0
0 0 0 7
2 0 0 0
0 2 0 0
(z4) A=
0 0
.
7 0
0 0 0 7
2 1 0 0
0 3 0 0
(z5) A=
0 0
.
4 1
0 0 0 6
206 SECTION 15
2 0 0 0
0 3 0 0
(z6) A = .
0 0 4 1
0 0 0 6
2 0 0 0
0 3 0 0
(z7) A =
0
.
0 4 0
0 0 0 6