2 Solving systems of differential equations

The Laplace transform method is also well suited to solving systems of differential equations. A simple example will illustrate the technique.

Let x ( t ) , y ( t ) be two independent functions which satisfy the coupled differential equations

d x d t + y = e t d y d t x = 3 e t x ( 0 ) = 0 , y ( 0 ) = 1

Now, using a traditional approach, we could try to eliminate one of the unknown functions from this system: for example, from the first:

d y d t = e t d 2 x d t 2  (taking the derivative and rearranging)

This can then be substituted in the second equation: d y d t x = 3 e t , to give:

d 2 x d t 2 x = 4 e t

which can then be solved in the normal way (either using the complementary function/particular integral approach or else the Laplace transform approach.) However, this approach is not workable if we have large numbers of first order differential equations to deal with. Let us instead use the Laplace transform directly.

If we use the notation that

L { x ( t ) } = X ( s )  and L { y ( t ) } = Y ( s )

then, by taking the Laplace transform of every term in the given differential equations, we obtain:

x ( 0 ) + s X ( s ) + Y ( s ) = 1 s + 1 y ( 0 ) + s Y ( s ) X ( s ) = 3 s + 1

which, using the initial conditions and rearranging gives

s X ( s ) + Y ( s ) = 1 s + 1 X ( s ) + s Y ( s ) = s + 4 s + 1
Key Point 13

Taking the Laplace transform converts a system of differential equations

into a system of algebraic simultaneous equations.

We can solve these algebraic equations (in X ( s ) and Y ( s ) ) using a variety of techniques (inverse matrix; Cramer’s determinant method etc.) Here we will use Cramer’s method.

X ( s ) = 1 s + 1 1 s + 4 s + 1 s s 1 1 s = s s + 1 s + 4 s + 1 s 2 + 1 = 4 ( s 2 + 1 ) ( s + 1 ) = 2 ( s 1 ) s 2 + 1 2 s + 1

and

Y ( s ) = s 1 s + 1 1 s + 4 s + 1 s 1 1 s = s ( s + 4 ) s + 1 + 1 s + 1 s 2 + 1 = s 2 + 4 s + 1 ( s 2 + 1 ) ( s + 1 ) = 1 s + 1 + 2 ( s + 1 ) s 2 + 1

The last lines in each case having been obtained using partial fractions. We can now invert X ( s ) , Y ( s ) to find x ( t ) , y ( t ) :

x ( t ) = L 1 { X ( s ) } = 2 L 1 { s s 2 + 1 } 2 L 1 { 1 s 2 + 1 } 2 L 1 { 1 s + 1 } = ( 2 cos t 2 sin t 2 e t ) u ( t ) y ( t ) = L 1 { Y ( s ) } = L 1 { 1 s + 1 } + 2 L 1 { s s 2 + 1 } + 2 L 1 { 1 s 2 + 1 } = ( e t + 2 cos t + 2 sin t ) u ( t )

(Note that once the solution for x ( t ) is found the solution for y ( t ) may be easier to obtain by substituting in the differential equation: y = e t d x d t rather than using Laplace transforms.)

Task!

Use the Laplace transform to solve the coupled differential equations:

d y d t x = 0 , d x d t + y = 1 , x ( 0 ) = 1 , y ( 0 ) = 1

Begin by obtaining a system of algebraic equations for X ( s ) and Y ( s ) :

Writing L { x ( t ) } = X ( s ) and L { y ( t ) } = Y ( s ) you should obtain the set of transformed equations

1 + s Y ( s ) X ( s ) = 0

1 + s X ( s ) + Y ( s ) = 1 s

which, when re-arranged, are

X ( s ) + s Y ( s ) = 1

s X ( s ) + Y ( s ) = 1 s s

Now solve these equations for X ( s ) and Y ( s ) :

X ( s ) = s 1 + s 2 Y ( s ) = 1 s 1 1 + s 2

Now find the required solution by obtaining the inverse Laplace transforms:

You should obtain x ( t ) = cos t . u ( t ) and y ( t ) = ( 1 sin t ) . u ( t ) . This follows since

L 1 { s 1 + s 2 } = cos t . u ( t ) L 1 { 1 s } = u ( t ) L 1 { 1 1 + s 2 } = sin t . u ( t )

Exercises
  1. Solve the given system of differential equations for the initial conditions specified.
    1. d x d t = y d y d t = x x ( 0 ) = 1 y ( 0 ) = 0
    2. d x d t = 4 x 2 y d y d t = 5 x + 2 y x ( 0 ) = 2 y ( 0 ) = 2
  2. The Laplace transform can also be used to solve a pair of coupled second order differential equations.

    Solve, for the given initial conditions,

    d 2 x d t 2 = y + sin t x ( 0 ) = 1 x ( 0 ) = 0

    d 2 y d t 2 = d x d t + cos t y ( 0 ) = 1 y ( 0 ) = 1

    (Note that the initial conditions on each of x ( t ) and y ( t ) are needed in the second order situation.)

    1. x = cosh t , y = sinh t
    2. x = e 3 t ( 2 cos 3 t + 2 sin 3 t ) , y = e 3 t ( 2 cos 3 t + 4 sin 3 t )
  1. x = cos t , y = cos t sin t