arxiv: v1 [cs.sy] 16 Nov 2014

Transcription

1 Estimation of transient process for singularly perturbed synchronization system with distributed parameters Vera Smirnova 1, Anton V. Proskurnikov and Natalia V. Utina 3 arxiv: v1 [cs.sy] 16 Nov 14 Abstract Many systems, arising in electrical and electronic engineering are based on controlled phase synchronization of several periodic processes ( phase synchronization systems, or PSS). ypically such systems are featured by the gradientlike behavior, i.e. the system has infinite sequence of equilibria points, and any solution converges to one of them. his property however says nothing about the transient behavior of the system, whose important qualitative index is the maximal phase error. he synchronous regime of gradient-like system may be preceded by cycle slipping, i.e. the increase of the absolute phase error. Since the cycle slipping is considered to be undesired behavior of PSSs, it is important to find efficient estimates for the number of slipped cycles. In the present paper, we address the problem of cycle-slipping for phase synchronization systems described by integro-differential Volterra equations with a small parameter at the higher derivative. New effective estimates for a number of slipped cycles are obtained by means of Popov s method of a priori integral indices. he estimates are uniform with respect to the small parameter. Index erms singularly perturbed systems, Asymptotic properties, Popov-type stability of feedback systems, frequencyresponse methods. I. INRODUCION A vast range of physical and mechanical systems are described by ordinary differential and integro-differential equations with a small parameter at the higher derivative. Such equations are usually called singularly perturbed, since the order of unperturbed equation is lower than the order of the perturbed one. So in the absence of special simplifying assumptions, the electron generator is described by a third order differential equation with a small parameter [], Van der Pol equation being the special case when the small parameter is equal to zero. he asymptotic properties of singularly perturbed equations may differ from these of unperturbed ones. So the problems of stability and oscillations for various singularly perturbed integro-differential equations became the subject for special research [4]. he armamentarium of singularly perturbed systems is a backbone of the mathematical of *he work was supported by the European Research Council (ERCStG- 377) and RFBR, grant he example (Section V) was supported solely by Russian Scientific Foundation, grant Vera Smirnova is with the Department of Mathematics, St.Petersburg State University of Architecture and Civil Engineering. She is also with the Department of Mathematics and Mechanics, St.Petersburg State University, St.Petersburg, Russia, root@al189.spb.edu Anton Proskurnikov is with the Research Institute of echnology and Management, University of Groningen. He is also with St.Petersburg State University, IMO University and Institute for Problems of Mechanical Engineering RAS, St.Petersburg, Russia avp198@gmail.com 3 Natalia Utina is with the Department of Mathematics, St.Petersburg State University of Architecture and Civil Engineering. unv74@mail.ru theory of systems with different time scales, arising in different areas of engineering and natural sciences, from mechanics to mathematical biology [5]. In this paper we examine asymptotic behavior of singularly perturbed phase synchronization systems (PSS). he PSSs are based on the principle of phase synchronization [7]. hese systems, sometimes referred to as synchronous or pendulum-like control systems, involve periodic nonlinearities and typically have infinite sequence of equilibria points. An important class of PSS is constituted by phase-locked systems, which are based on the seminal idea of phaselocked loop (PLL) and widely used in telecommunications and electronics [9], [8]. Recent decades a vast literature examining asymptotic behavior and other dynamical properties of PSSs has been published, motivated by numerous applications of these systems in mechanical, electric, electronic and telecommunication engineering. Most of these papers address the problem of gradient-like behavior, aiming at obtaining conditions which guarantee convergence of the solutions to equilibria, which means that the generators of the system are synchronized for any initial state. For details and bibliography see e.g. [7], [1] and references therein. But as a rule the synchronous regime of gradient-like system is preceded by cycle slipping, i.e. the increase of the absolute phase error. Its amplitude depends on the initial state of PSS and is an important characteristic of the transient process of the system. he phenomenon of cycle slipping was set forth in the book [16], for mathematical pendulum with viscous friction proportional to the square of angular velocity. For mathematical pendulum the number of full rotations around the point of suspension was called the number of slipped cycles. he extension of this notion to PSSs is as follows. Suppose that a gradient-like phase synchronization system has a - periodic input and let σ(t) be its phase error. hey say that the output function σ(t) has slipped k N } cycles if there exists such a moment t that however for all t one has σ( t) σ() = k, (1) σ(t) σ() < (k +1). () So to give the adequate description for behavior of PSSs one must establish possibly close estimates for the number of slipped cycles. And since large number of slipped cycles is undesirable for PSSs the problem of its estimation is important.

