Problem 36 - Advanced Engineering Mathematics

Question

Solving a nonhomogeneous linear system \(\mathbf{X}^{\prime}=\mathbf{A} \mathbf{X}+\mathbf{F}(t)\) by variation of parameters when \(\mathbf{A}\) is a \(3 \times 3\) (or larger) matrix is almost an impossible task to do by hand. Consider the system \(\mathbf{X}^{\prime}=\left(\begin{array}{rrrr}2 & -2 & 2 & 1 \\ -1 & 3 & 0 & 3 \\ 0 & 0 & 4 & -2 \\ 0 & 0 & 2 & -1\end{array}\right) \mathbf{X}+\left(\begin{array}{c}t e^{t} \\ e^{-t} \\ e^{2 t} \\\ 1\end{array}\right)\) (a) Use a CAS or linear algebra software to find the eigenvalues and eigenvectors of the coefficient matrix. (b) Form a fundamental matrix \(\boldsymbol{\Phi}(t)\) and use the computer to find \(\boldsymbol{\Phi}^{-1}(t)\) (c) Use the computer to carry out the computations of \(\mathbf{\Phi}^{-1}(t) \mathbf{F}(t), \int \mathbf{\Phi}^{-1}(t) \mathbf{F}(t) d t, \mathbf{\Phi}(t) \int \mathbf{\Phi}^{-1}(t) \mathbf{F}(t) d t, \mathbf{\Phi}(t) \mathbf{C}\) and \(\boldsymbol{\Phi}(t) \mathbf{C}+\int \boldsymbol{\Phi}^{-1}(t) \mathbf{F}(t) d t\), where \(\mathbf{C}\) is a column matrix of constants \(c_{1}, c_{2}, c_{3}\), and \(c_{4}\). (d) Rewrite the computer output for the general solution of the system in the form \(\mathbf{X}=\mathbf{X}_{c}+\mathbf{X}_{p}\), where \(\mathbf{X}_{c}=c_{1} \mathbf{X}_{1}+c_{2} \mathbf{X}_{2}+c_{3} \mathbf{X}_{3}+c_{4} \mathbf{x}_{4}\)

Step-by-Step Solution

Answer

The general solution is \( \mathbf{X} = \mathbf{X}_c + \mathbf{X}_p \), using \( \boldsymbol{\Phi}(t) \) and \( \integral \boldsymbol{\Phi}^{-1}(t) \mathbf{F}(t) dt \).

Step 1: Find Eigenvalues and Eigenvectors

To start, use a computer algebra system (CAS) to find the eigenvalues and eigenvectors of the matrix \( \mathbf{A} = \begin{pmatrix} 2 & -2 & 2 & 1 \ -1 & 3 & 0 & 3 \ 0 & 0 & 4 & -2 \ 0 & 0 & 2 & -1 \end{pmatrix} \). Compute the characteristic equation \( \det(\mathbf{A} - \lambda \mathbf{I}) = 0 \) and solve for \( \lambda \). Find corresponding eigenvectors for each eigenvalue.

Step 2: Form the Fundamental Matrix \( \boldsymbol{\Phi}(t) \)

Using the eigenvectors found in Step 1, form the fundamental matrix \( \boldsymbol{\Phi}(t) \). Each column of \( \boldsymbol{\Phi}(t) \) is built from the solution of the homogeneous part of the system formed by the eigenvectors. Compute \( \boldsymbol{\Phi}(t) \) using \( e^{\lambda t} \) terms.

Step 3: Compute \( \boldsymbol{\Phi}^{-1}(t) \)

Find the inverse of the fundamental matrix \( \boldsymbol{\Phi}^{-1}(t) \) using a CAS. The inverse is required to apply the method of variation of parameters.

Step 4: Compute \( \boldsymbol{\Phi}^{-1}(t) \mathbf{F}(t) \)

Using \( \boldsymbol{\Phi}^{-1}(t) \), calculate the product with the vector \( \mathbf{F}(t) = \begin{pmatrix} t e^{t} \ e^{-t} \ e^{2t} \ 1 \end{pmatrix} \). This step involves matrix multiplication, which can be done using a computer.

Step 5: Integrate \( \int \boldsymbol{\Phi}^{-1}(t) \mathbf{F}(t) dt \)

Compute the integral of \( \boldsymbol{\Phi}^{-1}(t) \mathbf{F}(t) \) with respect to \( t \) to find the particular solution integral component. Use a CAS to handle potentially complicated integrals.

Step 6: Construct the General Solution

