๐ Introduction to the Dirac Delta Function
The Dirac delta function, denoted as $\delta(t)$, is not a function in the traditional sense but rather a distribution. It's often described as an idealized impulse. Informally, it's zero everywhere except at $t = 0$, where it's infinite, and its integral over the entire real line is equal to one:
- ๐ Definition: $\delta(t) = 0$ for $t \neq 0$, and $\int_{-\infty}^{\infty} \delta(t) dt = 1$.
- ๐ก Sifting Property: The key property we'll use is the sifting property: $\int_{-\infty}^{\infty} f(t) \delta(t - a) dt = f(a)$, where $f(t)$ is a continuous function at $t = a$. This property allows us to "sift out" the value of $f(t)$ at a specific point.
- ๐ Units: Note that $\delta(t)$ has units of 1/(time), so that the integral has units of 1 (dimensionless).
๐ History and Background
The Dirac delta function was introduced by physicist Paul Dirac in his work on quantum mechanics. While it initially lacked a rigorous mathematical foundation, it was later formalized by Laurent Schwartz in the theory of distributions.
- โ๏ธ Quantum Mechanics: Dirac used it to describe the probability density of a particle at a specific point in space.
- ๐ Signal Processing: It is also used to represent impulses in signal processing.
- ๐งโ๐ซ Mathematical Formalization: Laurent Schwartz developed the theory of distributions to provide a rigorous mathematical basis for the Dirac delta function and similar concepts.
๐ Key Principles for Solving Differential Equations with Dirac Delta
When solving differential equations involving the Dirac delta function, keep these principles in mind:
- ๐ Homogeneous Solution: First, find the homogeneous solution to the differential equation (i.e., the solution when the Dirac delta term is zero).
- ๐ฅ Impulse Effect: The Dirac delta function represents an impulse at a specific time. The solution will have a discontinuity in its derivative (or higher derivatives, depending on the order of the differential equation) at that time.
- ๐งฉ Matching Conditions: Use the sifting property to determine the jump in the derivative at the point where the Dirac delta function is applied.
- โ๏ธ Continuity: The solution itself (without derivatives) should be continuous at the point where $\delta(t)$ is applied.
๐งช Worked Problems
Example 1: A Simple First-Order Equation
Solve the differential equation $y' + 2y = \delta(t)$, with the initial condition $y(0) = 0$.
- ๐ Homogeneous Solution: The homogeneous equation is $y' + 2y = 0$, which has the solution $y_h(t) = Ce^{-2t}$.
- ๐ฅ Effect of Delta Function: Integrate the equation from $0^-$ to $0^+$: $\int_{0^-}^{0^+} (y' + 2y) dt = \int_{0^-}^{0^+} \delta(t) dt$. This gives $[y(0^+) - y(0^-)] + 2\int_{0^-}^{0^+} y dt = 1$. Since $y$ is continuous, the integral term is zero. Thus, $y(0^+) - y(0^-) = 1$. Since $y(0^-) = y(0) = 0$, we have $y(0^+) = 1$.
- ๐งฉ Complete Solution: For $t > 0$, the equation is $y' + 2y = 0$, with the initial condition $y(0^+) = 1$. The solution is $y(t) = e^{-2t}$ for $t > 0$. Combining this with the initial condition, we have
$y(t) = \begin{cases}
0, & t < 0 \\
e^{-2t}, & t > 0
\end{cases}$.
Example 2: A Second-Order Equation
Solve the differential equation $y'' + 4y = \delta(t)$, with initial conditions $y(0) = 0$ and $y'(0) = 0$.
- ๐ Homogeneous Solution: The homogeneous equation is $y'' + 4y = 0$, which has the general solution $y_h(t) = A\cos(2t) + B\sin(2t)$.
- ๐ฅ Effect of Delta Function: Integrate the equation from $0^-$ to $0^+$: $\int_{0^-}^{0^+} (y'' + 4y) dt = \int_{0^-}^{0^+} \delta(t) dt$. This gives $[y'(0^+) - y'(0^-)] + 4\int_{0^-}^{0^+} y dt = 1$. Since $y$ is continuous, the integral term is zero. Thus, $y'(0^+) - y'(0^-) = 1$. Since $y'(0^-) = y'(0) = 0$, we have $y'(0^+) = 1$. Also, $y(0^+) = y(0^-) = 0$.
- ๐งฉ Complete Solution: For $t > 0$, the equation is $y'' + 4y = 0$, with initial conditions $y(0) = 0$ and $y'(0) = 1$. Applying these conditions to the homogeneous solution, we get $A = 0$ and $2B = 1$, so $B = 1/2$. Therefore, for $t > 0$, $y(t) = (1/2)\sin(2t)$. Combining with the initial conditions, we have
$y(t) = \begin{cases}
0, & t < 0 \\
(1/2)\sin(2t), & t > 0
\end{cases}$.
Example 3: Forced Harmonic Oscillator
Solve the differential equation $y'' + 2y' + 5y = \delta(t-1)$, with initial conditions $y(0) = 0$ and $y'(0) = 0$.
- ๐ Homogeneous Solution: The homogeneous equation is $y'' + 2y' + 5y = 0$. The characteristic equation is $r^2 + 2r + 5 = 0$, with roots $r = -1 \pm 2i$. The general solution is $y_h(t) = e^{-t}(A\cos(2t) + B\sin(2t))$.
- ๐ฅ Effect of Delta Function: Integrate the equation from $1^-$ to $1^+$: $\int_{1^-}^{1^+} (y'' + 2y' + 5y) dt = \int_{1^-}^{1^+} \delta(t-1) dt$. This gives $[y'(1^+) - y'(1^-)] + 2[y(1^+) - y(1^-)] + 5\int_{1^-}^{1^+} y dt = 1$. Since $y$ is continuous at $t=1$, the second and third terms are zero. Thus, $y'(1^+) - y'(1^-) = 1$.
- ๐งฉ Complete Solution: For $0 < t < 1$, the solution is $y(t) = 0$ because of the initial conditions and no forcing function before $t=1$. For $t > 1$, we have $y'' + 2y' + 5y = 0$, and $y(1) = 0$. We need to find $y'(1)$. Let $y'(1^-) = 0$ (since $y(t) = 0$ before $t=1$). Then $y'(1^+) = 1$. Now, $y(t) = e^{-(t-1)}(A\cos(2(t-1)) + B\sin(2(t-1)))$. Because $y(1) = 0$, we have $A = 0$. Then, $y'(t) = e^{-(t-1)}(2B\cos(2(t-1)) - B\sin(2(t-1))) - e^{-(t-1)}B\sin(2(t-1))$. So, $y'(1) = 2B = 1$, and $B = \frac{1}{2}$. Thus, for $t > 1$, $y(t) = \frac{1}{2}e^{-(t-1)}\sin(2(t-1))$. Therefore, the entire solution is
$y(t) = \begin{cases}
0, & t < 1 \\
\frac{1}{2}e^{-(t-1)}\sin(2(t-1)), & t > 1
\end{cases}$.
๐ Real-World Examples
The Dirac delta function appears in various fields:
- ๐ก Signal Processing: Modeling impulsive noise in signals.
- ๐งฒ Physics: Describing point charges or point masses.
- ๐จ Engineering: Analyzing the response of a system to a sudden impact.
๐ก Conclusion
Understanding how to work with the Dirac delta function is crucial in many areas of science and engineering. By mastering the sifting property and understanding its effect on differential equations, you can solve a wide range of problems involving impulsive forces or signals. Practice is key! Keep solving those equations, and you'll become a Dirac delta pro in no time! ๐