Systems of Differential Equations

6.1 Review: Matrices

Linear algebra, particularly the study of matrices, is fundamental in understanding and solving systems of differential equations. This section provides a focused overview of the key concepts in matrix theory that are essential for this purpose.

A. Matrix Definition and Notation

A matrix is a rectangular array of numbers, symbols, or expressions, arranged in rows and columns. The individual items in a matrix are called its elements or entries. A matrix is typically denoted by a capital letter (e.g., A, B, C). The element in the [asciimath]i[/asciimath]-th row and [asciimath]j[/asciimath]-th column of a matrix A is denoted as [asciimath]a_(ij)[/asciimath]. The dimensions of a matrix are given as [asciimath]"rows" xx "columns"[/asciimath]. For example, a matrix A with [asciimath]m[/asciimath] rows and [asciimath]n[/asciimath] columns is an [asciimath]mxxn[/asciimath] matrix.

 [asciimath]A = [[a_11, a_12 ,..., a_(1n)],[a_21,a_22, ..., a_(2n)], [vdots,vdots ,ddots, vdots] ,[a_(m1), a_(m2) , ..., a_(mn)] ]_(mxxn)[/asciimath]

B. Special Matrices

A row matrix has only one row and multiple columns, while a column matrix has one column and multiple rows. These are also known as row vectors and column vectors, respectively.

 [asciimath]x = [[a_11, a_12,...,a_(1n)]]_(1xxn)[/asciimath]                     [asciimath]y = [[a_11],[a_(21)], [vdots] ,[a_(m1)] ]_(mxx1)[/asciimath]

A matrix with the same number of rows and columns is called a square matrix. For example, matrix B is an [asciimath]nxxn[/asciimath] square matrix.

 [asciimath]B = [[a_11, a_12 ,..., a_(1n)],[a_21,a_22, ..., a_(2n)], [vdots,vdots ,ddots, vdots] ,[a_(n1), a_(n2) , ..., a_(n\n)] ]_(nxxn)[/asciimath]

In a diagonal matrix, the elements outside the main diagonal are all zero. The main diagonal is the set of elements [asciimath]a_(ij)[/asciimath] where [asciimath]i=j[/asciimath]. For example, matrix C is an [asciimath]nxxn[/asciimath] diagonal matrix.

 [asciimath]C = [[a_11, 0 ,..., 0],[0,a_22, ..., 0], [vdots, vdots,ddots, vdots] ,[0, 0 , ..., a_(n\n)] ]_(nxxn)[/asciimath]

The identity matrix is a special type of diagonal matrix where all the elements on the main diagonal are 1. It is denoted as [asciimath]I[/asciimath] or [asciimath]I_n[/asciimath] to indicate its size ([asciimath]nxxn[/asciimath]).

 [asciimath]I_n = [[1, 0 ,..., 0],[0,1, ..., 0], [vdots, vdots,ddots, vdots] ,[0, 0 , ..., 1] ]_(nxxn)[/asciimath]

The zero matrix is a matrix in which all elements are zero. It is denoted by [asciimath]0_(mxxn)[/asciimath]to indicate its dimensions.

 [asciimath]0_(mxxn) = [[0, 0 ,..., 0],[0,0, ..., 0], [vdots, vdots,ddots, vdots] ,[0, 0 , ..., 0] ]_(mxxn)[/asciimath]

C. Matrix Operation

Matrix Addition and Subtraction

Matrix addition and subtraction are elementary operations where matrices of the same dimension are added or subtracted element by element. If [asciimath]A=[a_(ij)][/asciimath] and [asciimath]B=[b_(ij)][/asciimath] are matrices of the same size, their sum [asciimath]C=A+B[/asciimath] is a matrix where each element [asciimath]c_(ij)=a_(ij)+b_(ij)[/asciimath].

These operations are commutative (i.e., [asciimath]A+B=B+A[/asciimath]) and associative (i.e., [asciimath](A+B)+C=A+(B+C)[/asciimath]).

Scalar Multiplication

Scalar multiplication involves multiplying every element of a matrix by a scalar (a constant number). If [asciimath]k[/asciimath] is a scalar and [asciimath]A=[a_(ij)][/asciimath], then [asciimath]kA[/asciimath] is a matrix where each element is [asciimath]ka_(ij)[/asciimath].

