# Digital linear time-invariant systems

Series | Geophysical References Series |
---|---|

Title | Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing |

Author | Enders A. Robinson and Sven Treitel |

Chapter | 15 |

DOI | http://dx.doi.org/10.1190/1.9781560801610 |

ISBN | 9781560801481 |

Store | SEG Online Store |

A discrete-time signal is a signal of (real or complex) numbers defined for every integer *k*. The index *k* represents the set of discrete, equally spaced time points. A discrete time signal often is called a *digital signal*. It is customary to denote time spacing between the discrete time points as . Actual time *t* is related to the index *k* by the equation . Usually, time is measured in seconds.

An important type of digital signal is the *impulse function*. The discrete impulse function is defined as the signal made up entirely of zeros except for a value of one at time index . That is, is given by

**(**)

The impulse function is also called the *Kronecker delta function*. A more popular name for the delta function is the *unit spike*.

A digital system is a discrete time system that can be represented by a rule that transforms a signal into another signal . The signal is called the *input*, and the signal is the *output*. This relationship is indicated by the notation

**(**)

where *S* denotes the digital system. A system *L* is linear if

**(**)

for any constants *a* and *b* and for any digital signals and . In particular, if is the output of a linear system for input , then is the output for input . That is, if we amplify the input of a linear system by a constant factor, the output also is amplified by this same constant factor. In addition, if and are the respective outputs of a linear system for the inputs and , then is the output for input . That is, if we add the inputs of a linear system, the outputs also are added.

A digital system is said to be *time invariant* if for any integer *n*. That is, a system is time invariant if a shift in the input produces the same shift in the output.

Let us give some examples. We begin with the delay line , where *n* is a constant. The delay line is linear and time invariant. The rectifier is nonlinear and time invariant. The amplifier is linear and *time varying*. It is time varying because (1) a shift in the input produces the output a shift in the output produces ; and (3) the two outputs are not the same.

A signal is *one-sided*, or *causal*, if for . A signal is *anticausal* if for . Finally, a signal is *noncausal* if it has an anticausal component.

A system is said to be causal if a causal input yields a causal output. That is, a causal system is a system with this property: If for , then for . If the system represents a physical phenomenon operating in real time, then the system must be causal. However, if the system represents an analysis of past recordings (such as a seismic record) so that *k* represents nominal time (i.e., time marks on the record), then the system need not be causal.

A linear time-invariant system is characterized by its impulse response. The *impulse response* is defined as the output resulting from a spike input ; that is,

**(**)

Because the impulse is causal, we see that a causal system has a causal impulse response, whereas a noncausal system has a noncausal impulse response. Let us now express, in terms of , the output of a linear time-invariant system that has an arbitrary input . Because of time invariance, the output resulting from the input is for any *n*; that is,

**(**)

Because of linearity, the output resulting from the input is . The input can be written as

**(**)

The output resulting from this input is

**(**)

which is

**(**)

The expression on the right is a *convolution*. It represents the convolution of with and is denoted by . Convolution is commutative, so we can write the input-output relationship as

**(**)

The most important class of linear time-invariant digital systems is the class represented by a *finite-difference equation* (with constant coefficients) of the form

**(**)

This difference equation represents an ARMA(*p,q*) system. In the case of *q* = 0 and , it reduces to the AR(*p*) system

**(**)

On the other hand, if *p* = 0, the ARMA(*p,q*) system reduces to the MA(*q*) system

**(**)

Most of our work will deal with the three difference equations **10**, **11**, and **12**. Let us now make some remarks about these three equations. In the form in which they are written above, we want index *k* to represent the present instant of time, so would be the present input and the present output. As a result, the other input values , represent past values, and the other output values represent past values. Thus, the above difference equations, so interpreted, involve only present and past values of the input and output and no future values. Under this interpretation, the systems represented by equations **10**, **11**, and **12** are causal — that is, they can be implemented in real time. Conceptually, the MA system is the easiest one to visualize.

We now want to introduce the *backward-shift* (or *unit-delay*) *operator Z*, which is defined as for any signal . In terms of the backward-shift operator, the MA(*q*) equation is

**(**)

or

**(**)

A diagram for this system is shown in Figure 1. In such diagrams, a branch point is indicated by a solid circle. At a branch point, the signal that goes out on each branch is the same as the signal that enters the branch point. A small open circle represents a summation point.

At a summation point, all the incoming signals are added to produce the outgoing signal.

In Figure 1, the forward direction (i.e., the direction from input to output) is from left to right. Because all of the arrows associated with the constant multipliers , are in the forward direction, an MA system is seen to be a feedforward system. In a feedforward system, the present output depends on the present and past values , , , ... of the input. In the case of an MA system, only a finite number *q* of past values are involved.

Next, let us construct a diagram for an AR(*p*) system, which we can write as

**(**)

A diagram for this system is shown in Figure 2. Because all the arrows associated with the multipliers , are in the backward direction, an AR system is seen to be a (pure) feedback system. In a feedback system, the present input can be expressed in terms of the present and past values , , , ... of the output. In the case of an AR system, only a finite number *p*of past values are involved.

Finally, let us construct a diagram of the ARMA(*p*, *q*) system

**(**)

An ARMA(*p*, *q*) system can be regarded as the combination in series of an AR(*p*) system and an MA(*q*) system. A diagram for an ARMA system is shown in Figure 3.

## Continue reading

Previous section | Next section |
---|---|

Introduction | Analog linear time-invariant systems |

Previous chapter | Next chapter |

Absorption | none |

## Also in this chapter

- Introduction
- Analog linear time-invariant systems
- Digital transfer functions
- Analog transfer functions
- Causality and stability of digital systems
- Causality and stability of analog systems
- Frequency response of a digital system
- Digital prediction
- Digital prediction error
- Analog prediction error