Frequency response of a digital system
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 |
We recall that impulse response is the system’s response to a unit spike. Thus, the impulse response represents a transient that dies out over time for a stable system. On the other hand, the frequency response can be interpreted as a steady-state response of the system to signals, each with a pure frequency. Let the real variable denote frequency (in radians per second), and let be the discrete time spacing (in seconds). The function is a steady-state sinusoidal wave of pure frequency If we let this steady-state signal be the input to a unit-delay system Z, the steady-state output is the delayed signal . The frequency response is the ratio of the steady-state output to the steady-state input ; that is,
( )
Thus, the pure-delay operator Z can be represented as the point on the unit circle. More generally, if we let the steady-state input be the input to the system
( )
the steady-state output is
( )
The frequency response is defined as the ratio of the steady-state output to the steady-state input; that is, the frequency response is
( )
The frequency response written thus conforms to the standard electrical engineering convention, with negative signs in the exponents. Therefore, although we use the geophysics Z-transform (i.e., the generating function) instead of the electrical engineering z-transform, we use the same convention regarding frequency that electrical engineers use. The transfer function
( )
of a stable digital system exists in a region that includes the unit circle. If we substitute into the transfer function , we obtain the frequency response . That is, the values of H(Z) on the unit circle make up the frequency response of the system. Instead of the cumbersome notation , the frequency response is written more conveniently as . Thus, we have
( )
The frequency response is in the form of a discrete Fourier transform. The magnitude of the frequency response is called the magnitude spectrum, and the angle is called the phase spectrum — both of which we first encountered in Chapter 6. The phase-lag spectrum is defined as the negative of the phase spectrum; that is, . The slope of the phase-lag spectrum is called the group delay; that is, the group delay is . Because
( )
it follows that the frequency response is periodic in , with a period of . As a result, we have to examine only the frequency response over one period, which usually is taken from the point to the point . The frequency is known as the Nyquist frequency, which we discussed in Chapter 4. For a causal system, the frequency response is the one-sided discrete Fourier transform
( )
The frequency response of an ARMA system is
( )
which in factored form is
( )
An all-pass system is defined as a stable causal system with unit magnitude spectrum; that is,
( )
All-pass systems can be broken down into four types: trivial, pure-delay, dispersive, and impure-delay systems. The transfer function of a trivial all-pass system is merely a constant of magnitude one. In the case of real systems, this constant must be either +1 or -1, so a (real) trivial all-pass system either produces no phase shift or a 180° phase shift for all frequencies. A pure-delay all-pass system delays a signal by a given positive constant N; that is, it transforms the input into . As a result, the transfer function of a pure-delay all-pass system is
( )
so the frequency response is
( )
Thus, a pure-delay all-pass system produces a phase lag of
( )
which is a linear function of . A dispersive all-pass system is an ARMA system of the form
( )
In other words, an all-pass system has a special arrangement of poles and zeros: The poles and the zeros are symmetric with respect to the unit circle. If we write , then the pole is and the zero is . Both are on the same radial line, with the pole outside the unit circle and the zero inside the unit circle.
Therefore, the frequency spectrum is the product of q factors of the form
( )
The frequency spectrum 82 has a magnitude of unity because the complex conjugate of the numerator is equal to times the denominator. Hence, the numerator has the same magnitude as the denominator.
An impure-delay all-pass system never is used in engineering applications, so we will omit its consideration here.
Let us say a word about notation. We use and as the transfer function and the impulse response, respectively, of an all-pass system. This use of p should be kept distinct from the constant p appearing in the designation of an ARMA system.
The group-delay theorem (Bode, 1945^{[1]}; Robinson, 1963^{[2]}, theorem 5, p. 204) deals with group delay of a digital all-pass system. The group delay of a nontrivial digital all-pass system is positive for all frequencies.
We can prove the group-delay theorem as follows: Often it is convenient to change our time scale from a scale in seconds to a scale in which the discrete time-spacing is one unit. Then we can let so . The Nyquist frequency becomes simply , and the Nyquist range becomes . The pure-delay all-pass system, as we have seen, has phase lag where time delay , so it satisfies the theorem. The dispersive all-pass system is made up of q factors of the form . In terms of frequency, this factor becomes
( )
The constant increases the phase lag by an additive constant, so we can neglect it as far as this proof is concerned. The phase lag of the factor is the difference of the angles of the numerator and denominator; that is,
( )
In Figure 10, we show the angle for various points on the unit circle. We trace out one period of the frequency function if we travel counterclockwise around the unit circle from point A to B, C, D, E, and back to A. Over this journey, we see that the angle increases monotonically from 0 to . Thus, the given factor of a dispersive all-pass system produces an increase in phase lag of over one complete period of frequency Therefore, a dispersive all-pass system of q such factors would have a phase-lag spectrum that increases monotonically by an amount over one complete period; that is,
( )
Because the phase-lag spectrum of an all-pass system is increasing monotonically, it follows that its derivative is positive. This derivative is the group delay. Thus, a nontrivial all-pass system has positive group delay for all frequencies .
By definition, an all-pass system is a causal system with the property that the magnitude spectrum of the output signal is the same as the magnitude spectrum of the input signal. The energy-delay theorem (Robinson, 1963, theorem 7, p. 204) establishes that an all-pass system has the property that the time structure of output energy is delayed with respect to the time structure of input energy. In other words, an all-pass system differentially delays the return of energy in passing a signal from input to output. An all-pass system is comparable to a bureaucrat who day by day receives written requests on his desk and who holds them and later passes them on in some other order to a subordinate. Nothing in content has changed, but certain requests are delayed with respect to others. In seismic processing, an all-pass system always gives trouble.
The energy-delay theorem has the formal statement: Let a stable causal signal be the input to a nontrivial digital all-pass system , and let the stable causal signal be the resulting output. The front-end energy of the input is greater than or equal to the front-end energy of the output for all time indices n; that is,
( )
The strict inequality (i.e., greater than) holds for some . For , the equality holds (i.e., the energy of the input is equal to the energy of the output).
We can prove the energy-delay theorem as follows: Bessel’s equality gives
( )
Because , the above equation becomes
( )
which says that output and input have the same total energy. The truncated input (for time index n) is defined as
( )
In other words, the truncated input is the front end of the input signal. Let designate the output of the all-pass system to the truncated input . Two statements can be made immediately about the output : (1) because the system is all pass, the total energy of the truncated input is equal to the total energy of ; that is,
( )
and (2) because the system is causal, the front-end values of and are the same; that is,
( )
Putting these two things together, we have
( )
Because the all-pass system is nontrivial, it can be shown that
( )
for some n. Hence, it follows that the front-end energy of the signal exceeds the front-end energy of the all-pass filtered signal . In other words, an all-pass system delays the energy from input to output.
Let us consider a numerical example. The all-pass filter is , with impulse response (Table 1, first row below the column headings). The input signal to the all-pass filter is the wavelet (Table 1, second row) with three coefficients. The resulting output signal is the wavelet (Table 1, third row). For the cutoff n = 0, the front-end signal (i.e., the front end of the input ) is for other k (Table 1, fourth row). In other words, the front-end signal has only one nonzero coefficient, and it occurs at time 0. The resulting output of the all-pass system to the front-end signal is the signal (Table 1, fifth row). Note that the front ends of both the inputs u and f are the same (namely, –0.75), and the front ends of both the outputs y and g are the same (namely, –0.375).
Time k | 0 | 1 | 2 | 3 | 4 |
---|---|---|---|---|---|
All-pass filter | 0.500 | 0.750 | -0.375 | 0.188 | -0.094 |
Input signal | -0.750 | 2.750 | 1.000 | 0 | 0 |
Output for input signal | -0.375 | 0.813 | 2.844 | -0.422 | 0.211 |
Front end of input signal | –0.750 | 0 | 0 | 0 | 0 |
Output for front end | -0.375 | -0.563 | 0.282 | -0.141 | 0.070 |
Time k | 0 | 1 | 2 | 3 | 4 |
---|---|---|---|---|---|
Front end of | –0.375 | 0 | 0 | 0 | 0 |
Tailgate of | 0 | –0.563 | 0.282 | –0.141 | 0.070 |
Front-end energy of | 0 | 0 | 0 | 0 | |
Tailgate energy of | 0 |
The output is separated into two parts: (1) the front end (Table 2, first row below the column headings) and (2) the tailgate (Table 2, second row extended to infinity). The front-end energy of (Table 2, third row) is . The tailgate energy of is the sum of squares of its elements (Table 2, fourth row extended to infinity), which is 0.4219. Because the filter is an all-pass filter, the sum of front-end energy of and the tailgate energy of , namely 0.1406 + 0.4219 = 0.5625, is equal to the input front-end energy, namely . Thus, the front-end energy of the input to the all-pass filter exceeds the front-end energy of the output. This demonstration is shown graphically in Figure 11.
The concept of minimum phase is attributable to Bode (1945). Let be the poles and let be the zeros of a system. A minimum-phase (or minimum-delay) system is a stable causal system that has its zeros as well as its poles lying outside of the unit circle. That is, stability requires that , whereas minimum delay requires in addition that . The canonical representation theorem (Robinson, 1963, theorem 1, p. 203) can be stated as follows:
The signal is stable and causal if and only if
( )
where is the minimum-delay signal (with the same magnitude spectrum as that for ) and is an all-pass signal. In other words, a minimum-delay system is comparable to an executive who day by day receives written requests, acts on them as soon as possible, and immediately passes the results to a bureaucrat (an all-pass system). The bureaucrat holds the results, changes nothing, and then passes them on with various delays.
The inverse of the system is defined as the system 1/H(Z). Thus, the inverse of the ARMA(p,q) system B(Z)/A(Z) is the ARMA(q,p) system A(Z)/B(Z). The poles and zeros of the original system become, respectively, the zeros and poles of the inverse system. If is the impulse response for H(Z), then we write for the impulse response for . However, is not , even as is not .
Because an AR(p) system is an all-pole system, it follows that a stable causal AR(p) system is necessarily a minimum-delay system. The inverse of a stable AR(p) system is an MA(p) system, in which the poles of the AR(p) system become the zeros of the MA(p) system. Thus we see that the inverse of an AR(p) system is a minimum-delay MA(p) system.
In our first such example, the given system is the stable causal AR(1) system , where . It has the stable causal impulse response for k = 0, 1, 2, …. The inverse system is the minimum-delay MA(1) system , with the impulse response .
In a second example, the given system is the minimum-delay ARMA(2,1) system
( )
so the impulse response is the minimum-delay signal
( )
The inverse system is the minimum-delay ARMA(1,2) system
( )
so the impulse response of the inverse system is the causal signal
( )
In conclusion, let us state the (digital) minimum-delay theorem (Robinson, 1963, theorem 10, p. 205). Let the signal be a stable causal signal in the class of all stable causal signals with the same magnitude spectrum. Then each of the following conditions is necessary and sufficient for to be minimum delay:
1) The group delay of is minimum for all frequencies.
2) The front-end energy is a maximum for all time indices .
3) The inverse system is causal.
References
- ↑ Bode, H. W., 1945, Network analysis and feedback amplifier design: Van Nostrand.
- ↑ Robinson, E. A., 1963, Nekotorye svoystva razlozheniya vol’da statsionarnykh sluchaynykh protsessov [Properties of the Wold Decomposition of stationary stochastic processes; in Russian]: Teoriya Veroyatnostei i ee Primememiya Akademiya Nauk SSSR, 7, no. 2, 201–211.
Continue reading
Previous section | Next section |
---|---|
Causality and stability of analog systems | Digital prediction |
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 digital systems
- Causality and stability of analog systems
- Digital prediction
- Digital prediction error
- Analog prediction error