Calculate the Laplace Transform of a Function

The Laplace transform is an integral transform used in solving differential equations of constant coefficients. This transform is also extremely useful in physics and engineering.

While tables of Laplace transforms are widely available, it is important to understand the properties of the Laplace transform so that you can construct your own table.

Preliminaries

  • Let <math>f(t)</math> be a function defined for <math>t \geq 0.</math> Then we define the Laplace transform of <math>f(t)</math> as the following function for every value of <math>s</math> where the integral converges.
    • <math>F(s) = \mathcal{L}\{f(t)\} = \int_{0}^{\infty}f(t)e^{-st}\mathrm{d}t</math>
  • By applying a Laplace transform to a function, we are transforming a function from the t-domain (or time domain) to the s-domain (or Laplace domain), where <math>F(s)</math> is a complex function of a complex variable. In doing so, we are transforming the problem into a domain that is hopefully easier to solve in.
  • Obviously, the Laplace transform is a linear operator, so we can consider the transform of a sum of terms by doing each integral separately.
    • <math>\int_{0}^{\infty}[af(t) + bg(t)]e^{-st}\mathrm{d}t = a\int_{0}^{\infty}f(t)e^{-st}\mathrm{d}t + b\int_{0}^{\infty}g(t)e^{-st}\mathrm{d}t</math>
  • Remember that the Laplace transform only exists if the integral converges. If the function <math>f(t)</math> is discontinuous anywhere, we must be very careful to ensure that we split the boundaries of the integral to avoid blowup.

Steps

The Basics

  1. Substitute the function into the definition of the Laplace transform. Conceptually, calculating a Laplace transform of a function is extremely easy. We will use the example function <math>f(t) = e^{at}</math> where <math>a</math> is a (complex) constant such that <math>\operatorname{Re}(s) < \operatorname{Re}(a).</math>
    • <math>\mathcal{L}\{e^{at}\} = \int_{0}^{\infty}e^{at}e^{-st}\mathrm{d}t</math>
  2. Evaluate the integral using any means possible. In our example, our evaluation is extremely simple, and we need only use the fundamental theorem of calculus. In other more complicated cases, techniques like integration of parts or differentiation under the integral may be used. Our constraint that <math>\operatorname{Re}(s) < \operatorname{Re}(a)</math> means that the integrand converges, i.e. goes to 0 as <math>t\to\infty.</math>
    • <math>\begin{align}\mathcal{L}\{e^{at}\} &= \int_{0}^{\infty}e^{(a-s)t}\mathrm{d}t \\ &= \frac{e^{(a-s)t}}{a-s}\Bigg|_{0}^{\infty} \\ &= \frac{1}{s-a}\end{align}</math>
    • Notice that this gives us two Laplace transforms for "free": the sine and cosine functions, if we consider the related function <math>e^{iat}</math> via Euler's formula. Then in the denominator, we would have <math>s - ia,</math> and all that remains is to take the real and imaginary parts of this result. We could also just evaluate directly, but that would take a bit more work.
      • <math>\mathcal{L}\{\cos at\} = \operatorname{Re}\left(\frac{1}{s-ia}\right) = \frac{s}{s^{2} + a^{2}}</math>
      • <math>\mathcal{L}\{\sin at\} = \operatorname{Im}\left(\frac{1}{s-ia}\right) = \frac{a}{s^{2} + a^{2}}</math>
  3. Evaluate the Laplace transform of the power function. Before moving on, we must determine the transform of the power function, for the property of linearity allows us to determine the transform for all polynomials. The power function is the function <math>t^{n},</math> where <math>n</math> is any positive integer. We can use integration by parts to determine a recursive rule.
    • <math>\mathcal{L}\{t^{n}\} = \int_{0}^{\infty}t^{n}e^{-st}\mathrm{d}t = \frac{n}{s}\mathcal{L}\{t^{n-1}\}</math>
    • Our result is not written explicitly, but from substituting a few values of <math>n,</math> a clear pattern emerges (try it yourself), from which we can determine the following result.
      • <math>\mathcal{L}\{t^{n}\} = \frac{n!}{s^{n+1}}</math>
    • We can also determine Laplace transforms of fractional powers by using the Gamma function. This allows us to find transforms of functions like <math>f(t) = \sqrt{t}.</math>
      • <math>\mathcal{L}\{t^{n}\} = \frac{\Gamma(n+1)}{s^{n+1}}</math>
      • <math>\mathcal{L}\{t^{1/2}\} = \frac{\Gamma(3/2)}{s^{3/2}} = \frac{\sqrt{\pi}}{2s\sqrt{s}}</math>
    • Although functions with fractional powers must contain branch cuts (recall that for any complex numbers <math>z</math> and <math>\alpha,</math> we rewrite <math>z^{\alpha}</math> as <math>e^{\alpha\operatorname{Log}z}</math>), we can always define them such that the branch cuts lie in the left half-plane in order to avoid analyticity issues.

