Mathematical Communications (997), 47 57 47 On behaviour of solutions of system of linear differential equations Božo Vrdoljak Abstract. The present paper deals with the existence, behaviour and approximation of some solutions of the nonhomogeneous system of linear differential equations with functional coefficients. The behaviour of solutions in the neighbourhood of an arbitrary or respective curve is considered. The qualitative analysis theory of differential equations and the topological retraction method are used. Key words: system of linear differential equations, qualitative analysis, behaviour of solutions, approximation of solutions, asymptotic solutions. Sažetak. O ponašanju rješenja sustava linearnih diferencijalnih jednadžbi. Rad se odnosi na egzistenciju, ponašanje i aproksimaciju nekih rješenja nehomogenog sistema linearnih diferencijalnih jednadžbi s funkcionalnim koeficijentima. Promatra se ponašanje rješenja u okolini proizvoljne ili odredjene krivulje. Koristi se teorija kvalitativne analize diferencijalnih jednadžbi i topološka metoda retrakcije. Ključne riječi: sustav linearnih diferencijalnih jednadžbi, kvalitativna analiza, ponašanje rješenja, aproksimacija rješenja, asimptotska rješenja. AMS subject classifications: 34C05 Received May 5, 997, Accepted June 9, 997. Introduction Let us consider a system of equations x = A (t) x f (t), () Part of this paper is presented at the Mathematical Colloquium in Osijek organized by Croatian Mathematical Society - Division Osijek, January 7, 997. Faculty of Civil Engineering, University of Split, Matice hrvatske 5, HR- 000 Split, Croatia, e-mail: vrdoljak@gradst.hr
48 B. Vrdoljak where t I = (a, ), a R; x R n, n ; the matrix-function A (t) = (a ij (t)) i,j=,...n is real and continuous on interval I; the vector-function f (t) = (f (t), f (t),..., f n (t)) T is real and continuous on interval I. Let Γ = (x, t) Ω : x = ϕ (t), where Ω = R n I, ϕ (t) = (ϕ (t), ϕ (t),..., ϕ n (t)), ϕ i C (I, R), be an arbitrary curve. The system () has been studied by various authors, e.g., [, 4, 6, 8]. In [8] the system () is considered with respect to the set σ = (x, t) Ω : x i ϕ i (t) < r i (t), i =,..., n, where r i C (I, R ), i =,..., n, R = (0, ). In this paper we shall consider the behaviour of integral curves of () with respect to the set ω = (x, t) Ω : αi (t) [x i ϕ i (t)] < (the neighbourhood of curve Γ), where α i C (I, R ), i =,..., n. The obtained results are new and independent of known results.. Preliminaries To obtain our results we need the following results concerning the applicability of the qualitative analysis theory and the topological retraction method of T. Wažewski (see, e.g., [,3,5,8,9]) We shall consider the behaviour of the integral curves (x (t), t), t I of the system () with respect to the set ω. The boundary surface of ω is W = (x, t) Ω : H (x, t) := αi (t) [x i ϕ i (t)] = 0. Let us denote the tangent vector field to an integral curve (x (t), t) of () by X. The vector H is the external normal on surface W. We have X (x, t) = a j (t) x j f (t),..., a nj (t) x j f n (t),, H (x, t) = j= j= α (t) [x ϕ (t)],..., αn (t) [x n ϕ n (t)], α i (t) [x i ϕ i (t)] αi (t) ϕ i (t) [x. i ϕ i (t)] By means of scalar product P = ( H, X) on W we shall establish the behaviour of integral curves of () with respect to the set ω. Let us designate with S n (I) a class of solutions defined on I which depends on n parameters. We shall simply say that the class of solutions S n (I) belongs to the
System of linear differential equations 49 set ω if graphs of functions in S n (I) are contained in ω. In that case we shall write S n (I) ω. The results of this paper are based on the following Lemmas. Lemma. If, for the system (), the scalar product P < 0 on W, then the system () has a class of solutions S n (I) belonging to the set ω for every t I, i.e., S n (I) ω. Proof. The condition P < 0 on W means that integral curves of the system () pass through points of the set W entering ω, i.e., the set W is a set of points of strict entrance of integral curves of () with respect to the sets Ω and ω (see, e.g., [3, 5]). Hence, all integral curves passing through points of the set ω remain in ω for all further values of variable t. Consequently, an n-parameter family of solutions of () S n (I) ω. Lemma. If, for the system (), the scalar product P > 0 on W, then the system () has at least one solution on I whose graph belongs to the set ω. Proof. The condition P > 0 on W means that integral curves of the system () pass through points of the set W exiting from ω, i.e., the set W is a set of points of strict exit of integral curves of () with respect to the sets Ω and ω (see, e.g., [3, 5]). Let Z be a set defined by the intersection of the set ω W and the plane t = t I, i.e., Z = (x, t) ω W : t = t. The set Z W is not a retract of Z, but is a retract of W. Hence, according to the retraction method, there is a point P Z\W, such that the integral curve of (), passing through P, remains in ω for every t I, i.e., K (P ) ω. Another intersection of ω W with plane t = t I, t t, should yield point P and K (P ) ω. However, those two integral curves can coincide, so we can only grant the existence of at least one solution of () on I whose graph belongs to the set ω. 3. Main Results Theorem. Let u C (I, R 0 ), R 0 = [0, ), such that a ij (t) ϕ j (t) f i (t) ϕ i (t) u (t), t I. () (i) If j=,j i j= αi (t) a ij (t) αj (t) a ji (t) u (t) (3) α j (t) < a ii (t) α i (t), i =,..., n, t I, then the system () has a class of solutions S n (I) ω. (ii) If j=,j i α α j (t) i (t) a ij (t) αj (t) a ji (t) u (t) (4)
50 B. Vrdoljak < a ii (t) α i (t), i =,..., n, t I, then the system () has at least one solution on I whose graph belongs to the set ω. Proof. For the scalar product P we have [ n ] P (x, t) = αi (t) [x i ϕ i (t)] a ij (t) x j f i (t) α i (t) [x i ϕ i (t)] αi (t) ϕ i (t) [x i ϕ i (t)] = [x i ϕ i (t)] a ij (t) [x j ϕ j (t)] j= a ij (t) ϕ j (t) f i (t) j= = = α i (t) [x i ϕ i (t)] αi (t) ϕ i (t) [x i ϕ i (t)] αi (t) a ii (t) [x i ϕ i (t)] n j=i [ α i (t) a ij (t) α j (t) a ji (t) ] [x i ϕ i (t)] [x j ϕ j (t)] [x i ϕ i (t)] a ij (t) ϕ j (t) f i (t) j= α i (t) [x i ϕ i (t)] αi (t) ϕ i (t) [x i ϕ i (t)] [ α i (t) a ii (t) α i (t) ] [x i ϕ i (t)] n j=i [ α i (t) a ij (t) α j (t) a ji (t) ] [x i ϕ i (t)] [x j ϕ j (t)] [x i ϕ i (t)] a ij (t) ϕ j (t) f i (t) ϕ i (t) j= = n [ a ii (t) α i (t) ] αi (t) [x i ϕ i (t)] j=i αi (t) a ij (t) αj (t) a ji (t) [x i ϕ i (t)] α j (t) [x j ϕ j (t)] α j (t)
System of linear differential equations 5 [x i ϕ i (t)] a ij (t) ϕ j (t) f i (t) ϕ i (t). If we introduce a notation A ij (t) = we have R (x, t) = P (x, t) = j= X i (x, t) = [x i ϕ i (t)], A ii (t) = a ii (t) α i (t), i =,..., n, (5) [ α α j (t) i (t) a ij (t) αj (t) a ji (t) ], i j (A ij = A ji ), (6) X i (x, t) a ij (t) ϕ j (t) f i (t) ϕ i (t), j= n A ii (t) Xi (x, t) j=i A ij (t) X i (x, t) X j (x, t) R (x, t). Now, according to the conditions (), (3) and (4) and the inequality ab a, b R, the following estimates for P on W are valid: (i) ( a b ) P (x, t) R (x, t) = < n A ii (t) Xi (x, t) A ii (t) j=,j i j=i A ij (t) [ X i (x, t) X j (x, t) ] A ij (t) Xi (x, t) R (x, t) [ u (t)] Xi (x, t) R (x, t) = u (t) R (x, t) 0, (ii) P (x, t) R (x, t) = > n A ii (t) Xi (x, t) A ii (t) j=,j i j=i A ij (t) [ X i (x, t) X j (x, t) ] A ij (t) Xi (x, t) R (x, t) u (t) Xi (x, t) R (x, t) = u (t) R (x, t) 0.
5 B. Vrdoljak The statements obtained for the function P, according to the given Lemmas, grant the presented statement. Theorem. Let u, v i C (I, R 0 ), i =,..., n, such that a ij (t) ϕ j (t) f i (t) ϕ i (t) v i (t), vi (t) u (t), (7) j= i =,..., n. (i) Under condition (3) the system () has a class of solutions S n (I) ω. (ii) Under condition (4) the system () has at least one solution on I whose graph is contained in the set ω. Proof. Using the proof of Theorem it is sufficient to note that in view of (7), R (x, t) u (t) on W. Indeed, using the concept of the relative extremum of the function R (x, t) on W, we have R (x, t) a ij (t) ϕ j (t) f i (t) ϕ i (t) X i (x, t) 4. Applications j= ( n v i (t) X i (x, t) vi (t) ) u (t). Let us note that conditions of Theorems and are simplified if curve Γ is some particular curve and if we take the function α (t) = (α (t),..., α n (t)) in an appropriate special form, e.g., = β (t) or = q = const, i =,..., n. Theorem 3. Let Γ be the curve of stationary points of () and let u C (I, R 0 ) such that [ ϕ i (t)] u (t), t I. (8) (i) If (3) holds true, the system () has a class of solutions S n (I) ω. (ii) If (4) holds true, the system () has at least one solution on I whose graph belongs to the set ω. Proof. It is sufficient to note that the condition (8) follows from (7) when Γ is the curve of stationary points of the system (). Example. The system of equations x = ( t 3) x t 3 x sin t sin t, (9) x = ( t 3 sin t ) x t 3 x t 3 ( sin t ) sin t has a two-parameter family of solutions x (t) = (x (t), x (t)) satisfying the condition [x (t) sin t] [x (t) ] < t 4 t >. (0)
System of linear differential equations 53 Here we consider the behaviour of solutions of (9) with respect to the set ω = (x, x, t) R 3 : (x sin t) (x ) < t 4, t >, i.e., in δ neighbourhood, δ (t) = t, of the curve of stationary points (sin t,, t), t >, of the system (9). Here we have α (t) = ( t, t ) and u (t) = t. Theorem 4. Let Γ, be an integral curve of the system (). (i) If j=,j i α α j (t) i (t) a ij (t) αj (t) a ji (t) () < a ii (t) α i (t), i =,..., n, t I, then the system () has a class of solutions S n (I) ω. (ii) If j=,j i α i (t) a ij (t) αj (t) a ji (t) () α j (t) < a ii (t) α i (t), i =,..., n, t I, then the system () has at least one solution on I whose graph belongs to the set ω. Proof. Here is a ij (t) ϕ j (t) f i (t) ϕ i (t) 0 on I, i =,..., n, j= and the proof of this Theorem follows from the proof of Theorem for u (t) = 0. Theorem 5. Let Γ be an integral curve of the system () and let terms A ij (t) be defined by (5) and (6). (i) If A (t)... A k (t) ( ) k det. A k (t)... A kk (t) > 0, k =,..., n, t I, (3) then the system () has a class of solutions S n (I) ω. (ii) If A (t)... A k (t) det.. A k (t)... A kk (t) > 0, k =,..., n, t I, (4) then the system () has at least one solution on I whose graph is contained in the set ω.
54 B. Vrdoljak Proof. In this case the function of a scalar product P on the corresponding set W is P (x, t) = = n A ii (t) Xi (x, t) j=i A ij (t) X i (x, t) X j (x, t), i,j= A ij (t) X i (x, t) X j (x, t) i.e., P is a quadratic symmetric form. According to the Sylvester criterion condition (3) grants that P < 0, and condition (4) grants that P > 0 on the corresponding set W. This, according to given Lemmas, confirms the statements of the Theorem 5. Now, let us consider the system () in which coefficients a ij (t) 0 for i j >, i.e., system x = a (t) x a (t) x f (t), x i = a i,i (t) x i a ii (t) x i a i,i (t) x i f i (t), i =,..., n, (5) x n = a n,n (t) x n a nn (t) x n f n (t), where coefficients satisfy the condition (for every i,..., n) or a i,i (t) a i,i (t) < 0, t I (6) a i,i (t) = a i,i (t) = 0, t I. (7) Let us define positive functions η k, k =,..., n, η (t) =, ηk (t) = ηk (t) a k,k (t) a k,k (t), k =,..., n. (8) If condition (7) is valid for a certain i, then we should take η i (t) =. Corollary. Let µ C (I, R ) and = µ (t) η i (t), i =,..., n, where functions η i are defined by (8) and let Γ be an integral curve of the system (5). (i) If a ii (t) α i (t) < 0, i =,..., n, t I, then the system (5) has a class of solutions S n (I) ω. (ii) If a ii (t) α i (t) > 0, i =,..., n, t I, then the system (5) has at least one solution on I whose graph belongs to the set ω.
System of linear differential equations 55 Proof. It is sufficient to note that α k (t) a k,k (t) α k (t) a k,k (t) = 0, k =,..., n, t I and conditions of the Theorem 4 are evidently satisfied. Example. The system of equations x = x t sin t g (t) x f (t), x = g (t) x t x t x 3 f (t), (9) x 3 = 4 (t ) x t sin t x 3 f 3 (t), where g, f i C (I, R), has a tree-parameter family of solutions x(t)=(x (t), x (t), x 3 (t)) satisfying the condition ( 4 ) [x (t) ϕ (t)] [x (t) ϕ (t)] [x 3 (t) ϕ 3 (t)] t ( < C ) t > t where ϕ (t) = (ϕ (t), ϕ (t), ϕ 3 (t)) is an arbitrary solution of the system (9), C > 0 is an arbitrary constant. In the case of the system (9) we have t η (t) = η (t) =, η 3 (t) = (t ) and µ (t) = ( ). C t Remark. The obtained results also contain an answer to the question on approximation of solutions x (t) whose existence is established. The errors of the approximation are defined by the function δ (t) = α(t) which can be arbitrarily small t I. In case when the function δ (t) tends to zero as t we have asymptotic solutions of the corresponding system of differential equations. In Example we have δ (t) = ( t, t ) and, for a sufficiently large t, each curve K ω can be taken as an approximate integral curve of the system (5), for example x (t) = (sin t, ) is an approximate solution and the inequality (0) presents an estimate of the accuracy of this approximate ( solution. ( ) ( ) ( C In Example we have δ (t) = t, C t, C t) ), where C > 0 can be arbitrarily small. Remark. The obtained results give the possibility that we can talk about the stability of the system (). For example, in view of the Theorem 4, in case (i) the system () is stable, and in case (ii) the system () is unstable according to Ljapunov, only if the conditions () and () are valid for α (t) = q = const. The system (9) is stable, and the system (9) is asymptotically stable. Remark 3. Let us have a quasilinear system x = A (x, t) x f (x, t), (0)
56 B. Vrdoljak where a matrix-function A (x, t) and a vector-function f (x, t) satisfy the conditions sufficient for the existence and uniqueness of solutions which are determined by initial data in Ω. The obtained results for the system () are also valid for the system (0) only if A (x, t) and f (x, t) satisfy the corresponding conditions on ω. Example 3. Consider the system of Volterra type x = [p i (x, t) f i (x, t)] x i, i =,..., n, () where p i, f i C (Ω, R). Corollary. Let p R 0, q R and ω 3 = (x, t) Ω : x i < q e pt. If or p i (x, t) f i (x, t) < p 0, i =,..., n on ω 3 () p i (x, t) < p 0, x i f i (x, t) 0, i =,..., n on ω 3, (3) then the system () has S n (I) ω 3 and the system () is exponential asymptotically stable. Proof. Here we consider the behaviour of solutions of () with respect to the set ω 3. The system () is exponential asymptotically stable because the conditions () and (3) are valid with the function qe pt for every q > 0. Takeuchi [7] considers the stability of the solution x = 0 of the system References x = [p i f i (x)] x i, i =,..., n, p i R. [] B. P. Demidowitsch, Lectures on Mathematical Theory of Stability (in Russian), Moscow, 967. [] J. Diblík, On existence of δ-bounded solutions of a nonhomogeneous linear system of differential equations, Funkcialaj Ekvacioj 34(99), 8. [3] J. Diblík, On asymptotic behaviour of solutions of certain classes of ordinary differential equations, J. Differential Equations 95(99), 03 7. [4] W. A. Harris, D. A. Lutz, Recent results in the asymptotic integration of linear differential system, International Conference on Differential Equations, Academic Press (975), 345 349. [5] P. Hartman, Ordinary Differential Equations, Wiley, New York, 98. [6] K. J. Palmer, Exponential dichotomy, integral separation and diagonalizability of linear systems of ordinary differential equations, J. Differential Equations 43(98), 84 03.
System of linear differential equations 57 [7] Y. Takeuchi, Existence of stability equilibrium point for dynamical systems of Volterra type, J. Math. Anal. Appl. 79(98), 4 6. [8] B. Vrdoljak, On parameter classes of solutions for system of linear differential equations, Glasnik Mat. 0(40)(985), 6 69. [9] T. Wažewski, Sur un principe topologique de l examen de l allure asymptotique des intégrales des équations différentielles ordinaires, Ann. Soc. Polon. Math. 0 (947), 79 33.