Szidarovszky, F, Bahill, A..Stability Analysis The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton CRC Press llc. 2000
Szidarovszky, F., Bahill, A.T. “Stability Analysis” The Electrical Engineering Handbook Ed. Richard C. Dorf Boca Raton: CRC Press LLC, 2000
2 Stability analysis 12.2 Using the State of the System to Determine Stability 12.3 Lyapunov Stability Theory Ferenc Szidarovszky 12.4 Stability of Time-Invariant Linear Systems Stability Analysis with State-Space Notation. The Transfer A. Terry Bahill 12.5 BIBO Stability University of Arizona 12.6 Physical Examples 12.1 Introduction In this chapter, which is based on Szidarovszky and Bahill [ 1992], we first discuss stability in general and then present four techniques for assessing the stability of a system:(1)Lyapunov functions, (2)finding the eigenvalues for state-space notation, (3)finding the location in the complex frequency plane of the poles of the closed loop transfer function, and(4) proving bounded outputs for all bounded inputs. Proving stability with Lyapunov functions is very general: it works for nonlinear and time-varying systems. It is also good for doing proofs. Proving the stability of a system with Lyapunov functions is difficult, however, and failure to find a Lyapunov function that proves a system is stable does not prove that the system is unstable. The next techniques we present, finding the eigenvalues or the poles of the transfer function, are sometimes difficult, because they require factoring high-degree polynomials. Many commercial software packages are now available for this task, however. We think most engineers would benefit by having one of these computer programs. Jamshidi et al. [1992] and advertisements in technical publications such as the IEEE Control Systems Magazine and IEEE Spectrum describe many appropriate software packages. The last technique we present, bounded-input, bounded-output stability, is also quite general. Let us begin our discussion of stability and instability of systems informally. In an unstable system the state an have large variations, and small inputs or small changes in the initial state may produce large variations in the output. A common example of an unstable system is illustrated by someone pointing the microphone of a public address(PA)system at a speaker; a loud high-pitched tone results. Often instabilities are caused by oo much gain, so to quiet the PA system, decrease the gain by pointing the microphone away from the speaker Discrete systems can also be unstable. A friend of ours once provided an example. She was sitting in a chair she went over and turned up the thermostat on the heater. The house warme She got hot, so she got up and turned down the thermostat. The house cooled off. She got cold and turned up the thermostat. This process continued until someone finally suggested that she put on a sweater(reducing the gain of her heat loss system). She did, and was much more comfortable. We modeled this as a discrete system, because she seemed to sample the environment and produce outputs at discrete intervals about 15 c 2000 by CRC Press LLC
© 2000 by CRC Press LLC 12 Stability Analysis 12.1 Introduction 12.2 Using the State of the System to Determine Stability 12.3 Lyapunov Stability Theory 12.4 Stability of Time-Invariant Linear Systems Stability Analysis with State-Space Notation • The Transfer Function Approach 12.5 BIBO Stability 12.6 Physical Examples 12.1 Introduction In this chapter, which is based on Szidarovszky and Bahill [1992], we first discuss stability in general and then present four techniques for assessing the stability of a system: (1) Lyapunov functions, (2) finding the eigenvalues for state-space notation, (3) finding the location in the complex frequency plane of the poles of the closedloop transfer function, and (4) proving bounded outputs for all bounded inputs. Proving stability with Lyapunov functions is very general: it works for nonlinear and time-varying systems. It is also good for doing proofs. Proving the stability of a system with Lyapunov functions is difficult, however, and failure to find a Lyapunov function that proves a system is stable does not prove that the system is unstable. The next techniques we present, finding the eigenvalues or the poles of the transfer function, are sometimes difficult, because they require factoring high-degree polynomials. Many commercial software packages are now available for this task, however. We think most engineers would benefit by having one of these computer programs. Jamshidi et al. [1992] and advertisements in technical publications such as the IEEE Control Systems Magazine and IEEE Spectrum describe many appropriate software packages. The last technique we present, bounded-input, bounded-output stability, is also quite general. Let us begin our discussion of stability and instability of systems informally. In an unstable system the state can have large variations, and small inputs or small changes in the initial state may produce large variations in the output. A common example of an unstable system is illustrated by someone pointing the microphone of a public address (PA) system at a speaker; a loud high-pitched tone results. Often instabilities are caused by too much gain, so to quiet the PA system, decrease the gain by pointing the microphone away from the speaker. Discrete systems can also be unstable. A friend of ours once provided an example. She was sitting in a chair reading and she got cold. So she went over and turned up the thermostat on the heater. The house warmed up. She got hot, so she got up and turned down the thermostat. The house cooled off. She got cold and turned up the thermostat. This process continued until someone finally suggested that she put on a sweater (reducing the gain of her heat loss system). She did, and was much more comfortable. We modeled this as a discrete system, because she seemed to sample the environment and produce outputs at discrete intervals about 15 minutes apart. Ferenc Szidarovszky University of Arizona A. Terry Bahill University of Arizona
12.2 Using the State of the System to Determine Stability The stability of a system is defined with respect to a given equilibrium point in state space. If the initial state x is selected at an equilibrium state x of the system, then the state will remain at x for all future time. w the initial state is selected close to an equilibrium state, the system might remain close to the equilibrium state or it might move away. In this section we introduce conditions that guarantee that whenever the system starts near an equilibrium state, it remains near it, perhaps even converging to the equilibrium state as time increases For simplicity, only time-invariant systems are considered in this section. Time-variant systems are discussed in Section 12.5 Continuous, time-invariant systems have the for x(t)=f(x(t)) and discrete, time-invariant systems are modeled by the difference equation x(t+1)=f(x(t) (12.2) Here we assume that f: X-R", where XsR" is the state space. We also assume that function f is continuous; furthermore, for arbitrary initial state xo E X, there is a unique solution of the corresponding initial value roblem x(to)=Xo, and the entire trajectory x(o) is in X. Assume furthermore that to denotes the initial time It is also known that a vector x E X is an equilibrium state of the continuous system, Eq (12.1), if and only if f(x)=0, and it is an equilibrium state of the discrete system, Eq (12.2), if and only if x=f(x). In this chapter the equilibrium of a system will always mean the equilibrium state, if it is not specified otherwise. In analyzing the dependence of the state trajectory x(t) on the selection of the initial state xo nearby the equilibrium, the ollowing stability types are considered. Definition 12.1 1. An equilibrium state x is stable if there is an eo >0 with the following property: For all E1,00 such that if llx-xolt 2. An equilibrium state x is asymptotically stable if it is stable and there is an e >0 such that whenever ‖x-x‖<e, then x(t)→xast→ 3. An equilibrium state x is globally asymptotically stable if it is stable and with arbitrary initial state xo ∈Xx(t)→xast→∞. Unstable The first definition says an equilibrium state x is stable if the entire trajectory x(t) is closer to the equilibrium state than any small E, if the initial state xo is selected close enough to the equilibrium state. For asymptotic stability, in addition, x(t)has to converge to the equilibrium state as I-o0. If an equilibrium state is globally asymptotically stable, then x(n) converges to the equilibrium state regard less of how the initial state x, is selected These stability concepts are called internal, because they stable roperties of the state of th illustrated in Fig. 1 FIGURE 12.1 Stability concepts. Source: F. Szi In the electrical engineering literature, sometimes our darovszky and AT Bahill, Linear Systems Theory Boca tability definition is called marginal stability, and our Raton, Fla. CRC Press, 1992, P 168. With permission. asymptotic stability is called stability. e 2000 by CRC Press LLC
© 2000 by CRC Press LLC 12.2 Using the State of the System to Determine Stability The stability of a system is defined with respect to a given equilibrium point in state space. If the initial state x0 is selected at an equilibrium state x of the system, then the state will remain at x for all future time. When the initial state is selected close to an equilibrium state, the system might remain close to the equilibrium state or it might move away. In this section we introduce conditions that guarantee that whenever the system starts near an equilibrium state, it remains near it, perhaps even converging to the equilibrium state as time increases. For simplicity, only time-invariant systems are considered in this section. Time-variant systems are discussed in Section 12.5. Continuous, time-invariant systems have the form (12.1) and discrete, time-invariant systems are modeled by the difference equation (12.2) Here we assume that f: X Æ Rn, where X Õ Rn is the state space. We also assume that function f is continuous; furthermore, for arbitrary initial state x0 Œ X, there is a unique solution of the corresponding initial value problem x(t0) = x0, and the entire trajectory x(t) is in X. Assume furthermore that t0 denotes the initial time period of the system. It is also known that a vector x Œ X is an equilibrium state of the continuous system, Eq. (12.1), if and only if f(x) = 0, and it is an equilibrium state of the discrete system, Eq. (12.2), if and only if x = f(x). In this chapter the equilibrium of a system will always mean the equilibrium state, if it is not specified otherwise. In analyzing the dependence of the state trajectory x(t) on the selection of the initial state x0 nearby the equilibrium, the following stability types are considered. Definition 12.1 1. An equilibrium state x is stable if there is an e0 > 0 with the following property: For all e1, 0 0 such that if || x – x0 || t0. 2. An equilibrium state x is asymptotically stable if it is stable and there is an e > 0 such that whenever || x – x0 || < e, then x(t) Æ x as t Æ •. 3. An equilibrium state x is globally asymptotically stable if it is stable and with arbitrary initial state x0 Œ X, x(t) Æ x as t Æ •. The first definition says an equilibrium state x is stable if the entire trajectory x(t) is closer to the equilibrium state than any small e1, if the initial state x0 is selected close enough to the equilibrium state. For asymptotic stability, in addition, x(t) has to converge to the equilibrium state as t Æ •. If an equilibrium state is globally asymptotically stable, then x(t) converges to the equilibrium state regardless of how the initial state x0 is selected. These stability concepts are called internal, because they represent properties of the state of the system. They are illustrated in Fig. 12.1. In the electrical engineering literature, sometimes our stability definition is called marginal stability, and our asymptotic stability is called stability. x˙( )t = f(x(t)) x(t + 1) = f(x(t)) FIGURE 12.1 Stability concepts. (Source: F. Szidarovszky and A.T. Bahill, Linear Systems Theory, Boca Raton, Fla.: CRC Press, 1992, p. 168. With permission.)
2.3 Lyapunov Stability Theory Assume that x is an equilibrium state of a continuous or discrete system, and let Q2 denote a subset of the state pace X such that x∈ Definition 12.2 A real-valued function V defined on Q2 is called a Lyapunov function, if 2. V has a unique global minimum at x with respect to all other points in $2 for any state trajectory x(t) contained in Q2, v(x(n)) is nonincreasing in t The Lyapunov function can be interpreted as the generalization of the energy function in electrical systems The first requirement simply means that the graph of V has no discontinuities. The second requirement mear that the graph of V has its lowest point at the equilibrium, and the third requirement generalizes the well- known fact of electrical systems, that the energy in a free electrical system with resistance always decreases, unless the system is at rest. Theorem 12.1 Assume that there exists a Lyapunov function V on the spherical region x||x-x‖0 is given; furthermore Q2S X. Then the equilibrium state is stable Theorem 12. 2 Assume that in addition to the conditions of Theorem 12. 1, the Lyapunov function V(x(o)is strictly decreasing in t, unless x(t)=x. Then the equilibrium state is asymptotically stable Theorem 12.3 Assume that the Lyapunov function is defined on the entire state space X, v(x( o))is strictly decreasing in t unless x(r)=x; furthermore, V(x) tends to infinity as any component of x gets arbitrarily large in magnitude Then the equilibrium state is globally asymptotically stable Consider the differential equation The stability of the equilibrium state(1/@, 0)can be verified directly by using Theorem 12. 1 without computing the solution. Select the Lyapunov function V(x)=(x-x)(x-x)=‖ where the euclidian norm is used This is continuous in x; furthermore, it has its minimal(zero)value at x=x. Therefore, to establish the stability the equilibrium state we have to show only that v(x(n)) is decreasing. Simple differentiation shows that (x(t)=2(x-x)2·文=2(x-x)(Ax+b) c 2000 by CRC Press LLC
© 2000 by CRC Press LLC 12.3 Lyapunov Stability Theory Assume that x is an equilibrium state of a continuous or discrete system, and let W denote a subset of the state space X such that x Œ W. Definition 12.2 A real-valued function V defined on W is called a Lyapunov function, if 1. V is continuous; 2. V has a unique global minimum at x with respect to all other points in W; 3. for any state trajectory x(t) contained in W, V(x(t)) is nonincreasing in t. The Lyapunov function can be interpreted as the generalization of the energy function in electrical systems. The first requirement simply means that the graph of V has no discontinuities. The second requirement means that the graph of V has its lowest point at the equilibrium, and the third requirement generalizes the wellknown fact of electrical systems, that the energy in a free electrical system with resistance always decreases, unless the system is at rest. Theorem 12.1 Assume that there exists a Lyapunov function V on the spherical region (12.3) where e0 > 0 is given; furthermore W Õ X. Then the equilibrium state is stable. Theorem 12.2 Assume that in addition to the conditions of Theorem 12.1, the Lyapunov function V(x(t)) is strictly decreasing in t, unless x(t) = x. Then the equilibrium state is asymptotically stable. Theorem 12.3 Assume that the Lyapunov function is defined on the entire state space X, V(x(t)) is strictly decreasing in t unless x(t) = x; furthermore, V(x) tends to infinity as any component of x gets arbitrarily large in magnitude. Then the equilibrium state is globally asymptotically stable. Example 12.1 Consider the differential equation The stability of the equilibrium state (1/w, 0)T can be verified directly by using Theorem 12.1 without computing the solution. Select the Lyapunov function where the Euclidian norm is used. This is continuous in x; furthermore, it has its minimal (zero) value at x = x . Therefore, to establish the stability of the equilibrium state we have to show only that V(x(t)) is decreasing. Simple differentiation shows that W = {x * ** x - x ** < e } 0 ˙ x = x - Ê Ë Á ˆ ¯ ˜ + Ê Ë Á ˆ ¯ ˜ 0 0 0 1 w w V T (x) = - (x x) (x - x) = ** x x - ** 2 2 d dt V t T T (x( )) = 2 2 (x - x) × x˙ = (x - x) (Ax + b)
That is, with=(x, x,)T 2 =2(Ox1x2-x2-Ox1x2+x2)=0 Therefore, function V(x(n))is a constant, which is a(not strictly)decreasing function. That is, all conditions of Theorem 12.1 are satisfied, which implies the stability of the equilibrium state. Theorems 12.1, 12.2, and 12.3 guarantee, respectively, the stability, asymptotic stability, and global asymptotic stability of the equilibrium state, if a Lyapunov function is found. Failure to find such a Lyapunov function does not mean that the system is unstable or that the stability is not asymptotic or globally asymptotic. It only means that you were not clever enough to find a Lyapunov function that proved stability 12. 4 Stability of Time-Invariant Linear Systems This section is divided into two subsections. In the first subsection the stability of linear time-invariant systems given in state-space notation is analyzed. In the second subsection, methods based on transfer functions are discussed Stability analysis with State-Space Notation Consider the time-invariant continuous linear system Ax+b (12.4) and the time-invariant discrete linear syste x(t +1)=Ax(t)+b Assume that x is an equilibrium state, and let o(t, t,)denote the fundamental matrix. Theorem 12. 4 1. The equilibrium state x is stable if and only if o(t, to)is bounded for t2 to 2. The equilibrium state x is asymptotically stable if and only if o(t, t )is bounded and tends to zero as t-oo. e use the symbol s to denote complex frequency, i. e, s=0+ jo. For specific values of s, such as eigenvalues and poles, we use the symbol. Theorem 12.5 I.If for at least one eigenvalue of A, Re n >0 (or a>1 for discrete systems), then the system is unstable 2. Assume that for all eigenvalues a, of A, Re i, s 0 in the continuous case(or s 1 in the discrete case), and all eigenvalues with the property Re =0(or =1)have single multiplicity; then the equilibrium 3. The stability is asymptotic if and only if for all i, Re l, <0(or <1 e 2000 by CRC Press LLC
© 2000 by CRC Press LLC with That is, with x = (x1, x2)T , Therefore, function V(x(t)) is a constant, which is a (not strictly) decreasing function. That is, all conditions of Theorem 12.1 are satisfied, which implies the stability of the equilibrium state. Theorems 12.1, 12.2, and 12.3 guarantee,respectively, the stability, asymptotic stability, and global asymptotic stability of the equilibrium state, if a Lyapunov function is found. Failure to find such a Lyapunov function does not mean that the system is unstable or that the stability is not asymptotic or globally asymptotic. It only means that you were not clever enough to find a Lyapunov function that proved stability. 12.4 Stability of Time-Invariant Linear Systems This section is divided into two subsections. In the first subsection the stability of linear time-invariant systems given in state-space notation is analyzed. In the second subsection, methods based on transfer functions are discussed. Stability Analysis with State-Space Notation Consider the time-invariant continuous linear system (12.4) and the time-invariant discrete linear system (12.5) Assume that x is an equilibrium state, and let f(t,t0) denote the fundamental matrix. Theorem 12.4 1. The equilibrium state x is stable if and only if f(t,t0) is bounded for t ³ t0. 2. The equilibrium state x is asymptotically stable if and only if f(t,t0) is bounded and tends to zero as t Æ •. We use the symbol s to denote complex frequency, i.e., s = s + jw. For specific values of s, such as eigenvalues and poles, we use the symbol l. Theorem 12.5 1. If for at least one eigenvalue of A, Re li > 0 (or *li * > 1 for discrete systems), then the system is unstable. 2. Assume that for all eigenvalues li of A, Re li £ 0 in the continuous case (or *li * £ 1 in the discrete case), and all eigenvalues with the property Re li = 0 (or *li *= 1) have single multiplicity; then the equilibrium state is stable. 3. The stability is asymptotic if and only if for all i, Re li < 0 (or *li * < 1). A = b - Ê Ë Á ˆ ¯ ˜ = Ê Ë Á ˆ ¯ ˜ 0 0 0 1 w w and d dt V t x x x x x x x x x x ( ( )) , ( ) x = - Ê Ë Á ˆ ¯ ˜ - + Ê Ë Á ˆ ¯ ˜ = - - + = 2 1 1 2 0 1 2 2 1 1 2 2 1 2 2 w w w w w x˙ = + Ax b x(t + = 1) Ax(t) + b
emark 1. Note that Part 2 gives only sufficient conditions for the stability of the equilibrium state. As the following example shows, these conditions are not necessary. le12.2 Consider first the continuous system x=Ox, where O is the zero matrix. Note that all constant functions x(t) = x are solutions and also equilibrium states. Since φ(t,t0) is bounded(being independent of t), all equilibrium states are stable, but o has only one eigenvalue M=0 with zero real part and multiplicity n, where n is the order of the system Consider next the discrete systems x(t 1)=Ix(o), when all constant functions x(o)=x are also solutio and equilibrium states. Furthermore, (,t0)=A I which is obviously bounded. Therefore, all equilibrium states are stable, but the condition of Part 2 of the theorem is violated again, since M= l with unit absolute value having a multiplicity n. Remark 2. The following extension of Theorem 12.5 can be proven. The equilibrium state is stable if and only if for all eigenvalues of A,Reλ:≤0or内≤1), and ifλ M is a repeated eigenvalue of A such that Re= 0 or =1), then the size of each block containing 2 in the Jordan canonical form of A is 1 X1. if the same holds for the equilibrium states of the corres niog homogeneous equation ss Remark 3. The equilibrium states of inhomogeneous equations are stable or asymptotically stable if and only Example 12.3 0 the stability of which was analyzed earlier in Example 12.1 by using the Lyapunov function method. The characteristic polynomial of the coefficient matrix is therefore, the eigenvalues are M=jo and n2=-jo Both eigenvalues have single multiplicities, and Re M1= Re of Part 3 do not hold. Consequenty the system is not asymptotically stable m state is st n2=0. Hence, the conditions of Part 2 are satisfied, and therefore the equilibrium state is stable. The conditions If a time-invariant system is nonlinear, then the Lyapunov method is the most popular choice for stability alysis. If the system is linear, then the direct application of Theorem 12.5 is more attractive, since the igenvalues of the coefficient matrix A can be obtained by standard methods. In addition, several condition are known from the literature that guarantee the asymptotic stability of time-invariant discrete and continuou systems even without computing the eigenvalues For examining asymptotic stability, linearization is an alternative approach to the Lyapunov method as is shown here. Consider the time- invariant continuous and discrete syste x(t)=f(x(t)) c 2000 by CRC Press LLC
© 2000 by CRC Press LLC Remark 1. Note that Part 2 gives only sufficient conditions for the stability of the equilibrium state. As the following example shows, these conditions are not necessary. Example 12.2 Consider first the continuous system x· = Ox, where O is the zero matrix. Note that all constant functions x(t) º x are solutions and also equilibrium states. Since is bounded (being independent of t), all equilibrium states are stable, but O has only one eigenvalue l1 = 0 with zero real part and multiplicity n, where n is the order of the system. Consider next the discrete systems x(t + 1) = Ix(t), when all constant functions x(t) º x are also solutions and equilibrium states. Furthermore, which is obviously bounded. Therefore, all equilibrium states are stable, but the condition of Part 2 of the theorem is violated again, since l1= 1 with unit absolute value having a multiplicity n. Remark 2. The following extension of Theorem 12.5 can be proven. The equilibrium state is stable if and only if for all eigenvalues of A, Re li £ 0 (or *li * £ 1), and if li is a repeated eigenvalue of A such that Re li = 0 (or *li * = 1), then the size of each block containing li in the Jordan canonical form of A is 1 3 1. Remark 3. The equilibrium states of inhomogeneous equations are stable or asymptotically stable if and only if the same holds for the equilibrium states of the corresponding homogeneous equations. Example 12.3 Consider again the continuous system the stability of which was analyzed earlier in Example 12.1 by using the Lyapunov function method. The characteristic polynomial of the coefficient matrix is therefore, the eigenvalues are l1 = jw and l2 = –jw. Both eigenvalues have single multiplicities, and Re l1 = Re l2 = 0. Hence, the conditions of Part 2 are satisfied, and therefore the equilibrium state is stable. The conditions of Part 3 do not hold. Consequently, the system is not asymptotically stable. If a time-invariant system is nonlinear, then the Lyapunov method is the most popular choice for stability analysis. If the system is linear, then the direct application of Theorem 12.5 is more attractive, since the eigenvalues of the coefficient matrix A can be obtained by standard methods. In addition, several conditions are known from the literature that guarantee the asymptotic stability of time-invariant discrete and continuous systems even without computing the eigenvalues. For examining asymptotic stability, linearization is an alternative approach to the Lyapunov method as is shown here. Consider the time-invariant continuous and discrete systems f( , ) ( ) t t e t t 0 0 = = O - I f(t t, ) t t t t 0 0 0 = = = - - A I I ˙ x = - Ê Ë Á ˆ ¯ ˜ + Ê Ë Á ˆ ¯ ˜ 0 0 0 1 w w x j w w (s) w s s = s - - - Ê Ë Á ˆ ¯ ˜ det = + 2 2 x˙( )t = f x( (t))
x(t+1)=f(x(t) Let j(x)denote the Jacobian of f(x), and let x be an equilibrium state of the system. It is known that the method of linearization around the equilibrium state results in the time-invariant linear systems x6()=J(x)x(t) x6(t+1)=J(x)x(t) where x(t)=x(o)-x. It is also known from the theory of ordinary differential equations that the asymptotic tability of the zero vector in the linearized system implies the asymptotic stability of the equilib For continuous systems the following result has a special importance The equilibrium state of a continuous system [Eq.(12.4) is asymptotically stable if and only if equation A Q+QA has positive definite solution Q with some positive definite matrix M We note that in practical applications the identity matrix is almost always selected for M. An initial stability check is provided by the following result Theorem 12 Let o()=A"+P2-1+..+P,+ Po be the characteristic polynomial of matrix A Assume that all eigenvalues of matrix A have negative real parts. Then P: >0(i=0, 1,.,n-1) Corollary. If any of the coefficients p is negative or zero, the equilibrium state of the system with coefficient matrix A cannot be asymptotically stable. However, the conditions of the theorem do not imply that the eigenvalues of A have negative real part the characteristic polynominal is ((s)=s2+0. Since the coefficient of s is zero, the system of Example 12.3 tically stable The Transfer Function Approach The transfer function of the time invariant linear continuous system (12.7) e 2000 by CRC Press LLC
© 2000 by CRC Press LLC and Let J(x) denote the Jacobian of f(x), and let x be an equilibrium state of the system. It is known that the method of linearization around the equilibrium state results in the time-invariant linear systems and where xd(t) = x(t) – x. It is also known from the theory of ordinary differential equations that the asymptotic stability of the zero vector in the linearized system implies the asymptotic stability of the equilibrium state x in the original nonlinear system. For continuous systems the following result has a special importance. Theorem 12.6 The equilibrium state of a continuous system [Eq. (12.4)] is asymptotically stable if and only if equation (12.6) has positive definite solution Q with some positive definite matrix M. We note that in practical applications the identity matrix is almost always selected for M. An initial stability check is provided by the following result. Theorem 12.7 Let j(l) = l n + pn–1 ln–1 + . . . + p1l + p0 be the characteristic polynomial of matrix A.Assume that all eigenvalues of matrix A have negative real parts. Then pi > 0 (i = 0, 1,..., n – 1). Corollary. If any of the coefficients pi is negative or zero, the equilibrium state of the system with coefficient matrix A cannot be asymptotically stable. However, the conditions of the theorem do not imply that the eigenvalues of A have negative real parts. Example 12.4 For matrix the characteristic polynominal is j(s) = s 2 + w2 . Since the coefficient of s1 is zero, the system of Example 12.3 is not asymptotically stable. The Transfer Function Approach The transfer function of the time invariant linear continuous system (12.7) x(t + 1) = f(x(t)) x˙ ( ) Jxx ( ) ( ) d d t = t x Jxx d d (t + 1) = ( ) (t) A Q QA M T + = - A = - Ê Ë Á ˆ ¯ ˜ 0 0 w w x˙ Ax Bu y Cx = + =
and that of the time invariant linear discrete system x(t+1)= Ax(t)+ Bu(t) (12.8) (t )= Cx(t) have the common form TF(S)=C(Sl-A) If both the input and output are single, ther Y(s) or in the familiar electrical engineering notation KG(s) TF(s) (12.9) 1+ KG(S)H(s) where K is the gain term in the forward loop, G(s) represents the dynamics of the forward loop, or the plant, and H(s)models the dynamics in the feedback loop. We note that in the case of continuous systems s is the variable of the transfer function, and for discrete systems the variable is denoted by z. After the Second World War systems and control theory flourished. The transfer function representation was the most popular representation for systems. To determine the stability of a system we merely had to factor the denominator of the transfer function(12.9)and see if all of the poles were in the left half of the complex frequency plane. However, with manual techniques, factoring polynomials of large order is difficult. So engi leers, being naturally lazy people, developed several ways to determine the stability of a system without factoring the polynomials [Dorf, 1992]. First, we have the methods of Routh and Hurwitz, developed a century ago, that looked at the coefficients of the characteristic polynomial. These methods showed whether the system was stable or not, but they did not show how close the system was to being stable What we want to know is for what value of gain, K, and at what frequency, a, will the denominator of the transfer function(12.9)become zero Or, when will KGH =-1, meaning, when will the magnitude of KGH equal I with a phase angle of -180 degrees? These parameters can be determined easily with a Bode diagram. Construct a Bode diagram for KGH of the system, look at the frequency where the phase angle equals -180 degrees, and look up at the magnitude plot. If it is smaller than 1.0, then the system is stable If it is larger than 1.0, then the system is unstable. Bode diagram techniques are discussed in Chapter 11 would be encountered by a signal in one loon P transfer function of the system, because it is the effect that The quantity KG(s)H(s) is called the open p around the system if the feedback loop were artificially opened Bahill, 1981. To gain some intuition, think of a closed-loop negative feedback system. Apply a small sinusoid at frequency o to the input. Assume that the gain around the loop, KGH, is 1 or more, and that the phase lag is 180 degrees. The summing junction will flip over the fed back signal and add it to the original signal. The result is a signal that is bigger than what came in. This signal will circulate around this loop, getting bigger and bigger until the real system no longer matches the model. This is what we call instability. The question of stability can also be answered with Nyquist diagrams. They are related to Bode diagram but they give more information. A simple way to construct a Nyquist diagram is to make a polar plot on the complex frequency plane of the Bode diagram. Simply stated, if this contour encircles the -l point in the omplex frequency plane, then the system is unstable. The two advantages of the nyquist technique are(1)in c 2000 by CRC Press LLC
© 2000 by CRC Press LLC and that of the time invariant linear discrete system (12.8) have the common form If both the input and output are single, then or in the familiar electrical engineering notation (12.9) where K is the gain term in the forward loop, G(s) represents the dynamics of the forward loop, or the plant, and H(s) models the dynamics in the feedback loop. We note that in the case of continuous systems s is the variable of the transfer function, and for discrete systems the variable is denoted by z. After the Second World War systems and control theory flourished. The transfer function representation was the most popular representation for systems. To determine the stability of a system we merely had to factor the denominator of the transfer function (12.9) and see if all of the poles were in the left half of the complex frequency plane. However, with manual techniques, factoring polynomials of large order is difficult. So engineers, being naturally lazy people, developed several ways to determine the stability of a system without factoring the polynomials [Dorf, 1992]. First, we have the methods of Routh and Hurwitz, developed a century ago, that looked at the coefficients of the characteristic polynomial. These methods showed whether the system was stable or not, but they did not show how close the system was to being stable. What we want to know is for what value of gain, K, and at what frequency, w, will the denominator of the transfer function (12.9) become zero. Or, when will KGH = –1, meaning, when will the magnitude of KGH equal 1 with a phase angle of –180 degrees? These parameters can be determined easily with a Bode diagram. Construct a Bode diagram for KGH of the system, look at the frequency where the phase angle equals –180 degrees, and look up at the magnitude plot. If it is smaller than 1.0, then the system is stable. If it is larger than 1.0, then the system is unstable. Bode diagram techniques are discussed in Chapter 11. The quantity KG(s)H(s) is called the open-loop transfer function of the system, because it is the effect that would be encountered by a signal in one loop around the system if the feedback loop were artificially opened [Bahill, 1981]. To gain some intuition, think of a closed-loop negative feedback system. Apply a small sinusoid at frequency w to the input. Assume that the gain around the loop, KGH, is 1 or more, and that the phase lag is 180 degrees. The summing junction will flip over the fed back signal and add it to the original signal. The result is a signal that is bigger than what came in. This signal will circulate around this loop, getting bigger and bigger until the real system no longer matches the model. This is what we call instability. The question of stability can also be answered with Nyquist diagrams. They are related to Bode diagrams, but they give more information. A simple way to construct a Nyquist diagram is to make a polar plot on the complex frequency plane of the Bode diagram. Simply stated, if this contour encircles the –1 point in the complex frequency plane, then the system is unstable. The two advantages of the Nyquist technique are (1) in x Ax Bu y Cx ( ) ( ) ( ) ( ) ( ) t t t t t + = + = 1 TF( )s = - C(sI A) B -1 TF Y U ( ) ( ) ( ) s s s = TF G G H ( ) ( ) ( ) ( ) s K s K s s = 1 +
addition to the information on Bode diagrams, there are about a dozen rules that can be used to help construct Nyquist diagrams, and(2) Nyquist diagrams handle bizarre systems better, as is shown in the following rigorous statement of the Nyquist stability criterion. The number of clockwise encirclements minus the number of counterclockwise encirclements of the point s=-1+jo by the Nyquist plot of KG(s)H(s) is equal to the number of poles of Y(s)/U(s)minus the number of poles of KG(s)H(s) in the right half of the s-plane The root-locus technique was another popular technique for assessing stability. It furthermore allowed the engineer to see the effects of small changes in the gain, K, on the stability of the system. The root-locus diagram shows the location in the s-plane of the poles of the closed-loop transfer function, Y(S)/U(s). All branches of the root-locus diagram start on poles of the open-loop transfer function, KGH, and end either on zeros of the open-loop transfer function, KGH, or at infinity. There are about a dozen rules to help draw these trajectories. The root-locus technique is discussed in Chapter 93.4 We consider all these techniques to be old fashioned. They were developed to help answer the question of stability without factoring the characteristic polynomial. However, many computer programs are currently available that factor polynomials. We recommend that engineers merely buy one of these computer packages and find the roots of the closed-loop transfer function to assess the stability of a system. The poles of a system are defined as all values of s such that sI-A is singular. The poles of a closed-loop transfer function are exactly the same as the eigenvalues of the system: engineers prefer the term poles and the symbol s, and mathematicians prefer the term eigenvalues and the symbol n. We will use s for complex frequency Sometimes, some poles could be canceled in the rational function form of TF(s)so that they would not be explicitly shown. However, even if some poles could be canceled by zeros, we still have to consider all poles in the following criteria which is the statement of Theorem 12.5. The equilibrium state of the continuous system [Eq(12.7)] with constant input is unstable if at least one pole has a positive real part, and is stable if all poles of TF(s)have nonpositive real parts and all poles with zero real parts are single. The equilibrium state is asymptotically stable if and only if all poles of TF(s) have negative real parts; that is, all poles are in the left half of the s-plane. Similarly, the equilibrium state of the discrete system [Eq (12. 8] with constant input is unstable if the absolute value of at least one pole is greater than one, and is stable if all poles of TF(z) have absolute values less than or equal to one and all poles with unit absolute values are single. The equilibrium state is asymptotically stable if and only if all poles of TF(z) have absolute values less than one; that is, the poles all inside the unit circle of the z-plane Consider again the system which was discussed earlier. Assume that the output equation has the form y=(1,1)x Then TF(s)=5+o The poles are jo and jo, which have zero real parts; that is, they are on the imaginary axis of the e equilibrium state is stable but not asymptotically stable. A system such would produce constant amplitude sinusoids at frequency @. So it seems natural to assume that such would be used to build sinusoidal signal generators and to model oscillating systems. However, this is not the case, because(1)zero resistance circuits are hard to make; therefore, most function generators use other e 2000 by CRC Press LLC
© 2000 by CRC Press LLC addition to the information on Bode diagrams, there are about a dozen rules that can be used to help construct Nyquist diagrams, and (2) Nyquist diagrams handle bizarre systems better, as is shown in the following rigorous statement of the Nyquist stability criterion. The number of clockwise encirclements minus the number of counterclockwise encirclements of the point s = –1 + j 0 by the Nyquist plot of KG(s)H(s) is equal to the number of poles of Y(s)/U(s) minus the number of poles of KG(s)H(s) in the right half of the s-plane. The root-locus technique was another popular technique for assessing stability. It furthermore allowed the engineer to see the effects of small changes in the gain, K, on the stability of the system. The root-locus diagram shows the location in the s-plane of the poles of the closed-loop transfer function, Y(s)/U(s). All branches of the root-locus diagram start on poles of the open-loop transfer function, KGH, and end either on zeros of the open-loop transfer function, KGH, or at infinity. There are about a dozen rules to help draw these trajectories. The root-locus technique is discussed in Chapter 93.4. We consider all these techniques to be old fashioned. They were developed to help answer the question of stability without factoring the characteristic polynomial. However, many computer programs are currently available that factor polynomials. We recommend that engineers merely buy one of these computer packages and find the roots of the closed-loop transfer function to assess the stability of a system. The poles of a system are defined as all values of s such that sI – A is singular. The poles of a closed-loop transfer function are exactly the same as the eigenvalues of the system: engineers prefer the term poles and the symbol s, and mathematicians prefer the term eigenvalues and the symbol l.We will use s for complex frequency and l for specific values of s. Sometimes, some poles could be canceled in the rational function form of TF(s) so that they would not be explicitly shown. However, even if some poles could be canceled by zeros, we still have to consider all poles in the following criteria which is the statement of Theorem 12.5. The equilibrium state of the continuous system [Eq. (12.7)] with constant input is unstable if at least one pole has a positive real part, and is stable if all poles of TF(s) have nonpositive real parts and all poles with zero real parts are single. The equilibrium state is asymptotically stable if and only if all poles of TF(s) have negative real parts; that is, all poles are in the left half of the s-plane. Similarly, the equilibrium state of the discrete system [Eq. (12.8)] with constant input is unstable if the absolute value of at least one pole is greater than one, and is stable if all poles of TF(z) have absolute values less than or equal to one and all poles with unit absolute values are single. The equilibrium state is asymptotically stable if and only if all poles of TF(z) have absolute values less than one; that is, the poles are all inside the unit circle of the z-plane. Example 12.5 Consider again the system which was discussed earlier. Assume that the output equation has the form Then The poles are jw and –jw, which have zero real parts; that is, they are on the imaginary axis of the s-plane. Since they are single poles, the equilibrium state is stable but not asymptotically stable. A system such as this would produce constant amplitude sinusoids at frequency w. So it seems natural to assume that such systems would be used to build sinusoidal signal generators and to model oscillating systems. However, this is not the case, because (1) zero resistance circuits are hard to make; therefore, most function generators use other ˙ x = x - Ê Ë Á ˆ ¯ ˜ + Ê Ë Á ˆ ¯ ˜ 0 0 0 1 w w y = (1, 1)x TF(s) s s = + + w w 2 2
chniques to produce sinusoids; and(2) such systems are not good models for oscillating systems, because most real-world oscillating systems (i. e, biological systems) have energy dissipation elements in them More generally, real-world function generators are seldom made from closed-loop feedback control systems with 180 degrees of phase shift, because(1)it would be difficult to get a broad range of frequencies and several waveforms from such systems,(2) precise frequency selection would require expensive high-precision compo lents, and(3)it would be difficult to maintain a constant frequency in such circuits in the face of changing mperatures and power supply variations. Likewise, closed-loop feedback control systems with 180 degrees of phase shift are not good models for oscillating biological systems, because most biological systems oscillate because of nonlinear network properties A special stability criterion for single-input, single-output time-invariant continuous systems will be intro- duced next. Consider the system Ax+b d (12.10) where A is an nx n constant matrix and b and c are constant n-dimensional vectors The transfer function of this system is TF(s)=C(SI-A which is obviously a rational function of s. Now let us add negative feedback around this system so that u ky, where k is a constant. The resulting system can be described by the differential equation X= Ax +kbcX =(A + kbc x (12.11) The transfer function of this feedback system is TE(S) TFS) To help show the connection between the asymptotic stability of systems(1210)and(12.11), we introduce the Definition 12.3 Let r(s) be a rational function of s. Then the locus of points L(r)=(a+jbla= Re(r(iv), b= Im(r(v),V E RE is called the response diagram of r. Note that L(r)is the image of the imaginary line Re(s)=o under the mapping r We shall assume that L(r) is bounded, which is the case if and only if the degree of the denominator is not less than that of the numerator and r has no poles on the line re(s)=o Theorem 12. 8 The Nyquist stability criterion. Assume that TF, has a bounded response diagram L(TF). If TF, has v poles in the right half of the s-plane, where Re(s)>0, then H has p +v poles in the right half of the s-plane where Re(s)>0 if the point 1/k+.0 is not on L(TF,), and I(TE) encircles 1/k+. 0 p times in the clockwise sense Corollary. Assume that system(12.10)is asymptotically stable with constant input and that I(TFi)is bounded and traversed in the direction of increasing v and has the point 1/k +j.0 on its left. Then the feedback system (12.11)is also asymptotically stable c 2000 by CRC Press LLC
© 2000 by CRC Press LLC techniques to produce sinusoids; and (2) such systems are not good models for oscillating systems, because most real-world oscillating systems (i.e., biological systems) have energy dissipation elements in them. More generally, real-world function generators are seldom made from closed-loop feedback control systems with 180 degrees of phase shift, because (1) it would be difficult to get a broad range of frequencies and several waveforms from such systems, (2) precise frequency selection would require expensive high-precision components, and (3) it would be difficult to maintain a constant frequency in such circuits in the face of changing temperatures and power supply variations. Likewise, closed-loop feedback control systems with 180 degrees of phase shift are not good models for oscillating biological systems, because most biological systems oscillate because of nonlinear network properties. A special stability criterion for single-input, single-output time-invariant continuous systems will be introduced next. Consider the system (12.10) where A is an n ¥ n constant matrix, and b and c are constant n-dimensional vectors. The transfer function of this system is which is obviously a rational function of s. Now let us add negative feedback around this system so that u = ky, where k is a constant. The resulting system can be described by the differential equation (12.11) The transfer function of this feedback system is (12.12) To help show the connection between the asymptotic stability of systems (12.10) and (12.11), we introduce the following definition. Definition 12.3 Let r(s) be a rational function of s. Then the locus of points is called the response diagram of r. Note that L(r) is the image of the imaginary line Re(s) = 0 under the mapping r. We shall assume that L(r) is bounded, which is the case if and only if the degree of the denominator is not less than that of the numerator and r has no poles on the line Re(s) = 0. Theorem 12.8 The Nyquist stability criterion. Assume that TF1 has a bounded response diagram L(TF1). If TF1 has n poles in the right half of the s-plane, where Re(s) > 0, then H has r + n poles in the right half of the s-plane where Re(s) > 0 if the point 1/k + j · 0 is not on L(TF1), and L(TF1) encircles 1/k + j · 0 r times in the clockwise sense. Corollary. Assume that system (12.10) is asymptotically stable with constant input and that L(TF1) is bounded and traversed in the direction of increasing n and has the point 1/k + j · 0 on its left. Then the feedback system (12.11) is also asymptotically stable. x˙ = + Ax bu and y = c x T TF s s 1 1 ( ) = - ( ) - c I A b T x˙ = + Ax k k bc x = (A + bc )x T T TF TF TF ( ) ( ) ( ) s s k s = - 1 1 1 L(r) = {a + jb *a = Re(r(jv)), b = Im(r(jv)), v Œ R}