Abstract. Linear time-invariant discrete digital plant, with no restrictions on the form of the state equation and the number of controls, is consider...

12 downloads 6 Views 209KB Size

DISCRETE-TIME VARIABLE STRUCTURE CONTROL SYSTEMS - MULTIVARIABLE LINEAR PLANT CASE UDC 62-52

Zoran M. Bučevac Mechanical Engineering Department, University of Belgrade

Abstract. Linear time-invariant discrete digital plant, with no restrictions on the form of the state equation and the number of controls, is considered. The stabilizing state feedback control algorithm is developed by Lyapunov's second method leading to the variable structure system with sliding modes. Essentially the control algorithm drives the system from an arbitrary initial state to a prescribed so-called sliding subspace S, in finite time. Inside the sliding subspace S the system is switched to the sliding mode regime and stay in it for ever. Key words: Discrete-time digital system, variable structure systems, Lyapunov's second method, sliding mode

1. INTRODUCTION In the sixties, S.V. Emel'yanov originated the variable structure control systems (VSCS) theory, although this type of ideas could have been recognized in the early works of discontinuous systems in the West in the fifties; see, for example, Flüge-Lotz (1953). Enormous number of papers appeared as a result of early research in the area of VSCS. The results were mostly gathered in two monographs (Emel'yanov, 1967; Emel'yanov et. al., 1970). In Utkin's research in this area we meet more modern and general approach. He oriented himself especially to the VSCS with sliding modes (Utkin, 1974; Utkin, 1981) where did generalization observing no restrictions on the form of the state equations and the number of controls. Research was continued and many papers on the subject have appeared (see, for example, Young, 1977 and 1978; Itkis, 1983; Ryan, 1983; Slotine and Received November 20, 1996; in revised form October 20, 1998 This work was supported in part by the National Science Foundation under project No.11E08PT1

984

Z.M. BUČEVAC

Sastry, 1983; Ambrosino et al., 1984; Gough et al., 1984; Khurana et al., 1986; Geering and Guzzella, 1986; Bartolini and Zolezzi, 1986; Slotine et al., 1986; Xing-Huo and Chun-Bo, 1987; Li et al., 1994; El-Khazali et al., 1995; Karakasoglu et al., 1995; Zheung et al., 1995 and their bibliographies). Almost all published in the area is related to the VSCS with sliding modes where the main feature of the systems is existence of their so - called sliding mode regime inside the so - called sliding subspace. The popularity of the systems stems from their well known advantages compared to the other types of VSCS or other types of control systems at all: 1. In sliding mode regime, the systems' behavior is independent of plant parameters' variations, in some cases, 2. Realizability of desired dinamical properties of the systems in their sliding mode regime, 3. Realization simplicity of compensator of the systems, 4. Sliding mode regime invariance to disturbances, in some cases. The first motivation for the paper to be related to VSCS with sliding modes goes from the all previously mentioned advantages of the systems. The second and the most important motivation for this work stems from the fact that results to date mostly concern continuous time type of these systems. There were small number of papers related to the discrete time version of the problem (see, Itkis, 1976; Milosavljevic, 1985; Sarpturk, 1987; Kotta, 1989; Furuta, 1990; Sira-Ramirez, 1991; Bartolini et al., 1995; Gao et al., 1995). The author of this paper, treated discrete digital type of the problem (see, Bucevac, 1985, 1988), clearly to make place of microprocessor compensator application in the systems. By this paper the author continues the series of his discrete digital type of the problem papers. In Itkis (1976), the author considers amplitude as well as width modulation of a scalar control applied to a continuous linear plant, in companion canonical form, of a VSS. Since the resulting control is piecewise constant, discrete nature of the special case system is evident. Also, in the work of Milosavljevic (1985), the linear plants in companion canonical form with scalar control were considered. He observed all combinations of discretization of compensator input variables. Sarpturk considered discrete equivalent of linear time invariant continuous system and gave necessary and sufficient convergence and sliding conditions. Kotta claimed that Sarpturk's conditions led to the upper and lower bounded control. Furuta treated the scalar control case and used Lyapunov function which was analogous to the one used for continuous type of the problem. Sira-Ramirez treated single input single output nonlinear time discrete system. His approach led to chattering around sliding manifold. Bartolini et al. gave control algorithm guaranteeing finite reaching time to the sliding manifold as well as taking in account boundness in value. Also, the "equivalent control" was introduced providing motion in sliding manifold. Gao treated scalar control case of time discrete VSCS. Bucevac discovered phenomenon of discrete digital VSCS: during sliding mode regime inside so-called sliding subspace actual so-called sliding control must be applied to the system. Also, one step control was discovered. In the work linear time invariant discrete digital system is under consideration. A state feedback discrete-time variable structure control algorithm is designed. The control algorithm is stabilizing but through the art of VS with sliding modes. Lyapunov's second method, as an elegant technique, is used. This alines the work among the small number of papers where the method was primarily used for analysis or design of VSCS, see, Gough

Discrete-Time Variable Structure Control Systems Multivariable Linear Plant Case

985

et al. (1984), Xing-Huo and Chun-Bo (1987). Essentially, the objective is to push the system state from an arbitrary initial position to the so-called sliding subspace S, in finite time. Once the system reaches the sliding subspace S it stays there for ever working in socalled sliding mode regime. The main feature of the regime is the system state assymptotical approaching to the zero equilibrium state. The solution of the problem makes place for making of commercial microprocessor compensator of the type. So far the problem related to the most general case of multivariable linear time invariant plant was not solved by Lyapunov's approach as proposed in the work. 2. PROBLEM STATEMENT, NOTATION AND SOME DEFINITIONS Linear time - invariant discrete digital system is considered, which is described by its state equation (2.1) x ( k + 1) = Ax ( k ) + Bu ( k ) n

m

where k ∈ N0, N0 = {0,1,2,...,}, x(k)∈ℜ is state vector at time k, u(k)∈ℜ is control n n n m vector at time k, m < n; A∈ℜ × and B∈ℜ × are real constant matrices and (A,B) is assumed to be controllable pair. A set of hyperplanes described by

ci x = 0 ; i = 1,2,..., m , n

where ci∈ℜ is a constant non zero row vector, is taken and jointed to the system (2.1). The intersection of these hyperplanes is the so-called sliding subspace S (in further text only S). Clearly, S is described by S = {x: Cx = 0} m n

where C∈ℜ × has i-th row equal to c j. Also, invertibility of matrix CB is assumed. The objective of this paper is to develop variable structure type of state feedback control law (2.2) u = u(x ) which guarantees that the state χ(k;x(0);u(⋅)) of the system (2.1) reaches S in finite time. Once S has been reached the controller is required to keep the state within it thereafter which means positive invariance of S relative to the system motion and what is denoted as sliding mode regime. During this regime, inside S, convergence to zero equilibrium with appropriately prescribed modes λi, can be guaranteed if the c j has been appropriately chosen. For c j choice see Utkin (1981). More rigorously , S is positive invariant relative to the system (2.1) motion iff (if and only if) x(0) ∈ S implies χ(k;x(0);u(⋅))∈S, ∀k∈N0. Further, some other notations and definitions are given for the reason of their usage in theorems which are main results and by means of which formulated problem is solved. n Real n-dimensional state space ℜ is with Euclidean norm denoted by || ⋅ ||. For S, d(x, S) = inf (||x − y||: y∈S) n

is the distance between S and the point x∈ ∈ℜ . Also for S

986

Z.M. BUČEVAC

d N 0 = (χ(k;x(0);u(⋅)),S) = sup[d (x(k),S): k∈N0] = sup{inf [||χ(k;x(0);u(⋅)) − y||: y∈S]: k∈N0} is distance between solution χ(k;x(0);u(⋅)) of system (2.1) and S on N 0 . Definition 2.1: The state x = 0 of the system (2.1), (2.2) is stable in S (with respect to S) iff ∀ε∈ℜ+(ℜ+ = [0,+∞[), ∃δ = δ(ε)∈ℜ+ such that x(0)∈S and ||x(0)|| < δ(ε) implies that χ(k;x(0);u(⋅)) exists, ||χ(k;x(0);u(⋅))|| < ε and χ(k;x(0);u(⋅))∈S, ∀k ≥ 0. Definition 2.2: The state x = 0 of the system (2.1), (2.2) is attractive (globally) in S (with respect to S) iff ∃∆ > 0 (∆ = +∞) such that x(0) ∈ S and ||x(0)|| < ∆ implies that χ(k;x(0);u(⋅)) exists, lim {||χ(k;x(0);u(⋅))||: k→+ ∞} = 0 and χ(k;x(0);u(⋅))∈S, ∀k∈N0. Definition 2.3: The state x = 0 of the system (2.1), (2.2) is (globally) asymptotically stable in S (with respect to S) iff it is both stable and (globally) attractive in S. Definition 2.4: The system (2.1), (2.2) is stable in S (with respect to S) iff its state x = 0 is globally asymptotically stable in S. Definition 2.5: S is stable, relative to the system (2.1), (2.2), iff ∀ε ∈ ℜ+, ∃δ = δ(ε)∈ℜ+ such that distance d[x(0),S] < δ implies that χ(k;x(0);u(⋅)) exists and d N 0 [χ(k ; x(0); u(⋅)), S ] < ε, ∀k ≥ 0. Definition 2.6: S is attractive (globally) relative to the system (2.1), (2.2) iff ∃∆∈]0,+∞[ (∆ = +∞) such that d[x(0), S] < ∆ implies lim{d[χ(k;x(0);u(⋅)),S]: k → +∞} = lim{d[x(k ), S]: k → +∞} = 0. Definition 2.7: S is (globally) asymptotically stable relative to the system (2.1), (2.2) iff S is stable and (globally) attractive at the same time. The problem stated in the objective of the paper can be reformulated as following: the objective of the paper is to develop the control law (2.2) such that S is globally asymptotically stable relative to the system (2.1), (2.2), and the system (2.1), (2.2) is stable in the sliding mode regime with appropriately prescribed modes. 3. MAIN RESULTS In this section solution of the problem stated in section 2 is provided by the following theorems. Theorem 3.1: Application of actual so-called sliding control usl(k) = −(CB)−1CAx(k) to the system (2.1) inside S is necessary and sufficient for S to be positive invariant relative to the system (2.1) solution χ(k; x(0); usl(k)). The application of actual sliding control usl(k) (in further text only usl(k)) to the system (2.1) inside S is phenomenon of discrete digital VSCS with sliding modes. Using this fact and intention for developing of variable structure type of controller for the system (2.1), control (2.2) components could be more precisely specified as

Discrete-Time Variable Structure Control Systems Multivariable Linear Plant Case

u i+ [x(k )], ci x(k )>0 u i [x(k )] = u isl [x(k )], ci x(k ) = 0 − i u i [x(k )], c x(k )<0

987

(3.1)

Evidently, the control is discontinuous which makes the right hand side of the system (2.1) state equation also discontinuous. Opposite to the usl(k), which is applied to the system (2.1) inside S, let us call control, which is applied to the system (2.1) outside of S, as outer control, uot(k). The following theorems are basic for solving of the stated problem by Lyapunov's second method. Theorem 3.2: Let the system (2.1), (3.1) be considered. S is globally asymptotically stable relative to the system (2.1), (3.1) if there exists scalar function v such that: n a) v(x)∈C(ℜ ) n b) v(x) > 0, ∀x∈ℜ \S c) v(x) = 0, ∀x∈S d) v(x) → ∞ as d(x, S) → +∞ n e) ∆v(x) < 0, ∀x∈ℜ \S f) ∆v(x) = 0, ∀x∈S where ∆v(x) is the first forward finite difference of v(x). Scalar function v defined by

v(x)=(signCx)T⋅Cx

as Lyapunov candidate and outer control law uot(k) = −(CB)−1[CAx(k) − Ψ(k)signCx(k)], Ψ = diag[ψ11,ψ22,…, ψmm, ], ψii(k) ≥ 0, ∀x∈N0 are chosen relative to the system (2.1). Evidently, such sort of uot together with usl, which is unique and already stated, represents control of (3.1) type. Using all these facts the following theorem was obtained as corollary to the previous theorem. Theorem 3.3: S is globally asymptotically stable relative to the system (2.1) if the following control is applied: u ot (k ) = −(CB) −1[CAx(k ) − Ψ (k )signCx(k )] , ψ (k ) = α | c i x(k ) | , α ∈ [0,1[ . u (k ) = ii sl −1 u (k ) = −(CB) CAx(k )

To be guaranteed for x(k) to reach S in finite time ∆v(x(k)) = −v(x(0))/P is taken with "P" as desired number of reaching steps. Taking the fact into account the following theorem is come to. Theorem 3.4: S is globally asymptotically stable relative to the system (2.1) if the following control is applied:

988

Z.M. BUČEVAC

u ot (k ) = −(CB) −1[CAx(k ) − Ψ (k )signCx(k )] m 1 j j ψ ii (k ) = ρii (k )∑ [| c x(k ) | − P | c x(0) |], ∀i = 1,, m ; ∀k = 0,, P − 1 j =1 m u (k ) = ρii (k ) ∈ [0,1[∧ ∑ ρ ii (k ) | signc i x(k ) |= 1 i =1 u sl (k ) = −(CB) −1 CAx(k )

4. PROOFS OF THEOREMS In this section proofs of theorems from the preceding section are provided. Proof of theorem 3.1: (Necessity) Positive invariance of S relative to the system (2.1) solution χ(k; x(0); usl(k)) is assumed. So, every time x(k)∈S, x(k+1)∈S also. Last fact, description of S and state equation of the system (2.1), give equation

CAx ( k ) + CBu sl ( k ) = 0 and further

u sl ( k ) = −(CB ) −1CAx ( k ) . (Sufficiency) Let assume, relative to system (2.1), x(k)∈S and control

u sl ( k ) = −(CB ) −1CAx ( k ) is applied. Then, state equation of system (2.1) and expression for S give

Cx ( k + 1) = CAx ( k ) − CB(CB ) −1CAx ( k ) = 0 which implies that x(k+1) ∈ S . Proof of theorem 3.2: The proof starts with assumption that sufficient conditions are fulfilled while stability and attraction of S should be shown. n The following notation is introduced, too. O(S)⊆ℜ is neighbourhood of S in meaning O(S) is superset of an open set containing S. O(.)(S) = {x: d(x,S) < (⋅)} is (⋅) − neighbourhood of S, (⋅) > 0. V(.) = {x: 0 ≤ ν(x) < (⋅)} is the largest connected neighbourhood of S with indicated property. ϕ(.)∈K is comparison function from class K, with property ϕ(ν1) < ϕ(ν2), ∀ν1,ν2∈[0, +∞[, ν1 < ν2. n Stability: As function v is positive definite on ℜ with respect to S, there exists open and connected set A = O(S) such that Vξ1 ⊂ Vξ2 ∧ ∂Vξ1 ∩ ∂Vξ2 = ∅ when ξ1 < ξ2 ∧ Vξ⊆A, ∀i = 1,2. Let ξm = max{ξ: Vξ⊆A} and ε* = max{ε: Oε(S)⊆Vξm}. Let ε∈]0, ε*]. As ν(x) is n n positive definite on ℜ relative to S, there exists ϕ ∈K such that ν(x) ≥ ϕ[d(x, S)] ∀x∈ℜ . Last statement and d(x,S) = ε, ∀x∈∂Oε(S) imply ν(x) ∂Oε ( S ) ≥ ϕ ( ε) . Now δ∈]0,+∞[ is chosen in the way Oδ(S) is the largest δ-neighbourhood inside Vϕ(ε) which finally implies

Discrete-Time Variable Structure Control Systems Multivariable Linear Plant Case

989

Oδ(S)⊆Vϕ(ε)⊆Oε(S)⊆Oε*(S)⊆Vξm⊆A. For initial state x(0)∈Oδ(S), ν[x(0)] < ϕ(ε) which means that hypersurface defined by ν[x(0)] = constans is inside Vϕ(ε). This fact, conditions d) and e) of the theorem and expression ∆ν[x(k)] = ν[χ(k+1;k;x(k))]− ν[χ(k;k;x(k))], x(k) = χ(k;0;x(0)) imply v[χ(k;0;x(0))] < ν[x(0)] < ϕ(ε) ∀k∈N0, ∀x(0) ∈Oδ(S)\S ∨ v[χ(k;0;x(0))] = ν[x(0)] < ϕ(ε) ∀k∈N0, ∀x(0) ∈ S. On the other hand v[χ(k;0;x(0))] ≥ ϕ{d[χ(k;0;x(0)),S]} ∀k∈N0. Last three inequalities give ϕ(ε) > ϕ{d[χ(k;0;x(0)),S]} ∀k∈N0, ∀x(0)∈Oδ(S) and further ε > d[χ(k;0;x(0)),S] ∀k∈N0, ∀x(0)∈Oδ(S), ∀ε∈]0,ε*]. For ε > ε*, δ(ε) = δ(ε*) is chosen. Previous analysis gives d[χ(k;0;x(0)),S] < ε* ∀k∈N0 and hence d[χ(k;0;x(0)),S] < ε ∀k∈N0. In that way it is shown that all what is stated by definition 2.5 is fulfilled, that is, that stability of S is proved. Attraction: S is assumed not to be globally attractive while the conditions of the theorem are fulfilled. The consequence is that lim d[χ(k;0;x(0)),S] = v, v∈]0,∞[ for some k →∞

x(0), x(0) is such that d[x(0),S] < ∞. In that case there exists α∈]0, v] such that α = inf{d[χ(k;0;x(0)),S]: k∈N0}. Let ξ1 = inf[ν(x): x∈∂Oα(S)]. Previous assumptions n imply that inf{v[χ(1;0;x(0))]: k∈N0} ≥ ξ1 for x(0)∈ℜ , that is, working point of (2.1), n (3.1) system never comes into Vξ1, χ(k;0;x(0)) ∈ℜ \Vξ1, ∀k∈N0. Condition e) of the n

theorem implies that ∆ν(x) < 0 ∀x∈ℜ \Vξ1 so that there exists supremum of function ∆ν n

n

on set ℜ \Vξ1, sup[∆ν(x): x∈ℜ \Vξ1] = η< 0. This further gives v[χ(1;0;x(0))] − ν[x(0)] ≤ η v[χ(2;1;x(1))] − ν[x(1)] ≤ η

v[χ(k;k − 1;x(k −1))] − ν[x(k − 1)] ≤ η. By summing all these relations, the following expression is got: ν[x(k)] − ν[x(0)] ≤ kη, that is, ν[x(k)] ≤ ν[x(0)] + kη. It implies that ν[x(k)]< 0 ∀k ≥ [0,ν[x(0)]/(−η)]N + 1 where [0, ν[x(0)]/(−η)]N denotes the biggest integer of the indicated segment. The last statement is contradictory to the b) condition of the theorem. This is caused by failed assumption that S is not globally attractive while conditions of the theorem are fulfilled, that is, S is globally attractive. Proof of theorem 3.3: Previously introduced Lyapunov's function candidate fulfills evidently conditions a) − d) of theorem 3.2. Let be shown that conditions e) and f) are fulfilled. First, state x(k) of the system (2.1) is adopted to be out of S. Then ∆ν[x(k)] = ν[x(k+1)]-ν[x(k)] = = {sign[Ψ(k)signCx(k)]}TΨ(k)signCx(k) − [signCx(k)]TCx(k) = m

1− 1T | Cx( k )| = (α−1) 1T | Cx(k ) |< 0

= ∑ | ψ ii ( k )||⋅ signci x ( k )|− 1T | Cx ( k )| ≤ 1T Ψ( k ) i =1

m

where |Cx(k)| = [|c1x(k)|, |c2x(k)|, …, |cmx(k)|]T, 1 = (1 1 …1)T, 1∈ℜ . Second, state x(k) of the system (2.1) is adopted to be inside of S. Then ∆ν[x(k)] = {sign[CAx(k) − CB(CB)-1CAx(k)]}T[CAx(k) − B(CB)−1CAx(k)] −1T|Cx(k)| = 0 .

990

Z.M. BUČEVAC

Proof of theorem 3.4: Differently from previous theorem, concerning condition e) of theorem 3.2, in this theorem decreasing of ν[x(k)] along motion of (2.1) with prescribed "velocity" in section 3 should be shown. State of system (2.1) at time k is adopted to be out of S. Then T ∆ν[x(k)] = 1 | Ψ (k )signCx(k ) | − 1T | Cx ( k )| = −

T

1

Cx(k) P

that is, m

∑ ψ ii (k ) | sign c i x(k ) | = i =1

, k = 0, 1, …, P −1,

j | c j x (0 ) | | ( ) | k − c x . ∑ P j =1 m

Further, including expression for ψii(k) left-hand side of the last relation becomes: i =1 m

m

j =1

1

m

j =1

1

∑ ρii (k ) ∑ | c j x(k ) | − P | c j x(0) | ⋅ | signci x(k ) |= ∑ | c j x(k ) | − P | c j x(0) | ⋅ m

m

j =1

1

∀k = 0, 1, 2,..., P − 1

∑ ρii (k ) | signci x(k ) |=∑ | c j x(k ) | − P | c j x(0) | , i =1

in which way the proof is finished. 5. EXAMPLE In this section an example is provided to illustrate presented results. Example 4.1: Let us consider unstable system 1 1 − 1 1 2 x(k + 1) = − 2 1 − 1 x(k ) + 2 − 1u(k ) − 2 − 1 1 − 2 1

(4.1)

with controllable pair (A,B). S, which guarantees prescribed mode λ = 0.9 in sliding mode regime, is designed by analogical procedure to the known one in literature; see [16]. This sliding subspace S is defined by 1 1.9 −2.1 C= . -1 -1.1 0.9 Let be chosen, in presented control algorithm 1 , l (k ) ≠ m ρii (k ) = ρ(k ) = m − l (k ) 0, l (k ) = m where l(k) is number of zero components in vector Cx(k). Further, "P" is adopted to be 10. System (4.1) with developed control algorithm included becomes:

Discrete-Time Variable Structure Control Systems Multivariable Linear Plant Case

991

−2 0.9 − 0.6 0.9 −1 sign[σ1 (k )] 0.5 ⋅ x(k + 1) = − 1.8 12 − 1.8 x(k ) + ρ(k ){Ω(k ) − 0.1Ω(0)} 0.5 sign[σ 2 (k )] − 1.2 0.8 − 1.2 − 0.5 − 0.5 where, σ1 (⋅) = x1 (⋅) + 1.9 x2 (⋅) − 2.1x3 (⋅) , σ 2 (⋅) = − x1 (⋅) − 1.1x2 (k ) + 0.9 x3 (⋅) , Ω(⋅) =| σ1 (⋅) | + | σ 2 (⋅) | . Solving of the discrete digital system was done by digital computer. As illustration, simulation results from the initial state x(0) = (5 5 10)T , controls, c1x(k) and c2x(k) are given in figures 4.1, 4.2, 4.3 and 4.4 respectively. Evidently at time k = 10 working point reaches S, staying afterwards in it for ever. x1 x2 x3

10

x1(k), x2(k), x3(k)

5 0 -5 -10 -15 -20 0

5

10

15

20

25

k

Fig.4.1 u1 u2

10 0 -10 -20

u1(k), u2(k)

-30 -40 -50 -60 -70 -80 -90 -100 0

5

10

15

k

20

25

Fig.4.2

992

Z.M. BUČEVAC 4

3

c1x(k)

2

1

0 0

5

10

15

20

25

k

Fig.4.3

0

c2x(k)

-1

-2

-3

0

5

10

15

20

25

k

Fig.4.4 6. CONCLUSION In the paper, discrete digital type of VSCS with sliding mode is considered but this time related to the most general case of linear time invariant plant with no restrictions on the form of state equation and number of controls. The following phenomenon, characteristic only for this type of VSCS with sliding modes, is discovered: it is necessary to apply actual control so-called sliding control usl to the plant during sliding mode regime, what is not characteristic for continuous time VSCS with sliding modes. Solution of this type of problem makes place for microprocessor compensator application as compensator in the systems. Acknowledgment. The author would like to express deep gratitude to Prof. Lj.T.Grujic for all sorts of support in research the paper resulted from.

Discrete-Time Variable Structure Control Systems Multivariable Linear Plant Case

993

REFERENCES 1. Ambrosino, G., Gelentano, G., Garofalo, F.: Variable Structure Model Reference Adaptive Control Systems, International Journal of Control, Vol. 39, pp.1339-1349, 1984. 2. Bartolini, G., Zolezzi, T.: Some Developments of Variable Structure Control for Nonlinear Systems, Preprints of 25 IEEE Conference on Decision and Control, Athens, Greece, 1986. 3. Bartolini, G., Ferrara, A., Utkin, V. I.: Adaptive Sliding Mode Control in Discrete Time Systems, Automatica, Vol. 31, No. 5, pp. 769-773, 1995. 4. Bucevac, Z.: Design of Digital Discrete Control Systems with Sliding Mode - Ph. D. Dissertation (in Serbian), Mechanical Engineering Department, Belgrade University, Belgrade, Yugoslavia, 1985. 5. Bucevac, Z.: Lyapunov's Method Approach in a Stabilizing Control Algorithm Design for Digital Discrete VS Systems with Sliding Modes - Linear Plant Case, Proceedings of the 7. IASTED International Symposium Modelling, Identification and Control, Grindelwald, Switzerland, l988. 6. El-Khazali, R., DeCarlo, R.: Output Feedback VSCS Design, Automatica, Vol.31, No.6, 805-816, 1995. 7. Emel'yanov, S.V.: Variable Structure Automatic Control Systems (in Russian), Nauka, Moscow, 1967. 8. Emel'yanov, S.V., Utkin, V.I.,Taran,V.A., Kostyleva,N.E., Subladze, A.M., Ezerov, V.B., Dubrovskij, E.N.: Theory of Variable Structure Systems (in Russian), Nauka, Moscow, 1970. 9. Fluge-Lotz, I.: Discontinuous Automatic Control Systems, Academic Press, Princeton, New Jersey, 1953. 10. Furuta, K.: Sliding Mode Control of a Discrete System, System and Control Letters, Vol. 14, pp. 145-152. 1990. 11. Gao, W., Wang, Y., Homaifa, A.: Discrete Time Variable Structure Control Systems, IEEE Transactionon on Industrial Electronics, Vol.42, No.2, 117-122, 1995. 12. Geering, H.P., Guzzella, L.: Model Following Variable Structure Control for a Classof Uncertain Mechanicalth Systems, Preprints of 25th IEEE Conferenceon Decision and Control, Athens, Grecce, 1986. 13. Gough, E.N., Ismail,M.Z., King,E.R.: Analysis of Variable Structure Systems with Sliding Modes, International Journal of Systems Science, Vol. 15, pp. 401- 409, 1984. 14. Grujic, Lj., Ribar, Z.: Nonstationary Lourie Variable Structure Control System Design, II SAUM conference, Belgrade, 1986. 15. Itkis, Y.: Control Systems of Variable Structure, New York, Wiley, 1976. 16. Itkis,Y.: Dynamic Switching of Type-I/Type-II / Structures in Tracking Servosystems, IEEE Transactionon Automatic Control, AC-28, pp. 531-534, 1983. 17. Karakasoglu, A., Sundareshau, M.K.: A Recurrent Neural Network-based Adaptive VS Model Following Control of Robotic Manipulators, Automatica, Vol.31, No.10, 1995. 18. Khurana, H., Ashon, S.I., Lamba, S.S.: On Stabilization of Large-Scale Control Systems Using Variable Structure Systems Theory, IEEE Transactionon Automatic Control, AC-31, pp. 176-178, 1986. 19. Kotta, U.: Comments on "On the Stability of Discrete Time Sliding Mode Control Systems", IEEE Transactionon Automatic Control, Vol.-34, No.9, pp. 1021-1022, 1989. 20. Milosavljevic, C.: General Conditions for the Existence of a Quasisliding Mode on the Switching Hyperplanein Discrete Variable Structure Systems, Automation and Remote Control, No. 3, pp. 36-44, 1985. 21. Li, K-C., Leung, T-P., Hu, Y-M.: Sliding Mode Control of Distributed Parameter Systems, Automatica, Vol.30, No.12, 1994. 22. Ryan, E.P.: A Variable Structure Approach Regulation of Uncertain Dynamical Systems, International Journal of Control, Vol.38, pp 1121-1134, 1983. 23. Sarpturk, S. Z., Istefanopulos, Y., Kaynak, O.: On the Stability of Discrete Time Sliding Mode Control Systems, IEEE Transactionon Automatic Control, AC-32, No.10, pp. 930-932, 1987. 24. Sira-Ramirez, S.: Non linear discrete variable structure systems in quasi-sliding mode, Int. J. Control, Vol.54, pp.1171-1187, 1991. 25. Slotine ,J.J., Sastry, S.S.: Tracking Control of Nonlinear Systems Using Sliding Surfaces with Application to Robot Manipulators, International Journal of Control, Vol.38, pp.465-492, 1983. 26. Slotine, J.J.E., Hedrick, J.K., Misawa, E.: On Sliding Observers for Nonlinear Systems, Preprints of 25th IEEE, Conference on Decision and Control, Athens, Greece, 1986. 27. Utkin, V.I.: Sliding Modes and Their Applicationsin Variable Structure Systems (in Russian), Nauka, Moscow, 1974. 28. Utkin, V.I.: Sliding Modes in Optimization and Control Tasks (in Russian), Nauka, Moscow, 1981. 29. Xing-Huo, Y., Chun-Bo, F.: A Lyapunov-function Approach for the Design of Variable-Systems, Preprints of 10th World Congress for Structure Automatic Control, Munich,W.Germany,1987.

