A fast matrix iterative technique for the WLS design of 2-D quadrantally symmetic FIR filters

High computational complexity is a major problem encountered in the optimal design of two-dimensional (2-D) finite impulse response (FIR) filters. In this paper, we present an iterative matrix solution with very low complexity to the weighted least square (WLS) design of 2-D quadrantally symmetric FIR filters with two-valued weighting functions. Firstly, a necessary and sufficient condition for the WLS design of 2-D quadrantally symmetric filters with general nonnegative weighting functions is obtained. Then, based on this optimality condition, a novel iterative algorithm is derived for the WLS design problem with a two-valued weighting function. Because the filter parameters are arranged in their natural 2-D form and the transition band is not sampled, the computation amount of the proposed algorithm is reduced significantly, especially for high-order filters. The exponential convergence of the algorithm is established, and its computational complexity is estimated. Design examples demonstrating the convergence rate and solution accuracy of the algorithm, as well as the relation between the iteration number of the algorithm and the size and transition-band width of the filter are given.