Polynomial division to separate FIR and IIR in cascade

Question

I have a FIR filter $L(z)$ and an IIR filter $H(z) = \frac{B(z)}{A(z)}$ which are cascaded together as $L(z)H(z) = L(z)\frac{B(z)}{A(z)}$.

For specific reasons I want to calculate a parallel topology that is equivalent to the original $L(z)H(z)$. I can achieve this using long division to divide $L(z)B(z)$ by $A(z)$ to obtain:

$L(z)H(z) = \frac{L(z)B(z)}{A(z)} = Q(z) + \frac{R(z)}{A(z)}$

Where $Q(z)$ is the quotient and $R(z)$ is the remainder of the long division of $L(z)B(z)$ by $A(z)$.

However, the resulting IIR, $\frac{R(z)}{A(z)}$, has a delay equal to the number of taps of the original FIR $L(z)$

Is it possible to do some clever re-arranging before the division so that the resulting IIR part of the long division is not delayed?

Of course one approach is to flip the filters so the magnitude of the poles are greater than one and perform long division using an unstable filter and flip back (which is a simple change of variable $z \to \frac{1}{z}$ and back again), but this isn't a useful solution.

Answer 1

The way that you currently split the FIR from the IIR part results in an FIR filter that exactly matches the first $N$ samples of the impulse response of the original cascade. Consequently, the parallel IIR filter must have $N$ initial zeros in its impulse response, otherwise it would mess with that first part of the original response which is already taken care of by the FIR filter. If you choose another FIR filter, then the corresponding IIR filter won't have as many initial zeros (i.e., delays without multipliers in its implementation). The question is why you would want this, because in that case the number of additions and multiplications necessary for implementing the IIR filter would increase. You would get the same number of delays, but instead of just delaying the input you would get non-zero filter coefficients.

Note that if $Q(z)$ denotes the transfer function of the FIR filter, and $R(z)$ is the numerator of the IIR filter, you just need to satisfy the following equation:

$$R(z)=L(z)B(z)-A(z)Q(z)\tag{1}$$

There are infinitely many solutions. You can choose any arbitrary FIR filter $Q(z)$, and the numerator $R(z)$ of the corresponding IIR filter follows from $(1)$. The most efficient solution, however, is obtained by choosing the FIR filter that is just a truncated version of the impulse response of the original cascade. In that case, the IIR filter $R(z)/A(z)$ delays the input by the length of the FIR filter, and takes care of the remaining tail of the original impulse response. That solution is just the one obtained by long division.

Also take a look at this related question and its answer.