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

книги из ГПНТБ / Васильева А.Б. Асимптотические разложения решений сингулярно возмущенных уравнений

.pdf
Скачиваний:
20
Добавлен:
24.10.2023
Размер:
8.96 Mб
Скачать

190

КРАЕВЫЕ

ЗАДАЧИ

 

[ГЛ. 4

поводу f), имеет

вид ( ^ = 1 / 4 ,

t2 =

3/4)

 

 

 

t/2

при

0

1/4,

 

Уоі (0

1/4 —1/2

при

1/4 < * <

3/4,

 

 

t/2 — \/2

при 3 / 4 < * <

1,

(4.353)

 

—t

при

0

1/4,

 

 

У02 (0

—1/4

при

1 / 4 < і < 3 / 4 ,

 

 

г! —1

при

3 / 4 < / < 1 .

 

Используя выражение для точного решения задачи (4.351), (4.352), которое нетрудно найти, можно непосредственно

убедиться

в

том, что

«/-компонента решения

при р,—*-0

стремится

к

ломаной

(4.353)

 

*

 

lim y(t,

p) = ï U 0

п р и 0 < / < 1.

(4.354)

 

Ц->-0

 

 

 

 

У п р а ж н е н и е .

Найти точное

решение задачи"(4.351), (4.352)

ипроверить равенства (4.353), (4.354).

4.Переходы между корнями одинаковой степени устой­ чивости. В этом пункте будет показано, что явление пе-

рехода с корня ср на корень ф может иметь место не <*) 1

только в том случае, когда фѴотвечает большее число собственных значений Я(- с отрицательной действительной частью, а также и тогда, когда число собственных зна­ чений Я,- с положительной действительной частью и число

собственных

значений Я,- с

отрицательной

действитель-

ной

частью,

(/)

(к)

т. е. корни

отвечающих ф и ф, совпадают,

(/")

<*)

 

 

 

Фи ф имеют одинаковую степень устойчивости.

Этот случай будет детально исследован для системы,

эквивалентной одному

уравнению

второго

порядка

Р2

&=Р(г,

t),

 

 

 

 

 

 

^

= F(zt,

t),

P ^

= Zi

( 0 < * < 1 ) ,

(4.355)

с

краевыми

условиями

 

 

 

 

 

 

г2 (0,

р) = 0,

z 2 ( l , [І) = 0.

(4.356)

Для изучения указанного явления рассмотрим сначала автономную систему, получающуюся из (4.355), если в

 

 

ВНУТРЕННИЙ

ПОГРАНИЧНЫЙ

слой

191

правой

части

(4.355)

положить

t

равным

постоянной Т

и сделать замену t = тр:

 

 

 

 

 

 

 

 

~

= F(z2,

T) =

F(z2),

 

 

*£ = гх

(4.357)

(это присоединенная

система (см.

главу

1),

отвечающая

(4.355)

при t = T). Пусть

функция

F (z2)

имеет три про­

стых

нуля г2 = ф,-

( і = 1 ,

2,

3),

причем

Z^., (ф,) > О

Рис. 5. График функции F (z2 ).

( / = 1 ,

3), a Fz%

2 ) < 0

(рис. 5).

Введем фазовую пло­

скость

(Zj, z2 ). Из (4.357),

исключая

т, имеем

 

 

 

2

гл

 

 

 

 

 

 

 

 

-J-

= J

F(z)dz

+ c =

Q>(z2) + c,

 

откуда

 

 

 

 

 

 

 

 

 

 

 

 

г1=±Ѵг2~ѴФ(гі)

+ с.

(4.358)

Из этого

выражения

непосредственно

ясно, что в зависи­

мости

от

характера

графика Ф (z2)

семейство

фазовых

кривых

(4.358)

выглядит

по-разному.

Можно

выделить

три случая, для каждого из которых на рис. 6 пред­

ставлен график

Ф (z2) и соответствующее

семейство фазо­

вых траекторий:

 

 

 

 

Фг

Фг

 

 

а)

J F{z)dz>\

 

F{z)dz,

 

 

Фі

Фз

 

 

т . е . площадь / больше площади

II, или

Ф ( ф х ) < Ф ( ф 3 ) ;

 

Фг

Фг

 

 

б)

