# Causality and stability of digital 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 |

As we have seen, the transfer function *H*(*Z*) is the *Z*-transform of the impulse-response signal . A causal system is one for which is one sided; that is, for . The reason is that a causal system cannot produce a response to the spike input before the spike occurs. Because the spike occurs at time zero, the response thus must be zero before time zero. Certainly, all physical systems operating in real time must be causal.

Suppose now that we take a causal system and watch its impulse-response signal as the time index *k* becomes larger and larger. If the infinite sum

**(**)

is finite, then the impulse response damps out quite rapidly, and we say that the system is stable. However, if this infinite sum diverges, we say the system is unstable. The problem of modifying the elements of an unstable system to make it stable is one of the major problems of control theory.

Let us look at two examples, which we call Example A and Example B. Let *a* be a complex constant. In Example A, the prototype causal digital signal is the causal signal given by the geometric signal

**(**)

Its *Z*-transform is the power series

**(**)

This power series converges for , which is . That is, the region of convergence is the interior of the circle of radius .

Two cases are important. Case A1 is for (Figure 6a). In such a case, the region of convergence includes the unit circle, and the causal signal is stable. Case A2 is for (Figure 6b). In that case, the region of convergence does not include the unit circle, and the causal signal is unstable.

In Example B, the prototype anticausal digital signal is the anticausal signal given by the geometric signal

**(**)

Its *Z*-transform is the Laurent series that involves only negative powers of *Z *given by

**(**)

This Laurent series converges for , which is . That is, the region of convergence is the exterior of the circle of radius .

Two cases are important. Case B1 is for (Figure 7a). In such a case, the region of convergence does not include the unit circle, and the anticausal signal is unstable. Case B2 is for (Figure 7b). In that case, the region of convergence includes the unit circle, and the anticausal signal is stable.

The above two examples show that the causal signal given by equation **50** and the anticausal signal given by equation **52** have the same expression, , for their *Z*-transforms. However, the regions of convergence are different. The region of convergence for the causal signal is inside the circle of radius , and the region of convergence for the anticausal signal is outside the circle of radius . Thus, the complete specification of the *Z*-transform requires giving the region of convergence as well as the algebraic expression.

Let us summarize. Given the *Z*-transform , where , we have two choices. If we choose the power-series representation (equation **50**), we obtain the stable causal signal

**(**)

where the entry 1 occurs at time index . Because , this signal damps geometrically in the positive time direction. On the other hand, if we choose the Laurent-series representation (equation **52**), we obtain the unstable anticausal signal

**(**)

where the first 0 occurs at time index . Because , the entries in this signal grow geometrically in the negative time direction. Thus, the choice of Laurent series leads not only to an anticausal digital signal but also to an unstable one. Clearly, such a choice is to be avoided.

Given the *Z*-transform , where , again we have two choices. If we choose the power-series representation (equation **50**), we obtain the unstable causal signal

**(**)

where the entry 1 occurs at time index . Because , this signal grows geometrically in the positive time direction. On the other hand, if we choose the Laurent-series representation (equation **52**), we obtain the stable anticausal signal

**(**)

where the first 0 occurs at time index . Because , this signal damps geometrically in the negative time direction. Thus we have no clear-cut choice, as we did when . We cannot have both causality and stability; we must choose one or the other. The usual choice is equation **56**, which gives us stability at the expense of causality.

We now wish to present the method of partial fractions to find the impulse response of a digital ARMA system from its transfer function

**(**)

The right-hand expression is the factored form of the transfer function. The are its poles and the are its zeros. For simplicity of presentation, we will assume that , so that is a proper fraction in *Z*. In addition, we will suppose that all the poles , ,..., are simple ones. We now write as the sum of partial fractions

**(**)

The constants are determined by means of the equation

**(**)

We now divide the poles into two sets: those with magnitude greater than one, which we designate as , , , and those with magnitude less than one, which we designate as , , ,. For the first set, we pick the regions as the regions of convergence. For the second set, we pick the regions as the regions of convergence. Thus, all the regions include the unit circle. It follows that the stable impulse response is the two-sided signal given by

**(**)

Unless otherwise stated, it is assumed that an ARMA system is both stable and causal. The above expression for the impulse response shows that the system is causal if and only if *the set of poles of magnitude less than one is empty*. Thus we have proved the following theorem: All the poles , ,..., of a stable causal digital ARMA system have magnitude greater than one. In the case in which and in which the poles are distinct, the stable causal impulse response has the form

**(**)

where the coefficients are given by equation **59**.

## Continue reading

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

Analog transfer functions | Causality and stability of analog systems |

Previous chapter | Next chapter |

Absorption | none |

## Also in this chapter

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