Nyquist Plot

From ControlTheoryPro.com

Revision as of 11:12, 18 April 2013 by Gabe Spradlin (Talk | contribs) (1 revision)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search
Nyquist Plot
Green carrot left.gif
Classical Control Nichols Plot
Green carrot.jpg
In order to prevent spam, users must register before they can edit or create articles.

1 Introduction to Nyquist Stability Criteria

The Nyquist Stability Criteria is a test for system stability, just like the Routh-Hurwitz test, or the Root Locus. However, the Nyquist Criteria can also give us additional information about a system. The Nyquist Criteria can provide information about the frequency characteristics of the system. For instance, some systems with constant gain might be stable for low-frequency inputs, but become unstable for high-frequency inputs. Finally, the Nyquist Criteria can tell us things about the phase of the input signals and the time-shift of the system.

2 Contours

A contour is a complicated mathematical construct, but luckily we only need to worry ourselves with a few points about them. We will denote contours with the greek letter Γ (gamma). Contours are lines, drawn on a graph, that follow certain rules:

  1. The contour must close (it must form a complete loop).
  2. The contour may not cross directly through a pole of the system.
  3. Contours must have a direction (clockwise or counterclockwise, generally).
  4. A contour is called "simple" if it has no self-intersections. We only consider simple contours here.

Once we have such a contour, we can develop some important theorems about them, and finally use these theorems to derive the Nyquist stability criterion.

3 Argument Principal

If we have a contour, Γ, drawn in one plane (say the complex laplace plane, for instance), we can map that contour into another plane, the F(s) plane, by transforming the contour with the function F(s). The resultant contour, LaTeX: \Gamma_{F(s)} will circle the origin point of the F(s) plane N times, where N is equal to the difference between Z and P (the number of zeros and poles of the function F(s), respectively).

When we have our contour, Γ, we transform it into LaTeX: \Gamma_{F(s)} by plugging every point of the contour into the function F(s), and taking the resultant value to be a point on the transformed contour.

3.1 Example: First Order System

Let's say, for instance, that Γ is a unit square contour in the complex s plane. The verticies of the square are located at points E, F, G, H, as follows:

LaTeX: I = 1 + j1
LaTeX: J = 1 - j1
LaTeX: K = -1 - j1
LaTeX: L = -1 + j1

we must also specify the direction of our contour, and we will say (arbitrarily) that it is a clockwise contour (travels from I to J to K to L). We will also define our tranform function, F(s), to be the following:

LaTeX: F(s) = s2 + 1

We can factor the denominator of F(s), and we can show that there is one zero at s → -0.5, and no poles. Plotting this root on the same graph as our contour, we see clearly that it lies within the contour. Since s is a complex variable, defined with real and imagnary parts as:

LaTeX: s = \sigma + j \omega

We know that F(s) must also be complex. We will say, for reasons of simplicy, that the axes in the F(s) plane are u and v, and are related as such:

LaTeX: F(s) = u + jv = 2(\sigma + j\omega) + 1

From this relationship, we can define u and v in terms of σ and ω:

LaTeX: u = 2 \sigma + 1
LaTeX: v = 2 \omega

Now, to transform Γ, we will plug every point of the contour into F(s), and the resultant values will be the points of LaTeX: \Gamma_{F(s)}. We will solve for complex values u and v, and we will start with the verticies, because they are the simplest examples:

LaTeX: u + jv = F(I) = 3 + j2
LaTeX: u + jv = F(J) = 3 - j2
LaTeX: u + jv = F(K) = -1 + j2
LaTeX: u + jv = F(L) = -1 - j2

We can take the lines in between the vertices as a function of s, and plug the entire function into the transform. Luckily, because we are using straight lines, we can simplify very much:

  • Line from I to J: LaTeX: \sigma = 1, u = 3, v = \omega
  • Line from J to K: LaTeX: \omega = -1, u = 2 \sigma + 1, v = -1
  • Line from K to L: LaTeX: \sigma = -1, u = -1, v = \omega
  • Line from L to I: LaTeX: \omega = 1, u = 2 \sigma + 1, v = 1

And when we graph these functions, from virtex to virtex, we see that the resultant contour in the F(s) plane is a square, but not centered at the origin, and larger in size. Notice how the contour encircles the origin of the F(s) plane one time. This will be important later on.

4 Nyquist Contour

The Nyquist contour, the contour that makes the entire nyquist criterion work, must encircle the entire unstable region of the complex plane. For analog systems, this is the right half of the complex s plane. For digital systems, this is the entire plane outside the unit circle. Remember that if a pole to the closed-loop transfer function (or equivalently a zero of the characteristic equation) lies in the unstable region of the complex plane, the system is an unstable system.

Analog Systems
The Nyquist contour for analog systems is an infinite semi-circle that encircles the entire right-half of the s plane. The semicircle travels up the imaginary axis from negative infinity to positive infinity. From positive infinity, the contour breaks away from the imaginary axis, in the clock-wise direction, and forms a giant semicircle.
Digital Systems
The Nyquist contour in digital systems is a counter-clockwise encirclement of the unit circle.

5 Nyquist Criteria

Let us first introduce the most important equation when dealing with the Nyquist criterion:

LaTeX: N = Z - P


LaTeX: N| is the number of encirclements of the (-1, 0) point.
LaTeX: Z| is the number of zeros of the characteristic equation.
LaTeX: P| is the number of poles of the open-loop characteristic equation.

With this equation stated, we can now state the Nyquist Stability Criterion: A feedback control system is stable, if and only if the contour LaTeX: \Gamma_{F(s)} in the F(s) plane does not encircle the (-1, 0) point when P is 0.

A feedback control system is stable, if and only if the contour LaTeX: \Gamma_{F(s)} in the F(s) plane encircles the (-1, 0) point a number of times equal to the number of poles of F(s) enclosed by Γ.

In other words, if P is zero then N must equal zero. Otherwise, N must equal P. Essentially, we are saying that Z must always equal zero, because Z is the number of zeros of the characteristic equation (and therefore the number of poles of the closed-loop transfer function) that are in the right-half of the s plane.

Keep in mind that we don't necessarily know the locations of all the zeros of the characteristic equation. So if we find, using the nyquist criterion, that the number of poles is not equal to N, then we know that there must be a zero in the right-half plane, and that therefore the system is unstable.

6 Nyquist ↔ Bode

A careful inspection of the Nyquist plot will reveal a surprising relationship to the Bode plots of the system. If we use the Bode phase plot as the angle θ, and the Bode magnitude plot as the distance r, then it becomes apparent that the Nyquist plot of a system is simply the polar representation of the Bode plots.

To obtain the Nyquist plot from the Bode plots, we take the phase angle and the magnitude value at each frequency ω. We convert the magnitude value from decibels back into gain ratios. Then, we plot the ordered pairs (r, θ) on a polar graph.

7 Nyquist in the Z Domain

The Nyquist Criteria can be utilized in the digital domain in a similar manner as it is used with analog systems. The primary difference in using the criteria is that the shape of the Nyquist contour must change to encompass the unstable region of the Z plane. Therefore, instead of an infinitesimal semi-circle, the Nyquist contour for digital systems is a counter-clockwise unit circle. By changing the shape of the contour, the same N = P - Z equation holds true, and the resulting Nyquist graph will typically look identical to one from an analog system, and can be interpreted in the same way.