Scalar multiplication is distributive over matrix addition or subtraction (i.e., [asciimath]k(A+B)=kA+kB[/asciimath] ) and associative with respect to the multiplication of scalars (i.e., [asciimath]k(lA)=(kl)A[/asciimath]).

 

Example 6.1.1: Matrix Subtraction and Scalar Multiplication

Find matric [asciimath]C[/asciimath] where [asciimath]C=3A-B[/asciimath] given matrices [asciimath]A[/asciimath] and [asciimath]B[/asciimath] .

 [asciimath]A=[[-1,3],[0,7]][/asciimath]  and [asciimath]B=[[2,-4],[-2,5]][/asciimath]

Show/Hide Solution

 

Matrices [asciimath]A[/asciimath] and [asciimath]B[/asciimath] are the same size and thus can be subtracted.

 [asciimath]C=3A-B[/asciimath]

 [asciimath]C=3[[-1,3],[0,7]] -[[2,-4],[-2,5]][/asciimath]

We first multiply all entries of matrix A by 3.

 [asciimath]=[[-3,9],[0,21]] -[[2,-4],[-2,5]][/asciimath]

We then subtract the corresponding entries.

 [asciimath]=[[-3-2,9-(-4)],[0-(-2),21-5]][/asciimath]

 [asciimath]C=[[-5,13],[2,16]][/asciimath]

 

Try an Example

 

Matrix Multiplication

Matrix multiplication is only possible when the number of columns in the first matrix matches the number of rows in the second matrix. Consider two matrices [asciimath]A_(mxxn)[/asciimath]  and [asciimath]B_(nxxp)[/asciimath]. The product of these matrices is a new matrix [asciimath]C_(mxxp)[/asciimath], where the dimension of C is [asciimath]mxxp[/asciimath]. Each element of C is computed by taking the dot product of a corresponding row from A and a column from B. This computation for each element in the [asciimath]i[/asciimath]-th row and [asciimath]j[/asciimath]-th column of C is given by the formula

 [asciimath]c_{ij} = a_{i1}b_{1j} + a_{i2}b_{2j} + \cdots + a_{i\n}b_{nj}[/asciimath] (6.1.1)

Matrix multiplication is associative, meaning [asciimath](AB)C = A(BC)[/asciimath]. It is also distributive over addition, which implies [asciimath]A(B + C) = AB + AC[/asciimath]. However, it is not commutative, meaning [asciimath]AB[/asciimath] may not equal [asciimath]BA[/asciimath]).

Special cases in matrix multiplication include interactions with identity and zero matrices. Multiplying any matrix by an identity matrix of appropriate size leaves the matrix unchanged (i.e., [asciimath]AI = IA = A[/asciimath]). Any matrix multiplied by a zero matrix results in a zero matrix of appropriate dimensions.

 

Example 6.1.2: Matrix Multiplication

Compute matrix [asciimath]C=AB[/asciimath] given

 [asciimath]A=[[1,4,-1],[2,0,-5]][/asciimath] and [asciimath]B=[[-7,3,-1,0],[-5,1,4,3],[0,-2,1,2]].[/asciimath]

Show/Hide Solution

 

To compute the product of matrices A and B, [asciimath]AB[/asciimath], we first verify that multiplication is possible. Matrix A has dimensions [asciimath]2xx3[/asciimath], and matrix B has dimensions [asciimath]3xx4[/asciimath]. Since the number of columns in A (3) matches the number of rows in B (3), multiplication can be performed. The resulting matrix C will have dimensions [asciimath]2xx4[/asciimath].

We compute each entry of matrix C using Equation 6.1.1:

 [asciimath]c_{11} = (1)(-7)+(4)(-5)+(-1)(0)=-27[/asciimath]

 [asciimath]c_{12} = (1)(3)+(4)(1)+(-1)(-2)=9[/asciimath]

 [asciimath]c_{13} = (1)(-1)+(4)(4)+(-1)(1)=14[/asciimath]

 [asciimath]c_{14} = (1)(0)+(4)(3)+(-1)(2)=10[/asciimath]

 [asciimath]c_{21} = (2)(-7)+(0)(-5)+(-5)(0)=-14[/asciimath]

 [asciimath]c_{22} = (2)(3)+(0)(1)+(-5)(-2)=16[/asciimath]

 [asciimath]c_{23} = (2)(-1)+(0)(4)+(-5)(1)=-7[/asciimath]

 [asciimath]c_{24} = (2)(0)+(0)(3)+(-5)(2)=-10[/asciimath]

