Then the transfer function of the overall system is
The closed loop poles are the roots of the denominator polynomial (called the characteristic polynomial).
Definition 30
All points on the root locus must satisfy Equation 15.
Root Locus Rules
Theorem 6
Next, because we are dealing with real systems, complex poles must have a corresponding conjugate pole.
Theorem 7
The root locus is symmetric about the real axis.
Going back to Equation 15, we can alternatively express the angle criteria as
Theorem 8
The real axis segments of the root locus are to the left of an odd number of open loop poles and zeros.
Theorem 9
The root locus begins at the open loop poles and ends at the open loop zeros
Theorem 10
In the limit, poles will asymptotically approach
and the real axis intercept of these asymptotes is
If there is a gap between real-axis segments, in order to end at an open loop zero, poles must sometimes break away from the real axis and then re-enter.
Theorem 11
The break-in and break-away points satisfy the equation
Since the angles can travel asymptotically, they sometimes cross the imaginary axis.
Theorem 12
The root locus intersects the imaginary axis at points where
Similarly, if the poles begin at complex locations, then we can find their angle of departure.
Theorem 13
Theorem 14
Generalized Root Locus
Thus if we apply the root locus rules to the open loop system
then we can capture the behavior of the closed loop poles of the original system as we vary the location of the open loop pole we control.
Bode Plots
Definition 31
A Bode plot is a plot of the magnitude and phase of the frequency response with the magnitude on a log-log scale and the phase on a semi-log scale.
If we write the frequency response in polar form,
If we take the convert this to decibels, we get
One useful way to use bode plots is to approximate the gain and phase margin because they can easily be seen visually from the plots themselves.
Nyquist Criteria
Consider the basic feedback system in Figure 8 and suppose that
Then the feedback transfer function is
If we focus specifically on the poles of the system, then we see
Definition 32
A contour is a closed loop set of points in the complex plane.
Definition 33
A mapping is a function that takes a point in the complex plane and transforms it into another point in the complex plane.
Definition 34
The contour which is easiest to find the mapping for is the one which starts at the origin, travels up the imaginary axis, encapsulates the right half plane, and then travels back up the imaginary axis back to the origin in a counter-clockwise fasion. This is the easiest because while the contour is on the imaginary axis, the mapping is just the frequency response of the system, and we can use the Bode plot of the system in order to draw the contour because each point on the mapping is a complex vector, and the bode plot can give us both the magnitude and angle of that vector.
Suppose we choose to control the plant by scaling our error signal by K and then put the controller on the feedback path like in Figure 7. It would be helpful for us to understand how the closed loop poles of the feedback system change as K is varied over the range [0,∞). First, suppose
Clearly, no matter what K is, the poles must satisfy two criteria.
kDGDHNGNH=1∠KDGDHNGNH=(2r+1)π,r∈Z(15)
The root locus is the set of all s∈C such that ∃K where Δ(s)=0.
First, notice that the roots of Δ(s) are the closed loop poles of the system.
The number of branches in the root locus is equal to the number of closed loop poles where a branch is the path traveled by a single pole as Kis varied.
where zi are open loop zeros and pi are open loop poles. If we restrict ourselves to the real axis, then given a closed loop pole s, each zi>s will contribute −180˚and each pi>s will contribute 180˚while the zi,pi<s will contribute 0˚.
When K is small, then the poles look like the open loop poles. As K grows very large, then the poles look like the open loop zeros.
If there are more poles (n) then zeros (m), then not all of the poles will end up at a zero in the limit. This means that n−m poles must branch off to infinity.
Poles beginning at complex locations will depart at an angle θ where
∑i=1m(p+zi)−∑i=1n(p+pi)=(2l+1)π
Finally, since ∣KG(s)H(s)∣=1, we can determine K if we know a particular pole location.
Given a pole location p,
K=G(p)H(p)1.
Because the Root Locus rules are derived from the characteristic polynomial Δ(s) of the closed-loop system, they can be used not just to find how the closed loop poles vary with a gain, but also to find how the closed loop poles vary with an open loop pole. Suppose that G(s)=(s+k)∏i(s+pi)NG and H(s)=1. Then
Each r is the magnitude of a factor 1+ωnjω where ωn is either a zero or a pole, zi,pk are the phases of each factor, and Nz0,Np0 are the number of zeros and poles at 0. By writing G(ω) this way, it is clear that
Likewise, the exponential form of G(ω) tells us that
∠G(ω)=2π(Nz0−Np0)+(∑i=0nzi−∑k=0mpk).
Each pk and zi are of the form 1+ωnjω. If ω>10ωn, then pk,zi≈ωn. Likewise, if ω<10ωn, pk,zi≈1. This means we can approximate bode plots using piece-wise linear segments using the following rules.
Each zero ωz contributes 20dB/dec to magnitude starting at ωz.
Each pole ωp contributes −20dB/dec to magnitude starting at ωp.
Each zero ωz contributes 45˚/dec to phase starting at 10ωz and ending at 10ωz.
Each pole ωp contributes −45˚/dec to phase starting at 10ωp and ending at 10ωp.
G(s)=DGNGH(s)=DHNH.
R(s)Y(s)=1+G(s)H(s)G=DGDH+NGNHNGDH.
1+GH=1+DGDHNGNH=DGDHNGNH+DGDH
From here, we can see that the poles of 1+GH are the poles of the open loop system whereas the zeros of 1+GH are the poles of the closed loop system.
The Nyquist Criterion says that if N is the number of counter-clockwise encirclements of zero of a contour mapped by a transfer function F(s), P is the number of poles in the contour, and Z is the number of zeros in the contour, then
N=P−Z.(16)
Thus, given an open loop transfer function GH, we can determine its stability. We already know P from the poles of the open loop system, and we can find N by defining a contour which encapsulates the right half plane and use Equation 16 to find Z. However, remember that we need to find the RHP poles of 1+GH. This shifts our mapping to the right by 1, so we can instead just let N be the number of encirclements of −1. Once we have Z, we know how many RHP poles the closed loop transfer function will have because they are the same as the RHP zeros of the 1+GH. We can extend the Nyquist Criterion to finding a range of gains that would make the open-loop system kG(s)H(s) stable by looking for the encirclements of K−1.