J F(z)dz<

J

F(z)dz,

 

 

Фі

Фз

 

 

т. е. площадь / меньше площади /7, или Ф ( ф і ) > Ф ( ф 3 ) ;

192

КРАЕВЫЕ ЗАДАЧИ

[ГЛ. 4

Рис.

6.

1 — график Ф (г2 ).

Точки

Ах с координатами 1

=

0, г2 = фх )

и А3

с

координатами (Zi =

0, г2 =

ф3 )—седла; точка Аг

с

координа­

тами

1

= 0, г2 = ф2 ) — центр. 2—сепаратрисы;

3—одна

из замкнутых

траекторий, окружающих

центр.

Стрелки

указывают

направление

возрастания т.

 

 

 

 

 

В Н У Т Р Е Н Н И Й

П О Г Р А Н И Ч Н Ы Й СЛО Й

 

 

193

 

 

 

 

 

 

Ф 2

 

 

 

 

фг

 

 

 

 

 

 

 

 

в)

 

 

 

 

J F{z)dz=

 

J

F(z)dz,

 

 

 

(4.359)

т. е. площадь

/

равна

площади

/ / ,

или

Ф(фх ) = Ф(ср3).

В этом случае

 

седла

Ах

 

и

А3

(рис. 6, б))

оказываются

соединенными

сепаратрисами, т. е.

 

имеет

место

так на­

зываемая

ячейка.

 

 

 

 

 

 

 

 

 

 

 

 

 

Обратимся теперь

к

(4.355).

 

 

 

 

 

 

 

 

I .

Пусть

функция

F (z2,

t) имеет

 

непрерывные

частные

производные

до

 

второго

 

порядка

 

включительно

в

некото­

рой

области

D

 

на плоскости

(z2,

t).

 

 

 

 

 

I

I .

Пусть

уравнение

 

F (z2,

t) = 0 имеет

три корня

2 2 = Фі(0

( ' = 1 .

2, 3)

такие,

что:

 

 

 

 

 

 

1)

Фі (0 <

Ф2

(t)<

Фз (0

при

0 < / < 1 ,

 

 

 

 

2)

область

{(z2, t)\qt

 

(t)

s?; z2

<

ф3

(t), О <

t <

1 } с D,

3)

 

^ ( Ф . - ( 0 ) > 0

 

 

np«

0 < * < 1

( / = 1 ,

 

3),

 

 

 

^ Л Ф Л О Х О

 

 

при

 

0 < f < l .

 

 

 

Из условия

I I

следует,

что

график

/''(z,,,

^)

при

каждом

фиксированном

 

t € [0,

1]

имеет

вид

такой,

как

на рис. 5.

С изменением t этот график может как-то деформироваться,

в частности, может меняться соотношение между

площа­

дями

In

II.

Оказывается, явление

перехода

решения

за­

дачи

(4.325), (4.326)

с корня

ф 3 ( 0

на

корень ^(t)

 

(оба

эти

корня

условно

устойчивы:

%li2 (t)

= ±

V FZi

(ф,- (t)),

i = l , 3) может наблюдаться, если при некотором

значе­

нии

t = T

площади I

м I I становятся одинаковыми, т. е.

на фазовой

плоскости присоединенной

системы, отвечаю­

щей

значению

t = T,

имеет

место

ячейка

(случай

в)).

Образно

говоря,

для

перехода

с

седла на

седло

между

ними должен быть как бы

«мост» в

виде

сепаратрисы

(см.

рис.

6,

б)).

 

 

 

 

 

 

 

 

 

К этому выводу можно прийти, пользуясь асимпто­

тическими формулами § 14.

Покажем,

что

существует

решение такого типа, как изображенное

на рис. 7, с точ­

кой перехода Т (так условно

назовем значение

t,

принад­

лежащее (0, 1), при котором

г2 = ф2 ),

которой

соответ­

ствует

рис. 6, б). Для этого поступим следующим

образом.

Считая,

что Т—некоторое,

не

известное

пока

значение

7 А, Б. Васильева, В. Ф. Бутузов

