Control System Analysis, The Nyquist Criteria

Stability via the Nyquist Diagram – The Nyquist Criteria

The Nyquist criterion can tell us if the system is stable or unstable by determining how many closed-loop poles are in the right half-plane of the closed-loop system of Figure 1:

null

Figure 1

Consider the contour A defined in s-plane of Figure 2:

null

Figure 2

If a contour, A, that encircles the entire right half-plane of the root-locus of the system determined by the characteristic equation 1+ G(s)H(s), is mapped through G(s)H(s), then the number of closed-loop poles, Z, in the right half-plane equals the number of open-loop poles, P, that are in the right half-plane minus the number of counterclockwise revolutions, N, around -1 of the mapping; that is, Z:

null

Thus, to reach stability, Z must be equal to zero.

This mapping is called the Nyquist diagram, or Nyquist plot, of G(s)H(s).

To understand the Nyquist criteria for stability, we must first establish four important concepts that will be used during its application:

(1) the relationship between the poles of 1+ G(s)H(s) and the poles of G(s)H(s); (2) the relationship between the zeros of 1+ G(s)H(s) and the poles of the closed-loop transfer function (3) the concept of mapping points; and (4) the concept of mapping contours.

We could demonstrate that the poles of 1+ G(s)H(s) are the same as the
poles of G(s)H(s), the open-loop system, and (2) the zeros of  1+ G(s)H(s) are the
same as the poles of closed-loop transfer function of the system.

Example:

null

Step 1. Find the open-loop transfer function G(s)H(s) of the system.

Consider the closed-loop control system as follows:

null

The system characteristic equation is as follows:

null

The factor form of this characteristic equation is:

null

To determine the previous factor form:

null

Where the open-loop transfer function G(s)H(s) of the system is:

null

Step 2. Use Command Window of Matlab to draw the Nyquist Diagram, applying the following commands:

>> s=tf(‘s’);

>> G=10/(s^3+2*s^2+5*s);

>> nyquist(G);

null

null

We can see at the previous Diagram that for:

null

To reach stability, Z must be equal to zero:

null

Recalling that the poles of 1+ G(s)H(s), are the same as the poles of G(s)H(s), the open-loop system, we can determine P, the number of open-loop poles enclosed by the contour A from:

null

null

A detour around the poles on the contour is required:

null

In the Nyquist Diagram obtained for the system of Task 2, the point -1+j0 is highlighted in red:

null

We can see that N=0, so:

null

However, the Nyquist diagram intersects the real axis at -1+j0. Hence, according to the Nyquist Criteria, the system is marginally stable.

Sources:

  1. Modern_Control_Engineering, Ogata 4t
  2. Control Systems Engineering, Nise
  3. Sistemas de Control Automatico, Kuo

Literature review by:

Prof. Larry Francis Obando – Technical Specialist – Educational Content Writer

Resolving problems!!

WhatsApp:  +34633129287  Immediate Attention!!

Twitter: @dademuch

Copywriting, Content Marketing, Tesis, Monografías, Paper Académicos, White Papers (Español – Inglés)

Escuela de Ingeniería Electrónica de la Universidad Simón Bolívar, USB Valle de Sartenejas.

Escuela de Ingeniería Eléctrica de la Universidad Central de Venezuela, UCV CCs

Escuela de Turismo de la Universidad Simón Bolívar, Núcleo Litoral.

Contacto: España. +34633129287

Caracas, Quito, Guayaquil, Cuenca. 

WhatsApp:  +34633129287   +593998524011  

Twitter: @dademuch

FACEBOOK: DademuchConnection

email: dademuchconnection@gmail.com

Control System Analysis, The Nyquist Criteria

Sketching the Nyquist Diagram

The Nyquist diagram is also known as “The Polar Trace” of a transfer function G(jω), is a graph of the magnitude of G(jω) against the phase angle of G(jω) in polar coordinates, according to ω variables of zero to infinity. Therefore, The Nyquist diagram is the geometric place of vectors:

null
conforming ω modified from zero to infinity. Note that, in polar graphs, the phase angles are positive (negative) if they are measured counterclockwise (clockwise) from the positive real axis.

The following Figure shows an example of a Nyquist Diagram:

null

All points of the polar trace of G(jω) represent the end point of a vector at a given value of ω. The projections of G(jω) on the real and imaginary axes are their real and imaginary components. The magnitude and phase angle of G(jω) must be calculated directly for each frequency ω in order to construct polar traces.

