Information about Transfer Matrix

The transfer matrix is a formulation in terms of a block-Toeplitz matrix of the two-scale equation, which characterizes refinable functions. Refinable functions play an important role in wavelet theory and finite element theory.

For the mask , which is a vector with component indexes from to , the transfer matrix of , we call it here, is defined as
.
More verbosely
The effect of can be expressed in terms of the downsampling operator "":
.

Properties

  • .
  • If you drop the first and the last column and move the odd indexed columns to the left and the even indexed columns to the right, then you obtain a transposed Sylvester matrix.
  • The determinant of a transfer matrix is essentially a resultant.
More precisely:
Let be the even indexed coefficients of () and let be the odd indexed coefficients of ().
Then , where is the resultant.
This connection allows for fast computation using the Euclidean algorithm.
  • For the determinant of the transfer matrix of convolved mask holds
where denotes the mask with alternating signs, i.e. .
  • If , then .
This is a concretion of the determinant property above. From the determinant property one knows that is singular whenever is singular. This property also tells, how vectors from the null space of can be converted to null space vectors of .
  • If is an eigenvector of with respect to the eigenvalue , i.e.
,
then is an eigenvector of with respect to the same eigenvalue, i.e.
.
  • Let be the eigenvalues of , which implies and more generally . This sum is useful for estimating the spectral radius of . There is an alternative possibility for computing the sum of eigenvalue powers, which is faster for small .
Let be the periodization of with respect to period . That is is a circular filter, which means that the component indexes are residue classes with respect to the modulus . Then with the upsampling operator it holds
Actually not convolutions are necessary, but only ones, when applying the strategy of efficient computation of powers. Even more the approach can be further sped up using the Fast Fourier transform.
  • From the previous statement we can derive an estimate of the spectral radius of . It holds
where is the size of the filter and if all eigenvalues are real, it is also true that
,
where .

References

  • Gilbert Strang: Eigenvalues of and convergence of the cascade algorithm. IEEE Transactions on Signal Processing, 44:233-238, 1996.
  • Henning Thielemann: Optimally matched wavelets, 2006 (contains proofs of the above properties)
In mathematics, in the area of wavelet analysis, a refinable function is a function which fulfills some kind of self-similarity. A function is called refinable with respect to the mask if
This condition is called refinement equation,
..... Click the link for more information.
A wavelet is a kind of mathematical function used to divide a given function into different frequency components and study each component with a resolution that matches its scale. A wavelet transform is the representation of a function by wavelets.
..... Click the link for more information.
Finite element analysis (FEA) is a computer simulation technique used in engineering analysis. It uses a numerical technique called the finite element method (FEM). There are many finite element software packages, both free and proprietary.
..... Click the link for more information.
In signal processing, downsampling (or "subsampling") is the process of reducing the sampling rate of a signal. This is usually done to reduce the data rate or the size of the data.
..... Click the link for more information.
In mathematics, a Sylvester matrix is a matrix associated to two polynomials that gives us some information about those polynomials. It is named for James Joseph Sylvester.
..... Click the link for more information.
resultant of two monic polynomials and over a field is defined as the product



of the differences of their roots, where and take on values in the algebraic closure of .
..... Click the link for more information.
Euclidean algorithm (also called Euclid's algorithm) is an algorithm to determine the greatest common divisor (GCD) of two elements of any Euclidean domain (for example, the integers).
..... Click the link for more information.
In linear algebra, the trace of an n-by-n square matrix A is defined to be the sum of the elements on the main diagonal (the diagonal from the upper left to the lower right) of A, i.e.
..... Click the link for more information.
convolution is a mathematical operator which takes two functions f and g and produces a third function that in a sense represents the amount of overlap between f and a reversed and translated version of g.
..... Click the link for more information.
In algebra, a determinant is a function depending on n that associates a scalar, det(A), to every n×n square matrix A. The fundamental geometric meaning of a determinant is as the scale factor for volume when A
..... Click the link for more information.
invertible or non-singular if there exists an n-by-n matrix such that



where denotes the n-by-n identity matrix and the multiplication used is ordinary matrix multiplication.
..... Click the link for more information.
Null space may refer to:
  • kernel (matrix) In linear algebra, the null space of a matrix.
  • kernel (linear operator) In advanced linear algebra and functional analysis, the null space of a linear operator.

..... Click the link for more information.
In mathematics, the spectral radius of a matrix or a bounded linear operator is the supremum among the moduli of the elements in its spectrum, which is sometimes denoted by ρ(·).
..... Click the link for more information.
Upsampling is the process of increasing the sampling rate of a signal.

The upsampling factor (commonly denoted by L) is usually an integer or a rational fraction greater than unity.
..... Click the link for more information.
fast Fourier transform (FFT) is an efficient algorithm to compute the discrete Fourier transform (DFT) and its inverse. FFTs are of great importance to a wide variety of applications, from digital signal processing and solving partial differential equations to algorithms for
..... Click the link for more information.
In mathematics, the spectral radius of a matrix or a bounded linear operator is the supremum among the moduli of the elements in its spectrum, which is sometimes denoted by ρ(·).
..... Click the link for more information.


This article is copied from an article on Wikipedia.org - the free encyclopedia created and edited by online user community. The text was not checked or edited by anyone on our staff. Although the vast majority of the wikipedia encyclopedia articles provide accurate and timely information please do not assume the accuracy of any particular article. This article is distributed under the terms of GNU Free Documentation License.
Herod_Archelaus


page counter