Use the convolution integral closed in the left half-plane. The z-transform is the sum of residues of X(s)z−esTz at the poles of X(s). This is valid because X(s) has 4 poles and no zeros.
X(z)=poles of X(s)∑Res[(s+1)2(s+2)(s+3)1⋅z−esTz]
Poles: a double pole at s=−1, and simple poles at s=−2 and s=−3.
Residue at s=−2
R2=[(s+1)2(s+3)1⋅z−esTz]s=−2=(1)(1)1⋅z−e−2Tz=z−e−2Tz
Residue at s=−3
R3=[(s+1)2(s+2)1⋅z−esTz]s=−3=(4)(−1)1⋅z−e−3Tz=−41⋅z−e−3Tz
Residue at the double pole s=−1
Let F(s)=(s+2)(s+3)1, so F(−1)=21 and
F′(s)=−(s+2)2(s+3)22s+5⇒F′(−1)=−43
R1=dsd[F(s)z−esTz]s=−1=F′(−1)z−e−Tz+F(−1)(z−e−T)2Tze−T=−43⋅z−e−Tz+2Te−T⋅(z−e−T)2z
General result
X(z)=2Te−T⋅(z−e−T)2z−43⋅z−e−Tz+z−e−2Tz−41⋅z−e−3Tz
Substitute T=1
e−1=0.36788, e−2=0.13534, e−3=0.04979, and 2Te−T=0.18394:
X(z)=(z−0.36788)20.18394z−z−0.367880.75z+z−0.13534z−z−0.049790.25z
Combined into a single ratio:
X(z)=z4−0.92088z3+0.27828z2−0.030011z+0.0009120.030919z3+0.023057z2+0.000939z
The denominator is (z−0.36788)2(z−0.13534)(z−0.04979).
Check
The inverse Laplace transform is x(t)=41[(2t−3)e−t+4e−2t−e−3t]. Its samples are:
| k | 0 | 1 | 2 | 3 |
|---|
| x(k) | 0 | 0.03092 | 0.05153 | 0.03979 |
Long division of X(z) gives 0.03092z−1+0.05153z−2+0.03979z−3+…, which matches.
Answer: X(z)=(z−0.3679)20.18394z−z−0.36790.75z+z−0.1353z−z−0.04980.25z.