Singularity Functions

From PrattWiki
Revision as of 00:09, 8 January 2010 by DukeEgr93 (talk | contribs)
Jump to navigation Jump to search

Singularity Functions are a class of functions that - you guessed it - contain singularities. One notation for a singularity function is given as adapted from p. 77 Chapra[1]

\( \begin{align} <x-x_0>^n&= \left\{ \begin{array}{cc} 0, & x\leq x_0\\ (x-a)^n, & x>x_0 \end{array} \right. \end{align} \)

where \(x\)is the independent variable and \(n\) and \(x_0\) are constants. Note that this definition is only useful when \(n\geq 0\). Negative values of \(n\) are used to denote the impulse function (when \(n=-1\)) or its derivatives. Also note that it is common to leave the singularity undefined when \(n=0\) and \(x=x_0\). That is to say, a fuller definition of the singularity function might be:

\( \begin{align} <x-x_0>^n&= \left\{ \begin{array}{ccc} \delta^{(1-n)}(x-x_0), & ~ & n<0\\ 0, & x<x_0 & ~\\ \mbox{Undefined}, & x=x_0 & n=0\\ 1, & x>x_0 & n=0\\ (x-x_0)^n, & x\geq x_0 & n>0\\ \end{array} \right. \end{align} \)

where the delta function \(\delta(x)\) and its derivatives will be defined below.

Basic Definitions and Derivations

Alternate Names for \(n\geq 0\)

Graph of Unit Singularity Functions

In courses such as ECE 54 and EGR 119, alternate names are given to different non-negative powers of the singularity function. Furthermore, singularity functions are generally written as functions of time rather than space, so:

\( \begin{align} \begin{array}{|c|c|c|}\hline \mbox{Singularity Function} & \mbox{Name} & \mbox{Alternate Symbol} \\ \hline \hline <t-a>^0 & \mbox{unit step function} & u(t) \\ \hline <t-a>^1 & \mbox{unit ramp function} & r(t)=t~u(t) \\ \hline \frac{1}{2}<t-a>^2 & \mbox{unit parabola function} & p(t)=\frac{1}{2}t^2~u(t) \\ \hline \end{array} \end{align} \)

For the last entry, you may be tempted to ask, "How is that the unit parabola?" The simplest explanation is that the unit ramp is the integral of the unit step, and the unit parabola is the integral of the unit parabola.

Derivation of Unit Impulse Functions

You can also take derivatives of the singularity functions. For \(n>0\), this is quite easy as the unit ramp and above are continuous. The difficulty comes in taking the derivative of the \(<t-a>^0\) case. Mathematically, call the derivative of the unit step function \(\delta(t)\); you can then find

\( \delta(t)=\begin{align} \frac{d}{dt}u(t)&= \left\{ \begin{array}{cc} 0, & t\neq 0\\ \infty, & t=0 \end{array} \right. \end{align} \)

But how large is that infinity? The answer comes in integrating this "delta function" (also known as the "unit impulse function" or just "impulse function"):

\( \begin{align} u(t)&=\int_{-\infty}^t\delta(\tau)~d\tau \end{align} \)

and by definition,

\( \begin{align} u(t)&= \left\{ \begin{array}{cc} 0, & t< 0\\ 1, & t>0 \end{array} \right. \end{align} \)

meaning

\( \begin{align} \int_{-\infty}^t\delta(\tau)~d\tau &= \left\{ \begin{array}{cc} 0, & t< 0\\ 1, & t>0 \end{array} \right. \end{align} \)

In other words (well...in words) the total area is 0 when integrating between negative infinity and just before 0 and the total area is 1 when integrating between negative infinity and anything positive. That must mean, at exactly \(t=0\), \(\delta(t)\) has an area of 1 while for all other times, it has an area of 0. That is to say,

\( \begin{align} \delta(t)&= \left\{ \begin{array}{cc} 0, & t\neq 0\\ \mbox{Area of 1}, & t=0 \end{array} \right. \end{align} \)

Variable Transformations for Unit Step and Unit Impulse

Three common transformations of the independent variable (in the examples, \(t\)) are scaling, shifting, and reversing. The following pairs show how scaling (multiplying the independent variable by some positive constant \(a\)), shifting (subtracting some value \(t_0\) from the independent variable), and reversing (multiplying the independent variable by -1) may be written in alternate ways. First, for the unit step function:

\( \begin{align} u(at)&=u(t)\\ u(-t)&=1-u(t)\\ u(t-t_0)&=1-u(t_0-t)\\ u(a(t-t_0))&=u(t-t_0)\\ u(-a(t-t_0))&=u(a(t_0-t))=u(t_0-t)=1-u(t-t_0) \end{align} \)

Note that a positive time scaling does not have an effect on the unit step function. Also note in the final two cases above that the independent variable is first shifted, then scaled or scaled and reversed.

For the unit impulse function, time scaling is a bit more complex. Fortunately, time reversal is much simpler:

\( \begin{align} \delta(at)&=\frac{1}{a}\delta(t)\\ \delta(-t)&=\delta(t)\\ \delta(t-t_0)&=\delta(t_0-t)\\ \delta(a(t-t_0))&=\frac{1}{a}\delta(t-t_0)\\ \delta(-a(t-t_0))&=\delta(a(t-t_0))=\frac{1}{a}\delta(t-t_0)\\ \end{align} \)

Convolution Integral Simplification with Step Function Product as Part of Integrand

The convolution integral

\( \begin{align} y(t)=\int_{-\infty}^{\infty}x(\tau)h(t-\tau)~d\tau&= \int_{-\infty}^{\infty}x(t-\tau)h(\tau)~d\tau \end{align} \)

will generate products of step functions in the integrand if both \(x(t)\) and \(h(t)\) are defined using step functions. There are two main ways to get a product of step functions as an integrand - either they are both "pointing" in the same direction or they are pointing in opposite directions.

Same Directions in Integrand

If the step functions are "pointing" in the same direction in \(\tau\) space - that is

\( \begin{align} \int_{-\infty}^{\infty} u(\tau-a)~u(b-t+\tau)~f(g(\tau))d\tau \end{align} \)

the key is to recognize that, for the integrand to be non-zero, you must satisfy the inequalities

\( \begin{align} \tau-a&>0 & b-t+\tau&>0,\mbox{ or}\\ a&<\tau & t-b&<\tau \end{align} \)

You will thus need to determine which of these presents the stricter condition. Then:

\( \begin{align} \int_{-\infty}^{\infty} u(\tau-a)~u(b-t+\tau)~f(g(\tau))d\tau&= \begin{cases} a<t-b & \int_{t-b}^{\infty} f(g(\tau))d\tau\\ a>t-b & \int_{a}^{\infty} f(g(\tau))d\tau \end{cases} \end{align} \)

which can be re-written using unit step functions as:

\( \begin{align} \left(u(t-b-a) \int_{t-b}^{\infty} f(g(\tau))d\tau\right) &+ \left(u(a+b-t) \int_{a}^{\infty} f(g(\tau))d\tau\right) & \end{align} \)

Opposite Directions in Integrand

This is the case that most often comes into play during convolution:

\( \begin{align} \int_{-\infty}^{\infty} u(\tau-a)~u(t-\tau-b)~f(g(\tau))d\tau \end{align} \)

where the first step function is right-sided in \(\tau\) and the second function is left-sided. For the integrand to be nonzero,

\( \begin{align} \tau-a&>0 & t-\tau-b&>0,\mbox{ or}\\ a&<\tau & \tau&<t-b \end{align} \)

This clearly indicates the limits on the integral - but it also indicates a restriction on time. That is, for both inequalities to be satisfied,

\( \begin{align} a&<t-b,\mbox{ or}\\ t&>a+b \end{align} \)

which gives rise to a unit step function that can be written outside the integral since it is not a function of \(\tau\). Specifically,

\( \begin{align} \int_{-\infty}^{\infty} u(\tau-a)~u(t-\tau-b)~f(g(\tau))d\tau&= u(t-b-a)\int_{a}^{t-b} f(g(\tau))d\tau \end{align} \)

Questions

Post your questions by editing the discussion page of this article. Edit the page, then scroll to the bottom and add a question by putting in the characters *{{Q}}, followed by your question and finally your signature (with four tildes, i.e. ~~~~). Using the {{Q}} will automatically put the page in the category of pages with questions - other editors hoping to help out can then go to that category page to see where the questions are. See the page for Template:Q for details and examples.

External Links

References