Matrix representation of a linear map

Similar documents
Matrix representation of a linear map

Matrix Arithmetic. j=1

Elementary Row Operations on Matrices

Phys 201. Matrices and Determinants

CS100: DISCRETE STRUCTURES. Lecture 3 Matrices Ch 3 Pages:

Linear Algebra and Matrix Inversion

. a m1 a mn. a 1 a 2 a = a n

1 Matrices and matrix algebra

Kevin James. MTHSC 3110 Section 2.1 Matrix Operations

n n matrices The system of m linear equations in n variables x 1, x 2,..., x n can be written as a matrix equation by Ax = b, or in full

Matrices: 2.1 Operations with Matrices

Elementary maths for GMT

10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )

A FIRST COURSE IN LINEAR ALGEBRA. An Open Text by Ken Kuttler. Matrix Arithmetic

Section 9.2: Matrices.. a m1 a m2 a mn

We could express the left side as a sum of vectors and obtain the Vector Form of a Linear System: a 12 a x n. a m2

Matrix Multiplication

Matrices and Linear Algebra

Section 9.2: Matrices. Definition: A matrix A consists of a rectangular array of numbers, or elements, arranged in m rows and n columns.

Matrices Gaussian elimination Determinants. Graphics 2009/2010, period 1. Lecture 4: matrices

Prepared by: M. S. KumarSwamy, TGT(Maths) Page

Matrix Operations. Linear Combination Vector Algebra Angle Between Vectors Projections and Reflections Equality of matrices, Augmented Matrix

Linear Algebra March 16, 2019

Chapter 2. Ma 322 Fall Ma 322. Sept 23-27

What is the Matrix? Linear control of finite-dimensional spaces. November 28, 2010

Math 360 Linear Algebra Fall Class Notes. a a a a a a. a a a

Matrices BUSINESS MATHEMATICS

Kevin James. MTHSC 3110 Section 2.2 Inverses of Matrices

Matrices and systems of linear equations

CLASS 12 ALGEBRA OF MATRICES

Chapter 1: Systems of linear equations and matrices. Section 1.1: Introduction to systems of linear equations

Matrix Algebra & Elementary Matrices

MATRICES. a m,1 a m,n A =

A matrix over a field F is a rectangular array of elements from F. The symbol

ICS 6N Computational Linear Algebra Matrix Algebra

Matrices and Determinants

. =. a i1 x 1 + a i2 x 2 + a in x n = b i. a 11 a 12 a 1n a 21 a 22 a 1n. i1 a i2 a in

Review of linear algebra

Matrix & Linear Algebra

PH1105 Lecture Notes on Linear Algebra.

Knowledge Discovery and Data Mining 1 (VO) ( )

Matrices. Chapter Definitions and Notations

Math 3C Lecture 20. John Douglas Moore

Chapter 4 - MATRIX ALGEBRA. ... a 2j... a 2n. a i1 a i2... a ij... a in

Mathematics 13: Lecture 10

Math 4377/6308 Advanced Linear Algebra

CS 246 Review of Linear Algebra 01/17/19

MAT 2037 LINEAR ALGEBRA I web:

MAC Module 2 Systems of Linear Equations and Matrices II. Learning Objectives. Upon completing this module, you should be able to :

Linear Algebra Notes. Lecture Notes, University of Toronto, Fall 2016

Math 3108: Linear Algebra

Fall Inverse of a matrix. Institute: UC San Diego. Authors: Alexander Knop

Announcements Monday, October 02

Announcements Wednesday, October 10

Introduction. Vectors and Matrices. Vectors [1] Vectors [2]

Chapter 2 Notes, Linear Algebra 5e Lay

Topics. Vectors (column matrices): Vector addition and scalar multiplication The matrix of a linear function y Ax The elements of a matrix A : A ij

Matrix Algebra. Matrix Algebra. Chapter 8 - S&B

Foundations of Cryptography

MTH 464: Computational Linear Algebra