Properties of the Laplace Transform

  1. Determine the Laplace transform of a function multiplied by <math>e^{at}</math>. The results in the previous section have allowed us to take a glimpse at some interesting properties of the Laplace transform. The Laplace transform of functions like cosine, sine, and the exponential function seem to be simpler than the transform of the power function. We will see that multiplication by <math>e^{at}</math> in the t-domain corresponds to a shift in the s-domain.
    • <math>\mathcal{L}\{e^{at}f(t)\} = \int_{0}^{\infty}f(t)e^{-(s-a)t}\mathrm{d}t = F(s-a)</math>
    • This property immediately allows us to find transforms of functions like <math>f(t) = e^{3t}\sin 2t</math> without having to directly evaluate the integral.
      • <math>\mathcal{L}\{e^{3t}\sin 2t\} = \frac{2}{(s-3)^{2} + 4}</math>
  2. Determine the Laplace transform of a function multiplied by <math>t^{n}</math>. Let's consider multiplying by <math>t</math> first. Then from the definition, we can differentiate under the integral to obtain a surprisingly clean result.
    • <math>\begin{align}\mathcal{L}\{tf(t)\} &= \int_{0}^{\infty}tf(t)e^{-st}\mathrm{d}t \\ &= -\int_{0}^{\infty}f(t)\frac{\partial}{\partial s}e^{-st}\mathrm{d}t \\ &= -\frac{\mathrm{d}}{\mathrm{d}s}\int_{0}^{\infty}f(t)e^{-st}\mathrm{d}t \\ &= -\frac{\mathrm{d}F}{\mathrm{d}s}\end{align}</math>
    • By repeating this process, we arrive at the general result.
      • <math>\mathcal{L}\{t^{n}f(t)\} = (-1)^{n}\frac{\mathrm{d}^{n}F}{\mathrm{d}s^{n}}</math>
    • The exchange of the integral and the differentiation operators takes a bit of justification as far as rigor is concerned, but we will not justify it here except to note that the operation is allowed as long as our final answer makes sense. A bit of comfort can be sought in the fact that <math>s</math> and <math>t</math> are variables that are independent of one another.
    • Of course, using this property, Laplace transforms of functions like <math>t^{2}\cos 2t</math> are easily found without having to repeatedly use integration by parts.
      • <math>\mathcal{L}\{t^{2}\cos 2t\} = \frac{\mathrm{d}^{2}}{\mathrm{d}s^{2}}\frac{s}{s^{2}+4} = \frac{2s^{3}-24s}{(s^{2}+4)^{3}}</math>
  3. Determine the Laplace transform of a stretched function <math>f(at)</math>. Using the definition, we can also easily determine this transform using a u-substitution.
    • <math>\begin{align}\mathcal{L}\{f(at)\} &= \int_{0}^{\infty}f(at)e^{-st}\mathrm{d}t,\ \ u = at \\ &= \frac{1}{a}\int_{0}^{\infty}f(u)e^{-su/a}\mathrm{d}u \\ &= \frac{1}{a}F\left(\frac{s}{a}\right)\end{align}</math>
    • Previously, we found the Laplace transforms of <math>\sin at</math> and <math>\cos at</math> from the exponential function directly. We can use this property to arrive at the same result, starting from finding the real and imaginary parts of <math>\mathcal{L}\{e^{it}\} = \frac{1}{s-i}</math>.
  4. Determine the Laplace transform of a derivative <math>f^{\prime}(t)</math>. Unlike our previous results that saved a bit of labor from integration by parts, we must use integration by parts here.
    • <math>\begin{align}\mathcal{L}\{f^{\prime}(t)\} &= \int_{0}^{\infty}f^{\prime}(t)e^{-st}\mathrm{d}t, \ \ u = e^{-st},\ \mathrm{d}v = f^{\prime}(t)\mathrm{d}t \\ &= f(t)e^{-st}\Big|_{0}^{\infty} + s\int_{0}^{\infty}f(t)e^{-st}\mathrm{d}t \\ &= sF(s) -f(0)\end{align}</math>
    • Because the second derivative comes up in many physical applications, we list the Laplace transform of a second derivative as well.
      • <math>\mathcal{L}\{f^{\prime\prime}(t)\} = s^{2}F(s) - sf(0) - f^{\prime}(0)</math>
    • In general, it turns out that the Laplace transform of the nth derivative is given by the following result. This result is important in solving differential equations via Laplace transforms.
      • <math>\mathcal{L}\{f^{(n)}(t)\} = s^{n}F(s) - \sum_{k=0}^{n-1}s^{n-k-1}f^{(k)}(0)</math>

