Skip to main content

Convolution and Z-transform via convolution integral

Short answer

Convolution combines an input with a system's impulse response to yield the output; in the s or z domain it becomes multiplication. For X(s)=1/s^2(s+1), the convolution integral method gives X(z)=Tz/(z-1)^2 - z/(z-1) + z/(z-e^{-T}).

What do you mean by convolution? Obtain Z-transform of X(s)=1/s²(s+1) by use of the convolution integral in the left half s-plane.

Convolution

Convolution combines two signals to give the response of a linear time-invariant system. If gg is the impulse response and xx is the input, the output is the convolution y=g∗xy = g * x.

  • Continuous time: y(t)=∫0tg(t−τ) x(τ) dτy(t) = \displaystyle\int_0^t g(t-\tau)\,x(\tau)\,d\tau
  • Discrete time (convolution sum): y(k)=∑h=0kg(k−h) x(h)y(k) = \displaystyle\sum_{h=0}^{k} g(k-h)\,x(h)

The key property is that convolution in time becomes multiplication in the transform domain: Y(s)=G(s)X(s)Y(s) = G(s)X(s) and Y(z)=G(z)X(z)Y(z) = G(z)X(z). This is why pulse transfer functions are useful.

Convolution integral for the z-transform

The impulse-sampled signal is x∗(t)=x(t)∑kδ(t−kT)x^*(t) = x(t)\sum_k\delta(t-kT), which is a product in time. Its Laplace transform is therefore a complex convolution of X(s)X(s) with the transform of the impulse train, 11−e−Ts\dfrac{1}{1 - e^{-Ts}}:

X∗(s)=12πj∫c−j∞c+j∞X(p) 11−e−T(s−p) dpX^*(s) = \frac{1}{2\pi j}\int_{c-j\infty}^{c+j\infty}X(p)\,\frac{1}{1 - e^{-T(s-p)}}\,dp

The contour is closed in the left half of the p-plane, which encloses the poles of X(p)X(p). This is valid when X(s)X(s) has at least two more poles than zeros. Setting eTs=ze^{Ts} = z gives:

X(z)=∑poles of X(s)Res[X(s)zz−esT]X(z) = \sum_{\text{poles of } X(s)}\text{Res}\left[X(s)\frac{z}{z-e^{sT}}\right]

Apply to X(s)=1s2(s+1)X(s) = \dfrac{1}{s^2(s+1)}

Poles: a double pole at s=0s = 0 and a simple pole at s=−1s = -1.

Simple pole s=−1s = -1:

R1=[1s2⋅zz−esT]s=−1=zz−e−TR_1 = \left[\frac{1}{s^2}\cdot\frac{z}{z - e^{sT}}\right]_{s=-1} = \frac{z}{z - e^{-T}}

Double pole s=0s = 0:

R2=dds[1s+1⋅zz−esT]s=0=[−1(s+1)2⋅zz−esT+1s+1⋅zTesT(z−esT)2]s=0=−zz−1+Tz(z−1)2\begin{aligned} R_2 &= \frac{d}{ds}\left[\frac{1}{s+1}\cdot\frac{z}{z - e^{sT}}\right]_{s=0} \\ &= \left[-\frac{1}{(s+1)^2}\cdot\frac{z}{z-e^{sT}} + \frac{1}{s+1}\cdot\frac{zTe^{sT}}{(z-e^{sT})^2}\right]_{s=0} \\ &= -\frac{z}{z-1} + \frac{Tz}{(z-1)^2} \end{aligned}

Sum of residues

X(z)=Tz(z−1)2−zz−1+zz−e−TX(z) = \frac{Tz}{(z-1)^2} - \frac{z}{z-1} + \frac{z}{z - e^{-T}}

Over the common denominator:

X(z)=z[(T−1+e−T)z+(1−e−T−Te−T)](z−1)2(z−e−T)X(z) = \frac{z\left[(T - 1 + e^{-T})z + (1 - e^{-T} - Te^{-T})\right]}{(z-1)^2(z - e^{-T})}

Check

x(t)=t−1+e−tx(t) = t - 1 + e^{-t} (from partial fractions 1s2−1s+1s+1\frac{1}{s^2} - \frac{1}{s} + \frac{1}{s+1}). Its z-transform is Tz(z−1)2−zz−1+zz−e−T\frac{Tz}{(z-1)^2} - \frac{z}{z-1} + \frac{z}{z-e^{-T}}, which is the same.

For T=1T = 1 s: X(z)=z(0.3679z+0.2642)(z−1)2(z−0.3679)X(z) = \dfrac{z(0.3679z + 0.2642)}{(z-1)^2(z-0.3679)}.

Answer: X(z)=z[(T−1+e−T)z+(1−e−T−Te−T)](z−1)2(z−e−T)X(z) = \dfrac{z[(T-1+e^{-T})z + (1-e^{-T}-Te^{-T})]}{(z-1)^2(z-e^{-T})}.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Digital Control System chapter 2 questions →