Form the general solution \( \mathbf{X}(t) = \boldsymbol{\Phi}(t)\mathbf{C} + \boldsymbol{\Phi}(t)\int \boldsymbol{\Phi}^{-1}(t) \mathbf{F}(t) dt \). Here, \( \mathbf{C} = \begin{pmatrix} c_1 \ c_2 \ c_3 \ c_4 \end{pmatrix} \) contains arbitrary constants.

Step 7: Rewrite in Standard Form

Rewrite the solution in the form \( \mathbf{X} = \mathbf{X}_c + \mathbf{X}_p \), where \( \mathbf{X}_c = c_1 \mathbf{X}_1 + c_2 \mathbf{X}_2 + c_3 \mathbf{X}_3 + c_4 \mathbf{X}_4 \) is the complementary solution, and \( \mathbf{X}_p \) is the particular solution from the integral.

Key Concepts

Nonhomogeneous Linear System

When dealing with a nonhomogeneous linear system of differential equations, we are looking at an equation of the form \( \mathbf{X}^{\prime} = \mathbf{A} \mathbf{X} + \mathbf{F}(t) \). This equation involves a matrix \( \mathbf{A} \) and a vector function \( \mathbf{F}(t) \) that adds a non-zero component, differentiating it from a homogeneous system.
The challenge arises with the solution of these systems as the nonhomogeneous part \( \mathbf{F}(t) \) interacts with the system's intrinsic dynamics defined by \( \mathbf{A} \). Variation of parameters is a technique often used to solve such systems. Here, the solution is typically expressed as the sum of the complementary solution, which solves \( \mathbf{X}^{\prime} = \mathbf{A} \mathbf{X} \), and a particular solution that addresses the nonhomogeneous part.
Variation of parameters utilizes the inverse of a fundamental matrix to construct the particular solution. This can often be complex and involves robust calculations, making the use of computational tools like a computer algebra system essential for solving larger systems, such as those with a 3x3 matrix or bigger.

Eigenvalues and Eigenvectors

The process of finding eigenvalues and eigenvectors is crucial for analyzing the system defined by the matrix \( \mathbf{A} \). Eigenvalues are special numbers associated with a matrix that give us insight into the system's behavior over time.
To find them, we solve the characteristic equation \( \det(\mathbf{A} - \lambda \mathbf{I}) = 0 \). Each solution \( \lambda \) is an eigenvalue, and for each eigenvalue, there is an associated eigenvector that satisfies \( (\mathbf{A} - \lambda \mathbf{I})\mathbf{v} = 0 \). This vector remains a scaled version of itself under the transformation \( \mathbf{A} \).
These eigenvectors form the fundamental building blocks of the system's solutions. Together, they help form the fundamental matrix, \( \boldsymbol{\Phi}(t) \), used extensively in constructing both the complementary and particular solutions for the nonhomogeneous systems.

Fundamental Matrix

The fundamental matrix, denoted \( \boldsymbol{\Phi}(t) \), is a central element in solving linear differential systems. It is constructed using the eigenvectors of the coefficient matrix and encompasses the behavior of the homogeneous system.
The columns of this matrix are solutions derived from the eigenvectors \( \mathbf{v}_1, \mathbf{v}_2, ..., \mathbf{v}_n \) and their corresponding eigenvalues \( \lambda_1, \lambda_2, ..., \lambda_n \). Each column in the fundamental matrix is tackled using terms like \( e^{\lambda_i t} \), forming solutions that account for how the system evolves over time.
This matrix not only gives us the structure for the complementary solution but also plays a crucial role in calculating the inverse \( \boldsymbol{\Phi}^{-1}(t) \). The inverse is crucial for applying the variation of parameters method, used to obtain the particular solution of the nonhomogeneous system.

Computer Algebra System

A Computer Algebra System (CAS) is a powerful tool when solving complex mathematical problems, especially those involving large matrices or intricate functions. For differential equations like those in the nonhomogeneous linear systems, a CAS can handle operations that would otherwise be tremendously labor-intensive.
Using a CAS, students can quickly calculate the eigenvalues and eigenvectors of matrices, perform matrix inversion, and integrate complex functions. These tasks, though possible by hand, can be prone to error and exceedingly time-consuming
In our scenario, the CAS aids in managing the heavy computational demand of the variation of parameters. It provides a means to find the eigenstructure, build the fundamental matrix, compute its inverse, and evaluate the integrals essential for deriving the particular solution, all with much more reliability and efficiency than manual calculations would allow. This simplifies the process and ensures the solution's robustness, making it an invaluable resource for students tackling advanced mathematical problems.