Chapter 1 Matrices and Systems of Equations

MATH2210 Notebook 2 Spring 2018

Matrices. In this chapter: matrices, determinants. inverse matrix

Finite Mathematics Chapter 2. where a, b, c, d, h, and k are real numbers and neither a and b nor c and d are both zero.

INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRICES

Linear Algebra. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey

Matrices. Math 240 Calculus III. Wednesday, July 10, Summer 2013, Session II. Matrices. Math 240. Definitions and Notation.

Instruction: Operations with Matrices ( ) ( ) log 8 log 25 = If the corresponding elements do not equal, then the matrices are not equal.

Determinants - Uniqueness and Properties

Graduate Mathematical Economics Lecture 1

Definition 2.3. We define addition and multiplication of matrices as follows.

1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )

Matrix operations Linear Algebra with Computer Science Application

ICS141: Discrete Mathematics for Computer Science I

Matrix Algebra Determinant, Inverse matrix. Matrices. A. Fabretti. Mathematics 2 A.Y. 2015/2016. A. Fabretti Matrices

Calculus II - Basic Matrix Operations

Vectors and matrices: matrices (Version 2) This is a very brief summary of my lecture notes.

1 Matrices and Systems of Linear Equations. a 1n a 2n

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

is a 3 4 matrix. It has 3 rows and 4 columns. The first row is the horizontal row [ ]

Lecture 3: Matrix and Matrix Operations

Matrix Algebra 2.1 MATRIX OPERATIONS Pearson Education, Inc.

