A recursive method for the approximation of LTI systems using subband processing