To understand that many systems of equations cannot be solved analytically but can be solved using numerical techniques such as Eulerβs Method or Runge-Kutta methods to find approximate solutions for the system.
If we are unable to find an analytic solution to a first-order differential equation , there is no reason to expect that it will be any easier to solve a system of equations. However, many of the numerical techniques that are used to solve a first-order equation can be extended to solve a system.
Large mobile cranes can reach up to 700 feet (Figure 2.3.1). This is about the same height as a 50 story building. At such heights, cranes are susceptible to winds and the end of the boom can move back and forth several feet even in moderate winds. We can model the motion as a harmonic oscillator if the side-to-side motion is not too great. However, if the wind becomes too strong and the crane moves too much, gravity will become a factor and the crane may topple due to its weight.
Before we can construct a model that might best describe the swaying of a large crane, we should remind ourselves of how harmonic oscillators work (see Section 1.1). A simple mass-spring system can be modeled by the equation
where is the displacement of the mass at time . Given such a system, the mass will oscillate forever with constant amplitude. If we want to be a bit more realistic, we can introduce a damping force,
where . Our equation now becomes
If we let and , we can rewrite this last equation as
If we set , we can rewrite this second-order equation as a system of first-order equations,
In the case of our swaying crane, we will let be the displacement of the end of the boom from the perfect vertical position. When , the boom is bent, and the structure of the boom supplies a strong restorative force to bring everything back to a true vertical position. Thus, the swaying of our craneβs boom might be described by an equation such as
or the equivalent system,
We will learn how to solve such systems later, but it is easy to check that
is a solution to our equation. The Sage code for solving this system is given below.
The constants and can be determined if we know the initial position and initial velocity of the end of the craneβs boom. We show several solutions to our equation in Figure 2.3.2.
Modeling the swaying motion of a crane with a harmonic oscillator might work only if the side-to-side motion is small. If the motion is larger, we must account for the effect of gravity in our model. When is large, part of the crane will not be above any other part of the crane. Thus, gravity will pull downward on that part of the crane and will cause the crane to bend even further. We can model this effect by setting the equation for our harmonic oscillator equal to a factor of ,
(2.3.1)
When is small, this forcing term will not contribute much to the motion of the crane. However, when is large, the term will contribute a great deal. The equation is an example of Duffingβs equation.
The system of equations (2.3.2)β(2.3.3) is nonlinear due the term. There is little hope to finding an analytic solution to such a system. We can use software to plot the phase plane in order to learn something about the solutions (Figure 2.3.3). However, we still need to numerically generate solutions in order to plot the phase plane.
Let us see how we might find a numerical solution for a system. Consider the system
with initial conditions and . We can rewrite our system in vector form
where ,,, and . We wish to find approximate values and for the solution and at the points
where is the step size. If we let and , then Eulerβs method now becomes
or
The initial conditions are used to determine , which is the tangent vector to the graph of the solution in the -plane (Figure 2.3.3). We can move in the direction of this tangent vector for time in order to find the next point . We then calculate a new tangent vector and then move along this new vector for a time step to find . We can repeat this technique to generate an approximate solution curve in the phase plane.
The surprising thing about Duffingβs equation is that it is extremely sensitive to initial conditions. A slight change in the initial conditions can yield dramatically different solutions. If we change the initial velocity to , we obtain a very different graph of the solution (Figure 2.3.5 ). For small initial velocities, solutions spiral towards the origin. However, a larger initial velocity will send the solution in the phase plane away from the origin. If the crane sways too violently, we will have a disaster.
Notice that our system is not autonomous and depends on time. Therefore, we cannot graph the phase plane of this system; however, we can graph a solution curve in three dimensions (Figure 2.3.9).
Just as in the case of a single first-order differential equation, we can think of Eulerβs approximation as the first two terms of the Taylor series expansion,
To get a more accurate approximation, we can take the first two terms of the Taylor series. In this case, we must compute . If we return to our example,
we can see how this is done. First, note that
Our algorithm now becomes clear,
Of course, this algorithm requires us to compute some derivatives.
Subsection2.3.4A Word about Existence and Uniqueness
Theorem2.3.10.
Let
be a system of differential equations such that is continuously differentiable. If is the initial value at time , there exists a unique solution for the initial value problem on some open interval for some .
If we have an autonomous system , our solution does not depend on time. Thus, we obtain the same solution curve if we start at the same point even though we might start at different times.
Many systems of equations cannot be solved analytically. However, we can use numerical techniques such as Eulerβs Method to find approximate solutions for the system. Given the system
with initial condition and step size , we can approximate a solution to the system by
We can extend this method to the improved Eulerβs method, Taylor series methods, or Runge-Kutta methods.
Provided certain conditions are met, we are guaranteed unique local solutions to any system of first-order differential equations. Some of the following are consequences of existence and uniqueness.
If we have an autonomous system , our differential equation does not depend on time. Thus, we obtain the same solution curve if we start at the same point even though we might start at different times.