Acessibilidade / Reportar erro

Simple answers to usual questions about unusual forms of the Evans' root locus plot

Abstracts

After a first contact with Evans' root locus plots, in an introductory course about classical control theory, students usually pose questions for which the answers are not trivially found in the usual textbooks. Examples of such questions are: Can a branch intersect itself? Can two or more branches be coincident? Can a branch intersect its asymptote? In this paper devoted to helping teaching, numerical examples and an incremental property are used for answering some questions about unusual forms of root loci related to closed-loop control systems. In some cases, answering to these questions can be fundamental for making a correct sketch of a root locus.

feedback control; root locus method


Após um primeiro contato com gráficos de Evans de lugar das raízes, num curso introdutório sobre teoria de controle clássico, os estudantes normalmente fazem perguntas para as quais as respostas não são encontradas trivialmente nos livros-textos usuais. Exemplos de tais questões são: Um ramo pode se auto-interseccionar? Dois ou mais ramos podem ser coincidentes? Um ramo pode interseccionar sua assíntota? Neste artigo com fins didáticos, usam-se exemplos numéricos e uma propriedade incremental para responder algumas questões sobre formas incomuns de lugar das raízes relacionadas a sistemas de controle em malha fechada. Em alguns casos, responder a essas questões pode ser fundamental para fazer um esboç o correto de um lugar das raízes.

controle com realimentação; método do lugar das raízes


SISTEMAS DE CONTROLE

Simple answers to usual questions about unusual forms of the Evans' root locus plot

L.H.A. MonteiroI; J.J. Da CruzII

IEscola de Engenharia da Universidade Presbiteriana Mackenzie e Escola Politécnica da Universidade de São Paulo - luizm@mackenzie.br, luizm@usp.br

IIEscola Politécnica da Universidade de São Paulo - jaime@lac.usp.br

ABSTRACT

After a first contact with Evans' root locus plots, in an introductory course about classical control theory, students usually pose questions for which the answers are not trivially found in the usual textbooks. Examples of such questions are: Can a branch intersect itself? Can two or more branches be coincident? Can a branch intersect its asymptote? In this paper devoted to helping teaching, numerical examples and an incremental property are used for answering some questions about unusual forms of root loci related to closed-loop control systems. In some cases, answering to these questions can be fundamental for making a correct sketch of a root locus.

Keywords: feedback control, root locus method.

RESUMO

Após um primeiro contato com gráficos de Evans de lugar das raízes, num curso introdutório sobre teoria de controle clássico, os estudantes normalmente fazem perguntas para as quais as respostas não são encontradas trivialmente nos livros-textos usuais. Exemplos de tais questões são: Um ramo pode se auto-interseccionar? Dois ou mais ramos podem ser coincidentes? Um ramo pode interseccionar sua assíntota? Neste artigo com fins didáticos, usam-se exemplos numéricos e uma propriedade incremental para responder algumas questões sobre formas incomuns de lugar das raízes relacionadas a sistemas de controle em malha fechada. Em alguns casos, responder a essas questões pode ser fundamental para fazer um esboç o correto de um lugar das raízes.

Palavras-chave: controle com realimentação, método do lugar das raízes.

1 INTRODUCTION

Analysis and design of linear time-invariant, single-input/single-output, feedback control systems can be performed by using the root locus method. For instance, this method can be employed for determining the lock-in range of phase-locked loops (e.g. Piqueira and Monteiro, 2006). Root locus method was introduced in 1948 by W.R. Evans (1948, 1950) in order to make possible sketching the path (in the complex plane) of the closed-loop poles of control systems, as the value of a single parameter is varied. This parameter is normally the open-loop gain k > 0. If G(s) is the open-loop transfer function, the closed-loop poles are the roots of the characteristic equation (e.g. Ogata, 2001):

By taking G(s) ≡ b(s)/a(s), where a(s) and b(s) are real coprime polynomials of degrees n and m, respectively, then (1) can be written as:

