Exercise 6.2.1.
Verify Table 6.2.
Section 6.2 Transforms of derivatives and ODEs
Subsection 6.2.1 Transforms of derivatives
Let us see how the Laplace transform is used for differential equations. First we find the Laplace transform of the derivative of a function. Suppose \(g(t)\) is a differentiable function of exponential order, that is, \(\lvert g(t) \rvert \leq M e^{ct}\) for some \(M\) and \(c\text{.}\) So \(\mathcal{L} \bigl\{ g(t) \bigr\}\) exists, and what is more, \(\lim_{t\to\infty} e^{-st}g(t) = 0\) when \(s > c\text{.}\) Then
\begin{equation*}
\mathcal{L} \bigl\{ g'(t) \bigr\}
=
\int_0^\infty
e^{-st}
g'(t) \,dt
=
\Bigl[e^{-st} g(t) \Bigr]_{t=0}^\infty
-
\int_0^\infty
(-s)\,
e^{-st}
g(t) \,dt
=
-g(0) + s \mathcal{L} \bigl\{ g(t) \bigr\} .
\end{equation*}
We repeat this procedure for higher derivatives. The results are listed in Table 6.2. The procedure also works for continuous piecewise smooth functions, that is, functions that are continuous with a piecewise continuous derivative.
| \(f(t)\) | \(\mathcal{L} \bigl\{ f(t) \bigr\} = F(s)\) |
| \(g'(t)\) | \(sG(s)-g(0)\) |
| \(g''(t)\) | \(s^2G(s)-sg(0)-g'(0)\) |
| \(g'''(t)\) | \(s^3G(s)-s^2g(0)-sg'(0)-g''(0)\) |
Subsection 6.2.2 Solving ODEs with the Laplace transform
Notice that the Laplace transform turns differentiation into multiplication by \(s\text{.}\) It is this property that makes it useful to apply the transform to differential equations.
Example 6.2.1.
Take the equation
\begin{equation*}
x''(t) + x(t) = \cos (2t), \quad x(0) = 0, \quad x'(0) = 1 .
\end{equation*}
We will take the Laplace transform of both sides. By \(X(s)\) we will, as usual, denote the Laplace transform of \(x(t)\text{.}\)
\begin{equation*}
\begin{aligned}
\mathcal{L} \bigl\{ x''(t) + x(t) \bigr\} & = \mathcal{L} \bigl\{ \cos (2t) \bigr\} , \\
s^2 X(s) -sx(0)-x'(0) + X(s) & = \frac{s}{s^2 + 4} .
\end{aligned}
\end{equation*}
We plug in the initial conditions now—making computations more streamlined—to find
\begin{equation*}
s^2 X(s) -1 + X(s) = \frac{s}{s^2 + 4} .
\end{equation*}
We solve for \(X(s)\text{,}\)
\begin{equation*}
X(s) = \frac{s}{(s^2+1)(s^2 + 4)} + \frac{1}{s^2+1} .
\end{equation*}
We use partial fractions (exercise) to write
\begin{equation*}
X(s) =\frac{1}{3} \, \frac{s}{s^2+1} -
\frac{1}{3}\, \frac{s}{s^2+4} + \frac{1}{s^2+1} .
\end{equation*}
Now take the inverse Laplace transform to obtain
\begin{equation*}
x(t) =\frac{1}{3} \cos (t) -
\frac{1}{3} \cos (2t) + \sin (t) .
\end{equation*}
The procedure for linear constant coefficient equations is as follows: Take an ordinary differential equation in the time variable \(t\text{.}\) Apply the Laplace transform to transform the equation into an algebraic (non differential) equation in the frequency domain. All the \(x(t)\text{,}\) \(x'(t)\text{,}\) \(x''(t)\text{,}\) and so on, will be converted to \(X(s)\text{,}\) \(sX(s) - x(0)\text{,}\) \(s^2X(s) - sx(0) - x'(0)\text{,}\) and so on. Solve the equation for \(X(s)\text{.}\) Then taking the inverse transform, if possible, find \(x(t)\text{.}\)
It should be noted that since not every function has a Laplace transform, not every equation can be solved in this manner. Also if the equation is not a linear constant coefficient ODE, then by applying the Laplace transform we may not obtain an algebraic equation.
Subsection 6.2.3 Using the Heaviside function
Before we move on to more general equations than those we could solve before, we want to consider the Heaviside function. See Figure 6.1 for the graph.
\begin{equation*}
u(t) =
\begin{cases}
0 & \text{if } \; t < 0 , \\
1 & \text{if } \; t \geq 0 .
\end{cases}
\end{equation*}
The Heaviside function is useful for putting other functions together or cutting functions off. Usually we use \(u(t-a)\) for some constant \(a\text{;}\) we just shift the graph to the right by \(a\text{.}\) That is, \(u(t-a) = 0\) when \(t < a\) and \(u(t-a) = 1\) when \(t \geq a\text{.}\) For example, suppose \(f(t)\) is a “signal”, and we started receiving \(\sin t\) at time \(t=\pi\text{.}\) The function \(f(t)\) is then defined as
\begin{equation*}
f(t) =
\begin{cases}
0 & \text{if } \; t < \pi , \\
\sin t & \text{if } \; t \geq \pi .
\end{cases}
\end{equation*}
Using the Heaviside function, \(f(t)\) can be written as
\begin{equation*}
f(t) = u(t - \pi) \, \sin t .
\end{equation*}
Similarly, the step function that is \(1\) on the interval \([1,2)\) and \(0\) elsewhere is written as
\begin{equation*}
u(t - 1) - u(t-2) .
\end{equation*}
With the Heaviside function we can express functions defined piecewise. If \(f(t) = t\) when \(t\) is in \([0,1]\text{,}\) \(f(t) = -t+2\) when \(t\) is in \([1,2]\text{,}\) and \(f(t) = 0\) otherwise, then you write
\begin{equation*}
f(t) = t \, \bigl( u(t) - u(t-1) \bigr) +
(-t+2) \, \bigl( u(t - 1) - u(t-2) \bigr) .
\end{equation*}
How does the Heaviside function interact with the Laplace transform? We saw that
\begin{equation*}
\mathcal{L} \bigl\{ u(t-a) \bigr\} = \frac{e^{-as}}{s} .
\end{equation*}
This computation can be generalized into a shifting property or second shifting property.
\begin{equation}
\mybxbg{~~
\mathcal{L} \bigl\{ f(t-a) \, u(t-a) \bigr\}
= e^{-as} \mathcal{L} \bigl\{ f(t) \bigr\} .
~~}\tag{6.1}
\end{equation}
For example,
\begin{equation*}
\mathcal{L} \bigl\{ t \, u(t-1) \bigr\}
=
\mathcal{L} \bigl\{ \bigl( (t-1) + 1\bigr) \, u(t-1) \bigr\}
= e^{-as} \mathcal{L} \bigl\{ t+1 \bigr\}
= e^{-as} \left(\frac{1}{s^2} + \frac{1}{s} \right)
.
\end{equation*}
Example 6.2.2.
Consider the mass-spring system
\begin{equation*}
x''(t) + x(t) = f(t) , \quad x(0) = 0, \quad x'(0) = 0,
\end{equation*}
where \(f(t) = 1\) if \(1 \leq t < 5\) and zero otherwise. The \(f(t)\) is not periodic and defined piecewise; no problem for Laplace. Imagine a rocket attached to the mass is fired for 4 seconds starting at \(t=1\text{.}\) Or perhaps imagine an RLC circuit, where the voltage is raised at a constant rate for 4 seconds starting at \(t=1\) and then held steady again starting at \(t=5\text{.}\)
We write \(f(t) = u(t-1) - u(t-5)\text{.}\) We transform the equation and we plug in the initial conditions as before
\begin{equation*}
s^2 X(s) + X(s) = \frac{e^{-s}}{s} - \frac{e^{-5s}}{s} .
\end{equation*}
We solve for \(X(s)\text{,}\)
\begin{equation*}
X(s) = \frac{e^{-s}}{s(s^2+1)} - \frac{e^{-5s}}{s(s^2+1)} .
\end{equation*}
We leave it as an exercise to the reader to show that
\begin{equation*}
{\mathcal{L}}^{-1} \left\{ \frac{1}{s(s^2+1)} \right\}
= 1 - \cos t .
\end{equation*}
In other words, \(\mathcal{L} \{ 1 - \cos t \} =
\frac{1}{s(s^2+1)}\text{.}\) Using the second shifting property (6.1), we find
\begin{equation*}
{\mathcal{L}}^{-1} \left\{ \frac{e^{-s}}{s(s^2+1)} \right\}
=
{\mathcal{L}}^{-1} \bigl\{
e^{-s}
\mathcal{L} \{ 1 - \cos t \}
\bigr\}
=
\bigl( 1 - \cos (t-1) \bigr) \, u(t-1) .
\end{equation*}
Similarly,
\begin{equation*}
{\mathcal{L}}^{-1} \left\{ \frac{e^{-5s}}{s(s^2+1)} \right\}
=
{\mathcal{L}}^{-1} \bigl\{
e^{-5s}
\mathcal{L} \{ 1 - \cos t \}
\bigr\}
=
\bigl( 1 - \cos (t-5) \bigr) \, u(t-5) .
\end{equation*}
Hence, the solution is
\begin{equation*}
x(t) =
\bigl( 1 - \cos (t-1) \bigr) \, u(t-1) -
\bigl( 1 - \cos (t-5) \bigr) \, u(t-5) .
\end{equation*}
The plot of this solution is given in Figure 6.2.
Subsection 6.2.4 Transfer functions
The Laplace transform leads to the following useful concept for studying the steady state behavior of a linear system. Consider an equation of the form
\begin{equation*}
L x = f(t) ,
\end{equation*}
where \(L\) is a linear constant coefficient differential operator. Then \(f(t)\) is usually thought of as input of the system and \(x(t)\) is thought of as the output of the system. For example, for a mass-spring system the input is the forcing function and the output is the behavior of the mass. We would like to have a convenient way to study the behavior of the system for different inputs.
Let us suppose that all the initial conditions are zero. We take the Laplace transform of the equation to obtain the equation
\begin{equation*}
A(s) X(s) = F(s) .
\end{equation*}
Solving for the ratio \(\nicefrac{X(s)}{F(s)}\text{,}\) we obtain the so-called transfer function \(H(s) = \nicefrac{1}{A(s)}\text{,}\) that is,
\begin{equation*}
H(s) = \frac{X(s)}{F(s)} .
\end{equation*}
In other words, \(X(s) = H(s) F(s)\text{.}\) We obtain an algebraic dependence of the output of the system based on the input. We can now easily study the steady state behavior of the system given different inputs by simply multiplying by the transfer function. Moreover, it is possible to compute the \(H(s)\) without knowing exactly what the equation is by observing the output \(X(s)\) for a given input \(F(s)\text{.}\) Once \(H(s)\) is known, you can find the output for any input.
Example 6.2.3.
Given \(x'' + \omega_0^2 x = f(t)\) (assume the initial conditions are zero), let us find the transfer function.
First, we take the Laplace transform of the equation,
\begin{equation*}
s^2 X(s) + \omega_0^2 X(s) = F(s) .
\end{equation*}
Now we solve for the transfer function \(\nicefrac{X(s)}{F(s)}\text{,}\)
\begin{equation*}
H(s) = \frac{X(s)}{F(s)} = \frac{1}{s^2 + \omega_0^2} .
\end{equation*}
Let us see how to use the transfer function. Suppose we have the constant input \(f(t) = 1\text{.}\) Hence \(F(s) = \nicefrac{1}{s}\text{,}\) and
\begin{equation*}
X(s) = H(s) F(s) = \frac{1}{s^2+\omega_0^2} \, \frac{1}{s} .
\end{equation*}
Taking the inverse Laplace transform of \(X(s)\text{,}\) we obtain
\begin{equation*}
x(t) = \frac{1-\cos(\omega_0 t)}{\omega_0^2} .
\end{equation*}
Similarly, for any other input \(F(s)\text{,}\) the output is \(X(s) = H(s) F(s) = \frac{1}{s^2+\omega_0^2} \, F(s)\text{.}\)
Subsection 6.2.5 Transforms of integrals
A feature of Laplace transforms is that it is also able to easily deal with integral equations. That is, equations in which integrals rather than derivatives of functions appear. The basic property, which can be proved by applying the definition and doing integration by parts, is
\begin{equation*}
\mybxbg{~~
\mathcal{L} \left\{
\int_0^t f(\tau) \, d\tau
\right\} = \frac{1}{s} \, F(s) .
~~}
\end{equation*}
It is sometimes useful (e.g. for computing the inverse transform) to write this as
\begin{equation*}
\int_0^t f(\tau) \, d\tau
=
{\mathcal{L}}^{-1} \left\{
\frac{1}{s} \, F(s) \right\} .
\end{equation*}
Example 6.2.4.
To compute \({\mathcal{L}}^{-1} \left\{\frac{1}{s(s^2+1)}\right\}\) we could apply this integration rule:
\begin{equation*}
{\mathcal{L}}^{-1} \left\{
\frac{1}{s} \, \frac{1}{s^2+1} \right\}
=
\int_0^t
{\mathcal{L}}^{-1} \left\{
\frac{1}{s^2+1} \right\} \, d\tau
=
\int_0^t
\sin \tau \, d\tau
=
1 - \cos t .
\end{equation*}
Example 6.2.5.
An equation containing an integral of the unknown function is called an integral equation. Consider
\begin{equation*}
x(t) - t = \int_0^t x(\tau) \, d\tau ,
\end{equation*}
where we wish to solve for \(x(t)\text{.}\) We apply the Laplace transform to get
\begin{equation*}
X(s) - \frac{1}{s^2} = \frac{1}{s} X(s) ,
\end{equation*}
where \(X(s) = \mathcal{L} \bigl\{ x(t) \bigr\}\text{.}\) Thus
\begin{equation*}
X(s) = \frac{1}{s(s-1)} = \frac{1}{s-1} - \frac{1}{s} .
\end{equation*}
The inverse Laplace transform gives
\begin{equation*}
x(t) = e^t - 1 .
\end{equation*}
Subsection 6.2.6 Periodic functions
The reader might ask: What about periodic functions as our input \(f(t)\text{?}\) That is, a function \(f(t)\) where \(f(t) = f(t+P)\) for some constant \(P\) (the period). Well, let us compute \(F(s)\text{:}\)
\begin{equation*}
\begin{split}
F(s)
& = \int_0^\infty e^{-st} f(t) \, dt
= \int_0^P e^{-st} f(t) \, dt
+ \int_P^\infty e^{-st} f(t) \, dt
\\
& = \int_0^P e^{-st} f(t) \, dt
+ \int_0^\infty e^{-s(t+P)} f(t+P) \, dt
= \int_0^P e^{-st} f(t) \, dt
+ e^{-Ps} \int_0^\infty e^{-st} f(t) \, dt
\\
& = \int_0^P e^{-st} f(t) \, dt
+ e^{-Ps} F(s) .
\end{split}
\end{equation*}
Solving for \(F(s)\text{,}\) we get
\begin{equation*}
F(s)
=
\frac{1}{1-e^{-Ps}}
\int_0^P e^{-st} f(t) \, dt .
\end{equation*}
As before, computing the inverse would be more complex and possibly involve consulting a table. Let us not worry about computing the inverse here.
Example 6.2.6.
Suppose our function \(f(t)\) is a version of the sawtooth, that is, let \(f(t) = t\) for \(0 \leq t < 1\) and use \(f(t)=f(t+1)\) to extend it periodically. So \(f(t) = t-1\) for \(1 \leq t < 2\text{,}\) \(f(t) = t-2\) for \(2 \leq t < 3\text{,}\) etc. Then \(P=1\) and a short computation with integration by parts gets
\begin{equation*}
F(s) =
\frac{1}{1-e^{-s}}
\int_0^1 e^{-st} t \, dt
=
\frac{1}{1-e^{-s}}
\left(
\frac{-e^{-s}}{s}-\frac{e^{-s}}{s^2} + \frac{1}{s^2}
\right)
=
\frac{-e^{-s}}{(1-e^{-s})s} + \frac{1}{s^2} .
\end{equation*}
Exercises 6.2.7 Exercises
6.2.2.
Using the Heaviside function write down the piecewise function that is 0 for \(t < 0\text{,}\) \(t^2\) for \(t\) in \([0,1]\) and \(t\) for \(t > 1\text{.}\)
6.2.3.
Using the Laplace transform solve
\begin{equation*}
m x'' + c x' + k x = 0 , \quad x(0) = a, \quad x'(0) = b ,
\end{equation*}
where \(m > 0\text{,}\) \(c > 0\text{,}\) \(k > 0\text{,}\) and \(c^2 - 4km > 0\) (system is overdamped).
6.2.4.
Using the Laplace transform solve
\begin{equation*}
m x'' + c x' + k x = 0 , \quad x(0) = a, \quad x'(0) = b ,
\end{equation*}
where \(m > 0\text{,}\) \(c > 0\text{,}\) \(k > 0\text{,}\) and \(c^2 - 4km < 0\) (system is underdamped).
6.2.5.
Using the Laplace transform solve
\begin{equation*}
m x'' + c x' + k x = 0 , \quad x(0) = a, \quad x'(0) = b ,
\end{equation*}
where \(m > 0\text{,}\) \(c > 0\text{,}\) \(k > 0\text{,}\) and \(c^2 = 4km\) (system is critically damped).
6.2.6.
Solve \(x'' + x = u(t-1)\) for initial conditions \(x(0) = 0\) and \(x'(0) = 0\text{.}\)
6.2.7.
Show the ``differentiation of the transform’’ property. Suppose \(\mathcal{L} \bigl\{ f(t) \bigr\} = F(s)\text{,}\) then show
\begin{equation*}
\mathcal{L} \bigl\{ -t f(t) \bigr\} = F'(s) .
\end{equation*}
Hint: Differentiate under the integral sign.
6.2.8.
Solve \(x''' + x = t^3 u(t-1)\) for initial conditions \(x(0) = 1\) and \(x'(0) =
0\text{,}\) \(x''(0) = 0\text{.}\)
6.2.9.
Show the second shifting property: \(\mathcal{L} \bigl\{ f(t-a) \, u(t-a) \bigr\}
= e^{-as} \mathcal{L} \bigl\{ f(t) \bigr\}\text{.}\)
6.2.10.
Let us think of the mass-spring system with a rocket from Example 6.2.2. We noticed that the solution kept oscillating after the rocket stopped running. The amplitude of the oscillation depends on the time that the rocket was fired (for 4 seconds in the example).
- Find a formula for the amplitude of the resulting oscillation in terms of the amount of time the rocket is fired.
- Is there a nonzero time (if so what is it?) for which the rocket fires and the resulting oscillation has amplitude 0 (the mass is not moving)?
6.2.11.
Define
\begin{equation*}
f(t) =
\begin{cases}
{(t-1)}^2 & \text{if } \; 1 \leq t < 2, \\
3-t & \text{if } \; 2 \leq t < 3, \\
0 & \text{otherwise} .
\end{cases}
\end{equation*}
- Sketch the graph of \(f(t)\text{.}\)
- Write down \(f(t)\) using the Heaviside function.
- Solve \(x''+x=f(t)\text{,}\) \(x(0)=0\text{,}\) \(x'(0) = 0\) using Laplace transform.
6.2.12.
Find the transfer function for \(m x'' + c x' + kx = f(t)\text{,}\) \(x(0)=0\text{,}\) \(x'(0)=0\text{.}\)
6.2.13.
Suppose \(Lx = f(t)\text{,}\) \(x(0)=0\text{,}\) \(x'(0)=0\) for the input function \(f(t) = t\) has the output \(x(t) = 2e^{-t} + te^{-t} + (t-2)\text{.}\)
- Find \(F(s)\text{,}\) \(X(s)\text{,}\) and the transfer function \(H(s)\text{.}\)
- If the input is instead \(f(t) = \sin(t)\) instead, find the new output \(x(t)\text{.}\)
6.2.14.
Suppose \(f(t) = 1\) if \(0 \leq t < 1\) and \(f(t)=0\) if \(1 \leq t < 2\text{,}\) and then extend periodically for all \(t \geq 0\) so that \(f(t)=f(t+2)\text{.}\) Compute the Laplace transform \(F(s)\text{.}\)
6.2.101.
Using the Heaviside function \(u(t)\text{,}\) write down the function
\begin{equation*}
f(t) =
\begin{cases}
0 & \text{if } \; \phantom{1 \leq {}} t < 1 , \\
t-1 & \text{if } \; 1 \leq t < 2 , \\
1 & \text{if } \; 2 \leq t .
\end{cases}
\end{equation*}
Answer.
\(f(t) =
(t-1)\bigl(u(t-1) - u(t-2)\bigr) +
u(t-2)\)
6.2.102.
Solve \(x''-x = (t^2-1) u(t-1)\) for initial conditions \(x(0)=1\text{,}\) \(x'(0) = 2\) using the Laplace transform.
Answer.
\(x(t) = (2e^{t-1}-t^2-1) u(t-1) -\frac{1}{2}e^{-t}+\frac{3}{2}e^t\)
6.2.103.
Find the transfer function for \(x' + x = f(t)\text{,}\) \(x(0)=0\text{.}\)
Answer.
\(H(s) = \frac{1}{s+1}\)
6.2.104.
Suppose \(Lx = f(t)\text{,}\) \(x(0)=0\text{,}\) \(x'(0)=0\) for the input function \(f(t) = 4\) has the output \(x(t) = 1 - \cos(2 t)\text{.}\) Find \(F(s)\text{,}\) \(X(s)\text{,}\) and the transfer function \(H(s)\text{.}\)
Answer.
\(F(s) = \frac{4}{s}\text{,}\) \(X(s) = \frac{1}{s} - \frac{s}{s^2+4}\text{,}\) \(H(s) = \frac{1}{s^2+4}\text{.}\)