Therefore, the resulting matrix C is

 [asciimath]C=[[-27,9,14,10],[-14,16,-7,-10]][/asciimath]

 

Try an Example

 

D. Matrix Determinant

The determinant is a scalar value that is associated with every square matrix. It provides critical information about the matrix, such as its invertibility. The determinant of a matrix A is denoted as

 [asciimath]det(A)=|A|[/asciimath]

For a [asciimath]2xx2[/asciimath] matrix, the determinant is calculated as

 [asciimath]|[a_11,a_12],[a_21,a_22]|=a_11a_22-a_21a_12[/asciimath](6.1.2)

For larger square matrices, the determinant is typically calculated using the method of cofactor expansion. For instance, the determinant of a [asciimath]3xx3[/asciimath] matrix can be computed by expanding along any row or column. Expanding along the first row, the formula is

 [asciimath]|[a_11,a_12,a_13],[a_21,a_22,a_23] ,[a_31,a_32,a_33] |=[/asciimath] [asciimath]a_11|[a_22,a_23],[a_32,a_33]|-[/asciimath] [asciimath]a_12|[a_21,a_23],[a_31,a_33]|+[/asciimath] [asciimath]a_13|[a_21,a_22],[a_31,a_32]|[/asciimath](6.1.3)

Another approach to compute determinants, especially for large matrices, is to use row reduction to transform the matrix into an upper triangular form. The determinant is then the product of the diagonal elements.

 

Example 6.1.3: Find Determinant

Find the determinant of the given matrices.

 [asciimath]A=[[-5,-1],[3,2]][/asciimath] and [asciimath]B= [[2,4,7],[5,-3,8],[0,-1,3] ][/asciimath]

Show/Hide Solution

 

To find the determinant of matrix A, we use Formula 6.1.2.

 [asciimath]|A|=[/asciimath] [asciimath]|[-5,-1],[3,2]|=(-5)(2)-(3)(-1)=-7[/asciimath]

To find the determinant of matrix B, we use Formula 6.1.3.

 [asciimath]|B|=|[2,4,7],[5,-3,8] ,[0,-1,3] |=[/asciimath] [asciimath]2|[-3,8],[-1,3]|-[/asciimath] [asciimath]4|[5,8],[0,3]|+[/asciimath] [asciimath]7|[5,-3],[0,-1]|[/asciimath]

[asciimath]=2((-3)(3)-(-1)(8))-4((5)(3)-(0)(8))+[/asciimath] [asciimath]7((5)(-1)-(0)(-3))[/asciimath]

[asciimath]=2(-1)-4(15)+7(-5)[/asciimath]

[asciimath]=-97[/asciimath]

 

Try an Example

 

E. Matrix Inverse

The inverse of a square matrix [asciimath]A[/asciimath], denoted as [asciimath]A^(-1)[/asciimath], is a matrix that, when multiplied with [asciimath]A[/asciimath], yields the identity matrix.

 [asciimath]A A^-1=A^-1 A=I_n[/asciimath]

One common method to find a matrix inverse is to use the adjugate and determinant. The formula is

 [asciimath]A^{-1} = \frac{1}{det(A)} \text{adj}(A)[/asciimath]

where [asciimath]"adj"(A)[/asciimath] is the adjugate of [asciimath]A[/asciimath], calculated from the cofactors of [asciimath]A[/asciimath]. This method involves computing the determinant and then the cofactor matrix, which is then transposed to get the adjugate matrix. For a [asciimath]2xx2[/asciimath] matrix [asciimath]A=[[a_11,a_12],[a_21,a_22]][/asciimath], the inverse is given by

 [asciimath]A^-1=1/(det(A))[[a_22,-a_12],[-a_21,a_11]][/asciimath](6.1.4)

