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

0387986995Basic TheoryC

.pdf
Скачиваний:
23
Добавлен:
10.06.2015
Размер:
11.75 Mб
Скачать

298 IX. AUTONOMOUS SYSTEMS

Example IX-5-16. Assume that f ( =

fl (pl'Y2) satisfies the condition f (arty)

 

f2(y142)

= aP f (rte, where A is a real variable and p is an integer. Set y = r os0 I. Then,

the autonomous system

= f(it) can be written in the form

dr

 

rPFI (0),

rdO

rPF2(0),

dt =

 

dt

 

where

F1(0) = f 1(cos 8, sin o) cos 0 + f2(cos 0, sin 0) sin 0,

Sl F2(0) = -f1(cos0,sin6)sin0 + f2(cos0,sin0)cos0

(cf. Exercise VIII-3). If a ray to is defined by F2(0) = 0, then to is an orbit of the

system dff = f (r))). Thus, the entire y"-plane can be divided into sectorial regions

by those orbits QB determined by equation F2(0) = 0. For example, in the case of system (IX.5.5), we obtain Fl (0) = - sin 0 and F2(0) = cos 0. Observe that

 

>0

for

-r<0<0,

>o

for

-2 <0<2,

 

 

 

 

 

dr

= 0

for

0=-r and 0, and

=0

f o r

0 = -

and .

 

dt

<0

for 0<0<r

dt

 

2

2 '

 

< 0

for

2 < 0 < 31r

 

 

 

 

Hence, E = 2 and H = 0 (cf. Figure 19). Therefore, I f-{0) = 2 (cf. Example

IX-5-15).

FIGURE 18-1.

FIGURE 18-2.

FIGURE 18-3.

FIGURE 19.

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

398-402] and [Har2, Chapter VII, §§2 and 3, pp. 146-151, and §§8 and 9, pp.

161-1741. For a geometric and topological treatment of indices, see (Mil.

EXERCISES IX

IX-1. For each of the following three systems, using the given function V (x, y), show that all orbits are bounded as t -. +oo .

 

 

 

EXERCISES IX

299

(1)

_ -xy2 - 4y,

_ -yx2 + 3x,

V(x,y) = 3x2 +4Y2;

(2)

= y,

 

= -x3 - y,

V(x,y) = x4 +2Y2;

 

dy

= - (xs - 3x4 + 2x3 + 120x2 - 23x + 5) - (1 + x2)y,

 

= y,

 

V(x,y) =

+ J (ss - 3x4 + 2x3 + 120x2 - 23s + 5)ds.

 

2

 

0

 

 

 

 

 

IX-2. Show that every solution y(t) and its derivative L(t) of the differential

equation j2 + dt + y3 =0 tend to 0 as t -i +oo.

IX-3. Consider an autonomous system

(S1)

yj

__

fj(11j,112)

 

[y2]

 

[f2(yi,y2)

 

Wt

 

'

where f j and f2 are continuously differentiable on the entire (yj, y2)-plane. Assume that

(1)

fl(yl,y2)>0

for y2>0 and

-oo<yl<+oo,

fi(yi,y2) < 0

for y2<0 and

- oo < yl < +oo,

(2)fj(yi,1) > 1 and fi(yi, -1) < -1 for -oo < yl < +00,

(3)f2(yi,1) = 0 and f2(yi, -1) = 0 for -oo < y1 < +oo,

(4)

1

f2(11j,y2) <_

YI(112 - 1)

for

jy2j <

1

 

 

 

 

 

 

f2(Y1,Y2) ? Y1 (Y22 - 1)

for

jy21 <

1

 

 

 

 

and 0 < yi < +oo,

and - oo < yi < 0,

Find G+((n1,n2)) for (771,772) such that In2t < 1.

Hint. There are two possibilities:

Case 1. The solution (yi(t),y2(t)) of (Si) that satisfies the initial condition

(C)

yi(0) = 771

and

112(0) = 172

is bounded as t -+ +oo.

Case 2. The solution (y1(t),y2(t)) of (SI) that satisfies the initial condition (C) is unbounded as t +oo.

In Case 1, G+((nj, n2)) is either {(0, 0)} or a periodic orbit. In Case 2, G+((nj, 92))

is{(x,±1):-oo<x<+oo}.

300

IX. AUTONOMOUS SYSTEMS

Examples.

(a)Every orbit in Iy21 < 1 of the system with fl(yl,y2) = y2 and f2(yi,y2) _ yl(y2 - 1) is periodic.

(b)The stationary point (0,0) is a stable spiral point if fl(yl,y2) = y2(2 - sin(yly2)) and f2(yj,y2) = 2yi(y2 - 1).

(c)The stationary point (0, 0) is an unstable spiral point if ff(yi, y2) = y2(2 +

sin(yiy2)) and f2(yi,y2) = 2y, (y22 - 1).

Verify these statements by using the function V(yl, y2) = yi - ln(1 - y2).

IX-4. Let us consider a system

(52)

dt -

W= [y2 ] ,

f(U _ [f2(Y-)

where f, (y-) and f2(y-) are continuously differentiable with respect toy in a domain

Do C R2. Assume that

(i)system (S2) has a periodic orbit p-(t, ik) of period 1 which is contained in the domain Do,

(ii) RA t' v) 6,

(iii) the integral (1 ( i (pjt, qo)) +

(5t, o))) dt is negative.

,9y

0

Show that the periodic orbit p(t, jo) is orbitally asymptotically stable as t -' +00.

Hint. This is called Poincnre's criterion. Look at

dzV 49f

dt =

Then, I det Y(1)I < 1 for the fundamental matrix solution Y(t) of this system suc that Y(0) = 12 (cf. (4) of Remark IV-2-7). Hence, an eigenvalue p of Y(1) must satisfy the condition dpi < 1. (The other eigenvalue of Y(1) is 1.) Now, use Theorem

IX-3-4.

IX-5. Assume that two functions f (x) and g(x) are continuously differentiable in x E R. Assume also that the differential equation dt2 + f (x) drt + g(x) = 0 has

a nontrivial periodic solution x(t) of period I such that f f (x(t))dt > 0. Show

0

that (x(t), x'(t)) is an orbitally asymptotically stable orbit as t

+00 in the (x, x')

phase plane.

 

Hint. Use Exercise IX-4-

 

IX-6. Show that there exists a nontrivial periodic orbit of the system

dI7

y=

 

W)

=f (y),

Ly2J,

Lf2h (Y-)

 

 

 

where

(1)the entries of the 1R2-valued function f (y-) is continuously differentiable on the entire y-plane,

EXERCISES IX

301

(2) f(6)

and f(y)

ifil96 ,

 

 

(3)

of (O) = i,

eft (0) = 8, 2Z-2 (0) = -2 and

(p) = 1,

 

ay1

 

OW

ay1

 

 

(4)

Ivll+lliM

 

 

+oo(ft(Y1,Y2) + yr) and

1Y11+liM

+oo(f2(yt,N2) + yz) exist.

 

 

 

 

IX-7. Given that

 

 

 

 

 

 

 

 

Y1 - y2

 

y2

find the total number E of elliptic sectorial regions, the total number H of hyper- bolic sectorial regions, and the total number of parabolic sectorial regions in the

neighborhood of the isolated stationary point 0 of the system

= f (y, E) in the

 

 

dt

following two cases: (i) E = 2 and (ii) E _ 2. Also, find I y(0) for e 96 0.

IX-8. Let us consider a system

 

 

(S3)

dt = f (y),

 

where the entries of the l 3-valued functionf is continuously differentiable on the entire y-space R3. Assume that (S3) has a periodic orbit pit,ip) of period 1 such that f (r'(t, o)) 0 0. Assume also that the first variation of system (S3) with respect to the solution p(t,' o), i.e.,

d9 Of

(r(t,io))u,

dt

has three multipliers p1 = 1, p2i and p3 satisfying the condition: 1P21 < 1 and jp3j > 1, respectively. Construct the general orbits pit, r) of (53) such that distance(At,17),

C(s'7o)) tends to0ast-+too.

IX-9. Show that the differential equation d-t2 + (x + 3)(x + 2) drt + x(x + 1) = 0 does not have nontrivial periodic solutions.

Hint. Two stationary points are a node (0,0) and a saddle (-1,0). Furthermore,

x2

 

setting f(x1,x2) = I -(xt + 3)(x1 + 2)x2 - xt(xt + 1) J , we obtain divf(x1,x2) _

-(x1 + 3)(x1 + 2) < 0 if x1 > -1. Also, use index in §IX-5.

IX-10. For the system

 

= f(x,y) = x(1 - x2 - y2) - 3y,

= g(x,y) = y(1 - x2 - y2) + 3x,

dt

 

(1) find and classify all critical points,

 

(2) find

 

dt = f (x, y)

+ 9(x, y) 5

 

.-

302

IX. AUTONOMOUS SYSTEMS

for the function V (x, y) = x2 + y2,

(3)find the set S = {(x): dVdt = 01,

(4)examine if S is an invariant set,

(5)find the phase-portrait of orbits.

IX-11. For the system

dx f(x,y,x) = - x(1 - x2 _ y2)2 + 3xz + y, dt =

dt = 9(x,y,z) = - y(1 - x2 - y2)2 + 3yz - x,

dz

h(x,y,z) = - z - 3(x2 + y2),

dt =

(1)find all critical points and determine if they are asymptotically stable,

(2)find

dV

= f (x, y, x) 8V + g(x, y, z)

8V + h(x, y, x) 8V

dt

8y

8z

8x

for the function V(x, y, z) = x2 + y21+ z2,

(3)find thesetS=((x,y,z): dV =0},

(4)find the maximal invariant set M in S,

(5)find the phase portrait of orbits.

IX-12. Find the a phase portrait of orbits of the system

dt = y,

dt = y + (1 - x2)x2(4 - x2).

 

IX-13. Let f (x, y) and g(x, y) be real-valued, continuous, and continuously differentiable functions of two real variables (x, y) in an open, connected, and sim-

ply connected set D in the (x, y)-plane such that 8f (x, y) + Lf (x, y) j4 0 for all

(x, y) E D. Show that the system dxt = f (x, y), d = g(x, y) does not have any nontrivial periodic orbit that is contained entirely in D.

IX-14. Let p(z) be a polynomial in a complex variable z and deg p(z) > 0. Set x = 3t(z] and y = Q [z]. Verify the following statements.

(i) In the neighborhood of each of stationary points of system

(A)

d = R[p(z)],

dt = `3`[6 (x)],

 

 

there are no hyperbolic sectors. Also, system (A) does not have any isolated nontrivial periodic orbit.

(ii) In the neighborhood of each of stationary points of system

(B) dt = 3R[p(z)], L = -`3'[p(x)],

there are no elliptic sectors. Also, system (B) does not have any nontrivial periodic orbits.

1Rk and let T2 (j = 1, 2, ...) be analytic mappings of Do to l

EXERCISES IX

303

IX-15. Find the phase portrait of orbits of the system

d = x2 - y2 -3x+2,

d _ -2xy + 3y.

 

Hint. See Exercise IX-14 with p(z) = (z - 1)(z - 2).

 

 

dx

dy

IX-16. Find explicitly a two-dimensional system dt = f (x, y),

d t= 9(x, y) so

that it has exactly five stationary points and all of them are centers.

IX-17. Consider a system

(S4)

dy

= f (y1,

dt

 

 

where the entries of the R'-valued function f are analytic with respect toy in a domain Do C R'. Assume that

(i)system (S2) has a periodic orbit pit, ijo) of period 1 which is contained in the domain Do,

(ii)AP-Tt,no)) 0

(iii)for any open subset V of Do which contains the periodic orbits p(t, i"p), there

exists an open subset U of V which also contains p'(t, o) such that for any point i in U, the orbits p'(t, i of (S4) is contained in V and periodic in t.

Show that if U is sufficiently small, for any point i in U, we can fix a positive period of p1t, y) so that is bounded and analytic with respect to i in any simply

connected bounded open subset of U.

Hint. Apply the following observation.

Observation. Let Do be a connected, simply connected, open, and bounded set in

. Suppose that, for any pointy E Do, there exists a j such that T3 [yl = y, where j may depends on Y.

Then, there exists a jo such that Tjo [yj = y for all y E Do.

Proof.

Set

E j _ {y E D o : T j [ y 1 = Y - ) .

j = 1, 2, ... .

Then,

(1) E. is closed in Do,

+00

(2) Do = U E,,

3=1

(3) Do is of the second category in the sense of Baire.

Hence, for some jo, the set Ego contains a nonempty open subset (cf. Baire's Theorem). Since Tea is analytic, we obtain Tao [ 7 = y" for all y" E Do. O

For Baire's Theorem, see, for example, [Bar, pp. 91-921.

CHAPTER X

THE SECOND-ORDER DIFFERENTIAL

EQUATION dt2 + h(x) + g(x) = 0

In this chapter, we explain the basic results concerning the behavior of solutions of a system

112

[-h(yi)112 - 9(111)

as t -- +oo. In §X-2, using results given in §IX-2, we show the boundedness of solutions and apply these results to the van der Pol equation

(E)

x + E(x2 - 1)

d

+ x = 0

 

 

 

(cf. Example X-2-5). The boundedness of solutions and the instability of the unique stationary point imply that the van der Pol equation has a nontrivial periodic solution. This is a consequence of the Poincar&$endixson Theorem (cf. Theorem IX-4-1). In §X-3, we prove the uniqueness of periodic orbits in such a way that it can be applied to equation (E). In §X-4, we show that the absolute value of one of the two multipliers of the unique periodic solution of (E) is less than 1. The argument in §X-4 gives another proof of the uniqueness of periodic orbit of (E). In §X-5, we explain how to approximate the unique periodic solution of (E) in the case when a is positive and small. This is a typical problem of regular perturbations. In

§X-6, we explain how to locate the unique periodic solution of (E) geometrically as e - +oo. In §X-8, we explain how to find an approximation of the periodic solution of (E) analytically as a +oc. This is a typical problem of singular perturbations.

Concerning singular perturbations, we also explain a basic result due to M. Nagumo [Na6] in §X-7. In §X-1, we look at a boundary-value problem

y = F (t, y, d i ) , y(a) = o, y(b) = 8.

Using the Kneser Theorems (cf. Theorems 1II-2-4 and III-2-5), we show the exis- tence of solutions for this problem in the case when F(t, y, u) is bounded on the entire (y, u)-space. Also, we explain a basic theorem due to M. Nagumo [Na4] (cf.

Theorem X-1-3) which we can use in more general situations including singular perturbation problems (cf. [How]).

For more singular perturbation problems, see, for example, [Levi2], [LeL], (FL],

[HabL], [Si5], [How], [Wasl], and [O'M].

304

1. TWO-POINT BOUNDARY-VALUE PROBLEMS

305

X-1. Two-point boundary-value problems

In this section, first as an application of Theorems 111-2-4 and 111-2-5 (cf. [Kn]), we prove the following theorem concerning a boundary-value problem

(X.1.1)

d 2 = F t, y, dt) ,

y(a) = a, y(b) = Q

Theorem X-1-1. If the function F(t, yt, y2) is continuous and bounded on a region 11 = {(t, Y1, y2) : a < t < b, lyt I < +oo, Iy2I < +oo}, then problem (X.1.1) has a solution (or solutions).

Proof.

For any positive number K, the set Aa = {(a, or, y2) : Iy21 < K} is a compact and connected subset of ft We shall show that A0 satisfies Assumptions 1 and 2 of §111-2 for every positive number K. In fact, writing the second-order equation (X.1.1) as a system

(X.1.2)

dyt

_

dye

= F(t, yt, y2),

dt

- Y2,

dt

 

 

 

we derive

 

 

 

 

y1(t) = y1(a) + J y2(s)ds, a

y2(t) = y2(a) + JF(s,yi(s),y2(s))ds` .

Hence, if (a, y1(a), Y2 (a)) E A0, we obtain

Jy2(t)I < K + M(b - a),

ly1(t)I < lal + [K + M(b - a)](b - a),

where I F(t, yi, Y2)1 < Al on Q. Therefore, Ao satisfies Assumptions 1 and 2 of

§111-2. Thus, Theorem 111-2-5 implies that SS is also compact and connected for every c on the interval Te = {t : a < t < b}.

We shall prove that if K > 0 is sufficiently large, the set Sb contains two points

(771,772) and (t;t, (2) such that

(X.1.3)

n, < 0 < (1.

 

In fact, by using the Taylor series at t = a, write yt (b) in the form

 

yt(b) = a + y2(a)(b - a) + 2 dt (c) (b -

a)2,

 

 

where c is a certain point in the interval 10. Since I ddt2 (c)

M, the quantity

yI(b)I can be made as large as we wish by choosing Iy2(a)I sufficiently large. Thus, there are two points (nt, n2) and ((1, (2) in Sb such that (X.1.3) is satisfied.

Since the set Sb is compact and connected, there must be a point (Q, () in the set Sb. This implies the existence of a solution of problem (X.1.1). 0

306

X. THE SECOND-ORDER DIFFERENTIAL EQUATION

Example X-1-2. Theorem X-1-1 applies to the following two problems:

(X.1.4)

dt2 + sin y = 0,

y(a) = a,

y(b)

and

 

 

 

(X.1.5)

0,

y(a) = a,

y(b) = Q

 

dt2 + p-+-, =

 

 

However, Theorem X-1-1 does not apply to

 

(X.1.6)

d 22 + y = 0,

y(a) = a,

y(b) = Q.

For more general cases, the following theorem due to M. Nagumo (Na4] is useful.

Theorem X-1-3. Assume that

(i) a real-valued function f (t, x, y) and its derivatives az and

09 f

are continuous

in a region V = {(t, x, y) : (t, x) E A, -oo < y < +oo}, where 0 is a bounded and closed set in the (t, x) -space;

(ii) in the region D, the function f satisfies the condition

(I) 1f(t,x,y)1 : 0(IyI),

where 0(u) is a positive-valued function on the interval 0 < u < +oo such that

(II)

+°° udu

+00;

(iii)two real-valued functions wi (t) and w2(t) are twice continuously differentiable on an interval a < t < b and satisfy the conditions

wi(t) < w2(t) for a < t < b,

x<w2(t)}co,

Ao

= {(t x)-a<t<b wi

 

 

_(t) <

and

d2wi (t)

2

(IV)

d2dt2(t)

> f

< f

1 t, wi (t), d dt t) )

t, W2 (t), ddt t)) , for a < t <b;

(iv) two real numbers A and B satisfy the condition

(V)

wi(a) < A < w2(a),

and wi(b) < B < w2(b).

Then, the boundary-value problem

 

 

d 22 = f I t, x, ') ,

x(a) = A, x(b) = B,

has a solution x(t) such that (t, x(t)) E Ao for a < t < b, i.e.,

 

wi(t) < x(t) < w2(t),

for a < t < b.

Proof.

The main tools are the following two lemmas.

1. TWO-POINT BOUNDARY-VALUE PROBLEMS

307

Lemma X-1-4. Let x(t, to,rl) be the solution of the initial-value problem

d2x = f (t, x, d .) ,

x(to)

x'(to) = n,

where a < to < b, (to, l;) E Ao. Then, for any given positive number M, there exists a positive number a(M) such that l x'(t, to, l;, il) j < a(M) if

(X.1.7)

jr71 < M and (T, x(r, to, e, rl)) E Do for to < r < t or t:5 r < to.

Proof.

 

 

 

Letting L be a positive number such that

 

(X.1.8)

w2(t) - w1(t) < L

for

a < t < b,

choose a(M) > 0 for any given positive number M in such a way that a(M) > 111 and

(X.1.9)

a(M)

udu > L.

 

 

+!

q5(u)

 

Suppose that there exist ri and r2 such that to < Ti < r2 < t and that

x'(r1) = M < x'(r) < x'(r2) = a(M)

for r1 < r < T2,

where x(T) = x(r, to, t ,17). Then, since x'(r) > 0 for ri < r < r2, it follows that

xO(X'x

 

x/ (7-)

for

T1 < -r < -r2.

(T))) <

 

 

 

Hence,

 

 

 

 

a(M)

u/du

 

 

f

0u) -

 

 

 

 

 

This contradicts the choice of L by (X.1.9). Therefore, Lenuna X-1-4 is true for to < r < t. We can treat the case t < r < to similarly, since if we change t by -t, the differential equation x" = f (t, x, x') becomes x" = f (-t, x, -x'). 0

Lemma X-1-5. Set

7 = rim {a(ba)I wi(a)+e, -w2' (a)-eI,

where e is an arbitrarily fixed positive number. Let also x(t, c) be the solution to the initial-value problem

d2x = f (t, x,

)

,

x(a) = A, i (a) =