Differential Equations - Systems: Examples

In the introduction to this section we briefly discussed how a system of differential equations can arise from a population problem in which we keep track of the population of both the prey and the predator. It makes sense that the number of prey present will affect the number of the predator present. Likewise, the number of predator present will affect the number of prey present. Therefore the differential equation that governs the population of either the prey or the predator should in some way depend on the population of the other. This will lead to two differential equations that must be solved simultaneously in order to determine the population of the prey and the predator.

The whole point of this is to notice that systems of differential equations can arise quite easily from naturally occurring situations. Developing an effective predator-prey system of differential equations is not the subject of this chapter. However, systems can arise from

n^{th}

order linear differential equations as well. Before we get into this however, let’s write down a system and get some terminology out of the way.

We are going to be looking at first order, linear systems of differential equations. These terms mean the same thing that they have meant up to this point. The largest derivative anywhere in the system will be a first derivative and all unknown functions and their derivatives will only occur to the first power and will not be multiplied by other unknown functions. Here is an example of a system of first order, linear differential equations.

\begin{aligned} {x^{'}}_{1} & = x_{1} + 2 x_{2} \\ {x^{'}}_{2} & = 3 x_{1} + 2 x_{2} \end{aligned}

We call this kind of system a coupled system since knowledge of

x_{2}

is required in order to find

x_{1}

and likewise knowledge of

x_{1}

is required to find

x_{2}

. We will worry about how to go about solving these later. At this point we are only interested in becoming familiar with some of the basics of systems.

Now, as mentioned earlier, we can write an

n^{th}

order linear differential equation as a system. Let’s see how that can be done.

Example 1 Write the following 2^nd order differential equation as a system of first order, linear differential equations.

2 y^{″} - 5 y^{'} + y = 0 y (3) = 6 y^{'} (3) = - 1

We can write higher order differential equations as a system with a very simple change of variable. We’ll start by defining the following two new functions.

