Matrices A brief introduction Basilio Bona DAUIN – Politecnico di Torino Semester 1, 2014-15 B. Bona (DAUIN) Matrices Semester 1, 2014-15 1 / 44
Definitions Definition A matrix is a set of N real or complex numbers organized in m rows and n columns, with N = mn a 11 a 12 · · · a 1 n · · · a 21 a 22 a 2 n � a ij � A = ≡ i = 1 , . . . , m j = 1 , . . . , n · · · · · · · · · a ij · · · a m 1 a m 2 a mn A matrix is always written as a boldface capital letter viene as in A . To indicate matrix dimensions we use the following symbols A m × n A ∈ F m × n A m × n A ∈ F m × n where F = R for real elements and F = C for complex elements. B. Bona (DAUIN) Matrices Semester 1, 2014-15 2 / 44
Transpose matrix Given a matrix A m × n we define a transpose matrix the matrix obtained exchanging rows and columns · · · a 11 a 21 a m 1 a 12 a 22 · · · a m 2 A T n × m = . . . ... . . . . . . · · · a 1 n a 2 n a mn The following property holds ( A T ) T = A B. Bona (DAUIN) Matrices Semester 1, 2014-15 3 / 44
Square matrix A matrix is said to be square when m = n A square n × n matrix is upper triangular when a ij = 0 , ∀ i > j · · · a 11 a 12 a 1 n 0 a 22 · · · a 2 n A n × n = . . . ... . . . . . . · · · 0 0 a nn If a square matrix is upper triangular its transpose is lower triangular and viceversa 0 · · · 0 a 11 · · · a 12 a 22 0 A T n × n = . . . ... . . . . . . · · · a 1 n a 2 n a nn B. Bona (DAUIN) Matrices Semester 1, 2014-15 4 / 44
Symmetric matrix A real square matrix is said to be symmetric if A = A T , or A − A T = O In a real symmetric matrix there are at least n ( n + 1) independent 2 elements. If a matrix K has complex elements k ij = a ij + j b ij (where j = √− 1) its conjugate is K with elements k ij = a ij − j b ij . Given a complex matrix K , an adjoint matrix K ∗ is defined, as the T = K T conjugate transpose K ∗ = K A complex matrix is called self-adjoint or hermitian when K = K ∗ . Some textbooks indicate this matrix as K † or K H B. Bona (DAUIN) Matrices Semester 1, 2014-15 5 / 44
Diagonal matrix A square matrix is diagonal if a ij = 0 for i � = j 0 · · · 0 a 1 0 · · · 0 a 2 A n × n = diag( a i ) = . . . ... . . . . . . 0 0 · · · a n A diagonal matrix is always symmetric. B. Bona (DAUIN) Matrices Semester 1, 2014-15 6 / 44
Skew-symmetric matrix Skew-symmetric matrix A square matrix is skew-symmetric or antisymmetric if A + A T = 0 A = − A T → Given the constraints of the above relation, a generic skew-symmetric matrix has the following structure · · · 0 a 12 a 1 n − a 12 · · · 0 a 2 n A n × n = . . . ... . . . . . . − a 1 n − a 2 n · · · 0 In a skew-symmetric matrix there are at most n ( n − 1) non zero 2 independent elements. We will see in the following some important properties of the skew-symmetric 3 × 3 matrices. B. Bona (DAUIN) Matrices Semester 1, 2014-15 7 / 44
Block matrix It is possible to represent a matrix with blocks as · · · A 11 A 1 n · · · · · · A = A ij A m 1 · · · A mn where the blocks A ij have suitable dimensions. Given the following matrices A 11 · · · A 1 n A 11 O O A 11 O O A 2 = A 3 = A 1 = · · · · · · O A ij A ij O O A ij O · · · O O A mn A m 1 A mn O O A mn A 1 is upper block triangular, A 2 is lower block triangular, and A 3 is block diagonal B. Bona (DAUIN) Matrices Semester 1, 2014-15 8 / 44
Matrix algebra Matrices are elements of an algebra , i.e., a vector space together with a product operator. The main operations of this algebra are: product by a scalar , sum , and matrix product Product by a scalar · · · · · · a 11 a 12 a 1 n α a 11 α a 12 α a 1 n a 21 a 22 · · · a 2 n α a 21 α a 22 · · · α a 2 n α A = α = . . . . . . ... ... . . . . . . . . . . . . · · · · · · a m 1 a m 2 a mn α a m 1 α a m 2 α a mn Sum a 11 + b 11 a 12 + b 12 · · · a 1 n + b 1 n · · · a 21 + b 21 a 22 + b 22 a 2 n + b 2 n A + B = . . . ... . . . . . . a m 1 + b m 1 a m 2 + b m 2 · · · a mn + b mn B. Bona (DAUIN) Matrices Semester 1, 2014-15 9 / 44
Matrix sum Sum properties A + O = A A + B = B + A ( A + B ) + C = A + ( B + C ) A T + B T ( A + B ) T = The null (neutral, zero) element O takes the name of null matrix. The subtraction (difference) operation is defined using the scalar α = − 1: A − B = A + ( − 1) B B. Bona (DAUIN) Matrices Semester 1, 2014-15 10 / 44
Matrix product Matrix product The operation is performed using the well-known rule “ rows by columns ”: the generic element c ij of the matrix product C m × p = A m × n · B n × p is n � c ij = a ik b kj k =1 The bi-linearity of the matrix product is guaranteed, since it is immediate to verify that, given a generic scalar α , the following identity holds: α ( A · B ) = ( α A ) · B = A · ( α B ) B. Bona (DAUIN) Matrices Semester 1, 2014-15 11 / 44
Product Product properties A · B · C = ( A · B ) · C = A · ( B · C ) A · ( B + C ) = A · B + A · C ( A + B ) · C = A · C + B · C ( A · B ) T = B T · A T In general: the matrix product is non-commutative: A · B � = B · A , apart from particular cases; A · B = A · C does not imply B = C , apart from particular cases; A · B = O does not imply A = O or B = O , apart from particular cases. B. Bona (DAUIN) Matrices Semester 1, 2014-15 12 / 44
Identity matrix A neutral element wrt product exists and is called identity matrix, written as I n or simply I when no ambiguity arises; given a rectangular matrix A m × n the following identities hold A m × n = I m A m × n = A m × n I n Identity matrix · · · 1 0 0 · · · · · · 0 0 I = . . . ... . . . . . . 0 0 · · · 1 B. Bona (DAUIN) Matrices Semester 1, 2014-15 13 / 44
Trace Trace The trace of a square matrix A n × n is the sum of its diagonal elements n � tr ( A ) = a kk k =1 The matrix traces satisfies the following properties tr ( α A + β B ) = α tr ( A ) + β tr ( B ) tr ( AB ) = tr ( BA ) tr ( A ) = tr ( A T ) tr ( A ) = tr ( T − 1 AT ) for non singular T (see below) B. Bona (DAUIN) Matrices Semester 1, 2014-15 14 / 44
Determinant Once defined the cofactor, the determinant of a square matrix A can be defined “by row”, i.e., choosing a generic row i , n n a ik ( − 1) i + k det ( A ( ik ) ) = � � det ( A ) = a ik A ik k =1 k =1 or, choosing a generic column j , we have the definition “by column”: n n a kj ( − 1) k + j det ( A ( kj ) ) = � � det ( A ) = a kj A kj k =1 k =1 Since these definition are recursive and assume the computation of determinants of smaller order minors, it is necessary to define the determinant of a matrix 1 × 1 (scalar), that is simply det ( a ij ) = a ij . B. Bona (DAUIN) Matrices Semester 1, 2014-15 15 / 44
Properties of determinant det( A · B ) = det( A ) det( B ) det( A T ) = det( A ) det( k A ) = k n det( A ) if one makes a number of s exchanges between rows or columns of A , obtaining a new matrix A s , we have det( A s ) = ( − 1) s det( A ) if A has two equal or proportional rows/columns, we have det( A ) = 0 if A has a row or a column that is a linear combination of other rows or columns, we have det( A ) = 0 e upper or lower triangular, we have det( A ) = � n if A ` i =1 a ii if A is block triangular, with p blocks A ii on the diagonal, we have det( A ) = � p i =1 det A ii B. Bona (DAUIN) Matrices Semester 1, 2014-15 16 / 44
Singular matrix and rank A matrix A is singular if det( A ) = 0. We define the rank of matrix A m × n , the number ρ ( A m × n ), computed as the maximum integer such that at least a non zero minor D p exists. The following properties hold: ρ ( A ) ≤ min { m , n } if ρ ( A ) = min { m , n } , A is said to have full rank if ρ ( A ) < min { m , n } , the matrix does not have full rank and one says that there is a fall of rank ρ ( AB ) ≤ min { ρ ( A ) , ρ ( B ) } ρ ( A ) = ρ ( A T ) ρ ( AA T ) = ρ ( A T A ) = ρ ( A ) if A n × n and det A < n then it has no full rank B. Bona (DAUIN) Matrices Semester 1, 2014-15 17 / 44
Invertible matrix Given a square matrix A ∈ R n × n , it is invertible of nonsingular if an inverse matrix A − 1 n × n exists, such that AA − 1 = A − 1 A = I n The matrix is invertible iff ρ ( A ) = n , or rather it has full rank; this implies det( A ) � = 0. The inverse matrix can be computed as 1 A − 1 = det( A ) Adj ( A ) The following properties hold: ( A − 1 ) − 1 = A ; ( A T ) − 1 = ( A − 1 ) T . The inverse matrix, if exists, allows to compute the following matrix equation y = Ax obtaining the unknown x as x = A − 1 y . B. Bona (DAUIN) Matrices Semester 1, 2014-15 18 / 44
Recommend
More recommend