In this control engineering and control, we explain how to compute an inverse Laplace transform of transfer functions and expressions in the s-domain. We consider the case when the poles of complex expressions and transfer functions have distinct real poles. To compute the inverse Laplace transform, we use the partial fraction expansion and cover-up methods. The YouTube tutorial accompanying this webpage is given below.
Compute the Partial Fraction Expansion by Using the High-School-Level Mathematics
Let us consider the following transfer function
(1)
and let us assume that our goal is to compute the inverse Laplace transform of this function. If we can somehow decompose this transfer function into the following form:
(2)
where and are constants that need to be determined, then our problem would be solved since we can easily compute the inverse Laplace transforms of the terms and .
The equation (2) is referred to as the partial fraction expansion. To compute this decomposition, let us first use a high-school-level approach. This approach originates from computing integrals of rational functions. From the last equation, we have:
(3)
From the last equation, we obtain the following equation
(4)
Our goal is to determine the constants and , such that (4) and (1) are identical. By equating the numerators of these two fractions, we obtain
(5)
The expressions on both sides of this equation are actually polynomials in . Two polynomials are equal if their coefficients multiplying the corresponding powers of are equal. Consequently, we have
(6)
We obtain and . In order to verify this solution, we can simply substitute these values of and in (2) and investigate if we can get the original expression, given by (1). We can also use MATLAB to verify our computations:
syms s
F=1/((s-2)*(s-3))
partfrac(F)
The previous code produces the following result
F =
1/((s - 2)*(s - 3))
Factorized =
1/(s - 3) - 1/(s - 2)
Obviously, we have correctly computed the solution. Next, let us explain the MATLAB code. With “syms s”, we define a symbolic variable . On code line 3 we can define a symbolic expression for our function. Code line 5 is used to compute the expansion.
The above-presented approach can also be applied to more complex functions with a large number of poles. However, it becomes tedious to solve a system of linear equations to compute every coefficient of the partial fraction expansion. Consequently, we need to search for an alternative method to compute the partial fraction expansion.
Cover-Up Method for Computing the Inverse Laplace Transform
Here, we present a more elegant and easily applicable approach for computing the partial fraction expansion and for computing the inverse Laplace transform. The approach presented below is applicable to expressions with (real) distinct poles.
First, let us explain what are the poles and zeros of a function. Consider a general form of a complex function
(7)
The last expression can be written as follows:
(8)
where is a constant. The constant is referred to as a zero of the function. The constant is referred to as a pole of the function. We can see that when then the function becomes zero, and when the function value approaches infinity (for real ). The poles and zeros can be real or complex. Now, assuming that the poles are distinct, we can write the function as follows:
(9)
Our goal is to develop the procedure for computing the constants . Let us multiply the last expression with . As a result, we obtain
(10)
Now by substituting in the expression , we obtain
(11)
That is, we have obtained an expression for computing the constant . To compute the constant corresponding to the -th pole, we multiply the unfactorized expression by and substitute the value for .
Why this partial fraction expansion is important? Namely, after we have computed (9), we can simply compute the inverse Laplace transform by noting that:
(12)
where is a Heaviside step function. We can do this for all other fractions to obtain the time-domain transformation.
Let us do an example.
Example 1: Compute the step response of the function
(13)
Solution: The step response of the function in the complex domain is obtained by multiplying the transfer function by the term . This term corresponds to the Laplace transform of the Heaviside step function. That is, we want to compute the inverse Laplace transform of
(14)
Let us now apply the previously explained approach. We want to expand as follows:
(15)
By using the formula (11), we obtain:
(16)
(17)
(18)
Consequently, our factorized expression becomes
(19)
By computing the inverse Laplace transform of the last expression we obtain
(20)
The following MATLAB code is used to compute the partial fraction expansion and the inverse Laplace transform.
clc
syms s
F=(s+1)*(s+3)/(s*(s+2)*(s+4))
Factorized=partfrac(F)
ilaplace(Factorized)
The following output is generated.
F =
((s + 1)*(s + 3))/(s*(s + 2)*(s + 4))
Factorized =
1/(4*(s + 2)) + 3/(8*(s + 4)) + 3/(8*s)
ans =
exp(-2*t)/4 + (3*exp(-4*t))/8 + 3/8
Notice that we used the MATLAB function “ilaplace()” to compute the inverse Laplace transform.