MHB System of differential equations involving matrices

Guest2
Messages
192
Reaction score
0
Let $\mathbf{A} = \begin{pmatrix} 9 & 2 \\ 1 & 8 \end{pmatrix}.$ Obtain the general solution $\mathbf{y}(t)$ of the system of differential equations $\displaystyle \frac{d\mathbf{y}}{dt} = \mathbf{Ay}$:

$\begin{cases} \frac{dy_1}{dt} = 9y_1+2y_2 \\ \frac{dy_2}{dt} = y_1+8y_2 \end{cases}$

and find the unique solution satisfying $y_1(0) = 1$ and $y_2(0) = 5$.​

I've found that $\mathbf{y} = ke^{\mathbf{A}t}$, where $k \in \mathbb{R}$. That also eigenvalues of $\mathbf{A}$ are $\lambda = 7,10$. But I do not know how to proceed.
 
Last edited:
Physics news on Phys.org
Re: System of ifferential equations involving matrices

Guest said:
Let $\mathbf{A} = \begin{pmatrix} 9 & 2 \\ 1 & 8 \end{pmatrix}.$ Obtain the general solution $\mathbf{y}(t)$ of the system of differential equations $\displaystyle \frac{d\mathbf{y}}{dt} = \mathbf{Ay}$:

$\begin{cases} \frac{dy_1}{dt} = 9y_1+2y_2 \\ \frac{dy_2}{dt} = y_1+8y_2 \end{cases}$

and find the unique solution satisfying $y_1(0) = 1$ and $y_2(0) = 5$.​

I've found that $\mathbf{y} = ke^{\mathbf{A}t}$, where $k \in \mathbb{R}$. That also eigenvalues of $\mathbf{A}$ are $\lambda = 7,10$. But I do not know how to proceed.
Well, you have non-zero eigenvalues for A. That means you can diagonalize it. What happens to your two equations if you change the basis so that A is diagonal?

-Dan
 
Re: System of ifferential equations involving matrices

topsquark said:
Well, you have non-zero eigenvalues for A. That means you can diagonalize it. What happens to your two equations if you change the basis so that A is diagonal?

-Dan
I haven't obtained the two equations yet. They're asking me to show it. I'm at the step where I've $\displaystyle \frac{d\mathbf{y}}{dt} = \mathbf{Ay} \implies \mathbf{y} = ke^{\mathbf{A}t}$.
 
Guest said:
Let $\mathbf{A} = \begin{pmatrix} 9 & 2 \\ 1 & 8 \end{pmatrix}.$ Obtain the general solution $\mathbf{y}(t)$ of the system of differential equations $\displaystyle \frac{d\mathbf{y}}{dt} = \mathbf{Ay}$:

$\begin{cases} \frac{dy_1}{dt} = 9y_1+2y_2 \\ \frac{dy_2}{dt} = y_1+8y_2 \end{cases}$

and find the unique solution satisfying $y_1(0) = 1$ and $y_2(0) = 5$.​
Ignore the y = ke^{At} thing for a moment.

This is the long way around for this level of problem, but it generalizes nicely.

What I want you to do is diagonalize A. You have non-zero eigenvalues, so we can do the following. Find a transformation matrix S to act on you original equations as so:
[math]S \vec{y}' = SA \vec{y} = (S A S^{-1}) (S \vec{y} )[/math]

such that
[math]S A S^{-1} = \left ( \begin{matrix} 7 & 0 \\ 0 & 10 \end{matrix} \right )[/math]

One possibility is [math]S = \left ( \begin{matrix} 1 & -2 \\ 1 & 1 \end{matrix} \right )[/math]

Thus
[math]S \vec{y}' = \left ( \begin{matrix} 7 & 0 \\ 0 & 10 \end{matrix} \right ) S(\vec{y})[/math]

Call [math]S \vec{y} = Y[/math] and thus we also have [math]S \vec{y}' = Y'[/math]

So you need to solve
[math]\left ( \begin{matrix} Y_1' \\ Y_2 ' \end{matrix} \right ) = \left ( \begin{matrix} 7 & 0 \\ 0 & 10 \end{matrix} \right ) \left ( \begin{matrix} Y_1 \\ Y_2 \end{matrix} \right )[/math]

This decouples nicely. Once you have Y_1 and Y_2 then you can solve the simultaneous equations problem to find y_1 and y_2.

-Dan
 
topsquark said:
...
I was actually trying to get the system of equations in the cases, i.e. the one I was supposed solve (I completely misread the problem)! (Rofl)

Thank you. I'll study your post.
 
There is the following linear Volterra equation of the second kind $$ y(x)+\int_{0}^{x} K(x-s) y(s)\,{\rm d}s = 1 $$ with kernel $$ K(x-s) = 1 - 4 \sum_{n=1}^{\infty} \dfrac{1}{\lambda_n^2} e^{-\beta \lambda_n^2 (x-s)} $$ where $y(0)=1$, $\beta>0$ and $\lambda_n$ is the $n$-th positive root of the equation $J_0(x)=0$ (here $n$ is a natural number that numbers these positive roots in the order of increasing their values), $J_0(x)$ is the Bessel function of the first kind of zero order. I...
Are there any good visualization tutorials, written or video, that show graphically how separation of variables works? I particularly have the time-independent Schrodinger Equation in mind. There are hundreds of demonstrations out there which essentially distill to copies of one another. However I am trying to visualize in my mind how this process looks graphically - for example plotting t on one axis and x on the other for f(x,t). I have seen other good visual representations of...
Back
Top