[ Here 21 is the dot product of (3, 1, 2, 5) with (2, 3, 1, 2), and 31 is the dot product of

ELEMENTARY LINEAR ALGEBRA

Matrix Arithmetic. a 11 a. A + B = + a m1 a mn. + b. a 11 + b 11 a 1n + b 1n = a m1. b m1 b mn. and scalar multiplication for matrices via.

Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition

Introduction to Matrix Algebra

Stage-structured Populations

SUPPLEMENTARY NOTES: CHAPTER 1

DM559 Linear and Integer Programming. Lecture 3 Matrix Operations. Marco Chiarandini

Systems of Linear Equations. By: Tri Atmojo Kusmayadi and Mardiyana Mathematics Education Sebelas Maret University

Multiplying matrices by diagonal matrices is faster than usual matrix multiplication.

Jim Lambers MAT 610 Summer Session Lecture 1 Notes

Linear Systems and Matrices

Matrices and Vectors

* is a row matrix * An mxn matrix is a square matrix if and only if m=n * A= is a diagonal matrix if = 0 i

MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants.

MA 106 Linear Algebra

Undergraduate Mathematical Economics Lecture 1

Chapter 4. Matrices and Matrix Rings

Transcription:

Matrix representation of a linear map As before, let e i = (0,..., 0, 1, 0,..., 0) T, with 1 in the i th place and 0 elsewhere, be standard basis vectors. Given linear map f : R n R m we get n column vectors (of size m) viz., f (e 1 ),..., f (e n ). Place them side by side., Thus if f (e j ) = (f 1j, f 2j,..., f mj ) T, then we obtain an array of the form M f = f 11 f 12... f 1n f 21 f 22... f 2n.. f m1 f m2... f mn. Such an array is called a matrix with m rows and n columns or a matrix of size m n. In short, we write M f = (f ij ). We call M f to be the matrix associated with the linear map f : R n R m. 12/45

13/45 If we let Mat m n (R), or simply, M m,n, denote the set of all m n matrices with entries in R, then we obtain: M : L(n, m) M m,n f M f The map M is a one-to-one correspondence (i.e., one-one and onto map) between the set of all linear maps from R n to R m and the set of all m n matrices with entries in R. Note: Two matrices are equal if their sizes are the same and corresponding entries are the same. M is one-one (or injective), means M f = M g for some f, g L(n, m) f = g, whereas M is onto (or surjective) means that A Mat m n (R) A = M f for some f L(n, m).

13/45 If we let Mat m n (R), or simply, M m,n, denote the set of all m n matrices with entries in R, then we obtain: M : L(n, m) M m,n f M f The map M is a one-to-one correspondence (i.e., one-one and onto map) between the set of all linear maps from R n to R m and the set of all m n matrices with entries in R. Note: Two matrices are equal if their sizes are the same and corresponding entries are the same. M is one-one (or injective), means M f = M g for some f, g L(n, m) f = g, whereas M is onto (or surjective) means that A Mat m n (R) A = M f for some f L(n, m).

13/45 If we let Mat m n (R), or simply, M m,n, denote the set of all m n matrices with entries in R, then we obtain: M : L(n, m) M m,n f M f The map M is a one-to-one correspondence (i.e., one-one and onto map) between the set of all linear maps from R n to R m and the set of all m n matrices with entries in R. Note: Two matrices are equal if their sizes are the same and corresponding entries are the same. M is one-one (or injective), means M f = M g for some f, g L(n, m) f = g, whereas M is onto (or surjective) means that A Mat m n (R) A = M f for some f L(n, m).

13/45 If we let Mat m n (R), or simply, M m,n, denote the set of all m n matrices with entries in R, then we obtain: M : L(n, m) M m,n f M f The map M is a one-to-one correspondence (i.e., one-one and onto map) between the set of all linear maps from R n to R m and the set of all m n matrices with entries in R. Note: Two matrices are equal if their sizes are the same and corresponding entries are the same. M is one-one (or injective), means M f = M g for some f, g L(n, m) f = g, whereas M is onto (or surjective) means that A Mat m n (R) A = M f for some f L(n, m).

If we let Mat m n (R), or simply, M m,n, denote the set of all m n matrices with entries in R, then we obtain: M : L(n, m) M m,n f M f The map M is a one-to-one correspondence (i.e., one-one and onto map) between the set of all linear maps from R n to R m and the set of all m n matrices with entries in R. Note: Two matrices are equal if their sizes are the same and corresponding entries are the same. M is one-one (or injective), means M f = M g for some f, g L(n, m) f = g, whereas M is onto (or surjective) means that A Mat m n (R) A = M f for some f L(n, m). M preserves addition and scalar multiplication, i.e., M f +g = M f +M g and M αf = αm f, f, g L(n, m), α R. 13/45

Examples: 14/45

14/45 Examples: 1. For the identity map id : R n R n, 1 0... 0 0 1... 0 M id = I = I n :=.. 0 0... 1

14/45 Examples: 1. For the identity map id : R n R n, 1 0... 0 0 1... 0 M id = I = I n :=.. = (δ ij) 0 0... 1 wehre δ ij is the Kronecker delta defined by δ ij = 1 if i = j and δ ij = 0, if i j.

14/45 Examples: 1. For the identity map id : R n R n, 1 0... 0 0 1... 0 M id = I = I n :=.. = (δ ij) 0 0... 1 wehre δ ij is the Kronecker delta defined by δ ij = 1 if i = j and δ ij = 0, if i j. 2. The matrix of the linear map T : R 2 R 2 which interchanges coordinates, is

14/45 Examples: 1. For the identity map id : R n R n, 1 0... 0 0 1... 0 M id = I = I n :=.. = (δ ij) 0 0... 1 wehre δ ij is the Kronecker delta defined by δ ij = 1 if i = j and δ ij = 0, if i j. 2. The matrix of the linear map T : R 2 ( R 2 which ) 0 1 interchanges coordinates, is M T =. 1 0

14/45 Examples: 1. For the identity map id : R n R n, 1 0... 0 0 1... 0 M id = I = I n :=.. = (δ ij) 0 0... 1 wehre δ ij is the Kronecker delta defined by δ ij = 1 if i = j and δ ij = 0, if i j. 2. The matrix of the linear map T : R 2 ( R 2 which ) 0 1 interchanges coordinates, is M T =. 1 0 3. For the linear map µ α : R n R n corresponding to multiplication by a fixed scalar α, the associated matrix is

14/45 Examples: 1. For the identity map id : R n R n, 1 0... 0 0 1... 0 M id = I = I n :=.. = (δ ij) 0 0... 1 wehre δ ij is the Kronecker delta defined by δ ij = 1 if i = j and δ ij = 0, if i j. 2. The matrix of the linear map T : R 2 ( R 2 which ) 0 1 interchanges coordinates, is M T =. 1 0 3. For the linear map µ α : R n R n corresponding to multiplication by a fixed scalar α, the associated matrix is the scalar matrix diag(α,..., α) := (αδ ij )

15/45 4. Rotation of the plane by an angle θ defines a linear map R θ : R 2 R 2 (Linearity of R θ can be shown by the law of congruent triangles or by finding a formula for R θ (x, y) using polar coordinates and trigonometric identities). The matrix of R θ is about the ( ) cos θ sin θ M Rθ = sin θ cos θ

15/45 4. Rotation of the plane by an angle θ defines a linear map R θ : R 2 R 2 (Linearity of R θ can be shown by the law of congruent triangles or by finding a formula for R θ (x, y) using polar coordinates and trigonometric identities). The matrix of R θ is about the ( ) cos θ sin θ M Rθ = sin θ cos θ 5. The inclusion map R 2 R 3 that sends (x, y) R 2 to (x, y, 0) R 3 is clearly a linear map and its matrix is

15/45 4. Rotation of the plane by an angle θ defines a linear map R θ : R 2 R 2 (Linearity of R θ can be shown by the law of congruent triangles or by finding a formula for R θ (x, y) using polar coordinates and trigonometric identities). The matrix of R θ is about the ( ) cos θ sin θ M Rθ = sin θ cos θ 5. The inclusion map R 2 R 3 that sends (x, y) R 2 to (x, y, 0) R 3 is clearly a linear map and its matrix is 1 0 0 1 0 0 In a similar way, we can write down the (n + t) n matrix corresponding to the inclusion map R n R n+t.

Operations on matrices 16/45 Set of all matrices of size m n is denoted by M m,n

Operations on matrices 16/45 Set of all matrices of size m n is denoted by M m,n As before, for A = (a ij ), B = (b ij ) M m,n and α R, A + B := (a ij + b ij ); αa := (αa ij ) Matrix addition and scalar multiplicaion satisfies the usual properties. The zero element is the zero matrix 0, all of whose entries are 0 and the negative of a matrix A = (a ij ) is the matrix A = ( a ij ).

Operations on matrices 16/45 Set of all matrices of size m n is denoted by M m,n As before, for A = (a ij ), B = (b ij ) M m,n and α R, A + B := (a ij + b ij ); αa := (αa ij ) Matrix addition and scalar multiplicaion satisfies the usual properties. The zero element is the zero matrix 0, all of whose entries are 0 and the negative of a matrix A = (a ij ) is the matrix A = ( a ij ).

Operations on matrices 16/45 Set of all matrices of size m n is denoted by M m,n As before, for A = (a ij ), B = (b ij ) M m,n and α R, A + B := (a ij + b ij ); αa := (αa ij ) Matrix addition and scalar multiplicaion satisfies the usual properties. The zero element is the zero matrix 0, all of whose entries are 0 and the negative of a matrix A = (a ij ) is the matrix A = ( a ij ). Transpose operation introduced earlier can be extended to all matrices A T := (b ij ) M n,m where b ij := a ji Note that (αa + βb) T = αa T + βb T for any α, β R.

Operations on matrices 16/45 Set of all matrices of size m n is denoted by M m,n As before, for A = (a ij ), B = (b ij ) M m,n and α R, A + B := (a ij + b ij ); αa := (αa ij ) Matrix addition and scalar multiplicaion satisfies the usual properties. The zero element is the zero matrix 0, all of whose entries are 0 and the negative of a matrix A = (a ij ) is the matrix A = ( a ij ). Transpose operation introduced earlier can be extended to all matrices A T := (b ij ) M n,m where b ij := a ji Note that (αa + βb) T = αa T + βb T for any α, β R.

Operations on matrices Set of all matrices of size m n is denoted by M m,n As before, for A = (a ij ), B = (b ij ) M m,n and α R, A + B := (a ij + b ij ); αa := (αa ij ) Matrix addition and scalar multiplicaion satisfies the usual properties. The zero element is the zero matrix 0, all of whose entries are 0 and the negative of a matrix A = (a ij ) is the matrix A = ( a ij ). Transpose operation introduced earlier can be extended to all matrices A T := (b ij ) M n,m where b ij := a ji Note that (αa + βb) T = αa T + βb T for any α, β R. Products of linear maps is not linear (in fact not even defined), in general. But composites of linear maps are linear. How does this relate to matrices? 16/45

17/45 Composition of Linear Maps and Matrix Multiplication Question: Let f : R n R m, g : R p R n be linear maps. If A := M f and B := M g then M f g =?

Composition of Linear Maps and Matrix Multiplication Question: Let f : R n R m, g : R p R n be linear maps. If A := M f and B := M g then M f g =? f g is a map of R p to R m. For 1 j p, ( n ) (f g)(e j ) = f (g(e j )) = f b kj e k = = n b kj f (e k ) = k=1 k=1 ( n m ) b kj a ik e i k=1 ( m n ) a ik b kj e i i=1 k=1 Thus M f g = C, where C = (c ij ) is the m p matrix whose (i, j) th entry is c ij = n k=1 a ikb kj. We call C the product of the m n matrix A and the n p matrix B, and write C = AB. i=1 17/45

Basic Properties of Matrix Multiplication 18/45 The following properties are readily derived from the corresponding properties of linear maps and the one-to -one correspondence between linear maps and matrices. Assume that the matrices A, B, C are of appropriate sizes so that all relevant products are defined. 1. Associativity: A(BC) = (AB)C.

Basic Properties of Matrix Multiplication 18/45 The following properties are readily derived from the corresponding properties of linear maps and the one-to -one correspondence between linear maps and matrices. Assume that the matrices A, B, C are of appropriate sizes so that all relevant products are defined. 1. Associativity: A(BC) = (AB)C. 2. Right and Left Distributivity: A(B + C) = AB + AC, (B + C)A = BA + CA

Basic Properties of Matrix Multiplication The following properties are readily derived from the corresponding properties of linear maps and the one-to -one correspondence between linear maps and matrices. Assume that the matrices A, B, C are of appropriate sizes so that all relevant products are defined. 1. Associativity: A(BC) = (AB)C. 2. Right and Left Distributivity: A(B + C) = AB + AC, (B + C)A = BA + CA 3. Multiplicative identity: For any A M m,n, AI n = A and I m A = A In particular, if A is a n n matrix, then AI n = A = I n A. Remark: Matrix multiplication is not commutative in general, i.e., AB need not equal BA even when both the products are defined. Exercise: Find examples! 18/45

Basic Properties of Matrix Multiplication The following properties are readily derived from the corresponding properties of linear maps and the one-to -one correspondence between linear maps and matrices. Assume that the matrices A, B, C are of appropriate sizes so that all relevant products are defined. 1. Associativity: A(BC) = (AB)C. 2. Right and Left Distributivity: A(B + C) = AB + AC, (B + C)A = BA + CA 3. Multiplicative identity: For any A M m,n, AI n = A and I m A = A In particular, if A is a n n matrix, then AI n = A = I n A. 4. Transpose of Product: (AB) T = B T A T. Remark: Matrix multiplication is not commutative in general, i.e., AB need not equal BA even when both the products are defined. Exercise: Find examples! 18/45

Invertible Maps 19/45 Recall the set-theoretic definition of invertible functions. Let X, Y be sets. A function f : X Y is said to be invertible if there exists g : Y X satisfying g f = id X and f g = id Y. In this case the function g is unique (check!) and it is called the inverse of f and denoted by f 1.

Invertible Maps Recall the set-theoretic definition of invertible functions. Let X, Y be sets. A function f : X Y is said to be invertible if there exists g : Y X satisfying g f = id X and f g = id Y. In this case the function g is unique (check!) and it is called the inverse of f and denoted by f 1. For any function f : X Y, f is invertible f is one-one and onto. If f : R n R n is an invertible linear map, then f 1 is linear. [To see this, let y 1, y 2 R n. Then there are x 1, x 2 R n such that f (x i ) = y i, i.e., x i = f 1 (y i ) for i = 1, 2. Thus f 1 (y 1 + y 2 ) = f 1 (f (x 1 ) + f (x 2 )) and using the linearity of f, we find that this is equal to f 1 (f (x 1 + x 2 )) = x 1 + x 2 = f 1 (y 1 ) + f 1 (y 2 ). Similarly f 1 preserves scalar multiplication. 19/45

Invertible Matrices 20/45 A square matrix A of size n n is said to be invertible if there exists an n n matrix B such that AB = BA = I n. If such a matrix B exists, then it is unique.

Invertible Matrices 20/45 A square matrix A of size n n is said to be invertible if there exists an n n matrix B such that AB = BA = I n. If such a matrix B exists, then it is unique. Indeed, if C is another matrix such that CA = AC = I n, then C = CI n = C(AB) = (CA)B = I n B = B. We thus call B to be the inverse of A and denote it by A 1. Note that if f A, f B : R n R n are the linear maps associated with A, B resp., then B = A 1 implies that f B = f 1 A.

Invertible Matrices 20/45 A square matrix A of size n n is said to be invertible if there exists an n n matrix B such that AB = BA = I n. If such a matrix B exists, then it is unique. Indeed, if C is another matrix such that CA = AC = I n, then C = CI n = C(AB) = (CA)B = I n B = B. We thus call B to be the inverse of A and denote it by A 1. Note that if f A, f B : R n R n are the linear maps associated with A, B resp., then B = A 1 implies that f B = f 1 A. Basic Properties/Examples: If A 1, A 2 are invertible, then so is A 1 A 2.

Invertible Matrices 20/45 A square matrix A of size n n is said to be invertible if there exists an n n matrix B such that AB = BA = I n. If such a matrix B exists, then it is unique. Indeed, if C is another matrix such that CA = AC = I n, then C = CI n = C(AB) = (CA)B = I n B = B. We thus call B to be the inverse of A and denote it by A 1. Note that if f A, f B : R n R n are the linear maps associated with A, B resp., then B = A 1 implies that f B = f 1 A. Basic Properties/Examples: If A 1, A 2 are invertible, then so is A 1 A 2. What is its inverse?

Invertible Matrices 20/45 A square matrix A of size n n is said to be invertible if there exists an n n matrix B such that AB = BA = I n. If such a matrix B exists, then it is unique. Indeed, if C is another matrix such that CA = AC = I n, then C = CI n = C(AB) = (CA)B = I n B = B. We thus call B to be the inverse of A and denote it by A 1. Note that if f A, f B : R n R n are the linear maps associated with A, B resp., then B = A 1 implies that f B = f 1 A. Basic Properties/Examples: If A 1, A 2 are invertible, then so is A 1 A 2. What is its inverse? Ans: (A 1 A 2 ) 1 = A 1 2 A 1 1.

Invertible Matrices A square matrix A of size n n is said to be invertible if there exists an n n matrix B such that AB = BA = I n. If such a matrix B exists, then it is unique. Indeed, if C is another matrix such that CA = AC = I n, then C = CI n = C(AB) = (CA)B = I n B = B. We thus call B to be the inverse of A and denote it by A 1. Note that if f A, f B : R n R n are the linear maps associated with A, B resp., then B = A 1 implies that f B = f 1 A. Basic Properties/Examples: If A 1, A 2 are invertible, then so is A 1 A 2. What is its inverse? Ans: (A 1 A 2 ) 1 = A 1 2 A 1 1. Clearly I n and more generally, any diagonal matrix diag(a 1, a 2,..., a n ) with a i 0 for all i is invertible; diag(a 1,..., a n ) 1 = diag(1/a 1,..., 1/a n ). 20/45

Elementary Matrices 21/45 Consider a m m matrix E ij whose (i, j) th entry is 1 and all other entries are 0.

Elementary Matrices 21/45 Consider a m m matrix E ij whose (i, j) th entry is 1 and all other entries are 0. If we multiply a m n matrix A by E ij on the left, then what we get is a m n matrix whose i th row is equal to the j th row of A and all other rows are zero.

Elementary Matrices 21/45 Consider a m m matrix E ij whose (i, j) th entry is 1 and all other entries are 0. If we multiply a m n matrix A by E ij on the left, then what we get is a m n matrix whose i th row is equal to the j th row of A and all other rows are zero. In particular, if i j, then E ij E ij = 0 and so for any α, (I + αe ij )(I αe ij ) = I + αe ij αe ij α 2 E ij E ij = I.

Elementary Matrices 21/45 Consider a m m matrix E ij whose (i, j) th entry is 1 and all other entries are 0. If we multiply a m n matrix A by E ij on the left, then what we get is a m n matrix whose i th row is equal to the j th row of A and all other rows are zero. In particular, if i j, then E ij E ij = 0 and so for any α, (I + αe ij )(I αe ij ) = I + αe ij αe ij α 2 E ij E ij = I. On the other hand, if i = j, then (I + αe ii ) = diag(1,..., 1, 1 + α, 1,..., 1) and this is invertible if α 1.

Elementary Matrices Consider a m m matrix E ij whose (i, j) th entry is 1 and all other entries are 0. If we multiply a m n matrix A by E ij on the left, then what we get is a m n matrix whose i th row is equal to the j th row of A and all other rows are zero. In particular, if i j, then E ij E ij = 0 and so for any α, (I + αe ij )(I αe ij ) = I + αe ij αe ij α 2 E ij E ij = I. On the other hand, if i = j, then (I + αe ii ) = diag(1,..., 1, 1 + α, 1,..., 1) and this is invertible if α 1. Matrices of the type I + αe ij, with α R and i j or of the type I + αe ii with α 1 provide simple examples of invertible matrices whose inverse is of a similar type. These are two among 3 possible types of elementary matrices. 21/45

The third remaining type of elementary matrix and its basic propery is described in the following easy exercise. 22/45 Exercise: Given any i j, show that the square matrix T ij := I + E ij + E ji E ii E jj is precisely the matrix obtained from the identity matrix by interchanging the i th and j th rows. Also show that T ij T ij = I and deduce that T ij is invertible and T 1 ij = T ij. A square matrix is said to be elementary if it of the type T ij (i j) or I + αe ii (α 1) or I + αe ij (i j). Thanks to the discussion in the last slide and the exercise above, we have the following result.

The third remaining type of elementary matrix and its basic propery is described in the following easy exercise. 22/45 Exercise: Given any i j, show that the square matrix T ij := I + E ij + E ji E ii E jj is precisely the matrix obtained from the identity matrix by interchanging the i th and j th rows. Also show that T ij T ij = I and deduce that T ij is invertible and T 1 ij = T ij. A square matrix is said to be elementary if it of the type T ij (i j) or I + αe ii (α 1) or I + αe ij (i j). Thanks to the discussion in the last slide and the exercise above, we have the following result. Theorem Every elementary matrix is invertible and its inverse is an elementary matrix of the same type.