Line 68: Line 68:
 
<math>L[f+g]=\int_{0}^{\infty} (f(t)+g(t)) e^{-st} dt=\int_{0}^{\infty} f(t) e^{-st} dt + \int_{0}^{\infty} g(t) e^{-st} dt = L[f]+L[g]</math>,
 
<math>L[f+g]=\int_{0}^{\infty} (f(t)+g(t)) e^{-st} dt=\int_{0}^{\infty} f(t) e^{-st} dt + \int_{0}^{\infty} g(t) e^{-st} dt = L[f]+L[g]</math>,
  
L[cf]=\int_{0}^{\infty} cf(t) e^{-st} dt = c \int_{0}^{\infty} f(t) e^{-st} dt = cL[f]</math>.
+
<math>L[cf]=\int_{0}^{\infty} cf(t) e^{-st} dt = c \int_{0}^{\infty} f(t) e^{-st} dt = cL[f]</math>.
  
 
Therefore, it is linear.
 
Therefore, it is linear.
 +
  
 
We can do the same test for the inverse transform, which is also linear:
 
We can do the same test for the inverse transform, which is also linear:

Revision as of 03:49, 22 November 2017

Laplace Transforms

A slecture by Yijia Wen

7.0 Abstract

Sometimes an ODE can be really complex and hard to solve by any basic methods we have looked at in previous tutorials. For example, it may involve exponential functions, trigonometric functions, Heaviside functions, and anything else you can imagine. Hence, a famous French physicist Pierre-Simon Laplace found a transform method, which converts the functions in "time-domain" to "complex number-domain", to overcome the problem. It is great for transforming the calculation of differentiation and integration to the simple algebraic calculation.


7.1 Concept

Here is a basic concept map from my course note when first learning Laplace transform:

$ \begin{bmatrix} ODE & in & y(t) \\ "Time-domain" \end{bmatrix} → \begin{bmatrix} Laplace & Transform \\ L=[y(t)] \end{bmatrix} → \begin{bmatrix} Algebraic & equation & in & Y(s) \\ "Complex & number-domain" \end{bmatrix} $;

$ \begin{bmatrix} Algebraic & equation & in & Y(s) \\ "Complex & number-domain" \end{bmatrix} → \begin{bmatrix} Inverse & Laplace & Transform \\ L^{-1}[Y(s)] \end{bmatrix} → \begin{bmatrix} ODE & in & y(t) \\ "Time-domain" \end{bmatrix} $.


That means for some complicated ODEs, Laplace Transform is used to convert them to some algebraic equations. Then inverse Laplace Transform is used to convert the solution back.


Laplace Transform is defined as $ Y(s)=L[y(t)]=\int_{0}^{\infty} y(t) e^{-st} dt $, which is also time function $ y(t) $ expressed in the "complex frequency" domain.

Let's do a simple example. Consider the time function $ y(t)=e^{2t} $. Follow the definition above, we have $ Y(s)=L[e^{2t}]=\int_{0}^{\infty} e^{2t} e^{-st} dt=\int_{0}^{\infty} e^{(2-s)t} dt=\lim_{b \to \infty} {\int_{0}^{b} e^{2-s}t dt}=\lim_{b \to \infty} {\frac{e^{2-s}t}{2-s} |_{0}^{b}}=\lim_{b \to \infty} {\frac{1}{2-s} (e^{(2-s)b}-e^0)} $

$ =\begin{cases} \frac{1}{s-2}, & if & s>2 \\ undefined, & if & s≤2 \end{cases} $.


If doing Laplace Transform to a derivative $ \frac{dy}{dt} $:

$ L[\frac{dy}{dt}]=\int_{0}^{\infty} \frac{dy}{dt} e^{-st} dt $, by the definition of Laplace Transform,

$ =[ye^{-st} - \int (-s e^{-st} y dt)]_{0}^{\infty} $, integration by part,

$ =[ye^{-st} + s \int y e^{-st} dt]_{0}^{\infty} $,

$ =[ye^{-st}]_{0}^{\infty} + sY(s) $, again by the definition of Laplace Transform,

$ =sY(s) - y(0) (= sL[y] - y(0)) $,


Now consider the second order derivative $ \frac{d^2y}{dt^2} $,

$ L[\frac{d^2y}{dt^2}]=L[\frac{d}{dt} (\frac{dy}{dt})]=s L[\frac{dy}{dt}] - \frac{dy}{dx}(0) = s(sL[y]-y(0)) - \frac{dy}{dt}(0) = s^2 L[y] - sy(0) - \frac{dy}{dt}(0) $.

Here we can see the regular pattern and infer the general form $ L[\frac{d^ny}{dt^n}]=s^n Y(s) - s^{n-1}y(0) - ... - \frac{d^{n-1}y}{dt^{n-1}(0)} $.


7.2 Linearity of Laplace Transform

By the definition of linear transformation, if a function $ f(x) $ is linear, then $ f(a+bx)=f(a)+bf(x) $. Test the linearity for the definition of Laplace Transform $ L[y(t)]=\int_{0}^{\infty} y(t) e^{-st} dt $ by this property,

$ L[f+g]=\int_{0}^{\infty} (f(t)+g(t)) e^{-st} dt=\int_{0}^{\infty} f(t) e^{-st} dt + \int_{0}^{\infty} g(t) e^{-st} dt = L[f]+L[g] $,

$ L[cf]=\int_{0}^{\infty} cf(t) e^{-st} dt = c \int_{0}^{\infty} f(t) e^{-st} dt = cL[f] $.

Therefore, it is linear.


We can do the same test for the inverse transform, which is also linear:

$ L^{-1} [X(s)+Y(s)]=L^{-1}[X(s)] + L^{-1}[Y(s)] $,

$ L^{-1} [cX(s)]=cL^{-1} [Y(s)] $.


7.3 Heaviside Unit Step Function for Discontinuous Functions


7.4 Exercises


7.5 References

Institute of Natural and Mathematical Science, Massey University. (2017). 160.204 Differential Equations I: Course materials. Auckland, New Zealand.

Robinson, J. C. (2003). An introduction to ordinary differential equations. New York, NY., USA: Cambridge University Press.

Alumni Liaison

Meet a recent graduate heading to Sweden for a Postdoctorate.

Christine Berkesch