Skip to main content

The Ordinary Differential Equations Project

Thomas W. Judson

Section 4.4 Forcing and Resonance

Resonance happens when the natural frequency and the forcing frequency of an undamped harmonic oscillator are the same. There are many real world examples of resonance. For example, it is possible for a trained singer to shatter a champagne glass by simply singing the correct note (see www.youtube.com/watch?v=IZD8ffPwXRo). Shattering a glass by singing the correct note can be modeled by an equation of an undamped harmonic oscillator with sinusoidal forcing. While a certain amount of damping occurs in all physical systems, the damping is often so small that it makes sense to use an undamped harmonic oscillator.

Subsection 4.4.1 Resonance

Resonance was responsible for the collapse of the Broughton suspension bridge near Manchester, England in 1831. The collapse occurred when a column of soldiers marched in cadence over the bridge, setting up a periodic force of rather large amplitude. The frequency of the force was approximately equal to the natural frequency of the bridge. Thus, the bridge collapsed when large oscillations occurred. For this reason soldiers are ordered to break cadence whenever they cross a bridge.
The Millennium Bridge, the first new bridge to span the Thames River in London in over 100 years, is a modern example of how resonance can effect a bridge (Figure 4.4.1). This pedestrian bridge, which opened to the public in June 2000, was quickly closed after the bridge experienced high amplitude horizontal oscillations during periods of high traffic. Studies by designers found that the bridge experienced high amplitude horizontal oscillations in response to horizontal forcing at a rate of one cycle per second. Typically, people walk at a rate of two steps per second, so the time between two successive steps of the left foot is about one second. Thus, if people were to walk in cadence, they could set up strong horizontal forcing that would place a destructive load on the bridge. The engineers did not envision this to be a problem since tourists do not generally march in time. However, a video of tourists crossing the bridge revealed the opposite. When the bridge began oscillating, people tended to walk in cadence in order to keep their balance (www.youtube.com/watch?v=gQK21572oSU).
a photo of the Millennium Bridge in London
Figure 4.4.1. The Millennium Bridge
Consider the case of an undamped harmonic oscillator with a sinusoidal forcing function. We know that we can model this situation with the second-order linear differential equation
\begin{equation*} x'' + \omega_0^2 x = A \cos \omega t. \end{equation*}
Since the homogeneous solution to this equation is
\begin{equation*} x_h(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t, \end{equation*}
we can see that the natural frequency of the harmonic oscillator is \(\omega_0\text{.}\) This is the frequency of the oscillations if there is no forcing term. Depending on whether or not the frequency of the forcing term \(\omega\) is equal to the natural frequency, we will have two very different situations. If the two frequencies are equal, then they will re-enforce one another and we will have resonance. If the two frequencies are close but not equal, we will see the phenomenon of beats. We wish to investigate what happens when the driving frequency \(\omega\) is equal to the natural frequency and when it is not.
Let us examine the case where the forcing frequency and the natural frequency of the oscillator are the same,
\begin{equation} x'' + \omega_0^2 x = A \cos \omega_0 t.\tag{4.4.1} \end{equation}
Since \(A \cos \omega_0 t\) is a solution to the homogeneous equation \(x'' + \omega_0^2 x = 0\text{,}\) we cannot assume that a particular solution to equation (4.4.1) has the form \(a \cos \omega_0 t + b \sin \omega_0 t\text{.}\) Equivalently, if we use the complex method, we cannot assume that our solution has the form \(ae^{i \omega_0t}\text{.}\) Therefore, we will look for a complex solution of the form
\begin{equation*} x_c = a t e^{i \omega_0 t} \end{equation*}
for the complex equation that corresponds to (4.4.1),
\begin{equation} x'' + \omega_0^2 x = A e^{i \omega_0 t}.\tag{4.4.2} \end{equation}
In this case,
\begin{align*} x_c' & = a ( 1 + i \omega_0 t )e^{i \omega_0 t} \\ x_c'' & = a (2 i \omega_0 - \omega_0^2 t) e^{i \omega_0 t}. \end{align*}
Substituting \(x_c\) and \(x_c''\) into the left-hand side of (4.4.2), we have
\begin{align*} x_c'' + \omega_0^2 x_c & = a (2 i \omega_0 - \omega_0^2 t )e^{i \omega_0 t} + \omega_0^2 a t e^{i \omega_0 t}\\ & = 2ai \omega_0 e^{i \omega t}. \end{align*}
In order for \(x_c = a t e^{i \omega_0 t}\) to be a solution to (4.4.2), we must have
\begin{equation*} a = \frac{A}{2 i \omega_0}. \end{equation*}
Thus, our solution to the complex form of the differential equation is
\begin{align*} x_c & = \frac{A}{2i \omega_0} t e^{i \omega_0 t}\\ & = - \frac{Ai}{2 \omega_0} t e^{i \omega_0 t}\\ & = - \frac{Ai}{2 \omega_0} t (\cos \omega_0 t + i \sin \omega_0 t)\\ & = \frac{A}{2 \omega_0} t \sin \omega_0 t - i \frac{A}{2 \omega_0} t \cos \omega_0 t. \end{align*}
The real part of \(x_c\text{,}\)
\begin{equation*} x_p(t) = \frac{A}{2 \omega_0} t \sin \omega_0 t, \end{equation*}
is a particular solution to (4.4.1). Thus, our general solution is
\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t + \frac{A}{2 \omega_0} t \sin \omega_0 t. \end{equation*}

Example 4.4.2.

Now let us consider the initial value problem
\begin{align*} x'' + 144 x & = 4 \cos 12t\\ x(0) & = 0\\ x'(0) & = 0. \end{align*}
where \(A = 4\) and \(\omega_0 = 12\text{.}\) The solution to the homogeneous equation \(x'' + 144 x = 0\) is
\begin{equation*} x_h(t) = c_1 \cos 12t + c_2 \sin 12t. \end{equation*}
To find a particular solution to \(x'' + 144 x = 4 \cos 12t\text{,}\) we will use the complex method and try to find a particular solution to \(x'' + 144 x = 4 e^{12it}\text{.}\) We must assume that the solution has the form \(x_c(t) = a t e^{12it}\text{,}\) since \(a e^{12it}\) is a solution to the homogeneous equation. As before, we have \(x_c'' = a (24i - 144 t) e^{12it}\text{.}\) If we substitute \(x_c\) and \(x_c''\) into the left-hand side of our differential equation, we have
\begin{equation*} x_c'' + 144 x_c = a (24i - 144 t) e^{12it} + 144 a t e^{12i t} = 24ai e^{12i t}. \end{equation*}
Therefore,
\begin{equation*} a = \frac{4}{24i} = - \frac{1}{6}i, \end{equation*}
and our complex solution is
\begin{equation*} - \frac{1}{6} i t e^{12it} = - \frac{1}{6} i t (\cos 12t + i \sin 12 t) = \frac{1}{6} t \sin 12 t - \frac{1}{6} i t \cos 12 t. \end{equation*}
Taking the real part of our complex solution, we have a particular solution
\begin{equation*} x_p(t) = \frac{1}{6} t \sin 12t. \end{equation*}
Thus, the general solution to \(x'' + 144 x = 4 \cos 12t\) is
\begin{equation*} x(t) = x_h(t) + x_p(t)= c_1 \cos 12t + c_2 \sin 12t + \frac{1}{6} t \sin 12t. \end{equation*}
Applying the initial conditions, both \(c_1 = 0\) and \(c_2 = 0\text{.}\) Consequently, the solution to the initial value problem is
\begin{equation*} x(t) = \frac{1}{6} t \sin 12t. \end{equation*}
The graph of this solution is given in Figure 4.4.3. Notice that our solution grows with time. This growth is due to the fact that the frequency of the forcing term is equal to the natural frequency of the oscillator. Since the force pulls and pushes at a frequency equal to the natural frequency of the oscillator, the amplitude increases with time. This type of behavior is called resonance.
an oscillating position curve with amplitudes increasing as time increases
Figure 4.4.3. The solution to \(x'' + 4 x = 4 \sin 2t\text{,}\) \(x(0) = 0\text{,}\) \(x'(0) = 0\text{.}\)

Activity 4.4.1. Resonance.

Solve each of the following initial value problems and plot the solution for \(t \geq 0\text{.}\) In each case, what is the natural (and forcing) frequency of the harmonic oscillator?
(a)
\begin{align*} x'' + 121 x & = 3 \cos 11t\\ x(0) & = 0\\ x'(0) & = 0. \end{align*}
(b)
\begin{align*} x'' + 169 x & = 5 \cos 13t\\ x(0) & = 0\\ x'(0) & = 0. \end{align*}

Subsection 4.4.2 Beats or the Case \(\omega \neq \omega_0\)

If the forcing frequency and natural frequency are not equal but close, then we have the phenomenon of beats. For example, a piano tuner uses this phenomenon to make certain that a particular string is correctly tuned. The tuner will strike a tuning fork which vibrates at the correct frequency. Next, the tuner hits the poorly tuned piano string. This modulation gives rise to beats in the tone that are readily audible. When the string is properly tuned, the beats will go away. To see another example of beats, watch the video www.youtube.com/watch?v=pRpN9uLiouI.
To understand the situation where the forcing and natural frequencies differ by a small amount, \(\omega \neq \omega_0\text{,}\) we will first consider the equation,
\begin{equation} x'' + \omega_0^2 x = A \cos \omega t.\tag{4.4.3} \end{equation}
The solution to the homogeneous equation \(x'' + \omega_0^2 x = 0\) is
\begin{equation*} x_h(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t. \end{equation*}
Since the forcing frequency is not equal to the natural frequency, we can look for a particular solution to the complex of the form of equation (4.4.3),
\begin{equation} x'' + \omega_0^2 = A e^{i \omega t}.\tag{4.4.4} \end{equation}
If we assume that our solution has the form \(x_c = a e^{i\omega t}\) and substitute \(x_c\) into the left-hand side of (4.4.4), we find
\begin{align*} x_c'' + \omega_0^2 x_c \amp = -a \omega^2 e^{i \omega t} + \omega_0^2 a e^{i \omega t}\\ \amp = a(\omega_0^2 - \omega^2)e^{i \omega t}. \end{align*}
Therefore, for \(x_c = a e^{i\omega t}\) to be a solution for
\begin{equation*} x'' + \omega_0^2 x = Ae^{i\omega t}, \end{equation*}
we must have
\begin{equation*} a = \frac{A}{\omega_0^2 - \omega^2}. \end{equation*}
The real part of \(x_c\) is a particular solution to (4.4.3),
\begin{equation*} x_p(t) = \frac{A}{\omega_0^2 - \omega^2} \cos \omega t. \end{equation*}
Thus, the general solution to (4.4.3) is
\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 \cos \omega_0 t + c_2 \sin \omega_0 t + \frac{A}{\omega_0^2 - \omega^2} \cos \omega t. \end{equation*}
Now let us examine the case where the motion starts at equilibrium. That is, we will investigate what happens when \(x(0) = 0\) and \(x'(0) = 0\text{.}\) In this case, we can easily determine that \(c_1 = -A/(\omega_0^2 - \omega^2)\) and \(c_2 = 0\text{.}\) Hence the solution to our initial value problem is
\begin{equation} x(t) = \frac{A}{\omega_0^2 - \omega^2} (\cos \omega t - \cos \omega_0 t).\tag{4.4.5} \end{equation}
Thus, we have a superposition to two oscillations of different frequencies.
To understand how the superposition of the two frequencies works, we will let
\begin{equation*} \overline{\omega} = \frac{\omega_0 + \omega}{2} \end{equation*}
be the mean frequency and
\begin{equation*} \delta = \frac{\omega_0 - \omega}{2} \end{equation*}
be the half difference. Since \(\omega = \overline{\omega} - \delta\) and \(\omega_0 = \overline{\omega} + \delta\text{,}\) we know that
\begin{align*} x(t) & = \frac{A}{\omega_0^2 - \omega^2} (\cos \omega t - \cos \omega_0 t)\\ & = \frac{A}{(\overline{\omega} + \delta)^2 - (\overline{\omega} - \delta)^2} (\cos [(\overline{\omega} - \delta) t] - \cos[(\overline{\omega} + \delta) t])\\ & = \frac{A}{4\overline{\omega} \delta} [\cos(\overline{\omega} t - \delta t) - \cos(\overline{\omega}t + \delta t)]\\ & = \frac{A}{4\overline{\omega} \delta} [\cos(\overline{\omega} t) \cos( \delta t) + \sin(\overline{\omega} t) \sin( \delta t) - \cos(\overline{\omega} t) \cos( \delta t) + \sin(\overline{\omega} t) \sin( \delta t)]\\ & = \frac{A}{4\overline{\omega} \delta} [2 \sin(\overline{\omega} t) \sin( \delta t)]\\ & = \frac{A \sin \delta t}{2 \overline{\omega} \delta} \sin \overline{\omega} t. \end{align*}
Thus, we have a sine function, \(\sin \overline{\omega} t\) whose amplitude varies according to \(|(A/2 \overline{\omega} \delta) \sin \delta t|\text{,}\) and
\begin{equation} x(t) = \frac{A \sin \delta t}{2 \overline{\omega} \delta} \sin \overline{\omega} t.\tag{4.4.6} \end{equation}

Example 4.4.4.

In the initial value problem
\begin{align} x'' + 144 x \amp = 4 \cos(13t)\tag{4.4.7}\\ x(0) \amp = 0\tag{4.4.8}\\ x'(0) \amp = 0,\tag{4.4.9} \end{align}
we have \(A= 4\text{,}\) \(\omega_0 = 12\text{,}\) and \(\omega = 13\text{.}\) The solution to this equation is
\begin{equation*} x(t) = \frac{4}{25}(\cos 12 t - \cos 13 t). \end{equation*}
The graph of this solution is given in Figure 4.4.5 and illustrate the phenomenon of beats. The two frequencies are almost equal and interfere with each other.
an oscillating solution curve with the amplitudes also oscillating
Figure 4.4.5. Solution to the initial value problem (4.4.7)
In our example, \(\delta = -1/2\) and \(\overline{\omega} = 25/2\text{.}\) The factor \(\sin \delta t = \sin( - t/2)\) in (4.4.6) oscillates very slowly in comparison to \(\sin(\overline{\omega}) = \sin(25t/2)\text{.}\) Thus, we can view the solution (4.4.6) as a fast oscillation with a frequency \(\overline{\omega}\) and amplitude
\begin{equation} \left| \frac{A \sin \delta t}{2 \overline{\omega} \delta} \right| = \left| \frac{8}{25} \sin\left( \frac{t}{2}\right)\right|,\tag{4.4.10} \end{equation}
which oscillates much more slowly. If we superimpose the graph of (4.4.10) onto the graph in Figure 4.4.5, we obtain an envelope of the faster oscillation.
an oscillating solution curve with the amplitudes also oscillating and an envelope curve that bounds the oscillating amplitudes
Figure 4.4.6. Envelope for the initial value problem (4.4.7)

Activity 4.4.2. Beats.

Consider the initial value problem
\begin{align*} x'' + 196 x \amp = 4 \cos(15t)\\ x(0) \amp = 0\\ x'(0) \amp = 0. \end{align*}
(a)
Solve the initial value problem and plot its solution.
(b)
Calculate the mean frequency, \(\overline{\omega} = (\omega_0 + \omega)/2\text{,}\) and the half difference, \(\delta = (\omega_0 - \omega)/2\text{,}\) where \(\omega_ 0\) is the natural frequency and \(\omega\) is the forcing frequency and rewrite the solution to the initial value problem in the form (4.4.6).
(c)
Plot the envelope, \(|(A \sin \delta t)/(2 \overline{\omega} \delta)|\text{,}\) over the plot of the solution to the initial value problem.
(d)
Discuss what is going on in this activity.

Subsection 4.4.3 Forced Damped Harmonic Motion

Since perpetual motion does not really exist, we should consider what happens when we add some damping to our harmonic oscillator,
\begin{equation} x'' + 2c x' + \omega_0^2 x = A \cos \omega t.\tag{4.4.11} \end{equation}

Example 4.4.7.

Let us consider the differential equation
\begin{equation} x'' + 2x' + 2x = 3 \cos 4t.\tag{4.4.12} \end{equation}
The corresponding homogeneous equation is
\begin{equation*} x'' + 2x' + 2x = 0, \end{equation*}
which has solution
\begin{equation*} x_h(t) = c_1 e^{-t} \cos t + c_2 e^{-t} \sin t. \end{equation*}
To find a particular solution for the equation (4.4.12), we will look for a solution of the form \(x_c = a e^{4 i t}\) to the equation
\begin{equation} x'' + 2x' + 2x = 3e^{4 it}\tag{4.4.13} \end{equation}
and then set \(x_p(t)\) to the real part of the solution. Substituting \(x_c\) into the lefthand side of (4.4.13), we obtain
\begin{equation*} x_c'' + 2x_c' + 2 x_c = (-14 + 8i) a e^{4it} = P(4i)a e^{4it} = 3e^{4 it}, \end{equation*}
where \(P(\lambda) = \lambda^2 + 2 \lambda + 2\) is the characteristic polynomial of \(x'' + 2x' + 2x = 0\text{.}\) Therefore,
\begin{equation*} a = \frac{3}{P(4i)} = \frac{3}{-14 + 8i} = - \frac{3}{130}(7 + 4i), \end{equation*}
and the solution to the complex equation is
\begin{align*} x_c \amp = - \frac{3}{130}(7 + 4i)e^{4 i t}\\ \amp = - \frac{3}{130}(7 + 4i)(\cos 4t + i \sin 4t)\\ \amp = \frac{3}{130}[(- 7 \cos 4t + 4 \sin 4t) + i(- 4 \cos 4t - 7 \sin 4t)]. \end{align*}
Thus, we can determine a particular solution by taking the real part of \(x_c\text{,}\)
\begin{equation*} x_p(t) = \frac{3}{130}(- 7 \cos 4t + 4 \sin 4t). \end{equation*}
The general solution for (4.4.12) is
\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 e^{-t} \cos t + c_2 e^{-t} \sin t + \frac{3}{130}(- 7 \cos 4t + 4 \sin 4t). \end{equation*}
It is easy to see that all solutions approach \(x_p(t)\) as \(t \to \infty\) no matter what the initial conditions may be. However, it would be much more useful if we would write the particular solution in the form
\begin{equation*} x_p(t) = K \cos(4 t - \phi ) \end{equation*}
for some constants \(K\) and \(\phi\) rather than a linear combination of \(\cos 4t\) and \(\sin 4t\text{.}\)
Returning to the general case, the associated homogeneous equation for (4.4.11) is
\begin{equation} x'' + 2c x' + \omega_0^2 x = 0.\tag{4.4.14} \end{equation}
The characteristic polynomial of (4.4.14),
\begin{equation} P(\lambda ) = \lambda^2 + 2c \lambda + \omega_0^2,\tag{4.4.15} \end{equation}
has roots
\begin{equation*} \lambda = - c \pm \sqrt{c^2 - \omega_0^2}. \end{equation*}
If our harmonic oscillator is undedamped (\(c \lt \omega_0\)), then the solution to the homogeneous equation (4.4.14) is
\begin{equation*} x_h = e^{-ct} (c_1 \cos \beta t + c_2 \sin \beta t), \end{equation*}
where \(\beta = \sqrt{\omega_0^2 - c^2}\text{.}\)
As in Example 4.4.7, we will look for a particular solution to (4.4.11) by examining the equivalent complex differential equation. That is, we shall look for a solution of the form \(x_c = a e^{i \omega t}\) to the equation
\begin{equation} x_c'' + 2c x_c' + \omega_0^2 x_c = Ae^{i \omega t},\tag{4.4.16} \end{equation}
and then set \(x_p\) equal to the real part of our solution. Substituting \(x_c\) into the right-hand side of (4.4.16), we obtain
\begin{equation*} x_c'' + 2c x_c' + \omega_0^2 x_c = [(i\omega)^2 + 2c(i \omega) + \omega_0^2] a e^{i \omega t} = P(i \omega) a e^{i \omega t}, \end{equation*}
where \(P\) is the characteristic polynomial (4.4.15). Thus, equation (4.4.16) becomes
\begin{equation*} P(i \omega) a e^{i \omega t} = Ae^{i \omega t}, \end{equation*}
and
\begin{equation*} x_c(t) = a e ^{i \omega t} = \frac{A}{P(i \omega)} e^{i \omega t} = H(i \omega) A e^{i \omega t}, \end{equation*}
where
\begin{equation*} H(i \omega) = \frac{1}{P(i \omega)}. \end{equation*}
We say that \(H(\lambda) = 1 / P(\lambda)\) is the transfer function.
Let us examine the transfer function more closely to see if we can use it to write our particular solution in a more useable form. First, let us write
\begin{equation*} P(i \omega ) = (i\omega)^2 + 2c(i \omega) + \omega_0^2 = (\omega_0^2 - \omega^2) + 2ic \omega. \end{equation*}
in polar form,
\begin{equation*} P(i \omega) = R e^{i \phi} = R(\cos \phi + i \sin \phi), \end{equation*}
where
\begin{equation*} R = \sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2} \end{equation*}
and \(\phi\) is the angle defined by the equations
\begin{align*} \cos \phi & = \frac{\omega_0^2 - \omega^2}{\sqrt{(\omega_0^2 - \omega^2)^2 - 4 c^2 \omega^2}},\\ \sin \phi & = \frac{2c\omega}{\sqrt{(\omega_0^2 - \omega^2)^2 - 4 c^2 \omega^2}}. \end{align*}
Since \(2 c \omega \gt 0\text{,}\) we know that \(\sin \phi \gt 0\text{.}\) Equivalently, \(0 \lt \phi \lt \pi\text{.}\) Thus,
\begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega} \right). \end{equation*}
Therefore, we can write the transfer function as
\begin{equation*} H(i \omega) = \frac{1}{P(i \omega)} = \frac{1}{R} e^{-i \phi}. \end{equation*}
We define the gain to be
\begin{equation*} G(\omega) = \frac{1}{R} = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}, \end{equation*}
and we will rewrite the transfer function as
\begin{equation*} H(i \omega) = G(\omega) e^{-i \phi(\omega)}. \end{equation*}
Thus, the solution to
\begin{equation*} x_c'' + 2c x_c' + \omega_0^2 x_c = Ae^{i \omega t} \end{equation*}
is
\begin{equation} x_c(t) = H(i \omega) A e^{i \omega t} = G(\omega) A e^{i(\omega t - \phi)}.\tag{4.4.17} \end{equation}
Taking the real part of (4.4.17), our particular solution is
\begin{equation*} x_p(t) = \real(x_c(t)) = G(\omega) A \cos(\omega t - \phi ). \end{equation*}
It is now clear that particular solution has the same frequency as the forcing term. In addition, \(x_p\) is out of phase with the driving force by the amount
\begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega}\right). \end{equation*}
The general solution to
\begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}
is
\begin{equation*} x(t) = x_h(t) + x_p(t) = e^{-ct} (c_1 \cos(\eta t) + c_2 \sin(\eta t)) + G(\omega) A \cos(\omega t - \phi). \end{equation*}
Since \(x_h\) has the factor \(e^{-ct}\text{,}\) the homogeneous part of the solution quickly decays to zero as \(t \to \infty\text{.}\) For this reason, \(x_h\) is called the transient term while \(x_p\) is called the steady-state term.