Another method for finding the inverse is the row reduction method, which involves augmenting the matrix [asciimath]A[/asciimath] with the identity matrix

 [asciimath][A|I][/asciimath]

and then performing row operations to transform [asciimath]A[/asciimath] into the identity matrix. The operations that transform A into [asciimath]I[/asciimath] will transform the augmented identity matrix into [asciimath]A^{-1}[/asciimath].

 [asciimath][I|A^-1][/asciimath]

This method is particularly useful for numerical calculations and for larger matrices.

If a matrix is invertible, its inverse is unique. A square matrix is invertible if and only if it is nonsingular, meaning its determinant is not zero. If the determinant of a matrix is zero, the matrix does not have an inverse, and it is referred to as a singular matrix.

 

Example 6.1.4: Find Inverse of 2 by 2 Matrix

Find the inverse of matrix A, provided it exists.

 [asciimath]A=[[4,6],[2,6]][/asciimath]

Show/Hide Solution

 

We first find the determinant of A to determine if it has an inverse.

[asciimath]|A|=(4)(6)-(2)(6)=12[/asciimath]

The determinant is nonzero, so the inverse exists. For a [asciimath]2xx2[/asciimath] matrix, the cofactor approach, Formula 6.1.4, is fairly simple.

 [asciimath]A^-1=1/12[[6,-6],[-2,4]][/asciimath]

[asciimath]=[[1/2,-1/2],[-1/6,1/3]][/asciimath]

 

Try an Example

 

 

Example 6.1.5: Find Inverse of 3 by 3 Matrix

Find the inverse of matrix A, provided it exists.

 [asciimath]A=[[1,0,3],[-3,1,-9],[0,2,1]][/asciimath]

Show/Hide Solution

 

To find the inverse of a [asciimath]3xx3[/asciimath] matrix, the row reduction method is more straightforward. To find the inverse of matrix [asciimath]A[/asciimath] using the row reduction method, we start by forming an augmented matrix with matrix [asciimath]A[/asciimath] and the [asciimath]3xx3[/asciimath] identity matrix [asciimath]I_3[/asciimath]. The goal is to use row operations to transform the left side of the augmented matrix (the first three columns) into the identity matrix.

 [asciimath][A|I]=[/asciimath] [asciimath][[1,0,3,|,1,0,0],[-3,1,-9,|,0,1,0],[0,2,1,|,0,0,1]][/asciimath]

Apply the row operation:

[asciimath]R2=R2+3R1[/asciimath] (Add 3 times the first row to the second row):

 [asciimath][[1,0,3,|,1,0,0],[0,1,0,|,3,1,0],[0,2,1,|,0,0,1]][/asciimath]

[asciimath]R3=R3-2R2[/asciimath] (Subtract 2 times the second row from the third row):

 [asciimath][[1,0,3,|,1,0,0],[0,1,0,|,3,1,0],[0,0,1,|,-6,-2,1]][/asciimath]

[asciimath]R1=R1-3R3[/asciimath] (Subtract 3 times the third row from the first row):

 [asciimath][[1,0,0,|,19,6,-3],[0,1,0,|,3,1,0],[0,0,1,|,-6,-2,1]][/asciimath]

Since we have successfully transformed the left side of the augmented matrix into the identity matrix, the inverse of matrix A exists and is given by the right side of the augmented matrix:

 [asciimath]A^-1=[[19,6,-3],[3,1,0],[-6,-2,1]][/asciimath]

 

Try an Example

 

F. Matrix Calculus

Differentiation and integration of matrices are important in the context of systems of linear differential equations, particularly in finding the solution to nonhomogeneous systems.

Matrix Differentiation

Differentiating a matrix with function entries involves taking the derivative of each element of the matrix individually. Consider matrix [asciimath]A(t)[/asciimath]  whose entries are a function of [asciimath]t[/asciimath].

 [asciimath]A(t) = [[a_11(t), a_12(t) ,..., a_(1n)(t)],[a_21(t),a_22(t), ..., a_(2n)(t)], [vdots,vdots ,ddots, vdots] ,[a_(m1)(t), a_(m2)(t) , ..., a_(mn)(t) ] ][/asciimath]

