How To Matrix
How To Matrix
Vable
a m1 a m2
(D.2b)
a m1 a m2
468
M. Vable
The order of multiplication is important when two matrices are multiplied. In Equation (D.4a), matrix [A] is said to pre-multiply matrix [B] and matrix [B] is said to post-multiply matrix [A]. [C] = [A][B] (D.4a) In Equation (D.4a) the number of columns of matrix [A] must equal to number of rows of matrix [B]. If matrix [A] of size (m x n) pre-multiplies matrix [B] of size (n x p), the result is a matrix [C] of size (m x p). The elements of matrix [C] can be found from
c ij =
aik bkj
k=1
i = 1, 2 m j = 1, 2 p
(D.4b)
a 11 a 12 a 1n
a 21 a 22 a 2n
a m1 a m2 a mn (D.5)
a m1 a m2
The element aij of matrix [A] becomes element aji in the transposed matrix [A]T. A square matrix (same number of rows and columns) is said to be symmetric if the transpose of the matrix is the same as the original matrix as shown in Equation (D.6) Symmetric Matrix [A] = [A] (D.6) Equation (D.7) lists the rules that apply to transpose of matrices during addition and multiplications. ([A] + [B]) = [A] + [B]
T T T T
(D.7)
469
M. Vable
The minor Mij associated with a element aij is the determinant of the matrix in which the ith row and jth column have been removed. The determinant of a matrix can be found using Equation (D.9) where i is any row in the matrix or it can be found using Equation (D.10) where j is any column in the matrix.
A =
( 1 )
k=1
i+k
a ik M ik
(D.9)
A =
( 1 )
k=1
k+j
a kj M kj
(D.10)
If the determinant of a matrix is zero i.e., A = 0 then the matrix [A] is said to be singular. In a singular matrix either all rows are not independent or all columns are not independent.
(D.11)
By Cramers rule the jth unknown xj can be found by rst replacing the jth column by the right
470
M. Vable
hand side vector, taking the determinant of the resulting matrix, and then dividing by the determinant of the matrix [A] as shown in Equation (D.12) a 11 a 21 a 12 a 22 r1 r2 rj a 1n a 2n (D.12)
a n1 a n2 rn a nn x j = ----------------------------------------------------------------------------------------A
j = 1, 2 n
(D.13)
[A]{x} = {r } (D.14a) where, {x} represents the unknown vector with components xj and {r} represents the right hand side vector with components rj. By pre-multiplying by [A]-1 to both sides of Equation (D.14a) and using Equation (D.13) we obtain the unknown vector as shown in Equation (D.14b) [A] [A]{x} = [A] {r }
1 1
or
1
[I]{x} = [A] {r }
or (D.14b)
471