194

 

 

 

 

 

 

КРАЕВЫЕ ЗАДАЧИ

 

 

 

 

[ГЛ. 4

между

t = 0 и t — 1, построим, пользуясь алгоритмом § 14,

асимптотическое

представление

решения

системы

(4.355)

на отрезке [0, Т],

удовлетворяющее условиям z2

(0, р,) — 0,

z2 (7\

 

р) = ф2 (Т),

 

и стремящееся

при р,—>0

к

(0,

Ф3 (0)-

Другими

словами, построим

асимптотику решения, отве­

чающего

левому

 

полуциклу

на рис. 7. Все условия су­

 

 

 

 

 

 

 

 

 

ществования такого решения,

 

 

 

 

 

 

 

 

 

указанные в § 14, за исключе­

 

 

 

 

 

 

 

 

 

нием

пока

условия

V,

вы­

 

 

 

 

 

 

 

 

 

полнены

в

силу

I и

I I . Об

 

 

 

 

 

 

 

 

 

условии

V

будет

сказано

 

 

 

 

 

 

 

 

 

ниже *).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что мы имеем не

 

 

 

 

 

 

 

 

 

совсем тот

случай, который

 

 

 

 

 

 

 

 

 

подробно описан в § 14, так

 

 

 

 

 

 

 

 

 

как здесь

на концах задается

 

 

 

 

 

 

 

 

 

одна и та же компонента z2.

Рис.

7.

1,

2, 3—графики

г2

=

Но

 

теория

справедлива

и

= ф,- (t)(i=

1, 2,

3),

/ — график

для

такого

случая

(см. п.

г2-компоненты решения задачи

11

§

14).

 

 

 

 

 

 

(4.355),

(4.356),

Г —точка пе­

 

Точно

так

же

построим

 

 

 

рехода.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на отрезке [T, 1] асимптотику

решения

системы

(4.355),

удовлетворяющего

условиям

22(7",

р) = ф2 (7'),

z 2 ( l , р) = 0, и стремящегося

п р и р , - ^ 0

к (0,

фх

(t)),

т. е. решения,

отвечающего правому полу­

циклу

на рис. 7. Условия существования такого решения,

кроме пока условия V, также выполнены в силу

I и I I .

Приравнивая

в точке Т выражения для ги

отвечающие

верхнему

и нижнему

полуциклам, получим уравнение для

*)

Собственные векторы

Ьх

(t) =

(

fy1^]

и b2

(t) = (

fy1

^ )

 

 

 

(0I

 

 

у

J

\

b12 (t)J

 

 

\

b22

(t)J

( i = l ,

3)

 

1 т

 

 

 

 

 

 

матрицы

 

 

 

,

соответствующие собственным

значениям кх (t) и Я2

(t),

обладают

тем свойством,

что bkJ- (t)

ф 0

при

0 < * < 1 ,

k,j=l,2.

 

Поэтому

условия {Det Вп

(0) ф 0,

Det ß 2 2

(0) Ф 0, содержащиеся в V, а) § 14, и аналогичные

условия,

содержащиеся в V, б) § 14, которые сводятся в нашем случае к усло-

 

(3)

 

(3)

 

 

 

(3)

 

 

(3)

 

 

 

 

 

 

виям

Ъ,х (0) ф 0, Ь22

(0) ф 0,

bu (Г) Ф 0, Ь22

(Г) Ф О, будут выпол­

нены. Такое же замечание относится к правому полуциклу ( і = 1 ) . Поэтому ниже будут рассмотрены только условия, связанные с 5+ и 5~.

ВНУТРЕННИЙ ПОГРАНИЧНЫЙ СЛОЙ

195

определения того значения Т, которое действительно бу­ дет точкой перехода, так как в ней кривая, отвечающая левому полуциклу, гладко (z1 с точностью до множителя р

dzA

совпадает с J сопрягается с кривой, отвечающей пра­ вому полуциклу. Заметим, что Т является, вообще говоря, функцией р, и его следует искать в виде Т = Т0 + р7 \ + +... + \ькТк + . . . . Величина Т0 определится из условия равенства членов нулевого приближения для гх. Мы огра­ ничимся рассмотрением только нулевого приближения, поэтому индекс 0 у Т0 будем опускать 0 = Т).

Решение, отвечающее левому полуциклу с точностью О(р), дается согласно п. 5 § 14 следующими формулами (верхний индекс у пограничных функций означает номер отрезка, на которые разбит [0, 1]):

zx(t,

d)

(D

 

 

р) = П 0 г 1 ( т 0 ) + С 0 г 1 ( т г ) + О ( р ) ,

 

 

d)

