Sunday, July 19, 2009

Pole-Zero plot - Theory/Equations

Pole-Zero plot and its relation to Frequency domain:

Pole-Zero plot is an important tool, which helps us to relate the Frequency domain and Z-domain representation of a system. Understanding this relation will help in interpreting results in either domain. It also helps in determining stability of a system, given its transfer function H(z).

The Fourier transform of a sequence is given as

fourier transform

and Z-transform is given as


z-transform

There is a close relationship between these equations. If we replace z with e jw, then the z-transform reduces to the Fourier transform. When it exists, the Fourier transform is simply X(z) with z=e jw.i.e for |z| = 1, the z-transform corresponds to the Fourier transform. If instead we express z in polar form as

Frequency & z domain

This is the Fourier transform of the product of the original sequence x[n] and the exponential sequence r -n. For, r = 1, this is the Fourier transform. So, it is possible for the z-transform to converge even if the Fourier transform does not.

On a similar line, the Fourier transform and z-transform of a system can be given as

system FT

system z-transform

Visualizing Pole-Zero plot:

Since the z-transform is a function of a complex variable, it is convenient to describe and interpret it using the complex z-plane. In the z-plane, the contour corresponding to |z| = 1 is a circle of unit radius. This contour is referred to as the Unit Circle. Also, the z-transform is most useful when the infinite sum can be expressed as a simple mathematical formula. One important form of representation is to represent it as a rational function inside the Region Of Convergence.,i.e

rational TF

where, the numerator and denominator are polynomials in z. The values of z for which H(z) = 0 are called the zeros of H(z), and the values of z for which H(z) is ¥ are referred to as the poles of H(z). In other words, the zeros are the roots of the numerator polynomial and the poles of H(z) for finite values of z are the roots of the denominator polynomial. A plot of Pole and Zeros of a system on the z-plane is called a Pole-Zero plot. Usually, a Zero is represented by a 'o'(small-circle) and a pole by a 'x'(cross). Since H(z) evaluated on the unit-circle gives the frequency response of a system, it is also shown for reference in a pole-zero plot. The pole-zero GUI also uses this convention.

The pole-zero plot gives us a convenient way of visualizing the relationship between the Frequency domain and Z-domain. The frequency response H(e jw ) is obtained from the transfer function H(z), by evaluating the transfer function at specific values of z = e jw. Since, the frequency response is periodic with period 2p, we need to evaluate it over one period, such as -p < w < p. If we substitute these values of w in z= e jw, values of z lie on the unit circle and range from z = -1 all the way around and back to the point z = -1. This is shown in Figure 1 below. From this the periodicity of 2p in frequency domain corresponds to moving through an angle of 2p on the unit circle.

Complex z-plane
Figure 1

Interpreting Pole-Zero plot:

In the case of FIR filters, the location of zeros of H(z) can be used to design filters to null out specific frequencies. This can be done by placing zeros on the unit circle at locations corresponding to the frequencies where the gain needs to be 0.

In the case of IIR filters, the Pole-Zero plot can be used to determine stability. This can be done with some knowledge about the region of convergence and frequency response. For a given sequence, the set of values of z for which the z-transform converges is called the Region of Convergence (ROC). Using the definition of Z-transform given earlier in this section, the condition for convergence is as given below.

convergence in zdomain

So, as stated earlier, the z-transform might converge even if the Fourier transform does not converge. Because from the equations above, depending on the value of r the z-transform might converge while the Fourier transform might not still converge. Since, for Fourier transform to converge, we need the following condition.

convergence in zdomain

Further, the convergence is dependent only on the magnitude of z i.e |z|. Thus, if z = z1 is in the ROC, then all values of z on the circle defined by |z| = |z1| will also be in the ROC. As a result, the ROC will consist of a ring in the z-plane centered about the origin. The ROC might extend outward to ¥ or inward to origin depending on the sequence. If the ROC includes the Unit circle, this implies convergence of z-transform for |z| = 1 or, equivalently, the Fourier transform converges. The ROC cannot contain a Pole, since at a pole H(z) is infinite by definition and hence does not converge.

For a causal system (impulse response h(n) is zero for n<>), the ROC is the exterior of a circle, including ¥. Further, for a system to be stable, its impulse response must be absolutely summable. This implies the Fourier transform converges, and consequently, the ROC of H(z) must include the unit circle. Combining these two properties, a causal LTI system with a rational system function H(z) is stable if and only if all of the poles of H(z) lie inside the unit circle.i.e., the magnitude of poles must be less than 1. Thus, stability of a IIR system can be determined by looking at the pole location.


No comments:

Post a Comment