Transient Response of a Prototype Second-Order Dynamical System with Detailed Derivation Using Inverse Laplace Transform (damping, natural frequency, and damping ratio) – Fusion of Engineering, Control, Coding, Machine Learning, and Science

In this post, we explain how to analytically compute a transient response of a prototype second-order linear dynamical system. We start from a generic (prototypical) form of the second-order system, and we explain how to compute its step response (transient response) by using the inverse Laplace transform. The derivation and analysis presented in this post are very useful for understanding the transient response of linear dynamical systems and how different parameters, such as damping, natural frequency, etc. influence the transient response behavior. The YouTube video accompanying this post is given here

The prototype second-order system has the following form:

(1) $\begin{align*}W(s)=\frac{\omega_{n}^{2}}{s^2+2\zeta \omega_{n} s +\omega_{n}^{2}}\end{align*}$

where

$s$ is the complex variable
$\omega_{n}$ is the natural undamped frequency
$\zeta$ is the damping ratio

In this post, we derive the expression for the step response of the system (1). The step response is given by the following equation

(2) $\begin{align*}y(t)& =1- \frac{e^{-\zeta \omega_{n} t}}{\sqrt{1-\zeta^2}} \sin((\omega_{n} \sqrt{1-\zeta^2}) t +\cos^{-1}(\zeta))\end{align*}$

This equation can also be written as follows

(3) $\begin{align*}y(t)& =1-e^{-\zeta \omega_{n} t}(\frac{\zeta}{\sqrt{1-\zeta^2}} \sin((\omega_{n} \sqrt{1-\zeta^2}) t)+\cos((\omega_{n} \sqrt{1-\zeta^2})t) ) \\\end{align*}$

The general form of the second-order transfer function given by (1) is often used in control systems textbooks and articles. Moreover, many properties of the dynamical systems are often expressed in terms of the undamped natural frequency $w_{n}$ and the damping ratio $\zeta$ . Consequently, before we proceed with the derivation of the step response of the second-order system (1), we relate the poles of this system with the natural undamped frequency and the damping ratio.

The poles of the dynamical system (1), are computed as follows:

(4) $\begin{align*} s^2+2\zeta \omega_{n} s +\omega_{n}^{2}=0 \end{align*}$

(5) $\begin{align*}s_{1,2}&=\frac{-2\zeta \omega_{n}\pm \sqrt{4\zeta^2 \omega_{n}^{2}-4\omega_{n}^{2}}}{2} \\s_{1,2}&=\frac{-2\zeta \omega_{n}\pm 2\omega_{n}\sqrt{\zeta^2-1}}{2} \\ s_{1,2}&=\frac{-2\zeta \omega_{n}\pm 2\omega_{n}\sqrt{-1(1-\zeta^2)}}{2} \\s_{1,2}&=-\zeta \omega_{n}\pm j\omega_{n}\sqrt{1-\zeta^2} \\s_{1,2}&=-\alpha \pm j \omega \\\end{align*}$

where $j$ is the imaginary unit and

(6) $\begin{align*}\alpha&=\zeta \omega_{n} \\\omega&=\omega_{n}\sqrt{1-\zeta^2}\end{align*}$

These poles have a very useful graphical interpretation given in the figure below.

Figure 1: Graphical interpretation of the poles of the prototype system (1).

Here is an interpretation of Fig. 1:

The natural undamped frequency $\omega_{n}$ is the distance from the pole to the origin.
The damping ratio $\zeta$ can be interpreted as the cosine of the angle $\theta$ .
The parameter $\alpha=\zeta \omega_{n}$ is called the damping factor or the damping constant.
When $\zeta=1$ , the system is critically damped. In this case, $\alpha_{c}=\zeta\omega_{n}=\omega_{n}$ . This enables us to define the damping ratio as:
(7) $\begin{align*}\zeta=\frac{\text{actual damping factor}}{\text{damping factor at critical damping}}=\frac{\alpha}{\omega_{n}}\end{align*}$
The parameter $\omega$ is called the damped frequency or the conditional frequency.

In the sequel, we derive an analytical expression for the step response of the system (1). The output of the system can be represented as follows:

(8) $\begin{align*}Y(s)=W(s)U(s)\end{align*}$

Since we are interested in the step response, we assume $U(s)=1/s$ . By substituting this expression in (8), we obtain

(9) $\begin{align*}Y(s)=\frac{\omega_{n}^{2}}{s^2+2\zeta \omega_{n} s +\omega_{n}^{2}}\cdot \frac{1}{s} \end{align*}$

We use partial fraction expansion to compute the response. The partial fraction expansion of (9) is given by

(10) $\begin{align*}Y(s)=\frac{K_{1}}{s}+\frac{K_{2}}{s+\alpha -\omega j}+\frac{K_{3}}{s+\alpha +\omega j}\end{align*}$

