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

0387986995Basic TheoryC

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

328X. THE SECOND-ORDER DIFFERENTIAL EQUATION

(iii)the function Y(x) is a solution of the differential equation

 

F (x, Y,

I = 0,

and d'Yx(x)

M for a positive constant M on the interval 0 < x < t.

Let y(x) be any solution of the differential equation

 

 

AL2 + f (x,y,LY) = 0

(A>0)

satisfying the conditions y(O) = Y(0) and ly'(0) - Y'(0)l l< p. Then,

ly(x) - Y(x)I < {h, + a(L + k)}exp[ Lx]

on the interval 0 < x < if c and A are sufficiently small.

Proof.

We prove this theorem in four steps.

Step 1. Setting

y=u+Y(x) and

 

-__v+dY),

change the equation

 

 

 

 

ad.z + f (x, y,

d-T

) = 0

(A > 0)

 

 

 

 

to the systern

 

 

 

 

du

= v,

dv

 

 

(X.7.1)

 

F(X, u, v, a),

dx

 

where

 

 

 

 

F (x, Y(x), v + d

) I + F(x, u, v,.)l

< (e + AM) + Kf ul,

as long as (x, u, v) is in the regi/on

 

l

 

Do = {(x, u, v) : 0 < x < t, Jul <- a(x), lv) < pe-`= + b(x)}.

Also,

F (riY(x)v + d ()) > Lv

for

F x, Y(x). v + d()) < Lv

for

in Do. Hence, in Do,

v > 0,

v < 0,

F(x, u, v, A) < -Lv + Klul + (e + AM)

(X.7.2)

F(x, u, v, A) > -Lv - Klul - (E + AM)

for v > 0,

for v 0.

8. A SINGULAR PERTURBATION PROBLEM

329

Step 2. Suppose that two functions wl (x, .A) and w2(x, A) satisfy the following conditions:

(X.7.3)

10 < wl(x,.A) < a(x), 0 < w2(x,A) < pe-' + b(x),

wl(0,A) > 0,

w2(0,.1) > p,

wi (x, A) > W&, A), Aw2(x, A) > -Lw2(x, A) + Kwl (x, A) + (E + AM)

on the interval 0 < x < e. Then, as long as (x, u, v) is in the region

Dl = {(x,u,v): 0 <_ x < e, IuI < wi(x,A), IvI < w2(x,A)},

it holds that

IvI < wi (x, A),

.F(x, u, w2(x, A), A) < Aw2(x, )1),

F(x. u, -w2(x, A), .A) > -aw2(x, I\).

Look at the right-hand side of (X.7.1) on the boundary of V. Then, it can be easily seen that if a solution of (X.7.1) starts from Dl, it will stay in Dl on the interval 0 < x < e.

Step 3. Show that two functions

 

Jwi(x,A) = J0 x

A)l

+ 62(1+1),

 

 

 

 

 

w2(x A) =

pe-Lx/A + Alex:/L

 

 

 

 

ApA

e + AM + K62(e + 1)

, satisfy the requirements (X.7.3) if f, A,

where bl = L2 +

L

 

and a positive constant 62 are sufficiently small. Observe that two roots of AX2 +

LX -K=0 are -L -(o and (o = L +O(A).

Step 4. Note that

wl (x, A) = (L / (1 - e-Lx/A) + bK I (eKx/L - 1) + 82(x + 1).

To complete the proof, look at Jul < wl (x, A) as 62 -4 0. 0

The inequalityI

v< wr2( a,)as 62

0,

yields the following estimate of d

 

VI

 

dy

dY(x) < pe-Lx/A +

f E

+

(h +

M1 ] ex'/L.

dx

dx

L

L

L

J

Note that

lim a-Lx/A = 0

if x > 0.

330

X. THE SECOND-ORDER DIFFERENTIAL EQUATION

X-8. A singular perturbation problem

In this section, we look at behavior of solutions of the van der Pol equation

(X.4.1) as a --+ +oo more closely. Setting t = Er and A = c2, let us change (X.4.1) to

 

d2x

(x

2

-

1)dx

+ x = 0.

 

J1dr2 +

 

dr

 

 

 

 

 

 

Set A = 0. Then, (X.8.1) becomes

 

 

 

 

 

(X.8.2)

(x2 - 1) d + x = 0.

Solving (X.8.2) with an initial value x(O) = xo < -1, we obtain x = 0(r), where

2 )2 - In I0(r) I = -r + 2 - In(-xo).

Observe that

d0(r)

dr

0(r)

> 0 if d(r) < -1.

0(7)2-1

 

2

Note also that setting ro = 2 - In(-xo) -- 2 > 0, we obtain ,(ro) = -1 and

45(7-)2 - 1 > 0 for 0:5 r < ro. The graph of 0(r) is given in Figure 22.

FIGURE 22.

(i) Behavior for 0 < r < 7-o -bo, where do > 0: Let us denote by x(r, A) the solution to the initial-value problem

(X.8.3)

d2z

(x2 -

dx

x = 0,

x(O, A) = xo, 40,A) = 7-7,

 

a dr2 +

1) dr- +

 

 

where the prime is d7and q is a fixed constant. Using Theorem X-7-1 (Na6J, we

derive

Ix(r, A) - O(T)I 5 AK,

I Ix'(r, A) - 4'(r)l 5 I7-7 - 0'(0)le-P'1x + AK,

8. A SINGULAR PERTURBATION PROBLEM

331

for 0 < T < To - 60, where K and a are suitable positive constants.

(11) Behavior for lx+1l< o: First note that

lim p0'(T) _ +oo. Set ¢'(ro-2bo) _

1 > 0. Then, there exists T1(A) for sufficiently small A > 0 such that

Pi

 

 

 

 

{ O<ri(A)<ro_So,

llm T1(A) _ 7.0 - 280,

 

 

 

X

 

x(T1(A),A) _ O(Tp - 26p),

x'(T, A) >

for r1(A) < T < TO - 60.

 

 

 

 

2P1

Set p =

1A)

Then,

 

 

 

 

 

(X.8.4)

 

ALP

= p2(x2 - 1) + p3x,

regarding p = p(x, A) as a function of x and A. Note that

0 < p(x, A) < 2p'

for ¢(ro -26o) < x < x(ro-b0, A) < -1, 0 < A < AO,

where A0 is a sufficiently small positive number. It is important to notice that if we make 60 small and if we make A0 also small accordingly, we can make p1 small.

Set

(X.8.5)

o = -1 - ¢(T0 - 26o)

and M0 = max Ix2 - 11.

 

 

I1+xI<a

Then,

 

 

(X.8.6)

o > 0,

lien o = 0,

and

 

lim Mp = 0.

 

 

 

 

 

 

0-+0+

Furthermore, 0 < p1 < Mo since P1 = 1 - ((TO - 25p)2

Therefore,

 

 

 

 

0(TO - 260)

 

 

0 < p(z, A) < 2Mo for

 

¢(r0-25o) < x < x(ro - 50, A), 0 < A < A0.

Now, we shall show that

 

 

 

 

 

 

(X.8.7) 0 < p(x, A) < 2Mo

for

 

Ii + xt < o,

0 < A < A0.

If (X.8.7) is not true, there must exist a

such that

 

 

 

p({, A) = 2M0, 0 < p(x, A) < 2Mo for -1-a:5 x<

Then, this implies that

 

 

 

 

 

 

5

A)2(Mo +

A)2(1 + g) < 0.

332 X. THE SECOND-ORDER DIFFERENTIAL EQUATION

This is impossible.

(iii) Behavior for -1 + or < x < 2 - al: To begin with, it should be remarked that

 

[..-_]2 i

=

8

-2- ( -3+1J = 0,

 

J(2_1)d

 

 

 

 

(X.8.8)

 

-x)-3<0, for -1<x<2.

 

1)d. = C3

Fix a positive number al so that

 

 

 

 

(X.8.9)

j(2 - 1)4 < - y

for -1+a< x < 2-al,

where y is a sufficiently small positive constant. Note that ry --+ 0 as do

0+ and

al -0+-

 

so thatx,,\-

 

2

Integrating (X.8.4), it follows that

 

 

 

 

 

p(z ) = P(

I, A) - J

(C2

1) - f 1

agd

 

 

 

 

 

(I)P(

 

 

 

x<0. Hence, if we choose \0>0

A

7

0 < A < A0,

where

K = f

 

 

 

 

for

 

 

y

< 4K

 

 

 

0

 

w e obtain 0 < p(x, A) < 9K for -1+a <x<0, 0<A <A0.

a

(II) p(x A) > y -K max A) for 0 < x < 2 - al. Suppose that

0 <

for

0 < < x, 0 < A < A0.

 

2K

 

 

Then, 0 < p(x, A) < 2

< L. Thus, we proved that

 

 

2K

 

 

(X.8.10) 0 < p(x,A) < 2K for

-1+a < x < 2 - al,

0 < A < A0.

(iv) Behavior for 2 - al < x < 2 + a2: First look at

 

p(2 - at, A)

 

1) +

p(, )ld

2&0), A) -

 

8. A SINGULAR PERTURBATION PROBLEM

333

 

Notice that p(2 - al, A) is independent of 5o. Then, it can be shown that

 

A

 

0+

as a, -+ 0+

and A

0+.

 

 

 

 

p(2 - a,, A)

 

 

 

 

 

 

For a fixed positive number a2, assume that

 

 

 

0 < p(x, A) < 1

for

2 - al < x < 2 + a2i

0 < A < Ao.

 

Then,

 

 

 

 

 

 

 

 

 

2+02 g2 - 1) + P( A)fl > Q.

 

P(2 + 2, A)

p(2 -- o , A)

12-,,

 

 

 

Hence,

 

 

 

 

 

 

2+0a g2 - 1) +

 

<

A

for 2 - al <x < 2 + a2.

 

J-o,

 

p(2 - al, A)

 

 

 

This is a contradiction if al and A are small. Therefore, if Ao > 0 is sufficiently small, there exists an x(A) such that

2 - al < x(A) < 2 + a2

and

p(x(A),A) = 1

for

0 < A < Ao.

It can be shown that Urn x(A) = 2.

ao*

Setting r(x(A)) = r(A), it can be shown that liM +r(A) = ro. Now again, apply

Theorem X-7-1 (Na6) to the initial-value problem

AT2 + (x2 - 1)d

+ x = 0,

x(r(A)) = x(A), x'(r(A)) = 1.

 

 

Figure 23 shows the general behavior of x(r, A) for small positive A.

TO

FIGURE 23.

334

X. THE SECOND-ORDER DIFFERENTIAL EQUATION

 

 

 

 

EXERCISES X

 

 

X-1. (1) Show that if F(t,yl,y2) is continuous and bounded on

 

Il = {(t, Y1, Y2) a < t < b, [yl[ < +oo, 1y2[ < +oo},

the boundary-value problem

 

 

 

 

dt2

=

F (t, y, &)

!LY (a) = a,

y(b)

 

 

 

 

 

 

has a solution.

 

 

 

 

 

 

(2) Does the boundary-value problem

 

 

 

dtY

=

F

t, y, dY

dy(a)

= a'

dY

 

dt2

 

 

()dt

dt

dt (b) _ 3

have a solution?

 

 

 

 

 

 

(3) Show that the boundary-value problem

 

 

 

d2x

= tx + x3,

x(-2) = A,

x(3) = B

 

dt2

 

 

 

 

 

 

has a solution for any real numbers A and B.

Hint. A counterexample for (2) is

d2 = 0,

(n) = 0,

(b) = 1.

 

For (3), apply Theorem X-1-3 [Na4j with w,(t) = a and w2(t)

where -a

and 3 are sufficiently large positive constants.

X-2. Find the global phase portrait of each of the following two differential equa- tions:

(i)

d2x + x2(1 - x)2 dt +

(x - 1)2(x + 1)x = 0;

 

2

dx

-

 

 

 

d2x +

1

= 0.

 

dx2

dt

 

l + x2

 

X-3. Suppose that

(i)f (x, y, t) is continuously differentiable everywhere in the (x, y, t)-space (i.e., in p3),

(ii)g(x) is continuously differentiable in -oo < x < +oo,

(iii)e(t) is continuous on 0 < t < +oo,

(iv)f (x, y, t) 0 for x2 + y2 > r2 and t > 0, where r is a positive number r,

r>

 

(v) G(x) = J g(t)dC --' +oo as jx[

+oo,

0

 

(vi) 1 je(r)Idr is bounded for 0 < t < +oo.

0 `