2 In the paper [3] the problem of cycle slipping was considered for multidimensional PSSs. By periodic Lyapunov-like functions and the Kalman-Yakubovich-Popov (KYP) lemma some frequency-algebraic estimates were obtained. he results of [3] were formulated in terms of LMI-solvability in [17]. he estimates of [3] were extended to discretetime and distributed parameter PSSs in the paper [15] and the monograph [6] respectively. Distributed parameter PSSs were investigated by the method of a priori integral indices with the help of Popov-like functionals of special type. In paper [11] the generalized Popov-like functionals from [13] are exploited for estimation of cycle-slipping for distributed parameter PSSs. In this paper we make frequency-algebraic estimates from [11] more exact and extend them for singularly perturbed distributed parameter PSSs. II. PROBLEM SEUP Consider a distributed parameter synchronization system which is described by the integro-differential Volterra equation with a small parameter at the higher derivative: µ σ µ (t)+ σ µ (t) = α(t)+ρϕ(σ µ (t)) t γ(t τ)ϕ((σ µ (τ))dτ (t ). Here µ >, h, ρ R, γ,α : [,+ ) R; ϕ : R R. he map ϕ is assumed C 1 -smooth and -periodic with two simple isolated roots on [, ). he kernel function γ( ) is piece-wise continuous, the function α( ) is continuous. For each µ the solution of (3) is defined by specifying initial condition (3) σ µ (t) t [,] = σ (t), (4) σ ( ) is continuous and σ(+) = σ (). We assume for definiteness that ϕ(σ)dσ. (5) We assume also that the linear part of (3) is stable: Let α(t) + γ(t) Me rt (M,r > ). (6) dϕ α 1 = inf σ [, ) dσ, α dϕ = sup σ [, ) dσ. α 1 dϕ dσ α, σ R. Notice that α 1 < < α. In paper [1] some sufficient frequency-algebraic conditions for gradient-like behavior of equation (3) were obtained, uniform with respect to the parameter µ. hey are formulated in terms of the transform function of the linear part of the unperturbed Volterra equation (µ = ) t σ(t) = α(t)+ρϕ(σ(t)) γ(t τ)ϕ((σ(τ))dτ. (7) he transfer function for (7) is defined as follows t K(p) = ρe p + γ(t)e pt dt (p C). (8) Our goal is to estimate maximal deviations sup σ µ (t) σ µ () for a gradient-like system. Precisely, t we are going to obtain estimates for the number of slipped cycles. III. FREQUENCY-ALGEBRAIC ESIMAES IN CASE OF HE NUMBER OF CYCLES SLIPPED FOR UNPERURBED EQUAION Consider a distributed parameter system with a periodic nonlinearity, described by the integro-differential equation (7). his equation is the unperturbed equation for the equation (3). Its transfer function from the input ϕ to the output ( σ) is defined by (8). he initial condition for (7) has the form σ(t) t [,] = σ (t). (9) Let us obtain certain estimates for the number of slipped cycles for the unperturbed equation. We start with the following technical lemma which is a cornerstone point in estimating the number of slipped cycles. Lemma 1. Suppose there exist such positive ϑ, ε, δ, τ that for all ω the frequency-domain inequality holds: ReϑK(iω) τ(k(iω)+α 1 1 iω) (K(iω)+α 1 iω)} ε K(iω) δ (i = 1). (1) the following quadratic functionals I = ϑ σ(t)ϕ(σ(t)) +ε σ (t)+δϕ (σ(t))+ +τ(α 1 1 ϕ(σ(t)) σ(t))(α 1 ϕ(σ(t)) σ(t)) } dt are uniformly bounded along the solution of (7): I Q, (11) Q does not depend on. Proof: Let σ(t) be the solution of (7), (9) and η(t) = ϕ(σ(t)). Let be an arbitrary positive number. If η() we determine v(t) = for t <, t for t [,1], 1 for t > 1 in case η() = we put v(t) = 1. We introduce auxiliary functions, σ (t) := α(t)+(1 v(t))ρη(t) ζ (t) := t (1 v(τ))γ(t τ)η(τ)dτ, η(t) t, η()e λ( t) t > (λ > ) η (t) := v(t)ζ (t), },

3 t σ (t) := ρη (t) For t [,] we have γ(t τ)η (τ)dτ. σ(t) = σ (t)+σ (t). (1) For any >, consider a functional ρ := + ϑσ (t)η (t)+δη (t)+εσ (t)+ +τ(σ (t) α1 1 η (t))(σ (t) α 1 η (t)) } dt. It is demonstrated in [6] by means of the Popov method of a priori integral indices [14] that the frequency-domain inequality (1) implies It can be easily shown that ρ >. (13) ρ = I +I 1 +I +I 4, (14) the addends I 1,I,I 4 are defined by l I 1 := +α 1 1 α 1 τ η (t) ϑ(1 v(t)) σ(t)η(t)+δ(1 µ )η (t)+ ( +τ ( ( α 1 1 +α 1 ) σ(t) ) v(t)η(t) α 1 1 α 1 τ+ )} η(t) v(t)η(t) dt, with l =, if < 1, and l = 1, if 1; I := ϑσ (t)η(t)v(t) (ε+τ) σ(t)σ (t)+ + (τ +ε) σ (t)+τ ( α 1 1 +α 1 ) σ (t) η(t) } dt; I 4 := ϑσ (t)η (t)+δη (t)+(ε+τ)σ (t) ( α 1 1 +α 1 ) τσ (t) η (t)+τα1 1 α 1 η (t)} dt. From (13) and (14) it is immediate that I I 1 I I 4. (15) As can be seen from (6), the functionals I 1 and I are uniformly bounded: I 1 q 1 (δ,ϑ,τ,α 1,α ), (16) I q (ϑ,ε,τ,α 1,α ), (17) q 1 and q do not depend on. By definition of η (t), I 4 = ϑσ (t)η()e λ( t) +δη ()e λ( t) + +(ε+τ)σ (t)+λ( α 1 1 +α 1 ) τσ (t)η()e λ( t) + +λ τα 1 1 α 1 η ()e λ( t)} dt. Introducing the constant δ α1 α λ :=, (18) τ it can be shown that I 4 W 4(ε+τ) η ()e λ( t) dt, W = ϑ+ τδ α 1 α ( α 1 1 +α 1 I 4 q 3 (δ,ϑ,τ,α 1,α ) = ), and hence τ W maxϕ (σ) 8(ε+τ) δ α 1 α. (19) he inequalities (15), (16), (17), and (19) yield that I q 1 +q +q 3 Q(ϑ,δ,ε,τ,α 1,α ), () which finishes the proof of Lemma 1. he value of Q from (11) may be found explicitly. Lemma. Let α 1 = α and ϕ(σ()) = ϕ(σ()) =. if the conditions of Lemma 1 are fulfilled the following estimate holds I q := 1 r (ϑmm+(ε+τ)mm(m r +ρ)+(ε+τ)m ), (1) m = supϕ(σ). Proof: Let us exploit Lemma 1 for ϕ(σ()) = and ϕ(σ()) =. In this case I 1 = and I 4. So it follows from (15) that Since α 1 = α I I I. () ϑα(t)ϕ(σ(t)) +(τ +ε)α (t)+(τ +ε)α(t) (ρϕ(σ(t τ)) (τ +ε) t γ(t τ)ϕ(σ(τ))dτ) dt Hence it follows that the estimate (1) is true. (3) he next two theorems give estimates for the number of slipped cycles. o start with, we introduce auxiliary functions Φ(σ) = (1 α 1 1 ϕ (σ))(1 α 1 ϕ (σ)), r 1j (k,ϑ,ε,τ,x) = P(ε,τ,σ) = ε+τφ (σ), ϕ(σ)dσ +( 1) j x ϑk ϕ(σ) P(ε, τ, σ)dσ (j = 1,) Y j (σ) = ϕ(σ) r 1j ϕ(σ) P(ε,τ,σ) (j = 1,).

4 heorem 1. Suppose there exist such positive ϑ, ε, δ, τ and natural k that the following conditions are fulfilled: 1) for all ω the frequency-domain inequality (1) ) 4δ > ϑ (r 1j (k,ϑ,ε,τ,q)) (j = 1,), (4) Q is given by (11). any solution of (7) slips less than k cycles, that is, the inequalities hold σ() σ(t) < k t. (5) he proof of heorem 1 is presented in [11]. o proceed with the next result, we introduce the following functions for j = 1, ϕ(σ)dσ +( 1) j x ϑk r j (k,ϑ,x) :=, ϕ(σ) dσ ϕ(σ)dσ +( 1) j x ϑk r j (k,ϑ,x) :=, Φ(σ) ϕ(σ) dσ and matrices j (k,ϑ,x) := ε := aϑr j (k,ϑ,x) F j (σ) = ϕ(σ) r j ϕ(σ), Ψ j (σ) = ϕ(σ) r j ϕ(σ) Φ(σ) aϑr j (k,ϑ,x) δ a ϑr j (k,ϑ,x) a ϑr j (k,ϑ,x), τ a [,1] and a := 1 a. heorem. Suppose there exist positive ϑ, ε, δ, τ, a [, 1] and natural k satisfying the conditions as follows: 1) for all ω the frequency-domain inequality (1) ) the matrices j (k,ϑ,q) (j = 1,) the value of Q is defined by (11), are positive definite. for the solution of (7) the inequality (5) holds. Proof: Let σ(t) be the solution of (7), (9). Let ε > be so small that matrices j (Q + ε ) are positive definite. We consider the functions F j (σ), Ψ j (σ) (j = 1,) with r j = r j (k,ϑ,q+ε ), r j = r j (k,ϑ,q+ε ) and the functionals I from Lemma 1. It is true that I = ϑa σ() σ() σ() F j (σ)dσ +ϑa σ() Ψ j (σ)+ + ϑ σ(t)ϕ(σ(t))+ε σ (t)+δϕ (σ(t)) ϑaf j (σ(t)) σ(t) ϑa Ψ j (σ(t)) σ(t)+ +τ σ (t)φ(σ(t)) } dt (j = 1,). (6) In virtue of condition ) of the theorem, the third term in the right hand part of (6) is the integral of positive definite quadratic form. So ( I ϑ a σ() σ() F ) σ() j(σ)dσ +a σ() Ψ j(σ)dσ (7) (j = 1,). Suppose that σ(t1) F 1 (σ)dσ = k σ() σ(t1) σ() Ψ 1 (σ)dσ = k σ(t 1 ) = σ()+k. F 1 (σ)dσ = 1 ϑ (Q+ε ), Ψ 1 (σ)dσ = 1 ϑ (Q+ε ). I t1 Q+ε > Q. (8) which contradicts with (). So our hypothesis is wrong. With the help of F (σ) and Ψ (σ) we prove that As a result for all t > σ(t) σ() k. σ() k < σ(t) < σ()+k. heorem is proved. heorem 3 Let α 1 = α and σ() = σ ϕ(σ ) =. Suppose there exist such positive ϑ, ε, δ, τ, a [,1] and natural k that the following conditions are fulfilled: 1) for all ω the frequency-domain inequality(1) ) the matrices j (k,ϑ,q) (j = 1,), with q defined by (1), are positive definite. for any solution of (7) the estimate (5) holds. Proof: Let us repeat the proof of heorem with r j = r j (k,ϑ,q +ε ), r j = r j (k,ϑ,q +ε ). I tj > q (j = 1,). (9) Since ϕ(σ()) = we conclude that ϕ(σ(t 1 )) = and ϕ(σ(t )) =. it follows from Lemma that which contradicts (9). I tj q (j = 1,),

5 IV. EXAMPLE Let us consider a phase-locked loop (PLL) with a proportional integral lowpass filter, a sine-shaped characteristic of phase frequency detector and a time-delay in the loop. Its mathematical description is borrowed from [1]: σ(t)+ 1 σ(t)+ϕ(σ(t))+s ϕ(σ(t)) =, (3) ϕ(σ) = sinσ β, s (,1), β (,1], h >, >. he differential equation (3) can be reduced to integrodifferential equation (7) with }, t < h, γ(t) =, (1 s)e t, t h α(t) = e t (b (1 s)j), b = σ()+sϕ(σ()) and t e λ+h ϕ(σ(λ))dλ, t h, J =. e λ+h ϕ(σ(λ))dλ, t > h he transfer function of the lowpass filter here has the form: K(p) = sp+1 p+1 e ph We suppose that ϕ(σ()) = and apply heorem 3. Let α = α 1 = 1, ϑ = 1, a = 1. he assumption 1) of heorem 3 shapes into Ω(ω) τ ω 4 +ω ( 3 scosωh 4 s (ε+τ)+ +τ δ ) (1 s)ωsinωh+ cosωh (ε+τ) δ ω; as condition ) may be rewritten as εδ > Notice that for all ω R one has (31) πβ +q k 1 4(βarcsinβ + 1 β ). (3) Ω(ω) Ω (ω) (τ 1 3 sh )ω 4 +( 3 s 4 s (ε+τ)+τ δ 1 h (1 s) h)ω + +( (ε+τ) δ), ω and Ω(ω) Ω (ω) when ωh << 1. We consider the case.9, h = h 1, since for small and small h the PLL is gradient-like for all β (,1] [1]. Let us choose ε = β, δ = α, τ = γ 3. As Ω() = Ω () it is necessary that α + β + γ 4 1. the optimal values for α and β are α = β = 1 (1 γ 4 ), whence εδ = 1 γ 4. For γ = max 1 sh, 1 (h +1 s) } the polynomial Ω (ω) is nonnegative, ω. It follows from (3) that the number k of cycles slipped satisfies the inequality k r := q (8 εδ(βarcsinβ + 1 β ) πβ) 1, x stands for the integer floor of x. Let us consider the PLL with b = K()β [3]. by estimating the functional I we conclude that the value of q can be defined by the formula q = (A+Bh +Ch ), (33) A = ( 7 β +3), B = 3(1 s)(1+β)(3β +1), C = 3 (1 s) (1+β). (34) It follows from (33), (34) that the number of slipped cycles increases together with, withβ or withh. Let for example h = 1, s =.4, =.1. r = 1 for β =.9, r = for β =.9, and r = 5 for β =.95. V. ESIMAES FOR HE NUMBER OF CYCLES SLIPPED FOR SINGULARLY PERURBED EQUAION Equation (3) can be reduced to integro-differential Volterra equation t σ µ (t) = α µ (t) γ µ (t τ)ϕ(σ µ (τ))dτ (t > ), (35) α µ (t) = σ()e t 1 µ + e λ t ρ µ α(λ)dλ+ µ µ J, (36) t e λ+h t µ ϕ(σ(λ))dλ, t h, J =. γ µ (t) = 1 µ Let t e λ+h t µ ϕ(σ(λ))dλ, t > h, e λ t h t ρ µ γ(λ)dλ e µ, t h, µ, t < h he transfer function for equation (35) is as follows }. (37) K µ (p) = K(p) 1+µp. (38) q = q+(ϑm+(ε+τ)m( M r +ρ))ρmh+(ε+τ)ρ m h, q is defined by (1). heorem 4 Let α = α 1 and σ() = σ ϕ(σ ) =. Suppose there exist such positive ϑ, ε, δ, τ, a [,1] and natural k that the following conditions are fulfilled: 1) for all ω the frequency-domain inequality(1) ) the matrices j (k,ϑ,q ) (j = 1,) are positive definite. there exists such value µ that for all µ (,µ ) the following assertion is true: for any solution of (7) the estimates σ µ () σ µ (t) < k t (39)

