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,
F(x, u, v, A) < -Lv + Klul + (e + AM)
(X.7.2)
F(x, u, v, A) > -Lv - Klul - (E + AM)
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
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
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
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.
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, |
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 |