<< Chapter < Page Chapter >> Page >
H R ( z ) = z - M m + 1 l = 0 M - 1 z l H l ( z - M ) = l = 0 M - 1 z - ( M - 1 - l ) ( z - m M + M H l ( z - M ) ) = l = 0 M - 1 z - l H M - 1 - l R ( z M ) .

Therefore

( H R ) l ( z ) = ( H M - 1 - l ) R ( z )

and for linear-phase H ( z ) , since H R ( z ) = ± H ( z ) ,

H l ( z ) = ± ( H R ) M - 1 - l ( z ) .

Lemma 2 For even M , H p ( z ) is of the form

  • W 0 ( z ) W 1 ( z ) J W 0 ( z ) V ( - 1 ) M / 2 J W 1 ( z ) V = I 0 0 J W 0 ( z ) W 1 ( z ) W 0 ( z ) V ( - 1 ) M / 2 W 1 ( z ) V
  • W 0 ( z ) W 1 ( z ) J J W 1 R ( z ) V ( - 1 ) M / 2 J W 0 R ( z ) J V = I 0 0 J W 0 ( z ) W 1 ( z ) J W 1 R ( z ) V ( - 1 ) M / 2 W 0 R ( z ) J V
  • W 0 ( z ) D 0 W 0 R ( z ) J W 1 ( z ) D 1 W 1 R ( z ) J = W 0 ( z ) D 0 W 0 R ( z ) W 1 ( z ) D 1 W 1 R ( z ) I 0 0 J
    or
    Q W 0 ( z ) W 0 R ( z ) J W 1 ( z ) - W 1 R ( z ) J = Q W 0 ( z ) W 0 R ( z ) W 1 ( z ) - W 1 R ( z ) I 0 0 J
  • W 0 ( z ) D W 0 R ( z ) J J D W 0 ( z ) V ( - 1 ) M / 2 J W 0 R ( z ) J V = I 0 0 J W 0 ( z ) D W 0 R ( z ) J D W 0 ( z ) V ( - 1 ) M / 2 W 0 R ( z ) J V
  • W 0 ( z ) D W 0 R ( z ) J J W 0 ( z ) V ( - 1 ) M / 2 J D W 0 R ( z ) J V = I 0 0 J D W 0 ( z ) D W 0 R ( z ) J D W 0 ( z ) V ( - 1 ) M / 2 W 0 R ( z ) J V

Thus in order to generate H p ( z ) for all symmetries other than P S , we need a mechanism that generates a pair of matrices and their reflection (i.e., W 0 ( z ) , W 1 ( z ) W 0 R ( z ) and W 1 R ( z ) ). In the scalar case, there are two well-known lattice structures,that generate such pairs. The first case is the orthogonal lattice [link] , while the second is the linear-prediction lattice [link] . A K t h order orthogonal lattice is generated by the product

i = 0 K a i z - 1 b i - b i z - 1 a i a 0 b 0 - b 0 a 0 = def Y 0 ( z ) Y 1 ( z ) - Y 1 R ( z ) Y 0 R ( z ) = def X ( z ) .

This lattice is always invertible (unless a i and b i are both zero!), and the inverse is anticausal since

a i z - 1 b i - b i z - 1 a i - 1 = 1 a 1 2 + b i 2 a i - b i z b i z a i .

As we have seen, this lattice plays a fundamental role in the theory of two-channel FIR unitary modulated filter banks.The hyperbolic lattice of order K generates the product

i = 0 K a i z - 1 b i b i z - 1 a i a 0 b 0 b 0 a 0 = def Y 0 ( z ) Y 1 ( z ) Y 1 R ( z ) Y 0 R ( z ) = def X ( z ) .

where Y 0 ( z ) and Y 1 ( z ) are of order K . This lattice is invertible only when a i 2 b i 2 (or equivalently ( a i + b i ) / 2 and ( a i - b i ) / 2 are nonzero) in which case the inverse is noncausal since

a i z - 1 b i b i z - 1 a i - 1 = 1 a 1 2 - b i 2 a i - b i - z b i z a i .

Since the matrix a i b i b i a i can be orthogonal iff a i , b i = ± 1 , 0 , or a i , b i = 0 , ± 1 , the ( 2 × 2 ) matrix generated by the lattice can never be unitary.

Formally, it is clear that if we replace the scalars a i and b i with square matrices of size M / 2 × M / 2 then we would be able to generate matrix versions of these two lattices which can then be used togenerate filter banks with the symmetries we have considered. We will shortly see that both the lattices can generate unitary matrices, and this willlead to a parameterization of FIR unitary H p ( z ) for PCS, linear-phase, and PCS plus linear-phase symmetries. We prefer to call the generalization ofthe orthogonal lattice, the antisymmetric lattice and to call the generalization of the hyperbolic lattice, the symmetric lattice, which should be obvious from the form of the product. The reason for this is that the antisymmetric lattice may not generate a unitary matrix transfer function (in the scalar case, the 2 × 2 transfer function generated is always unitary). The antisymmetric lattice is defined by the product

X ( z ) = def i = 1 K A i z - 1 B i - B i z - 1 A i A 0 B 0 - B 0 A 0

where A i and B i are constant square matrices of size M / 2 × M / 2 . It is readily verified that X ( z ) is of the form

X ( z ) = Y 0 ( z ) Y 1 ( z ) - Y 1 R ( z ) Y 0 R ( z )

Given X ( z ) , its invertibility is equivalent to the invertibility of the constant matrices,

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Wavelets and wavelet transforms. OpenStax CNX. Aug 06, 2015 Download for free at https://legacy.cnx.org/content/col11454/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Wavelets and wavelet transforms' conversation and receive update notifications?

Ask