On the stability of solutions of a system of differential equations

Similar documents
(Received March 30, 1953)

On the convergence of solutions of the non-linear differential equation

Complex Variables. Chapter 2. Analytic Functions Section Harmonic Functions Proofs of Theorems. March 19, 2017

Global Attractivity in a Nonlinear Difference Equation and Applications to a Biological Model

Leplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math EECE

Limit. Chapter Introduction

JUST THE MATHS UNIT NUMBER ORDINARY DIFFERENTIAL EQUATIONS 3 (First order equations (C)) A.J.Hobson

Leplace s Equations. Analyzing the Analyticity of Analytic Analysis DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING. Engineering Math 16.

Math 4263 Homework Set 1

("-1/' .. f/ L) I LOCAL BOUNDEDNESS OF NONLINEAR, MONOTONE OPERA TORS. R. T. Rockafellar. MICHIGAN MATHEMATICAL vol. 16 (1969) pp.

LEBESGUE INTEGRATION. Introduction

Lecture Notes on Metric Spaces

MATH 205C: STATIONARY PHASE LEMMA

Math 10C - Fall Final Exam

1 Definition of the Riemann integral

Variational Methods & Optimal Control

(7: ) minimize f (x) subject to g(x) E C, (P) (also see Pietrzykowski [5]). We shall establish an equivalence between the notion of

Math 140A - Fall Final Exam

2.10 Saddles, Nodes, Foci and Centers

There is a more global concept that is related to this circle of ideas that we discuss somewhat informally. Namely, a region R R n with a (smooth)

Global Attractivity of a Higher-Order Nonlinear Difference Equation

LECTURE 10: REVIEW OF POWER SERIES. 1. Motivation

DIFFERENTIAL EQUATIONS

The local equicontinuity of a maximal monotone operator

Math 5630: Iterative Methods for Systems of Equations Hung Phan, UMass Lowell March 22, 2018

Calculus. Mathematics FRDIS. Mendel university in Brno

Functions of Several Variables

Math 201 Handout 1. Rich Schwartz. September 6, 2006

MATH34032 Mid-term Test 10.00am 10.50am, 26th March 2010 Answer all six question [20% of the total mark for this course]

PRELIMINARY THEORY LINEAR EQUATIONS

Available online at Adv. Fixed Point Theory, 3 (2013), No. 4, ISSN:

arxiv: v1 [math.oc] 21 Mar 2015

CHAPTER 3. Analytic Functions. Dr. Pulak Sahoo

DIFFERENTIAL EQUATIONS

1 Implicit Differentiation

First order Partial Differential equations

Minimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality

1. Introduction. 2. Outlines

THE NEARLY ADDITIVE MAPS

Chapter 3 Second Order Linear Equations

Ordinary Differential Equations

William P. Thurston. The Geometry and Topology of Three-Manifolds

1 The Existence and Uniqueness Theorem for First-Order Differential Equations

2.12: Derivatives of Exp/Log (cont d) and 2.15: Antiderivatives and Initial Value Problems

FINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.

NONLINEAR PERTURBATION OF LINEAR PROGRAMS*

LINEAR ALGEBRA REVIEW

MATH 409 Advanced Calculus I Lecture 12: Uniform continuity. Exponential functions.

Sec. 14.3: Partial Derivatives. All of the following are ways of representing the derivative. y dx

Tangential vector bundle and Todd canonical systems of an algebraic variety

Permanency and Asymptotic Behavior of The Generalized Lotka-Volterra Food Chain System

MORE ON CONTINUOUS FUNCTIONS AND SETS

NEUTRIX CALCULUS I NEUTRICES AND DISTRIBUTIONS 1) J. G. VAN DER CORPUT. (Communicated at the meeting of January 30, 1960)

THEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)

Into a Hilbert Space

EVANESCENT SOLUTIONS FOR LINEAR ORDINARY DIFFERENTIAL EQUATIONS