Series Methods

  1. Determine the Laplace transform of a periodic function. A periodic function is a function that satisfies the property <math>f(t) = f(t + nT),</math> where <math>T</math> is the period of the function and <math>n</math> is a positive integer. Periodic functions show up in many applications in signal processing and electrical engineering. Using a bit of manipulation, we arrive at the following answer.
    • <math>\begin{align}\mathcal{L}\{f(t)\} &= \int_{0}^{\infty}f(t)e^{-st}\mathrm{d}t \\ &= \sum_{n=0}^{\infty}\int_{nT}^{(n+1)T}f(t)e^{-st}\mathrm{d}t \\ &= \sum_{n=0}^{\infty}\int_{0}^{T}f(t+nT)e^{-s(t+nT)}\mathrm{d}t \\ &= \sum_{n=0}^{\infty}e^{-snT}\int_{0}^{T}f(t)e^{-st}\mathrm{d}t \\ &= \frac{1}{1-e^{-sT}}\int_{0}^{T}f(t)e^{-st}\mathrm{d}t\end{align}</math>
    • We see that the Laplace transform of a periodic function is related to the Laplace transform of one cycle of the function.
  2. See the article on Calculate the Laplace Transform of the Natural Logarithm. This integral cannot be evaluated using the fundamental theorem of calculus because the antiderivative cannot be expressed in terms of elementary functions. The article discusses a technique using the Integrate Using the Gamma Function and its Integrate Using Expansions of the Gamma Function to evaluate the natural log and its higher powers. The presence of the Euler-Mascheroni constant <math>\gamma</math> is enough to hint that the integral must be evaluated using series methods.
    • <math>\mathcal{L}\{\ln t\} = -\frac{\gamma+\ln s}{s}</math>
  3. Evaluate the Laplace transform of the (unnormalized) sinc function. The sinc function <math>\operatorname{sinc}(t) = \frac{\sin t}{t}</math> is a function widely encountered in signal processing, and may be recognizable from differential equations as equivalent to the zeroth-order spherical Bessel function of the first kind <math>j_{0}(x).</math> The Laplace transform of this function also cannot be computed in the standard way. We resort to transforming term-by-term, permissible because the individual terms are power functions and therefore their transforms certainly converge on the prescribed interval.
    • We begin by writing out the Taylor series of this function.
      • <math>\frac{\sin t}{t} = \sum_{n=0}^{\infty}\frac{(-1)^{n}t^{2n}}{(2n+1)!}</math>
    • Now we simply transform using the Laplace transform of the power function that we know. The factorials cancel out, and after staring at our expression, we recognize the Taylor series of the inverse tangent, the alternating series that looks like the Taylor series for the sine function but without the factorials.
      • <math>\begin{align}\mathcal{L}\left\{\frac{\sin t}{t}\right\} &= \sum_{n=0}^{\infty}\frac{(-1)^{n}(2n)!}{(2n+1)!}\frac{1}{s^{2n+1}} \\

&= \sum_{n=0}^{\infty}\frac{(-1)^{n}}{2n+1}\frac{1}{s^{2n+1}} \\ &= \tan^{-1}\frac{1}{s}\end{align}</math>

Video

{{#ev:youtube|OiNh2DswFt4}}

Related Articles

  • Calculate the Laplace Transform of the Natural Logarithm