Conceptually, the Nyquist diagram is plotted by substituting the points
of the contour that encloses the right half-plane into the function G(s)H(s).  This process is called mapping. Next Figure shows the process of mapping:

null

Consider the closed-loop control system of Figure 1:

null

Figure 1

Thus, in the Nyquist diagram, the contour that encloses the right half-plane, shown in Figure 2, can be mapped through the function G(s)H(s), derived from Figure 1,  by substituting points along the contour into G(s)H(s):

null

Figure 2

If a contour, A, that encircles the entire right half-plane of the root-locus of the system determined by the characteristic equation 1+ G(s)H(s), is mapped through G(s)H(s), then the number of closed-loop poles, Z, in the right half-plane equals the number of open-loop poles, P, that are in the right half-plane minus the number of counterclockwise revolutions, N, around -1+j0 of the mapping; that is, Z:

null

Thus, Z must be equal to zero to reach stability.

This mapping is called the Nyquist diagram, or Nyquist plot, of G(s)H(s). For more information an examples, see: The Nyquist Criteria

 

Example

Consider the control system whose block diagram and diagram are shown in the following Figure 3:

null

Figure 3

Conceptually, the Nyquist diagram is plotted by substituting the points of the contour shown in Figure 4(a) into G(s)H(s):

null

Each Pole and Zero term of G(s) shown in Figure 3(b) is a vector in Figure 4(a) and 4(b). The resultant vector, , found at any point along the contour is in general the product of the Zero vectors divided by the product of the Pole vectors (see Figure 4(c)). Thus, the magnitude of the resultant is the product of the Zero lengths divided by the product of the Pole lengths, and the angle of the resultant is the sum of the Zero angles minus the sum of the Pole angles.

null

Figure 4

The mapping from point A to point C can also be explained analytically. From
A to C the collection of points along the contour is imaginary. Hence, from A to C,
G(s)H(s)=G(s)*1=G(s)=G(jω), or from Figure 3(b):

nullAt zero frequency:

null

Thus, the Nyquist diagram starts at 50/3 at an angle of 0°. As ω increases the real part remains positive, and the imaginary part remains negative.

At null the real part becomes negative. At null, the Nyquist diagram crosses the negative real axis since the imaginary term goes to zero. The real value at the axis crossing, point Q in Figure 4(c), is -0.874. Continuing toward , the real part is negative, and the imaginary part is positive. At infinite frequency:

nullor approximately zero at 90°.

Around the infinite semicircle from point C to point D shown in Figure 4(b), the vectors rotate clockwise, each by 180°. Hence, the resultant undergoes a counterclockwise rotation of 3×180, starting at point C’ and ending at point D’ of Figure 4(c).

Nyquist diagram with Matlab

Consider the following open loop transfer function:

null

To create the Nyquist Diagram of the system, use the following commands in the command window of Matlab:

>> s=tf(‘s’)

>> G=1/(s^2+0.8*s+1)

>> nyquist(G)

This line of commands yields:

null

We can obtain information of points of interest in the Nyquist Diagram by cliking once over the contour. This yields:

null

Sources:

  1. Modern_Control_Engineering, Ogata 4t
  2. Control Systems Engineering, Nise
  3. Sistemas de Control Automatico, Kuo

Literature review by:

Prof. Larry Francis Obando – Technical Specialist – Educational Content Writer

I solve problems!!

WhatsApp:  +34633129287  Immediate Attention!!

Twitter: @dademuch

Copywriting, Content Marketing, Tesis, Monografías, Paper Académicos, White Papers (Español – Inglés)

Escuela de Ingeniería Electrónica de la Universidad Simón Bolívar, USB Valle de Sartenejas.

Escuela de Ingeniería Eléctrica de la Universidad Central de Venezuela, UCV CCs

Escuela de Turismo de la Universidad Simón Bolívar, Núcleo Litoral.

Contacto: España. +34633129287

Caracas, Quito, Guayaquil, Cuenca. 

WhatsApp:  +34633129287   +593998524011  

Twitter: @dademuch

FACEBOOK: DademuchConnection

email: dademuchconnection@gmail.com

 

 

 

Some content on this page was disabled on 3 April, 2020 as a result of a DMCA takedown notice from Pearson Education, Inc.. You can learn more about the DMCA here:

https://en.support.wordpress.com/copyright-and-the-dmca/