Show that every solution (x(t), y(t)) of the system

dx _

dy =

-f(x,y,t)Y

- g(x) + e(t)

dt - Y '

 

df

 

 

 

is well defined and bounded for 0 < t < +oo.

EXERCISES X

335

Hint. Set V (x, y) = 2 + G(x). Then, d = - f (x, y, t)yz + e(t)y. There exists a

positive number ro > r such that V(x, y) > 4 for x2 + y2 > ro. Hence, if an orbit (x(t), y(t)) satisfies conditions that x(t)2 + y(t)2 > ro for to < t < ti, we obtain

d V(x(t),y(t)) < Ie(t)Iiy(t)I <2(e(t)j

V(x(t),y(t))

for

to <t <t,.

This yields

 

 

 

 

 

 

V(x(t),y(t)) <

 

V(x(to),y(to)) + ft., (e(r)!dr

for

to <t < t,.

X-4. Show that the differential equation

 

 

 

( E.1 )

d

yi

_

ill

 

 

dt

[y2 ]

[ -E(yi - 1)y2 - i/i - Eyi

 

 

 

 

 

has exactly one periodic orbit, where a is a positive parameter-

X-5. Find an approximation for the unique periodic orbit of (E.1) for small e > 0 and for large e.

X-6. Considering the system of two differential equations

(E.2)

d

yi

y2

dt

[ y2]

- [ -e(yi - 1)y2 + yi - byl

 

where a and 6 are positive numbers,

(a)show that every orbit is bounded for t > 0,

(b)determine whether each stationary point is stable or unstable, examining every

possibility,

(c) using the function

!y,

-y2 + E(-yi + 6YNY2 + J (-s + bs3J[1 + E2(si)]ds

0

and Theorem IX-2-2, show that if 0 < 6< 3, every orbit of system (E.2) tends to

one of the stationary points as t -' +oo,

(d) discuss the uniqueness of periodic orbits.

X-7. Show that the differential equation

d2x

d z + f (x) d+ 9(x) = 0

has at most one nontrivial periodic solution if fix) and g(x) are continuously differentiable in !R and satisfy the following conditions

 

 

< 0

for

1x3 < 1,

(i)

f (x) { > 0

for

IzI > 1,

 

9(x)

<0

for

-2<x<0,

(ii)

> 0

for (x + 2)x > 0,

 

 

(iii)

 

J

g(x)dx = 0.

336

X. THE SECOND-ORDER DIFFERENTIAL EQUATION

Hint. See [Sat].

The main ideas are as follows:

(1)If there are more than one nontrivial periodic orbits, then at least one of them should not be orbitally asymptotically stable.

(2)For a nontrivial periodic orbit (x(t), x'(t)) of period T > 0, the inequality fT

f (x(t))dt > 0 implies that this orbit is orbitally asymptotically stable

o

(cf. Exercise IX-5).

(3) Set

 

 

( x'(t))2

 

f

x(t)

 

a(t) = V(x(t), x'(t))=

+

 

 

2

 

 

Then,

da(t) _ -f(x(t))(x'(t))2.

 

 

 

 

dt

 

 

 

 

 

Therefore, we obtain

 

 

 

 

 

( I )

J f(x(t) )dt = -

r dA

 

 

 

 

 

J (x') 2

 

 

Now, investigate the quantity (x'(t))2 as a function of A along the periodic orbit

(x(t), x'(t)). For a fixed value of A, compare (x'(t))2 for different values of x(t), using the assumptions on f and g.

Step 1. First the following remarks are very important.

(a) If we set y = x', the given differential equation is reduced to the system

(S)

dt = y,

dt = -f (x)y -

A careful observation shows that the index of the critical point (0, 0) is 1, while the index of the critical point (-2,0) is -1. There are only two critical points of (S).

(c)The periodic orbit and any line {x = a constant) intersect each other at most twice.

(d)The critical point (-2,0) should not be contained in the domain bounded by the periodic orbit. To show this, use the fact that the index of any periodic orbit is 1.

These facts imply that the periodic orbit should be confined in the half-plane x >

-2.

2

r:

 

Step 2. A level curve of V (x, y) = 2 + J0

is a closed curve if the curve is

totally confined in the strip jxj < 1. In particular, (1, 0) and (-1, 0) are on the same

level curve V (x, y) = / 1

_ 1 g(C)d{. Since A(t) is increasing as long as

0

Io

jx(t)I < 1, the periodic orbit cannot intersect the level curve V(x,y) = / ,

0

Therefore, the periodic orbit must intersect the lines x = 1 (as well as the line x = -1) twice. Denote these four points by (1,s1(A)), (1,-rt(B)), (-1,-q(C)),

and (-1, 77(D)), where r)(A), n(B), q(C), and q(D) are positive numbers.

EXERCISES X

337

Step 3. Set

 

 

AA = V(1, rl(A)), AB = V (1,

Ac = V(-1, -n(C)), AD = V (-1, rl(D))

Then,

AC > AD,

AA > AD-

AA > AB, AC > AB,

Set

Z(A) = { A : AA > A > max(AB, AD) },

Z(B) = { A : AB < A < min(AA, )'C) },

Z(C) = { A : Ac > A > max(AB, AD) },

Z(D) = { A: AD < A < min(AA, Ac) }.

Note that the interval lo = JA: min(AA, Ac) > A > max(AB, AD)} is contained in

Z(A) n Z(B) n Z(C) n I(D).

Now,

(a)on the arc between (1, rl(A)) and (1, -rl(B)) of the periodic orbit, regard y as a function of A and denote it by yi (A), where AB < A 5 AA,

(Q) on the arc between (1, -rl(B)) and (-1, -r1(C)) of the periodic orbit, regard y as a function of A and denote it by y2 (A), where AB < A < AC,

(ry) on the are between (-1, -rl(C)) and (-1, rl(D)) of the periodic orbit, regard y as a function of A and denote it by yi (A), where AD < A < AC,

(6) on the arc between (-1, rl(D)) and (1, rl(A)) of the periodic orbit, regard y as a function of A and denote it by yz (A), where AD :5 A < AA.

Then, it can be shown that

yi (A)2 < y2 (A)2

on

I(A),

yi (A)2

< y2 (A)2

on

Z(B),

yi (A)2 < y2 (A)2

on

I(C),

yi (A)2

< y2 (A)2

on

Z(D).

Step 4. Now, fixing a AO E Z(A) n Z(B) n 1(C) n T(D), evaluate the integral (I) as follows:

rT

{JA8 dA

 

rc

dA

 

rAD

dA

AA

dA 1

J

f(s(t))dt = -

w yi (A)2

+ J

y2 (A)2 +'Ac yl (A)2 + IAa

y2 (A)2 T ,

where

 

 

 

 

 

r"

 

 

 

 

 

dA

 

 

dA

 

dA

 

 

 

J\A

yi (A)2

_

as yi (A)2

lao

yi (A)2,

 

 

 

rac

dA

ra°

dA

+

AC

dA

 

 

 

Joe y2 (A)2

=

Jaa y2 (A)2

Jib

 

 

 

 

 

 

 

 

 

I"

dA

 

 

ra° dA

 

rAc

dA

 

 

 

 

 

 

 

 

 

yi (A)2 = - JaD yi (A)2 - ao yi

(A)2.

ao ys (A)2

aD y2 (A) 2

Jao y2