d)

(4.360)

z2 (t,

р) = фз (t) + П 0 г 2

0 ) + Q0z2

Г )4-0 (р),

(D d )

где пограничные функции П0г(- (т0 ) и Q0z,- (тг ) определяются из следующих уравнений и дополнительных условий:

 

n U 2

(0) = - Ф з (0), г Ц (оо) = 0

 

( t = l ,

2),

dXT

= F (Фз ( Г ) + Q z 2

,

'

T),

^ = Q z

x

( Т г = Ц ^ < 0

* V f » ^ ' 1 ^0~*'

 

"

àlj

vu .

v

.

^

 

 

 

 

 

 

 

 

 

 

 

 

(4.362)

(D

 

 

 

(D

 

 

 

 

 

2). (4.363)

Q0z2

(0) = Ф2

( Г ) - Ф з ( Т ) ,

Q 0 Z / ( - ° ° )

= 0

( t = l ,

С другой стороны, асимптотика решения системы (4.355), отвечающего правому полуциклу, имеет вид

(2)(2)

zx{t,

p) =

n 0

z 1 ( x R ) + Q 0 z 1

( T 1 )

+

O(p),

(4.364)

 

 

 

(2)

 

(2)

z2 (t,

p) =

фх (t) + П 0 2 2 ( T 7 ) + Q 0

Z 2 ( T X ) + О

(p),

где

 

 

 

 

 

 

 

^ 1

= £ ( ф 1

( Л + П 0 2 2 ,

Г),

^

= П Л

( т г > 0 ) , (4.365)

196

 

 

КРАЕВЫЕ ЗАДАЧИ

 

 

[ГЛ. 4

І.і)

 

 

(2)

 

 

 

 

 

П о 2г (0) = Ф 2

( 7 , ) - ф г ( Л ,

1102,(00) = 0 ( і = 1,2),

(4.366)

( 2 )

 

( )

( 2 )

(2)

 

 

 

( Ф

І ( І ) +

С А , О, ^

=

Qo^x

( т І =

Ь І

< о ) ,

С0 г2 (0) = - Ф

і ( 1 ) .

<?U(—оо) = 0

(г =

1,2).

 

 

Приравнивая

(4.360)

и (4.364)

при t — T и

учитывая,

 

 

 

(1)

(2)

 

 

 

 

 

что в этой

 

точке П0 г(

и Q0zi экспоненционально малы,

в нулевом

приближении

получим

 

 

 

 

 

 

 

< Qo z1 (0) = no z1 (0),

 

 

(4.367)

 

 

Фз (Т)

+ Q0 23

(0) = Ф і

(T) + S U

(0).

 

(4.368)

Равенство (4.368) обеспечивается условиями (4.363) и (4.366), а равенство (4.367) представляет собой уравне-

(1)(2)

ние относительно T (Q0z1(xT)

И

 

 

Т) зависят от Т, так

как Т входит в правую часть уравнений (4.362) и (4.365)

и

дополнительные

условия (4.363)

и (4.366)).

(4.363) за-

 

С целью

определения Т сделаем в (4.362),

 

 

 

 

 

(1)

 

 

_

 

(D

мену переменных zx

(xT) = Q0z1

г),

z2

( т г ) = Ф з ( Т ) 4 - Q 0 z 2 г).

В

этих

переменных (4.362),

(4.363)

примут вид

 

 

 

§T

=

F ( Z 2 ,

Т ) ,

^

= 2,,

(4.369)

 

22 (0) = Ф 2 (Л - ^

( - « , ) = О,

2 f

( - o o ) = q>8(7).

Точно

так же в (4.365), (4.366)

сделаем замену г2 Г ) =

 

(2)

 

=

 

(2)

 

 

 

 

 

=

U0zx

г),

г2 г ) = Ф і

(Т) + Пц г2

т).

Получим

 

 

 

WT

= F & '

П

§Т

= 1 ,

 

 

(4.370)

 

г2 (0) = ф 2 ( Л ,

Zi (оо) = 0,

22 (оо) = ф1

( Т ) .

В

новых переменных

уравнение

 

(4.367) примет вид

( 0 ) ^ ( 0 ) . (4.371)

Нетрудно видеть, что как (4.369), так и (4.370), с точ­ ностью до обозначений, представляет собой систему (4.357).

§ 15]

ВНУТРЕННИЙ ПОГРАНИЧНЫЙ слой

197

Предположим, что ей отвечает фазовая картина рис. 6, а).

Тогда значению гг (0) отвечает точка А, а значению гх (0)— точка В, лежащие на вертикальной прямой / рис. 8, и условие (4.371) не выполняется. Точно так же (4.371) не может выполняться, если имеет место фазовая картина рис. 6, б). И лишь фазовая картина рис. 6, в) обеспечит

2 3

7

— ^

в

Рис. 8. 7 —прямая z2 = <p2

на фазовой плоскости ( г ь

г2 ), отвечающей

t = T,

2

и 3—два

различных

положения прямой

г2

= 0 на фазовой

 

 

 

плоскости

І Р

г2 ), отвечающей t =

0.

 

 

 

нужное

равенство. Итак,

мы приходим к выводу, что Т

должно

быть

тем значением,

при

котором

на

фазовой

плоскости

образуется

ячейка.

Это будет сформулировано

ниже

в

виде требования

IV,

после

того

как

уравнение

(4.371) будет преобразовано к надлежащему виду.

Кроме этого, для реализации каждого полуцикла фа­

зовые

картины, при ^ = 0 и t= \ должны

удовлетворять

определенным условиям,соответствующим требованиям V, а)

и V, б) § 14.

В точке

^ = 0 для левого полуцикла

нужно

обеспечить выполнение требования V, а), а в

точке ^ = 1

для правого полуцикла—выполнение

требования

V, б).

Заметим,

что

требование

V, б) для левого полуцикла и

V, а) для

правого полуцикла,

т. е. условия в точке t = T,

автоматически удовлетворяются вследствие самого требо­

вания

существования

ячейки

при t T.

Действительно,

точка

 

(0, ФаСО)

лежит

внутри ячейки и,

следовательно,

прямая

z2 = ф2 (Т)

пересекает

сепаратрису,

идущую из

седла

 

(0, ФзСГ)) в седло (0,

ух{Т)). Иначе

говоря,

зна­

чение

z2 = ф2 (Т)

принадлежит

области определения функ­

ции

z1

= f(z2),

описывающей

многообразие

S"

седла

(0, Фя(Т)), которое

совпадает

с многообразием S +

седла

198

КРАЕВЫЕ

ЗАДАЧИ

[ГЛ. 4

(О. Фі<Т)) и

является именно

той сепаратрисой,

которая

соединяет эти седла. То же самое относится и к г2 = ф г (Т). Выясним теперь, что нужно, чтобы удовлетворить требованию V, а) для левого полуцикла. На фазовой плоскости, отвечающей'присоединенной системе для ^ = 0

<і)

(она

 

получается

из

(4.361)

заменой

Z1(T0)

=

 

UOZ1(XQ),

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z2 0 ) = фз (0) + П 0 2 2

0 ) и

имеет вид (4.357), если опустить

значок

~

у

 

z (

( i = l ,

2), и положить

в

(4.357)

Т = 0,

т = т0 ),

имеется, вообще говоря, петля (рис. 6,

а) или 6,6)).

Например,

для случая

рис. 6,

а)

условие

V, а) означает,

что

прямая

z2 = 0 не должна

находиться

вне петли,

так

как

иначе

эта прямая

(прямая

2 на рис. 8)*) не пересе­

кала

 

бы многообразия

S+,

отвечающего (0, Ф3 (0)) и гео­

метрически

 

изображаемого

этой

петлей. Таким

образом,

для

выполнения

условия

V, а) требуется,

чтобы прямая

z3 = 0

пересекала петлю

(прямая

3 на

рис. 8). Аналогич­

ное

требование

должно

выполняться

и при t=\.

Итак,

I I I .

Пусть

прямая

z2

=

0 на фазовой

плоскости,

 

отве­

чающей

системе

^-

= F(z2,

0),

^~ =

zi>

пересекает

сепа­

ратрису,

входящую

в седло (0, Ф3 (0))

при

т0 «-сю.

 

 

Пусть

прямая

z2 = 0 на

фазовой

 

плоскости,

отвечаю­

щей

системе

d^- = F(z2,

1), ^

 

 

пересекает

сепарат­

рису,

 

входящую

в седло

(0, ф х (1)) при

хгу

— с о .

 

 

Напишем теперь уравнение (4.371) в другой форме.

Равенство

z1 = f(z2)

представляет

собой

уравнение

сепа­

ратрисы системы (4.369) (или (4.357)),

 

исходящей

из

точки

(0,

фР

(71))

при

т г =

с о .

Выражение

для / ( z 2 )

 

 

 

 

 

 

 

 

 

 

 

ф 2 ( Т)

 

 

 

 

 

 

 

 

 

 

получим из (4.358)

при с—

\

F (z,

 

 

T)dz,

 

 

 

 

~zx

 

 

 

 

 

 

=

- V

l

\

 

J F(z,

 

T)dz.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ф , (Л

 

 

 

 

 

 

 

 

 

 

*) Выше мы считали рис. 8 относящимся к присоединенной си­ стеме в точке t = T, а теперь — к присоединенной системе в точке / = 0 (чтобы не делать лишнего рисунка).

ВНУТРЕННИЙ ПОГРАНИЧНЫЙ СЛОЙ

199

Точно так же для сепаратрисы системы (4.370), входящей

в точку

(0,

Ф1 (Т))

при т г = оо, получим

 

 

 

г1 =

-

] /

2 "І // 5 \ (Т) F (г,

T)dz.

 

Уравнение

(4.371)

дает

 

 

 

_

/ " Ф І m

 

 

 

Г ф 2 (Г)

 

— / 2

у

J

F(z,

T)dz = — / 2 | /

J

F(z,T)dz.

Отсюда

 

Ф . ( Л

 

 

 

 

Ф , ( Г )

 

 

 

 

 

 

 

 

 

 

 

Фі(Г)

 

Фг ( Л

 

 

 

 

J

F(z,

=

J ^(z,

0^ г .

(4.372)

 

 

<Р.(Г)

 

Ф

І (T)

 

 

что совпадает с условием ячейки (4.359). Таким образом, для существования в некоторой точке отрезка [0, 1]

ячейки,

обеспечивающей

переход

с корня ф3 (t)

на корень

Фі (0. нужно,

чтобы выполнялось

условие

 

IV.

Пусть

уравнение

(4.372)

относительно

Т имеет

решение

Т = Т° и пусть

0 < Т° <

1.

 

Мы

рассмотрели

условия

формального

построения

асимптотики

решения

такого

типа, как на рис. 7. Вы­

ясним

теперь

достаточные условия его фактического су­

ществования.

Выражение

 

 

 

ß ^ ^ Q U ^ - f i U (0)

обращается в

нуль

при Т~Т".

Пусть

 

В(Т)

меняет знак

при

переходе

через

Т = Т°. Выражение

В (Т) с точностью

О (р,)

представляет

собой разность значений

гг для левого

и правого полуциклов в точке

t T,

в

которой z2 (и для

того и для другого

полуцикла)

равно

у2(Т).

Отсюда сле­

дует, что в окрестности Т° найдется Т(р) такое, в кото­ ром разность значений zx для левого и правого полуциклов обращается в нуль, т. е. левый полуцикл соединяется с правым с сохранением непрерывности как z2 , так и z\, другими словами, действительно существует решение за­ дачи (4.355) такого вида, как на рис. 7.

Величину В (Т), пользуясь преобразованиями, при­ водящими (4.371) к (4.372), можно записать в виде,

Соседние файлы в папке книги из ГПНТБ