A linear system of differential equations is an ODE (ordinary differential equation) of the type: $$$x'(t)=A(t)\cdot x+ b(t)$$$ Where, $$A(t)$$ is a matrix, $$n \times n$$, of functions of the variable $$t$$, $$b (t)$$ is a dimension $$n$$ vector of functions of the variable $$t$$, and $$x$$ is a vector of dimension $$n$$ that is the function that we want to find.
An example of system of linear ODE would be: $$$\begin{pmatrix}x \\ y \end{pmatrix}'= \begin{pmatrix} 0 & \ln t \\ -e^t & 3 \cos t \end{pmatrix} \cdot \begin{pmatrix} x \\ y \end{pmatrix}+\begin{pmatrix} e^t \\ 3e^t \end{pmatrix}$$$
A linear system of dimension $$n$$ has $$n$$ linearly independent solutions and solving the system means finding all of them. Any solution is a linear combination of these $$n$$ solutions. Therefore, given a vector of initial conditions ($$n$$), we will determine the $$n$$ constants and will find the exact solution.
When we solve a linear system we will place the $$n$$ vectors that are the solution (linearly independent) in the columns of a matrix, the so-called fundamental matrix of the system ($$n \times n$$). Therefore, solving the system is to find the fundamental matrix. By multiplying this matrix by a vector of arbitrary constants we will have the general solution.
An important property of the fundamental matrices is that if we multiply a fundamental matrix by a constant matrix with a determinant different from zero, the result is another fundamental matrix (it is important that the constant matrix multiplies from the right, if not this is not true).
There are no explicit methods to solve these types of equations, (only in dimension 1). Nevertheless, there are some particular cases that we will be able to solve: Homogeneous systems of ode's with constant coefficients, Non homogeneous systems of linear ode's with constant coefficients, and Triangular systems of differential equations.
In this unit we are going to explain the Triangular systems of differential equations.
We know how to solve linear systems with constant coefficients, and there are no methods for solving systems where the matrix $$A$$ is a matrix of functions unless we have a very special matrix, for example a triangular matrix.
Let's suppose that we have the following system: $$$\begin{pmatrix} x' \\ y' \\ z' \end{pmatrix} = \begin{pmatrix} a_{11}(t) & 0 & 0 \\ a_{21}(t) & a_{22}(t) & 0 \\ a_{31}(t) & a_{32}(t) & a{33}(t) \end{pmatrix} \cdot \begin{pmatrix} x \\ y \\ z \end{pmatrix} + \begin{pmatrix} b_1(t) \\ b_2(t) \\ b_3(t) \end{pmatrix}$$$ i.e., a non homogeneous triangular linear system.
The idea is to solve the system step by step. Notice the first equation:$$$x'(t)=a_{11}(t) \cdot x(t) +b_1(t)$$$
This is a linear equation like the ones we solved in the first unit of differential equations.
Therefore we can calculate its solution.
Now we concentrate on the second equation: $$$y'(t)=a_{21}(t)\cdot x(t)+a_22(t)\cdot y(t)+b_2(t)$$$ Now, $$x (t)$$, is already known.
Replacing the function obtained in the previous step we have it $$$y'(t)=a_{22}(t) \cdot y(t) +\widetilde{b}_2(t), \mbox{ where } \widetilde{b}_2(t) = a_{12}(t) \cdot x(t)+b_2(t)$$$ is a well-known function.
Therefore we obtain another linear equation that we can already solve.
Finally, we take the third equation and replace the values of $$x (t)$$ and $$y(t)$$, obtaining: $$$z'(t)=a_{33}(t) \cdot z(t) + \widetilde{b}_3 \mbox{ where } \widetilde{b}_3(t)=a_{31}(t) \cdot x(t)+a_{32} (t) \cdot y(t)+b_3(t) $$$ which is a well-known function. This is another linear equation that we can solve.
Let's notice that we have considered a low triangular counterfoil, but we could use the same procedure if we had an upper triangular matrix but starting from below.
Thus we already solved the system.
Let's see it more clearly with an example.
Let's consider the following system: $$$\displaystyle \begin{pmatrix} x' \\ y' \\ z' \end{pmatrix} = \begin{pmatrix} \frac{1}{t} & 0 & 0\\ t & \frac{2}{t} & 0 \\ -t^3 & t & \frac{3}{t} \end{pmatrix} \cdot \begin{pmatrix} x \\ y \\ z \end{pmatrix} + \begin{pmatrix} 2t^2 \\ -t^4 \\ t^6 \end{pmatrix}$$$ We are going to proceed as we described:
-
We take the first equation:$$$\displaystyle x'=\frac{1}{t} \cdot x +2 t^2$$$ Since this is a linear ODE, we solve it following the procedure described in other units:
-
Resolution of the homogeneous part: $$$\displaystyle x_h(t)=k_1\cdot e^{\displaystyle \int \frac{1}{t} \ dt }=k_1\cdot e^{\ln t}=k_1 \cdot t$$$
-
Search for a particular solution:$$x_p(t)=t\cdot u(t)$$, where $$u(t)$$ satisfies: $$\displaystyle u'=\frac{2t^2}{t}=2t$$. Therefore,$$u(t)=t^2$$ and thus $$x_p(t)=t^3$$
- General solution $$x(t)=k_1 \cdot t+t^3$$
-
-
We take the second equation, replacing the obtained value of x: $$\displaystyle y'=t \cdot x+ \displaystyle\frac{2}{t}\cdot y- t^4=t \cdot\Big(k_1\cdot t+t^3\Big)+\displaystyle\frac{2}{t}\cdot y-t^4=\displaystyle\frac{2}{t}\cdot y +k_1\cdot t^2$$
-
Resolution of the homogeneous part: $$y_h(t)=k_2\cdot e^{\displaystyle \int \frac{2}{t} \ dt }= k_2 \cdot e^{ \ln t^2}=k_2 \cdot t^2$$
-
Search for a particular solution:$$y_p(t)= t^2 \cdot u(t)$$, where $$u(t)$$ satisfies: $$u'= \displaystyle \frac{k_1 \cdot t^2}{t^2}=k_1$$. Therefore, $$u(t)=k_1\cdot t$$ and thus $$y_p(t)=k_1\cdot t^3$$
- General solution $$y(t)=k_2 \cdot t^2+k_1\cdot t^3$$
-
-
We take the third equation, replacing the obtained values of $$x$$ and $$y$$: $$$z'=-t^3\cdot x+t \cdot y +\displaystyle \frac{3}{t}\cdot z+t^6=$$$ $$$=-t^3\cdot \Big(k_1\cdot t+ t^3 \Big) +t\cdot \Big(k_2 \cdot t^2+k_1 \cdot t^3 \Big)+\displaystyle \frac{3}{t} \cdot z+t^6=$$$ $$$=\displaystyle \frac{3}{t} \cdot z+k_2 \cdot t^3$$$
-
Resolution of the homogeneous part: $$z_h(t)=k_3\cdot e^{\displaystyle \frac{3}{t} \ dt}=k_3 \cdot e^{ \ln t^3}=k_ 3\cdot t^3$$
-
Search for a particular solution: $$z_p(t)=t^3\cdot u(t)$$, where $$u(t)$$ satisfies: $$u'(t)=\displaystyle \frac{k_2 \cdot t^3}{t^3}=k_2$$. Therefore, $$u(t)=k_2\cdot t$$ and thus $$z_p(t)=k_2\cdot t^4$$
- General solution $$z(t)=k_3\cdot t^3+k_2 \cdot t^4$$
-
Therefore the solution of the system is:$$$\begin{pmatrix} x(t) \\ y(t) \\ z(t) \end{pmatrix}=\begin{pmatrix} k_1\cdot t+t^3 \\ k_2 \cdot t^2+k_1 \cdot t^3 \\ k_3 \cdot t^3+k_2 \cdot t^4\end{pmatrix}$$$