Matrix differential equation
A differential equation is a mathematical equation for an unknown function of one or several variables that relates the values of the function itself and of its derivatives of various orders. A matrix differential equation contains more than one function stacked into vector form with a matrix relating the functions to their derivatives.
For example, a simple matrix ordinary differential equation is
where x(t) is an n×1 vector of functions of an underlying variable t, x'(t) is the vector of first derivatives of these functions, and A is an matrix, of which all elements are constants.
In the case where A has n distinct eigenvalues, this differential equation has the following general solution,
By use of the Cayley-Hamilton theorem and Vandermonde-type matrices, this formal matrix exponential solution may be reduced to a simple form. Below, this solution is displayed in terms of Putzer's algorithm.
- 1 Stability and steady state of the matrix system
- 2 Solution in matrix form
- 3 Putzer Algorithm for computing eAt
- 4 Deconstructed example of a matrix ordinary differential equation
- 5 Solving deconstructed matrix ordinary differential equations
- 6 Solved example of a matrix ODE
- 7 See also
- 8 References
Stability and steady state of the matrix system
The matrix equation
assuming A is invertible.
Thus, the original equation can be written in homogeneous form in terms of deviations from the steady state,
An equivalent way of expressing this is that x* is a particular solution to the non-homogeneous equation, while all solutions are in the form
with a solution to the homogeneous equation (b=0).
Stability of the two-state-variable case
In the n = 2 case (with two state variables), the stability conditions that the two eigenvalues of the transition matrix A each have a negative real part are equivalent to the conditions that the trace of A be negative and its determinant be positive.
Solution in matrix form
The formal solution of is the celebrated matrix exponential,
evaluated in a multitude of techniques.
Putzer Algorithm for computing eAt
Given a matrix A with eigenvalues ,
The equations for are simple first order inhomogeneous ODEs.
Deconstructed example of a matrix ordinary differential equation
A first-order homogeneous matrix ordinary differential equation in two functions x(t) and y(t), when taken out of matrix form, has the following form:
where and may be any arbitrary scalars.
Higher order matrix ODE's may possess a much more complicated form.
Solving deconstructed matrix ordinary differential equations
The process of solving the above equations and finding the required functions, of this particular order and form, consists of 3 main steps. Brief descriptions of each of these steps are listed below:
The final, third, step in solving these sorts of ordinary differential equations is usually done by means of plugging in the values, calculated in the two previous steps into a specialized general form equation, mentioned later in this article.
Solved example of a matrix ODE
To solve a matrix ODE according to the three steps detailed above, using simple matrices in the process, let us find, say, a function x and a function y both in terms of the single independent variable t, in the following linear differential equation of the first order,
To solve this particular ordinary differential equation system, at some point of the solution process we shall need a set of two initial values (corresponding to the two state variables at the starting point). In this case, let us pick x(0)=y(0)=1.
The first step, already mentioned above, is finding the eigenvalues of A in
The derivative notation x' etc. seen in one of the vectors above is known as Lagrange's notation,(first introduced by Joseph Louis Lagrange. It is equivalent to the derivative notation dx/dt used in the previous equation, known as Leibniz's notation, honouring the name of Gottfried Leibniz.)
Once the coefficients of the two variables have been written in the matrix form A displayed above, one may evaluate the eigenvalues. To that end, one finds the determinant of the matrix that is formed when an identity matrix, , multiplied by some constant λ, is subtracted from the above coefficient matrix to yield the characteristic polynomial of it,
and solve for its zeroes.
Applying further simplification and basic rules of matrix addition yields
Applying the rules of finding the determinant of a single 2×2 matrix,yields the following elementary quadratic equation,
which may be reduced further to get a simpler version of the above,
The values and , calculated above are the required eigenvalues of A. In some cases, say other matrix ODE's, the eigenvalues may be complex, in which case the following step of the solving process, as well as the final form and the solution, may dramatically change.
As mentioned above, this step involves finding the eigenvectors of A from the information originally provided.
Simplifying the above expression by applying basic matrix multiplication rules yields
All of these calculations have been done only to obtain the last expression, which in our case is α=2β. Now taking some arbitrary value, presumably a small insignificant value, which is much easier to work with, for either α or β (in most cases it does not really matter), we substitute it into α=2β. Doing so produces a simple vector, which is the required eigenvector for this particular eigenvalue. In our case, we pick α=2, which, in turn determines that β=1 and, using the standard vector notation, our vector looks like
This final step actually finds the required functions that are 'hidden' behind the derivatives given to us originally. There are two functions, because our differential equations deal with two variables.
The equation, which involves all the pieces of information that we have previously found has the following form,
Applying further simplification,
Simplifying further and writing the equations for functions and separately,
The above equations are, in fact, the general functions sought, but they are in their general form; whilst we want to actually find their exact forms and solutions, so now is the time to look back at the information given to us, the so-called initial value problem. We chose , which plays the role of starting point for our ordinary differential equation; application of this condition specifies the constants, A and B. As we see from the condition, when t=0, the overall equation is equal to 1. Thus we may construct the following system of linear equations,
Solving these equations, we find that both constants A and B equal 1/3. Therefore substituting these values into the general form of these two functions specifies their exact forms,
the two functions sought.
- Nonhomogeneous equations
- Matrix difference equation
- Newton's law of cooling
- Fibonacci Sequence
- Difference equations
- Wave equation
- H. Moya-Cessa, F. Soto-Eguibar, DIFFERENTIAL EQUATIONS: AN OPERATIONAL APPROACH, RINTON PRESS, New Jersey, 2011. ISBN 978-1-58949-060-4
- E. J. Putzer (1966). "Avoiding the Jordan Canonical Form in the Discussion of Linear Systems with Constant Coefficients", The American Mathematical Monthly, 73, No. 1 (1966) 2-7. , Paper on JSTOR