Matrix free essay sample

In linear algebra, an n-by-n (square) matrix A is called invertible or nonsingular or nondegenerate if there exists an n-by-n matrix B such that where In denotes the n-by-n identity matrix and the multiplication used is ordinary matrix multiplication. If this is the case, then the matrix B is uniquely determined by A and is called the inverse of A, denoted by A-I . It follows from the theory of matrices that if for square matrices A and B, then also Non-square matrices (m-by-n matrices for which m * n) do not have an inverse. However, in some cases such a matrix may have a left inverse or right inverse. If A is m-by-n and the rank of A is equal to n, then A has a left inverse: an n-by-m matrix B such that BA = l. If A has rank m, then it has a right inverse: an n-by-m matrix B such that AB = l. We will write a custom essay sample on Matrix or any similar topic specifically for you Do Not WasteYour Time HIRE WRITER Only 13.90 / page While the most common case is that of matrices over the real or complex numbers, all these definitions can be given for matrices over any commutative ring. A square matrix that is not invertible is called singular or degenerate. A square matrix is singular if and only if its determinant is O. Singular matrices are rare in the sense hat if you pick a random square matrix, it will almost surely not be singular. Matrix inversion is the process of finding the matrix B that satisfies the prior equation for a given invertible matrix A. Properties of invertible matrices Let A be a square n by n matrix over a field K (for example the field R of real numbers). Then the following statements are equivalent: A is invertible. A is row-equivalent to the n-by-n identity matrix In. A is column-equivalent to the n-by-n identity matrix In. A has n pivot positions. detA rank A = n. The equation Ax = O has only the trivial solution x = O (i. e. Null A = {O}) The equation Ax = b has exactly one solution for each b in Kn. The columns of A are linearly independent. The columns of A span Kn (i. e. Col A = Kn). The columns of A form a basis of Kn. The linear transformation mapping x to Ax is a bijection from Kn to Kn. There is an n by n matrix B such that AB = In. The transpose AT is an invertible matrix. The number O is not an eigenvalue of A. The matrix A can expressed as a finite product of elementary matrices. In general, a unit in that ring. A matrix that is its own inverse, i. e. and , is called an involution. Proof for matrix product rule If A1, A2, , An are nonsingular square matrices over a field, then It becomes evident why this is the case if one attempts to find an inverse for the product of the Ais from first principles, that is, that we wish to determine B such that where B is the inverse matrix of the product. To remove A1 from the product, we can then write which would reduce the equation to Likewise, then, from which simplifies to If one repeats the process up to An, the equation becomes but B is the inverse matrix, i. . so the property is established Methods of matrix inversion Gaussian elimination Gauss-Jordan elimination is an algorithm that can be used to determine whether a given matrix is invertible and to find the inverse. An alternative is the LIJ decomposition which generates an upper and a lower triangular matrices which are easier to invert. For special purposes, it may be convenient to invert mat rices by treating mn-by-mn matrices as m-by-m matrices of n-by-n matrices, and applying one or another formula recursively (other sized matrices can be padded out with dummy rows and columns). For other purposes, a variant of Newtons method may be convenient (particularly when dealing with families of related matrices, so inverses of arlier matrices can be used to seed generating inverses of later matrices). Analytic solution Writing the transpose of the matrix of cofactors, known as an adJugate matrix, can also be an efficient way to calculate the inverse of small matrices, but this recursive method is inefficient for large matrices. To determine the inverse, we calculate a matrix of cofactors: where IAI is the determinant of A, CiJ is the matrix of cofactors, and AT represents matrix to solve a system of linear equations; however, for a unique solution, it is necessary that the matrix involved be invertible. Decomposition techniques like LIJ decomposition are much faster than inversion, and various fast algorithms for special classes of linear systems have also been developed. Inversion of 22 matrices The cofactor equation listed above yields the following result for 22 matrices. Inversion of these matrices can be done easily as follows: [2] This is possible because l/(ad-bc) is the reciprocal of the determinant of the matrix in question, and the same strategy could be used for other matrix sizes. Blockwise inversion Matrices can also be inverted blockwise by using the following analytic inversion ormula: where A, B, C and D are matrix sub-blocks of arbitrary size. (A and D must, of course, be square, so that they can be inverted. Furthermore, this is true if and only if A and D-CA-I B are nonsingular [3] ).