11.1: Applications of Sinusoids
- Page ID
- 80826
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In the same way exponential functions can be used to model a wide variety of phenomena in nature,1 the cosine and sine functions can be used to model their fair share of natural behaviors. In section 10.5, we introduced the concept of a sinusoid as a function which can be written either in the form \(C(x)=A \cos (\omega x+\phi)+B \text { for } \omega>0\) or equivalently, in the form \(S(x)=A \sin (\omega x+\phi)+B \text { for } \omega>0\). At the time, we remained undecided as to which form we preferred, but the time for such indecision is over. For clarity of exposition we focus on the sine function2 in this section and switch to the independent variable \(t\), since the applications in this section are time-dependent. We reintroduce and summarize all of the important facts and definitions about this form of the sinusoid below.
Properties of the Sinusoid \(S(t)=A \sin (\omega t+\phi)+B\)
- The amplitude is |A|
- The angular frequency is ω and the ordinary frequency is \(f=\frac{\omega}{2 \pi}\)
- The period is \(T=\frac{1}{f}=\frac{2 \pi}{\omega}\)
- The phase is φ and the phase shift is \(-\frac{\phi}{\omega}\)
- The vertical shift or baseline is B
Along with knowing these formulas, it is helpful to remember what these quantities mean in context. The amplitude measures the maximum displacement of the sine wave from its baseline (determined by the vertical shift), the period is the length of time it takes to complete one cycle of the sinusoid, the angular frequency tells how many cycles are completed over an interval of length \(2 \pi\), and the ordinary frequency measures how many cycles occur per unit of time. The phase indicates what angle \(\phi\) corresponds to \(t = 0\), and the phase shift represents how much of a ‘head start’ the sinusoid has over the un-shifted sine function. The figure below is repeated from Section 10.5.
In Section 10.1.1, we introduced the concept of circular motion and in Section 10.2.1, we developed formulas for circular motion. Our first foray into sinusoidal motion puts these notions to good use.
Recall from Exercise 55 in Section 10.1 that The Giant Wheel at Cedar Point is a circle with diameter 128 feet which sits on an 8 foot tall platform making its overall height 136 feet. It completes two revolutions in 2 minutes and 7 seconds. Assuming that the riders are at the edge of the circle, find a sinusoid which describes the height of the passengers above the ground \(t\) seconds after they pass the point on the wheel closest to the ground.
Solution
We sketch the problem situation below and assume a counter-clockwise rotation.3
We know from the equations given on page 732 in Section 10.2.1 that the y-coordinate for counterclockwise motion on a circle of radius \(r\) centered at the origin with constant angular velocity (frequency) \(\omega\) is given by \(y=r \sin (\omega t)\). Here, \(t = 0\) corresponds to the point \((r, 0)\) so that \(\theta\), the angle measuring the amount of rotation, is in standard position. In our case, the diameter of the wheel is 128 feet, so the radius is \(r\) = 64 feet. Since the wheel completes two revolutions in 2 minutes and 7 seconds (which is 127 seconds) the period \(T=\frac{1}{2}(127)=\frac{127}{2}\) seconds. Hence, the angular frequency is \(\omega=\frac{2 \pi}{T}=\frac{4 \pi}{127}\) radians per second. Putting these two pieces of information together, we have that \(y=64 \sin \left(\frac{4 \pi}{127} t\right)\) describes the \(y\)-coordinate on the Giant Wheel after \(t\) seconds, assuming it is centered at (0, 0) with \(t = 0\) corresponding to the point \(Q\). In order to find an expression for \(h\), we take the point \(O\) in the figure as the origin. Since the base of the Giant Wheel ride is 8 feet above the ground and the Giant Wheel itself has a radius of 64 feet, its center is 72 feet above the ground. To account for this vertical shift upward,4 we add 72 to our formula for \(y\) to obtain the new formula \(h=y+72=64 \sin \left(\frac{4 \pi}{127} t\right)+72\). Next, we need to adjust things so that \(t = 0\) corresponds to the point \(P\) instead of the point \(Q\). This is where the phase comes into play. Geometrically, we need to shift the angle \(\theta\) θ in the figure back \(\frac{\pi}{2}\) radians. From Section 10.2.1, we know \(\theta=\omega t=\frac{4 \pi}{127} t\), so we (temporarily) write the height in terms of \(\theta\) as \(h=64 \sin (\theta)+72\). Subtracting \(\frac{\pi}{2}\) from \(\theta\) gives the final answer \(h(t)=64 \sin \left(\theta-\frac{\pi}{2}\right)+72=64 \sin \left(\frac{4 \pi}{127} t-\frac{\pi}{2}\right)+72\). We can check the reasonableness of our answer by graphing \(y = h(t)\) over the inteval \(\left[0, \frac{127}{2}\right]\).
A few remarks about Example 11.1.1 are in order. First, note that the amplitude of 64 in our answer corresponds to the radius of the Giant Wheel. This means that passengers on the Giant Wheel never stray more than 64 feet vertically from the center of the Wheel, which makes sense. Second, the phase shift of our answer works out to be \(\frac{\pi / 2}{4 \pi / 127}=\frac{127}{8}=15.875\). This represents the ‘time delay’ (in seconds) we introduce by starting the motion at the point \(P\) as opposed to the point \(Q\). Said differently, passengers which ‘start’ at \(P\) take 15.875 seconds to ‘catch up’ to the point \(Q\).
Our next example revisits the daylight data first introduced in Section 2.5, Exercise 6b.
According to the U.S. Naval Observatory website, the number of hours \(H\) of daylight that Fairbanks, Alaska received on the 21st day of the \(n\)th month of 2009 is given below. Here \(t = 1\) represents January 21, 2009, \(t = 2\) represents February 21, 2009, and so on.
\[ \begin{array}{|l|r|r|r|r|r|r|r|r|r|r|r|r|} \hline \begin{array}{l} \text { Month } \\ \text { Number } \end{array} & 1 & 2 & 3 & 4 & 5 & 6 & 7 & 8 & 9 & 10 & 11 & 12 \\ \hline \begin{array}{l} \text { Hours of } \\ \text { Daylight } \end{array} & 5.8 & 9.3 & 12.4 & 15.9 & 19.4 & 21.8 & 19.4 & 15.6 & 12.4 & 9.1 & 5.6 & 3.3 \\ \hline \end{array}\nonumber\]
- Find a sinusoid which models these data and use a graphing utility to graph your answer along with the data.
- Compare your answer to part 1 to one obtained using the regression feature of a calculator.
Solution
- To get a feel for the data, we plot it below.
The data certainly appear sinusoidal,5 but when it comes down to it, fitting a sinusoid to data manually is not an exact science. We do our best to find the constants \(A, \omega, \phi\) and \(B\) so that the function \(H(t)=A \sin (\omega t+\phi)+B\) closely matches the data. We first go after the vertical shift \(B\) whose value determines the baseline. In a typical sinusoid, the value of \(B\) is the average of the maximum and minimum values. So here we take \(B=\frac{3.3+21.8}{2}=12.55\). Next is the amplitude \(A\) which is the displacement from the baseline to the maximum (and minimum) values. We find \(A=21.8-12.55=12.55-3.3=9.25\). At this point, we have \(H(t)=9.25 \sin (\omega t+\phi)+12.55\). Next, we go after the angular frequency \(\omega\). Since the data collected is over the span of a year (12 months), we take the period \(T=12\) months.6 This means \(\omega=\frac{2 \pi}{T}=\frac{2 \pi}{12}=\frac{\pi}{6}\). The last quantity to find is the phase \(\phi\). Unlike the previous example, it is easier in this case to find the phase shift \(-\frac{\phi}{\omega}\). Since we picked \(A>0\), the phase shift corresponds to the first value of \(t\) with \(H(t) = 12.55\) (the baseline value).7 Here, we choose \(t = 3\), since its corresponding \(H\) value of 12.4 is closer to 12.55 than the next value, 15.9, which corresponds to \(t = 4\). Hence, \(-\frac{\phi}{\omega}=3\), so \(\phi=-3 \omega=-3\left(\frac{\pi}{6}\right)=-\frac{\pi}{2}\). We have \(H(t)=9.25 \sin \left(\frac{\pi}{6} t-\frac{\pi}{2}\right)+12.55\). Below is a graph of our data with the curve \(y = H(t)\).
- Using the ‘SinReg’ command, we graph the calculator’s regression below.
While both models seem to be reasonable fits to the data, the calculator model is possibly the better fit. The calculator does not give us \(r^{2}\) value like it did for linear regressions in Section 2.5, nor does it give us \(R^{2}\) value like it did for quadratic, cubic and quartic regressions as in Section 3.1. The reason for this, much like the reason for the absence of \(R^{2}\) for the logistic model in Section 6.5, is beyond the scope of this course. We’ll just have to use our own good judgment when choosing the best sinusoid model.
11.1.1 Harmonic Motion
One of the major applications of sinusoids in Science and Engineering is the study of harmonic motion. The equations for harmonic motion can be used to describe a wide range of phenomena, from the motion of an object on a spring, to the response of an electronic circuit. In this subsection, we restrict our attention to modeling a simple spring system. Before we jump into the Mathematics, there are some Physics terms and concepts we need to discuss. In Physics, ‘mass’ is defined as a measure of an object’s resistance to straight-line motion whereas ‘weight’ is the amount of force (pull) gravity exerts on an object. An object’s mass cannot change,8 while its weight could change.
An object which weighs 6 pounds on the surface of the Earth would weigh 1 pound on the surface of the Moon, but its mass is the same in both places. In the English system of units, ‘pounds’ (lbs.) is a measure of force (weight), and the corresponding unit of mass is the ‘slug’. In the SI system, the unit of force is ‘Newtons’ (N) and the associated unit of mass is the ‘kilogram’ (kg). We convert between mass and weight using the formula9 \(w=m g\). Here, \(w\) is the weight of the object, \(m\) is the mass and \(g\) is the acceleration due to gravity. In the English system, \(g=32 \frac{\text { feet }}{\text { second}^{2}}\), and in the SI system, \(g=9.8 \frac{\text { meters }}{\text { second}^{2}}\). Hence, on Earth a \(mass\) of 1 slug \(weighs\) 32 lbs. and a \(mass\) of 1 kg \(weighs\) 9.8 N.10 Suppose we attach an object with mass \(m\) to a spring as depicted below. The weight of the object will stretch the spring. The system is said to be in ‘equilibrium’ when the weight of the object is perfectly balanced with the restorative force of the spring. How far the spring stretches to reach equilibrium depends on the spring’s ‘spring constant’. Usually denoted by the letter \(k\), the spring constant relates the force \(F\) applied to the spring to the amount \(d\) the spring stretches in accordance with Hooke’s Law11 \(F=k d\). If the object is released above or below the equilibrium position, or if the object is released with an upward or downward velocity, the object will bounce up and down on the end of the spring until some external force stops it. If we let \(x(t)\) denote the object’s displacement from the equilibrium position at time \(t\), then \(x(t) = 0\) means the object is at the equilibrium position, \(x(t)<0\) means the object is above the equilibrium position, and \(x(t)>0\) means the object is below the equilibrium position. The function \(x(t)\) is called the ‘equation of motion’ of the object.12
If we ignore all other influences on the system except gravity and the spring force, then Physics tells us that gravity and the spring force will battle each other forever and the object will oscillate indefinitely. In this case, we describe the motion as ‘free’ (meaning there is no external force causing the motion) and ‘undamped’ (meaning we ignore friction caused by surrounding medium, which in our case is air). The following theorem, which comes from Differential Equations, gives \(x(t)\) as a function of the mass \(m\) of the object, the spring constant \(k\), the initial displacement \(x_{0}\) of the object and initial velocity \(v_{0}\) of the object. As with \(x(t)\), \(x_{0}=0\) means the object is released from the equilibrium position, \(x_{0}<0\) means the object is released above the equilibrium position and \(x_{0}>0\) means the object is released below the equilibrium position. As far as the initial velocity \(v_{0}\) is concerned, \(v_{0}=0\) means the object is released ‘from rest,’ \(v_{0}<0\) means the object is heading upwards and \(v_{0}>0\) means the object is heading downwards.13
Suppose an object of mass \(m\) is suspended from a spring with spring constant \(k\). If the initial displacement from the equilibrium position is \(x_{0}\) and the initial velocity of the object is \(v_{0}\), then the displacement \(x\) from the equilibrium position at time \(t\) is given by \(x(t)=A \sin (\omega t+\phi)\) where
- \(\omega=\sqrt{\frac{k}{m}} \text { and } A=\sqrt{x_{0}^{2}+\left(\frac{v_{0}}{\omega}\right)^{2}}\)
- \(A \sin (\phi)=x_{0} \text { and } A \omega \cos (\phi)=v_{0}\).
It is a great exercise in ‘dimensional analysis’ to verify that the formulas given in Theorem 11.1 work out so that \(\omega\) has units \(\frac{1}{s}\) and \(A\) has units ft. or m, depending on which system we choose.
Suppose an object weighing 64 pounds stretches a spring 8 feet.
- If the object is attached to the spring and released 3 feet below the equilibrium position from rest, find the equation of motion of the object, \(x(t)\). When does the object first pass through the equilibrium position? Is the object heading upwards or downwards at this instant?
- If the object is attached to the spring and released 3 feet below the equilibrium position with an upward velocity of 8 feet per second, find the equation of motion of the object, \(x(t)\). What is the longest distance the object travels above the equilibrium position? When does this first happen? Confirm your result using a graphing utility.
Solution
In order to use the formulas in Theorem 11.1, we first need to determine the spring constant \(k\) and the mass of the object \(m\). To find \(k\), we use Hooke’s Law \(F = kd\). We know the object weighs 64 lbs. and stretches the spring 8 ft.. Using \(F = 64\) and \(d = 8\), we get \(64=k \cdot 8\), or \(k=8 \frac{\mathrm{lbs.}}{\mathrm{ft.}}\). To find \(m\), we use \(w = mg\) with \(w = 64\) lbs. and \(g=32 \frac{\mathrm{ft.}}{s^{2}}\). We get \(m = 2\) slugs. We can now proceed to apply Theorem 11.1.
- With \(k = 8\) and \(m = 2\), we get \(\omega=\sqrt{\frac{k}{m}}=\sqrt{\frac{8}{2}}=2\). We are told that the object is released 3 feet below the equilibrium position ‘from rest.’ This means \(x_{0}=3\) and \(v_{0}=0\). Therefore, \(A=\sqrt{x_{0}^{2}+\left(\frac{v_{0}}{\omega}\right)^{2}}=\sqrt{3^{2}+0^{2}}=3\). To determine the phase \(\phi\), we have \(A \sin (\phi)=x_{0}\), which in this case gives \(3 \sin (\phi)=3\) so \(\sin (\phi)=1\). Only \(\phi=\frac{\pi}{2}\) and angles coterminal to it satisfy this condition, so we pick14 the phase to be \(\phi=\frac{\pi}{2}\). Hence, the equation of motion is \(x(t)=3 \sin \left(2 t+\frac{\pi}{2}\right)\). To find when the object passes through the equilibrium position we solve \(x(t)=3 \sin \left(2 t+\frac{\pi}{2}\right)=0\). Going through the usual analysis we find \(t=-\frac{\pi}{4}+\frac{\pi}{2} k\) for integers \(k\). Since we are interested in the first time the object passes through the equilibrium position, we look for the smallest positive \(t\) value which in this case is \(t=\frac{\pi}{4} \approx 0.78\) 8 seconds after the start of the motion. Common sense suggests that if we release the object below the equilibrium position, the object should be traveling upwards when it first passes through it. To check this answer, we graph one cycle of \(x(t)\). Since our applied domain in this situation is \(t \geq 0\), and the period of \(x(t)\) is \(T=\frac{2 \pi}{\omega}=\frac{2 \pi}{2}=\pi\), we graph \(x(t)\) over the interval \([0, \pi]\). Remembering that \(x(t)>0\) means the object is below the equilibrium position and \(x(t)<0\) means the object is above the equilibrium position, the fact our graph is crossing through the \(t\)-axis from positive x to negative \(x\) at \(t=\frac{\pi}{4}\) confirms our answer.
- The only difference between this problem and the previous problem is that we now release the object with an upward velocity of \(8 \frac{\mathrm{ft}}{s}\). We still have \(\omega=2\) and \(x_{0}=3\), but now we have \(v_{0}=-8\), the negative indicating the velocity is directed upwards. Here, we get \(A=\sqrt{x_{0}^{2}+\left(\frac{v_{0}}{\omega}\right)^{2}}=\sqrt{3^{2}+(-4)^{2}}=5\). 5. From \(A \sin (\phi)=x_{0}\), we get \(5 \sin (\phi)=3\) which gives \(\sin (\phi)=\frac{3}{5}\). From \(A \omega \cos (\phi)=v_{0}\), we get \(10 \cos (\phi)=-8\), or \(\cos (\phi)=-\frac{4}{5}\). This means that \(\phi\) is a Quadrant II angle which we can describe in terms of either arcsine or arccosine. Since \(x(t)\) is expressed in terms of sine, we choose to express \(\phi=\pi-\arcsin \left(\frac{3}{5}\right)\). Hence, \(x(t)=5 \sin \left(2 t+\left[\pi-\arcsin \left(\frac{3}{5}\right)\right]\right)\).Since the amplitude of \(x(t)\) is 5, the object will travel at most 5 feet above the equilibrium position. To find when this happens, we solve the equation \(x(t)=5 \sin \left(2 t+\left[\pi-\arcsin \left(\frac{3}{5}\right)\right]\right)=-5\), the negative once again signifying that the object is above the equilibrium position. Going through the usual machinations, we get \(t=\frac{1}{2} \arcsin \left(\frac{3}{5}\right)+\frac{\pi}{4}+\pi k\) for integers \(k\). The smallest of these values occurs when \(k = 0\), that is, \(t=\frac{1}{2} \arcsin \left(\frac{3}{5}\right)+\frac{\pi}{4} \approx 1.107\) seconds after the start of the motion. To check our answer using the calculator, we graph \(y=5 \sin \left(2 x+\left[\pi-\arcsin \left(\frac{3}{5}\right)\right]\right)\) on a graphing utility and confirm the coordinates of the first relative minimum to be approximately (1.107,-5).
It is possible, though beyond the scope of this course, to model the effects of friction and other external forces acting on the system.15 While we may not have the Physics and Calculus background to derive equations of motion for these scenarios, we can certainly analyze them. We examine three cases in the following example.
- Write \(x(t)=5 e^{-t / 5} \cos (t)+5 e^{-t / 5} \sqrt{3} \sin (t)\) in the form \(x(t)=A(t) \sin (\omega t+\phi)\). Graph \(x(t)\) using a graphing utility.
- Write \(x(t)=(t+3) \sqrt{2} \cos (2 t)+(t+3) \sqrt{2} \sin (2 t)\) in the form \(x(t)=A(t) \sin (\omega t+\phi)\). Graph \(x(t)\) using a graphing utility.
- Find the period of \(x(t) = 5 sin(6t) − 5 sin (8t)\). Graph \(x(t)\) using a graphing utility.
Solution
- We start rewriting \(x(t)=5 e^{-t / 5} \cos (t)+5 e^{-t / 5} \sqrt{3} \sin (t)\) by factoring out \(5 e^{-t / 5}\) from both terms to get \(x(t)=5 e^{-t / 5}(\cos (t)+\sqrt{3} \sin (t))\). We convert what’s left in parentheses to the required form using the formulas introduced in Exercise 36 from Section 10.5. We find \((\cos (t)+\sqrt{3} \sin (t))=2 \sin \left(t+\frac{\pi}{3}\right)\) so that \(x(t)=10 e^{-t / 5} \sin \left(t+\frac{\pi}{3}\right)\). Graphing this on the calculator as \(y=10 e^{-x / 5} \sin \left(x+\frac{\pi}{3}\right)\) reveals some interesting behavior. The sinusoidal nature continues indefinitely, but it is being attenuated. In the sinusoid \(A \sin (\omega x+\phi)\), the coefficient \(A\) of the sine function is the amplitude. In the case of \(y=10 e^{-x / 5} \sin \left(x+\frac{\pi}{3}\right)\), we can think of the function \(A(x)=10 e^{-x / 5}\) as the amplitude. As \(x \rightarrow \infty, 10 e^{-x / 5} \rightarrow 0\) which means the amplitude continues to shrink towards zero. Indeed, if we graph \(y=\pm 10 e^{-x / 5}\) along with \(y=10 e^{-x / 5} \sin \left(x+\frac{\pi}{3}\right)\), we see this attenuation taking place. This equation corresponds to the motion of an object on a spring where there is a slight force which acts to ‘damp’, or slow the motion. An example of this kind of force would be the friction of the object against the air. In this model, the object oscillates forever, but with smaller and smaller amplitude.
- Proceeding as in the first example, we factor out \((t+3) \sqrt{2}\) from each term in the function \(x(t)=(t+3) \sqrt{2} \cos (2 t)+(t+3) \sqrt{2} \sin (2 t)\) to get \(x(t)=(t+3) \sqrt{2}(\cos (2 t)+\sin (2 t))\). We find \((\cos (2 t)+\sin (2 t))=\sqrt{2} \sin \left(2 t+\frac{\pi}{4}\right)\), so \(x(t)=2(t+3) \sin \left(2 t+\frac{\pi}{4}\right)\). Graphing this on the caculation as \(y=2(x+3) \sin \left(2 x+\frac{\pi}{4}\right)\), we find the sinusoid’s amplitude growing. Since our amplitude function here is \(A(x)=2(x+3)=2 x+6\), which continues to grow without bound as \(x \rightarrow \infty\), this is hardly surprising. The phenomenon illustrated here is ‘forced’ motion. That is, we imagine that the entire apparatus on which the spring is attached is oscillating as well. In this case, we are witnessing a ‘resonance’ effect – the frequency of the external oscillation matches the frequency of the motion of the object on the spring.16
- Last, but not least, we come to \(x(t) = 5 sin(6t)−5 sin(8t)\). To find the period of this function, we need to determine the length of the smallest interval on which both \(f(t) = 5 sin(6t)\) and \(g(t) = 5 sin(8t)\) complete a whole number of cycles. To do this, we take the ratio of their frequencies and reduce to lowest terms: \(\frac{6}{8}=\frac{3}{4}\). This tells us that for every 3 cycles \(f\) makes, \(g\) makes 4. In other words, the period of \(x(t)\) is three times the period of \(f(t)\) (which is four times the period of \(g(t)\)), or \(\pi\). We graph \(y = 5 sin(6x) − 5 sin(8x)\) over \([0, \pi]\) on the calculator to check this. This equation of motion also results from ‘forced’ motion, but here the frequency of the external oscillation is different than that of the object on the spring. Since the sinusoids here have different frequencies, they are ‘out of sync’ and do not amplify each other as in the previous example. Taking things a step further, we can use a sum to product identity to rewrite \(x(t) = 5 sin(6t) − 5 sin(8t) as x(t) = −10 sin(t) cos(7t)\). The lower frequency factor in this expression, \(−10 sin(t)\), plays an interesting role in the graph of \(x(t)\). Below we graph \(y = 5 sin(6x) − 5 sin(8x)\) and \(y=\pm 10 \sin (x)\) over \([0,2 \pi]\). This is an example of the ‘beat’ phenomena, and the curious reader is invited to explore this concept as well.17
11.1.2 Exercises
- The sounds we hear are made up of mechanical waves. The note ‘A’ above the note ‘middle C’ is a sound wave with ordinary frequency \(f=440 \text { Hertz }=440 \frac{\text { cycles }}{\text { second }}\). Find a sinusoid which models this note, assuming that the amplitude is 1 and the phase shift is 0.
- The voltage \(V\) in an alternating current source has amplitude \(220 \sqrt{2}\) and ordinary frequency \(f=60 \text { Hertz }\). Find a sinusoid which models this voltage. Assume that the phase is 0.
- The London Eye is a popular tourist attraction in London, England and is one of the largest Ferris Wheels in the world. It has a diameter of 135 meters and makes one revolution (counterclockwise) every 30 minutes. It is constructed so that the lowest part of the Eye reaches ground level, enabling passengers to simply walk on to, and off of, the ride. Find a sinsuoid which models the height \(h\) of the passenger above the ground in meters \(t\) minutes after they board the Eye at ground level.
- On page 732 in Section 10.2.1, we found the \(x\)-coordinate of counter-clockwise motion on a circle of radius \(r\) with angular frequency \(\omega\) to be \(x=r \cos (\omega t)\), where \(t = 0\) corresponds to the point \((r, 0)\). Suppose we are in the situation of Exercise 3 above. Find a sinsusoid which models the horizontal displacement \(x\) of the passenger from the center of the Eye in meters \(t\) minutes after they board the Eye. Here we take \(x(t)>0\) to mean the passenger is to the right of the center, while \(x(t)<0\) means the passenger is to the left of the center.
- In Exercise 52 in Section 10.1, we introduced the yo-yo trick ‘Around the World’ in which a yo-yo is thrown so it sweeps out a vertical circle. As in that exercise, suppose the yo-yo string is 28 inches and it completes one revolution in 3 seconds. If the closest the yo-yo ever gets to the ground is 2 inches, find a sinsuoid which models the height \(h\) of the yo-yo above the ground in inches \(t\) seconds after it leaves its lowest point.
- Suppose an object weighing 10 pounds is suspended from the ceiling by a spring which stretches 2 feet to its equilibrium position when the object is attached.
- Find the spring constant \(k\) in \(\frac{\mathrm{lbs} .}{\mathrm{ft} .}\) and the mass of the object in slugs.
- Find the equation of motion of the object if it is released from 1 foot below the equilibrium position from rest. When is the first time the object passes through the equilibrium position? In which direction is it heading?
- Find the equation of motion of the object if it is released from 6 inches above the equilibrium position with a downward velocity of 2 feet per second. Find when the object passes through the equilibrium position heading downwards for the third time.
-
\[T=2 \pi \sqrt{\frac{l}{g}}\nonumber\]
where \(l\) is the length of the pendulum and \(g\) is the acceleration due to gravity.
- Find a sinusoid which gives the angular displacement \(\theta\) as a function of time, \(t\). Arrange things so \(\theta(0)=\theta_{0}\).
- In Exercise 40 section 5.3, you found the length of the pendulum needed in Jeff’s antique Seth-Thomas clock to ensure the period of the pendulum is \(\frac{1}{2}\) of a second. Assuming the initial displacement of the pendulum is \(15^{\circ}\), find a sinusoid which models the displacement of the pendulum \(\theta\) as a function of time, \(t\), in seconds.
- The table below lists the average temperature of Lake Erie as measured in Cleveland, Ohio on the first of the month for each month during the years 1971 – 2000.19 For example, \(t = 3\) represents the average of the temperatures recorded for Lake Erie on every March 1 for the years 1971 through 2000.
\[\begin{array}{|l|r|r|r|r|r|r|r|r|r|r|r|r|} \hline \begin{array}{l} \text { Month } \\ \text { Number, } t \end{array} & 1 & 2 & 3 & 4 & 5 & 6 & 7 & 8 & 9 & 10 & 11 & 12 \\
\hline \begin{array}{l} \text { Temperature } \\ \left({ }^{\circ} \mathrm{F}\right), T \end{array} & 36 & 33 & 34 & 38 & 47 & 57 & 67 & 74 & 73 & 67 & 56 & 46 \\ \hline
\end{array}\nonumber\]- Using the techniques discussed in Example 11.1.2, fit a sinusoid to these data.
- Using a graphing utility, graph your model along with the data set to judge the reasonableness of the fit.
- Use the model you found in part 8a to predict the average temperature recorded for Lake Erie on April 15th and September 15th during the years 1971–2000.20
- Compare your results to those obtained using a graphing utility.
- The fraction of the moon illuminated at midnight Eastern Standard Time on the \(t^{\text {th }}\) day of June, 2009 is given in the table below.21
- Using the techniques discussed in Example 11.1.2, fit a sinusoid to these data.22
- Using a graphing utility, graph your model along with the data set to judge the reasonableness of the fit.
- Use the model you found in part 9a to predict the fraction of the moon illuminated on June 1, 2009. 23
- Compare your results to those obtained using a graphing utility.
- With the help of your classmates, research the phenomena mentioned in Example 11.1.4, namely resonance and beats.
- With the help of your classmates, research Amplitude Modulation and Frequency Modulation.
- What other things in the world might be roughly sinusoidal? Look to see what models you can find for them and share your results with your class.
11.1.3 Answers
- \(S(t)=\sin (880 \pi t)\)
- \(V(t)=220 \sqrt{2} \sin (120 \pi t)\)
- \(h(t)=28 \sin \left(\frac{2 \pi}{3} t-\frac{\pi}{2}\right)+30\)
- \(x(t)=67.5 \cos \left(\frac{\pi}{15} t-\frac{\pi}{2}\right)=67.5 \sin \left(\frac{\pi}{15} t\right)\)
- \(h(t)=28 \sin \left(\frac{2 \pi}{3} t-\frac{\pi}{2}\right)+30\)
-
- \(k=5 \frac{\text { lbs. }}{\text { ft. }} \text { and } m=\frac{5}{16} \text { slugs }\)
- \(x(t)=\sin \left(4 t+\frac{\pi}{2}\right)\). The object first passes through the equilibrium point when \(t=\frac{\pi}{8} \approx 0.39\) 0.39 seconds after the motion starts. At this time, the object is heading upwards.
- \(x(t)=\frac{\sqrt{2}}{2} \sin \left(4 t+\frac{7 \pi}{4}\right)\). The object passes through the equilibrium point heading downwards for the third time when \(t=\frac{17 \pi}{16} \approx 3.34\) seconds.
-
- \(\theta(t)=\theta_{0} \sin \left(\sqrt{\frac{g}{l}} t+\frac{\pi}{2}\right)\)
- \(\theta(t)=\frac{\pi}{12} \sin \left(4 \pi t+\frac{\pi}{2}\right)\)
-
- \(T(t)=20.5 \sin \left(\frac{\pi}{6} t-\pi\right)+53.5\)
- Our function and the data set are graphed below. The sinusoid seems to be shifted to the right of our data.
- The average temperature April \(15^{\text {th }}\) is approximately \(T(4.5) \approx 39.00^{\circ} \mathrm{F}\) and the average temperature on September \(15^{\text {th }}\) is approximately \(T(9.5) \approx 73.38^{\circ} \mathrm{F}\).
- Using a graphing calculator, we get the following
This model predicts the average temperature for April \(15^{\mathrm{th}}\) to be approximately \(42.43^{\circ} \mathrm{F}\) and the average temperature on September \(15^{\text {th }}\) to be approximately \(70.05^{\circ} \mathrm{F}\). This model appears to be more accurate.
-
- Based on the shape of the data, we either choose \(A<0\) or we find the second value of \(t\) which closely approximates the ‘baseline’ value, \(F=0.505\). We choose the latter to obtain \(F(t)=0.475 \sin \left(\frac{\pi}{15} t-2 \pi\right)+0.505=0.475 \sin \left(\frac{\pi}{15} t\right)+0.505\).
- Our function and the data set are graphed below. It’s a pretty good fit.
- The fraction of the moon illuminated on June 1st, 2009 is approximately \(F(1) \approx 0.60\)
- Using a graphing calculator, we get the following.
This model predicts that the fraction of the moon illuminated on June 1st, 2009 is approximately 0.59. This appears to be a better fit to the data than our first model.
Reference
1 See Section 6.5.
2 Sine haters can use the co-function identity \(\cos \left(\frac{\pi}{2}-\theta\right)=\sin (\theta)\) to turn all of the sines into cosines.
3 Otherwise, we could just observe the motion of the wheel from the other side.
4 We are readjusting our ‘baseline’ from \(y = 0\) to \(y =72\).
5 Okay, it appears to be the '\(\wedge\)' shape we saw in some of the graphs in Section 2.2. Just humor us.
6 Even though the data collected lies in the interval [1, 12], which has a length of 11, we need to think of the data point at t = 1 as a representative sample of the amount of daylight for every day in January. That is, it represents \(H(t)\) over the interval [0, 1]. Similarly, \(t = 2\) is a sample of \(H(t)\) over [1, 2], and so forth.
7 See the figure on page 882.
8 Well, assuming the object isn’t subjected to relativistic speeds . . .
9 This is a consequence of Newton’s Second Law of Motion \(F=m a\) where \(F\) is force, \(m\) is mass and \(a\) is acceleration. In our present setting, the force involved is weight which is caused by the acceleration due to gravity.
10 Note that \(1 \text { pound }=1 \frac{\text { slug foot }}{\text { second }^2}\) and \(1 \text { Newton }=1 \frac{\mathrm{kg} \text { meter }}{\text { second }^2}\).
11 Look familiar? We saw Hooke’s Law in Section 4.3.1.
12 To keep units compatible, if we are using the English system, we use feet (ft.) to measure displacement. If we are in the SI system, we measure displacement in meters (m). Time is always measured in seconds (s).
13 The sign conventions here are carried over from Physics. If not for the spring, the object would fall towards the ground, which is the ‘natural’ or ‘positive’ direction. Since the spring force acts in direct opposition to gravity, any movement upwards is considered ‘negative’.
14 For confirmation, we not that \(A \omega \cos (\phi)=v_{0}\), which in this case reduces to \(6 \cos (\phi)=0\).
15 Take a good Differential Equations class to see this!
16 The reader is invited to investigate the destructive implications of resonance.
17 A good place to start is this article on beats.
18 Provided \(\theta\) is kept ‘small.’ Carl remembers the ‘Rule of Thumb’ as being \(20^{\circ}\) or less. Check with your friendly neighborhood physicist to make sure.
19 See this website: http://www.erh.noaa.gov/cle/climate/cle/normals/laketempcle.html.
20 The computed average is \(41^{\circ} \mathrm{F}\) for April \(15^{\text {th }}\) and \(71^{\circ} \mathrm{F}\) for September \(15^{\text {th }}\).
21 See this website: http://www.usno.navy.mil/USNO/astronomical-applications/data-services/frac-moon-ill.
22 You may want to plot the data before you find the phase shift.
23 The listed fraction is 0.62.