home

PDF (letter size)

Mapping the system function from the s-plane to the z-plane in the presence of multiple order poles.

Nasser M. Abbasi

April 22, 2010   Compiled on January 29, 2024 at 2:53am

Given H(s) of order N with all its poles pi being distinct, it can be expressed in terms of partial fraction expansion in the form of H(s)=k=1NAkspk and the resulting H(z) can be found to be k=1NzAkzepkT where T is the sampling period.

In the case when H(s) contains a pole q of order 2, then H(s) can be written as (k=1N2Akspk)+Aq(sq)2 and the resulting H(z) can be found to be (k=1N2zAkzepkT)+TzeqT(eqTz)2.

In the case when H(s) contains a pole q of order 3, then H(s) can be written as(k=1N3Akspk)+Aq(sq)3 and the resulting H(z) can be found to be (k=1N3zAkzepkT)+(e2qTT2z+eqTT2z22(eqTz)3).

The following table was generated in order to obtain the general formula. This table below shows only the part of H(z) due to the multiple order pole.

n pole order H(z)
2 TzeqT(eqTz)2
3 e2qTT2z+eqTT2z22(eqTz)3
4 e3qTT3z+4e2qTT3z2+eqTT3z36(eqTz)4
5 e4qTT4z11e3qTT4z211e2qTT4z3eqTT4z424(eqTz)5
6 e5qTT5z+26e4qTT5z2+66e3qTT5z3+26e2qTT5z4+eqTT5z5120(eqTz)6

It is easy to see that the denominator of H(z) has the general form (n1)!(eqTz)n where n is the pole order, the hard part is to find the general formula for the numerator. The following table is a rewrite of the above table, where only the numerator is show, and eqT was written as A to make it easier to see the general pattern

n pole order numerator of H(z)
2 (1)n(AT)z
3 (1)n[(AT)2zA(Tz)2]
4 (1)n[(AT)3z+4A2T3z2+A(Tz)3]
5 (1)n[(AT)4z11A3T4z211A2T4z3A(Tz)4]
6 (1)n[(AT)5z+26A4T5z2+66A3T5z3+26A2T5z4+A(Tz)5]

I am trying to determine the general formula to generate the above. This seems to involve some combination of binomial coefficient. But so far, I did not find the general formula.

1 References

  1. Digital signal processing, by Oppenheim and Scafer, page 201
  2. Mathematica software version 7