Inverse Z-transform of a rational function

Question

Let $a \in (-1,1)$ and $n=0,1,2,\cdots$. By decomposing a rational function into simple fractions we have derived the following identity: \begin{equation} \left.\frac{1}{n!} \frac{d^n}{d z^n} \left( \frac{A z + B}{z^2-(1+a) z + 1}\right) \right|_{z=0} = \frac{A \sin(\phi \cdot n) + B \sin(\phi \cdot (n+1))}{\sin(\phi)} \end{equation} where $\phi := \arccos((1+a)/2)$. In other words we have shown that the inverse Z-transform of the rational function in parentheses on the left hand side is a waveform with both time independent frequency and amplitude. Now the question would be what happens if the second order polynomial in the denominator is replaced by some integer power of it. What kind of waveform are we going to get in this case? In other words we ask for the following: \begin{equation} \left.\frac{1}{n!} \frac{d^n}{d z^n} \left( \frac{A z + B}{(z^2-(1+a) z + 1)^p}\right) \right|_{z=0} = ? \end{equation} where $p=1,2,3,\cdots$.

Answer 1

We start from the following identity: \begin{equation} \frac{1}{\left[z^2-(1+a) z+1\right]^p} = 2 Re\left[ \sum\limits_{l=0}^{p-1} \binom{2p-2-l}{p-1} \frac{(-1)^p}{(z-\exp(\imath \phi))^{l+1} (-2 \imath \sin(\phi))^{2 p-l-1}} \right] \end{equation} which is just a simple fraction decomposition of the left hand side. Now since $A z + B = A\cdot (z - \exp(\imath \phi)) + A \exp(\imath \phi) + B$ we can express the desired rational function through inverse powers of monomials of $z$ only, i.e. through linear combinations of terms of the kind $(z-\exp(\imath \phi))^{-l}$ and their complex conjugates only. Those monomials are very simple to differentiate and therefore we arrive at the result without any substantial problems. Having said all this we just state the final result: \begin{eqnarray} &&\left.\frac{1}{n!} \frac{d^n}{d z^n} \left(\frac{A z+ B}{[z^2-(1+a) z+1]^p}\right)\right|_{z=0}=\\ && \sum\limits_{l=0}^{p-1} \frac{2(-1)^{l+p+1}}{ (2\sin (\phi ))^{-l+2 p-1}}\binom{l+n}{n} \binom{2 p-l-2}{p-1} \frac{1}{l+n}\\ && \left(A n \cos \left(\phi (l+n)-\frac{1}{2} \pi (-l+2 p-1)\right)+B (l+n) \cos \left(\phi (l+n+1)-\frac{1}{2} \pi (-l+2 p-1)\right)\right) \end{eqnarray}