\begin{aligned} x_{1} (t) & = y (t) \\ x_{2} (t) & = y^{'} (t) \end{aligned}

Now notice that if we differentiate both sides of these we get,

\begin{aligned} {x^{'}}_{1} & = y^{'} = x_{2} \\ {x^{'}}_{2} & = y^{″} = - \frac{1}{2} y + \frac{5}{2} y^{'} = - \frac{1}{2} x_{1} + \frac{5}{2} x_{2} \end{aligned}

Note the use of the differential equation in the second equation. We can also convert the initial conditions over to the new functions.

\begin{aligned} x_{1} (3) & = y (3) = 6 \\ x_{2} (3) & = y^{'} (3) = - 1 \end{aligned}

Putting all of this together gives the following system of differential equations.

\begin{aligned} {x^{'}}_{1} & = x_{2} & x_{1} (3) & = 6 \\ {x^{'}}_{2} & = - \frac{1}{2} x_{1} + \frac{5}{2} x_{2} & x_{2} (3) & = - 1 \end{aligned}

We will call the system in the above example an Initial Value Problem just as we did for differential equations with initial conditions.

Let’s take a look at another example.

Example 2 Write the following 4^th order differential equation as a system of first order, linear differential equations.

y^{(4)} + 3 y^{″} - \sin (t) y^{'} + 8 y = t^{2} y (0) = 1 y^{'} (0) = 2 y^{″} (0) = 3 y^{‴} (0) = 4

Just as we did in the last example we’ll need to define some new functions. This time we’ll need 4 new functions.

\begin{aligned} x_{1} & = y & \Rightarrow {x^{'}}_{1} & = y^{'} = x_{2} \\ x_{2} & = y^{'} & \Rightarrow {x^{'}}_{2} & = y^{″} = x_{3} \\ x_{3} & = y^{″} & \Rightarrow {x^{'}}_{3} & = y^{‴} = x_{4} \\ x_{4} & = y^{‴} & \Rightarrow {x^{'}}_{4} & = y^{(4)} = - 8 y + \sin (t) y^{'} - 3 y^{″} + t^{2} = - 8 x_{1} + \sin (t) x_{2} - 3 x_{3} + t^{2} \end{aligned}

The system along with the initial conditions is then,

\begin{aligned} {x^{'}}_{1} & = x_{2} & x_{1} (0) & = 1 \\ {x^{'}}_{2} & = x_{3} & x_{2} (0) & = 2 \\ {x^{'}}_{3} & = x_{4} & x_{3} (0) & = 3 \\ {x^{'}}_{4} & = - 8 x_{1} + \sin (t) x_{2} - 3 x_{3} + t^{2} & x_{4} (0) & = 4 \end{aligned}

Now, when we finally get around to solving these we will see that we generally don’t solve systems in the form that we’ve given them in this section. Systems of differential equations can be converted to matrix form and this is the form that we usually use in solving systems.

Example 3 Convert the following system to matrix from.

\begin{aligned} {x^{'}}_{1} & = 4 x_{1} + 7 x_{2} \\ {x^{'}}_{2} & = - 2 x_{1} - 5 x_{2} \end{aligned}

First write the system so that each side is a vector.

(\begin{matrix} {x^{'}}_{1} \\ {x^{'}}_{2} \end{matrix}) = (\begin{matrix} 4 x_{1} + 7 x_{2} \\ - 2 x_{1} - 5 x_{2} \end{matrix})

Now the right side can be written as a matrix multiplication,

(\begin{matrix} {x^{'}}_{1} \\ {x^{'}}_{2} \end{matrix}) = (\begin{matrix} 4 & 7 \\ - 2 & - 5 \end{matrix}) (\begin{matrix} x_{1} \\ x_{2} \end{matrix})

Now, if we define,

\vec{x} = (\begin{matrix} x_{1} \\ x_{2} \end{matrix})

then,

{\vec{x}}^{'} = (\begin{matrix} {x^{'}}_{1} \\ {x^{'}}_{2} \end{matrix})

The system can then be written in the matrix form,

{\vec{x}}^{'} = (\begin{matrix} 4 & 7 \\ - 2 & - 5 \end{matrix}) \vec{x}

Example 4 Convert the systems from Examples 1 and 2 into matrix form.

We’ll start with the system from Example 1.

\begin{aligned} {x^{'}}_{1} & = x_{2} x_{1} (3) = 6 \\ {x^{'}}_{2} & = - \frac{1}{2} x_{1} + \frac{5}{2} x_{2} x_{2} (3) = - 1 \end{aligned}

First define,

\vec{x} = (\begin{matrix} x_{1} \\ x_{2} \end{matrix})

The system is then,

{\vec{x}}^{'} = (\begin{matrix} 0 & 1 \\ - \frac{1}{2} & \frac{5}{2} \end{matrix}) \vec{x} \vec{x} (3) = (\begin{matrix} x_{1} (3) \\ x_{2} (3) \end{matrix}) = (\begin{matrix} 6 \\ - 1 \end{matrix})

Now, let’s do the system from Example 2.

\begin{aligned} {x^{'}}_{1} & = x_{2} x_{1} (0) = 1 \\ {x^{'}}_{2} & = x_{3} x_{2} (0) = 2 \\ {x^{'}}_{3} & = x_{4} x_{3} (0) = 3 \\ {x^{'}}_{4} & = - 8 x_{1} + \sin (t) x_{2} - 3 x_{3} + t^{2} x_{4} (0) = 4 \end{aligned}

In this case we need to be careful with the t² in the last equation. We’ll start by writing the system as a vector again and then break it up into two vectors, one vector that contains the unknown functions and the other that contains any known functions.

(\begin{matrix} {x^{'}}_{1} \\ {x^{'}}_{2} \\ {x^{'}}_{3} \\ {x^{'}}_{4} \end{matrix}) = (\begin{matrix} x_{2} \\ x_{3} \\ x_{4} \\ - 8 x_{1} + \sin (t) x_{2} - 3 x_{3} + t^{2} \end{matrix}) = (\begin{matrix} x_{2} \\ x_{3} \\ x_{4} \\ - 8 x_{1} + \sin (t) x_{2} - 3 x_{3} \end{matrix}) + (\begin{matrix} 0 \\ 0 \\ 0 \\ t^{2} \end{matrix})

Now, the first vector can now be written as a matrix multiplication and we’ll leave the second vector alone.

{\vec{x}}^{'} = (\begin{matrix} 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ - 8 & \sin (t) & - 3 & 0 \end{matrix}) \vec{x} + (\begin{matrix} 0 \\ 0 \\ 0 \\ t^{2} \end{matrix}) \vec{x} (0) = (\begin{matrix} 1 \\ 2 \\ 3 \\ 4 \end{matrix})

where,

\vec{x} (t) = (\begin{matrix} x_{1} (t) \\ x_{2} (t) \\ x_{3} (t) \\ x_{4} (t) \end{matrix})

Note that occasionally for “large” systems such as this we will go one step farther and write the system as,

{\vec{x}}^{'} = A \vec{x} + \vec{g} (t)

The last thing that we need to do in this section is get a bit of terminology out of the way. Starting with

{\vec{x}}^{'} = A \vec{x} + \vec{g} (t)

we say that the system is homogeneous if

\vec{g} (t) = \vec{0}

and we say the system is nonhomogeneous if

\vec{g} (t) \neq \vec{0}

Comments

Digital Signal Processing - Basic Continuous Time Signals

To test a system, generally, standard or basic signals are used. These signals are the basic building blocks for many complex signals. Hence, they play a very important role in the study of signals and systems. Unit Impulse or Delta Function A signal, which satisfies the condition, δ ( t ) = lim ϵ → ∞ x ( t ) δ ( t ) = lim ϵ → ∞ x ( t ) is known as unit impulse signal. This signal tends to infinity when t = 0 and tends to zero when t ≠ 0 such that the area under its curve is always equals to one. The delta function has zero amplitude everywhere except at t = 0. Properties of Unit Impulse Signal δ(t) is an even signal. δ(t) is an example of neither energy nor power (NENP) signal. Area of unit impulse signal can be written as; A = ∫ ∞ − ∞ δ ( t ) d t = ∫ ∞ − ∞ lim ϵ → 0 x ( t ) d t = lim ϵ → 0 ∫ ∞ − ∞ [ x ( t ) d t ] = 1 Weight or strength of the signal can be written as; y ( t ) = A δ ( t ) y ( t ) = A δ ( t ) Area of the weighted impulse s...

Differential Equations - First Order: Bernoulli

In this section we are going to take a look at differential equations in the form, y ′ + p ( x ) y = q ( x ) y n y ′ + p ( x ) y = q ( x ) y n where p ( x ) p ( x ) and q ( x ) q ( x ) are continuous functions on the interval we’re working on and n n is a real number. Differential equations in this form are called Bernoulli Equations . First notice that if n = 0 n = 0 or n = 1 n = 1 then the equation is linear and we already know how to solve it in these cases. Therefore, in this section we’re going to be looking at solutions for values of n n other than these two. In order to solve these we’ll first divide the differential equation by y n y n to get, y − n y ′ + p ( x ) y 1 − n = q ( x ) y − n y ′ + p ( x ) y 1 − n = q ( x ) We are now going to use the substitution v = y 1 − n v = y 1 − n to convert this into a differential equation in terms of v v . As we’ll see th...

Differential Equations - Systems: Solutions

Now that we’ve got some of the basics out of the way for systems of differential equations it’s time to start thinking about how to solve a system of differential equations. We will start with the homogeneous system written in matrix form, → x ′ = A → x (1) (1) x → ′ = A x → where, A A is an n × n n × n matrix and → x x → is a vector whose components are the unknown functions in the system. Now, if we start with n = 1 n = 1 then the system reduces to a fairly simple linear (or separable) first order differential equation. x ′ = a x x ′ = a x and this has the following solution, x ( t ) = c e a t x ( t ) = c e a t So, let’s use this as a guide and for a general n n let’s see if → x ( t ) = → η e r t (2) (2) x → ( t ) = η → e r t will be a solution. Note that the only real difference here is that we let the constant in front of the exponential be a vector. All we need to do then is plug this into the d...

Calculus III - 3-Dimensional Space: Equations of Lines

In this section we need to take a look at the equation of a line in R 3 R 3 . As we saw in the previous section the equation y = m x + b y = m x + b does not describe a line in R 3 R 3 , instead it describes a plane. This doesn’t mean however that we can’t write down an equation for a line in 3-D space. We’re just going to need a new way of writing down the equation of a curve. So, before we get into the equations of lines we first need to briefly look at vector functions. We’re going to take a more in depth look at vector functions later. At this point all that we need to worry about is notational issues and how they can be used to give the equation of a curve. The best way to get an idea of what a vector function is and what its graph looks like is to look at an example. So, consider the following vector function. → r ( t ) = ⟨ t , 1 ⟩ r → ( t ) = ⟨ t , 1 ⟩ A vector function is a function that takes one or more variables, one in this case, and returns a...

Differential Equations - First Order: Modeling - i

We now move into one of the main applications of differential equations both in this class and in general. Modeling is the process of writing a differential equation to describe a physical situation. Almost all of the differential equations that you will use in your job (for the engineers out there in the audience) are there because somebody, at some time, modeled a situation to come up with the differential equation that you are using. This section is not intended to completely teach you how to go about modeling all physical situations. A whole course could be devoted to the subject of modeling and still not cover everything! This section is designed to introduce you to the process of modeling and show you what is involved in modeling. We will look at three different situations in this section : Mixing Problems, Population Problems, and Falling Objects. In all of these situations we will be forced to make assumptions that do not accurately depict reality in most cases, but wi...

Digital Signal Processing - Miscellaneous Signals

There are other signals, which are a result of operation performed on them. Some common type of signals are discussed below. Conjugate Signals Signals, which satisfies the condition x ( t ) = x ∗ ( − t ) are called conjugate signals. Let x ( t ) = a ( t ) + j b ( t ) So, x ( − t ) = a ( − t ) + j b ( − t ) And x ∗ ( − t ) = a ( − t ) − j b ( − t ) By Condition, x ( t ) = x ∗ ( − t ) If we compare both the derived equations 1 and 2, we can see that the real part is even, whereas the imaginary part is odd. This is the condition for a signal to be a conjugate type. Conjugate Anti-Symmetric Signals Signals, which satisfy the condition x ( t ) = − x ∗ ( − t ) are called conjugate anti-symmetric signal Let x ( t ) = a ( t ) + j b ( t ) So x ( − t ) = a ( − t ) + j b ( − t ) And x ∗ ( − t ) = a ( − t ) − j b ( − t ) − x ∗ ( − t ) = − a ( − t ) + j b ( − t ) By Condition x ( t ) = − x ∗ ( − t ) ...

Differential Equations - Systems: Repeated Eigenvalues - i

This is the final case that we need to take a look at. In this section we are going to look at solutions to the system, → x ′ = A → x x → ′ = A x → where the eigenvalues are repeated eigenvalues. Since we are going to be working with systems in which A A is a 2 × 2 2 × 2 matrix we will make that assumption from the start. So, the system will have a double eigenvalue, λ λ . This presents us with a problem. We want two linearly independent solutions so that we can form a general solution. However, with a double eigenvalue we will have only one, → x 1 = → η e λ t x → 1 = η → e λ t So, we need to come up with a second solution. Recall that when we looked at the double root case with the second order differential equations we ran into a similar problem. In that section we simply added a t t to the solution and were able to get a second solution. Let’s see if the same thing will work in this case as well. We’ll see if → x = t e...

Differential Equations - Systems: Repeated Eigenvalues - ii

Example 3 Solve the following IVP. → x ′ = ( − 1 3 2 − 1 6 − 2 ) → x → x ( 2 ) = ( 1 0 ) x → ′ = ( − 1 3 2 − 1 6 − 2 ) x → x → ( 2 ) = ( 1 0 ) First the eigenvalue for the system. det ( A − λ I ) = ∣ ∣ ∣ ∣ − 1 − λ 3 2 − 1 6 − 2 − λ ∣ ∣ ∣ ∣ = λ 2 + 3 λ + 9 4 = ( λ + 3 2 ) 2 ⇒ λ 1 , 2 = − 3 2 det ( A − λ I ) = | − 1 − λ 3 2 − 1 6 − 2 − λ | = λ 2 + 3 λ + 9 4 = ( λ + 3 2 ) 2 ⇒ λ 1 , 2 = − 3 2 Now let’s get the eigenvector. ( 1 2 3 2 − 1 6 − 1 2 ) ( η 1 η 2 ) = ( 0 0 ) ⇒ 1 2 η 1 + 3 2 η 2 = 0 η 1 = − 3 η 2 ( 1 2 3 2 − 1 6 − 1 2 ) ( η 1 η 2 ) = ( 0 0 ) ⇒ 1 2 η 1 + 3 2 η 2 = 0 η 1 = − 3 η 2 → η = ( − 3 η 2 η 2 ) η 2 ≠ 0 → η ( 1 ) = ( − 3 1 ) η 2 = 1 η → = ( − 3 η 2 η 2 ) η 2 ≠ 0 η → ( 1 ) = ( − 3 1 ) η 2 = 1 Now find → ρ ρ → , ( 1 2 3 2 − 1 6 − 1 2 ) ( ρ 1 ρ 2 ) = ( − 3 1 ) ⇒ 1 2 ρ 1 + 3 2 ρ 2 = − 3 ρ 1 = − 6 − 3 ρ 2 ( 1 2 3 2 − 1 6 − 1 2 ) ( ρ 1 ρ 2 ) = ( − 3 1 ) ⇒ 1 2 ρ 1 + 3 2 ρ 2 = − 3 ρ 1 = − 6 − 3 ρ 2 → ρ = ( − 6 − 3 ρ 2 ρ 2 ) ⇒ → ρ = ( − 6 0 ) if ρ 2 = 0 ρ → ...

Differential Equations - Second Order: Repeated Roots

In this section we will be looking at the last case for the constant coefficient, linear, homogeneous second order differential equations. In this case we want solutions to a y ′′ + b y ′ + c y = 0 a y ″ + b y ′ + c y = 0 where solutions to the characteristic equation a r 2 + b r + c = 0 a r 2 + b r + c = 0 are double roots r 1 = r 2 = r r 1 = r 2 = r . This leads to a problem however. Recall that the solutions are y 1 ( t ) = e r 1 t = e r t y 2 ( t ) = e r 2 t = e r t y 1 ( t ) = e r 1 t = e r t y 2 ( t ) = e r 2 t = e r t These are the same solution and will NOT be “nice enough” to form a general solution. We do promise that we’ll define “nice enough” eventually! So, we can use the first solution, but we’re going to need a second solution. Before finding this second solution let’s take a little side trip. The reason for the side trip will be clear eventually. From the quadratic formula we know that the roots to the characteristic equation are, r 1 , 2 = ...

Differential Equations - Laplace Transforms: Table

f ( t ) = L − 1 { F ( s ) } f ( t ) = L − 1 { F ( s ) } F ( s ) = L { f ( t ) } F ( s ) = L { f ( t ) } 1 1 s 1 s e a t e a t 1 s − a 1 s − a t n , n = 1 , 2 , 3 , … t n , n = 1 , 2 , 3 , … n ! s n + 1 n ! s n + 1 t p t p , p > − 1 p > − 1 Γ ( p + 1 ) s p + 1 Γ ( p + 1 ) s p + 1 √ t t √ π 2 s 3 2 π 2 s 3 2 t n − 1 2 , n = 1 , 2 , 3 , … t n − 1 2 , n = 1 , 2 , 3 , … 1 ⋅ 3 ⋅ 5 ⋯ ( 2 n − 1 ) √ π 2 n s n + 1 2 1 ⋅ 3 ⋅ 5 ⋯ ( 2 n − 1 ) π 2 n s n + 1 2 sin ( a t ) sin ⁡ ( a t ) a s 2 + a 2 a s 2 + a 2 cos ( a t ) cos ⁡ ( a t ) s s 2 + a 2 s s 2 + a 2 t sin ( a t ) t sin ⁡ ( a t ) 2 a s ( s 2 + a 2 ) 2 2 a s ( s 2 + a 2 ) 2 t cos ( a t ) t cos ⁡ ( a t ) s 2 − a 2 ( s 2 + a 2 ) 2 s 2 − a 2 ( s 2 + a 2 ) 2 sin ( a t ) − a t cos ( a t ) sin ⁡ ( a t ) − a t cos ⁡ ( a t ) 2 a 3 ( s 2 + a 2 ) 2 2 a 3 ( s 2 + a 2 ) 2 sin ( a t ) + a t cos ( a t ) sin ⁡ ( a t ) + a t cos ⁡ ( a t ) 2 a s 2 ( s 2 + a 2 ) 2 2 a s 2 ( s 2 + a 2 ) 2 cos ( a t ) − a t sin ( a t ) cos ⁡ (...

Robert Karamagi

Search This Blog

Calculus III - 3-Dimensional Space: Equations of Lines