Projective Root-Locus: An Extension of Root-Locus Plot to the Projective Plane Francisco Mota Departamento de Engenharia de Computação e Automação Universidade Federal do Rio Grande do Norte Brasil e-mail:mota@dca.ufrn.br September 10, 2018 arxiv:1409.4476v2 [cs.sy] 18 Sep 2014 Abstract In this paper we present an extension of the classical Root-Locus (RL) method where the points are calculated in the real projective plane instead of the conventional affine real plane; we denominate this extension of the Root-Locus as Projective Root-Locus (PjRL). To plot the PjRL we use the concept of Gnomonic Projection in order to have a representation of the projective real plane as a semi-sphere of radius one in R 3. We will see that the PjRL reduces to the RL in the affine XY plane, but also we can plot the RL onto another affine component of the projective plane, like ZY affine plane for instance, to obtain what we denominate complementary plots of the conventional RL. We also show that with the PjRL the points at infinity of the RL can be computed as solutions of a set algebraic equations. Index terms Root-Locus, Projective Plane, Gnomonic Projection, Algebraic Geometry, Affine Algebraic Variety, Projective Algebraic Geometry, Ideal of Polynomials, Grobner Basis. 1 Introduction The Root-Locus (RL) method is a classical tool that has been used extensively in the feedback control literature for studying the stability and performance of a closed loop linear feedback system. It consists of a parametric plot of the roots of the polynomial p(s) = d(s) + kn(s) in in the complex plane, as the parameter k spans R; d and n are fixed coprime polynomials, and d is monic with degree, in general, greater than the degree of n. In fact, the polynomial p represents the denominator of the transfer function of a closed-loop feedback system that has the (irreducible) proper rational function G(s) = n(s)/d(s) as a linear time invariant plant model and k as a (proportional type) controller (see Figure 1) and that is why we use the terminology the RL for G(s). To plot the RL for a given G(s), most control theory textbooks presents a set of rules that allow us to make an approximate sketch of the plot ([1]), but a detailed plot, nowadays, in general, is obtained using a computer software that evaluates the roots of p, using numerical techniques, for a given range of the parameter k in R (e.g. Scilab ([2])). In Figure 2 we show the plot of the RL for the plant G(s) = (s + 1)/s 2, for some range of k R. The motivation to use the projective plane to analyze the RL method is that the RL plot for a given G(s) can have points at infinity: the parameter k itself has to reach an infinite value in order we can obtain the terminal points of the RL, that can, in turn, be finite (zeros of G(s)) or to be located at infinity. In this way, using the concepts of projective real line and projective plane we can account for these infinite points, and also obtain complementary plots of the RL where points at infinity can be plotted at a finite position onto an affine plane. We denominate this extension of the RL to the projective plane as Projective Root-locus (PjRL) and, in spite of its abstract definition, we will show that it can be relatively easy to Figure 1: Control Feedback Loop with a Proportional Controller 1
Figure 2: Root-Locus for G(s) = (s + 1)/s 2 obtain the PjRL for G(s) using a computer algebra software. Below we introduce the definitions and notation to be used along the paper: R, C and R[x 1, x 2,..., x n ]: Represents the field of real numbers, the field of complex numbers and the ring of polynomials with coefficient s in R and with indeterminates (x 1, x 2,..., x n ), respectively. Projective (real) line: The projective line P 1 (R) is the set of slopes y/x, (x, y) (0, 0) R 2 and 1/0 =. So, if k P 1 (R), then k = k n /k d, and k = corresponds to k n = 1 and k d = 0 (k n = k d = 0 is not allowed). We note that P 1 (R) = R { }. Projective (real) plane: The projective plane P 2 (R) is the set of equivalence classes of all nonzero triples (x, y, z) R 3 under the equivalence relation: (α 1, α 2, α 3 ) (β 1, β 2, β 3 ) if α i = λβ i, for some λ 0. We represent the equivalence class of (x, y, z) by (x : y : z), that is denominated homogeneous coordinate of (x, y, z). We note that P 2 (R) = R 2 H, where H represents the plane at infinity and it is disjoint from R 2. Mathematically, R 2 = {(x : y : 1) P 2 (R)}, the XY plane, and H = {(x : y : 0) P 2 (R)}. In fact, H has two types of points, namely, (1 : m : 0) and (0 : 1 : 0), where (1 : m : 0) represents the point of intersection of all (XY ) lines with finite slope m and (0 : 1 : 0) represents the intersection of all lines with infinite slope (vertical lines). Homogeneous polynomial A polynomial (in several variables) is homogeneous when all of its nonzero terms (monomials) have the same total degree. One important fact about a homogeneous polynomial p R[x 1, x 2,..., x n ] is that p(λx 1, λx 2,..., λx n ) = λ d p(x 1, x 2,..., x n ), where d is the total degree of p. We always can turn a non-homogeneous polynomial (q) into a homogeneous one (q h ) by adding a new variable (x n+1 ), with the following procedure: q h (x 1,..., x n, x n+1 ) = x d n+1 q(x 1 /x n+1, x 2 /x n+1,..., x n /x n+1 ), where d is the total degree of q; this process is denominated homogenization of q. We can always de-homogenize q h by setting x n+1 = 1 and recover back q. Affine Algebraic Variety An affine (real) algebraic variety V generated by a set of m polynomials, p i R[x 1, x 2,..., x n ], is a subset of the affine plane R n composed by the coordinates (x 1, x 2,..., x n ) that are simultaneous real roots of the m generating polynomials, that is p i (x 1, x 2,..., x n ) = 0, i = 1,..., m. Projective Algebraic Variety A projective (real) algebraic variety W generated by a set of m homogeneous polynomials, p i R[x 1, x 2,..., x n, x n+1 ], is a subset of the projective plane P n (R) composed by the homogeneous coordinates (x 1 : x 2 :... : x n : x n+1 ) such that (x 1, x 2,..., x n, x n+1 ) is a simultaneous real root of the m generating homogeneous polynomials, that is p i (x 1, x 2,..., x n, x n+1 ) = 0, i = 1,..., m. We note that if (x 1,..., x n, x n+1) is any member of the equivalence class (x 1 : x 2 :... : x n : x n+1 ), it is also a simultaneous root of p i, i = 1,..., m, since p i is homogeneous of degree d: x j = λx j, then p i (x 1,..., x n, x n+1) = λ d p i (x 1, x 2,..., x n, x n+1 ) = 0. Ideal of Polynomials: Let be {p 1, p 2,..., p t } a set of polynomials in R[x 1, x 2,..., x n ]. The set of polynomials I R[x 1, x 2,..., x n ] defined by I = t h i p i, h i R[x 1, x 2,..., x n ] i=1 2
is an ideal of R[x 1, x 2,..., x n ], and {p 1, p 2,..., p t } is denominated a generating set for I; in this case we write I = p 1, p 2,..., p t. A Grobner Basis for the ideal I is a particular kind of generating set that allows many important properties of the ideal to be deduced easily. Given a generating set {p 1, p 2,..., p t } for I, we can obtain a Grobner basis {g 1, g 2,..., g s } for I algorithmically (see [3, Ch. 2]). For more details about the concepts above see ([3], [4]). 2 The Projective Root-Locus - PjRL As discussed in Introduction, the conventional RL for an irreducible proper rational function G(s) = n(s)/d(s) is a plot of the roots of the polynomial p(s) = d(s) + kn(s), when k R; that is, we solve d(s) + kn(s) = 0 (1) for each k R and plot its roots in the affine plane R 2. But, since the parameter k belongs to R, to analyze the situation where k ±, we will modify Equation (1) slightly by considering k P 1 (R). So, following the definition of P 1 (R), we set k = k n /k d in (1) and clear the denominator to obtain: k d d(s) + k n n(s) = 0. (2) We note that setting k d = 1 in Equation (2) we recover Equation (1) and setting k d = 0, that is k =, corresponds to n(s) = 0 in Equation (2), or the finite terminal points of the RL (zeros of G(s)). We then see that the effect of passing from k R to k P 1 (R) is just that of including the roots of n(s), the finite terminal points, into the RL. As we will see in the next sections, the infinite terminal points of the RL will only appear when we extrapolate from R 2 to P 2 (R). We also note that we can treat the case where the degree of d is less than the degree of n in the same fashion we treat the case where the degree of of d is greater than the degree of n by just exchanging the positions of k d and k n in Equation 2. The case where the degree of d is equal the degree of n also can be treated by our approach, as shown in examples of Section 3. Since Equation (2) may admit complex solutions, if we write s = x + iy we have: d(x + iy) = q d (x, y) + ir d (x, y) and n(x + iy) = q n (x, y) + ir n (x, y) (3) where q d, r d, q n and r n are polynomials in R[x, y]. So we may rewrite (2) as: [k d q d (x, y) + k n q n (x, y)] + i [k d r d (x, y) + k n r n (x, y)] = 0 and finding a complex solution for (2), for given pair (k d, k n ), is equivalent of finding a solution in R 2 for the system of polynomial equations: k d q d (x, y) + k n q n (x, y) = 0 (4) k d r d (x, y) + k n r n (x, y) = 0 (5) for each k n /k d P 1 (R). It is important to stress the fact that any solution for the system (4 5) must be invariant when we pass from pair (k n, k d ) to (λk n, λk d ), λ 0, since they represent the same point in P 1 (R). This, in fact, is true because it is equivalent to multiply Equations (4) and (5) by λ 0. To obtain the Projective Root-Locus (PjRL) we need to extend the solutions of Equations (4 5), defined above, from the affine plane R 2, to the projective plane P 2 (R). To achieve this goal, we first need to interpret the solutions of Equation (4 5) as a real algebraic variety V generated by the set of two polynomials q and r defined as: q(x, y, k d, k n ) = k d q d (x, y) + k n q n (x, y) (6) r(x, y, k d, k n ) = k d r d (x, y) + k n r n (x, y). (7) Since the polynomials q and r are defined in R[x, y, k d, k n ] we would have a variety in R 4, i.e. V R 4 ; but, since k n /k d is defined in P 1 (R), in fact, we have V R 2 P 1 (R). Based on this, we could abstractly interpret the RL as the projection (represented by V k ) of V onto R 2, since each point of the RL is a solution of (4 5) for a fixed k = k n /k d P 1 (R). We note that, for each k P 1 (R), V k is an (finite) affine real variety defined in R 2, by the solutions of Equations (4 5), or equivalently, by the roots of Equation (2). Now we proceed with the question of extrapolating the RL from the affine plane (R 2 ) to the projective plane (P 2 (R)). Our approach will follow the two steps bellow: (1) Extrapolate the algebraic variety V R 2 P 1 (R), defined above, to obtain a projective algebraic variety W P 2 (R) P 1 (R); 3
(2) Obtain the projection of W onto P 2 (R). This projection, represented by W k, k P 1 (R), it will be what we denominate PjRL. To obtain W from V, we could simply homogenize the polynomials q and r, presented in Equations (6 7), and obtain a projective variety W in P 2 (R) P 1 (R), now generated by the homogenized polynomials q h (x, y, z, k d, k n ) = z d q(x/z, y/z, k d /z, k n /z) and r h (x, y, z, k d, k n ) = z e r(x/z, y/z, k d /z, k n /z), as defined in Introduction. The projective variety obtained this way reduces to V in R 2 P 1 (R), since the process of de-homogenization of q h and r h will restore back the polynomials q and r. The flaw with this approach is that the process of simply homogenizing the generating polynomials for V, in general, creates a projective variety that is greater than the necessary, in the sense that it may add points at infinity to the original variety, other than the existing ones (see [3, Ch. 8]). Then, in fact, W must be the projective closure of V, that is, a minimal projective variety in P 2 (R) P 1 (R) that reduces to V in R 2 P 1 (R). To compute the closure of V, instead of directly homogenizing the polynomials q and r that generates V, we need first to compute a Grobner basis, with respect a graded monomial order, for the ideal I = q, r (see [3, Ch. 8]). The projective closure of V will be the projective variety W generated by the homogenized polynomials of the obtained Grobner basis. For the sake of completeness we present the following definition for the PjRL: Definition 2.1. (PjRL) Let be an irreducible rational function G(s) = n(s)/d(s) and consider the polynomials q and r as defined in Equations (6 7). We call the PjRL of G(s) the projection onto P 2 (R) of the projective algebraic variety W P 2 (R) P 1 (R), where W is generated by the set of homogenized polynomials of the Grobner basis {g 1, g 2,..., g s }, with respect a graded monomial order, for the ideal q, r. We will denote this projection by W k, where k P 1 (R). Remark 2.1. We will denote the set homogenized polynomials of the Grobner basis for q, r by {g1 h,..., gs h }, where gi h R[x, y, z, k d, k n ], and the variable z comes from the homogenization process, as defined in Introduction. Since we analyze W k in P 2 (R), we consider k = k n /k d P 1 (R) as a parameter, and the homogeneous polynomials gi h can been seen as defined in R[x, y, z]. Based on the fact that k P 1 (R) = R { }, we define: Initial points of the PjRL (W W 0 ): k = 0/1 = 0; that is, W 0 is generated by the polynomials {g1 h, g2 h,..., gs h } setting k n = 0 and k d = 1. W Terminal points of the PjRL (WW ): k = 1/0 = ; that is, W {g1 h, g2 h,..., gs h } setting k n = 1 and k d = 0. is generated by the polynomials W λ Intermediary points of the PjRL (WW λ ): k = λ/1, λ 0; that is W λ is generated by the polynomials {g1 h, g2 h,..., gs h } setting k n = λ R\{0} and k d = 1. We have the following comments regarding the results presented above: Calculating the Grobner basis for the ideal q, r is a relatively easy procedure using an algebra software available such as Macaulay2 ([6]), since we have only two polynomials that depends on four indeterminates, namely x, y, k d and k n. In the classical RL method there is a procedure for calculating the asymptotes based on the difference between the number of poles and zeros of G(s). In our case, the direction of these asymptotes will appear as the solution of the algebraic equations that defines W k and it will represent points at infinity. 2.1 PjRL plot in Projective Real Plane In order to plot the PjRL we can use the concept of gnomonic projection ([5]) to obtain a geometric representation the projective real plane. In this representation, P 2 (R) is identified with a semi-sphere of radius one in R 3, as shown in Figure 3. We note that the points P (on the plane) and P (on the sphere surface) in Figure 3 have the same homogeneous coordinates, since they belong to the same line in R 3. The points at infinity in P 2 (R) are identified with the equatorial great circle, remembering that antipodal points (opposite relative to the sphere center) have the same homogeneous coordinates. Also, we note that any line in the plane z = 1 corresponds to a great semicircle on the semi-sphere and the the left (right) z-semi-plane corresponds to the left (right) half of the semi-sphere. So, the PjRL plot is made onto this semi-sphere, and its gnomonic projection onto the plane z = 1 coincides with the conventional RL (see examples in Section 3). 4
Figure 3: Gnomonic Projection of half unit sphere onto plane z = 1 2.2 Complementary Root-Locus plot in ZY affine plane The equations for the PjRL will reduce to the equations for the RL when we set z = 1 in the set of homogeneous polynomials {g h 1,... g h s } that defines W. This means that when this set of polynomials is dehomogenized with respect the variable z we obtain the the RL, that is the intersection of the PjRL with the affine XY plane. But since the projective plane contains three sets that are copies of the affine planes XY, ZY and XZ, the PjRL also can give another view of the RL plot, when we analyze the intersection of it with the affine plane ZY, for instance. In this situation we de-homogenize the set of polynomials {g h 1,... g h s } with respect to the variable x, instead of z, and obtain a set of polynomials that defines a new affine variety in ZY plane that we will denominate it complementary RL. So, the complementary RL can been as a gnomonic projection onto the plane x = 1, instead of onto the plane z = 1 as shown in Figure 3. Geometrically, the switch of the role of variables x and z in the complementary RL have the effect of moving all the points over the line x = 0 (in XY plane) to the infinite and bringing the points at infinity (z = 0) to a finite position. Intuitively we could state that the conventional RL is a plot as seen from the beginning (k = 0) while the complementary RL is a plot as seen from the end (k = ). Also there exists an interesting relation between Y crossing points in the conventional RL and asymptotes in complementary RL. More specifically, if the RL crosses the Y axis at a point, say, (0 : y : 1) in XY plane for a given value of k, when we translate this point to the ZY plane it will become (1 : y : 0), that is a point a infinity, in fact an asymptote with slope y/1 = y in ZY plane. So, we conclude that the Y axis crossing points by the RL will become asymptotes in complementary RL, and the absolute value of the variable z will explode to infinity for the corresponding value of k. We also can make a similar analysis, now considering the intersection of the PjRL with the affine plane XZ. In examples presented in Section 3 we explore the concept of complementary RL with concrete computations. 3 Examples As a matter of fixing ideas, we present a series of examples below. Example 3.1. Let be G(s) = s/(s 2 + 1). In this case, using notation introduced in Equation (3), we easily obtain: q d = x 2 y 2 + 1, r d = 2xy, q n = x r n = y and using the definition of q and r in (6 7), we have: q(x, y, k d, k n ) = k d (x 2 y 2 + 1) + k n x, r(x, y, k d, k n ) = k d (2xy) + k n y (8) 5
Now we compute the Grobner basis for the ideal q, r using the graded reversed lexicographic order ([3, pp. 56]), with x > y > k d > k n. We used the software Macaulay2 ([6]) to make the computations and obtained the Grobner basis {g 1, g 2, g 3, g 4 }, where: g 1 (x, y, k d, k n ) = 2xyk d + yk n (= r) (9) g 2 (x, y, k d, k n ) = x 2 k d y 2 k d + xk n + k d (= q) (10) g 3 (x, y, k d, k n ) = x 2 yk n + y 3 k n yk n (11) g 4 (x, y, k d, k n ) = 2y 3 k d xyk n 2yk d (12) and the homogenized polynomials gi h Introduction: 1 of the Grobner basis are obtained using the procedure indicated in the g h 1 = z 3 g 1 (x/z, y/z, k d /z, k n /z) = 2xyk d + yzk n (13) g h 2 = z 3 g 2 (x/z, y/z, k d /z, k n /z) = x 2 k d y 2 k d + xzk n + z 2 k d (14) g h 3 = z 4 g 3 (x/z, y/z, k d /z, k n /z) = x 2 yk n + y 3 k n yz 2 k n (15) g h 4 = z 4 g 4 (x/z, y/z, k d /z, k n /z) = 2y 3 k d xyzk n 2yz 2 k d (16) The PjRL is the projection onto P 2 (R) of the projective algebraic variety W defined by the four polynomials gi h presented in Eqs. (13 16) above. We will represent this projection by W k, where k P 1 (R). Initial points of the PjRL (W 0 ): Setting k n = 0 and k d = 1 in Eqs. (13 16), we obtain: g1 h = 2xy, g2 h = x 2 y 2 + z 2, g3 h = 0 and g4 h = 2y(y 2 z 2 ), and we see that g1 h = 0 implies x = 0 or y = 0. If x = 0, by g2 h = 0 we obtain y 2 = z 2. In this case we cannot have y = 0 or z = 0, since (0, 0, 0) is not valid as a solution. So we can set z = 1 (XY plane) and obtain y = ±1. Then the initial points onto plane z = 1 are W 0 = {(0 : 1 : 1), (0 : 1 : 1)} = {(0, 1), (0, 1)}. Or, onto the semi-sphere of radius one: {( W 0 = 0, 1/ 2, 1/ ) ( 2, 0, 1/ 2, 1/ )} 2. Terminal points of the PjRL (W ): Setting k n = 1 and k d = 0 in Eqs. (13 16), we obtain: g1 h = yz, g2 h = xz, g3 h = y(x 2 + y 2 z 2 ) and g4 h = xyz. To solve gi h = 0, we can simplify the set of equations calculating a Grobner basis with this set of polynomials. The new Grobner basis has three polynomials: g1 h = yz, g2 h = xz, g3 h = y(x 2 + y 2 ). Then we have two kinds of points: 1. Points at finite position (z 0): In this case, solving g h i = 0, we get y = 0 and x = 0. So, the homogeneous coordinate of the point is (0 : 0 : z) = (0 : 0 : 1) which corresponds to (0, 0) onto affine plane z = 1. 2. Points at the infinite plane H (z = 0): In this case we obtain: g h 1 = 0, g h 2 = 0, and g h 3 = y(x 2 +y 2 ). We must have y = 0 and the unique possible nonzero solution is (x, 0, 0), x 0 whose homogeneous coordinate is (x : 0 : 0) = (1 : 0 : 0)}, which corresponds to the intersection point of the horizontal lines in the plane XY or a pair of points (±1, 0, 0) onto the equatorial great circle over the half sphere of radius one. We then have: W = {(0 : 0 : 1), (1 : 0 : 0)}. Intermediary points of the PjRL (W λ ): In this case we set k d = 1 and k n = λ in the polynomials shown in Eqs. (13 16) and recalculate the Grobner basis for the resulting set of polynomials to obtain: g h 1 = y(2x + zλ), g h 2 = x 2 y 2 + z 2 + xzλ, g h 3 = y(x 2 + y 2 z 2 ) (17) We see that we must have z 0 in (17), since z = 0 will imply x = y = 0 which is not valid as a solution; so all intermediary points are at finite positions. 1 in fact, since the polynomials g i are already homogeneous relative to k d and k n, we can homogenize them relative only to x and y, and the resulting gi h will be the same. For example, g 1 could be homogenized as g1 h = z2 g 1 (x/z, y/z, k d, k n). 6
Figure 4: PjRL plot for G(s) = s/(s 2 + 1) Figure 5: Conventional and complementary RL for G(s) = s/(s 2 + 1) To plot the PjRL over the unit semi-sphere with radius one we need to add the equation x 2 +y 2 +z 2 = 1, with z 0 to the set of equations (17). The sketch of the PjRL plot is shown in Figure 4. To obtain the conventional RL plot, we set z = 1 in (17) and get g 1 = y(2x+λ), g 2 = x 2 y 2 +xλ+1, and g 3 = y(x 2 +y 2 1). We note that that if y 0, g 3 = 0 will require x 2 + y 2 = 1. The complete plot for g i = 0 is show in Figure 5. Now we will analyze the complementary RL (in plane ZY ): Switching the roles of the x and z axis in the projective plane, the coordinate (x : y : z) will become (z : y : x). Then, re-analyzing the initial and terminal points calculated above, we have: W 0 = {(1 : 1 : 0), (1, 1 : 0) and W = {(1 : 0 : 0), (0 : 0 : 1)} We note that the initial points now are at infinity, that is they are asymptotes with rates ±1. Regarding the terminal points, there is one at infinity, that is (1 : 0 : 0), or a horizontal asymptote; and other at origin (0 : 0 : 1). To obtain the intermediary points, we set x = 1 in the polynomials shown in (17) above and we easily see that, if y 0, we have the hyperbola z 2 y 2 = 1. The plot for both conventional and complementary RL are shown in Figure 5. Remark 3.2. We note, in this example, that if we directly homogenize the polynomials q and r, defined in Equation (8) (instead of the Grobner basis polynomials shown in Equations (9 12)), we are only left with two equations, namely: q h = x 2 k d y 2 k d + xzk n + z 2 k d and r h = 2xyk d + yzk n, as opposed to the four equations (13 16). So, if we evaluate the terminal points W using only q h and r h above, we would have, after setting k d = 0 and k n = 1, q h = xz and r h = yz. It is easy to see, by these two last equations, that (1 : 1 : 0) would be a possible terminal point at infinity, and this point, as we know, is 7
spurious. This confirms the point discussed above, that we need to evaluate the Grobner basis for the ideal q, r. Example 3.3. Let be G(s) = (s + 1)/s 2, so we have q d (x, y) = x 2 y 2, r d (x, y) = 2xy, q n (x, y) = x + 1, r n (x, y) = y and using the definition of q and r in (6 7), we have: q(x, y, k d, k n ) = k d (x 2 y 2 ) + k n (x + 1), r(x, y, k d, k n ) = 2k d xy + k n y. Now we compute the Grobner basis for the ideal q, r using the graded reversed lexicographic order with x > y > k d > k n and obtain {g 1, g 2, g 3, g 4 }, where: g 1 (x, y, k d, k n ) = 2xyk d + yk n (= r) g 2 (x, y, k d, k n ) = x 2 k d y 2 k d + xk n + k n (= q) g 3 (x, y, k d, k n ) = x 2 yk n + y 3 k n + 2xyk n g 4 (x, y, k d, k n ) = 2y 3 k d xyk n 2yk n Now we homogenize of the polynomials g i, using the procedure indicated in the Introduction: g h 1 = z 3 g 1 (x/z, y/z, k d /z, k n /z) = 2xyk d + yzk n (18) g h 2 = z 3 g 2 (x/z, y/z, k d /z, k n /z) = x 2 k d y 2 k d + xzk n + z 2 k n (19) g h 3 = z 4 g 3 (x/z, y/z, k d /z, k n /z) = x 2 yk n + y 3 k n + 2xyzk n (20) g h 4 = z 4 g 4 (x/z, y/z, k d /z, k n /z) = 2y 3 k d xyzk n 2yz 2 k n (21) The PjRL is the set of projective algebraic varieties W k generated by the four polynomials gi h presented in Eqs. (18 21) above, for each k = k n /k d P 1 (R). Initial points of the PjRL (W 0 ): Setting k n = 0 and k d = 1 in Eqs. (18 21), we obtain: g1 h = 2xy, g2 h = x 2 y 2, g3 h = 0 and g4 h = 2y 3, and we see that the simultaneous solution for gi h = 0 is (0, 0, z) where z R. The homogeneous coordinate for this point is (0 : 0 : z), but since (0 : 0 : 0) is undefined we need z 0 and the unique homogeneous coordinate possible is (0 : 0 : 1). Then we have: W 0 = {(0 : 0 : 1)}, which represents the point (0, 0) in the affine plane (XY ). Terminal points of the PjRL (W ): Setting k n = 1 and k d = 0 in Eqs. (18 21), we obtain: g1 h = yz, g2 h = xz + z 2, g3 h = x 2 y + y 3 + 2xyz and g4 h = xyz 2yz 2. To solve gi h = 0, we can simplify the set of equations by recalculating a new Grobner basis with this set of polynomials. The new Grobner basis has three polynomials: g1 h = yz, g2 h = xz + z 2, g3 h = x 2 y + y 3. Then we have two kinds of points: 1. Points at the affine plane XY (z = 1): In this case we get: g 1 = y, g 2 = x+1, and g 3 = y(x 2 +y 2 ). We see that the unique solution possible is ( 1, 0, 1) whose homogeneous coordinate is ( 1 : 0 : 1), and this (homogeneous) point corresponds to ( 1, 0) in the affine plane XY. 2. Points at the infinite plane H (z = 0): In this case we obtain: g 1 = 0, g 2 = 0, and g 3 = x 2 y +y 3 = y(x 2 + y 2 ). We must have y = 0 and the unique possible nonzero solution is (x, 0, 0), x 0 whose homogeneous coordinate is (x : 0 : 0) = (1 : 0 : 0)}, which corresponds to the intersection point of the horizontal lines in the plane XY. By the RL plot for this G(s) (Figure 6) we see that this is the direction of the asymptotes when k ±. We then have: W = {( 1 : 0 : 1), (1 : 0 : 0)}. Intermediary points of the PjRL (W λ ): Setting k n = λ 0 and k d = 1 in Eqs. (18 21) we obtain: g h 1 = 2xy + yzλ, g h 2 = x 2 y 2 + xzλ + z 2 λ, g h 3 = λ(x 2 y + y 3 + 2xyz), g h 4 = 2y 3 xyzλ 2yz 2 λ. Computing a new Grobner basis with this set of polynomials we obtain three polynomials: g h 1 = 2xy + yzλ, g h 2 = x 2 y 2 + xzλ + z 2 λ, g h 3 = x 2 y + y 3 + 2xyz (22) We easily see that z = 0 will imply x = y = 0 and, since we are interested in nonzero solutions, we must have z 0 (implying that doesn t exists intermediary points in the infinite plane H). To see the graph in the XY plane, we set z = 1 in these polynomials and obtain g 1 = y(2x + λ), g 2 = x 2 y 2 + λ(x + 1), and g 3 = y(x 2 + y 2 + 2x) = y[(x + 1) 2 + y 2 1]. It is easy to see that solving g i = 0 we obtain the same set of affine varieties V λ, λ R, that represents the RL for G(s) in the affine XY plane, as shown in Figure 6. 8
Figure 6: Conventional and complementary RL for G(s) = (s + 1)/s 2 We now analyze the question of obtaining the complementary RL, that corresponds to the intersection of the PjRL with ZY plane. We note that switching the roles of the x and z axis in the projective plane, the coordinate (x : y : z) will become (z : y : x). Then, re-analyzing the points calculated above, we have: Initial points: The projective point W 0 = {(0 : 0 : 1)} will become W 0 = {(1 : 0 : 0)} and this corresponds to a point at infinity in the affine plane ZY (intersection of all horizontal lines). Terminal points: The projective points W = {( 1 : 0 : 1), (1 : 0 : 0)} will become W = {(1 : 0 : 1) = ( 1 : 0 : 1), (0 : 0 : 1)}. This corresponds to the points ( 1, 0) and (0, 0) in the affine plane ZY. Intermediary Points: To obtain the intermediary points W λ (22) to obtain: in the affine plane ZY we set x = 1 in We only have two cases (y = 0 and y 0): g 1 (y, z) = 2y + yzλ = y(2 + zλ) g 2 (y, z) = 1 y 2 + zλ + z 2 λ g 3 (y, z) = y + y 3 + 2yz = y(1 + y 2 + 2z) 1. y = 0: In this case we are left only with g 2 (0, z) = 0 or λz 2 + λz + 1 = 0 or z 2 + z + 1/λ = 0, whose roots are z 1,2 = ( 1 ± 1 4/λ)/2. 2. y 0: In this case, from g 3 (y, z) = 0 we have y 2 +2z+1 = 0, which is the parabola z = y 2 /2 1/2. The plot for this set of equations in plane ZY is shown in Figure 6. Example 3.4. Let be G(s) = 1 s((s + 4) 2 + 4 2. In this case we have: ) q = k d (x 3 3xy 2 + 8x 2 8y 2 + 32x) + k n, r = k d ( y 3 + 3yx 2 + 16xy + 32y) and computing the Grobner basis using the graded reversed lexicographic order with x > y > z > k d > k n, we obtain the following set of homogenized Grobner polynomials: g h 1 = 3x 2 yk n y 3 k n + 16xyzk n + 32yz 2 k n (23) g h 2 = 3x 2 yk d y 3 k d + 16xyzk d + 32yz 2 k d (24) g h 3 = x 3 k d 3xy 2 k d + 8x 2 zk d 8y 2 zk d + 32xz 2 k d + k n z 3 (25) g h 4 = 24xy 3 k d + 64y 3 zk d 64xyz 2 k d + 256yz 3 k d 9yz 3 k n (26) g h 5 = 24y 5 k d 320y 3 z 2 k d + 1280xyz 3 k d 27xyz 3 k n + 4096yz 4 k d 72yz 4 k n (27) Initial points (W 0 ): In spite of knowing that the initial points are the roots of d(s) we calculate them here just as a matter of checking the theory. Setting k d = 1 and k n = 0 in the polynomials above we get: 9
g h 1 = 0 g h 2 = 3x 2 y y 3 + 16xyz + 32yz 2 g h 3 = x 3 3xy 2 + 8x 2 z 8y 2 z + 32xz 2 g h 4 = 24xy 3 + 64y 3 z 64xyz 2 + 256yz 3 g h 5 = 24y 5 320y 3 z 2 + 1280xyz 3 + 4096yz 4 We note that doesn t exists initial points at infinity, since if we choose z = 0, by g h 5 = 0 we necessarily have y = 0 and, by g h 3 = 0 we get x = 0, but (0 : 0 : 0) is not valid as solution. So, all the initial points are onto affine XY plane; then setting z = 1 in the homogeneous polynomials above we obtain the set of (non-homogeneous) polynomials: g 2 = 3x 2 y y 3 + 16xy + 32y = y(3x 2 y 2 + 16x + 32) g 3 = x 3 3xy 2 + 8x 2 8y 2 + 32x g 4 = 24xy 3 + 64y 3 64xy + 256y = y(24xy 2 + 64y 2 64x + 256) g 5 = 24y 5 320y 3 + 1280xy + 4096y = y(24y 4 320y 2 + 1280x + 4096) We have two cases, namely y = 0 or y 0; if we set y = 0 the only possible real value for x is obtained from g 3 = 0 and it is x = 0. So the first initial point is (0 : 0 : 1). Considering y 0 in the polynomials above we can rewrite them as: g 2 /y = 3x 2 y 2 + 16x + 32 g 3 = x 3 3xy 2 + 8x 2 8y 2 + 32x g 4 /y = 24xy 2 + 64y 2 64x + 256 g 5 /y = 24y 4 320y 2 + 1280x + 4096 and if we compute a Grobner basis for the set of polynomials above we obtain: h 1 = 3x 2 y 2 + 16x + 32 h 3 = 3xy 2 + 8y 2 8x + 32 h 4 = 3y 4 40y 2 + 160x + 512 To solve the set of equations h 1 = h 2 = h 3 = 0, we can eliminate y 2 from h 1 and h 2 to obtain the equation x 3 + 8x 2 + 24x + 32 = 0 whose only real solution is x = 4, and this implies y = ±4. Finally we have the following set of initial points: W 0 = {(0 : 0 : 1), ( 4 : 4 : 1), ( 4 : 4 : 1)} Terminal points (W ): Setting k d = 0 and k n = 1 in Eqs. (23 27) we obtain: 1. Points onto affine plane XY (z = 1): This will make g h 3 = 1, what means that doesn t exist terminal points onto affine plane XY. 2. Points at infinity (z = 0): We are left, from (23), only with 3x 2 y y 3 = y(3x 2 y 2 ) = 0. If y = 0, we will have the solution (x : 0 : 0) = (1 : 0 : 0); if y 0, we will have y 2 = 3x 2, that results in points (1 : ± 3 : 0). So, the terminal points are: W = {(1 : 0 : 0), (1 : 3 : 0), (1 : 3 : 0)} Intermediary points (W λ ): Setting k d = 1, k n = λ 0 in Eqs. (23 27), and computing a new Grobner basis we have the set below with only two polynomials: g h 1 = 3x 2 y y 3 + 16xyz + 32yz 2, and g h 2 = z 3 λ + x 3 3xy 2 + 8x 2 z 8y 2 z + 32xz 2 (28) In this example we will only analyze the complementary RL (in plane ZY ), and for that we will recalculate the initial and terminal points, just by switching the position of x and z in W 0 and W. So the new initial and terminal points are: W 0 = {(1 : 0 : 0), ( 1/4 : 1 : 1), ( 1/4 : 1 : 1)} 10
1 Figure 7: Conventional and complementary RL for G(s) = s((s + 4) 2 + 4) and W = {(0 : 0 : 1), (0 : 3 : 1), (0 : 3 : 1)}. To evaluate the intermediary points we set x = 1 in the polynomials shown in (28) to obtain: h 1 = 3y y 3 + 16yz + 32yz 2 = y(3 y 2 + 16z + 32z 2 ) h 2 = z 3 λ 8y 2 z 3y 2 + 32z 2 + 8z + 1 And clearly we have two cases to consider, namely y = 0 and y 0: y = 0 will imply h 1 = 0 and h 2 = z 3 λ + 32z 2 + 8z + 1. We easily see that the cubic polynomial h 2 always have one real and two complexes roots. The real root varies with λ as shown in Figure 7. y 0 will imply ( h 1 /y) = y 2 32z 16z 3, so we have the polynomials: h 1 /y = y 2 32z 2 16z 3 h 2 = z 3 λ 8y 2 z 3y 2 + 32z 2 + 8z + 1 Now if we compute again a Grobner basis for this set we will have: l 1 = y 2 32z 2 16z 3 l 2 = (λ 256)z 3 192z 2 64z 8 We note that l 1 = 0 represents the hyperbole shown in Figure 7, while l 2 = 0 determines as z depends on λ. We note that when λ 256, l 2 is a cubic polynomial with one real and two complex roots; moreover the real root explodes to infinity when λ 256. We also note that 256 is the value of the gain k where the RL crosses the Y axis (in plane XY ). Example 3.5. We now consider a case where the degree of d is equal the degree of n. The main point here is that the polynomial d(s) + kn(s) may decrease its degree for some (finite) value of k. To simplify the PjRL plot we will consider a simple rational function G defined as: G(s) = 1 s2 1 + s 2 In this example we easily see that (1 + s 2 ) + (1 s 2 ) = 2, so the degree of d(s) + kn(s) is zero for k = 1, and we have no finite roots. To analyze the PjRL we evaluate the polynomials q and r, defined in (6 7), which in this case are: q = k d (1 + x 2 y 2 ) + k n (1 x 2 + y 2 ), and r = k d (2xy) + k n ( 2xy) 11
Computing the Grobner basis using the graded reversed lexicographic order with x > y > z > k d > k n, we obtain the following set of homogenized Grobner polynomials: g h 1 = xyk n (29) g h 2 = xyk d (30) g h 3 = k d (z 2 + x 2 y 2 ) + k n (z 2 + y 2 x 2 ) (31) g h 4 = y(y 2 k d y 2 k n z 2 k d z 2 k n ) (32) Initial points (W 0 ): Using k d = 1 and k n = 0 in equations (29 32) above we obtain: g h 2 = xy, g h 3 = z 2 + x 2 y 2, and g h 4 = y(y 2 z 2 ) and the unique (non-null) possible solution is x = 0, y 0 and y 2 = z 2 ; so the homogeneous coordinates for W 0 is: W 0 = {(0 : 1 : 1), (0 : 1 : 1)} On the semi-sphere of radius one, we have: W 0 = {(0, 1/ 2, 1/ 2), (0, 1/ 2, 1/ 2)} Terminal points (W ): Using k d = 0 and k n = 1 in equations (29 32) above we obtain: g h 1 = xy, g h 3 = z 2 x 2 + y 2, and g h 4 = y(y 2 + z 2 ) and the unique (non-null) possible solution is x 0, y = 0 and x 2 = z 2 ; so the homogeneous coordinates for W is: W = {(1 : 0 : 1), ( 1 : 0 : 1)} On the semi-sphere of radius one, we have: W = {(1/ 2, 0, 1/ 2), ( 1/ 2, 0, 1/ 2)} Intermediary points (W λ ):Using k d = 1 and k n = λ in equations (29 32) and recalculating the Grobner basis we get: g h 1 = xy (33) g h 2 = (λ 1)x 2 + (1 λ)y 2 (1 + λ)z 2 (34) g h 3 = y[ (1 λ)y 2 + (λ + 1)z 2 )] (35) We have two cases to consider, namely, x = 0, y 0 and x 0, y = 0 (the case x = 0 and y = 0 will imply z = 0 which is not valid as a solution). Also, to plot the PjRL over the semi-sphere of radius one we have to consider the restriction x 2 + y 2 + z 2 = 1 with z 0. Then we have: 1. x = 0, y 0 and y 2 + z 2 = 1. Using equations (33 35) defined above, we necessarily have (1 λ)y 2 (1 + λ)z 2 = 0, and then we get: 1 + λ 1 λ y = ±, and z = 2 2 So, considering λ > 0, we must have 0 < λ 1 and for λ = 1 we have the point at infinity (0, 1, 0) = (0, 1, 0). 2. x 0, y = 0 and x 2 + z 2 = 1. Using equations (33 35) defined above, we necessarily have (λ 1)x 2 (1 + λ)z 2 = 0, and then we get: λ 1 λ 1 x = ± 2λ, and z = 2λ Again, considering λ > 0, we must have λ 1, and for λ = 1 we get the point at infinity (1, 0, 0) = ( 1, 0, 0). In Figure 8 we show the plot for the PjRL obtained in this example; as we can note, the PjRL has a discontinuity at infinity. 12
Figure 8: PjRL plot for G(s) = 1 s2 1 + s 2 4 Conclusions We have presented in this paper an extension of the classical Root-Locus (RL) method, denominated Projective Root-Locus (PjRL), where the coordinates of the points of the RL for an irreducible rational function G(s) = n(s)/d(s) are represented in the projective real plane P 2 (R) instead of the affine plane R 2. To obtain the PjRL we used results from algebraic geometry, representing the RL as an affine algebraic variety and extrapolating it to the projective plane. With this approach we could obtain the RL points at infinity as solutions of a set of algebraic equations. Also, we have shown how to plot the PjRL onto a semi-sphere of radius one that is a representation of the projective plane P 2 (R). Since the real projective plane contains three copies of real affine planes, we can plot the RL onto an affine real plane, other than the original XY one; we denominated this new plot complementary RL, and we have shown that the points where the RL crosses the Y in XY plane turns into asymptotes of the complementary RL in ZY affine plane, and vice-versa. Several examples were worked out in order to show that the PjRL can be relatively easily obtained using a computer algebra software. References [1] J. D Azzo and C. Houpis. Linear Control System Analysis and Design. Second Edition. MacGraw-Hill Kogakusha, Ltd., 1981. [2] Scilab Enterprises. Scilab: Free and Open Source Software for Numerical Computation. Orsay, France, 2012. Available at http://www.scilab.org. [3] D. Cox, J. Litlle and D. O Shea.Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra. Second Edition. Springer-Verlag New York Inc., 1997. [4] I. Shafarevich. Basic Algebraic Geometry 1: Varieties in Projective Space. Second Edition. Springer- Verlag Berlin Heildelberg, 1994. [5] H. Coxeter and S. Greitzer. Geometry Revisited. Mathematical Association of America, 1967. [6] D. Grayson and M. Stillman. Macaulay2, A Software System for Research in Algebraic Geometry. Available at http://www.math.uiuc.edu/macaulay2. 13