Loading [MathJax]/jax/output/HTML-CSS/jax.js

A second order differential equation is one containing the second derivative. These are in general quite complicated, but one fairly simple type is useful: the second order linear equation with constant coefficients.

Example 17.5.1 Consider the intial value problem ¨y˙y2y=0, y(0)=5, ˙y(0)=0. We make an inspired guess: might there be a solution of the form ert? This seems at least plausible, since in this case ¨y, ˙y, and y all involve ert.

If such a function is a solution then r2ertrert2ert=0ert(r2r2)=0(r2r2)=0(r2)(r+1)=0, so r is 2 or 1. Not only are f=e2t and g=et solutions, but notice that y=Af+Bg is also, for any constants A and B: (Af+Bg)(Af+Bg)2(Af+Bg)=Af+BgAfBg2Af2Bg=A(ff2f)+B(gg2g)=A(0)+B(0)=0. Can we find A and B so that this is a solution to the initial value problem? Let's substitute: 5=y(0)=Af(0)+Bg(0)=Ae0+Be0=A+B and 0=˙y(0)=Af(0)+Bg(0)=A2e0+B(1)e0=2AB. So we need to find A and B that make both 5=A+B and 0=2AB true. This is a simple set of simultaneous equations: solve B=2A, substitute to get 5=A+2A=3A. Then A=5/3 and B=10/3, and the desired solution is (5/3)e2t+(10/3)et. You now see why the initial condition in this case included both y(0) and ˙y(0): we needed two equations in the two unknowns A and B

You should of course wonder whether there might be other solutions; the answer is no. We will not prove this, but here is the theorem that tells us what we need to know:

Theorem 17.5.2 Given the differential equation a¨y+b˙y+cy=0, a0, consider the quadratic polynomial ax2+bx+c, called the characteristic polynomial. Using the quadratic formula, this polynomial always has one or two roots, call them r and s. The general solution of the differential equation is:

    (a) y=Aert+Best, if the roots r and s are real numbers and rs.

    (b) y=Aert+Btert, if r=s is real.

    (c) y=Acos(βt)eαt+Bsin(βt)eαt, if the roots r and s are complex numbers α+βi and αβi.

Example 17.5.3 Suppose a mass m is hung on a spring with spring constant k. If the spring is compressed or stretched and then released, the mass will oscillate up and down. Because of friction, the oscillation will be damped: eventually the motion will cease. The damping will depend on the amount of friction; for example, if the system is suspended in oil the motion will cease sooner than if the system is in air. Using some simple physics, it is not hard to see that the position of the mass is described by this differential equation: m¨y+b˙y+ky=0. Using m=1, b=4, and k=5 we find the motion of the mass. The characteristic polynomial is x2+4x+5 with roots (4±1620)/2=2±i. Thus the general solution is y=Acos(t)e2t+Bsin(t)e2t. Suppose we know that y(0)=1 and ˙y(0)=2. Then as before we form two simultaneous equations: from y(0)=1 we get 1=Acos(0)e0+Bsin(0)e0=A. For the second we compute ˙y=2Ae2tcos(t)+Ae2t(sin(t))2Be2tsin(t)+Be2tcos(t), and then 2=2Ae0cos(0)Ae0sin(0)2Be0sin(0)+Be0cos(0)=2A+B. So we get A=1, B=4, and y=cos(t)e2t+4sin(t)e2t.

Here is a useful trick that makes this easier to understand: We have y=(cost+4sint)e2t. The expression cost+4sint is a bit reminiscent of the trigonometric formula cos(αβ)=cos(α)cos(β)+sin(α)sin(β) with α=t. Let's rewrite it a bit as 17(117cost+417sint). Note that (1/17)2+(4/17)2=1, which means that there is an angle β with cosβ=1/17 and sinβ=4/17 (of course, β may not be a "nice'' angle). Then cost+4sint=17(costcosβ+sinβsint)=17cos(tβ). Thus, the solution may also be written y=17e2tcos(tβ). This is a cosine curve that has been shifted β to the right; the 17e2t has the effect of diminishing the amplitude of the cosine as t increases; see figure 17.5.1. The oscillation is damped very quickly, so in the first graph it is not clear that this is an oscillation. The second graph shows a restricted range for t.

