Funciones digitales de transferencia

ADVERTISEMENT
From SEG Wiki
Jump to: navigation, search
This page is a translated version of the page Digital transfer functions and the translation is 51% complete.

Other languages:
English • ‎español
Digital Imaging and Deconvolution: The ABCs of Seismic Exploration and Processing
DigitalImaging.png
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

In one key place, the conventions used in signal processing by geophysicists differ from those used by electrical engineers. At the outset, it is good to point out the differences so as to make the literature more accessible. Geophysicists and electrical engineers have different conventions with respect to the z-transform. Let , , , ... be the impulse response of a causal time-invariant linear filter. The engineering z-transform is


(24)

whereas the geophysics Z-transform is


(25)

To distinguish between the two, we use a capital Z in the geophysics Z-transform and a lowercase z in the engineering z-transform. Throughout this book, we use mostly the geophysics Z-transform. The two transforms are related thusly: . Whereas the engineering z represents a unit advance operator, the geophysics Z represents a unit-delay operator. The common mathematical terminology for the geophysics Z-transform is generating function.

The transfer function (or system function) of a digital system is defined as the Z-transform of the impulse response. If is the impulse response, the transfer function is


(26)
Figure 4.  Annular region of convergence of a Laurent series.

The Z-transform establishes a correspondence between the signal and the transfer function H(Z). This expression for H(Z) is in the form of a Laurent series (i.e., a power series that can involve negative as well as positive powers of Z). As we know from the theory of functions, the region of convergence of the Laurent series is an annular region (i.e., a ring-shaped region) , whose inner radius and outer radius depend on the behavior of the signal as and as , respectively. Within this ring, the Laurent series has no poles or other singularities (Figure 4). The Z-transform given above is called a two-sided Z-transform because it involves a signal with both negative and positive values of the time index k.

Causal signals are zero for negative values of k, and for such signals, we can use the one-sided Z-transform


(27)

In this case, the radius of the inner ring is zero, so the region of convergence is the interior of a circle. The one-sided Z-transform is in the form of a series that involves only nonnegative powers of Z. The symbol Z denotes the parameter in the Z-transform, whereas the symbol Z denotes the one-sided Z-transform itself.

The digital ARMA(p, q) model


(28)

represents a recursion relationship between the input signal and the output signal . If the input signal is causal, then the output signal also must be causal. The recursion relations are


(29)

Given the input signal, the first equation determines . The next equation can be used to find . The third equation then can be solved for . By continuing in this way, equation by equation, we can find all values of the output signal. We recall that the impulse response is the output resulting from a spike input. Thus we let the input be the Kronecker delta signal . The recursion relations (in the case in which p > q) are


(30)

Solving equation by equation, we can obtain the impulse response .

Alternatively, let us multiply the first equation by , the second by , the third by , and so on. We obtain


(31)

We now add these equations, column by column, to obtain


(32)

where


(33)

We recognize H(Z) and B(Z) as the Z-transforms of the impulse response and the feedforward coefficients , respectively. If we factor out H(Z), we have


(34)

We recognize the expression in parentheses as the Z-transform A(Z) of the feedback coefficients ; that is,


(35)

Thus, we have


(36)

The function H(Z), which is the Z-transform of the impulse response, is the transfer function. The Z-transforms A(Z) and B(Z) are polynomials. We can write H(Z) in factored form as


(37)

The constants , , ..., are the poles, and the constants , , ..., are the zeros of H(Z). In the case in which the and coefficients are real, all of the complex and values occur in complex-conjugate pairs. Because an AR(p) system has only poles, it also is called an all-pole system. Likewise, an MA(q) system is called an all-zero system.


Sigue leyendo

Sección previa Siguiente sección
Sistemas analógicos lineales invariantes en el tiempo Funciones análogicas de transferencia
Capítulo previo Siguiente capítulo
Absorción nada

Tabla de contenido


También en este capítulo


Vínculos externos

find literature about
Digital transfer functions/es
SEG button search.png Datapages button.png GeoScienceWorld button.png OnePetro button.png Schlumberger button.png Google button.png AGI button.png