Microcontroller based applied digital control (D. Ibrahim, 2006)
.pdf
208 System Stability
Example 8.14
The loop transfer function of a unity feedback sampled data system is given by
G(z) |
= z2 |
0.368z + 0.264 |
. |
|
− 1.368z + 0.368 |
||
Draw the Bode diagram and analyse the stability of the system. Assume that T = 1 s.
Solution |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Using the transformation |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
z |
1 + (T /2)w |
|
1 + 0.5w |
, |
|
|
|
|||||||
|
|
|
|
= 1 |
− |
(T /2)w = |
1 |
− |
0.5w |
||||||||
|
|
|
|
|
|
|
|||||||||||
we get |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G(w) |
|
|
|
0.368(1 + 0.5w/1 − 0.5w) + 0.264 |
|
||||||||||||
= (1 |
+ 0.5w/1 |
− 0.5w)2 − 1.368(1 + 0.5w/1 − 0.5w) + 0.368 |
|||||||||||||||
or |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
G(w) |
= − |
0.0381(w − 2)(w + 12.14) |
. |
|
||||||||||
|
|
|
|
|
w(w |
+ |
0.924) |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
To obtain the frequency response, we can replace w with j v, giving |
|||||||||||||||||
|
|
|
G( j v) |
= − |
0.0381( j v − 2)( j v + 12.14) |
. |
|||||||||||
|
|
|
|
||||||||||||||
|
|
|
|
|
|
j v( j v |
+ |
0.924) |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||
The magnitude is then
|
|
|
|
|
|
|
|
√ |
|
|
|
|
√ |
|
|
|
|
|
|
|
| |
G( j v) |
| = |
0.0381 |
v2 + 22 |
|
v2 + 12.142 |
|
|
||||||||||
|
|
|
|
|
|
v√ |
|
|
|
|
|
|
|||||||
and |
|
|
|
|
|
v2 + 0.9242 |
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G( j v) |
= |
tan−1 |
v |
+ |
tan−1 |
v |
|
− |
90 |
− |
tan−1 |
v |
. |
|||||
|
|
|
|||||||||||||||||
|
|
|
|
2 |
12.14 |
|
|
|
0.924 |
|
|||||||||
The Bode diagram is shown in Figure 8.13. The system is stable with a gain margin of 5 dB and a phase margin of 26◦.
8.7 EXERCISES
1.Given below are the characteristic equations of some sampled data systems. Using Jury’s test, determine if the systems are stable.
(a)z2 − 1.8z + 0.72 = 0
(b)z2 − 0.5z + 1.2 = 0
(c)z3 − 2.1z2 + 2.0z − 0.5 = 0
(d)z3 − 2.3z2 + 1.61z − 0.32 = 0
2.The characteristic equation of a sampled data system is given by
(z − 0.5)(z2 − 0.5z + 1.2) = 0.
Determine the stability of the system.
210 System Stability
7. Check the stability of the transfer function
G(z) =
1
z3 + 2.7z2 + 1.5z + 0.2
using:
(a)Jury’s test;
(b)the Routh–Hurwitz criterion;
(c)the oot locus
8.Repeat Exercise 7 using the Bode diagram.
9.A process with the transfer function
K
s(1 + as)
is preceded by a zero-order hold and is connected to form a unity feedback sampled data system.
(a)Assuming the sampling time is T , derive an expression for the closed-loop transfer function of the system.
(b)Draw the root locus of the system and hence determine the value of K for which the system becomes marginally stable.
10.The open-loop transfer function of a sampled data system is given by
1
G(z) = . z3 − 3.1z2 + 3.1z − 1.1
The closed-loop system is formed by using a unity gain feedback. Use Jury’s criterion to determine the stability of the system.
11. Use the Bode diagram to determine the stability of the sampled data system given by
G(z) = |
|
|
z |
|
|
. |
(z |
− |
1)(z |
− |
0.6) |
||
|
|
|
|
|
12.Repeat Exercise 11 using the Nyquist criterion.
13.The open-loop transfer function of a sampled data system is given by
G(z) |
|
|
K (z − 0.6) |
. |
|||
= (z |
|
||||||
|
− |
0.8)(z |
− |
0.4) |
|
||
|
|
|
|
||||
(a)Plot the Bode diagram by calculating the frequency response, assuming K = 1.
(b)From the Bode diagram determine the phase margin and the gain margin.
(c)Find the value of K for marginal stability.
(d)If the system is marginally stable, determine the frequency of oscillation.
14.The block diagram of a closed-loop sampled data system is shown in Figure 8.15. Determine the range of K for stability by:
(a)finding the roots of the characteristic equation;
(b)using Jury’s test;
(c)using the Routh–Hurwitz criterion;
FURTHER READING |
211 |
+ |
|
|
|
|
K ( z + 0.4) |
|
y(s) |
|||
r(s) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(z − 0.6)(z − 0.8) |
||||||
|
|
− |
|
|
||||||
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
Figure 8.15 Closed-loop system for Exercise 14
(d)using the root locus;
(e)drawing the Bode diagram;
(f)drawing the Nyquist diagram.
Which method would you prefer in this exercise and why?
15.Explain the mapping between the s-plane and the simple w-plane. How are the frequency points mapped?
FURTHER READING
[Bode, 1945] |
Bode, H.W. Network Analysis and Feedback Amplifier Design, Van Nostrand, |
|
Princeton, NJ, 1945. |
[D’Azzo and Houpis, 1966] |
D’Azzo, J.J. and Houpis, C.H. Feedback Control System Analysis and Synthesis, |
|
2nd edn., McGraw-Hill. New York, 1966. |
[Dorf, 1992] |
Dorf, R.C. Modern Control Systems, 6th edn., Addison-Wesley, Reading, MA, 1992. |
[Evans, 1948] |
Evans, W.R. Graphical analysis of control systems. Trans. AIEE, 67, 1948, pp. 547– |
|
551. |
[Evans, 1954] |
Evans, W.R. Control System Dynamics. McGraw-Hill. New York, 1954. |
[Houpis and Lamont, 1962] |
Houpis, C.H. and Lamont, G.B. Digital Control Systems: Theory, Hardware, Soft- |
|
ware, 2nd edn., McGraw-Hill, New York, 1962. |
[Jury, 1958] |
Jury, E.I. Sampled-Data Control Systems. John Wiley & Sons, Inc., New York, 1958. |
[Katz, 1981] |
Katz, P. Digital Control Using Microprocessors. Prentice Hall, Englewood Cliffs, |
|
NJ, 1981. |
[Kuo, 1962] |
Kuo, B.C. Automatic Control Systems. Prentice Hall, Englewood Cliffs, NJ, 1962. |
[Kuo, 1963] |
Kuo, B.C. Analysis and Synthesis of Sampled-Data Control Systems. Prentice Hall, |
|
Englewood Cliffs, NJ, 1963. |
[Lindorff, 1965] |
Lindorff, D.P. Theory of Sampled-Data Control Systems. John Wiley & Sons, Inc., |
|
New York, 1965. |
[Phillips and Harbor, 1988] |
Phillips, C.L. and Harbor, R.D. Feedback Control Systems. Englewood Cliffs, NJ, |
|
Prentice Hall, 1988. |
[Soliman and Srinath, 1990] Soliman, S.S. and Srinath, M.D. Continuous and Discrete Signals and Systems. |
|
|
Prentice Hall, Englewood Cliffs, NJ, 1990. |
[Strum and Kirk, 1988] |
Strum, R.D. and Kirk, D.E. First Principles of Discrete Systems and Digital Signal |
|
Processing. Addison-Wesley, Reading, MA, 1988. |
DIGITAL CONTROLLERS |
215 |
Equation (9.3) states that the required controller D(z) can be designed if we know the model of the process. The controller D(z) must be chosen so that it is stable and can be realized. One of the restrictions affecting realizability is that D(z) must not have a numerator whose order exceeds that of the numerator. Some common controllers based on (9.3) are described below.
9.1.1 Dead-Beat Controller
The dead-beat controller is one in which a step input is followed by the system but delayed by one or more sampling periods, i.e. the system response is required to be equal to unity at every sampling instant after the application of a unit step input.
The required closed-loop transfer function is then
T (z) |
= |
z−k , where k |
≥ |
1. |
(9.4) |
|
|
|
|
From (9.3), the required digital controller transfer function is
D(z) = HG(z) 1 |
|
T (z) |
= HG(z) |
1 |
z−z |
k . |
(9.5) |
|||
1 |
|
|
T (z) |
1 |
|
|
k |
|
||
|
|
|
− |
|
|
|
|
− − |
|
|
An example design of a controller using the dead-beat algorithm is given below.
Example 9.1
The open-loop transfer function of a plant is given by
e−2s
G(s) = .
1 + 10s
Design a dead-beat digital controller for the system. Assume that T = 1 s.
Solution
The transfer function of the system with a zero-order hold is given by
HG(z) = Z 1 −s − |
|
|
|
G(s) = (1 − z−1)Z s(1e− |
10s) |
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
e |
s T |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2s |
|
|
|
|
|
|
|
||||
or |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
10s) = (1 − z−1)z−2 Z |
s(s 1/ |
1/10) |
. |
||||||||||||||||||||||||||||
HG(z) = (1 − z−1)z−2 Z s(1 |
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
|
|||
From z-transform tables we obtain |
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
HG(z) (1 |
|
|
z−1)z−2 |
|
|
|
|
z(1 − e−0.1) |
|
|
|
z−3 |
(1 − e−0.1) |
|
|
|||||||||||||||||||||||||||
|
= |
|
− |
|
|
|
(z |
− |
1)(z |
− |
e |
− |
0.1) = |
|
1 |
− |
e |
− |
0.1z |
− |
1 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
or |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
HG(z) = |
|
|
0.095z−3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
1 − 0.904z−1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
From Equations (9.3) and (9.5), |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
D(z) |
|
|
|
1 |
|
T (z) |
|
|
|
|
|
1 − 0.904z−1 |
|
|
|
|
z−k |
|
|
. |
|
|
|
|
|
|
||||||||||||||||
= HG(z) 1 |
− |
|
T (z) = |
|
|
|
|
− |
|
|
− |
k |
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
0.095z |
− |
3 |
1 |
z |
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
