Control Exp 03
Control Exp 03
Control Exp 03
PRELAB WORK:
Read this laboratory manual carefully before coming to the
laboratory class, so that you know what is required.
Try to follow the lecture notes of EEE 363.
Familiarize yourself with relevant MATLAB functions and codes necessary for
this experiment.
Do not bring any prepared MATLAB code in the lab with any portable device.
DONOT copy others blindly!!!
Submit your lab report before the roll call.
Key MATLAB commands used in this experiment: feedback, rlocfind, rlocus, sgrid,
step
THEORY WITH EXAMPLES:
Closed-Loop Poles
The root locus of an (open-loop) transfer function H(s) is a plot of the locations (locus) of all
possible closed loop poles with proportional gain k and unity feedback:
and thus the poles of the closed loop system are values of s such that 1 + K H(s) = 0.
If we write H(s) = b(s)/a(s), then this equation have the form:
Let n = order of a(s) and m = order of b(s) [the order of a polynomial is the highest power of
s that appears in it].
We will consider all positive values of k. In the limit as k -> 0, the poles of the closed-loop
system are a(s) = 0 or the poles of H(s). In the limit as k -> infinity, the poles of the closedloop system are b(s) = 0 or the zeros of H(s).
No matter what we pick k to be, the closed-loop system must always have n poles, where n
is the number of poles of H(s). The root locus must have n branches, each branch starts at a
pole of H(s) and goes to a zero of H(s). If H(s) has more poles than zeros (as is often the
case), m < n and we say that H(s) has zeros at infinity. In this case, the limit of H(s) as s ->
infinity is zero. The number of zeros at infinity is n-m, the number of poles minus the number
of zeros, and is the number of branches of the root locus that go to infinity (asymptotes).
Since the root locus is actually the locations of all possible closed loop poles, from the root
locus we can select a gain such that our closed-loop system will perform the way we want. If
any of the selected poles are on the right half plane, the closed-loop system will be unstable.
The poles that are closest to the imaginary axis have the greatest influence on the closed-loop
response, so even though the system has three or four poles, it may still act like a second or
even first order system depending on the location(s) of the dominant pole(s).
Plotting the root locus of a transfer function
Consider an open loop system which has a transfer function of
How do we design a feed-back controller for the system by using the root locus method? Say
our design criteria are 5% overshoot and 1 second rise time. Make a MATLAB file called
rl.m. Enter the transfer function, and the command to plot the root locus:
num=[1 7];
den=conv(conv([1 0],[1 5]),conv([1 15],[1 20]));
sys=tf(num,den);
rlocus(sys)
axis([-22 3 -15 15])
On the plot above, the two dotted lines at about a 45 degree angle indicate pole locations with
Zeta = 0.7; in between these lines, the poles will have Zeta > 0.7 and outside of the lines Zeta
< 0.7. The semicircle indicates pole locations with a natural frequency Wn = 1.8; inside the
Note that since the root locus may has more than one branch, when you select a pole, you
may want to find out where the other pole (poles) are. Remember they will affect the
response too. From the plot above we see that all the poles selected (all the "+" signs) are at
reasonable positions. We can go ahead and use the chosen k as our proportional controller.
Closed-loop response
In order to find the step response, you need to know the closed-loop transfer function. You
could compute this using the rules of block diagrams, or let MATLAB do it for you:
sys_cl= feedback(k*sys,1)
The two arguments to the function feedback are the numerator and denominator of the openloop system. You need to include the proportional gain that you have chosen. Unity feedback
is assumed.
If you have a non-unity feedback situation, look at the help file for the MATLAB function
feedback, which can find the closed-loop transfer function with a gain in the feedback loop.
As we expected, this response has an overshoot less than 5% and a rise time less than 1
second.
A Realtime Example: Cruise Control Problem Using Root Locus Method
The open-loop transfer function for this problem is:
where
m=1000
b=50
U(s)=10
Y(s)=velocity output
proportional controller (Kp) will be considered to solve this problem. The closed-loop
transfer function becomes:
Also, from the Root-Locus Tutorial, we know that the MATLAB command called sgrid
should be used to find an acceptable region of the root-locus plot. To use the sgrid, both the
damping ratio (zeta) and the natural frequency (Wn) need to be determined first. The
following two equations will be used to find the damping ratio and the natural frequency:
where
Wn=Natural frequency
zeta=Damping ratio
Tr=Rise time
Mp=Maximum overshoot
One of our design criteria is to have a rise time of less than 5 seconds. From the first
equation, we see that the natural frequency must be greater than 0.36. Also using the second
equation, we see that the damping ratio must be greater than 0.6, since the maximum
overshoot must be less than 10%.
Now, we are ready to generate a root-locus plot and use the sgrid to find an acceptable region
on the root-locus. Create a new m-file and enter the following commands.
m=1000;
b=50;
u=10;
num=[1];
den=[m b];
cruise=tf(num,den);
rlocus(cruise)
axis([-0.6 0 -0.6 0.6]);
sgrid(0.6, 0.36)
[Kp, poles]=rlocfind(cruise)
sys_cl=feedback(Kp*cruise,1);
t=0:0.1:20;
step(u*sys_cl,t)
axis ([0 20 0 10])
Running this m-file should give you the following root-locus plot.
The two dotted lines in an angle indicate the locations of constant damping ratio (zeta=0.6);
the damping ratio is greater than 0.6 in between these lines and less than 0.6 outside the lines.
The semi-ellipse indicates the locations of constant natural frequency (Wn=0.36); the natural
frequency is greater than 0.36 outside the semi-ellipse, and smaller than 0.36 inside.
If you look at the MATLAB command window, you should see a prompt asking you to pick a
point on the root-locus plot. Since you want to pick a point in between dotted lines (zeta>0.6)
and outside the semi-ellipse (Wn>0.36), click on the real axis just outside the semi-ellipse
(around -0.4).
You should see the gain value (Kp) and pole locations in the MATLAB command window.
Also you should see the closed-loop step response similar to the one shown below.
With the specified gain Kp (the one you just picked), the rise time and the overshoot criteria
have been met; however, a steady-state error of more than 10% remains.
Lag controller
To reduce the steady-state error, a lag controller will be added to the system. The transfer
function of the lag controller is:
If you read the "Lag or Phase-Lag Compensator using Root-Locus" section in Lead and Lag
Compensator page, the pole and the zero of a lag controller need to be placed close together.
Also, it states that the steady-state error will be reduce by a factor of Zo/Po. For these
reasons, let Zo equal -0.3 and Po equal -0.03.
Create a new m-file, and enter the following commands.
m = 1000;
b = 50;
u = 10;
Zo=0.3;
Po=0.03;
num=[1];
den=[m b];
cruise=tf(num,den);
contr=tf([1 Zo],[1 Po]);
rlocus(contr*cruise);
axis([-0.6 0 -0.4 0.4])
sgrid(0.6,0.36);
[Kp, poles]=rlocfind(contr*cruise);
sys_cl=feedback(Kp*contr*cruise,1);
t=0:0.1:20;
step(u*sys_cl,t)
axis ([0 20 0 12])
Running this m-file should give you a root-locus plot similar to the following:
In the MATLAB command window, you should see the prompt asking you to select a point
on the root-locus plot. Once again, click on the real axis around -0.4. You should have the
following response.
As you can see, the steady-state error has been reduced to near zero. The slight overshoot is a
result of the zero added in the lag controller.
Now all of the design criteria have been met and no further iterations are needed.
Home Task
The dynamic equations and the open-loop transfer function of the DC Motor are: