Frequency Response Analysis



We have already discussed time response analysis of the control systems and the time domain specifications of the second order control systems. In this chapter, let us discuss the frequency response analysis of the control systems and the frequency domain specifications of the second order control systems.

What is Frequency Response?

The response of a system can be partitioned into both the transient response and the steady state response. We can find the transient response by using Fourier integrals. The steady state response of a system for an input sinusoidal signal is known as the frequency response. In this chapter, we will focus only on the steady state response.

If a sinusoidal signal is applied as an input to a Linear Time-Invariant (LTI) system, then it produces the steady state output, which is also a sinusoidal signal. The input and output sinusoidal signals have the same frequency, but different amplitudes and phase angles.

Let the input signal be −

$$r(t)=A\sin(\omega_0t)$$

The open loop transfer function will be −

$$G(s)=G(j\omega)$$

We can represent $G(j\omega)$ in terms of magnitude and phase as shown below.

$$G(j\omega)=|G(j\omega)| \angle G(j\omega)$$

Substitute, $\omega = \omega_0$ in the above equation.

$$G(j\omega_0)=|G(j\omega_0)| \angle G(j\omega_0)$$

The output signal is

$$c(t)=A|G(j\omega_0)|\sin(\omega_0t + \angle G(j\omega_0))$$

  • The amplitude of the output sinusoidal signal is obtained by multiplying the amplitude of the input sinusoidal signal and the magnitude of $G(j\omega)$ at $\omega = \omega_0$.

  • The phase of the output sinusoidal signal is obtained by adding the phase of the input sinusoidal signal and the phase of $G(j\omega)$ at $\omega = \omega_0$.

Where,

  • A is the amplitude of the input sinusoidal signal.

  • ω0 is angular frequency of the input sinusoidal signal.

We can write, angular frequency $\omega_0$ as shown below.

$$\omega_0=2\pi f_0$$

Here, $f_0$ is the frequency of the input sinusoidal signal. Similarly, you can follow the same procedure for closed loop control system.

Frequency Domain Specifications

The frequency domain specifications are resonant peak, resonant frequency and bandwidth.

Consider the transfer function of the second order closed loop control system as,

$$T(s)=\frac{C(s)}{R(s)}=\frac{\omega_n^2}{s^2+2\delta\omega_ns+\omega_n^2}$$

Substitute, $s = j\omega$ in the above equation.

$$T(j\omega)=\frac{\omega_n^2}{(j\omega)^2+2\delta\omega_n(j\omega)+\omega_n^2}$$

$$\Rightarrow T(j\omega)=\frac{\omega_n^2}{-\omega^2+2j\delta\omega\omega_n+\omega_n^2}=\frac{\omega_n^2}{\omega_n^2\left ( 1-\frac{\omega^2}{\omega_n^2}+\frac{2j\delta\omega}{\omega_n} \right )}$$

$$\Rightarrow T(j\omega)=\frac{1}{\left ( 1-\frac{\omega^2}{\omega_n^2} \right )+j\left ( \frac{2\delta\omega}{\omega_n} \right )}$$

Let, $\frac{\omega}{\omega_n}=u$ Substitute this value in the above equation.

$$T(j\omega)=\frac{1}{(1-u^2)+j(2\delta u)}$$

Magnitude of $T(j\omega)$ is -

$$M=|T(j\omega)|=\frac{1}{\sqrt {(1-u^2)^2+(2\delta u)^2}}$$

Phase of $T(j\omega)$ is -

$$\angle T(j\omega)=-tan^{-1}\left( \frac{2\delta u}{1-u^2} \right )$$

Resonant Frequency

It is the frequency at which the magnitude of the frequency response has peak value for the first time. It is denoted by $\omega_r$. At $\omega = \omega_r$, the first derivate of the magnitude of $T(j\omega)$ is zero.

Differentiate $M$ with respect to $u$.

$$\frac{\text{d}M}{\text{d}u}=-\frac{1}{2}\left [ (1-u^2)^2+(2\delta u)^2 \right ]^{\frac{-3}{2}} \left [2(1-u^2)(-2u)+2(2\delta u)(2\delta) \right ]$$

