Coupled differential equations

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
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.
 
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?
 
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]