Site Network: Home | Blogcrowds | Gecko and Fly | Free Online Libraries | Free Antivirus | Free CAD Software | Free Software | Telegraph TV | About

Advertisement

Your Ad Here

UTeM BMFR Google Search Engine

Custom Search

Laplace Transform

The Laplace transform is an integral transform perhaps second only to the Fourier transform in its utility in solving physical problems. The Laplace transform is particularly useful in solving linear ordinary differential equations such as those arising in the analysis of electronic circuits.

The (unilateral) Laplace transform L (not to be confused with the Lie derivative, also commonly denoted L) is defined by

 L_t[f(t)](s)=int_0^inftyf(t)e^(-st)dt,
(1)

where f(t) is defined for =0" width="27" border="0" height="14"> (Abramowitz and Stegun 1972). The unilateral Laplace transform is almost always what is meant by "the" Laplace transform, although a bilateral Laplace transform is sometimes also defined as

 L_t^((2))[f(t)](s)=int_(-infty)^inftyf(t)e^(-st)dt
(2)

(Oppenheim et al. 1997). The unilateral Laplace transform L_t[f(t)](s) is implemented in Mathematica as LaplaceTransform[f[t], t, s].

The inverse Laplace transform is known as the Bromwich integral, sometimes known as the Fourier-Mellin integral (see also the related Duhamel's convolution principle).

A table of several important one-sided Laplace transforms is given below.

fL_t[f(t)](s)conditions
11/s
t1/(s^2)
t^n(n!)/(s^(n+1))=0" width="57" border="0" height="14">
t^a(Gamma(a+1))/(s^(a+1))-1" width="62" border="0" height="14">
e^(at)1/(s-a)
cos(omegat)s/(s^2+omega^2)omega in R
sin(omegat)a/(s^2+omega^2)|I[omega]|" width="54" border="0" height="14">
cosh(omegat)s/(s^2-omega^2)|R[omega]|" width="59" border="0" height="14">
sinh(omegat)a/(s^2-omega^2)|I[omega]|" width="54" border="0" height="14">
e^(at)sin(bt)b/((s-a)^2+b^2)a+|I[b]|" width="73" border="0" height="14">
e^(at)cos(bt)(s-a)/((s-a)^2+b^2)b in R
delta(t-c)e^(-cs)
H_c(t)0" width="111" border="0" height="59">
J_0(t)1/(sqrt(s^2+1))
J_n(at)((sqrt(s^2+a^2)-s)^n)/(a^nsqrt(s^2+a^2))=0" width="57" border="0" height="14">

In the above table, J_0(t) is the zeroth-order Bessel function of the first kind, delta(t) is the delta function, and H_c(t) is the Heaviside step function.

The Laplace transform has many important properties. The Laplace transform existence theorem states that, if f(t) is piecewise continuous on every finite interval in [0,infty) satisfying

 |f(t)|<=Me^(at)
(3)

for all t in [0,infty), then L_t[f(t)](s) exists for all a" width="29" border="0" height="14">. The Laplace transform is also unique, in the sense that, given two functions F_1(t) and F_2(t) with the same transform so that

 L_t[F_1(t)](s)=L_t[F_2(t)](s)=f(s),
(4)

then Lerch's theorem guarantees that the integral

 int_0^aN(t)dt=0
(5)

vanishes for all 0" width="30" border="0" height="14"> for a null function defined by

 N(t)=F_1(t)-F_2(t).
(6)

The Laplace transform is linear since

L_t[af(t)+bg(t)]=int_0^infty[af(t)+bg(t)]e^(-st)dt
(7)
=aint_0^inftyfe^(-st)dt+bint_0^inftyge^(-st)dt
(8)
=aL_t[f(t)]+bL_t[g(t)].
(9)

The Laplace transform of a convolution is given by

L_t[f(t)*g(t)]=L_t[f(t)]L_t[g(t)]
(10)
L_t^(-1)[FG]=L_t^(-1)[F]*L_t^(-1)[G].
(11)

Now consider differentiation. Let f(t) be continuously differentiable n-1 times in [0,infty). If |f(t)|<=Me^(at), then

 L_t[f^((n))(t)](s)=s^nL_t[f(t)]-s^(n-1)f(0)-s^(n-2)f^'(0)-...-f^((n-1))(0).
(12)

This can be proved by integration by parts,

L_t[f^'(t)](s)=infty)int_0^ae^(-st)f^'(t)dt" width="112" border="0" height="35">
(13)
=infty){[e^(-st)f(t)]_0^a+sint_0^ae^(-st)f(t)dt}" width="209" border="0" height="35">
(14)

(15)
=infty)[e^(-sa)f(a)-f(0)+sint_0^ae^(-st)f(t)dt]" width="235" border="0" height="35">
(16)

(17)
=sL_t[f(t)]-f(0).
(18)

Continuing for higher-order derivatives then gives

 L_t[f^('')(t)](s)=s^2L_t[f(t)](s)-sf(0)-f^'(0).
(19)

This property can be used to transform differential equations into algebraic equations, a procedure known as the Heaviside calculus, which can then be inverse transformed to obtain the solution. For example, applying the Laplace transform to the equation

 f^('')(t)+a_1f^'(t)+a_0f(t)=0
(20)

gives

 {s^2L_t[f(t)](s)-sf(0)-f^'(0)}+a_1{sL_t[f(t)](s)-f(0)}+a_0L_t[f(t)](s)=0
(21)
 L_t[f(t)](s)(s^2+a_1s+a_0)-sf(0)-f^'(0)-a_1f(0)=0,
(22)

which can be rearranged to

 L_t[f(t)](s)=(sf(0)+f^'(0)+a_1f(0))/(s^2+a_1s+a_0).
(23)

If this equation can be inverse Laplace transformed, then the original differential equation is solved.

The Laplace transform satisfied a number of useful properties. Consider exponentiation. If L_t[f(t)](s)=F(s) for alpha" width="30" border="0" height="14"> (i.e., F(s) is the Laplace transform of f), then L_t[e^(at)f](s)=F(s-a) for a+alpha" width="52" border="0" height="14">. This follows from

F(s-a)=int_0^inftyfe^(-(s-a)t)dt
(24)
=int_0^infty[f(t)e^(at)]e^(-st)dt
(25)
=L_t[e^(at)f(t)](s).
(26)

The Laplace transform also has nice properties when applied to integrals of functions. If f(t) is piecewise continuous and |f(t)|<=Me^(at), then

 L_t[int_0^tf(t^')dt^']=1/sL_t[f(t)](s).
(27)

SEE ALSO: Bilateral Laplace Transform, Bromwich Integral, Fourier-Mellin Integral, Fourier Transform, Integral Transform, Laplace-Stieltjes Transform, Operational Mathematics, Unilateral Laplace Transform

0 Comments:

Post a Comment



UTeM BMFR Blog Directories

UTeM BMFR MyBlogLog