Example 4.4.8.

Let us examine the steady-state solution of the harmonic oscillator in Example 4.4.7. We have already determined that the second-order linear differential equation
\begin{equation*} x'' + 2x' + 2x = 3 \cos 4t \end{equation*}
has solution
\begin{equation*} x(t) = x_h(t) + x_p(t) = c_1 e^{-t} \cos t + c_2 e^{-t} \sin t + \frac{3}{130}(- 7 \cos 4t + 4 \sin 4t). \end{equation*}
The natural frequency is \(\omega_0 = \sqrt{2} \approx 1.4142\) If we choose the initial conditions so that \(c_1 = 0\) and \(c_2 = 1\text{,}\) our transient term is
\begin{equation*} x_h(t) = e^{-ct} \sin \beta t = e^{-t} \sin t. \end{equation*}
Since \(c = 1\text{,}\) \(\omega = 4\text{,}\) and \(A = 3\text{,}\) the gain is \(G(4) \approx 0.0620\text{,}\) the amplitude of the steady state solution is \(G(4) \cdot 4 \approx 0.2481\text{.}\) The phase is \(\phi = \phi(4) \approx 2.6224.\) The solution and the steady-state solution are given in Figure 4.4.9. As we can see, transient solutions can be quite large when compared to steady-state solutions.
an oscillating steady-state solution and a solution curve to a forced damped harmonic oscillator that converges to the steady-state solution
Figure 4.4.9. A forced damped harmonic oscillator and the steady-state solution.
Large transient currents can be destructive in electrical circuits. For example, a light bulb usually burns out when a large transient flows through an already weakened bulb. This, usually occurs at the time the light bulb is turned on or turned off. Large transient currents are particularly harmful to the hard drive in a computer. For this reason, computers should be left on except in the case that they are not being used for a period of several days.