$$\Rightarrow \frac{\text{d}M}{\text{d}u}=-\frac{1}{2}\left [ (1-u^2)^2+(2\delta u)^2 \right ]^{\frac{-3}{2}} \left [4u(u^2-1 +2\delta^2) \right ]$$

Substitute, $u=u_r$ and $\frac{\text{d}M}{\text{d}u}==0$ in the above equation.

$$0=-\frac{1}{2}\left [ (1-u_r^2)^2+(2\delta u_r)^2 \right ]^{-\frac{3}{2}}\left [ 4u_r(u_r^2-1 +2\delta^2) \right ]$$

$$\Rightarrow 4u_r(u_r^2-1 +2\delta^2)=0$$

$$\Rightarrow u_r^2-1+2\delta^2=0$$

$$\Rightarrow u_r^2=1-2\delta^2$$

$$\Rightarrow u_r=\sqrt{1-2\delta^2}$$

Substitute, $u_r=\frac{\omega_r}{\omega_n}$ in the above equation.

$$\frac{\omega_r}{\omega_n}=\sqrt{1-2\delta^2}$$

$$\Rightarrow \omega_r=\omega_n \sqrt{1-2\delta^2}$$

Resonant Peak

It is the peak (maximum) value of the magnitude of $T(j\omega)$. It is denoted by $M_r$.

At $u = u_r$, the Magnitude of $T(j\omega)$ is -

$$M_r=\frac{1}{\sqrt{(1-u_r^2)^2+(2\delta u_r)^2}}$$

Substitute, $u_r = \sqrt{1 − 2\delta^2}$ and $1 − u_r^2 = 2\delta^2$ in the above equation.

$$M_r=\frac{1}{\sqrt{(2\delta^2)^2+(2\delta \sqrt{1-2\delta^2})^2}}$$

$$\Rightarrow M_r=\frac{1}{2\delta \sqrt {1-\delta^2}}$$

Resonant peak in frequency response corresponds to the peak overshoot in the time domain transient response for certain values of damping ratio $\delta$. So, the resonant peak and peak overshoot are correlated to each other.

Bandwidth

It is the range of frequencies over which, the magnitude of $T(j\omega)$ drops to 70.7% from its zero frequency value.

At $\omega = 0$, the value of $u$ will be zero.

Substitute, $u = 0$ in M.

$$M=\frac{1}{\sqrt {(1-0^2)^2+(2\delta(0))^2}}=1$$

Therefore, the magnitude of $T(j\omega)$ is one at $\omega = 0$.

At 3-dB frequency, the magnitude of $T(j\omega)$ will be 70.7% of magnitude of $T(j\omega)$ at $\omega = 0$.

i.e., at $\omega = \omega_B, M = 0.707(1) = \frac{1}{\sqrt{2}}$

$$\Rightarrow M=\frac{1}{\sqrt{2}}=\frac{1}{\sqrt{(1-u_b^2)^2+(2\delta u_b)^2}}$$

$$\Rightarrow 2=(1-u_b^2)^2+(2\delta)^2 u_b^2$$

Let, $u_b^2=x$

$$\Rightarrow 2=(1-x)^2+(2\delta)^2 x$$

$$\Rightarrow x^2+(4\delta^2-2)x-1=0$$

$$\Rightarrow x=\frac{-(4\delta^2 -2)\pm \sqrt{(4\delta^2-2)^2+4}}{2}$$

Consider only the positive value of x.

$$x=1-2\delta^2+\sqrt {(2\delta^2-1)^2+1}$$

$$\Rightarrow x=1-2\delta^2+\sqrt {(2-4\delta^2+4\delta^4)}$$

Substitute, $x=u_b^2=\frac{\omega_b^2}{\omega_n^2}$

$$\frac{\omega_b^2}{\omega_n^2}=1-2\delta^2+\sqrt {(2-4\delta^2+4\delta^4)}$$

$$\Rightarrow \omega_b=\omega_n \sqrt {1-2\delta^2+\sqrt {(2-4\delta^2+4\delta^4)}}$$

Bandwidth $\omega_b$ in the frequency response is inversely proportional to the rise time $t_r$ in the time domain transient response.

Advertisements