Other physical systems that oscillate can also be described by such differential equations. Some electric circuits, for example, generate oscillating current.

1
2
3
4
5
0
1
y
3
4
5
0
0.01
0.02
y
Figure 17.5.1. Graph of a damped oscillation.

Example 17.5.4 Find the solution to the intial value problem ¨y4˙y+4y=0, y(0)=3, ˙y(0)=1. The characteristic polynomial is x24x+4=(x2)2, so there is one root, r=2, and the general solution is Ae2t+Bte2t. Substituting t=0 we get 3=A+0=A. The first derivative is 2Ae2t+2Bte2t+Be2t; substituting t=0 gives 1=2A+0+B=2A+B=2(3)+B=6+B, so B=7. The solution is 3e2t+7te2t.

Exercises 17.5

Ex 17.5.1 Verify that the function in part (a) of theorem 17.5.2 is a solution to the differential equation a¨y+b˙y+cy=0.

Ex 17.5.2 Verify that the function in part (b) of theorem 17.5.2 is a solution to the differential equation a¨y+b˙y+cy=0.

Ex 17.5.3 Verify that the function in part (c) of theorem 17.5.2 is a solution to the differential equation a¨y+b˙y+cy=0.

Ex 17.5.4 Solve the initial value problem ¨yω2y=0, y(0)=1, ˙y(0)=1, assuming ω0. (answer)

Ex 17.5.5 Solve the initial value problem 2¨y+18y=0, y(0)=2, ˙y(0)=15. (answer)

Ex 17.5.6 Solve the initial value problem ¨y+6˙y+5y=0, y(0)=1, ˙y(0)=0. (answer)

Ex 17.5.7 Solve the initial value problem ¨y˙y12y=0, y(0)=0, ˙y(0)=14. (answer)

Ex 17.5.8 Solve the initial value problem ¨y+12˙y+36y=0, y(0)=5, ˙y(0)=10. (answer)

Ex 17.5.9 Solve the initial value problem ¨y8˙y+16y=0, y(0)=3, ˙y(0)=4. (answer)

Ex 17.5.10 Solve the initial value problem ¨y+5y=0, y(0)=2, ˙y(0)=5. (answer)

Ex 17.5.11 Solve the initial value problem ¨y+y=0, y(π/4)=0, ˙y(π/4)=2. (answer)

Ex 17.5.12 Solve the initial value problem ¨y+12˙y+37y=0, y(0)=4, ˙y(0)=0. (answer)

Ex 17.5.13 Solve the initial value problem ¨y+6˙y+18y=0, y(0)=0, ˙y(0)=6. (answer)

Ex 17.5.14 Solve the initial value problem ¨y+4y=0, y(0)=3, ˙y(0)=2. Put your answer in the form developed at the end of example 17.5.3. (answer)

Ex 17.5.15 Solve the initial value problem ¨y+100y=0, y(0)=5, ˙y(0)=50. Put your answer in the form developed at the end of example 17.5.3. (answer)

Ex 17.5.16 Solve the initial value problem ¨y+4˙y+13y=0, y(0)=1, ˙y(0)=1. Put your answer in the form developed at the end of example 17.5.3. (answer)

Ex 17.5.17 Solve the initial value problem ¨y8˙y+25y=0, y(0)=3, ˙y(0)=0. Put your answer in the form developed at the end of example 17.5.3. (answer)

Ex 17.5.18 A mass-spring system m¨y+b˙y+ky has k=29, b=4, and m=1. At time t=0 the position is y(0)=2 and the velocity is ˙y(0)=1. Find y(t). (answer)

Ex 17.5.19 A mass-spring system m¨y+b˙y+ky has k=24, b=12, and m=3. At time t=0 the position is y(0)=0 and the velocity is ˙y(0)=1. Find y(t). (answer)

Ex 17.5.20 Consider the differential equation a¨y+b˙y=0, with a and b both non-zero. Find the general solution by the method of this section. Now let g=˙y; the equation may be written as a˙g+bg=0, a first order linear homogeneous equation. Solve this for g, then use the relationship g=˙y to find y.

Ex 17.5.21 Suppose that y(t) is a solution to a¨y+b˙y+cy=0, y(t0)=0, ˙y(t0)=0. Show that y(t)=0.