Matrix Extension with Symmetry and Its Application to Symmetric Orthonormal Multiwavelets

Let P be an r×s matrix of Laurent polynomials with symmetry such that P(z)P∗(z) = Ir for all z ∈ C\{0} and the symmetry of P is compatible. The matrix extension problem with symmetry is to find an s × s square matrix Pe of Laurent polynomials with symmetry such that [Ir,0]Pe = P (that is, the submatrix of the first r rows of Pe is the given matrix P), Pe is paraunitary satisfying Pe(z)P∗e(z) = Is for all z ∈ C\{0}, and the symmetry of Pe is compatible. Moreover, it is highly desirable in many applications that the support of the coefficient sequence of Pe can be controlled by that of P. In this paper, we completely solve the matrix extension problem with symmetry on deriving such a desired matrix Pe from a given matrix P. Furthermore, using a cascade structure, we obtain a complete representation of any r × s paraunitary matrix P having compatible symmetry, which in turn leads to an algorithm for deriving a desired matrix Pe from a given matrix P. Matrix extension plays an important role in many areas such as wavelet analysis, electronic engineering, system sciences, and so on. As an application of our general results on matrix extension with symmetry, we obtain a satisfactory algorithm for constructing symmetric orthonormal multiwavelets by deriving high-pass filters with symmetry from any given low-pass filters with symmetry. Several examples are provided to illustrate the proposed algorithms and results in this paper.