28.02.2013 Views

Introduction to Acoustics

Introduction to Acoustics

Introduction to Acoustics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

524 Part D Hearing and Signal Processing<br />

Part D 14.12<br />

limits,<br />

X(ω)<br />

TD/2 �<br />

dt<br />

TD<br />

−TD/2<br />

′ e −iωt′<br />

= 1<br />

×<br />

(N−1)/2<br />

�<br />

n=−(N−1)/2<br />

14.12 The z-Transform<br />

inωot ′<br />

X(nωo) e e −iωTD/2 inωoTD/2<br />

e ,<br />

(14.138)<br />

Like the discrete Fourier transform, the z-transform is<br />

well suited <strong>to</strong> describing sampled signals. We consider<br />

x(t) <strong>to</strong> be a sampled signal so that it is defined at discrete<br />

time points t = tk = kTs,whereTs is the sampling period.<br />

Then the time dependence of x can be described by an<br />

index, xk = x(tk). The z-transform of x is<br />

X(z) =<br />

∞�<br />

k=−∞<br />

xkz −k . (14.140)<br />

The quantity z is complex, with amplitude A and<br />

phase ϕ,<br />

z = A e iϕ = A e iωTs , (14.141)<br />

where ϕ is the phase advance in radians per sample.<br />

In the special case where A = 1, all values of z lie<br />

on a circle of radius 1 (the unit circle) inthecomplex<br />

z plane. In that case the z-transform is equivalent <strong>to</strong> the<br />

discrete Fourier transform. An often-overlooked alternative<br />

view is that the z-transform is an extension of<br />

the Fourier transform wherein the angular frequency ω<br />

becomes complex,<br />

so that<br />

ω = ωR + iωI , (14.142)<br />

z = e −ωITs e iωRTs . (14.143)<br />

The extended Fourier transform will not be pursued<br />

further in this chapter.<br />

A well-defined z-transform naturally includes a function<br />

of variable z, but the function itself is not enough.<br />

In order for the inverse transform <strong>to</strong> be unique, the definition<br />

also requires that the region of the complex plane<br />

in which the transform converges must also be specified.<br />

which reduces <strong>to</strong><br />

X(ω) =<br />

(N−1)/2<br />

�<br />

sin[(ω − nωo)TD/2]<br />

X(nωo)<br />

(ω − nωo)TD/2<br />

n=−(N−1)/2<br />

×e −iωTD/2 iπn<br />

e . (14.139)<br />

Table 14.2 z-Transform pairs<br />

xk X(z) Radius of<br />

convergence<br />

δk,k 0 z −k 0 all z<br />

a k uk z/(z − a) |z| > a<br />

ka k uk az/(z − a) 2 |z| > a<br />

a k cos(ωoTsk)uk<br />

a k sin(ωoTsk)uk<br />

z 2 −az cos(ωoTs)<br />

z 2 −2az cos(ωoTs)+a 2<br />

az sin(ωoTs)<br />

z 2 −2az cos(ωoTs)+a 2<br />

|z| > a<br />

|z| > a<br />

To illustrate that point, one can consider two different<br />

functions xk that have the same z-transform function,<br />

but different regions of convergence.<br />

Consider first the function<br />

xk = 2 k<br />

xk = 0 for k < 0 .<br />

for k ≥ 0 , (14.144)<br />

This two-line function can be written as a single line by<br />

using the discrete Heaviside function uk. The function<br />

uk is defined as zero when k is a negative integer and as<br />

+1whenk is any other integer, including zero. Then xk<br />

above becomes<br />

xk = 2 k uk . (14.145)<br />

The z-transform of xk is<br />

X(z) =<br />

∞�<br />

(2/z) k . (14.146)<br />

k=0<br />

The sum is a geometric series, which converges <strong>to</strong><br />

X(z) =<br />

1<br />

= z/(z − 2) (14.147)<br />

1 − 2/z

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!