2D Z-transform

Last updated

The 2D Z-transform, similar to the Z-transform, is used in multidimensional signal processing to relate a two-dimensional discrete-time signal to the complex frequency domain in which the 2D surface in 4D space that the Fourier transform lies on is known as the unit surface or unit bicircle. [1] The 2D Z-transform is defined by

Contents

where are integers and are represented by the complex numbers:

The 2D Z-transform is a generalized version of the 2D Fourier transform. It converges for a much wider class of sequences, and is a helpful tool in allowing one to draw conclusions on system characteristics such as BIBO stability. It is also used to determine the connection between the input and output of a linear shift-invariant system, such as manipulating a difference equation to determine the system's transfer function.

Region of Convergence (ROC)

The Region of Convergence is the set of points in complex space where:

In the 1D case this is represented by an annulus, and the 2D representation of an annulus is known as the Reinhardt domain. [2] From this one can conclude that only the magnitude and not the phase of a point at will determine whether or not it lies within the ROC. In order for a 2D Z-transform to fully define the system in which it means to describe, the associated ROC must also be know. Conclusions can be drawn on the Region of Convergence based on Region of Support of the original sequence .

Finite-support sequences

A sequence with a region of support that is bounded by an area within the plane can be represented in the z-domain as:

Because the bounds on the summation are finite, as long as z1 and z2 are finite, the 2D Z-transform will converge for all values of z1 and z2, except in some cases where z1 = 0 or z2 = 0 depending on .

First-quadrant and wedge sequences

Sequences with a region of support in the first quadrant of the plane have the following 2D Z-transform:

From the transform if a point lies within the ROC then any point with a magnitude

also lie within the ROC. Due to these condition, the boundary of the ROC must have a negative slope or a slope of 0. This can be assumed because if the slope was positive there would be points that meet the previous condition, but also lie outside the ROC. [2] For example, the sequence:

has the transform

It is obvious that this only converges for

So the boundary of the ROC is simply a line with a slope of -1 in the plane. [2]

In the case of a wedge sequence where the region of support is less than that of a half plane. Suppose such a sequence has a region of support over the first quadrant and the region in the second quadrant where . If is defined as the new 2D Z-Transform becomes:

Sequence with Region of support over a wedge and its corresponding ROC Wedge ROC.jpg
Sequence with Region of support over a wedge and its corresponding ROC

This converges if:

These conditions can then be used to determine constraints on the slope of the boundary of the ROC in a similar manner to that of a first quadrant sequence. [2] By doing this one gets:

and

Sequences with region of support in all quadrants

A sequence with an unbounded Region of Support can have an ROC in any shape, and must be determined based on the sequence . A few examples are listed below:

will converge for all . While:

will not converge for any value of . However, These are the extreme cases, and usually, the Z-transform will converge over a finite area. [2]

A sequence with support over the entire can be written as a sum of each quadrant sequence:

Now suppose:

and also have similar definitions over their respective quadrants. Then the Region of convergence is simply the intersection between the four 2D Z-transforms in each quadrant.

Using the 2D Z-transform to solve difference equations

A 2D difference equation relates the input to the output of a Linear Shift-Invariant (LSI) System in the following manner:

Due to the finite limits of computation, it can be assumed that both a and b are sequences of finite extent. After using the z transform, the equation becomes:

This gives:

Thus we have defined the relation between the input and output of the LSI system.

Using the 2D Z-transform to determine stability

Shanks' Theorem I

For a first quadrant recursive filter in which . The filter is stable iff: [3]

for all points such that or .

Shanks' Theorem II

For a first quadrant recursive filter in which . The filter is stable iff: [3]

Huang's Theorem

For a first quadrant recursive filter in which . The filter is stable iff: [3]

for any such that

Decarlo and Strintzis' Theorem

For a first quadrant recursive filter in which . The filter is stable iff: [3]

for any such that

for any such that

Calculation of 2D Z-transforms

Approach 1: Finite sequences

For finite sequences, the 2D Z-transform is simply the sum of magnitude of each point multiplied by raised to the inverse power of the location of the corresponding point. For example, the sequence:

has the Z-transform:

As this is a finite sequence the ROC is for all .

Approach 2: Sequences with values along only or

For a sequence with a region of support on only or , the sequence can be treated as a 1D signal and the 1D Z-transform can be used to solve for the 2D Z-transform. For example, the sequence:

Is clearly given by .

Therefore, its Z-transform is given by:

As this is a finite sequence the ROC is for all .

Approach 3: Separable sequences

A separable sequence is defined as

For a separable sequence, finding the 2D Z-transform is as simple as separating the sequence and taking the product of the 1D Z-transform of each signal and . For example, consider the sequence

.

Its Z-transform is given by

.

The ROC is given by

 ; .

Related Research Articles

In complex analysis, an entire function, also called an integral function, is a complex-valued function that is holomorphic on the whole complex plane. Typical examples of entire functions are polynomials and the exponential function, and any finite sums, products and compositions of these, such as the trigonometric functions sine and cosine and their hyperbolic counterparts sinh and cosh, as well as derivatives and integrals of entire functions such as the error function. If an entire function has a root at , then , taking the limit value at , is an entire function. On the other hand, the natural logarithm, the reciprocal function, and the square root are all not entire functions, nor can they be continued analytically to an entire function.

In mathematics, a series is, roughly speaking, an addition of infinitely many terms, one after the other. The study of series is a major part of calculus and its generalization, mathematical analysis. Series are used in most areas of mathematics, even for studying finite structures in combinatorics through generating functions. The mathematical properties of infinite series make them widely applicable in other quantitative disciplines such as physics, computer science, statistics and finance.

