THE ABSOLUTE STABILITY OF NONLINEAR SYSTEMS by TE-LUNG CHANG B.A.Sc. in E.E., Cheng Kung University, 1966 A THESIS SUBMITTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in the Department of Electrical Engineering We accept this thesis as conforming to the required standard Research Supervisor Members of Committee Acting Head of Department Members of the Department of Electrical Engineering THE UNIVERSITY OF BRITISH COLUMBIA February, 1970 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree tha permission for extensive copying of this thesis for scholarly purposes may be granted by the Head of my Department or by his representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of The University of British Columbia Vancouver 8, Canada ABSTRACT This thesis is in two parts, both considering the absolute stability of nonlinear systems. In the first two chapters the stability of certain classes of nonlinear time invariant systems involving several nonlinearities is considered. A number of graphical methods are given for testing the stability of these systems. The graphical tests are equivalent to a weakened form of the Popov criterion. The third chapter derives a stability condition for nonlinear systems involving a linear time-varying gain. The time-varying gain is assumed to satisfy conditions on its magnitude and rate of change. ii "TABLE OF CONTENTS Page LIST OF ILLUSTRATIONS v ACKNOWLEDGEMENT '. . vii 1. TIME INVARIANT NONLINEAR FEEDBACK SYSTEMS AND TIME VARYING NONLINEAR FEEDBACK'SYSTEMS 1 1.1 Introduction1.2 Lyapunov Methods 7 1.3 The Popov Criterion 8 1.4 L2 Stability 10 2. THE ABSOLUTE STABILITY OF A TIME INVARIANT NONLINEAR FEEDBACK SYSTEM 1 2.1 Introduction. 12.2 A Graphical Test of the Absolute Stability of a Series System with Nonlinearities and Identical Transfer Functions 11 2.2.1 Main Method-Common Popov Line 13 2.2.2 Simplification in Particular Cases 16 2.2.3 Use of the Modified Nichols Chart to Obtain Q 19 2.2.4 An Analogue-Computer Technique 23 2.3 The Absolute Stability of the Series Nonlinear System with Different Transfer Functions 27 2.4 The Absolute Stability of a Parallel Nonlinear System 30 3. THE ABSOLUTE STABILITY OF A TIME VARYING FEEDBACK SYSTEM WITH MONOTONIC NONLINEARITIES 34 3.1 The Absolute Stability of the Single-Loop Time Varying Nonlinear Feedback System 34 3.2 The Main Results 35 3.2.1 Theorem 13.2.2 A Speical Case of Theorem 1 36 3.3 Proof of Main Results •3.4 Examples 41 *in-i Page 3.5 The Absolute Stability of a System with Multiple Nonlinearities and Time Varying Gains 46 3.6 Theorem 2 49 3.7 Proof of Theorem 2 50 3.8 Example 51 4. CONCLUSIONS 3 APPENDIX 1 5 APPENDIX 2 , 57 APPENDIX 3 58 APPENDIX 4 9 REFERENCES 60 iv LIST OF ILLUSTRATIONS Page 1.1 General time invariant nonlinear feedback system 2 1.2 General time varying nonlinear feedback system1.3 Series system with m nonlinearities 4 1.4 Parallel system with m nonlinearities... 5 1.5 Internal feedback system with m nonlinearities 5 1.5 Multi-circuit system with m nonlinearities 6 1.7 Popov criterion 9 2.1 Graphical criterion for the time invariant series system with 5 identical linear transfer functions. 2.2 Modified Nyquist loci and the common Popov line for the series system, o nf \ 32(s+0.25) . _ • ' . Nm=3' G(s) = (s+l)(s+2)(s+4)' FOR EXAMPLE 2>1•• 15 2.3 Modified Nyquist locus for i=l, and the corresponding Popov line for the series system, m=3, G(s) = (g+i)(s+5)(s+20) ' ExamPle 2.2.1. 17 2.4 Nyquist plot of (1+0.788S)G(s), where G(s) = (s+1)(S+5)(s+20)' AND a graphical testing the stability of the series system, m=3, for Example 2.2.1 17 2.5 Modified Nyquist locus for i=l, and the corresponding Popov line for the series system, m=4, G(s) = , 1. , ,rw ,9n\ > for Example Example 2.2.2 \...\ \..".\ 18 2.6 Nyquist plot of (1+0.766s)G(s) , where G(s) = (s+]_) (3+5) (S+20) ' AN^ a graphical testing the stability of the series system, m=4, for Example 2.2.2 12.7 Gain-phase plot for G(s) = (S+JQ) (S+5Q) an<^ a ram:i-ly or boundaries for the series system, m=3, for Example 2.3 21 6 30 2.8 Gain-phase plot for G(s) =(S+2Q) (S+50) ANC^ A ^am^y °f .boundaries for the series system, m=4, for Example 2.3 22 580 2.9 Gain-phase plot for G(s) = (s+IQ)(S+50) and a family of boundaries for the series system, m=5, for Example 2.3 22 v Page 2.10A Computer program for R-I generator 24 2.10B Computer program for X-Y producer 5 2.11 The loci plotted by analogue computer 26 2.12 Extension of Popov criterion for the parallel system with m iden tical nonlinearities and m linear transfer functions 33 3.1 Single loop time varying nonlinear feedback system with zero-input 42 3.2 Bode diagram for a compensator l+y(s) = 43.3 The plots of G(s) = T4TV7^\ and G'(s) = (|^)G(s) for (s+l)(s+2) 4+s Example 3.1 44 3.4 Modified Nyquist plot of G(s) = (s-f]j (S+D) ExamP-'-e ^.2 47 3.5 General time varying nonlinear feedback system with zero-input.. 47 vi ACKNOWLEDGEMENT I wish to express my sincere gratitude to Dr. M.S. Davies, the supervisor of this project, for his invaluable guidance, and Dr. E.V. Bohn for reading the manuscript. Grateful acknowledgement is given to the National Research Council for financial support received under NRC Grant A-4148, and to the University of British Columbia for U.B.C. Fellowship awarded since 1968. I would also like to thank Miss Beverly Harasymchuk for typing the manuscript. VH 1 CHAPTER 1 TIME INVARIANT NONLINEAR FEEDBACK SYSTEM AND TIME VARYING NONLINEAR FEEDBACK SYSTEM §1.1 Introduction This thesis considers the absolute stability of the equilibrium position, X = 0, of feedback systems defined by X = AX + BY Y = 0(a) (1.1.1) T a = C X, or, in the time varying case, X = AX + BY Y = 0(a,t) (1.1.2) T a = C X, where X is an n-vector, Y is an m-vector, 0 is an m-vector, A is an nxn con-T stant matrix, B is annxm constant matrix, and C is an mxn constant matrix. In (1.1.1), each element of 0(a) is a nonlinear function of a alone, so that the system (1) is time invariant. In (1.1.2), each element of 0(a,t) is a time varying nonlinear function of both a and t. It is assumed that the i-th element can be separated into a nonlinear part 0 . (a . ) and the time ix varying gain k.(t), where 0.(a.,t) = 0.(a.)k.(t). The system (2) is thus a l ii iii time varying nonlinear feedback system. The transfer matrix TXs).of the linear part is T(s) = CT(sI-A)_1 B, (1.1.3) where I is an nxn unit matrix. The feedback system (1) may be depicted as shown in Fig. 1.1. The forward path consists of an mxm linear time invariant matrix r(s) and the nonlinearity matrix N.L. In Fig. 1.2, the forward path of the time varying nonlinear feedback system (2) consists of the transfer 2 Fig. 1.1 General time invariant nonlinear feedback system r(t) + <rtt) N N.U <P((Tl k(t) r(s) C(t) Fig. 1.2 General time varying nonlinear feedback system 3 matrix F(s), the mxm diagonal nonlinearity matrix N.L. with elements 0T(a, ), 0„(a„),. . • 0 (a ) and the mxm diagonal time-varying gain matrix k(t) 11 2 2 mm with elements k, (t),k (t) ,. . .k (t). 12m A wide variety of systems may be treated by choice of r(s) , some forms of particular interest are: (1) A series system. The forward path consists of the single-input single-output linear time invariant transfer functions separated by amnesic nonlinearities. This is shown in Fig. 1.3. (2) A parallel system. The forward path consists of m-parallel branches, each of which has one nonlinearity in series with one linear time invariant transfer function. This is shown in Fig. 1.4. (3) An internal feedback system. The forward path consists of m-single nonlinear feedback loops. This is shown in Fig. 1.5. (4) A multi-circuit system. Such systems do not fall into the previous classes. Such a system is shown in Fig. 1.6. The elements of the nonlinearity matrix are to be considered in 4 classes. In each it is assumed that 0(a) is a piece-wise continuous, angle valued function of a. (1) A sector nonlinearity; any function which satisfies the condition k. <^< k9, a i 0 la I 0(a) =0, a = 0 where k^ > 0 and k^ may be positive, negative or zero. (2) A first and third quadrant nonlinearity; this' is a special case of (1). 0 < < k < oo k > 0 , a — for all finite nonzero values of a. <ryt) <t>,«r,) Q,ls) <t>i(tTi) Gi(s) (Tmft) <Pm<F»t Fig. 1.3 Series system with m nonlinearities 5 Fig. 1.5 Internal feedback system with m nonlinearities 6 Fig. 1.6 Multi-circuit system with m nonlinearities (3) A monotonic nonlinearity; this belongs to a subclass of (2) , in which it is assumed that d|M>o. do (4) A monotonic odd nonlinearity; this belongs to a subclass of (3) , in which it is further assumed that 0(-a) = -0(a) for all a. Furthermore, the elements of the time varying gain matrix to be considered are assumed positive, bounded and continuous. §1.2 Lyapunov Method Investigations of the stability of such systems were initiated by Lure who proposed a Lyapunov function of the form: T m a. V(X) = X PX + Z 3. / 1 0.(z)dz, (1.2.1) i=l 1 0 1 T where P is an nxn symmetrical positive definite matrix (P = P > 0), all 3. > 0 and the upper limits of the integral terms are the elements of the matrix CTX. (1.2.2) It is obvious that V(X) is a positive function since 0 < 0(a)a < Ka2. (1.2.3) From (1.1.1), the derivative of (1.2.1) is V(X) = XT(APT+PA) +BTPX0(a)'+ XTPB0(a) - 3AX0(a) + BC0(a), (1.2.4) where g is a diagonal constant matrix with elements 3n,$0 >•••6 • I 2 m From the Lyapunov second method, in order to find the sufficient condition of absolute stabiltiy, it is necessary to determine the conditions under which V(X) is negative definite except at the null state, X = 0, where V(X) = 0. In [5], it was shown that the conditions of absolute stability of a time varying nonlinear system may be found from using the Lyapunov func tion (1.2.1) with all 3^ = 0. This method of using Lyapunov functions was further developed by Narendra and Taylor [6] using the modified Lyapunov func tion , viz. , T m °i V(X,t) = X PX + E g.k.(t) / 0.(z)dz, (1.2.5) i=l 1 1 0 1 where k^(t) is a time varying gain. In [4,5], it Is also proved that the sufficient condition under which the Lyapunov function is- valid is similar to the Popov criterion, dis cussed below. §1.3 The Popov Criterion The sufficient condition of absolute stability for a controllable and observable time invariant system with one nonlinearity satisfying (1.2.3) was established by V.M. Popov [1]. The Popov criterion takes the form Re[(l+qjco) G(jui)] +~ > 0 (1.3.1) for all to, where q is a nonnegative number. A convenient graphical method exists for testing (1.3.1). A sufficient condition for the absolute stability of a system with many nonlinearities was given by Jury and Lee [2]. The condition requires the Hermitian matrix: 2K_1 + H(jo>) + HT(-ju)) (1.3.2) to be positive definite for all OJ. Here K is a constant diagonal matrix of elements K... ,K , . . .K , all of which are positive numbers such that the inequality 1 I m 0 < 0.(o\)o\ < K.o.2 (1.3.3) iii ii is satisfied for i = l,2,...m, and where H(joi) = (I + JOJQ) TO), (1.3.4) unstable region Y 'Popov A 'me Stable region Fig. 1.7 Popov criterion 10 where Q is a constant diagonal matrix with elements q1,q„...q which are 12m nonnegative numbers. This forms a generalization of the Popov result. For the time varying nonlinear feedback system, the sufficient condition of absolute stability established by Rozenvasser is as below, Re G(jco) + ^ > 0 (1.3.5) K for all co. It is but the Popov criterion with q = 0. Other new criteria, similar to the Popov's, were introduced by Zames' and Falb [8,9], Yakubovich [4], Narendra and Taylor [6], Baker and Desoer [10,11], Bergen and Rault [12], and Anderson [13]. §1.4 L2 Stability The concept of the stability has been introduced by Sandberg [14]. It is closely related to asymptotic stability. L„ is the space of square integrable, valued functionjon [t ,°°), it is assumed that L -is-'a linear, inner-product, i/formed space; the inner product of x and y in is <x,y> = r x(t)-y(t)dt '< °°, (1.4.1) ' o and the norm of x is || x || 2 = /<x,x>. Suppose a(x) is in L^fO,00], a(t) is uniformly continuous, and <j(t) is bounded, then the state.o(t) approaches null state if the sufficient condi-tion of absolute stability of a time invariant nonlinear system, the Popov criterion, is satisfied. A further L2 bounded condition was introduced by Zames [8,9]. 11 CHAPTER 2 THE ABSOLUTE STABILITY OF A TIME INVARIANT NONLINEAR FEEDBACK SYSTEM §2.1 Introduction In the previous chapter, it was mentioned that the nonlinear systems can be considered in 4 classes according to the form of the transfer matrix res). Testing absolute stability of a single-loop time invariant nonlinear system using the modified Nyquist diagram was first initiated by Popov, basing the method on his criterion. Further developments using a graphical method to test the absolute stability of a nonlinear system have been furnished by Naumov [15], Meyer and Hsu [16], and Murphy [17]. A graphical method of testing the absolute stability of a time invariant series system with m-nonlinearities and m-identical linear transfer functions was introduced by Davies [18]. ' §2.2 A Graphical Test of the Absolute Stability of a Series System with Nonlinearities and Identical Transfer Functions Consider a series system with the linear time invariant transfer function matrix , 0 . 0v 0 G(s) r(s) = -G(s) 0 o 0 \ \ 0 -G(s>N \ N\ \ 0 \ \ -G(s) 0 (2.2.1) The input, a.(t), and output, 0^(oO> of the i-th nonlinear element satisfy the inequality 0 < a.0.(a.) < a. iii l (2.2.2) 12 Fig. 2.1 Graphical criterion for the time invariant series system with 5 ideal nonlinearities and 5 identical linear transfer functions. Such a system falls within the class for which a Popov-like stability criterion has been, established. In applying this result, the matrix where Q is an arbitrary, semi-positive, diagonal matrix of constants, is con sidered. A sufficient condition to establish the absolute stability of a nonlinear system is that the Hermitian matrix (1.3.2) must be positive definite for all co. If all nonlinear elements are assumed to have the same upper bound; that is, = K for i = l,2,...m, then, without loss of generality, K may be taken as the indentity matrix. If Q=0, the stability criterion is equivalent to the Nyquist plot of G(s) lying within a symmetric m-sided polygon. In the subsequent development it is not required that Q=0, but rather that all elements of Q are equal; that is, Q=ql, where q is a positive scalar constant. where r'(jco) is identical to r(jco) except that G(jco) has been replaced by G'(jco) = (l+jcoq)G(jco) • Thus it is possible to consider the case 0=0 by applying the earlier results for Q=0 to G'(jco) instead of G(jco) itself. If it can be shown that G'(jco) lies within the appropriate polygon, for any positive q, then stability has been established. H(jui) = (I+jcoQ)r(jco), (2.2.3) If Q is restricted in this manner, then (I+jcoQ)T(jco) = r'(jco), (2.2.4) §2.2.1 Main Method-Common Popov Line Let |G'(jco)| = *f and G'(jco) = 6. If G' (jco) lies within a polygon, then y cos (6 - a) < 1, (2.2.5) where a = 2iiT m and i = 1,2,... m, 14 and where i is one of m values each corresponding to the m sides of the polygon. '• Now where Thus and Define G'(jco) = (l+jwq)G(» = (l+jo)q)(R+jI), G(ju) = R(u) + jl(oo) Re G' (jcj) = R - jtol = ycosG , Ira G' (jco) = uqR + I = f sine , (R-coql)cosa + (ioqR+I)sina < 1, (Rcosa+Isina) + qco(-Icosa+Rsina) < 1. -X(co) = Rcosa + Isina, (2.2.6) (2.2.7A) (2.2.7B) (2.2.8A) (2.2.8B) (2.2.9A) (2.2.9B) Y((jj) = co(-Icosa + Rsina). (2.2.8B) gives -X + qY < 1. • (2.2.10) From (2.2.9A) and (2.2.9B), the m different modified Nyquist lcci are plotted each corresponding to one of the m sides of the polygon. To satisfy inequality (2.2.10), all these loci must be to the right side of a straight line, the Popov line, which passes the point (-1,0) having slope 1/q. If such a straight line exists, then the absolute stability of the system is established. Example.2.1 Consider a feedback system of the type shown in Fig. 1.3 (m=3). Let every nonlinearity satisfy the inequality 0 < 0.a. < a., l l l where i = 1,2,3 and let every linear block be represented by (2.2.11) 15 Fig. 2.2 Modified Nyquist loci and the common Popov line for the series system, „, . . 32(s+0.25) • m-3, G(s) = (s+1)(s+2)(s+4r ' f°r ExamPle 2-1' 16 n . , 32K(s+0,25) /, , ,^ l (s+1)(s+2)(s+4) . The m different modified Nyquist loci and common Popov line are shown in Fig. 2.2. Then, by setting q = 0.23, the absolute stability condition K < 0.476 (2.2.13) is obtained. §2.2.2 Simplification in Particular Cases Plotting the m different modified Nyquist loci is tedious if m >_ 4. A simpler and more direct approach is possible if G(s) is of the form G(s) = -p (2.2.14) n (S+D.) j=i 3 or n k n . , (s+N.) G(s) = —V — (2.2.15) P n (S+D.) j=i 2 where p > 1+m, N. and D. are real positive constants, and N. > D. for i = j. i 3 1 J If the transfer function G' (jco) = (l+jioq)G(jo)) having q^ and k^ which are found from the modified Nyquist locus for i = 1, satisfies the above conditions, then ^ (1) The phase of G'(jco) is decreasing as co increases. (2) For some co = co , giving Z.G'(jco ) > and |G'(jco )| < 1, |G'(jco)| c c — m c decreases when co increases for co > co . c And k^ is the maximum value satisfying the absolute stability conditon. 4 Y(u) 17 J -r 2.3 Modified Nyquist locus for i=l, and the corresponding Popov line for the series system, m=3, G(s) = (s+1) (s^s+2Q) ., for Example 2.2.1. 2.4 Nyquist plot of (1+0.788$)G(s), where G(s) = 100 , and a (s+l)(s+5)(s+20) graphical testing the stabilitv 18 a graphical testing the stability of the series system, m=4. for Example 2.2.2. - -19 Example 2.2 (1) Consider m = 3 with every linear element having a transfer function G(s) = 10°k (s+l)(s+5)(s+20). (2.2.16) From Fig. 2*3, q± = 0.788, (2.2.17) and k = 1.46. (2.2.18) Setting q = 0.788, k < 1.46 is the absolute stability condition since, as shown in Fig. 2.4, G'(s) staisfies the above conditions. (2) Using the same approach the absolute stability condition of the system with 4 nonlinearities is obtained as: k < 1.234 (2.2.19) by setting q = 0.766. The information required for the previous method may be obtained directly from the Nyquist Plot of G(s) by noting that the critical point is that having phase a^, and the corresponding value of q is given by tan(a^~3) where tanp is the slope of G(s) at a= a^. If the transfer function G(s) does not satisfy the special form of (2.2.14) or (2.2.15), the Nyquist locus of G'(s) haivng q^ and k^ found from inequality (2.2.10) for i=l may be tested by the polygon criterion. If this fails, the general approach must be adopted. §2.2.3 Use of the Modified Nichols Chart to Obtain 0 The use of the modified Nichols chart to obtain q and to test the absolute stability of the system with many nonlinearities is also possible. A polygon will be described by some relationship between the log amplitude and the phase of the form 20 M = 20 log,rt sec(- ^ + - - 0)db, (2.2.20) 10 mm > 2ii 2n where - (i-1) < 6 < i , m — — • m and i = 1,2,...m. The stable boundary G' can be represented in the modified Nichols chart. Let us consider _1 M = M - 20 log1()(l+N2)2, (2.2.21A) and where j = 0, tan"1 N. , (2.2.21B) J and N_. is an arbitrary, positive constant. (2.2.21A) and (2.2.21B) give a family of the stable boundaries as shown in Fig. 2.7, each corresponding to one of the constants N. Note that these curves are all of the same form and may easily be sketched. If the locus of G(jco) is sketched and q is chosen, such that each point co. on the locus of G(jco) on the modified Nichols chart is beneath the corresponding stable boundary N_. = co^ <J, then the absolute stability of the system is established. Also, if the position of the point co. is known, and is beneath a 2 N. family of stable boundaries N, then the corresponding values of q. = — j = 0,...°°, is known, so that every point has one corresponding stable q range. Then, any q in the interior q range is permitted to be chosen for establishing the stability of the system. This seemingly complex procedure is, in fact, quite straight-forward. If a range of q's is permissible, that value of q is chosen which gives the greatest possible value of k. Example 2.3 Consider the system in the Example 2.2 with every linear element G(s) where 700 Fig. 2.7 Gain-phase plot for. G(s) = -,—TTTTCT—, and a family of boundaries for (s+10)(s+50) the series system, m=3, for Example 2.3. 22 6db Fig. 2.8 630 Gain-phase plot for G(s) = (s+xo)(s+50) and a family of D°undaries for the series system, m=4, for Example 2.3. Fig. 2.9 Gain-phase plot for G(s) - 580 .____ (s+10)(s+50) and a family of boundaries for e series system, m=5, for Example 2.3. 23 G(S) = (s+10)(s+20T (2'2-22) and m = 3. Comparing the frequency response curve of G(juj) with a family of the stable boundaries as in Fig. 2.7, the value of q and stable condition are q = 0.067 (2.2.23A) k £ 700 . (2.2.23BSimilarly, from Fig. 2.8 and Fig. 2.9, by setting q = 0.067, the following absolute stability conditions k £ 630 for m = 4, . (2.2.24) and k £ 580 for m = 5, (2.2.25) are obtained. §2.2.4 An Analogue-Computer Technique This last method makes use of the analogue-computer to test the absolute stability of a nonlinear system. The computer arrangement, which is shown in Fig. 2.10A and Fig. 2.10B, is divided into two main parts. The first generates R and I, and depends on the particular transfer function being considered while the second gives the components of G(jco) which are set to remain unchanged for differing systems. The effect of varying q on the G'(jco) locus is easily obtained by adjusting a potentiometer. This method, however, suffers from difficulties in amplitude scaling. Besides, the analogue-computer set-up becomes more complex with increase in system order. It is also noted that to sketch the complete locus of G'(jco) from co = 0 to co = 00 by the computer is impossible because the describing time o M is proportional to . To improve the accuracy of the output it is necessary to use a three- (or two-) stage programming and rescalin'g technique, and also to limit the to range which is dependent on the computer, characteristics. 2.10B Computer program for X-Y producer Example 2.4 The system (m=3) has linear elements G^(s) where Gi(s) (s+0.5)(s+l) (2.2.26) The circuits of the two parts of the analogue computer set-up are shown in Fig. 2.10A and Fig. 2.10B, and the typical resultant loci for various values of q shown'in Fig. 2.11 lead to the choice q = 1.5. This therefore permits the choice of k = 0.685. max §2.3 Absolute Stability of the Series Nonlinear System with Different Transfer Functions Let us consider the case where the identical linear transfer functions in the previous section have been replaced by different linear transfer functions, Now the transfer matrix r(s) may be written / 0 r(s) = o . o -G-^s) 0 -G2(sJ\ G (s) N, m 0 I I 0 (2.3.1) -G ^(s) 0 m-1 Let us consider H(>) = (I+jwQ)r(ju)), where I is a unit matrix and Q is a diagonal constant matrix with elements q^,q2,...q^; thus H(s) may be considered as below, / 0 H(s) = G'(s)^ m G^(s) 0 -Gl(s) (2.3.2) 28 where G'(s) = (1+jwq.)G.(s), i=l,2,...m. i 11 (2.3.3) Let |G^(jco)I = Y£ and ZG^(JW) = 6i; then the Herraitian matrix 21 + H(jco) + H (-ju>) becomes / 2 -Yle~j6l 0 Y ej6m '1 \ 'm -Y2e 2 rJ69 ~Y0R 2 0 I _j9 ___0 < v e m ^ -j8 , \ -y e m-1 -Y ,e m-1^ 'm-1 2 (2.3.4) and must be positive definite for all a> to satisfy the stability condition. The first (m-1) principal minors of A^ are generated by the recurrence relation A. = 2A. l 1-1 " Yi-1 Ai-2 > 0, i = 3,4,...m-1, with A^ = 2 and = 4 - Yj_ A = m The last condition is ri1 •y.e ^1 0 V91 2 'Y9e 2 v e 2 m m -Y2E 2 e^m 'm Ym-1 \ ie \ 0 —Y _eJ m-1 2 'm-1 ! 6 e-J6m-2 > 0, (2.3.5) that is , 9 9 m m A = 2A ,-(Y , +Y )A „+(-!)( TT Y.)e(E 9.) m m-1 'm-1 m m-2 . .. 'i . , l 1=1 i=l m (-E 6.) m • -i 1 + (-l)m(Tf Y,)e1=1 i=l 2. 2, 2A-,-(Y /+Y")A „+(-Dm( Jl Y.)2 cos( E 0.) > 0. m-1 'm-1 'm m-2 . ., ' l . , l 1=1 1=1 m m (2.3.6) If all y. < 1 f°r aH w, i=l,2,...m, then 1 and A2 =-4-Y^ > 3, 2 2 A3~A2 = A2-Y2A-L > ^2~^1 = - 2 > 1, 2 A3 > A2+1 = 4~Y1+1 > 4: VA3 = VY3A2 > A3"A2 > X' > A3+l > 5, 2 A.-A. , = A. ., —Y. ,i. , > A. .-A. 0 > 1, l l-l l-l 1-2 l-l i-2 A. > A. ,+1 > i+1, l i-l ' 2 A = A -Y „A > A -A > 1, m-1 m-2 m-2 m-2 m-J m-2 m-3 A T > A „+l > m, m-1 m-2 m A > 2A ,-2A „+(-l)m2cos( E 6.) m m-l m-z . . l i=l m > 2+(-l)m2cos( E 6.) >_0. (2.3.8) 1=1 1 . Therefore, the absolute stability of a nonlinear system is assured if every locus of G^(joo) lies within'a circle of unit radius. Obviously, must be chosen zero so that G'(jw) = Gi(joJ) In order to test absolute stability of nonlinear system it is thus necessary to sketch the loci of G^(jio) and to observe whether all of the loci lie within the unit circle. 30 As m -> °°, this result coincides with the previous result in §2.2, but here the G_^(jaj)'s are not necessarily the same. §2.4 Absolute Stability of the Parallel Nonlinear System The transfer matrix of the parallel nonlinear system shown in Fig. 1.4 may be written Define and r(s) = , G^s) G (s) G (s) 1 / m Gl(s) ' j ! i G1(s) G2(s) m (2.4.1) then Suppose that G.(s) = G(s), i=l,2,...m, and G'(joo) = (l+jcoq)G(ju)) , ^G|(jw) G2(jco) G| (ja>) ^G^(jco) G2(ju) G' (ju>) ^ m. ,i I = G'(jco) 1 1 1 ' (2.4.2) G(jco) = R(co) + jl(co) , R' (jio) = R(co) + ql(w) , (2.4.3A) (2.4.3B) I1 (ju) = qtoR(w) + I(w) (2.4.3C) Let us suppose the i-th nonlinearity satisfies 0 < a.0.(a ) < a.. 1 x i I (2.2 The Hermitian matrix is 21 + 2R' ,1 1 1-I I I ^1 1 1/ /2(1+Rf) 2R' 2R' \ I 2R' 2(1+R')\ ' | \ 2R* s 2R' 2R' 2(1+R')' (2.4 The sufficient condition of absolute stability is that the Hermitian matrix must be positive definite for all w; consequently, > 0. Now 2(1+R') 2R' 2R' \ I 2R' 2(1+R'>s I \ \ \ I 2R1- -2R* \ I 2R' \ 2(1+R') 2(1+R') 2R1 2R* -2 2 0, -2 2 \ °N 1 ! \\Si I \ \ • 0 .0 2 32 = 2A. . + 21R' i-l = 2(2Ai_2 + 21_1R') + 21R' = 21~1A1 + (i-l)2V = 21(l+iR') > 0. (2.4.5) Substituting (2.4.3B) in (2.4.5), . R(u>) - q'ajl(u) + j > 0. (2.4.6) Define X(ui) = R(u>), (2.4.7A) and ^ Y(u>) = COI(OJ). (2.4.7B) (2.4.7A) and (2.4.7B) give X(u)) - qY(w) + j > 0. (2.4.8) The condition A = X(u>) - qY(u) + - > 0. (2.4.9) m m implies that A > 0, i = 1,2,...m-1. (2.4.5) Hence,the new Popov line,shown in Fig. 2.12, passes through the point (- ^, 0) with slope —. q 33 2.12 Extension of Popov criterion for the parallel system with m identical nonlinearities and m linear transfer functions 34 Chapter 3 ABSOLUTE STABILITY OF A TIME VARYING FEEDBACK SYSTEM WITH MONOTONIC NONLINEARITIES §3.1 The Absolute Stability of the Single-Loop Time Varying Nonlinear Feedback System In recent years some results concerning the absolute stability of a single-loop nonlinear system with a time varying gain have been obtained by Rozenvasser [5], Zames [22], Bergen and Rault [12]. The results to be presented here extend this previous work. Let us consider single-loop time varying nonlinear system shown in Fig. 3.1. In that system, a (t) = St g(t-r) e(i)dT for t > 0 e o ~~ (3.1.1) = 0 for t < 0 is the zero-state response of the linear time invariant part with transfer function G(jw) =^[g(t)]. The input il(t) represents the zero-input response of g(t). The complete response of g(t) is thus c(t) = o (t) + n(t). (3.1.2) e The input, a(t) = -c(t), and the output 0(c) of the amnesic nonlinearity N.L. are related in the following manner: (1) 0 < 00(a) < a2 for a i 0, 0(0) = 0, 0(a) - 0(a) (2) 0 < a _ a — — 1 for CTx * V The block, k(t), represents a linear time varying gain, thus e(t) = k(t)0[a(t)]. The instantaneous value of this gain is constrained so that (1) K± < k(t) < K2, where K2 > K± > 0, (2) bk(t)< k(t) < ak(t), where the number a> 0,and the number b (3.1.3) is finite. It is assumed that the linear part is stable, more specifically, (1) g(t) E L2(0,«0, g(t) e L1(0,K>)j (2) n(t) e ^(O.co), (3) n(t) is differentiable and f)(t) e L (0,»). Condition (1) above ensures that g(t) is bounded on (0,°°) and that g(t) 0 as t •+ °°; besides, conditions (2) and (3) ensure that n(t) behaves in the same manner. Denote n(t) = sup | n(t)|, t>0 g(t) = Sup |g(t)|. U t>0 The Fourier transforms of g(t), e(t), etc., are denoted by G(jco) E(jco), etc. The notation f| • II denotes norms in the space L^(0,°°). Thus ||n(t)|l = /~ |n(t)|dt. §3.2 • The Main Result The main result is the following theorem. §3.2.1 Theorem 1 Consider the system shown in Fig. 3.1 to which the assumptions made above apply. Let y(t) be any real function such that (1) y(t) =0 for t < 0, (2) y(t) <0 ' for t>0, Kl (3) ||y(t)|| < . K2 and let q be any nonnegative number. If Re( [l+qjw+Y(jco)] [G(jto)-4-]+aqG(jto')} K2 -lly(t)ll (f -f'-0 (Q1) Kl- K2 36 for all to, then (1) Sun |a(t) | < «, t>o (2) a(t) -v 0 as t <», (3) as |ln(t)|| + U n( t) 11 ->- 0, the corresponding a(t) has the property that Sup |o(t) J ->- 0. f >0 §3.2.2 A Special Case of the Theorem It should be noted that if the time varying gain k(t) is monotonically non-increasing and a <_ 0 V t >^ 0, then the condition (01) for absolute stability may be replaced by Re[l+qj<^Y(juO][G(ju))+^H - || y(t) II -K2 Kl K2 3 > 0 CQ1') for all a), 2K2 where again, q is any nonnegative real number. §3.3 Proof of Main Result The body of the proof of Theorem 1 will be given in a series of appendices; a brief summary is given below in this section. Define a(t) for t < T aT(- ' 0 for t > T, ,(t) ={ and Thus eT(t) = k(t)0[aT(t)], aeT(t) = g(t-x) eT(x)dT. (3.3.1) a T(t) = a (t) for t < T, eT e — and aeT(t) e L1(T,-) 37 Define The notation (x*y)(t) denotes convolution between x(t) and y(t); (x*y)(t) = x(T)y(t-T)dT. . (3.3.2) a - o + o*y, c = c + c*y. m- • J m J Then T e (t) /' [aAt) - -~—-]e(t)dt U m K.^. •= /J[a(t) - ^-]e(t)dt + /J[y*(a - |-)(t)e(t)dt. (3.3.3) Define k(t-x)0[aT(t-T.)] R(x) •= /J[aT(t-T) ^—- ]k(t)0[aT(t)]dt. (3.3.4) Now where and R(x) = R1(T) + R2(T.), (3.3.5) R1(t) = f0 toT(t-T)-0[aT(t-T)]]k(t)0[aT(t)]dt, (3.3.6) R2(T) = /Q[1 - k^ T) ]k(t)0[aT(t-x)]0[aT(t)]dt. (3.3.7) Frem Appendix 1. ro K k(t)0[a (t)] R1(T) l/0 K [aT(t) K ]k(t)0[aT(t)]dt, (3.3.8) and from Appendix 2, oo 1 1 2 R2(T) <_ /0(^- - K~)eT(t) dt. (3.3.9) Thus bo K2 eT(t) R(T) ±VK aT(t) ~ ~i ]ex(t)dt- (3.3.10) The first term of the right side of (3.3.3) is always positive. Let us now consider the second term of the right side of (3.3.3) ;0ty*(a " K^](T)DT 38 = /J /Jy(T)[o-(t-T) - e(^"T)]e(t)dT dt co co e (t-T) = fQ /Qy(x)[oT(t-T) - R ]eT(t)dT dt = /0y(T)R(x)dx. (3.3.11) But, employing (3.3. 10), it may be shown that fl[y*(o - f-)](t)e(t)dt 0 > /gyCOdx aT(t) - -|^-]eT(t)dt (3.3.12) - ' llyll/0[K^ °T(t) " K^~]eT(t)dt- (3.3.13) Substitution of (3.3.13) in (3,3.3) yields e (t) /S[0»(t) - -YT-^t + W " y" ;oeT(t)2dt S T 6T(T) > (1 - -^llyll/i[aT(t) - -| ]e_(t)dt (3.3.14) Kl 0 T K2 T and, from the assumptions of the Theorem 1, the right side of (3.3.14) is non-negative. Hence, e (t) - -f2-]e(t)dt + lly" 'SeT(t)2dt - °- (3-3-15) Consider the following integral . T EM(T) 1 1 I = /o[_CTem(t) ~ q°e(t) " ~T~ aqae(t) + (f- - £-) ||y|l e(t)]e(t)dt = I1 + I2, (3.3.16) where and T e (t) i l I = /'[-a (t) - + (f- - f-) ||yll e(t)]e(t)dt, (3.3.17) ± u em i^-^. 1^2 I0 = fl - q[a (t) + aa (t)]e(t)dt. (3.3.18) 2 0 e e Now 39 Ix = V°m(t) + ^-^)Hylle(t)]e(t)dt + nm(t)e(t)dt (3.3.19) giving, due to (3.3.15) I > i'l n (t)e(t)dt. (3.3.20) 1 — 0 m In a similar manner, I2 = /Qq[a(t)+aa(t)]e(t)dt+/Jq[n(t)+an(t)]e(t)dt. (3.3.21) Invoking the result of Appendix 3, I2 L q[k(T)$(T)-k(0)$(0)+/Qq[fi(t)+an(t)']e(t)dt. (3.3.22) Recall that /Jf(t)dt = /0 fT(t)dt, where f^(t) is the truncated version of f(t) to between 0 and T. Hence from (3.3.16), 1 = V^eTm^ " q6eT(t) " " qaCTeT(t) + (~ - jp) ||y I! eT(t)]eT(t)dt (3.3.23) and, from the conditions of Theorem 1 and since 1 CO 1 1 = " 27 L,Re{[l+qjco+Y(jw)][G(juO + + aqG(jw) -(^- " ~) llyll }ET(ju3)E*(jW)du), it follows that I 1 0. (3.3.24) Since I = I + I?, from Appendix 3 and after substituting (3.3.20) in (3.3.24), • q[k(T)*(T)-k(0)*(0)] < - /*[aqn(t)+n (t)+qn(t)]e(t)dt. (3.3.25) — Urn Define e = Sup |e(t)1 0 < t < T 40 and invoking the conditions already imposed upon n(t), fi(t) and y(t), the right side of (3.3.25) must be less than the quantity eM[(l+ llyll +aq) ||n|| +q llfill ] = MeM. (3.3.26) Furthermore, considering the first term of the'left side of (3.3.25), since 0 is monotonic, *(t) > |{0[a(t)]}2. (3.3.27) Using (3.3.26) and (3.3.27) in (3.3.25) yields ' | k(T){0[a(T)]}2 <_ qk(0)<D(0) + Me^. ' (3.3.28) From (3.3.28), since < 1, 2 ~ •L- e(T)2 £ qk(0)*(0) + MeM. (3.3.29) The inequality (3.3.29) holds .for any T >_0 and implies KM KM Sup |e(t)| ± + [(-~-)2 + 2K2k(0)$(0)]1/2. (3.3.30) Furthermore, since a (0) = 0 by (3.3.1), this bound on e(t) tends to zero with II hli + llfill • It remains to be shown that |e(t)| + 0 as t ->• °°. After substituting in (3.3.14) of(A 3.3) of Appendix 3 and using . (3.3.24) , K (1 - -^llyll )/J[a(t) - ^-]e(t)dt+/Qnm(t)e(t)dt -qk(0)$(0) + q/J[fi(t)+an(t)]e(t)dt <. 0. (3.3.31) Thus K (1 - r^llyll )fl[a(t) - ^-]e(t)dt < qk(0)*(0) .+ Me . (3.3.32) K^ 0 K2 — M From which it follows that T e(t) qk(0)*(0) + Me fQ[o(t) - ^-]e(t)dt i ^ . (3.3.33) 2 i-/llyll Kl 41 Since right side of (3.3.33) is independent of T, then letting l.-> «>, e(t) qk(0)*(0) + Me fAo(t) - ^-]e(t)dt < M 0LUVW K_ JC^'UI" - K , (3.3.34) 2 llyll 1 and this bound on the integral tends to zero with || nil + llnll . However, the bound already placed on e together with the conditions demanded of g(t) require 0" (t) to be bounded and to tend to zero. It may now be shown that the integral of (3.3.34) is infinite unless a (t) + 0 as t co, thus contradicting (3.3.34). Therefore it can be concluded that a(t) -> 0 as t 00. Since a (t) is uniformly continuous and if Sup |a(t)| does not go to zero as |lnll + llnll 0, then (3.3.34) does not tend to zero either. This is a contradiction. Hence, as llnll + II n II + 0, Sup |a(t)| -> 0. The proof is t>0 therefore completed. §3.4 Examples Example 3.1 In the single-loop nonlinear system with a linear time varying gain shown in Fig. 3.1, the linear part has a transfer function G(s) = (s+l)(s+2) ' the input and output of the nonlinear part, satisfy (3.1.2), and the time varying gain is such that (1) 1 <_ k(t) <_ 1.2, i.e. , K = 1, K = 1.2, (2) bk(t) <_ k(t) <_ ak(t) , where a,b are real numbers such that a > 0, and b is finite. Suppose that a is large enough, q must be chosen zero. Let us assume that y(t) is an exponential function such that 42 rct)=o p" o N.L. eft) Gcs) i|Ct7 c<t) -f 3.1 Single loop time varying nonlinear feedback system with zero-input 43 e ' as t > 0 yM-{ T 0 as t > 0, 2 where y > 0, and g > 1.2 = —. Then 1 lly(t)H =rQ l-ie'^l dt . Thus the conditions on y(t) are satisfied. After taking the Laplace transform of y(t), 1 + Y(s) = 1 - - 1 = + TS + Y S +1 3 + YS . (3.4.1) • Y The sufficient condition of absolute stability is Re{[Cp-l) + Y.V][G(ja)) +•_! ] _ I ^| o (3.4>2) 3 + YJW !•2 3 1-2 — for all 03. Let us define G1(Ja)) = GCjo,) + and From the relation (3.4.2), the locus of G|(joi) must lie on the right side of the vertical line passing through the point (r~L^~~ > 0) • The multiplier of 3 G^(jw) in G^(jo)) may be considered as a compensator which is shown on the Bode diagram Fig. 3.2. The function of the compensator is to improve the characteristics of G^(jo3). From the plot of G^(jw) in Fig. 3.3, it is a simple matter to choose the proper values of g and y. From the fact that the left-most point OJ^ = 2.5 • 3—1 3 on the locus of G.. (JOJ) lies between and — and that g satisfied the relation 1 y y li-l 1 -j 1 ., . . ft"Urz _ ir~J <<c the stability condition 3 Kx K2 K2 K < 23.75 (3.4.3) 44 G= (UJOJXUJCO) Kf =23. 75 Fig. 3.3 The plots of G(s) = (3+i)(g+2) and G'<s) = C^^J-)GCs) for Example 3.1. 45 has been found by running a suitably written programme on the digital computer. It is noted that the combination 3=2 and y = is not the best one because of the particular choice of y(t). That optimum y(t) which gives the best combination of 3 and y maY be determined by a digital computer technique. Example 3.2 Consider G(s) K (s+1)(s+D)' where D is any positive real constant, instead of the-G(s) in the previous problem and let 0 < k(t) <_ 1, . and bk(t) <_ R(t) <_ ak(t) . In this case, let us suppose y(t) = 0, then the' condition of absolute stability is Re{(l+qjio) [G(jco)+l]+aqG(jco)} >_ 0 C3 4.4) for all co. Rewritting, R(co) - qcol(co) + 1 + aqR(oj) ^ 0 (3.4.5) ' for all co, since G(jco) = R(co) + jl(co) . (3.4.6) Here „ /n 2. K^CD-co ) and Define -K (1+D) oi = (l+co2)(D2+.2) (3-4'7B) X(co) = R(co) (3.4.8A) and Y(OJ) = uiy - aR(oj). (3.4.8B) Substituting (3.4.8A) and (3.4.8B) in (3.4.5), X(OJ) - qY(a)) +1 >0. (3.4.9) To satisfy inequality (3.4.9), the locus must be on the right side of the straight line passing through the point (-1,0) having positive slope —. q From (3.4.8A), (.3.4. 8B). and (3.4.9), if a <_ l+D-e, where the small number e > 0 is arbitrarily chosen, and q is then the stability condition (Ql) is satisfied for any nonnegative real constant k. Besides, the modified Nyquist plot is on the right side of the Popov line. This is shown in Fig. 3.4. If a ^ 0, the sufficient condition of absolute stability is R(OJ) - qoKu) + 1 >_ 0. (3.4.10) This is the Popov criterion and is satisfied for any nonnegative real constant k 1 by choosing q = 1+D" §3.5 Absolute Stability of a System with Many Nonlinearities and Many Time Varying Gains In the previous secLions, the absoluLe stability of the system with one nonlinearity and one time varying gain is established. Now, let us consider i system with many nonlinearities and many time varying gains. Such a system is shown in Fig. 3.5. The input a. and the output 0.(o.) of the i-th nonlinearity l ii are related by the following: (1) 0 < a.0.(a.) < a2 for a. i 0, 0.(0) = 0, lii—i l l d0.(a.) (3.5.1) <2> 0<-<t^ l1-i and the instantaneous value of the i-th time varying gain is constrained so that (1) K, . < k.(t) < K. . , where K„. > K_ . > 0, li — I — zi 2i li (2) b.k(t) < k.(t) < a.k.(t), where the number a > 0 and the number 1111 i b. is finite, l 47 Fig. 3.5 General time varying nonlinear feedback system with zero-input 48 Besides, ° = fn r(t-T)e(x)dT e (J — is the zero state response of the linear time invariant transfer matrix r(jw) = J [F(t) ] , where /G11(JOJ) G12(jaj)-G21(ju) G (jo)—• r(ju) = Glm(^^ (3.5.2) I I kGml(U) G^Cjo.) The input vector n_(t) represents the zero-input response of r(t). The complete response of T(t) is thus c(t) = o^t) + _n(t). It is assumed that.all elements of the linear transfer matrix r(t) are stable, more specifically, (1) g±j(t) e L2(0,-),. g (t) £^(0,00), i=l,2,...m, j = l,2,...m, (2) n.(t) e ^(O.oo), (3) ni(t) is differentiable and f^Ct) c L (0,«). Condition (1) above ensures that each element of r(t) is bounded on (0,<») and that g_^.—»nj and conditions (2) and (3) ensure that ^(t) behaves in the same manner as g (t). ij Denote n = E n. = E Sup |n. (t)| i=l i=l t>0 m m m m E E Sup |g (t)| . r = E E g.... . . 1=1 j=l J i=l j=l t>0 J The notation 11-11 denotes norms in the space L (0,c°) such that, for example, 49 m H II = 'n E 'ni(t) ldt i=l Define and K^"*" = diagonal matrix O^-^l 'K] 2 '' ' * "Klm^ = diagonal matrix (K21'K22'"''*K2m^' A = diagonal matrix (a^,,....a^) §3.6 Theorem 2 . • Consider the system shown in Fig. 3.5 to which assumptions given above apply. Let Y(t) = diag{y (t), y (t) . . . .y (t) } be such that each element is a real function and that (1) y±(t)' = 0 for t < 0, i = 1,2,...m, (2) y.(t) < 0 for t > 0, l — — K (3) l|y.(t)ll < i KN Let Q be any positive semi-definite constant diagonal matrix. If there exists anmxm matrix H(jco) such that H(jco) = [I+JWQ+Y(jaj)] [r(jco)+K21] + AQTCjoj) - II YIt [K"1 - K"1], (02) and then T (1) H(jco) + H (-jco) is a positive semi-definite Hermitian matrix for all co, (2) H*(jw) = H(-jco) (3) Every element' of H(jco) is analytic for all co, (1) Sup | o.(t) | < «=, and Sup|a.(t)| < i=l,2,...m, t>0 t>0 1 50 (3) (2) o;(t) -> 0 as t -> 00, and o\ -> 0 as t -> oo, as Hull + Ilil II _> 0> tne corresponding _g_ has the property Sup|o_(t)| -K), and Sup|CT.(t)| ->0. t>0 t>0 1 §3.7 Proof of Theorem 2 This proof follows the same vein as that of the previous Theorem. The only difference here is that all vectors such as e(t), a (t) , n(t), n(t), — —e — — jc(t) , 0^.t) , cr(t) are m-vectors, and all matrices such as 0,Y are mxm matrices. All formulae and the proofs of which developed in §3.3 still hold, except that the proof of (3.3.24) must be performed in the following manner. where h(t) is the inverse Fourier transform of H(jco). From the condition of Theorem 2 and Newcomb's result [25], From (3.3.23), (3.3.23') = -/^eJ(t)/Qh(t-T)eT(T)dx dt, I' < 0. (3.3.24') This, however, is the same as (3.3.24). Hence, by the same argument used in the proof of Theorem 1, (1) Sup|a(t)| < t>0 CO (2) a(t) +0 as t -> CO Now, Sup|o.(t)| t>0 < 00 if and only if every component a.(t) satisfies Sup 1°" (t) I < 00. t>0 i Similarly, cr(t) -> 0 as t -> ~ if and only if every con\ponent ^(t) satisfies rj^(t) -> 0 as t -> co, and Sup|c[(t)| -> 0 if and only if every component a. (t) satisfies t>0 Sup\o. (t) I -> 0. t>0 1 §3.8 Example Example 3.3 Let us consider the parallel system, where m = 3, each branch of which has one nonlinearity and one time varying gain in series with one linear time invariant transfer function. Here these three parts of each the three branches are identical to the corresponding ones used in Example 3.2. Suppose that the matrix Y(t) = 0. Let us consider the matrix H(jw) = (I+ju)Q)(r(ju))+I) + a Qr(joi), (3.8.1) where r(jco) = G(ju) 1 1 IN 111 111'' and Q = ql-Rewritting, H(jco) = (l+jcoq+aq)G(jco) 111 + U+jcoq)!. Ill 111 Obviously, H*(jco) = H(-jco), and the elements of H(jco) are all analytic for all 60. -Invoking the proof of §2.4, the Hermatian matrix (3.8.2) 52 H(ju) + HT(-jw) = Re(l+ju)q+aq)G(joj) + -j. (3.8.3) Following the argument used in Example 3.2, if a <_ 1+D-E. and q is Y> the sufficient condition Q2 of absolute stability is satisfied for any positive real constant K. 53 Chapter 4 CONCLUSIONS A graphical method using the Popov line is possible for a particular class of time invariant nonlinear system. The method may be simplified in a number of cases. Two alternative approaches, one using the Nichols chart, the other the analogue computer, are mentioned briefly and illustrated. No simple graphical method exists to test the absolute stability of the parallel system with many different linear transfer functions, although a graphical method using the Popov line to test the absolute stability of the parallel system with many identical nonlinear transfer functions is possible. Neither is there any simple graphical method available to establish the criterion of absolute stability of a multi-circuit or an internal feedback system. However, work on the determination of the criterion of absolute stability for any one of the four classes mentioned in §1.1 by digital technique is underway. The digital technique [23] is in essence concerned with location of the optimum combination of matrices Q and K which will define the ^boundary of absolute stability region. It must be pointed out, however, that with the systems that have been discussed so far in this thesis the graphical method is so far simpler and less cumbersome in obtaining the requisite conditions for absolute stability. In chapter 3, Theorems 1 and 2 provide the sufficient, but not necessary, conditions for the absolute stability of a time varying nonlinear system in which the nonlinear part must be monotonically nonlinear. Of great importance in establishing the sufficient condition of absolute stability is the appropriate choice of A and y(t). 1 + Y(s) is identical to the function describing some RC passive network if y(t) is an exponential function. The optimum region of absolute stability may be found by a digital technique. 54 If the time varying gain is frozen, that is,k(t) = 1, the sufficient condition. (Ql) may be rewritten Re[l+qjw+Y(jui)] [G(jto) > 0. 2 This is the result of Baker and Desoer's [11]. For some classes of time varying nonlinear systems, it is possible to use the- graphical method discussed in chapter 2. APPENDIX 1 From (3.3.6) R^O-R^T) = /o([aT(t)-0[aT(t)]] - [aT(t-T)-0[aT(t-T) 1 ] } 0[aT(t)]k(t)dt Noting that 0(t) and [a(t)-0(t)] are monotonic, and 0 < < 1 "r°2 ' or or (tf1-a2)(01-02) - (01"02)2 > 0, [(a1-01)-(a2-02) 1(0^02) > 0, thus 0(t) is monotonic increasing in [a(t) - 0(t)]. Let us define A= ([cT(t)-0[aT(t)]-[aT(t-T)-0[aT(t-x)])0[aT(t)], and observe that A >_P(t) - P(t-r), where aT(t) P(t) =rQ 0[aT(x)]d[aT(T)-0[aT(x)]]. From (Al.l), (A1.4) and (A1.5), 00 R1(0)-R1(x) >_ /Q[P(t)-P(t-T)]k(t)dt and But R1(T) 1^(0) + /0[k(t+T)-k(t)]P(t)dt. V0) -C[0l(t) - ^taT(t)]]k(t)0[aT(t)]dt, and,since, P(t) < [oT(t)-0[aT(t)]]0[aT(t)l < taT(t) - -1 0taT(t)]] ^ 0[aT(t)], then /~[k(t+x)-k(t)]P(t)dt <_ (Z2~\) /QP(t)dt K —K J"o[aT(t) - ^ 0[aT(t)]]k(t)0[CTT(t)]dt. Substituting (A1.9) and (ALIO) in (A1.8) yields 57 APPENDIX 2 From (3.3.7) R2(T) = /Q[1 - ^^-]k(t)0[aT(t-T)]0[aT(t)]dt. (A2.1) Therefore, R2(T) - Qk^O ~ K^l l^-tW^tt-!)]] |k(t)0[aT(t)] |dt - " K"' lk(t-r)0[aT(t-T)] I |k(t)0[aT(t)] |dt (A2.2) giving r2(t) -I(K" " ^HyK(T~T)0[AT(T~T)]2)DT + /Q[k(t)0[aT(t)]2Jdt} (A2.3)' from which (3.3.9) follows. 58 clearly^ APPENDIX 3 Since ak(t) > k(t), ak(t) > 0, and a(t)0[a(t)] >_faQ(t) 0[a]da = *(t) > 0, /Jak(t)a(t)0[a(t)]dt >/Jk(t)$(t)dt, (A3.1) J = q/J(d(t)aa(t))k(t)0(t) >q [k(t)4>(t)+k(t)0[a(t)] d(t)]dt - q fl IF tk(t)*(t)]dt >_ q[k(T)$(T) - k(0)*(0)] . (A3.2) Therefore, I2 = i + q /J[f|(t)+an(t)]e(t)dt ^ q[k(T)$(T)-k(0)$(0) + /^[f|(t)+an(t)]e(t)dt. (A3.3) 59 APPENDIX 4 2 Since k(t) <_ ak(t) <_ 0 and 0 <_ M^i£(t)J_ ,< $-(t) f clearly, fl ^ e(t)2dt = /J ak(t) • ^t)0[a(t)]2 ^ ,T Ut)$(t)dtj (M<1) J' = q/J[6(t) + ~- e(t)]e(t)dt >_ q/J [k$ (t)+k (t) 0 [a (t) ]d (t) ] dt >_ q[k(T)$(T) - k(0)$(0)] . (A4.2) Therefore, 1\ £q[k(T)$(T) - k(0)$(0) + fi(t)e(t)dt]. ' (A4.3) 60 REFERENCES 1. V.M. Popov, "The Absolute Stability of Nonlinear Automatic-Control System", Avtomatika i Telemekhanika, Vol. 22, No. 8, pp. 857-875, August, 1961. 2. E.I. Jury and B.W. Lee, "The Absolute Stability of Systems with Many Nonlinearities", Avtomatika i Telemekhanika, Vol. 26, No. 6, pp. 945-965, July 1965. 3. V.A. Yakubovich, "The Solution of Some Matrix Inequalities Occurring in the Theory of Automatic Control", Dokl. AN SSSR, Vol. 143, No. 6, 1962. 4. V..A. Yakubovich, "The Matrix-Inequality Method in the Theory of the Stability of Nonlinear Control Systems", Avtomatika i Telemekhanika, Vol. 25, No. 7, pp. 1017-1029, July, 1964. 5. E.N. Rozenvasser, "The Absolute Stability of Nonlinear Systems", Avtomatika i Telemekhanika, Vol. 24, No. 3, pp. 304-313, March, 1963. 6. Kumpati S. Narendra and James H. Taylor, "Lyapunov Functions for Non linear Time-Varying Systems", Information and Control, Vol. 12, pp. 378-393, 1968. 7. Kumpati S. Narendra and Yo-sung Cho, "Stability Analysis of Nonlinear and Time-Varying Discrete Systems", SIAM J. Control, Vol. 6, No. 4, pp. 625-646, 1968. 8. G. Zames, "On the Input-Output Stability of Time Varying Nonlinear Feedback Systems", Part- I, II, IEEE Trans. Automatic Control Conference, Troy, New York, 1965, pp. 735-747. 9. G. Zames and P.L. Falb, "Stability Conditions for Systems with Monotone and Slope-Restricted Nonlinearities", SIAM J. Control, Vol. 6, No. 1, pp. 89-108, 1968. 10. C.A. Desoer, "A Generalization of the Popov Criterion", IEEE Trans. Auto matic Control, Vol. AC-10, pp. 182-18 , 1965. 11. R.A. Baker and C.A. Desoer, "Asympotic Stability in the Large of a Class of Single-Loop Feedback Systems", SIAM J. Control, Vol. 6, No. 1, pp. 1-8, 1968. 12. A.R. Bergen and A.J. Rault, "Absolute Input-Output Stability of Feedback Systems with a Single Time-Varying Gain", Journal of the Franklin Institute, Vol. 286, No. 4, October, 1968. 13'. B.D.O. Anderson, "Stability of Distributed-Parameter Dynamical Systems with Multiple Nonlinearities", Int. J. Control, Vol. 3, No. 6, pp. 535-540, 1966. 14. I.W. Sandberg, "On the L-Boundedness of Solutions of Nonlinear Functional Equations", BSIJ, Vol. 143, No. 4, 1964. 61 15. B.N. Naumov, "An Investigation of Absolute Stability of the Equilibrium State in Nonlinear Automatic Control Systems by Means of Logarithmic. Frequency Characteristics", Avtomatika i Telemekhanika, Vol. 26, No. 4, pp. 591-600, April, 1965. 16. A.U. Meyer and J.C. Hsu, "The Evaluation of Frequency Response Stability Criteria for Nonlinear System via Logarithmic Gain-Phase Plots", Proc. 1968 JACC Conf. pp. 791-811. 17. Gordon J. Murphy, "A Frequency-Domain Stability Chart for Nonlinear Feed back Systems", IEEE Trans. Automatic Control Vol. AC-12, pp. 740-743, December, 1967. 18. M.S. Davies., "Stability of a Class of Nonlinear Systems", Electronic Letters, Vol. 4, No. 5, pp. 322, July, 1968. 19. T.L. Chang and M.S. Davies, "Stability of a Class of Nonlinear Systems", Int. J. Control, Vol. 10, No. 2, 1970. 20. N.M. Trukhan, "Single-Loop Systems which are Absolute Stable in the Hurwitz Sector", Avtomatika i Telemekhanika, Vol. 27, No. 11, pp. 5-8, November, 1966. 21. N.G. Meadows, "New Analog-Computer Technique for Automatic Frequency-Response Locus Plotting", Proc. IEE, Vol. 114, No..2, December, 1967. 22. G. Zames, "Nonlinear Time-Varying Feedback Systems-Conditions for Loo-Boundedness Derived Using Conic Operators on Exponentially Weighted Spaces", Proc. 1965 Allerton Conf. pp. 460-471. 23. A.G. Dewey, "Frequency Domain Stability Criteria for Nonlinear Multi-Variable Systems", Int. J. Control, Vol. 5, No. 1, pp. 77-84, 1967. 24. T.L. Chang and M.S. Davies, "The Absolute Stability of Nonlinear Systems Involving a Time-Varying Gain", Proc. 1969 Allerton Conf. pp. 721-729. 25. R.W. Newcomb, "Linear Multiport Synthesis", McGraw-Hill, New York, 1966.
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- The absolute stability of nonlinear systems
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
The absolute stability of nonlinear systems Chang, Te-Lung 1970-12-31
pdf
Page Metadata
Item Metadata
Title | The absolute stability of nonlinear systems |
Creator |
Chang, Te-Lung |
Publisher | University of British Columbia |
Date | 1970 |
Date Issued | 2011-05-25T23:10:50Z |
Description | This thesis is in two parts, both considering the absolute stability of nonlinear systems. In the first two chapters the stability of certain classes of nonlinear time invariant systems involving several nonlinearities is considered. A number of graphical methods are given for testing the stability of these systems. The graphical tests are equivalent to a weakened form of the Popov criterion. The third chapter derives a stability condition for nonlinear systems involving a linear time-varying gain. The time-varying gain is assumed to satisfy conditions on its magnitude and rate of change. |
Subject |
Nonlinear mechanics Control theory Stability |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Collection |
Retrospective Theses and Dissertations, 1919-2007 |
Series | UBC Retrospective Theses Digitization Project |
Date Available | 2011-05-25 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0102180 |
URI | http://hdl.handle.net/2429/34874 |
Degree |
Master of Applied Science - MASc |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Electrical and Computer Engineering, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-UBC_1970_A7 C43.pdf [ 2.85MB ]
- Metadata
- JSON: 831-1.0102180.json
- JSON-LD: 831-1.0102180-ld.json
- RDF/XML (Pretty): 831-1.0102180-rdf.xml
- RDF/JSON: 831-1.0102180-rdf.json
- Turtle: 831-1.0102180-turtle.txt
- N-Triples: 831-1.0102180-rdf-ntriples.txt
- Original Record: 831-1.0102180-source.json
- Full Text
- 831-1.0102180-fulltext.txt
- Citation
- 831-1.0102180.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0102180/manifest