Chapter 5
Modeling with Differential Equations
5.4 Modeling With Circular Functions
-
- What is the frequency of the oscillation of a spring-mass system with mass \(m\) and spring constant \(k\)?
- What is the period?
- If the mass attached to a spring is thrown (vertically up or down) instead of let go, all that changes in our initial value problem is that the initial velocity is no longer zero.
- The initial velocity can be either positive or negative. Which sign corresponds to throwing the mass up and which down?
- Find values for the constants \(A\) and \(B\) such that the function
\(x(t)=A\,\sin\left(\sqrt{\frac{k}{m}}\,t\right)+B\,\cos\left(\sqrt{\frac{k}{m}}\,t\right)\)
satisfies the initial value problem\(\frac{d^2x}{dt^2}=-\frac{k}{m}x\), with \(x(0)=1\) and \(x'(0)=-0.5\).
- Interpret the initial value problem
\(\frac{d^2y}{dt^2}=-1.07y\), with \(y(0)=-0.5\) and \(y'(0)=0.5\)
as representing the motion of a spring-mass system. What are the mass and the spring constant? What is the initial state of the system? - A mass of 4 kilograms is suspended from a spring attached to the ceiling that has a spring constant of 9 kilograms per second squared. The mass is pulled downward 1 meter below its equilibrium position and then released from this point with an initial velocity upward of 1.5 meters per second.
- Write down the differential equation and initial conditions satisfied by the function \(x(t)\) that gives displacement of the mass from equilibrium as a function of time.
- Solve the initial value problem to find \(x(t)\).
- At what times will the mass be closest to the ceiling?
-
- Graph \(\cos(2t)\) and \(\cos^2t \) for \(0 \leq t \leq 2\pi\).
- Use the identities
\(\cos\,2t=\cos^2t-\sin^2t\) and \(\cos^2t+\sin^2t=1\)
to show that\(\cos^2t=\frac{1+\cos\,2t}{2}\)
for all \(t\). - Find all the points of inflection on the curve \(y=\cos^2t\), and note them on your graph in part (a). Why must these points also be points of inflection for the curve \(y=\sin^2t\)?
- Add the graph of \(\sin^2t\) for \(0 \leq t \leq 2\pi\) to your graphs in part (a). (You may have to restart the Graph tool.)
- Show that
\(\sin^2t=\frac{1-\cos\,2t}{2}\)
for all \(t\).
- A paper in the journal Statistics and Medicine reported on efforts to find baseline models for normal biological rhythms, such as the sleep-wake cycle, hormone levels, and core body temperature. The goal of such research is to be able to diagnose conditions such as depression by measuring disturbances in these rhythms. The principal example in the paper is based on temperature data gathered over a 5-day period from a healthy 28-year-old female. An embedded temperature probe recorded core body temperature every 2 minutes. The researchers used every tenth data point, so their time unit was 20 minutes. On that scale, they found the primary rhythm to be modeled quite well by a function of the form
\(f(t)=\mu+A\,\cos(2\pi \omega t)+B\,\sin(2\pi \omega t)\),
where \(\mu=37.5\), \(\omega=0.0138\), \(A=-0.103\), and \(B=-0.429\).- Graph the model function.
- What is the physiological significance of \(\mu\)?
- What is the physiological significance of \(\omega\)?
-
Figure P1 Approximate fit
to minutes-of-sunlight
data for Frederick, MD - Figure P2 repeats Figure 1 in this section, showing monthly averages of numbers of minutes after midnight for sunrise (blue) and sunset (red) in Frederick, MD.
Figure P2 Average monthly sunrise and sunset in Frederick, MD - Would you expect that either or both of these data curves could be modeled adequately by a function of the form \(f(t)=A+B\,\sin(C (t-D))\)? Why or why not?
- Use the Fit Sunrise Data tool to find as good a fit as you can to the sunrise data. Note the largest relative error for your fit to the sunrise data. (See Problem 7.)
- Repeat part (b) for the sunset data.
- How do your answers to parts (b) and (c) affect your answer to part (a), if at all?
-
Activity 2 was based on data from the US Naval Observatory, whose Astronomical ApplicationsWeb site provides data for sunrise/sunset, moonrise/moonset, or beginning/end of twilight for any location on the Earth. Select a location other than Frederick, MD, download one or more data sets, and explore whether your data can be modeled by a function of the form \(f(t)=A+B\,\sin(C(t-D))\). (Note: You should not expect this model to work well at extreme latitutdes. For example, in polar regions there are long periods of daylight and darkness.)
- In Activity 4, we found that all solutions of \(y''=-\omega^2y\) are of the form
\(y=A\,\sin(\omega t)+B\,\cos(\omega t)\).
- For each such function, explain why \(B\) is the value of the function at \(t=0\).
- How is \(A\) related to the value of the first derivative \(t=0\)?
- We discovered in Chapter 2 that solutions of \(y'=ky\) are all of the form \(y=Ae^{kt}\), and we discovered in this section that solutions of \(y''=-\omega^2y\) are all of the form
\(y=A\,\sin(\omega t)+B\,\cos(\omega t)\).
Now we seek solutions of \(y''=\omega^2y\).- Show that \(e^{\omega t}\) and \(e^{-\omega t}\) are both solutions of \(y''=\omega^2y\).
- Show that all functions of the form \(Ae^{\omega t}+Be^{-\omega t}\) are solutions.
- If \(\omega=2\), find a solution that satisfies the initial conditions \(y(0)=1\) and \(y'(0)=0.5\). Check your work by substituting your solution back into the differential equation and also checking the initial conditions.
- The solutions of \(y''=\omega^2y\) (see Problem 11) can be written in another form that is more suggestive of the parallel with trigonometric functions (also called circular functions). This other form uses so-called hyperbolic functions, the hyperbolic sine (abbreviated sinh) and hyperbolic cosine (abbreviated cosh). These functions may be defined from exponential functions by the following formulas:
\(\sinh(\omega t)=\frac{1}{2}\left(e^{\omega t}-e^{-\omega t}\right)\) and \(\cosh(\omega t)=\frac{1}{2}\left(e^{\omega t}+e^{-\omega t}\right)\)
- How are the first derivatives of these hyperbolic functions related to the functions themselves?
- Explain why you already know that \(\sinh(\omega t\) and \(\cosh(\omega t)\) are solutions of \(y''=\omega^2y\).
- Explain why every function of the form \(A\,\sinh(\omega t)+B\,\cosh(\omega t)\) is also a solution of the same differential equation.
- Why must \(B\) be the value of such a function at \(t=0\)?
- How is \(A\) related to the value of the first derivative at \(t=0\)?
- Compare what you did in parts (c)-(e) with what you did in Activity 4 and Exercise 15.
- If \(\omega=2\), find a solution expressed in hyperbolic functions that satisfies the initial conditions \(y(0)=1\) and \(y'(0)=0.5\). Is this the same solution you found in the preceding exercise? Why or why not?
- You showed in Activity 3 in Section 5.3 that the circular functions \(x=\cos(\omega t)\) and \(y=\sin(\omega t)\) satisfy the equation \(x^2+y^2=1\) — the equation of a circle — which is not surprising, given the origins of these functions in a circle.
- Find a similar equation satisfied by the hyperbolic functions, \(x=\cosh(\omega t)\) and \(y=\sinh(\omega t)\). (Hint: Calculate the squares of the hyperbolic sine and the hyperbolic cosine directly from their definitions in the preceding problem.)
- Why are these functions called hyperbolic?
-
- Make a sketch on paper, on a single set of coordinate axes, of the graphs of \(y=e^x\) and \(y=e^{-x}\).
- The function \(y=\cosh\,x\) is the average of the two functions whose graphs you just sketched. Sketch the average of the two graphs to find the graph of \(y=\cosh\,x\).
- Construct the graph of \(y=\sinh\,x\) in a similar manner, as the average of two exponential graphs.
- Check your work with the Graph tool (which knows the names sinh and cosh).