Docsity
Docsity

Prepare for your exams
Prepare for your exams

Study with the several resources on Docsity


Earn points to download
Earn points to download

Earn points by helping other students or get them with a premium plan


Guidelines and tips
Guidelines and tips

Inverse Laplace Transforms: Solving Differential Equations with Given Initial Conditions, Schemes and Mind Maps of Differential Equations

Solutions to various inverse Laplace transform problems, which are used to find the functions f(x) given their Laplace transforms F(s). different methods to solve these problems, including using properties of Laplace transforms and partial fractions. The applications of inverse Laplace transforms are demonstrated in the context of differential equations.

What you will learn

  • What are the properties of inverse Laplace transforms used in solving differential equations?
  • What is the role of partial fractions in solving inverse Laplace transform problems?
  • What is the procedure to solve a differential equation using inverse Laplace transforms?
  • How do you find the inverse Laplace transform of a given function?
  • How do you find the constants in the solution of a differential equation using inverse Laplace transforms?

Typology: Schemes and Mind Maps

2021/2022

Uploaded on 09/12/2022

maraiah
maraiah 🇺🇸

3.3

(3)

253 documents

1 / 12

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Math 334 6.3. INVERSE LAPLACE TRANSFORMS 90
theorem we get
L[cos ax] = 1
1e
2πs
aZ2π
a
0
esx cos ax dx =1
1e
2πs
aesx
s2+a2(asin ax scos ax)
2π
a
0
=1
1e
2πs
a(e
2πs
a(s)(s)
s2+a2)=s
s2+a2.
Example 6.23. Finf L[f(x)] for the function
f(x) = (1 0 < x 61
1 1 < x 62, f(x+ 2n) = f(x)nZ.
Solution
The function fis p eriodic with period 2, so we have
L[f(x)] = 1
1e2sZ2
0
esxf(x)dx =1
1e2sZ1
0
esx dx Z2
1
esx dx
=1
1e2se2s2es+ 1
s=(1 es)2
s(1 e2s)=1es
s(1 + es=es/2es/2
s(es/2+es/2)=1
stanh( 1
s).
6.3 Inverse Laplace Transforms
Recall the solution procedure outlined in Figure 6.1. The final stage in that solution pro cedure involves
calulating inverse Laplace transforms. In this section we look at the problem of finding inverse Laplace
transforms. In other words, given F(s), how do we find f(x) so that F(s) = L[f(x)].
We begin with a simple example which illustrates a small problem on finding inverse Laplace transforms.
Example 6.24. Consider the functions
f(x) = x2,and g(x) =
x2x6= 2,3
48 x= 2
π x = 3
.
Then L[f(x)] = L[g(x)] = 2
s3. Since an integral is not affected by the changing of its integrand at a few
isolated points, more than one function can have the same Laplace transform.
Example 6.24 illustrates that inverse Laplace transforms are not unique. However, it can b e shown that, if
several functions have the same Laplace transform, then at most one of them is continuous. This prompts
us to make the following definition.
Definition 6.25. The inverse Laplace transform of F(s), denoted L1[F(s)], is the function fdefined
on [0,) which has the fewest number of discontinuities and satisfies
L[f(x)] = F(s).
pf3
pf4
pf5
pf8
pf9
pfa

Partial preview of the text

Download Inverse Laplace Transforms: Solving Differential Equations with Given Initial Conditions and more Schemes and Mind Maps Differential Equations in PDF only on Docsity!

Math 334 6.3. INVERSE LAPLACE TRANSFORMS 90

theorem we get

L[cos ax] =

1 − e

− 2 πs a

∫ (^2) aπ

0

e−sx^ cos ax dx =

1 − e

− 2 πs a

e−sx s^2 + a^2

(a sin ax − s cos ax)

2 π a

0

=

1 − e − (^2) aπs

e

− 2 πs a (^) (−s) − (−s) s^2 + a^2

s s^2 + a^2

Example 6.23. Finf L[f (x)] for the function

f (x) =

1 0 < x 6 1 − 1 1 < x 6 2

, f (x + 2n) = f (x) ∀n ∈ Z.

Solution The function f is periodic with period 2, so we have

L[f (x)] =

1 − e−^2 s

0

e−sxf (x) dx =

1 − e−^2 s

0

e−sx^ dx −

1

e−sx^ dx

1 − e−^2 s

e−^2 s^ − 2 e−s^ + 1 s

(1 − e−s)^2 s(1 − e−^2 s)

1 − e−s s(1 + e−s^

es/^2 − e−s/^2 s(es/^2 + e−s/^2 )

s tanh(

s

6.3 Inverse Laplace Transforms

Recall the solution procedure outlined in Figure 6.1. The final stage in that solution procedure involves calulating inverse Laplace transforms. In this section we look at the problem of finding inverse Laplace transforms. In other words, given F (s), how do we find f (x) so that F (s) = L[f (x)].

We begin with a simple example which illustrates a small problem on finding inverse Laplace transforms.

Example 6.24. Consider the functions

f (x) = x^2 , and g(x) =

x^2 x 6 = 2, 3 48 x = 2 −π x = 3

Then L[f (x)] = L[g(x)] =

s^3

. Since an integral is not affected by the changing of its integrand at a few

isolated points, more than one function can have the same Laplace transform.

Example 6.24 illustrates that inverse Laplace transforms are not unique. However, it can be shown that, if several functions have the same Laplace transform, then at most one of them is continuous. This prompts us to make the following definition.

Definition 6.25. The inverse Laplace transform of F (s), denoted L−^1 [F (s)], is the function f defined on [0, ∞) which has the fewest number of discontinuities and satisfies

L[f (x)] = F (s). ◭

Math 334 6.3. INVERSE LAPLACE TRANSFORMS 91

Example 6.26.

  1. L−^1 [

s^3 ] = x^2.

2. L−^1 [

s s^2 + 9

] = cos 3x.

3. L−^1 [

s − 1 s^2 − 2 s + 5

] = L−^1 [

s − 1 (s − 1)^2 + 4 ] = ex^ L−^1 [

s s^2 + 4 ] = ex^ cos 2x. (using property 1 of Theorem 6.17 in reverse)

The inverse Laplace transform is a linear operator.

Theorem 6.27. If L−^1 [F (s)] and L−^1 [G(s)] exist, then L−^1 [αF (s) + βG(s)] = αL−^1 [F (s)] + βL−^1 [G(s)].

Proof Starting from the right hand side we have

L[αL−^1 [F (s)] + βL−^1 [G(s)]] = αL[L−^1 [F (s)]] + βL[L−^1 [G(s)]] = αF (s) + βG(s).

The result follows.

Most of the properties of the Laplace transform can be reversed for the inverse Laplace transform.

Theorem 6.28. If L−^1 [F (s)] = f (x), then the following hold:

  1. L−^1 [F (s + a)] = e−axf (x);
  2. L−^1 [sF (s)] = f ′(x), if f (0) = 0;
  3. L−^1 [

s F (s)] =

∫ (^) x

0

f (t) dt;

  1. L−^1 [e−asF (s)] = ua(x) f (x − a).

Proof

  1. L[e−axf (x)] = F (s + a) from Theorem 6.17, property 1. The result follows.
  2. L[f ′(x))] = −f (0) + sF (s) from Theorem 6.17, property 4. The result follows.

3. L[

∫ (^) x

0

f (t) dt] =

s

F (s), from Theorem 6.17, property 5. The result follows.

  1. L[ua(x)f (x − a)] = e−asL[f (x)] = e−asF (s), from Theorem 6.19. The result follows.

Example 6.29. Find L−^1 [

s(s^2 + 1)

].

Solution

We can write

s(s^2 + 1)

s

F (s), where F (s) =

s^2 + 1

. Then f (x) = L−^1 [F (s)] = sin x, so we get

L−^1 [

s(s^2 + 1)

] = L−^1 [

s F (s)] =

∫ (^) x

0

f (t) dt =

∫ (^) x

0

sin t dt = 1 − cos x.

Math 334 6.4. APPLICATIONS TO DIFFERENTIAL EQUATIONS 93

Example 6.35 (quadratic factors). Find L−^1 [

2 s^2 + 10s s^2 − 2 s + 5 (s + 1)].

Solution We write the expression in the form

L−^1 [

2 s^2 + 10s s^2 − 2 s + 5

(s + 1)] =

A(s − 1) + B (s − 1)^2 + 4

C

s + 1

Solving for the constants yields: A = 3, B = 8, and C = −1. Thus, we get

L−^1 [

2 s^2 + 10s s^2 − 2 s + 5

(s + 1)] = 3L−^1 [

s − 1 (s − 1)^2 + 4

] + 4L−^1 [

(s − 1)^2 + 4

] − L−^1 [

s + 1

]

= 3ex^ cos 2x + 4ex^ sin 2x − e−x.

6.4 Applications to Differential Equations

The easiest way to see how to apply Laplace transforms to differential equations is to work through some examples.

Example 6.36. Solve the following initial value problem:

y′′^ − y′^ = 2x, y(0) = 1, y′(0) = − 2.

Solution Method 1 (the old approach) First solve the homogeneous equation: y′′^ − y = 0.

y = erx^ =⇒ r^2 − r = 0 =⇒ r = 0, 1 =⇒ yh(x) = c 1 + c 2 ex.

Now look for a particular solution: yp(x) = x(Ax + B) = Ax^2 + Bx. Plug into the DE to get

y′′ p − y′ p = 0 =⇒

2 A − B = 0

− 2 A = 2

A = − 1

B = − 2

=⇒ yp(x) = −x^2 − 2 x.

Thus, we have y(x) = c 1 + c 2 ex^ − x^2 − 2 x, y′(x) = c 2 ex^ − 2 x − 2.

Apply the initial conditions:

y(0) = 1 y′(0) = − 2

c 1 + c 2 = 1 c 2 − 2 = − 2

c 1 = 1 c 2 = 0

=⇒ y(x) = 1 − x^2 − 2 x.

Method 2 (using Laplace transforms) Take Laplace transforms of the DE:

L[y′′] − L[y′] = 2L[x] =⇒

[

s^2 Y (s) − sy(0) − y′(0)

]

− [sY (s) − y(0)] =

s^2 =⇒ s^2 Y (s) − s + 2 − sY (s) + 1 =

s^2

=⇒ Y (s) =

s^3 − 3 s^2 + 2 s^3 (s − 1)

(s − 1)(s^2 − 2 s − 2) s^3 (s − 1)

s

s^2

s^3

Math 334 6.4. APPLICATIONS TO DIFFERENTIAL EQUATIONS 94

Finally, taking the inverse Laplace transform, we arrive at the final solution:

y(x) = L−^1 [Y (s)] = 1 − 2 x − x^2.

Example 6.37. Solve the following initial value problem:

y′′^ − 2 y′^ + 5y = − 8 e−x, y(0) = 2, y′(0) = 12.

Solution Take Laplace transforms of the DE:

L[y′′] − 2 L[y′] + 5L[y] = − 8 L[e−x] =⇒

[

s^2 Y (s) − sy(0) − y′(0)

]

− 2 [sY (s) − y(0)] + 5Y (s) =

s^2 =⇒ s^2 Y (s) − 2 s − 12 − 2 sY (s) + 4 + 5Y (s) =

s + 1

=⇒ Y (s) =

2 s^2 + 10s (s^2 − 2 s + 5)(s + 1)

=⇒ Y (s) = 3 (s − 1) (s − 1)^2 + 4

(s − 1)^2 + 4

s + 1

Finally, taking the inverse Laplace transform, we arrive at the final solution:

y(x) = 3ex^ cos 2x + 4ex^ sin 2x − e−x.

Example 6.38. Solve the following initial value problem:

y′′^ − 2 y′^ + 5y = − 8 e^7 −x, y(7) = 2, y′(7) = 12.

Solution It appears that we can not use Laplace transforms since L[y′] = sY (s) − y(0), and we don’t know y(0). But we can get around this by moving the initial point (in this case x 0 = 7) to the origin by means of a translation.

Let t = x − 7 and w(t) = y(x). Then we get

w′(t) = y′(x), w′′(t) = y′′(x), w(0) = y(7), w′(0) = y′(7),

so, the initial value problem becomes

w′′^ − 2 w′^ + 5w = − 8 e−t, w(0) = 2, w′(0) = 12.

This is just the initial value problem we had in Example 6.40. The solution is

w(t) = 3et^ cos 2t + 4et^ sin 2t − e−t.

The solution to the original problem is

y(x) = w(x − 7) = 3ex−^7 cos[2(x − 7)] + 4ex−^7 sin[2(x − 7)] − e−(x−7).

Next we consider an initial value problem with discontinuous forcing.

Math 334 6.5. CONVOLUTION 96

This is a linear ODE in Y (s). Look for an integrating factor μ:

μ′ μ

s

s 2

=⇒ ln μ = 3 ln s −

s^2 4

=⇒ μ = s^3 e−s

(^2) / 4 .

The ODE for Y becomes:

d ds [μY (s)] = μ

Y ′(s) +

μ′ μ Y (s)

μ 2 s^2

s 2 e−s

(^2) / 4 .

Integrating yields:

μY (s) = −

s 2 e−s

(^2) / 4 ds = e−s

(^2) / 4

  • C =⇒ Y (s) =

s^3

C

s^3 es

(^2) / 4 .

It remains to determin the value of the constant if integration C. There is no auxiliary condition inherited from the original ODE. To get an appropriate condition to enable us to determine C, we utilize Theorem 6.12 that states that Y (s) → 0 as s → ∞. Therfore

lim s→∞ Y (s) = 0 =⇒ C = 0 =⇒ Y (s) =

s^3

Finally, taking the inverse Laplace transform, we arrive at the final solution:

y(x) = L−^1 [

s^3

] =

x^2 2

We can summarize the application of Laplace transforms to differential equations as follows.

  1. For an ODE with constant coefficients, the equation for the Laplace transform is of the form AY = B.
  2. For an ODE with polynomial coefficients in x, the equation for the Laplace transform is an ODE with polynomial coefficients in s.
  3. If the coefficient functions of the ODE are liear in x, the ODE for Y is a first order ODE.
  4. The auxiliary condition to use when solving a first order ODE in Y is lim s→∞ Y (s) = 0.

6.5 Convolution

Consider the following initial value problem:

y′′^ + y = g(x), y(0) = y′(0) = 0.

Take the Laplace transform of the equation to get:

s^2 Y (s) − sy(0) − y′(0) + Y (s) = G(s) =⇒ Y (s) =

G(s) s^2 + 1

= F (s)G(s), where F (s) =

s^2 + 1

We would like to express the solution y(x) in terms of f (x) and g(x), i.e. we would like to express L−^1 [F (s)G(s)] in terms of L−^1 [F (s)] and L−^1 [G(s)]. To do this, we define a special type of product of functions. Let f, g ∈ P C(0, ∞).

Math 334 6.5. CONVOLUTION 97

Definition 6.41. The convolution of f and g, denoted f ∗ g, is defined as:

(f ∗ g)(x) :=

∫ (^) x

0

f (x − t)g(t) dt.

Example 6.42.

(1) 1 ∗ x =

∫ (^) x

0

1 · t dt =

x^2 2

(2) x ∗ x^2 =

∫ (^) x

0

(x − t) · t^2 dt =

x^4 12

Theorem 6.43. The convolution product satisfies the following properties:

  1. f ∗ g = g ∗ f ; (convolution product is commutative)
  2. f ∗ (g + h) = f ∗ g + f ∗ h; (convolution product is distributive over addition)
  3. f ∗ (g ∗ h) = (f ∗ g) ∗ h; (convolution product is associative)
  4. f ∗ 0 = 0.

Proof Exercise.

Remark. While the convolution product has many of the properties of ordinary multiplication of functions, it different in that it has no multiplicative identity element, i.e. there is no function g with the property that g ∗ f 6 = f for all funtions f. ◭

Theorem 6.44. If (i) f, g ∈ P C(0, ∞); (ii) F (s) = L[f (x)] and G(s) = L[g(x)], then

L[f ∗ g] = F (s)G(s), or equivalently L−^1 [F (s)G(s)] = (f ∗ g)(x).

Proof

L[f ∗ g] =

0

e−sx(f ∗ g)(x) dx =

0

e−sx

∫ (^) x

0

f (x − t)g(t) dt dx

0

t

e−sxf (x − t)g(t) dx dt =

0

0

e−s(ξ+t)f (ξ)g(t) dξ dt

0

e−sξf (ξ) dξ

0

e−stg(t) dt

= F (s)G(s).

Example 6.45. Find L−^1 [

s(s^2 + 4)

].

Solution Method 1 (the old approach)

L−^1 [

s(s^2 + 4)

] = L−^1 [

s

s s^2 + 4

] =

(1 − cos 2x).

Method 2 (using convolution)

Math 334 6.6. THE DELTA FUNCTION 99

The integral of the force acting over an interval of time ic called “impulse:” I =

∫ (^) t 1

t 0

F (t) dt. From Newton’s

Law, F = ma, we get

I =

∫ (^) t 1

t 0

F (t) dt =

∫ (^) t 1

t 0

m dv dt

dt = mv(t 1 ) − mv(t 0 ) (i.e. impulse = change in momentum).

Consider an impulse over shorter and shorter time intervals.

v

t (^0) t 1 t t (^0) t 1 t

F

v

t (^0) t 1 t t t t 0 1

F

v

t t 0

t 1

F

t 1 t^0 t

Figure 6.5: A plot of a localized force and impulse for decreasing time intervals.

The exact nature of the force in the interval [t 0 , t 1 ] is frequently unknown. What usually is known is the state of the system before and after the application of the force. These duration of the force is often so short, that it is convenient to think of the force as acting instantaneously.

Math 334 6.6. THE DELTA FUNCTION 100

v

t 0 t 1 t t t t 0 1

F

Figure 6.6: A plot of a localized force and impulse with simplified approximation.

We define the following functions:

δh(t) =

2 h |t| 6 h 0 |t| > h

for h > 0.

t -h h

Figure 6.7: A plot of a δh(t).

Thes functions have the following property:

∫ (^) ∞

−∞

δh(t) dt =

∫ (^) h

−h

2 h

dt =

2 h

(h − (−h)) = 1.

Definition 6.47. The Dirac Delta Function is defined implicitly by the following properties:

  1. δ(t) = 0 for t 6 = 0;

−∞

δ(t) = 1. ◭

Remark. The Dirac delta function is not a function in the usual sense. However, one can “think” of it as follows:

δ(t) = lim h→ 0

δh(t) =

∞ t = 0 0 t 6 = 0