994

Z.M. BUČEVAC

30. Young, K-K.D.: Asymptotic Stability of Model Reference Systems with Variable Structure Control, IEEE Transaction on Automatic Control, AC-22, pp. 279-281, 1977. 31. Young,K-K.D.: Design of Variable Structure Model Following Central Systems, IEEE Transaction on Automatic Control, AC-23, pp. 1079-1085, 1978. 32. Zheung, F., Chen, M., Gao, W.: VSC of Time Delay Systems with a Simulation Study on Stabilizing Combustion in Liquid Propellant Rocket Motors, Automatica, Vol.31, No.7, 1995.

DIGITALNI SISTEMI AUTOMATSKOG UPRAVLJANJA PROMENLJIVI STRUKTURE - SLUČAJ MULTIVARIJABILNOG LINEARNOG OBJEKTA Zoran M. Bučevac Razmatra se linearni stacionarni diskretni digitalni objekt bez ograničenja na oblik njegove jednačine stanja i broj upravljanja. Razvijen je stabilizirajući algoritam upravljanja sa povratnom spregom po stanju pomoću Ljapunovljeve Druge metode koji daje sistem promenljive strukture sa kliznim radnim režimom. Ovaj algoritam upravljanja obezbeđuje da se stanje sistema kreće iz proizvoljnog početnog položaja do propisanog takozvanog podprostora klizanja S, za konačno vreme. Unutar podprostora klizanja S sistem prelazi u klizni radni režim i u njemu ostaje zauvek.