<span class="mw-page-title-main">Natural logarithm</span> Logarithm to the base of the mathematical constant e

The natural logarithm of a number is its logarithm to the base of the mathematical constant e, which is an irrational and transcendental number approximately equal to 2.718281828459. The natural logarithm of x is generally written as ln x, logex, or sometimes, if the base e is implicit, simply log x. Parentheses are sometimes added for clarity, giving ln(x), loge(x), or log(x). This is done particularly when the argument to the logarithm is not a single symbol, so as to prevent ambiguity.

<span class="mw-page-title-main">Fourier series</span> Decomposition of periodic functions into sums of simpler sinusoidal forms

A Fourier series is an expansion of a periodic function into a sum of trigonometric functions. The Fourier series is an example of a trigonometric series. By expressing a function as a sum of sines and cosines, many problems involving the function become easier to analyze because trigonometric functions are well understood. For example, Fourier series were first used by Joseph Fourier to find solutions to the heat equation. This application is possible because the derivatives of trigonometric functions fall into simple patterns. Fourier series cannot be used to approximate arbitrary functions, because most functions have infinitely many terms in their Fourier series, and the series do not always converge. Well-behaved functions, for example smooth functions, have Fourier series that converge to the original function. The coefficients of the Fourier series are determined by integrals of the function multiplied by trigonometric functions, described in Fourier series§Definition.

In mathematics, an infinite series of numbers is said to converge absolutely if the sum of the absolute values of the summands is finite. More precisely, a real or complex series is said to converge absolutely if for some real number Similarly, an improper integral of a function, is said to converge absolutely if the integral of the absolute value of the integrand is finite—that is, if A convergent series that is not absolutely convergent is called conditionally convergent.

<span class="mw-page-title-main">Stirling's approximation</span> Approximation for factorials

In mathematics, Stirling's approximation is an asymptotic approximation for factorials. It is a good approximation, leading to accurate results even for small values of . It is named after James Stirling, though a related but less precise result was first stated by Abraham de Moivre.

In mathematics and signal processing, the Z-transform converts a discrete-time signal, which is a sequence of real or complex numbers, into a complex valued frequency-domain representation.

<span class="mw-page-title-main">Inverse trigonometric functions</span> Inverse functions of sin, cos, tan, etc.

In mathematics, the inverse trigonometric functions are the inverse functions of the trigonometric functions, under suitably restricted domains. Specifically, they are the inverses of the sine, cosine, tangent, cotangent, secant, and cosecant functions, and are used to obtain an angle from any of the angle's trigonometric ratios. Inverse trigonometric functions are widely used in engineering, navigation, physics, and geometry.

In mathematics, a Dirichlet series is any series of the form where s is complex, and is a complex sequence. It is a special case of general Dirichlet series.

<span class="mw-page-title-main">Digamma function</span> Mathematical function

In mathematics, the digamma function is defined as the logarithmic derivative of the gamma function:

In mathematics, the ratio test is a test for the convergence of a series

<span class="mw-page-title-main">Dirichlet eta function</span> Function in analytic number theory

In mathematics, in the area of analytic number theory, the Dirichlet eta function is defined by the following Dirichlet series, which converges for any complex number having real part > 0:

In mathematics, an asymptotic expansion, asymptotic series or Poincaré expansion is a formal series of functions which has the property that truncating the series after a finite number of terms provides an approximation to a given function as the argument of the function tends towards a particular, often infinite, point. Investigations by Dingle (1973) revealed that the divergent part of an asymptotic expansion is latently meaningful, i.e. contains information about the exact value of the expanded function.

<span class="mw-page-title-main">Series expansion</span> Expression of a function as an infinite sum of simpler functions

In mathematics, a series expansion is a technique that expresses a function as an infinite sum, or series, of simpler functions. It is a method for calculating a function that cannot be expressed by just elementary operators.

In mathematics, the Riemann series theorem, also called the Riemann rearrangement theorem, named after 19th-century German mathematician Bernhard Riemann, says that if an infinite series of real numbers is conditionally convergent, then its terms can be arranged in a permutation so that the new series converges to an arbitrary real number, and rearranged such that the new series diverges. This implies that a series of real numbers is absolutely convergent if and only if it is unconditionally convergent.

<span class="mw-page-title-main">Mercator series</span> Taylor series for the natural logarithm

In mathematics, the Mercator series or Newton–Mercator series is the Taylor series for the natural logarithm:

In mathematics, convergence tests are methods of testing for the convergence, conditional convergence, absolute convergence, interval of convergence or divergence of an infinite series .

This is a summary of differentiation rules, that is, rules for computing the derivative of a function in calculus.

In mathematical analysis and applications, multidimensional transforms are used to analyze the frequency content of signals in a domain of two or more dimensions.

References

  1. Siamak Khatibi, “Multidimensional Signal Processing: Lecture 11”, BLEKINGE INSTITUTE OF TECHNOLOGY, PowerPoint Presentation.
  2. 1 2 3 4 5 Dan E. Dudgeon, Russell M. Mersereau, “Multidimensional Digital Signal Processing”, Prentice-Hall Signal Processing Series, ISBN   0136049591, 1983.
  3. 1 2 3 4 Ed. Alexander D. Poularikas, “The Handbook of Formulas and Tables for Signal Processing”, Boca Raton: CRC Press LLC, 1999.