Differentiation by integration using orthogonal polynomials, a survey

This survey paper discusses the history of approximation formulas for n-th order derivatives by integrals involving orthogonal polynomials. There is a large but rather disconnected corpus of literature on such formulas. We give some results in greater generality than in the literature. Notably we unify the continuous and discrete case. We make many side remarks, for instance on wavelets, Mantica's Fourier-Bessel functions and Greville's minimum R_alpha formulas in connection with discrete smoothing.Comment: 35 pages, 3 figures; minor corrections, subsection 3.11 added; accepted by J. Approx. Theor

All-Pole Recursive Digital Filters Design Based on Ultraspherical Polynomials

A simple method for approximation of all-pole recursive digital filters, directly in digital domain, is described. Transfer function of these filters, referred to as Ultraspherical filters, is controlled by order of the Ultraspherical polynomial, nu. Parameter nu, restricted to be a nonnegative real number (nu ≥ 0), controls ripple peaks in the passband of the magnitude response and enables a trade-off between the passband loss and the group delay response of the resulting filter. Chebyshev filters of the first and of the second kind, and also Legendre and Butterworth filters are shown to be special cases of these allpole recursive digital filters. Closed form equations for the computation of the filter coefficients are provided. The design technique is illustrated with examples

Computation of the para-pseudoinverse for oversampled filter banks: Forward and backward Greville formulas

This is the author's accepted manuscript. The final published article is available from the link below. Copyright @ 2008 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other users, including reprinting/ republishing this material for advertising or promotional purposes, creating new collective works for resale or redistribution to servers or lists, or reuse of any copyrighted components of this work in other works.Frames and oversampled filter banks have been extensively studied over the past few years due to their increased design freedom and improved error resilience. In frame expansions, the least square signal reconstruction operator is called the dual frame, which can be obtained by choosing the synthesis filter bank as the para-pseudoinverse of the analysis bank. In this paper, we study the computation of the dual frame by exploiting the Greville formula, which was originally derived in 1960 to compute the pseudoinverse of a matrix when a new row is appended. Here, we first develop the backward Greville formula to handle the case of row deletion. Based on the forward Greville formula, we then study the computation of para-pseudoinverse for extended filter banks and Laplacian pyramids. Through the backward Greville formula, we investigate the frame-based error resilient transmission over erasure channels. The necessary and sufficient condition for an oversampled filter bank to be robust to one erasure channel is derived. A postfiltering structure is also presented to implement the para-pseudoinverse when the transform coefficients in one subband are completely lost

Oversampling in shift-invariant spaces with a rational sampling period

8 pages, no figures.It is well known that, under appropriate hypotheses, a sampling formula allows us to recover any function in a principal shift-invariant space from its samples taken with sampling period one. Whenever the generator of the shift-invariant space satisfies the Strang-Fix conditions of order r, this formula also provides an approximation scheme of order r valid for smooth functions. In this paper we obtain sampling formulas sharing the same features by using a rational sampling period less than one. With the use of this oversampling technique, there is not one but an infinite number of sampling formulas. Whenever the generator has compact support, among these formulas it is possible to find one whose associated reconstruction functions have also compact support.This work has been supported by the Grant MTM2009-08345 from the D.G.I. of the Spanish Ministerio de Ciencia y Tecnología

Construction of M - Band bandlimited wavelets for orthogonal decomposition

While bandlimited wavelets and associated IIR filters have shown serious potential in areas of pattern recognition and communications, the dyadic Meyer wavelet is the only known approach to construct bandlimited orthogonal decomposition. The sine scaling function and wavelet are a special case of the Meyer. Previous works have proposed a M - Band extension of the Meyer wavelet without solving the problem. One key contribution of this thesis is the derivation of the correct bandlimits for the scaling function and wavelets to guarantee an orthogonal basis. In addition, the actual construction of the wavelets based upon these bandlimits is developed. A composite wavelet will be derived based on the M scale relationships from which we will extract the wavelet functions. A proper solution to this task is proposed which will generate associated filters with the knowledge of the scaling function and the constraints for Mband orthogonality

Cyclic LTI systems in digital signal processing

Cyclic signal processing refers to situations where all the time indices are interpreted modulo some integer L. In such cases, the frequency domain is defined as a uniform discrete grid (as in L-point DFT). This offers more freedom in theoretical as well as design aspects. While circular convolution has been the centerpiece of many algorithms in signal processing for decades, such freedom, especially from the viewpoint of linear system theory, has not been studied in the past. In this paper, we introduce the fundamentals of cyclic multirate systems and filter banks, presenting several important differences between the cyclic and noncyclic cases. Cyclic systems with allpass and paraunitary properties are studied. The paraunitary interpolation problem is introduced, and it is shown that the interpolation does not always succeed. State-space descriptions of cyclic LTI systems are introduced, and the notions of reachability and observability of state equations are revisited. It is shown that unlike in traditional linear systems, these two notions are not related to the system minimality in a simple way. Throughout the paper, a number of open problems are pointed out from the perspective of the signal processor as well as the system theorist

Orthonormal and biorthonormal filter banks as convolvers, and convolutional coding gain

Convolution theorems for filter bank transformers are introduced. Both uniform and nonuniform decimation ratios are considered, and orthonormal as well as biorthonormal cases are addressed. All the theorems are such that the original convolution reduces to a sum of shorter, decoupled convolutions in the subbands. That is, there is no need to have cross convolution between subbands. For the orthonormal case, expressions for optimal bit allocation and the optimized coding gain are derived. The contribution to coding gain comes partly from the nonuniformity of the signal spectrum and partly from nonuniformity of the filter spectrum. With one of the convolved sequences taken to be the unit pulse function,,e coding gain expressions reduce to those for traditional subband and transform coding. The filter-bank convolver has about the same computational complexity as a traditional convolver, if the analysis bank has small complexity compared to the convolution itself