University:
Massachusetts Institute of Technology
Course:
2.161 | Signal Processing: Continuous and Discrete
Academic year:

2019
Views:

255

Pages:

6
Author:

aesirdidactical

MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING 2.161 Signal Processing – Continuous and Discrete The Dirac Delta and Unit-Step Functions 1 The Dirac Delta (Impulse) Function The Dirac delta function is a non-physical, singularity function with the following deﬁnition δ(t) = but with the requirement that 0 undeﬁned ∞ −∞ for t = 0 at t = 0 (1) δ(t)dt = 1, (2) that is, the function has unit area. Despite its name, the delta function is not truly a function. Rigorous treatment of the Dirac delta requires measure theory or the theory of distributions. 1 /T @ T ( t ) @ (t) 1 1 /T 2 1 /T 3 1 /T 4 0 T 1 T 2 T 3 T t 0 4 a ) U n it p u ls e s o f d iffe r e n t e x te n ts t b ) T h e im p u ls e fu n c tio n Figure 1: Unit pulses and the Dirac delta function. Figure 1 shows a unit pulse function δT (t), that is a brief rectangular pulse function of extent T , deﬁned to have a constant amplitude 1/T over its extent, so that the area T × 1/T under the pulse is unity: ⎧ ⎪ for t ≤ 0 ⎨ 0 1/T 0 T . The Dirac delta function (also known as the impulse function) can be deﬁned as the limiting form of the unit pulse δT (t) as the duration T approaches zero. As the extent T of δT (t) decreases, the amplitude of the pulse increases to maintain the requirement of unit area under the function, and δ(t) = lim δT (t). T →0 (4) The impulse is therefore deﬁned to exist only at time t = 0, and although its value is strictly undeﬁned at that time, it must tend toward inﬁnity so as to maintain the property of unit area in the limit. 1 In addition to the one-sided unit pulse described above, the delta function may be considered as the limit of several other functions δ(t) = lim δa (t) a→0 where δa (t) is sometimes called a nascent delta function. The following are some commonly used approximations: ⎧ ⎪ ⎪ ⎪ 1 ⎪ ⎪ ⎨ δa (t) = for |t| ≤ a ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ 0 ⎧ ⎪ ⎪ ⎪ 2 4 ⎪ ⎪ − 2 |t| ⎨ a a δa (t) = ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ 0 a 2 for |t| > a 2 for |t| ≤ a 2 δa (t) = δa (t) = 1 aπ 1 δa (t) = aπ triangular function a for |t| > 2 1 2 2 δa (t) = √ e−t /a a π δa (t) = “rect” function Gaussian function 1 a π a2 + t2 Cauchy-Lorentz distribution 1 −|t/a| e a Cauchy ϕ function sin (t/a) t/a sin (t/a) t/a sinc function 2 sinc2 function All of these functions have the property that as a → 0, the function becomes impulse-like in the neighborhood of t = 0, with unit area under the curve. The impulse function is used extensively in the study of linear systems, both spatial and tem poral. Although true impulse functions are not found in nature, they are approximated by short duration, high amplitude phenomena such as a hammer impact on a structure, or a lightning strike on a radio antenna. As we will see below, the response of a causal linear system to an impulse deﬁnes its response to all inputs. 1.1 1.1.1 Properties Shift An impulse occurring at t = a is δ(t − a). 1.1.2 The strength of an impulse Because the amplitude of an impulse is inﬁnite, it does not make sense to describe a scaled impulse by its amplitude. Instead, the strength of a scaled impulse Kδ(t) is deﬁned by its area K. 2 1.1.3 The “Sifting” Property of the Impulse When an impulse appears in a product within an integrand, it has the property of ”sifting” out the value of the integrand at the point of its occurrence: ∞ −∞ f (t)δ(t − a)dt = f (a) (5) This is easily seen by noting that δ(t − a) is zero except at t = a, and for its inﬁnitesimal duration f (t) may be considered a constant and taken outside the integral, so that ∞ −∞ f (t)δ(t − a)dt = f (a) ∞ −∞ δ(t − a)dt = f (a) (6) from the unit area property. 1.1.4 Scaling A helpful identity is the scaling property: ∞ −∞ δ(αt)dt = ∞ −∞ and so du 1 = |α| |α| 1 δ(t). |α| δ (αt) = 1.1.5 δ(u) Laplace Transform L {δ(t)} = ∞ δ(t)e−st dt = 1 0− by the sifting property. 1.1.6 Fourier Transform F {δ(t)} = ∞ −∞ δ(t)e−jΩt dt = 1 by the sifting property. 2 Practical Applications of the Dirac Delta Function • The most important application of δt in linear system theory is directly related to its Laplace transform property, L {δ(t)} = 1. Consider a SISO LTI system with transfer function H(s), with input u(t) and output y(t), so that in the Laplace domain Y (s) = H(s)U (s). If the input is u(t) = δ(t), so that U (s) = 1, then Y (s) = H(s).1, and through the inverse Laplace transform y(t) = h(t) = L−1 {H(s)} . where h(t) is deﬁned as the system’s impulse response. The impulse response completely characterizes the system, in the sense that it allows computation of the transfer function (and hence the diﬀerential equation). 3 • The impulse response h(t) is used in the convolution integral. • In signal processing the delta function is used to create a Dirac comb (also known as an impulse train, or Shah function): ∞ ΔT (t) = δ(t − nT ) n=−∞ is used in sampling theory. A continuous waveform f (t) is sampled by multiplication by the Dirac comb ∞ f ∗ (t) = f (t)ΔT (t) = f (t − nT )δ(t − nT ), n=−∞ where f ∗ (t) is the sampled waveform, producing a train of weighted impulses. 3 The Heaviside (Unit-Step) Function The unit-step function us (t), also known as the Heaviside function, is a discontinuous function with a value 0 for negative arguments, and unity for positive arguments, as shown in Fig. 2. K I J J Figure 2: The Heaviside, or unit-step function. The unit-step is an antiderivative of the Dirac delta in the sense that us (t) = t −∞ δ(u)du for t = 0 where this relation may not hold (or even make sense) for t = 0. The value of us (0) is subject to some debate, and is often given as us (0) = 0, us (0) = 1, or us (0) = 0.5. In practice the value rarely matters because us (t) usually appears within an integral. However, the value us (0) = 0.5 has some conceptual appeal, especially when using an analytic approximation to the unit-step. A useful representation of us (t) is us (t) = 1 + sgn(t) 2 where sgn() denotes the signum (sign) function ⎧ ⎪ ⎨ −1 sgn(t) = ⎪ ⎩ for t < 0 for t = 0 . for t > 0 0 1 Clearly this deﬁnition of us (t) supports the convention us (0) = 0.5. 4 3.1 Analytic Approximations The limiting form of many sigmoid type functions centered on t = 0 may be used as an approxi mation to the Heaviside function, for example: 1 1 (1 + tanh (kt)) = k→∞ 2 1 + e−2kt � � 1 1 us (t) = lim + arctan (kt) k→∞ 2 π 1 us (t) = lim (1 + erf (kt)) k→∞ 2 us (t) = 3.2 lim Integral representation An integral representation is often useful us (t) = lim − �→0+ 3.3 3.3.1 1 2π � ∞ −∞ 1 e−jtτ dτ τ + j� Properties Shift A unit-step occurring at t = a is simply us (t − a). 3.3.2 Derivative and Antiderivative While not strictly diﬀerentiable, the derivative of the unit-step is taken to be the Dirac delta d us (t) = δ(t). dt The integral of us (t) is � t −∞ � us (µ)dµ = 0 t for t < 0 for t ≥ 0 which is the unit ramp function r(t). 3.3.3 Laplace Transform L {us (t)} = 3.3.4 1 s Fourier Transform Although the Fourier integral does not converge directly for us (t), it is easy to show that 1 1 F {us (t)} = δ (jΩ) + 2 jΩ See the class handout on Fourier transforms for details. 5 3.4 Practical Applications of the Unit-Step Function • Perhaps the most common use of us (t) is as a multiplicative weighting function to create causal (one-sided) functions. For example, a causal exponential time function may be expressed as f (t) = us (t)e−σt . • As a test input signal for characterizing the response of linear systems throughout linear system theory and feedback control theory. • As the basis for synthesizing zero-order (staircase angular pulse ⎧ ⎪ for ⎨ 0 3 for f (t) = ⎪ ⎩ 0 for or boxcar ) functions. For example a rect t < 1.5 1.5 ≤ t < 5 t≥5 may be written f (t) = 3 (us (t − 1.5) − us (t − 5)) . Similarly, a zero-order approximation rT (t) to a continuous unit ramp r(t) = t for t > 0, and updated at intervals T , may be written as a sum of shifted step functions: rT (t) = ∞ � n=0 6 us (t − nT )

The Dirac Delta and Unit-Step Functions

Related Documents

Recommended Documents

Description

Related Documents

Free up your schedule!

Take 5 seconds to unlock