Abstract:Systems made up of arrays of sensors (seismographs, hydrophones, antennas, meteorological stations) receive geophysical data on a number of channels. Such multichannel data can be processed by means of a multichannel least‐squares (Wiener) digital filter; the filter is determined by requiring that the sum of the mean‐square errors between the desired outputs and the actual outputs for each channel be minimized. This requirement leads to a set of simultaneous linear equations, called the normal equations, where each element is itself a square matrix. The coefficients of the digital filter are given by the solution of these normal equations. An exact recursive procedure is derived for the numerical solution of the normal equations. In this recursive procedure advantage is taken of the fact that the matrix of the normal equations is an autocorrelation matrix, and hence is in the form of a Toeplitz matrix, namely one with equal elements along any diagonal. A set of auxiliary coefficients that are generated in the recursive procedure are the coefficients of matrix‐valued polynomials which are the counterparts of the classical (scalar‐valued) polynomials orthogonal on the unit circle. This set of auxiliary coefficients may also be used for a sideward recursion that corresponds to shifting the time index of the right‐hand side of the normal equations. This sideward recursion is valuable in applications because it represents a relatively inexpensive way of determining the filter that incorporates the optimum time delay between the multichannel input signal and the multichannel output signal. The machine time required to solve the normal equations for a filter with m (matrix‐valued) coefficients is proportional to m2 for the proposed recursive procedure, as compared with m3 for conventional techniques for the solution of simultaneous equations. Another advantage of using the recursive procedure is that it only requires computer storage space proportional to m, rather than to m2 as in the case of the conventional techniques.

Exact Least Squares Algorithm for Signal Matched Multirate Whitening Filter Bank: Part I

Exact Least Squares Algorithm for Signal Matched Synthesis Filter Bank: Part II

Iterative algorithm for matched-filter design in minimum variance spectral estimation

Special low-order IIR filter bank design

Bi-iterative least squares algorithms for blind channel identification and equalization with second-order statistics

The Analysis and Design of Two-Dimensional Nearly-Orthogonal Symmetric Wavelet Filter Banks

A unified framework for synthesis of cosine-modulated filter banks and corresponding wavelets

Adaptive Reduced-Rank Linear Equalization of Multiple Input Multiple Output Frequency Selective Channels Based on Multistage Wiener Filter

Properties for Perfect Reconstruction Using Multirate Linear Phase Two-Dimensional Nonseparable Filter Banks and Its Implementation

Design perfect reconstruction cosine-modulated filter bank by variable neighbourhood search-least-mean-square error

Design of two-channel filter bank using nature inspired optimization based fractional derivative constraints

The Analysis and Design of One-Dimensional Nearly-Orthogonal Symmetric Wavelet Filter Banks

On Lattice Factorization of Symmetric-Antisymmetric Multifilter Banks.

On Minimal Lattice Factorizations of Symmetric-Antisymmetric Multifilterbanks

Filter Bank Design Based on Minimization of Individual Aliasing Terms for Minimum Mutual Information Subband Adaptive Beamforming.

Synthesis of quasi-optimal fast filters by the least square criterion

A Reduced Gaussian Kernel Least-Mean-Square Algorithm for Nonlinear Adaptive Signal Processing

Recursive solution to the multichannel filtering problem

A simple design method of reduced-order filters and its applications to multirate filter bank design

Constrained Digital Matched Filter Method for Optimum Filter Synthesis

Reduced-Rank Adaptive Filtering Based on Joint Iterative Optimization of Adaptive Filters