EC9A0: Pre-sessional Advanced Mathematics Course. Lecture Notes: Unconstrained Optimisation By Pablo F. Beker 1

DIFFERENTIABILITY VIA DIRECTIONAL DERIVATIVES

(x x 0 ) 2 + (y y 0 ) 2 = ε 2, (2.11)

DIFFERENTIAL EQUATIONS

2. Second-order Linear Ordinary Differential Equations

MATH 452. SAMPLE 3 SOLUTIONS May 3, (10 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic.

A TWO PARAMETERS AMBROSETTI PRODI PROBLEM*

Math Exam III - Spring

arxiv: v2 [math.mg] 14 Apr 2014

MEAN VALUE METHODS IN ITERATION

Bessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics

(1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min.

Theory of Ordinary Differential Equations

Chain Rule. MATH 311, Calculus III. J. Robert Buchanan. Spring Department of Mathematics

COMMON FIXED POINT THEOREM OF THREE MAPPINGS IN COMPLETE METRIC SPACE

Math 5210, Definitions and Theorems on Metric Spaces

ON AN INFINITE SET OF NON-LINEAR DIFFERENTIAL EQUATIONS By J. B. McLEOD (Oxford)

Chapter 2: Differentiation

Math 222 Spring 2013 Exam 3 Review Problem Answers

(a) If om = [MI is a directed family in X, 9i = f(m) ], M E An, is a directed family

i c Robert C. Gunning

A NEIGHBORHOOD TREATMENT OF GENERAL TOPOLOGICAL SPACES*

Polynomials as Generators of Minimal Clones

Lesson 3: Linear differential equations of the first order Solve each of the following differential equations by two methods.

Lecture Introduction

. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,

s P = f(ξ n )(x i x i 1 ). i=1

Series Solution of Linear Ordinary Differential Equations

ON A DIFFERENCE EQUATION WITH MIN-MAX RESPONSE

Maths 212: Homework Solutions

Asymptotic Behavior of a Higher-Order Recursive Sequence

Sample Solutions of Assignment 10 for MAT3270B

3. Identify and find the general solution of each of the following first order differential equations.

MAT 570 REAL ANALYSIS LECTURE NOTES. Contents. 1. Sets Functions Countability Axiom of choice Equivalence relations 9

The Implicit Function and Inverse Function Theorems

On pseudomonotone variational inequalities

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck

Definition 5.1. A vector field v on a manifold M is map M T M such that for all x M, v(x) T x M.

MATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck

Asymptotic behaviour of solutions of third order nonlinear differential equations

The Mean Value Theorem. Oct

Mathematical Computing

Smarandachely Precontinuous maps. and Preopen Sets in Topological Vector Spaces

Lecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018

Transcription:

MEMOIRS O F TH E COLLEGE O F SCIENCE, UNIVERSITY OF KYOTO, SERIES A Vol. XXIX, Mathematics No. 1, 1955. On the stability of solutions of a system of differential equations By T aro YOSIIIZAWA (Received March 31, 1954) Introduction. About the nonlinear differential equation as well as about the linear one, various authors have discussed the bounded ness, the convergence" of solutions and the existence of a periodic solution ( 2 ) w h ic h w e have also discussed recently. Moreover we have obtained necessary and sufficient conditions for the boundedness of solutions of a certain type.' In these researches we have utilized the existence of a function of points which is characterised by the fact that it is nonincreasing o r nondecreasing along any solution of the given differential equation. Now by this idea, we will discuss the stability of solutions about which Liapounoff( 4 and various mathematicians and physicists have discussed actively and they have obtained many remarkable results in connection with the above mentioned problems. Now we consider a system of differential equations, (1) where y denotes a n ndimensional vector and F(x, y ) is a given vector field which is defined and continuous in the domain 0 < x < co, 1Y1< R (I.Y1='VY,'+ +Y:). Before the research of a solution, w e assume at first that y(x)=70 is a solution of (1), so th a t F(x, 0) is identically zero. A nd let y=y(x ; y o, x ) be any solution of ( (1) satisfying the initial condition Y=y0 when x Here we state the following definitions: ' 5 a) T he solutio n is said to be stable if, given any &> 0, there exists a a> 0 such that if ly o l 8, then ly(x ; y o, 0)1 & for all x > O.

28 T aro Yoshizawa b ) The solution y = 0 is said to be asymptotically stable if it is stable and there exists a (fixed) 80 > 0 such that if IY01 80, then lim y(x; yo, 0)=0. In this paper we will see the relations between these notions and the strong stability introduced by Okamure and one more, say equistability for the sake of convenience as compared with th e above defined stability. 1. Stability. At first, we will obtain a necessary and sufficient condition in order that is stable. Theorem 1. In order that the solution y 0 o f (1 ) is stable, it is necessary and sufficient that there exists a function 0 (x, y ) defined in 4 satisfy ing the following conditions ; namely 1 0 (x, y)> 0, provided ly 2 G (x, 0) =0 f or all x, 3 0 there ex ists a K () f or any (R 0 such that 60 (x, K (8) (> 0), whenever 0 x <co and 1Y1= 4 0 P(0 y ) is continuous at y=0, 5 for any solution of (1 ), y = y (x), the function 0 (x, y(x)) is a nonincreasing function of 1. P roof. A t first, we show that the condition is sufficient. Now we assume that there exists such a function 0(x, y). By 3 there is a K (6 ) for any 8 > 0 and by 4, we can choose a suitable positive number 8( < 8) for this K(8) such that (2) 0(0, y) < K (8 ), if IY1_ since 0(0, y ) is continuous at y = 0 a n d 0 (x, 0 )= 0. Now for a solution of (1) y=y(x ; y 0, 0 ) such as INI 8, we suppose that we have at some x, say x ly(x, ; 310, 0 )1 = 8. Thus we have by 3 and hence we have (X i, y (X1; Yo, 0)) K (6) (P(0, y(0 ; yo, 0))_/C(8) since 0(x, y(x)) is nonincreasing for x b y 5. O n the other hand we have (P (0, y(0 ; y0, 0)) G lc (8) b y (2 ) a n d therefore there arises a contradiction. Hence any solution of (1) y= y(x ; y(,, 0) such as y 0 ô satisfies for x < op

On the stability of solutions of a system of differential equations 29 I.Y(x ; Yo, 0) I < 8, that is to say, the solution y = 0 is stable. Nextly w e show th a t the condition is necessary. N ow w e suppose that y= 0 is sta b le. Here we consider the solutions of (1) going to the right from a point P in d and we represent this fam ily by Dp. W e will define a(p ) a s follows ; namely if there exists a solution of?),. which arrives a t ly =R, then we take a (P) =R and if otherwise, i.e. if all the solutions of (1) passing through P lie in the interior of 4, a( P ) = s u p IY(x) and when P lies on ly I a(p) = R. Then if P lies on the xaxis, we have a(p) 0, since for any point P on the xaxis the solution through P is unique to the right and it is y==_o b y the condition that 3170 is sta b le. W hen P does not lie on the xaxis, we have (P) > 0. Thus for any 6 (< R ) and P o n ly = E, w e have 8(P) >._.&whenever 0 <x < c o. This &is K (&) in the condition 3. Moreover if fo r any 8> 0 we choose a suitable a> 0, w e have for y o such a s Iy o < a and x 0 IY(x ; Y0, 0 ) I, since the solution y== 0 is stable. Therefore we have 0_<8(P) 8 fo r P lying in x=0, 3. Clearly if P and Q lie on the same solution of (1) and x, <4, w e have Here we put w ith P= (x, y) evi, ( P ) a (Q) (x, y) (P), and then it is clear from the above m entioned that this 0(x, y) satisfies the conditions 1, 2, 3, 4 and 5. Therefore this 0(x, y ) is the desired. 2. Equistability. In 1 w e have considered the case where y = 0 is stable for solutions starting from a neighborhood of y =0 at x = 0. Compared with this stability, here we consider the following case.

30 T aro Yoshizawa Definition. The solution y 0 is said to be equistable if, given any E>0, there exists a a(x 3 ) > O for every x such that if I y l < a(x0), then ly(x ; Yo, x9)1 for all x_> _ x.. Namely for every x is stable when we consider 31 0 as a solution going to th e right starting from x=x o. Often the equistability is said merely stable,'" specially in physical problems. A s Theorem 3 below shows, the stability and the equistability are equivalent if y===.0 is the unique solution starting from any point on the xaxis. If we can choose a independent of x, we shall say it uniformly equistable. Then we have in the same way as in the preceding theorem of the stability. Theorem 2. In order that the solution y = 0 o f (1 ) is equistable, it is necessary an d sufficient that there exists a similar function (A, y ) as o n e in Theorem 1, but only the condition 4 being replaced by 4 ' (x y ) is continuous at y=0 f o r every x o. Fo r th e uniform equistability, (x, y) is continuous at y = 0 with regard to y uniformly f o r x. When the solution y=. 0 is equistable, then it is of course stable ; b u t there is a case where it is stable and yet it is not equistable. For example, let us consider such an equation, (3) d y f(x, y), where 2(x1)y 1) 2(x1) (0 < x < 1, y) f (x, y)= 231 11 2 (0 < x 1, 0 < y (x 1) 2 ) 2 ( y) (0 _< x < 1, (x 1) '<y ( 0) 2(x1) (0 x<1, y.< (x1)=). Clearly y 0 is stable, being considered as the solution starting from the point x such as 0 x < 1, but it is not stable for A such as 1. The solution going to the left from x such as 1 on the xaxis is not unique in this case. T h is fact characterises the nonequistability. About this matter we have Theorem 3. W hen the solution 5, =_ 0 o f (1 ) is unique to the

On the stability of solutions of a system of diff erential equations 31 left for all points on the xaxis in Cauchyproblem, a if y 0 i s stable, ( ) it is also equistable. Namely in this case the stability and the equistability are the same. P roof. Suppose that y = 0 is not equistable. Then we assume that, even if we take any neighborhood of y = 0 at x=x, for some x o (> 0) and some E, there exists a solution going to the right from there such a s ly & at an x (> x(). For this & and a suitable 8> O we have ly (x ; y 0, 0)1<& if IY015 (3, since y===.0 is sta b le. By the uniqueness to th e left of ye=0, if we choose a suitable neighborhood U(x, 0) of y = 0 at x=x o, all the solutions of (1) going to the left from U arrive at the plane x=0 and moreover we have 13,1. at x = 0. Then there exists a solution starting from I I <a at x=0 such a s ly (x ; y o, 0)1 E and there arises a contradiction. Therefore y,= 0 is equistable. In the case where the righthand side of the system (1) does not depend on x, i.e. where (4) dy _ F ( y ), we can prove easily the following theorem, where we assume that F(0)=0. Theorem 4. If the solution is stable for the sy stem (4), it is also equistable (uniformly). The proof is omitted. 3. Strong stability. Okamura has introduced a notion, i.e. the strong stability. Definition. L et Y (x ) be an arbitrary function such that Y (0 )= 0, I Y(x) I < R (0 x< X ) which is of bounded variation in 0 x <X (X being arbitrary, 0 < X < 0 0 ). The solution y o f (1 ) is said to be strongly stable when for any given > 0 a n d fo r a suitable positive number E, we have I Y (X )I < 2, Provided, V [Y (x) F (x, Y (x))d x]< & ( q ) 0Șr;S X holds, where V denotes the total variation. Okamura has obtained the following theorems fo r th e strong stability. Theorem 5. A function 40(y) being defined in a neighborhood of y=0, if 59(y) is positive or zero according as I Iis positive or zero

32 T aro Yoshizawa and v (y ) satisfies the Lipschitz condition 14'00 so (30 I Klyi y21 and finally it is a nonincreasing function of x along any solution y (x) of (1), then the solution y=.0 of (1 ) is strongly stable. Theorem 6. In order that the solution y=0 of the system (4) is strongly stable, it is necessary and sufficient that there exists such a function F (y) as in Theorem 5 for F of the righthand side of (4). Moreover for an equation of the first order d Y f(y ) ( f ( 0 ) = 0 ), in order that y=0 is strongly stable, it is necessary and sufficient that, given any 6 > 0, both the measure of the set of y such as [0 <y < f ( y ) 0 ] and that of the set of y such as [ 6 <y < O, f (y )_ 0 ] are positive. Hence for the equation, (5) d y y sin' (y 0) 0 (Y=0), the solution y 0 is stable, but it is not strongly stable. Besides in (5) y =0 is uniformly equistable by Theorem 4. Hence there is a case where is uniformly equistable, but it is not strongly stable. By considering 90(y) in Theorems 5 and 6 as y ) in Theorem 2, we can see in these cases that if y =0 is strongly stable, then it is also uniformly equistable. And yet we have always Theorem 7. I f the solution y O o f (1 ) is strongly stable, it is also uniformly equistable. P roof. For every x o 0 ), given any )7 > 0 and for a suitable 6> 0, let Y (x ) be a function such that Y (x )= 0 for 0 < x < x o and it is a solution of (1) satisfying < 8 as x >x o + 0 for x e < x < X. Then by the strong stability of y = 0, we have I <ij, since for this Y(x) we have V [Y(x) IF(x, Y(x))] <6. O z X Therefore 0 is uniformly equistable. Thus the equation (5) shows that the strong stability is stronger than the uniform equistability, a fortiori, than the stability. Therefore in the case where y O is equistable simply (not uniformly),

On the stability of solutions o f a system of differential equations 33 y==.0 can not be strongly stable. For example, in d y y 2 e ( y 1 Y y=_0 is equistable, but it is not uniformly ; hence it shall not be strongly stable. ) Finally we add a theorem" due to Okamura showing a relation between the strong stability and the asymptotical stability. Theorem 8. If the solution of the system (4) is asymptotically stable, y O is strongly stable. REFERENCES (1) Yoshizawa ; " O n conv ergence o f solutions of the non linear differential equation," Mem. Coll. S c i. tyo Univ., A V ol. 28, Mathematics, No. 2 (1953), pp. 143151. (2) Yoshizawa ; " On the nonlinear deercntial equation," ibid., pp. 133 141 and "N ote on the existence theorem at' a periodic solution of the nonlinear differential equation," ibid., pp. 153 159. (3) Yoshizawa ; "N ote on the boundedness of solutions of a system of differential equations," these Memoirs, Vol. 28, No. 3 (1953), pp. 293298. (4) Liapounoff ; "Problèm e général de la stabilité du mouvement," Annals of Math. Studies, No. 17. (5) Massera ; " On L iapounoff's conditions of stability," Ann. of Math., Vol. 50 (1949), pp. 705 721. (6) Okamura ; " On t h e strong stability of the stationary Point in the current," Functional Equations (in Japanese), Vol. 40 (1943), pp. 619. (7) Lefschetz ; "L ectures on doerential equations," Annals of Math. Studies, No. 14. (8) Hayashi and Yoshizawa ; " New treatise of solutions of a system al ordinary differential equations an d its application to the uniqueness theorems," these Memoirs, Vol. 26 (1951), pp. 225 233. (9) Okamura ; " S u r une sorte de distance relative b système differentiel," Proc. of PhysicoMath. Soc. of Japan, 3rd Ser., Vol. 25 (1943), pp. 514523. (10) Okamura ; pp. 1617 in (6).