Subsection 4.4.4 Important Lessons

  • A forced damped harmonic oscillator,
    \begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}
    has homogeneous solution
    \begin{equation*} x_h = e^{-ct} (c_1 \cos \eta t + c_2 \sin \eta t), \end{equation*}
    where \(\eta = \sqrt{\omega_0^2 - c^2}\text{.}\) A particular solution is given by
    \begin{equation*} x_p(t) = \real(x_c(t)) = G(\omega) A \cos(\omega t - \phi ), \end{equation*}
    where the phase angle \(\phi\) is
    \begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega} \right) \end{equation*}
    and the gain \(G\) is given by
    \begin{equation*} G(\omega) = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}. \end{equation*}
  • The amplitude and phase of the steady-state solution,
    \begin{equation*} x_p(t) = G(\omega) A \cos(\omega t - \phi), \end{equation*}
    are determined by
    \begin{equation*} \omega_0^2G = \frac{1}{\sqrt{(1- s^2)^2 + D^2 s^2}}, \end{equation*}
    where \(\omega = s \omega_0\) and \(D = 2c/\omega_0\text{.}\)
  • In an undamped harmonic oscillators with a sinusoidal forcing term, the natural frequency of the solution interacts with the frequency of the forcing term. As the frequency of the forcing term approaches the natural frequency of the equation, we can observe a phenomenon called resonance. We can use the equation
    \begin{equation*} x'' + \omega_0^2 x = A \cos \omega t \end{equation*}
    to model an undamped harmonic oscillator with sinusoidal forcing.
  • If the driving frequency is not equal to the natural frequency (\(\omega \neq \omega_0\)), we have a particular solution
    \begin{equation*} x_p(t) = \frac{A}{\omega_0^2 - \omega^2} (\cos \omega t - \cos \omega_0 t). \end{equation*}
    If the two frequencies are close, we can see a phenomenon called beats.
  • If the driving frequency is equal to the natural frequency, we have a particular solution
    \begin{equation*} x_p(t) = \frac{A}{2 \omega_0} t \sin \omega_0 t. \end{equation*}
    In this case, we can see the resonance phenomenon.
  • A forced damped harmonic oscillator,
    \begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}
    has homogeneous solution
    \begin{equation*} x_h = e^{-ct} (c_1 \cos \eta t + c_2 \sin \eta t), \end{equation*}
    where \(\eta = \sqrt{\omega_0^2 - c^2}\text{.}\) A particular solution is given by
    \begin{equation*} x_p(t) = \real(x_c(t)) = G(\omega) A \cos(\omega t - \phi), \end{equation*}
    where the phase angle \(\phi\) is
    \begin{equation*} \phi = \phi(\omega) = \cot^{-1} \left( \frac{\omega_0^2 - \omega^2}{2c \omega} \right) \end{equation*}
    and the gain \(G\) is given by
    \begin{equation*} G(\omega) = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}. \end{equation*}

