Christoph Cremer

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

$\mathbf {x'} (t)=\mathbf {Ax} (t)$

where $\mathbf {x} (t)$ is an n×1 vector of functions of an underlying variable Template:Mvar, $\mathbf {x'} (t)$ is the vector of first derivatives of these functions, and $\mathbf {A}$ is an n×n matrix, of which all elements are constants.

Note that by using the Cayley-Hamilton theorem and Vandermonde-type matrices, a solution may be given in a simple form. Below the solution is displayed in terms of Putzer's algorithm.

This differential equation has the following general solution:

$\mathbf {x} (t)=c_{1}e^{\lambda _{1}t}\mathbf {u} _{1}+c_{2}e^{\lambda _{2}t}\mathbf {u} _{2}+\cdots +c_{n}e^{\lambda _{n}t}\mathbf {u} _{n}$

where $\lambda _{1}$, $\lambda _{2}$, ... and $\lambda _{n}$ are the eigenvalues of $\mathbf {A}$; $\mathbf {u} _{1}$, $\mathbf {u} _{2}$, ... and $\mathbf {u} _{n}$ are the respective eigenvectors of $\mathbf {A}$ and $c_{1}$, $c_{2}$, .... and $c_{n}$ are constants.

Stability and steady state of the matrix system

The matrix equation $\mathbf {x} '(t)=\mathbf {Ax} (t)+\mathbf {b}$ with n×1 parameter vector b is stable if and only if all eigenvalues of the matrix A have a negative real part. The steady state x* to which it converges if stable is found by setting $\mathbf {x} '(t)=\mathbf {0}$, yielding $\mathbf {x} ^{*}=-\mathbf {A} ^{-1}\mathbf {b}$, assuming A is invertible. Thus the original equation can be written in homogeneous form in terms of deviations from the steady state: $\mathbf {x} '(t)=\mathbf {A} [\mathbf {x} (t)-\mathbf {x} ^{*}]$. A different way of expressing this (closer to regular usage) is that x* is a particular solution to the in-homogenous equation, and all solutions are in the form $\mathbf {x} _{h}+\mathbf {x} ^{*}$, with $\mathbf {x} _{h}$ a solution to the homogenous equation (b=0).

Solution in matrix form

Matrix exponentials can be used to express the solution of $\mathbf {x} '(t)=\mathbf {A} [\mathbf {x} (t)-\mathbf {x} ^{*}]$

$\mathbf {x} (t)=\mathbf {x} ^{*}+e^{\mathbf {A} t}[\mathbf {x} (0)-\mathbf {x} ^{*}]$.

Putzer Algorithm for computing Buying, selling and renting HDB and personal residential properties in Singapore are simple and transparent transactions. Although you are not required to engage a real property salesperson (generally often known as a "public listed property developers In singapore agent") to complete these property transactions, chances are you'll think about partaking one if you are not accustomed to the processes concerned.

Professional agents are readily available once you need to discover an condominium for hire in singapore In some cases, landlords will take into account you more favourably in case your agent comes to them than for those who tried to method them by yourself. You need to be careful, nevertheless, as you resolve in your agent. Ensure that the agent you are contemplating working with is registered with the IEA – Institute of Estate Brokers. Whereas it might sound a hassle to you, will probably be worth it in the end. The IEA works by an ordinary algorithm and regulations, so you'll protect yourself in opposition to probably going with a rogue agent who prices you more than they should for his or her service in finding you an residence for lease in singapore.

There isn't any deal too small. Property agents who are keen to find time for any deal even if the commission is small are the ones you want on your aspect. Additionally they present humbleness and might relate with the typical Singaporean higher. Relentlessly pursuing any deal, calling prospects even without being prompted. Even if they get rejected a hundred times, they still come again for more. These are the property brokers who will find consumers what they need eventually, and who would be the most successful in what they do. 4. Honesty and Integrity

This feature is suitable for you who need to get the tax deductions out of your PIC scheme to your property agency firm. It's endorsed that you visit the correct site for filling this tax return software. This utility must be submitted at the very least yearly to report your whole tax and tax return that you're going to receive in the current accounting 12 months. There may be an official website for this tax filling procedure. Filling this tax return software shouldn't be a tough thing to do for all business homeowners in Singapore.

A wholly owned subsidiary of SLP Worldwide, SLP Realty houses 900 associates to service SLP's fast rising portfolio of residential tasks. Real estate is a human-centric trade. Apart from offering comprehensive coaching applications for our associates, SLP Realty puts equal emphasis on creating human capabilities and creating sturdy teamwork throughout all ranges of our organisational hierarchy. Worldwide Presence At SLP International, our staff of execs is pushed to make sure our shoppers meet their enterprise and investment targets. Under is an inventory of some notable shoppers from completely different industries and markets, who've entrusted their real estate must the expertise of SLP Worldwide.

If you're looking for a real estate or Singapore property agent online, you merely need to belief your instinct. It is because you don't know which agent is sweet and which agent will not be. Carry out research on a number of brokers by looking out the internet. As soon as if you find yourself certain that a selected agent is dependable and trustworthy, you'll be able to choose to utilize his partnerise find you a house in Singapore. More often than not, a property agent is considered to be good if she or he places the contact data on his web site. This is able to imply that the agent does not thoughts you calling them and asking them any questions regarding properties in Singapore. After chatting with them you too can see them of their office after taking an appointment.

Another method by way of which you could find out whether the agent is sweet is by checking the feedback, of the shoppers, on the website. There are various individuals would publish their comments on the web site of the Singapore property agent. You can take a look at these feedback and the see whether it will be clever to hire that specific Singapore property agent. You may even get in contact with the developer immediately. Many Singapore property brokers know the developers and you may confirm the goodwill of the agent by asking the developer.

Given a matrix $\mathbf {A}$ with eigenvalues $\lambda _{1},\lambda _{2},\dots ,\lambda _{n}$ then

$e^{\mathbf {A} t}=\sum _{j=0}^{n-1}r_{j+1}{\left(t\right)}\mathbf {P} _{j}$

Where

$\mathbf {P} _{0}=\mathbf {I}$

$\mathbf {P} _{j}=\prod _{k=1}^{j}\left(\mathbf {A} -\lambda _{k}\mathbf {I} \right)=\mathbf {P} _{j-1}\left(\mathbf {A} -\lambda _{j}\mathbf {I} \right),j=1,2,\dots ,n-1$

${\dot {r}}_{1}=\lambda _{1}r_{1}$

$r_{1}{\left(0\right)}=1$

${\dot {r}}_{j}=\lambda _{j}r_{j}+r_{j-1},j=2,3,\dots ,n$

$r_{j}{\left(0\right)}=0,j=2,3,\dots ,n$

The equations for $r_{i}{\left(t\right)}$ are simple first order nonhomogeneous ODEs.

Notice the algorithm does not require that the matrix $\mathbf {A}$ is diagonalizable and avoids the complexity of using Jordan canonical form when it is not needed.

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:

${\frac {dx}{dt}}=a_{1}x+b_{1}y,\quad {\frac {dy}{dt}}=a_{2}x+b_{2}y$

where $a_{1},a_{2},b_{1}\,\!$ and $b_{2}\,\!$ 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:

• Finding the eigenvalues
• Finding the eigenvectors
• Finding the needed functions

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 the matrix ODE's according to the three steps above, using simple matrices in the process, let us find, say, function $x\,\!$ and function $y\,\!$, both in terms of the single underlying variable t, in the following linear differential equation of the first order:

${\frac {dx}{dt}}=3x-4y,\quad {\frac {dy}{dt}}=4x-7y.$

To solve this particular ordinary differential equation, at some point of the solving process, we need an initial value, a starting point. In this case, we use $x(0)=y(0)=1.\,\!$

First step

The first step, that has already been mentioned above, is finding the eigenvalues. The process of finding the eigenvalues is not a very difficult process. Both eigenvalues and eigenvectors are useful in numerous branches of mathematics, including higher engineering mathematics/calculations(i.e. Applied Mathematics), mechanics, physical mathematics, mathematical economics, and linear algebra.

Therefore, the process consists of the following:

${\begin{pmatrix}x'\\y'\end{pmatrix}}$ = ${\begin{pmatrix}3&-4\\4&-7\end{pmatrix}}{\begin{pmatrix}x\\y\end{pmatrix}}.$

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 shown above, we may start the process of evaluating the eigenvalues. To do that we are going to have to find the determinant of the matrix that is formed when an identity matrix, $I_{n}\,\!$, multiplied by some constant lambda, symbol λ, is subtracted from our coefficient matrix in the following way:

$\det \left({\begin{bmatrix}3&-4\\4&-7\end{bmatrix}}-\lambda {\begin{bmatrix}1&0\\0&1\end{bmatrix}}\right)$.

Applying further simplification and basic rules of matrix addition we come up with the following:

$\det {\begin{bmatrix}3-\lambda &-4\\4&-7-\lambda \end{bmatrix}}$.

Applying the rules of finding the determinant of a single 2×2 matrix, we obtain the following elementary quadratic equation:

$\det {\begin{bmatrix}3-\lambda &-4\\4&-7-\lambda \end{bmatrix}}=0$

$-21-3\lambda +7\lambda +\lambda ^{2}+16=0\,\!$

which may be reduced further to get a simpler version of the above:

$\lambda ^{2}+4\lambda -5=0\,\!$.

Now finding the two roots, $\lambda _{1}\,\!$ and $\lambda _{2}\,\!$ of the given quadratic equation by applying the factorization method we get the following:

$\lambda ^{2}+5\lambda -\lambda -5=0\,\!$

$\lambda (\lambda +5)-1(\lambda +5)=0\,\!$

$(\lambda -1)(\lambda +5)=0\,\!$

$\lambda =1,-5\,\!$.

The values, $\lambda _{1}=1\,\!$ and $\lambda _{2}=-5\,\!$, that we have calculated above are the required eigenvalues. Once we find these two values, we proceed to the second step of the solution. We'll use the calculated eigenvalues later in the final solution. 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, dramatically change.

Second step

As it was already mentioned above, in a simple description, this step involves finding the eigenvectors by means of using the information originally given to us.

For each of the eigenvalues calculated we are going to have an individual eigenvector. For our first eigenvalue, which is $\lambda _{1}=1\,\!$, we have the following:

${\begin{pmatrix}3&-4\\4&-7\end{pmatrix}}{\begin{pmatrix}\alpha \\\beta \end{pmatrix}}=1{\begin{pmatrix}\alpha \\\beta \end{pmatrix}}.$

Simplifying the above expression by applying basic matrix multiplication rules we have:

$3\alpha -4\beta =\alpha$

$\alpha =2\beta .$.

All of these calculations have been done only to obtain the last expression, which in our case is $\alpha =2\beta \,\!$. Now taking some arbitrary value, presumably a small insignificant value, which is much easier to work with, for either $\alpha \,\!$ or $\beta \,\!$ (in most cases it does not really matter), we substitute it into $\alpha =2\beta \,\!$. Doing so produces a very simple vector, which is the required eigenvector for this particular eigenvalue. In our case, we pick $\alpha =2\,\!$, which, in turn determines that $\beta =1\,\!$ and, using the standard vector notation, our vector looks like this:

$\mathbf {\hat {v}} _{1}={\begin{pmatrix}2\\1\end{pmatrix}}.$

Performing the same operation using the second eigenvalue we calculated, which is $\,\!\,\lambda =-5\,,$, we obtain our second eigenvector. The process of working out this vector is not shown, but the final result is as follows:

$\mathbf {\hat {v}} _{2}={\begin{pmatrix}1\\2\end{pmatrix}}.$

Once we've found both needed vectors, we start the third and last step. Don't forget that we'll substitute the eigenvalues and eigenvectors determined above into a specialized equation (shown shortly).

Third (final) step

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:

${\begin{pmatrix}x\\y\end{pmatrix}}=Ae^{\lambda _{1}t}\mathbf {\hat {v}} _{1}+Be^{\lambda _{2}t}\mathbf {\hat {v}} _{2}.$

Substituting the values of eigenvalues and eigenvectors we get the following expression:

${\begin{pmatrix}x\\y\end{pmatrix}}=Ae^{t}{\begin{pmatrix}2\\1\end{pmatrix}}+Be^{-5t}{\begin{pmatrix}1\\2\end{pmatrix}}.$

Applying further simplification rules we have:

${\begin{pmatrix}2&1\\1&2\end{pmatrix}}{\begin{pmatrix}Ae^{t}\\Be^{-5t}\end{pmatrix}}={\begin{pmatrix}x\\y\end{pmatrix}}.$

Simplifying further and writing the equations for functions $x\,\!$ and $y\,\!$ separately:

$x=2Ae^{t}+Be^{-5t}\,\!$

$y=Ae^{t}+2Be^{-5t}.\,\!$

The above equations are in fact the functions that we needed to find, but they are in their general form and if we want to actually find their exact forms and solutions, now is the time to look back at the information given to us, the so-called initial value problem. At some point during solving these equations we have come across $x(0)=y(0)=1\,\!$, which plays the role of starting point for our ordinary differential equation. Now is the time to apply this condition, which lets us find the constants, A and B. As we see from the $x(0)=y(0)=1\,\!$ condition, when $t=0\,\!$, the overall equation is equal to 1. Thus we may construct the following system of linear equations:

$1=2A+B\,\!$

$1=A+2B.\,\!$

Solving these equations we find that both constants A and B are equal to 1/3. Therefore if we substitute these values into the general form of these two functions we have their exact forms:

$x={\frac {2}{3}}e^{t}+{\frac {1}{3}}e^{-5t}\,\!$

$y={\frac {1}{3}}e^{t}+{\frac {2}{3}}e^{-5t}.\,\!$

which is our final form of the two functions we were required to find.

See also

• Nonhomogeneous equations
• Matrix difference equation
• Newton's law of cooling
• Fibonacci Sequence
• Difference equations
• Wave equation

References

43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.