Laplace Transform

4.7 Impulse and Dirac Delta Function

In prior sections, we explored initial value problems for second-order differential equations with constant coefficients, focusing on cases where the forcing function, [asciimath]f(t)[/asciimath], is either continuous or piecewise continuous on the interval [asciimath][0, oo)[/asciimath].

 [asciimath]ay''+by'+cy=f(t)\ ;[/asciimath] [asciimath]y(0)=k_0, \ \ y'(0)=k_1[/asciimath]

Now, let’s turn our attention to a different type of forcing function: one that represents an impulsive force. Impulsive forces are characterized by very large magnitudes over extremely short time intervals, effectively appearing as a sudden “jolt” or “spike” in the system. Such impulses occur in various contexts, including electrical circuits during a switch-on event, mechanical systems during a collision, or any scenario where a sudden, significant force is applied for a brief period.

A. Dirac Delta Function

To mathematically model these impulsive forces, we use the Dirac Delta function, denoted as [asciimath]delta(t)[/asciimath]. The Dirac Delta function is not a function in the traditional sense but rather a generalized function or distribution with the following properties.

1. Zero everywhere except at zero:

  [asciimath]delta(t)={(0 \ if\ t!=0),(oo \ if \ t=0):}[/asciimath]

2. Integral equals one:

[asciimath]int_-oo^oo delta(t) = 1[/asciimath]

3. Sifting property:

[asciimath]int_-oo^oo f(t) delta(t) = f(0)[/asciimath]  for any [asciimath]f(t)[/asciimath] that is continuous on the interval that contains [asciimath]t=0[/asciimath]

By shifting the argument [asciimath]t[/asciimath] in [asciimath]delta(t)[/asciimath], we can model impulses that occur at times other than [asciimath]t=0[/asciimath]. The shifted Dirac Delta function, [asciimath]delta(t-a)[/asciimath], has a spike at [asciimath]t=a[/asciimath] and is defined as

 [asciimath]delta(t-a)={(0 \if \ t!=a),(oo \if \ t=a):}[/asciimath]

Thus the sifting property extends to

4. Sifting at [asciimath]t=a[/asciimath]:

 [asciimath]int_-oo^oo f(t) delta(t-a) = f(a)[/asciimath] for any [asciimath]f(t)[/asciimath] that is continuous on the interval that contains [asciimath]t=a[/asciimath]

B. Laplace Transform of the Dirac Delta Function

The Laplace Transform provides a convenient way to handle the Dirac Delta function in the context of solving differential equations. The transform of a shifted Dirac Delta function is given by

 [asciimath]\mathcal{L}{delta(t-a)}=e^(-as)[/asciimath](4.7.1)

Understanding the Dirac Delta function and its properties is crucial for modeling and analyzing systems subjected to impulsive forces.

 

Example 4.7.1: Solve IVP with Impulsive Forcing Function

Find the solution to the initial value problem

 [asciimath]y''+16y=4delta(t-pi),[/asciimath]   [asciimath]y(0)=1,[/asciimath]   [asciimath]y'(0)=0[/asciimath]

Show/Hide Solution

 

Taking the Laplace transform of the equation, applying Equation 4.7.1 with [asciimath]a=pi[/asciimath] to the Delta function, yields

 [asciimath]s^2 Y(s)-s+16Y(s)=4e^(-pis)[/asciimath]

Solving for [asciimath]Y(s)[/asciimath], we obtain

 [asciimath]Y(s)=(4e^(-pis) +s)/(s^2+16)[/asciimath]

 [asciimath]=e^(-pis) 4/(s^2+16)+s/(s^2+16)[/asciimath]

Computing the inverse Laplace transform gives

 [asciimath]y(t)=u(t-pi) sin(4(t-pi))[/asciimath] [asciimath]+cos(4t)[/asciimath]

Which equivalently is

 [asciimath]y(t)= {(cos(4t) if tltpi),(sin(4t)+cos(4t) if  t>=pi):}[/asciimath]

   [asciimath]={(cos(4t) if tltpi),(sqrt(2)sin(4t +pi/4) if  t>=pi):}[/asciimath]

The below figure shows [asciimath]y(t)[/asciimath]. The impulsive force is applied and adds momentum to the system at [asciimath]t=pi[/asciimath]. For comparison, the dotted line represents the undisturbed system.

 

Try an Example

 

Section 4.7 Exercises

  1. Solve the initial value problem

     [asciimath]y'' + 25 y = delta(t-4), \quad \ y(0) = y'(0) = 0.[/asciimath]

    Show/Hide Answer

     [asciimath]y(t)={(0 \ if t<4),(1/5 sin(5(t-4)) \if t>=4):}[/asciimath]

  2. Solve the initial value problem

     [asciimath]y'' + 4 y = 80 e^(4 t) + delta(t-7), \quad \ y(0) = 11,\ y'(0) = 32[/asciimath]

    Show/Hide Answer

     [asciimath]y(t)={(7cos(2 t) + 8sin(2 t) + 4 e^(4 t) \ if t<7),(7cos(2 t) + 8sin(2 t) + 4 e^(4 t) + 1/2 sin(2(t-7)) \if t>=7):}[/asciimath]

License

Icon for the Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License

Differential Equations Copyright © 2024 by Amir Tavangar is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License, except where otherwise noted.

Share This Book