Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

0387986995Basic TheoryC

.pdf
Скачиваний:
23
Добавлен:
10.06.2015
Размер:
12 Мб
Скачать
☆
b denotes the usual dot product. Since fl(t) is smooth, there exists a

288

IX. AUTONOMOUS SYSTEMS

Thus, the manifold M is transversal to the orbit C(fjo) at rp. Note that (t, t

At, ¢(0, 6).

Moreover, by continuity of the solutions of (IX.3.1) with respect to the initial conditions, there exist two positive numbers 61 and ro such that for i satisfying the condition

11 -401

there is a real number r(q such that 1r(q-)l < ro and p(r(g71,'7) E M. Therefore,

(IX.3.18) implies that

1pit+r(17),q - P-(t,il0){ <Koboe-°`

for t > 0.

Thus, C(ilo) is orbitally asymptotically stable as t --+ +oc.

Remark IX-3-5. If the entries of an R'-valued function f (t,1/) is continuously differentiable on the entire f.-space lRn and system (IX.3.1) has a periodic solution p ( t , i l o ) of period 1 such that f (p(t, i )) j4 0, there exist n - 1 vectors 4j(t) C Rn (j = 2, ... , n) such that

(i)the entries of vectors 9'J (t) (j = 2,... , n) are continuously differentiable and periodic in t of period 1,

(ii) n vectors dpig0), q"2(t), ...,

form an orthogonal system.

Proof.

If n = 2, it is easy to find 42(t). For n > 3, there is the recipe.

Note first that

dpdtilo}

- -

f(At, 40))

 

f(p1t,ilo))00. Set g1(t)=

 

 

 

 

f (pit, no)) f (pit, ilo))

where u"

constant vector i71 such that nl ill = 1 and that ili + qi (t) A 0 for 0 < t < 1.

Choosing an orthonormal set {ill.... , in I, where n, are constant vectors, define

Q'n(t) = T)n -

on(t)

(rJi + ji(t))

(1 = 2, ... , n),

 

+ai(t)

 

 

where ah(t) _'h - ji(t). 0

For this construction, see, for example, (U).

Example IX-3-6. The orthogonal system of Remark IX-3-5 is useful to study solutions of (IX.3.1) in a neighborhood of a periodic solution. To illustrate this application, let us look at a system

dydt = f

where f is an 1R3-valued function whose entries are continuously differentiable on the entire. y-space R3. Assume that (SI) has a periodic solution p(t, irjo) of period

I such that f(plt,r'b)) 54 0. In this case, there exist two vectors j2(t) E IR3 and

4fi(t) E L3 such that

3. ORBITAL STABILITY

289

(i) the entries of vectors 4"4(t) (j = 2,3) are continuously differentiable and peri-

odic in t of period 1,

 

 

 

 

 

 

 

(ii) three vectors d

Q , 4'2(t), and 43(t) form an orthogonal system.

Note that dp(ttrlo) = f (plt, no)).

 

 

 

 

 

 

For any fixed real non-negative number r, the set

 

 

 

P(T) = {l7-,m) + u142(r) + u2Q3(T)

(ul,u2) E R2}

is the plane which is perpendicular to the orbit

 

at piT,, o). Letting t, u1, and

u2 be three functions of r to be determined, set

 

 

 

 

(IX.3.19)

 

f(t) = p( r, no) + u1Q2(r) + u243(T)

 

From (IX.3.19) and the given system (S1), we derive

 

 

 

dt

{/,-/

 

du1

 

due _

 

 

dge(r) + u2 dQ3(r)

f(f(t))dT = f(P(r>io)) +

dT 42(r) +

dT 93(r)

+ u1 d-,

dT

This yields

 

y' dr

u1 -

 

dr

 

 

 

 

 

 

 

 

 

(IX.3.20)

 

 

dt

(

42(T) 1

 

 

d93(r)1

df11

 

 

 

 

W -

2(T) f(y(t))dr -

2(t)' dr l u1 - ((t)

dT J u2,

due

 

 

dt

((t)

dQ2(T)\

 

93(t)

dg3(r))

and

 

 

 

 

 

 

 

 

 

(IX.3.21)

= [f(PIT'770))

f(PIT,, )) + (fcvlr,no))

u1

 

 

 

d9d(T)) uz]

x [f(Ar, *lo}).f(N(t))} -',

 

 

 

 

 

 

+ (f(p1r,TIo))

 

where a" b denotes the usual dot product and we assumed that

q"1(r} = 1 (j =

2, 3).

 

 

 

 

 

 

 

 

 

Note that

 

 

 

 

 

 

 

 

 

AM)) = RAT, 7lo))+ul

(PIT,r, ))42(r)+ u2

f(YlT,r/o))93(T)+O(ItII +Iu2I)

 

 

89

 

 

09

 

 

 

Hence, from(IX.3.20) and (IX.3.21), we derive

dt

= 1 + 00U11 + Iu21)

77.

290 IX. AUTONOMOUS SYSTEMS

and

(IX.3.22)

dul

 

 

ul + [i(T).

(r))1 U2

42(r) '

l

 

dr

 

 

 

 

 

- (o) . _d_r) Jut -

(0) -

d9drr)) u2 + O(lul l + 1U21),

93(T) '

W

(FFIr, 1I0))42(r) 1 ul + [i(r).

r, 90))93(r) 11 U2

d7-

 

 

 

 

(i(t). ddrr)) ul -

(t(t). d (r) J U2 + O(lulI + 1U21)

Let Q(t) be the 3 x 3 matrix whose column vectors are (at, qo)), fi(t), and

,fi(t), i.e., Q(t) = [ f (p1t, rlo)) 6(t) q"3(t)! . The transformation w = Q(t)v changes the linear system

(S2)

dtr,

_

 

L(PIt,rl *6

to

 

 

 

 

dv'

=

0

/31(t)

32(t)

0

atl(t)

a12(t) 16,

 

 

0

a21(t)

a22(t)

Using these notations, write (IX.3.22) in the form

 

dul

= all(r)ul + a12(r)u2 + 91(T,ul,u2),

(IX.3.23)

d7-

 

 

1 dug

 

 

 

= a21(r)u1 + a22(r)u2 + 92(T,u1,u2),

 

dr

 

 

 

where

 

 

 

191(r,u1,u2) -

 

K(r)(lul-V11+1u2-V21)

(i = 1,2)

with K(r) > 0 and limK(r) = 0 for lull + lull < r and Ivll + 1v21 < r. Observe r-O

that the two multipliers of the linear system

du' =

[all(t)

a12(t)

U

j

a21(t)

a22(t)

 

are also multipliers of system (S2). Therefore, using system (IX.3.23), Theorem IX-3-4 can be proven. In general, we obtain more precise information concerning the behavior of solutions in a neighborhood of a periodic solution in this way.

The materials of this section are also found in 1CL, Chapter 13, §2, pp. 321-3271.

is a periodic
Lemma IX-4-4. If pit, qo) is bounded fort > 0 and if there

4. THE POINCARE-BENDIXSON THEOREM

291

IX-4. The Poincare-Bendixson theorem

In this section, we explain the structure of C+(i-) on the plane. Consider an

R2-valued function f (y) of y E R2 such that the entries of f (y7) are continuously differentiable on the entire g-plane R2. Denote again by pit, ill the unique solution

of the initial-value problem dg = f (y-), y(0) = i. The main result of this section

is the following theorem due to H. Poincare [Poll] and 1. Bendixson [Ben2].

Theorem IX-4-1. Suppose that the solution pit, qo) is bounded for t > 0 and that C+(ilo) contains only a finite number of stationary points. Then, there are the following three possibilities:

(i)C+(i-Io) is a periodic orbit,

(ii)C+(i o) consists of a stationary point,

(iii)C+(i)) consists of a finite number of stationary points and a set of orbits each of which tends to one of these stationary points as I ti tends to +oo.

To prove this theorem, we need some preparation.

Definition IX-4-2. A finite closed segment t of a straight line in R2 is called a transversal with respect to f if f (y1 0 0 at every point on t and if the vector f (y) is not parallel to t at every point on t (cf. Figure 3).

Observation IX-4-3. For every transversal t and every point i , the set tnC+( contains at most one point (cf. Figures 4-1 and 4-2).

I

PU,

FIGURE 3.

FIGURE 4-1.

FIGURE 4-2.

The following lemma is the main part of the proof of Theorem IX-4-1.

exists a point i7 E C+(7) such that C+ (11) contains a nonstationary point, then C+(i )

orbit.

Proof.

We prove this lemma in three steps. Since C+(qo) is invariant and rj1 E C+(7 ),

it follows that p1t,i7l) E C+(jo) for all t. Furthermore, C+(ijt) C C+(io) since

C+(ilo) is closed.

Step 1. Let rj be a nonstationary point on C+(t7 ). Also, let t be a transversal with respect to f that passes through I. Then, t n C+(10) = {77} since rj E G+(ijl) C

C+(no).

C+(7-) for

292 IX. AUTONOMOUS SYSTEMS

Step 2. Since tj E C+(r)l ), there exists a sequence {tk : k = 1,2.... } of real

numbers such that

lira

tk = +oo and p1tk, ill) E t (k = 1, 2, ... ). Note that

 

m+oo

p(tk, ill) E C+(ilo).

Therefore, p(tk, m) = #(k = 1,2.... ). This implies that

there exist two distinct real numbers rl and T2 such that it = p(rl, ill) = (r2, ill) and hence p(t, rl = p(t, p(ri , ili )) = p1t, p(r2i ill )). This, in turn, implies that p(t + T1, ill) = p(t + r2i ill) for t > 0. Therefore, the orbit C(ill) is periodic in t of period Irl - r21. Furthermore, C(ill) C C+(rjo). Note that there is no stationary point on C(qj ).

Step 3. Since C+(ilo) is connected, it follows that distance(C(ijl), C+(ilo)-C(ni)) = 0. Hence, if C(ijl) # L+ (10), there exists a sequence { k E C+(i o) : k = 1, 2.... }

such that tk

C(rj,) (k = 1, 2, ...) and

lim £k = { E C(ill ). Assume that

 

 

k-.+oo

there exists a transversal t such that E e and lk E t (k = 1, 2.... ). Note that l; E C(ill) C C+(ilo). Then, k = t c C(ill) (k = 1, 2, ... ). This is a contradiction.

Thus, it is concluded that C+(rlo) = C(iji).

Now, we complete the proof of Theorem IX-4-1 as follows.

Proof of Theorem IX-4-1-

If C+ (, ) does not contain any stationary points, then (i) follows (cf. Lemma

IX-4-4). If C+(no) consists of stationary points only, we obtain (ii), since C+(ilo) is connected. If f-+ (Q contains stationary and nonstationary points, then

any point it E C+(ilo) does not contain nonstationary points (cf. Lemma IX-4-4). This is true also for t < 0. Hence, (iii) follows.

Observation IX-4-5. In cases (i) and (iii), the set 1R2 - C+(ilo) is not connected.

Furthermore, if an orbit C(i)) is contained in C+(i o), two sides of the curve p1t, tl") belong to two different connected components of R2 - C+(t ). In fact, if we consider a simple Jordan curve C which intersects with the orbit C(tl at n transversally, then the curve p(t,,o) intersects with C in a neighborhood of two distinct points on C infinitely many times (cf. Figure 5).

Theorem IX-4-6. If C(ijo) n C+(ilo) 0 0, then C+(ijo) = C(rjo) and either no is a stationary point or the orbit C(ilo) is periodic.

Proof

In this case, C(ilo) C C+(ilo). If C+ (8o) contains nonstationary points, the orbit

C(i o) consists of nonstationary points. Choose a transversal a at ilo. Then, it can be shown that C(ilo) is periodic in a way similar to the proof of Lemma IX-4-4, since r1o E C+(tlo).

Observation IX-4-7. If C(ilo) n L+(no) = 0, then

lim distance(plt,ilo) and

,C+(Q) = 0.

 

t-+oo

 

It follows that if C+(ilo) consists of a stationary point

then

lim p(t, irjo)

If C+(ijo) is a periodic orbit, then C+(rjo) is called a limit

t+oo

 

 

 

cycle (cf. Figures 6-1 and 6-2).

5. INDICES OF JORDAN CURVES

293

FIGURE 5.

FIGURE 6-1.

FIGURE 6-2.

The materials of this section are also found in [CL, Chapter 16, §§1 and 2, pp.

391-3981 and (Har2, Chapter VII, §§4 and 5, pp. 151-1581. In [Hart[, the PoincareBendixson Theorem was proved without the uniqueness of solutions of initial-value problems.

IX-5. Indices of Jordan curves

In this section, we explain applications of index of a Jordan curve in the plane to the study of solutions of an autonomous system in 1R2. Consider again an R'-valued function f (y-) of y E 7,k2 whose entries are continuously differentiable on the entire

y-plane R. Denote also by pit, J) the unique solution of the initial-value problem

dg

= f (y1, yi(0)

dt

To begin with, let us introduce the concept of indices of Jordan curves. Let C be a Jordan curve y = il(s) (0 < s < 1) with the counterclockwise orientation (cf.

Figure 7). Assume that &(s)) # 0 for 0 < s < 1. Set u(s) = f ('l(s)) (0 <

If(s))1

s < 1). Then, il(s) is the unit vector in the direction of f'(il(s)). There exists a real-valued continuous function 0(s) defined on the interval 0 < s < I such that

a (s) _ Isin0(s)I'

Definition IX-5-1. The index of the Jordan curve C with respect to the vector

field f (y) is given by I f{C) = 0(1) - 8(0) .

2r,

This definition is independent of the choice of a parameterization rj(s) of C and the function 0(s). Let us denote by D the domain bounded by C (cf. Figure 7).

Observation IX-5-2. If the domain D is divided into two domains DI and D2 by a simple curve £, then 1 f{C) = l f{",) + If-(8D2) if f (rl) 0 on £ U C, where 8D? (j = 1, 2) denote the boundaries of domains DI and D2, respectively, as the portions of I f{8Di) and I j{d Dt) along £ canceled each other (cf. Figure 8).

Observation IX-5-3. If All) A 0 on C U D, then I f{C) = 0. To show this, divide

V into sufficiently small subdomains, use the fact that the vector field 19(s) has no change on a small subdomain, and apply Observation IX-5-2.

294

IX. AUTONOMOUS SYSTEMS

Observation IX-5-4. Assume that a point d E V is a stationary point, i.e., f(d) = 0. Assume also that f(,) 96 0 on C U V except at a. Then, If(C) depends only on aa. Therefore, we define the index of an isolated stationary point d with respect to f by Ilea) = I f(OV), where V is a neighborhood of as such that there are no stationary points in V other than a.

Observation IX-5-5. If D contains only a finite number of stationary points

N

dl, d2, ... , dN, then I j<C) = FI J(a"k). Proof of this result is similar to that of

k=1

Observation IX-5-3.

Observation IX-5-6. If f (y) # 0 on C and if f(,7(s)) is tangent to C at every point i(s) (0 < s < 1) of C, then I1{C) = 1 (cf. (CL, Chapter 16, §4, Theorem

4.31).

Proof.

It is evident that the index of C with respect to f and the index of C with respect

to the tangent vector to C are the same, i.e.,

If4C) = hn)(C)-

Assume without any loss of generalities that

# ( 0 )=# ( I ) ) = 0 and i (s) > 0 f o r 0 < s < 1,

where rig (s) (j = 1, 2) are the entries of the vector r"(s) (cf. Figure 9).

 

 

- __11

n2=0

 

 

(0.0)

 

FIGURE 7.

FIGURE 8.

FIGURE 9.

 

For 0<r<s<1,set

 

if

T = s,

V(T,s) =

if

T < s (mod 1),

 

1

(T, s) _ (0, 1).

Then, v(T, s) is continuous in (T, s) for 0 <_ r < s < 1. Hence, studying how 16(0, s) changes from s = 0 to s = 1 and how i6(T,1) changes from T = 0 to T = 1, we obtain Ij.(.)(C) = 1. 13

 

5. INDICES OF JORDAN CURVES

295

Observation IX-5-7. If C is a periodic orbit of the autonomous system

 

(IX.5.1)

di

f(+1+

 

dt

 

 

 

 

then 1 f-(C) = 1.

Observation IX-5-8. If C is a periodic orbit of autonomous system (IX.5.1), then the domain D contains at least a stationary point (cf. Observations IX-5-3 and

IX-5-7).

Observation IX-5-9. Suppose that f(it(s)) is tangent to C only at a finite number of points ij(rl ), #(7-2),. .. ,1t(TN) on C. Assume also that at each point it(rk), either

(1) p(t,it(rk)) V V for ItI < 5k, or

(2) p(t, it(rk)) E D for 0 < ItI < bk,

for some sufficiently small positive numbers 51, ... , SN. Let us call rj(rk) an exterior

(respectively interior) contact point in case (1) (respectively (2)) (cf. Figure 10).

Let us denote by E (respectively H) the total number of interior (respectively exterior) contact points among N points it(-r1), ij(r2), ... , #(T-N). Then E+H = N and

(IX.5.2)

1 f-(C) = E-H+2

 

2

Sketch of proof.

 

Let rj(r,) and fj(rk) be two consecutive contact points such that T. < rk. Then,

(i)if both of these two points are exterior contact points, then the tangent to C changes 7r in angle more than the vector field f does from the point #(T.) to the point fj(rk) (cf. Figure 11-1),

(ii)if both of these two points are interior contact points, then the vector field f changes it in angle more than the tangent to C does from the point #(T,) to the point it(rk) (cf. Figure 11-2),

(iii)if these two points are an exterior contact point and an interior contact point, then the tangent to C and the vector field f change in the same amount in angle (cf. Figures 12-1 and 12-2).

Since the total amount of change of the tangent to C in angle is 27r (cf. Observation

IX-5-6), we arrive at formula (IX.5.2). 0

q(

FIGURE 10.

FIGURE 11-1.

FIGURE 11-2.

296

IX. AUTONOMOUS SYSTEMS

Observation IX-5-10. Let d be an isolated stationary point. Assume that a neighborhood V of d is divided into a finite number of sectorial regions S1, S2,... , SN by a finite number of orbits C(rji ), C(il3),... , C(iv) in such a way that

(a)rli, Tt2, - , nN E 8V,

(Q)either p1r,rjk) E V fort > 0 and tends to 99 as t -+ +oo, or p5r,ilk) E V for t < 0 and tends to d as t - -cc (cf. Figures 13-1 and 13-2).

FIGURE 12-1.

FIGURE 12-2.

FIGURE 13-1. FIGURE 13-2.

Let us assume also that each sectorial region Sk satisfies one of the following four conditions:

(I)if it E Sk and {il - d] is sufficiently small, then C(r)) C Sk and pit, q-) tends to d as Itl +oo (cf. Figure 14),

(II)if i) E Sk and Jil - dl is sufficiently small, then pit, i)) E Sk only on a finite t-interval ak < t < Yk (cf. Figure 15),

FIGURE 14.

FIGURE 15.

(111) if it E Sk and ail - dl is sufficiently small, then p1t,i) tends to dd in Sk as t +oc, but plt, q") does not tend to dd in Sk as t -+ -oo (cf. Figure 16-1),

(IV) if ij E Sk and Jrj - al is sufficiently small, then p1t,71) tends to ad in Sk as t -. -oc, but plt, q+) does not tend to dd in Sk as t - +oo (cf. Figure 16-2).

The sectorial region Sk is said to be elliptic (respectively hyperbolic) in case (I)

(respectively (II)). In cases (III) and (IV), the sectorial region Sk is said to be parabolic. Let us denote by E (respectively H) the total number of elliptic sectorial regions (respectively hyperbolic sectorial regions) among St, ... , SN. Then, it is known that

E--H+2

(IX.5.3)

2

This result is similar to formula (1X.5.2) of Observation IX-5-9. For a proof in detail, see (Har2, Chapter VII, §9, pp. 166-172).

5. INDICES OF JORDAN CURVES

297

Observation IX-5-11. Suppose that the vector field f (y e) depends on a parameter e E 0 continuously, where 0 is a connected set. In this case, if f (if(s), e) # 0 at every point if(s) on the Jordan curve C for all e E A, then the index I f(.,()(C) is a constant independent of e. In fact, in this case, the integer I f( E)(C) depends on e continuously, and hence it is a constant.

Example IX-5-12. If an isolated stationary point d is a node, a spiral point, or a center, then E = H = 0. Hence It{a") = I.

Example IX-5-13. If an isolated stationary point d is a saddle point, then E = 0

and H = 4. Hence

-1 (cf. Figure 17).

FIGURE 16-1.

FIGURE 16-2.

FIGURE 17.

Example IX-5-14. More general cases of isolated stationary points a are shown

by Figures 18-1, 18-2, and 18-3. In fact,

1E = 0, H = 2 and hence II-(a") = 0

E = 2, H = 0 and hence It-(d) = 2

E = 1, H = 1 and hence l1-(d) = 1

in the case of Figure 18-1, in the case of Figure 18-2, in the case of Figure 18-3.

Example IX-5-15. Let p(z) be a polynomial in a complex variable z with complex coefficients. Set y' = t(zj and y2 = !3(zJ (i.e., z = yl + iy2) and regard the

differential equation dz = p(z) as a system of two differential equations

(IX.5.4)

ddt' _ 9(p(yl + iy2)],

Y2 = (p(yl + iy2)]

 

 

dt

on the y"-plane. It is easy to see that, if zo = 171 +ir12 is a zero of p(z) of multiplicity

m, then if#?) = m, where AM _

R(p(yl + iy2)]

and i) =

111

For example,

%P(YI + iy2)]

[112, .

if p(z) = iz2, system (IX.5.4) becomes

 

 

 

 

dtl

-2y1Y2,

ddt2 = y1 - y2

 

 

The point 0 is an isolated stationary point and II<O) = 2.

If we define a vector field f (y, e) by f (W, e) = [!1], then, 1- (0) = 2 for

e > 0 (cf. Observation IX-5-11).