We first need to determine the constants $K_{1}$ , $K_{2}$ , and $K_{3}$ , and after that, we can compute the inverse Laplace transform to compute the system response in the time domain. The partial fraction expansion and the cover method for computing it is explained in our previous post which can be found here. By applying the cover method, we obtain the value for the constant $K_{1}$

(11) $\begin{align*}K_{1}=Y(s)s\vert_{s=0}=\frac{\omega_{n}^{2}}{s^2+2\zeta \omega_{n} s +\omega_{n}^{2}}\vert_{s=0}=1 \end{align*}$

The constant $K_{2}$ is computed as follows

(12) $\begin{align*}K_{2}=Y(s)(s+\alpha-\omega j ) \vert_{s=-\alpha +\omega j} =\frac{\omega_{n}^{2}(s+\alpha-\omega j ) }{(s+\alpha-\omega j ) (s+\alpha+\omega j ) s} \vert_{s=-\alpha +\omega j} \end{align*}$

Here, we used the following expression

(13) $\begin{align*}s^2+2\zeta \omega_{n} s +\omega_{n}^{2}=(s+\alpha-\omega j ) (s+\alpha+\omega j )\end{align*}$

which originates from the following fact

(14) $\begin{align*}Y(s)=\frac{\omega_{n}^{2}}{(s^2+2\zeta \omega_{n} s +\omega_{n}^{2})s}=\frac{\omega_{n}^{2}}{(s+\alpha-\omega j ) (s+\alpha+\omega j )s}\end{align*}$

Going back to (12), we obtain

(15) $\begin{align*}K_{2}& =\frac{\omega_{n}^{2}}{(s+\alpha+\omega j ) s} \vert_{s=-\alpha +\omega j} = \frac{\omega_{n}^{2}}{ (-\alpha +\omega j +\alpha+\omega j ) (-\alpha +\omega j)} = \\& =\frac{\omega_{n}^{2}}{ 2\omega j (-\alpha +\omega j)}= \frac{\omega_{n}^{2}}{ 2 j (-\omega\alpha +\omega^2 j)}\end{align*}$

Next, we substitute the expressions for $\alpha$ and $\omega$ from (6) in (15), and as a result, we obtain

(16) $\begin{align*}K_{2} & =\frac{\omega_{n}^{2}}{2j (-\zeta \omega_{n}^{2}\sqrt{1-\zeta^2}+\omega_{n}^2(1-\zeta^2)j)} \\ & = \frac{\omega_{n}^{2}}{2j\omega_{n}^{2} (-\zeta \sqrt{1-\zeta^2}+(1-\zeta^2)j)} \\& = \frac{1}{2j\sqrt{1-\zeta^2} (-\zeta +j\sqrt{1-\zeta^2})} \end{align*}$

Next, let us write the complex number $(-\zeta +j\sqrt{1-\zeta^2})$ in the polar form

(17) $\begin{align*}-\zeta +j\sqrt{1-\zeta^2} & = \sqrt{(-\zeta)^{2}+(\sqrt{1-\zeta^2})^{2}}e^{j\arctan(\frac{\sqrt{1-\zeta^2}}{-\zeta})} \\& = 1\cdot e^{j\arctan(\frac{\sqrt{1-\zeta^2}}{-\zeta})} \\& = 1\cdot e^{j\arctan(\frac{\omega_{n}\sqrt{1-\zeta^2}}{-\omega_{n}\zeta})} \\& = 1\cdot e^{j\arctan(\frac{\omega }{-\alpha})}\end{align*}$

On the other hand, from Figure 1, we have that

(18) $\begin{align*}\arctan(\frac{\omega }{-\alpha})= \pi -\theta=\phi\end{align*}$

By substituting (18) in (17), we obtain

(19) $\begin{align*}-\zeta +j\sqrt{1-\zeta^2} =1\cdot e^{j \phi}, \;\; \text{where} \;\; \phi= \pi -\theta \end{align*}$

Substituting (19) in (16), we obtain

(20) $\begin{align*} K_{2}=\frac{1}{2j\sqrt{1-\zeta^2} ( e^{j \phi} )} = \frac{e^{- j \phi} }{2j\sqrt{1-\zeta^2}}\end{align*}$

Next, we compute the constant $K_{3}$ by using the same strategy

(21) $\begin{align*}K_{3}& =Y(s)(s+\alpha+\omega j ) \vert_{s=-\alpha -\omega j} =\frac{\omega_{n}^{2}(s+\alpha+\omega j ) }{(s+\alpha-\omega j ) (s+\alpha+\omega j ) s} \vert_{s=-\alpha -\omega j} \\& = \frac{\omega_{n}^{2}}{2j (\omega \alpha +\omega^{2}j)}\end{align*}$

Next, by substituting (6) in (21), we obtain

(22) $\begin{align*}K_{3}=\frac{1}{2j\sqrt{1-\zeta^2} (\zeta +j\sqrt{1-\zeta^2})} \end{align*}$

