Math/Control Systems Tutorial – Partial Fraction Expansion with MATLAB – Part III

This is the final lecture on partial fraction expansion. The first and second parts of this lecture can be found here and here. The video accompanying this post can be found here:

This lecture covers the case when the rational (transfer) function has repeated poles (poles are roots of the expression in the denominator). Let the poles of the transfer function be . We assume that the pole is repeated three times. Then, the partial fraction expansion of becomes:

(1)  

Our goal is to find the constants .
If we multiply the last expression with (s-p_{1})^{3}, we obtain

(2)  

From the last equation, we immediately obtain

(3)  

Taking the first derivative of (2), we obtain:

(4)  

If we evaluate the last expression for , the first term on the right-hand-side is zero. Also, the derivative of the expression in the brackets is zero becasue of the third power of that multiplies all the terms. Consequently, we obtain

(5)  

By taking the derivative of (4), we obtain

(6)  

Again, the second term on the right-hand-side of the last expression will be zero due to the term multiplying all the terms. From the last equation, we obtain

(7)  

This procedure can be generalized when a transfer function has repeated poles. The coefficients then become

(8)  

The remaining coefficients corresponding to distinct poles are computed using previously explained approaches that can be found here and here.
Let us now illustrate these derivations with an example.

Example 1: Find the inverse Laplace transform of the function

(9)  

The partial fraction expansion is

(10)  

The constant is

(11)  

The constant is

(12)  

The constant is

(13)  

The constant is

(14)  

Finally, we have

(15)  

The inverse Laplace transform is given by

(16)  

where is the Heaviside step function.

The MATLAB code for solving this problem is given below

The output is