Evans derived rules that allow to sketch the loci (the branches) of the roots of equation (2) as k is varied from zero to infinity, without solving such an equation (the case k < 0 is studied, for instance, by Eydgahi and Ghavamzadeh, 2001; Teixeira et alii, 2004). This sketch is based on the values of the poles pj of G(s) (j = 1,...,n), which are the n roots of the polynomial a(s); and on the zeros zj of G(s) (j = 1,...,m), which are the m roots of the polynomial b(s). For physically realizable systems, then n > m.

Evans' rules are usually presented in an introductory subject on control systems in undergraduate engineering courses and they are simple to be applied. Some well-known rules are (e.g. Kuo, 1967; Franklin et alii, 1986; Ogata, 2001): the number of branches of the root locus is equal to n; if n > m, then for k → ∞ there are v n – m branches approaching straight-line asymptotes; branches can leave from or arrive to the real axis, and so on.

However, after seeing some examples of root locus plot, students usually ask questions that can not be easily answered with the rules presented in the course and found in the textbooks. For instance: Which are the angles of departure (or arrival) of the branches leaving from (or arriving to) the real axis? Can a branch intersect itself? Can two or more branches have a common point? Can two or more branches be coincident for an open interval of values of k? Can a branch intersect its asymptote? Here, we answer such questions with numerical examples and/or arguments based on analytical expressions. These answers can be found in the literature (as cited below); however, they are barely used for sketching root loci. Here we present examples where the sketches of the corresponding root loci can not be appropriately made without such answers.

In the next section, an incremental property of the root locus plot is derived from a didactical block diagram. This property is employed in part of this study. The presented root loci were numerically plotted by using the software Matlab.

2 INCREMENTAL PROPERTY

Consider the feedback control system illustrated in Fig. 1 where 0 < q,k < ∞. The closed-loop poles are given by:

or:


When q = 0, the root locus plot of this system is sketched by the usual rules, because the characteristic equations (2) and (4) become the same. Now, suppose that q = q1 > 0. In this case, the root locus plot can be determined from the zeros of G(s) (the roots of b(s)) and the closed-loop poles of G(s) for q = q1 (the roots of a(s) + q1b(s)). This result represents an incremental property of the root locus plot: the closed-loop poles for k > k1q1 + k0 is a function of the open-loop zeros and of the closed-loop poles for k = k1. The usual situation occurs for k1 = 0, which implies that the roots of equation (2) depend on the roots of the polynomials a(s) and b(s), as Evans proved.

Incremental property, also denominated continuation property (e.g. Kuo, 1967; Pan and Chao, 1978; Franklin et alii, 1986; Paor, 2000), is an obvious consequence of the diagram shown in Fig. 1, which was not found in any control texts: the inner loop can be considered as the open-loop to draw the root locus for a varying k > 0.

3 ARRIVAL AND DEPARTURE ANGLES

The angle associated to the complex number s satisfying the equation (2) must be an odd multiple of 180º, because kG(s) = -1. Thus:

where l = 0, ± 1, ± 2, .... The angles of the vectors from the zeros zj and the poles pj to the number s are, respectively, øj≡ arg(s - zj) and qj ≡ arg(s - pj).

A breakaway point is a real root s = s1 of equation (2) occurring at k = k* where r branches (with r > 2) depart from the real axis at s = s1 for k > k*. A break-in point is a real root s = s1 of equation (2) occurring at k = k* where r branches (with r > 2) arrive to the real axis at s = s1 for k > k*. In both cases, when k = k*, then s1 is a root of equation (2) with multiplicity r.

The incremental property can be used for obtaining angles of arrival and angles of departure from break-in and breakaway points, respectively. The usual example is given by G(s) = 1/(s2 + 2s), where the closed-loop poles are the roots of s2 + 2s + k = 0. Notice that s = s1 = -1 is a root with multiplicity two for k = k* = 1; and for k > 1, both roots become complex numbers. The departure angles from the real axis at s = s1 can be determined from the closed loop poles at k = k*. The expression (5) reveals that these angles are:

with j = 0,1,..., r - 1. In this case, r = 2 and the departure angles are 90º and 270º, as shown in Fig. 2.