The derivative of [asciimath]A(t)[/asciimath] with respect to [asciimath]t[/asciimath], denoted as [asciimath]A'(t)[/asciimath] or [asciimath](dA)/(dt)[/asciimath], is a matrix of the same size where each entry is the derivative of the corresponding entry of [asciimath]A(t)[/asciimath].

 [asciimath]A'(t) = [[a'_11(t), a'_12(t) ,..., a'_(1n)(t)],[a'_21(t),a'_22(t), ..., a'_(2n)(t)], [vdots,vdots ,ddots, vdots] ,[a'_(m1)(t), a'_(m2)(t) , ..., a'_(mn)(t)] ][/asciimath]

The standard rules of differentiation, including the product rule, quotient rule, and chain rule, apply to each element of the matrix.

Matrix Integration

Integrating a matrix with function entries is similar to differentiation and is done element-wise. The integral of a matrix [asciimath]A(t)[/asciimath] over a variable [asciimath]t[/asciimath] is a matrix of the same size where each element is the integral of the corresponding element of [asciimath]A(t)[/asciimath].

 [asciimath]intA(t)dt = [[inta_11(t)dt, inta_12(t)dt ,..., inta_(1n)(t)dt],[inta_21(t)dt,inta_22(t)dt, ..., inta_(2n)(t)dt], [vdots,vdots ,ddots, vdots] ,[inta_(m1)(t)dt, inta_(m2)(t)dt , ..., inta_(mn)(t)dt ] ][/asciimath]

 

Example 6.1.6: Matrix Integration

Evaluate the integral of matrix [asciimath]A(t)[/asciimath] with respect to [asciimath]t[/asciimath].

 [asciimath]A(t) = [(4 e^(4t) , 3te^(-t^2)) , (t^2cos(-3t^3) , -5t^7)][/asciimath]

Show/Hide Solution

 

The integral of matrices is an element-wise operation.

[asciimath]intA(t)dt=[intA_(ij)(t)dt][/asciimath]

 [asciimath]=[(int4 e^(4t) dt, int3te^(-t^2)dt) , (intt^2cos(-3t^3)dt , int-5t^7dt)][/asciimath]

 [asciimath]=[(e^(4t) , -3/2e^(-t^2)) , (-1/9sin(-3t^3) , -5/8t^8)][/asciimath]

 

Try an Example

 

Section 6.1 Exercises

  1. Given [asciimath]A = [(0,-3,3),(3,4,-3)][/asciimath]  and [asciimath]B = [(-4,0,-4),(-3,0,0)][/asciimath], find matrix [asciimath]C=2 A -4 B[/asciimath].
    Show/Hide Answer

     [asciimath]C=[(16,-6,22),(18,8,-6)][/asciimath]

  2. Find the inverse of [asciimath]A = [(-8,3),(21,-8)][/asciimath].
    Show/Hide Answer

     [asciimath]A^-1=[(-8,-3),(-21,-8)][/asciimath]

  3. Find the inverse of [asciimath]A = [(3,2,0),(1,1,0),(-2,-4,1)].[/asciimath] 
    Show/Hide Answer

    [asciimath]A^-1=[(1,-2,0),(-1,3,0),(-2,8,1)][/asciimath]

  4. Given the matrices [asciimath]A=[(5,-5),(-2,-4)][/asciimath] and [asciimath]B=[(-2,1),(4,4)][/asciimath] , find their multiplication [asciimath]AB[/asciimath].
    Show/Hide Answer

     [asciimath]AB=[(-30,-15),(-12,-18)][/asciimath]

  5. Given the matrix

     [asciimath]A(t) = [(-2 e^(4t) , 3te^(-3t^2)) , (-5tsin(4t^2) , -7t^-5)][/asciimath]

    Evaluate the integral of [asciimath]A[/asciimath] with respect to [asciimath]t[/asciimath].

    Show/Hide Answer

    [asciimath]intAdt=[/asciimath] [asciimath][(-1/2 e^(4t) , -1/2e^(-3t^2)) , (5/8cos(4t^2) , 7/4t^-4)][/asciimath]

License

Icon for the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License

Differential Equations Copyright © 2024 by Amir Tavangar is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License, except where otherwise noted.

Share This Book