I am trying to solve a system of coupled differential equations to plot streamlines using Matlab.
The equations are these:
\begin{align} \frac{\mathrm dx}{\mathrm dt} &= -3x -5y \\ \frac{\mathrm dy}{\mathrm dt} &= 5x + 3y \end{align}
What method do you suggest for solving this system? I'd greatly appreciate any insight or suggestion. No need to solve the system :) as long as you tell me what literature I can refer to.
Thanks!
3 Answers
Another approach:
Consider the following IVP problem:
\begin{align} \frac{\mathrm{d}x}{\mathrm{d}t} &= -3x - 5y \\ \frac{\mathrm{d}y}{\mathrm{d}t} &= 5x + 3y \end{align}
with $x(0)=x_0$ and $y(0)=y_0$.
Then, Laplace-transform both sides of both equations to get:
\begin{align} s X(s) - x_0 & = - 3 X(s) -5 Y(s) \\ s Y(s) - y_0 & = 5 X(s) + 3Y(s) , \end{align}
which is an algebraic system for $X(s) = \mathcal{L}_sx(t)$ and $Y(s) = \mathcal{L}_sy(t)$. Solve for the unknowns using (for example) Gauss elimination and compute the inverse Laplace transfrom to get the solution.
Cheers!
The equations may be rewritten $$ \left[\begin{array}{cc} x'(t)\\ y'(t)\end{array} \right]=\left[\begin{array}{cc} -3&-5\\ 5&3 \end{array} \right]\left[\begin{array}{cc} x(t) \\ y(t) \end{array} \right] $$ or $$ X'(t)=AX(t)$$
Do you see the analogy with ordinary differential equations?
I think a simple finite difference scheme will do the trick here.
Put a mesh over your domain. Then a point $i$ has a certain position $(x_i,\ y_i)$. Then for the first equation you can write:
$\frac{dx_i}{dt} = -3x_i-5y_i$.
Similar for the second equation. Solving for all $i$ and combining both gives you the streamlines.