For G(s) = 1/(s + 1)3, the departure angles are 60º, 180º and 300º, as illustrated in Fig. 3. These angles are obtained from expression (6) with r = 3. Observe that s = s1 = -1 is the breakaway point at k = k* = 0, which is also the open-loop pole with multiplicity three.


For G(s) = 1/(s4 + 8s3 + 18s2 + 16s + 1), then s = s1 = -1 is a root with multiplicity three for k = k* = 4, as shown in Fig. 4. The departure angles obtained from (6) are the same of the previous example, which are 60º, 180º and 300º.


Similar plots were already presented by, for instance, Kuo (1967) and Krajewski and Viaro (2007).

4 INTERSECTION AMONG BRANCHES

Breakaway and break-in points are usual examples of common points among branches on the real axis. However, there can be such points also out of the real axis.

Two or more branches can intersect at a point that does not belong to the real axis.

Here intersection means having a common point (a different meaning is given by Lundberg (2007)). This situation occurs whenever there are multiple closed-loop poles outside the real axis. For instance, G(s) = 1/(s4 + 4s3 + 8s2 + 8s + 3) presents s = -1 ± i as closed-loop poles for k = 1, both poles with multiplicity two. In this example, the departure angles from these complex poles are 0º and 180º, as shown in Fig. 5 (found also in Kuo (1967)), because expression (6) reduces to:

with j = 0,1.


A branch can not intersect itself.

This property is immediate since, by construction, the gain k increases along the root locus branches. By contradiction, if such an intersection occurred there would be a point corresponding to two distinct values of k.

Two or more branches can not be coincident in an open interval of values of k.

This is because if branches intersect at a point of the complex plane for k = k', then they must separate following distinct angles for k = k' + q, with q → 0+ (these angles can be calculated by expression (6)). Consequently, these branches are no more coincident for k = k' + q.

5 INTERSECTION BETWEEN A BRANCH AND ITS ASYMPTOTE

In Fig. 2, both branches for k > 1 are coincident with the corresponding asymptotes. A similar behavior is observed in Fig. 3, where the three branches are coincident with their asymptotes for k > 0. Thus, it is easy to find examples where a branch and its asymptote coincide for an open interval of values of k.

However, it is hard to build a transfer function where this coincidence happens at a unique point. In fact, such a behavior is possible and it is illustrated in Fig. 6, obtained for G(s) = (s2 + 2s + 10)/(s5 + 9s4 + 40.25s3 + 106.5s2 + 141.75s + 67.5) (other examples can be found in Kuo (1967) and Ogata (2001)). Notice that two branches present a common point with the corresponding straight-line asymptotes (represented in the figure by dotted lines forming angles of 60º and 300º with respect to the real axis). The third asymptote (at 180º) is coincident with its corresponding real branch for k > 0.


6 HARD TO SKETCH

