A computational method for the LU decomposition of rectangular matrices and its implications to the realization of 2-D digital filters

The realization of 2-D digital filters based on the lower-upper triangular decomposition of the coefficient matrix is investigated. A numerical method based on the QA decomposition, which has some important characteristics, is proposed for reaching the LU structure. The coefficients in the final LU structure have values favorable to fixed-point arithmetic implementation. Furthermore, the QR structure can be used for the realization and possesses good numerical characteristics in terms of the approximate decomposition scheme. The symmetry in the impulse response coefficient matrix of an octagonally symmetric 2-D FIR filter is utilized to reduce the computational effort spent in the decomposition and the total number of multipliers in the final realization structure