Coupled differential equations

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
8 replies · 2K views
aaaa202
Messages
1,144
Reaction score
2
A bit related to my other topic but here goes:

I have the set of differential equations:

x1' = a-ax1-ax2
x2' = b-bx1-bx2

What is the solution to such a set? I could google systems of differential equations, but that turns up large texts on the theory. I just wonna know in a few lines what the solution to such a set is. I know how to diagonalize the coefficient matrix etc. if I need that.
 
Physics news on Phys.org
aaaa202 said:
A bit related to my other topic but here goes:

I have the set of differential equations:

x1' = a-ax1-ax2
x2' = b-bx1-bx2

What is the solution to such a set? I could google systems of differential equations, but that turns up large texts on the theory. I just wonna know in a few lines what the solution to such a set is. I know how to diagonalize the coefficient matrix etc. if I need that.

I did a search for "nonhomogeneous system of linear equations". Here is one link.
http://tutorial.math.lamar.edu/classes/de/NonhomogeneousSystems.aspx
 
aaaa202 said:
A bit related to my other topic but here goes:

I have the set of differential equations:

x1' = a-ax1-ax2
x2' = b-bx1-bx2

What is the solution to such a set? I could google systems of differential equations, but that turns up large texts on the theory. I just wonna know in a few lines what the solution to such a set is. I know how to diagonalize the coefficient matrix etc. if I need that.

Are the three 'a's in the first equation really supposed to be all the same (and ditto for the 'b's in the second equation)? If so, try looking at y = x1+x2.

BTW: what does the word "wonna" mean?
 
'wonna' is a variant of 'wanna' = 'want to'
 
SteamKing said:
'wonna' is a variant of 'wanna' = 'want to'

I know that of course; but it does no harm to try to encourage proper language.
 
There several different ways to handle
[itex]x_1'= a- ax_1- ax_2[/itex]
[itex]x_2'= b- bx_1- bx_2[/itex]

The more "elementary" method would be to differentiate the first equation again:
[itex]x_1''= -ax_1'- ax_2'= ax_1'- a(b- bx_1- bx_2)= (a+ ab)x_1+ abx_2- ab[itex] but, from the first equation, multiplied by b, [itex]abx_2= ab- bx_1'- abx_1[itex]so <br /> [itex]x_1''= (a+ ab)x_1+ (ab- bx_1'- abx_1)- ab[/itex]<br /> [itex]x_1''+ bx_1'- ax_1= 0[/itex].<br /> That has characteristic equation [itex]r^2+ br- a= 0[/itex] with roots [itex]r= \frac{-b\pm\sqrt{b^2+ 4a}}{2}[/itex] so that [itex]x_1(t)= e^{-bt}(C_1e^{\sqrt{b^2+ 4a}t}+ C_2e^{-\sqrt{b^2+ 4a}t}[/itex] if [itex]b^2+ 4a>0[/itex], [itex]x_1(t)= e^{-bt}(C_1 cos(\sqrt{-b^2- 4a}t)+ C_2 sin(\sqrt{-b^2- 4a}t)[/itex] if [itex]b^2+ 4a< 0[/itex], or [itex]x_1(t)= e^{-bt}(C_1+ C_2t)[/itex] if [itex]b^2+ 4a= 0[/itex].<br /> <br /> In any case, [itex]x_2= (ab- bx_1'- abx_1)/ab[/itex].<br /> <br /> A more "sophisticated" method would be to write the differential equations as a matrix equation:<br /> [tex]\begin{pmatrix}x_1 \\ x_2\end{pmatrix}'= \begin{pmatrix}-a & -a \\ -b & -b\end{pmatrix}\begin{pmatrix}x_1 \\ x_2 \end{pmatrix}+ \begin{pmatrix}a \\ b\end{pmatrix}[/tex]<br /> <br /> The characteristic equation for that matrix is <br /> [tex]\left|\begin{array}{cc}-a- \lambda & -a \\ -b & -b- \lambda\end{array}\right|= (-a- \lambda)(-b- \lambda)- ab= \lambda^2+ (a+ b)\lambda= 0[/tex]<br /> which has roots [itex]\lambda= 0[/itex] and [itex]\lambda= -a- b[itex].<br /> We can find the eigenvectors corresponding to those eigenvectors and so diagonalize the matrix,[/itex][/itex][/itex][/itex][/itex][/itex]
 
HallsofIvy said:
There several different ways to handle
[itex]x_1'= a- ax_1- ax_2[/itex]
[itex]x_2'= b- bx_1- bx_2[/itex]

The more "elementary" method would be to differentiate the first equation again:
[itex]x_1''= -ax_1'- ax_2'= ax_1'- a(b- bx_1- bx_2)= (a+ ab)x_1+ abx_2- ab[itex] but, from the first equation, multiplied by b, [itex]abx_2= ab- bx_1'- abx_1[itex]so <br /> [itex]x_1''= (a+ ab)x_1+ (ab- bx_1'- abx_1)- ab[/itex]<br /> [itex]x_1''+ bx_1'- ax_1= 0[/itex].<br /> That has characteristic equation [itex]r^2+ br- a= 0[/itex] with roots [itex]r= \frac{-b\pm\sqrt{b^2+ 4a}}{2}[/itex] so that [itex]x_1(t)= e^{-bt}(C_1e^{\sqrt{b^2+ 4a}t}+ C_2e^{-\sqrt{b^2+ 4a}t}[/itex] if [itex]b^2+ 4a>0[/itex], [itex]x_1(t)= e^{-bt}(C_1 cos(\sqrt{-b^2- 4a}t)+ C_2 sin(\sqrt{-b^2- 4a}t)[/itex] if [itex]b^2+ 4a< 0[/itex], or [itex]x_1(t)= e^{-bt}(C_1+ C_2t)[/itex] if [itex]b^2+ 4a= 0[/itex].<br /> <br /> In any case, [itex]x_2= (ab- bx_1'- abx_1)/ab[/itex].<br /> <br /> A more "sophisticated" method would be to write the differential equations as a matrix equation:<br /> [tex]\begin{pmatrix}x_1 \\ x_2\end{pmatrix}'= \begin{pmatrix}-a & -a \\ -b & -b\end{pmatrix}\begin{pmatrix}x_1 \\ x_2 \end{pmatrix}+ \begin{pmatrix}a \\ b\end{pmatrix}[/tex]<br /> <br /> The characteristic equation for that matrix is <br /> [tex]\left|\begin{array}{cc}-a- \lambda & -a \\ -b & -b- \lambda\end{array}\right|= (-a- \lambda)(-b- \lambda)- ab= \lambda^2+ (a+ b)\lambda= 0[/tex]<br /> which has roots [itex]\lambda= 0[/itex] and [itex]\lambda= -a- b[itex].<br /> We can find the eigenvectors corresponding to those eigenvectors and so diagonalize the matrix,[/itex][/itex][/itex][/itex][/itex][/itex]
[itex][itex][itex][itex][itex][itex] <br /> The two DEs added together give ##y' = c - cy## where ##c = a+b## and ##y = x_1 + x_2##. Once we know ##y## we can get ##x_1, x_2## from ##x_1' = a - ay##, etc.[/itex][/itex][/itex][/itex][/itex][/itex]