There are root loci hard of correctly sketching without investigating the existence of common points among branches out the real axis and without determining the arrival and departure angles related to such points. In some cases, intersections among branches and asymptotes should also be investigated in order to perform a right sketch. For instance, take G(s) = (s2 + (5/6)s +(1/3))/(s3(s2 + (6 + α)s + 16 + β). Fig. 7 corresponds to α = 0 and β = 0, Fig. 8 to α = 0 and β = 1, and Fig. 9 to α = 1 and β = 0. In all cases, the departure angles from the origin (a triple pole) are 60º, 180º and 300º, as in Fig. 4. In Fig. 7 there are two common points out of the real axis (in -1 ± i); in the other cases such points do not exist. In Figs. 7 and 8 there are neither breakaway point nor break-in point and the branches going to infinity in the right half-plane arrive from the complex poles. However, in Fig. 7 two branches intersect twice their asymptotes; in Fig. 8 such intersections do not occur. In Fig. 9 there are one breakaway point and one break-in point, the branches going to infinity in the right half-plane come from the origin, and intersections among branches and asymptotes do not happen. Observe that correct sketches of the these cases can be made only if the common points among branches out the real axis and the intersections among branches and asymptotes are calculated. For instance, as in Fig. 8 there is no intersection among branches and asymptotes, then the branches going to infinity in the right half-plane must arrive from the complex poles and the ones going to zeros must arrive from the origin. In Fig. 9 as such intersections do not occur either, then the branches going to infinity in the right half-plane must come from the origin (because, if not, then two branches coming from the origin should intersect the asymptotes in order to arrive in the break-in point situated at the left of the point where the asymptotes meet the real axis).




By comparing the three last figures, it is evident that a root locus plot can be drastically changed by slightly varying G(s) (in these cases, by slightly varying the positions of the open-loop poles and keeping unaltered the positions of the zeros).

As a "homework" we suggest to sketch the case with α = 0 and β = -1.

7 FINAL COMMENTS

By using an incremental property (here didactically explained) and/or numerical examples, we answered the questions posed at the Introduction about unusual forms of root loci. The investigated features can be really useful for sketching root loci, as shown in the last section. We hope that the examples as well as the reasoning presented can help root locus teaching.

The authors are partially supported by CNPq.

8 REFERENCES

Artigo submetido em 31/01/2008

1a. Revisão em 26/05/2008

2a. Revisão em 13/10/2008

Aceito sob recomendação do Editor Associado Prof. José Roberto Castilho Piqueira

  • Evans, W.R. (1948). Analysis of Control System, AIEE Transactions, vol. 67, pp. 547-551.
  • Evans, W.R. (1950). Control System Synthesis and Root Locus Method, AIEE Transactions, vol. 69, pp. 66-69.
  • Eydgahi, A.M., Ghavamzadeh, M. (2001). Complementary Root Locus Revisited, IEEE Transactions on Education, vol. 44, pp. 137-143.
  • Franklin, G.F., Powell, J.D., Emani-Naeini, A. (1986). Feedback Control of Dynamic Systems Addison-Wesley, Reading, USA.
  • Krajewski, W., Viaro, U. (2007). Root-Locus Invariance - Exploiting Alternative Arrival and Departure Points, IEEE Control Systems Magazine, vol. 27(1), pp. 36-43.
  • Kuo, B.C. (1967). Automatic Control Systems, Prentice-Hall, Englewood-Cliffs, USA.
  • Lundberg, K.H. (2007). Classical Crossroads, IEEE Control Systems Magazine, vol. 27(6), pp. 20-21.
  • Ogata, K. (2001). Modern Control Engineering, Prentice-Hall, New York, USA.
  • Pan, C.T., Chao, K.S. (1978). A Computer-Aided Root Locus Method, IEEE Transactions on Automatic Control, vol. 23, pp. 856-860.
  • Paor, A.M. (2000). The Root Locus Method: Famous Curves, Control Designs and Non-control Applications, International Journal of Electrical Engineering Education, vol. 37, pp. 344-356.
  • Piqueira, J.R.C., Monteiro, L.H.A. (2006). All-Pole Phase-Locked Loops: Calculating Lock-In Range by Using Evan's Root-Locus, International Journal of Control, vol. 79, pp. 822-829.
  • Teixeira, M.C.M., Assunção, E., Machado, E.R.M.D. (2004). A Method for Plotting the Complementary Root Locus Using the Root-Locus (Positive Gain) Rules, IEEE Transactions on Education, vol. 47, pp. 405-409.

Publication Dates

  • Publication in this collection
    20 Jan 2009
  • Date of issue
    Dec 2008

History

  • Reviewed
    13 Oct 2008
  • Received
    31 Jan 2008
Sociedade Brasileira de Automática Secretaria da SBA, FEEC - Unicamp, BLOCO B - LE51, Av. Albert Einstein, 400, Cidade Universitária Zeferino Vaz, Distrito de Barão Geraldo, 13083-852 - Campinas - SP - Brasil, Tel.: (55 19) 3521 3824, Fax: (55 19) 3521 3866 - Campinas - SP - Brazil
E-mail: revista_sba@fee.unicamp.br