ConceptComplete

Inverse Laplace Transform

The inverse Laplace transform recovers the original time-domain function from its transform, primarily using partial fraction decomposition and table lookups.


Definition and Uniqueness

Definition4.5Inverse Laplace transform

The inverse Laplace transform of F(s)F(s) is the function f(t)f(t) such that L{f}=F\mathcal{L}\{f\} = F. When it exists for piecewise continuous functions of exponential order, the inverse is unique (up to values at points of discontinuity). It is given by the Bromwich integral:

f(t)=12πiγiγ+iestF(s)dsf(t) = \frac{1}{2\pi i}\int_{\gamma - i\infty}^{\gamma + i\infty}e^{st}F(s)\,ds

where γ\gamma is greater than the real part of all singularities of FF.

Theorem4.5Lerch's theorem

If F(s)=G(s)F(s) = G(s) for all s>s0s > s_0, and both f,gf, g are piecewise continuous of exponential order with L{f}=F\mathcal{L}\{f\} = F and L{g}=G\mathcal{L}\{g\} = G, then f(t)=g(t)f(t) = g(t) at all points of continuity.


Partial Fraction Method

RemarkStrategy for rational transforms

Most inverse Laplace transforms in practice involve rational functions F(s)=P(s)/Q(s)F(s) = P(s)/Q(s) with degP<degQ\deg P < \deg Q. The method:

  1. Factor Q(s)Q(s) into linear and irreducible quadratic factors.
  2. Decompose FF into partial fractions.
  3. Invert each term using the table of standard transforms.
ExamplePartial fractions example

Find L1{3s+2(s1)(s2+4)}\mathcal{L}^{-1}\left\{\frac{3s+2}{(s-1)(s^2+4)}\right\}.

Partial fractions: 3s+2(s1)(s2+4)=As1+Bs+Cs2+4\frac{3s+2}{(s-1)(s^2+4)} = \frac{A}{s-1} + \frac{Bs+C}{s^2+4}.

3s+2=A(s2+4)+(Bs+C)(s1)3s+2 = A(s^2+4) + (Bs+C)(s-1). Setting s=1s=1: 5=5A5 = 5A, so A=1A = 1.

Comparing s2s^2: 0=A+B0 = A + B, so B=1B = -1. Comparing constants: 2=4AC=4C2 = 4A - C = 4 - C, so C=2C = 2.

F(s)=1s1+s+2s2+4=1s1ss2+4+2s2+4.F(s) = \frac{1}{s-1} + \frac{-s+2}{s^2+4} = \frac{1}{s-1} - \frac{s}{s^2+4} + \frac{2}{s^2+4}.

Inverting: f(t)=etcos2t+sin2tf(t) = e^t - \cos 2t + \sin 2t.


The Bromwich Integral and Residues

Theorem4.6Inverse via residues

If F(s)F(s) is rational with poles s1,,sns_1, \ldots, s_n, then

f(t)=k=1nRes(estF(s),sk).f(t) = \sum_{k=1}^n \text{Res}\left(e^{st}F(s), s_k\right).

ExampleInverse via residues

For F(s)=1(s1)2(s+2)F(s) = \frac{1}{(s-1)^2(s+2)}:

At s=1s = 1 (pole of order 2): Res=ddsests+2s=1=tet3et9=(3t1)et9\text{Res} = \frac{d}{ds}\frac{e^{st}}{s+2}\big|_{s=1} = \frac{te^t \cdot 3 - e^t}{9} = \frac{(3t-1)e^t}{9}.

At s=2s = -2 (simple pole): Res=e2t(21)2=e2t9\text{Res} = \frac{e^{-2t}}{(-2-1)^2} = \frac{e^{-2t}}{9}.

f(t)=(3t1)et+e2t9.f(t) = \frac{(3t-1)e^t + e^{-2t}}{9}.


Applications to Systems

RemarkSystems of ODEs via Laplace transforms

For a system x=Ax+g(t)\mathbf{x}' = A\mathbf{x} + \mathbf{g}(t) with x(0)=x0\mathbf{x}(0) = \mathbf{x}_0, the Laplace transform gives

sX(s)x0=AX(s)+G(s)    X(s)=(sIA)1(x0+G(s)).s\mathbf{X}(s) - \mathbf{x}_0 = A\mathbf{X}(s) + \mathbf{G}(s) \implies \mathbf{X}(s) = (sI - A)^{-1}(\mathbf{x}_0 + \mathbf{G}(s)).

The solution is x(t)=L1{(sIA)1}x0+L1{(sIA)1G(s)}\mathbf{x}(t) = \mathcal{L}^{-1}\{(sI-A)^{-1}\}\mathbf{x}_0 + \mathcal{L}^{-1}\{(sI-A)^{-1}\mathbf{G}(s)\}. The matrix L1{(sIA)1}=eAt\mathcal{L}^{-1}\{(sI-A)^{-1}\} = e^{At} is the matrix exponential.

Example2x2 system via Laplace

Solve x=3xyx' = 3x - y, y=x+yy' = x + y, x(0)=1x(0) = 1, y(0)=0y(0) = 0.

(sIA)X=x0(sI - A)\mathbf{X} = \mathbf{x}_0: (s311s1)(XY)=(10)\begin{pmatrix}s-3 & 1 \\ -1 & s-1\end{pmatrix}\begin{pmatrix}X\\Y\end{pmatrix} = \begin{pmatrix}1\\0\end{pmatrix}.

det(sIA)=(s3)(s1)+1=s24s+4=(s2)2\det(sI-A) = (s-3)(s-1)+1 = s^2-4s+4 = (s-2)^2. So X=s1(s2)2=1s2+1(s2)2X = \frac{s-1}{(s-2)^2} = \frac{1}{s-2} + \frac{1}{(s-2)^2}.

x(t)=e2t+te2t=(1+t)e2tx(t) = e^{2t} + te^{2t} = (1+t)e^{2t}. Similarly Y=1(s2)2Y = \frac{1}{(s-2)^2}, so y(t)=te2ty(t) = te^{2t}.