By using the previously explained polar transformation of complex numbers, we obtain

(23) $\begin{align*}K_{3}=\frac{e^{-j\arctan(\theta)}}{2j\sqrt{1-\zeta^2}} \end{align*}$

Now, taking into account that $\phi=\pi - \theta$ , from (23), we obtain

(24) $\begin{align*}K_{3} & =\frac{e^{j(-\pi+\phi)}}{2j\sqrt{1-\zeta^2}} = \frac{e^{-j\pi} e^{j\phi}}{2j\sqrt{1-\zeta^2}} \\& = -\frac{e^{j\phi}}{2j\sqrt{1-\zeta^2}}\end{align*}$

By using (11), (20), and (24), we can write (10) as follows

(25) $\begin{align*}Y(s)=\frac{1}{s}+\frac{e^{- j \phi} }{2j\sqrt{1-\zeta^2}}\cdot \frac{1}{s+\alpha -\omega j}-\frac{e^{j\phi}}{2j\sqrt{1-\zeta^2}}\cdot \frac{1}{s+\alpha +\omega j}\end{align*}$

Now, we need to find an inverse Laplace transform of (25). We need the following inverse Laplace transforms

(26) $\begin{align*}& \mathcal{L}^{-1}(\frac{1}{s} ) = 1 \\& \mathcal{L}^{-1}(\frac{1}{s+a} ) = e^{-at},\;\; t\ge 0\end{align*}$

By using these inverse transforms, from (25}), we obtain

(27) $\begin{align*}y(t)& =1 + \frac{e^{- j \phi} }{2j\sqrt{1-\zeta^2}} e^{-(\alpha -\omega j)t} -\frac{e^{j\phi}}{2j\sqrt{1-\zeta^2}} e^{-(\alpha +\omega j)t} \\y(t)& =1+\frac{e^{ j (\omega t -\phi)} }{2j\sqrt{1-\zeta^2}} \cdot e^{-\alpha t} - \frac{e^{-j(\omega t - \phi)} }{2j\sqrt{1-\zeta^2}} e^{-\alpha t} \\y(t) & = 1 + \frac{e^{-\alpha t}}{\sqrt{1-\zeta^2}} (\frac{e^{ j (\omega t -\phi)} - e^{-j(\omega t - \phi)}}{2j}) \end{align*}$

Next, we use the following trigonometric identity

(28) $\begin{align*}\frac{e^{j\beta}-e^{-j\beta}}{2j} & =\frac{1}{2j} (\cos(\beta)+j\sin(\beta)-(\cos(\beta)-j\sin(\beta))) \\& =\sin(\beta )\end{align*}$

for some angle $\beta$ . By using this expression in (27), we obtain

(29) $\begin{align*}y(t)& =1 + \frac{e^{-\alpha t}}{\sqrt{1-\zeta^2}} \sin(\omega t - \phi)\\y(t)& =1 + \frac{e^{-\alpha t}}{\sqrt{1-\zeta^2}} \sin(\omega t - \pi +\theta)\\\end{align*}$

Next, we use the following trigonometric expression

(30) $\begin{align*}\sin(\beta-\pi)=-\sin(\beta)\end{align*}$

By using this expression in (31), we obtain

(31) $\begin{align*}y(t)& =1 - \frac{e^{-\alpha t}}{\sqrt{1-\zeta^2}} \sin(\omega t +\theta)\\y(t)& =1- \frac{e^{-\alpha t}}{\sqrt{1-\zeta^2}} \sin(\omega t +\cos^{-1}(\zeta)) \\y(t)& =1- \frac{e^{-\zeta \omega_{n} t}}{\sqrt{1-\zeta^2}} \sin((\omega_{n} \sqrt{1-\zeta^2}) t +\cos^{-1}(\zeta))\end{align*}$

That is, the final expression is given by

(32) $\begin{align*}y(t)& =1- \frac{e^{-\zeta \omega_{n} t}}{\sqrt{1-\zeta^2}} \sin((\omega_{n} \sqrt{1-\zeta^2}) t +\cos^{-1}(\zeta))\end{align*}$

Using the procedure explained in the post that can be found here, we can also write the equation (32) as follows

(33) $\begin{align*}y(t)& =1-e^{-\alpha t}(\frac{\zeta}{\sqrt{1-\zeta^2}} \sin(\omega t)+\cos(\omega t) ) \\y(t)& =1-e^{-\alpha t}(\frac{\alpha}{\omega} \sin((\omega_{n} \sqrt{1-\zeta^2}) t)+\cos((\omega_{n} \sqrt{1-\zeta^2})t) ) \end{align*}$

admin

You might also like

Plot Phase Portraits and State-Space Trajectories of Dynamical Systems in MATLAB

Easy Explanation of Batch Gradient Descent Method with Python Implementation

Solve Constrained Optimization Problems in Python by Using SciPy Library, minimize() Function and Trust-Region Constrained Method