6 hold. Proof: Let α 1 = α = κ. In our case the inequality (1) takes the form Π(ω), (4) Π(ω) := τκ ω +ϑrek(iω)} (ε+τ) (K(iω) δ. (41) Consider δ < δ and Π(ω) = Π(ω)+δ δ. We have Π(ω) > ω. (4) Let us substitute in matrices j (k,ϑ,q ) the value δ by δ and denote the new matrices by j (k,ϑ,q ) Let δ δ be so small that j (k,ϑ,q ) are positive definite. For integro-differential Volterra equation (35) one can apply heorem 3. For transfer function (38) the frequencydomain inequality (1) with ϑ, ε, δ, τ takes the form or Π µ (ω) := ReϑK µ (iω)} (ε+τ) (K µ (iω) + +τκ ω δ. (43) Π(ω)+ ϑµωimk(iω)}+τµ κ ω 4 δµ ω. (44) Introduce a constant δ Ω := κ τ. (45) Let us fix a certain µ >. Since the value of K(iω) is bounded for ω R we can assert that there exists such number Ω > Ω that the inequality (44) is true for ω > Ω, µ µ. Let if δ 1 = inf Π(ω), ω [,Ω ] L 1 = sup ϑωimk(iω)}. ω [,Ω ] µ < min (46) } δ 1 δ1 τ,, µ, (47) L 1 κ δ the inequality (44) is true for any µ < µ and all ω. Consider now the functional I from Lemma 1 for the solution σ µ (t) of equation (35) and the nonlinear function ϕ(σ µ (t)). he estimate (3) for nonlinear function ϕ(σ µ (t)) and the term α µ (t)) takes then the form Hence I +(τ +ε) ϑα µ (t)ϕ(σ µ (t))+(τ +ε)α µ (t)+ t γ µ (t τ)ϕ(σ µ (τ)dτ) dt I q µ, (48) q µ = (ϑm+(ε+τ)m(ρ+ M r ))(µ σ() + M r +ρmh)+ M +(ε+τ)( µ 4µ σ ())+ (1 rµ) (µ 1 rµ + 1 r )+ +ρ m (h+µe µ µ)). (49) It is clear that limq µ = q (5) as µ. Since matrices j (k,ϑ,q ) are positive definite there exists a value ˆµ small enough that for µ ˆµ matrices j (k,ϑ,q µ ) are positive definite. hus if µ < µ := min µ, ˆµ} the frequency inequality of heorem 3 for the transfer function K µ (p) and the algebraic restrictions on the varying parameters are fulfilled. So for σ µ (t) the estimate (39) is true. VI. CONCLUSION he paper is devoted to the problem of cycle-slipping for singularly perturbed distributed parameter phase synchronization systems. he PSSs described by integro-differential Volterra equations with a small parameter at the higher derivative are addressed. he case of differentiable nonlinearities is considered. he problem is investigated with the help of the method of a priori integral indices. In the paper new effective multiparametric frequency-algebraic estimates for the number of slipped cycles of the output of the system are established. he estimates obtained are uniform with respect to a small parameter. REFERENCES [1] L.N. Belustina. Qualitative-numerical analysis for synchronization system with time-delay (in russian). In Proceedings of International Seminar Nonlinear circuits and systems, volume, pages 5 59, Wien, 199. Springer-Verlag. [] L. Caprioli. On triode oscillator described by differential equation of third order. In rudy Mezhdunarod. Simpoziuma po Nelineinym Kolebanijam, pages 8 87, Kiev, Naukova Dumka. [3] O. B. Ershova and G. A. Leonov. Frequency estimates of the number of cycle slippings in the phase control systems. Autom. Remote Control, 44(5):6 67, [4] M. Imanaliev. Oscillation and stability of solutins of singularlyperturbed integro-differential systems. ILIM, Frunze, [5] P. Kokotovic, H.K. Khalil, and J. OReilly. Singular Perturbation Methods in Control. 1999, Philadelphia, PA, [6] G. A. Leonov, V. Reitmann, and V. B. Smirnova. Non-local methods for pendulum-like feedback systems. eubner, Stuttgart-Leipzig, 199. [7] G.A. Leonov. Phase synchronization. theory and applications. Autom. Remote Control, 67(1): , 6. [8] G.A. Leonov and N.V. Kuznetsov. Nonlinear analysis of phaselocked loop (PLL): Global stability analysis, hidden oscillations and simulation problems. In A.K. Belyaev et al. (eds.), Mechanics and Model-Based Control of Advanced Engineering Systems, pages 199 7, Wien, 14. Springer-Verlag. [9] W. Margaris. heory of the Non-linear Analog Phase Locked Loop. Springer, New York, 4. [1] E. Pak, E. Perieva, and V. Smirnova. Stability of singularly perturbed integral and integro-differential volterra equations. In Proceedings of XVI International Conference Dinamical Systems Modelling and Stability Investigation, pages 47 48, Kiev, Ukraine, 13. [11] A. Perkin, A. Proskurnikov, and V. Smirnova. Asymptotic estimates for gradient-like distributed parameter systems with periodic nonlinearities. In Proceedings of 14 IEEE Multi-conference on Systems and Control (MSC 14), pages , Antibes, France, 14.

7 [1] A. Perkin, V. Smirnova, A. Shepeljavyi, and N. Utina. he employment of periodic lyapunov functions for asymptotic analysis of multidimensional phase control systems. In Proceedings of 4-th Conference on Physics and Control (PhysCon-9), Catania, Italy, 9. [13] A.A. Perkin, V.B. Smirnova, and A.I. Shepeljavyi. Frequencyalgebraic conditions for stability of phase systems with application to phase-locked loops and synchronization systems. Cybernetics and Physics, 1(3): , 1. [14] V. Rasvan. Four lectures on stability. Lecture 3. he frequency domain criterion for absolute stability. Control Engineering and Applied Informatics, 8():13, 6. [15] V. B. Smirnova, N. V. Utina, and A. I. Shepeljavyi. Asymptotic frequency-domain estimates for the amplitude of output in discrete phase systems. Vestnik SPb Univ. Math., 39(1):45 51, 6. [16] J.J. Stoker. Nonlinear vibrations in mechanical and electrical systems. Interscience, New York, 195. [17] Y. Yang and L. Huang. Cycle slipping in phase synchronization systems. Phys. Lett. A, 36: , 7.