Reading Questions 4.4.5 Reading Questions

1.

Describe what resonance means in your own words.

2.

Describe in words alone what the steady state solution of a damped harmonic oscillator is.

Exercises 4.4.6 Exercises

Beats and Resonance.

Solve the initial problems in Exercise Group 4.4.6.1–6. Do you find any beats or resonance?
1.
\(y'' + y = \cos t\text{,}\) \(y(0) = 0\text{,}\) \(y'(0) = 0\)
2.
\(y'' + y = \cos(11t/10)\text{,}\) \(y(0) = 0\text{,}\) \(y'(0) = 0\)
3.
\(y'' + 4y = \sin 2 t\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 0\)
4.
\(y'' + 4y = 3 \sin(13t/12)\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = -1\)
5.
\(y'' + 16y = 5\cos 4t\text{,}\) \(y(0) = 1\text{,}\) \(y'(0) = 0\)
6.
\(y'' + 16y = 5\cos(4t +1)\text{,}\) \(y(0) = 0\text{,}\) \(y'(0) = 0\)

Plotting the Envelope for Beats.

For each of the functions in Exercise Group 4.4.6.7–10 compute the envelope function. Then plot \(y(t)\) and the envelope function on the same set of axes.
7.
\(y = \cos 10 t - \cos 11t\text{.}\)
8.
\(y = \cos 11 t - \cos 12t\text{.}\)
9.
\(y = \sin 11 t - \sin 10t\text{.}\)
10.
\(y = \sin 12 t - \sin 11t\text{.}\)

11.

Consider the function \(y = \cos 10 t - \cos 11t\text{.}\)
  1. Find the mean frequency \(\overline{\omega} = (\omega_0 + \omega)/2\) and half difference \(\delta = (\omega_0 - \omega)/2\text{.}\)
  2. Use (a) to show that
    \begin{equation*} y = \cos 10 t - \cos 11t = 2 \sin\left( \frac{1}{2} t \right) \sin(\left( \frac{21}{2} t \right). \end{equation*}
    You will find the following trigonmetric identity very useful,
    \begin{equation*} 2 \sin A \sin B = \cos(A - B) - \cos(A + B). \end{equation*}
  3. Using (b), plot the graph of \(y\) and \(2 \sin(t/2)\) on the same axes.

12.

For large \(t\text{,}\) every solution of
\begin{equation*} y'' + py' + qy = \cos \omega t \end{equation*}
oscillates with angular frequency \(\omega\) and amplitude \(A\) given by
\begin{equation*} A( \omega, p, q) = \frac{1}{\sqrt{(q - \omega^2)^2 + p^2 \omega^2}}. \end{equation*}
That is, the amplitude of \(A\) is a function of the parameters \(\omega\text{,}\) \(p\text{,}\) and \(q\text{.}\)
  1. Compute \(\dfrac{\partial A}{\partial \omega}\text{.}\)
  2. For fixed \(p\) and \(q\text{,}\) let \(M(p, q)\) denote the maximum value of \(A( \omega, p, q)\) as a function of \(\omega\text{.}\) Compute an expression for \(M(p, q)\text{.}\)
  3. Set \(q = 1\) and plot \(M(p, q)\) as a function of \(p\text{.}\)
  4. Explain why \(M(p, q)\) is proportional to \(1/p\) as \(p \to 0\text{.}\)

13.

Let us examine the amplitude and phase of the steady-state solution,
\begin{equation*} x_p(t) = G(\omega) A \cos(\omega t - \phi), \end{equation*}
where
\begin{equation*} G(\omega) = \frac{1}{\sqrt{(\omega_0^2 - \omega^2)^2 + 4 c^2 \omega^2}}. \end{equation*}
Now let \(s = \omega / \omega_0\) and \(D = 2c/\omega_0\text{.}\) These new constants, \(s\) and \(D\text{,}\) measure the ratio of the driving frequency to the natural frequency and the effect of the damping force, respectively. Thus,
\begin{equation*} G = \frac{1}{\omega_0^2\sqrt{(1- s^2)^2 + D^2 s^2}} \end{equation*}
or
\begin{equation*} \omega_0^2G = \frac{1}{\sqrt{(1- s^2)^2 + D^2 s^2}}. \end{equation*}
This expression shows us how the gain varies as \(s = \omega/\omega_0\) varies. The natural frequency is fixed in
\begin{equation*} x'' + 2c x' + \omega_0^2 x = A \cos \omega t \end{equation*}
and \(D = 2c/\omega_0\) is proportional to the damping constant.
PrevTopNext