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

2 курс / Дифференциальные уравнения / Федорюк М., Обыкновенные дифф.уравнения, 1985

.pdf
Скачиваний:
844
Добавлен:
21.05.2015
Размер:
21.78 Mб
Скачать

§ 11] УРАВНЕНИЯ С ПЕРИОДИЧЕСКИМИ КОЭФФИЦИЕНТАМИ 221

Наметим доказательство (подробнее см. [281). Те же рассуждения, что и выше, приводят к тождеству (5), где С — постоянная невырожденная (п X п)-матрица. Поэтому существует матрица В такая, что е*в = С (§ 2). Положим

Р(х) = YU)e~xB, тогда

Р(х + со) - Пх)Се-»ве-*в = Р(х),

т. е. Р(х) есть периодическая с периодом со матрицафункция.

2. Зоны устойчивости и неустойчивости. Различные

физические и технические задачи приводят к уравнениям второго порядка, содержащими вещественный параметр

к>0, вида

 

-у"

+q(x)y

= k2y

(12)

или

вида

 

 

 

 

у"

+k2q(x)y

= 0.

(13)

Здесь q(x) —вещественная со-периодическая

функция.

К

уравнению вида (12) приводит, например, задача о дви-

жении квантовомеханической частицы в периодической кристаллической решетке, к уравнению вида (13) — задачи

ораспространении электромагнитных или звуковых волн

всредах с периодическими свойствами.

Физическая постановка задачи приводит к условию:

решение у{х) должно бытьограничено на всей оси —°о < <х<°°. Это условие не может выполняться на полуоси вида к > к0 (за исключением очевидного случая q(x) в 555 const, который мы далее не рассматриваем). По определению, точка к2 принадлежит зоне устойчивости, если все решения уравнения (12) (или (13)) ограничены при —оо < х < о©, и принадлежит зоне неустойчивости в противном случае. Из вида решений (2а) следует, что точка к2 принадлежит зоне устойчивости, если

 

 

 

(14)

Действительно, в этом случае, в силу (3),

 

 

|г/Д* + ш)1 = \р,\\у,(х)\

- \у,Ш,

1-U

2,

и потому решения у^х), у2(х)

и все их линейные комби-

нации ограничены на всей оси х. Если

же

Ip^l Ф1, то

точка к2 принадлежит зоне неустойчивости. Пусть IpJ >

>

1, для определенности (напомним, что piPz — D, тогда

из

(3) имеем

222

ЛИНЕЙНЫЕ

УРАВНЕНИЯ И СИСТЕМЫ

[ГЛ. 3

Фиксируем

точку

х0

такую,

что yi(x0:09

тогда

\yi(xQ + n(i))\-> оо при

/г-^«5> П

решение

уг(х)

не явля-

ется ограниченным.

 

 

 

 

 

К сожалению, единственный

пример

интегрируемого

уравнения

вида (12)

или

(13),

где q(x) — элементарная

функция, есть уравнение с кусочно постоянной функцией q(x). Но даже в этом случае анализ зон устойчивости и неустойчивости оказывается весьма сложным. Проще всего этот анализ проводится в примере, принадлежащем

П.Дираку [52]. Рассмотрим уравнение

'• — U Z 6 (* — mo) p = О,

(15

где U Ф 0 — постоянная и б есть дельта-функция Дирака (§ 12). Функция q(x), входящая в уравнение (15), называется диракоеской потенциальной гребенкой [52J. На интервале 0 < х < ш всякое решение уравнения (15) имеет вид

где А, В— постоянные. Будем искать решение у(х) такое, что у(х + ы) = ру(х) (см. (3)), так что

у{х) - р[АеЩх'а) + Ве~тх-а)]

при со< х < 2со. В точке х = со должны выполняться условия (§ 12)

г/(со + 0) = уЫ - 0), г/'(со + 0) - z/'(co - 0) + Uy(<* - 0),

и мы получаем систему уравнений относительно неизвестных А, В:

р(А+В) *= Aeih« + Ве~**,

ikp(A - В) = ik(Aeik« - Ве~гЫ) + tfUe** + Дв"**).

Система имеет нетривиальное решение тогда и только тогда, когда ее определитель равен нулю, и это дает уравнение для мультипликаторов

р2 — 2/ {к)р + 1 = 0, / (к) « cos ok + sin шА:.

При выводе этого уравнения следует помнить, что свободный член обязательно равен коэффициенту при р2# так как р^ 2 = 1. Решая уравнение, находим

l. (16)

§ 11] УРАВНЕНИЯ С ПЕРИОДИЧЕСКИМИ КОЭФФИЦИЕНТАМИ 223

Зоны устойчивости определяются из условия

f(ft)<l,

в этом случае

 

 

 

р,, a

 

 

и

Ipil = |ра1 = i. Если же f(ft)>

1, то pi > 1 при

/(fc)>0,

р 2

> 1 при /(ft) < О, и точка ft2 принадлежит зоне неустой-

чивости. Можно проверить, что

если /(/с) = ±1,

то точка

к2 принадлежит зоне неустойчивости. Итак, зоны устойчивости определяются неравенством

 

cos Arco + ^тг sin ft© < 1.

(17)

I

4ml Hi

I

Если построить график функции /(ft), то границами зон устойчивости будут абсциссы точек пересечения графика

с прямыми

г/= dbl. Поэтому зоны устойчивости и неус-

тойчивости

чередуются.

 

Если £/>0, то некоторый интервал вида

(0, ft0) будет

зоной неустойчивости, так как при малых к имеем /(ft) ~

~ 1 + (7со/2 > 1. Покажем, что существует

бесконечно

много зон устойчивости и неустойчивости. Точки экстремума функций /U) определяются из уравнения

Из сравнения графиков этих функций видно, что уравнение имеет бесконечно много корней, видаА*тг= -JJJ-+ Tn»

где Yn ~* 0 при п -* оо. Подставляя их в уравнение, получаем

Так как tg«r —х при х -> 0, то Yn ^ UH2nn)4 и х\1ы получаем асимптотику корней уравнения (1):

Из этой формулы и (18) находим

( 2 0 >

224

ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ •

1ГЛ. 3

так что 1/(/сп)1>1 при п > 1 . Поэтому существует бесконечно много зон устойчивости и неустойчивости. Так как \f(kn)\ -•1 при и-+•«>, то ширина п-й зоны неустойчивости (эти зоны называют также лакунами) стремится к нулю при п-> °°. Пусть /„«•(а„, bn) — лакуна, содержащая точку кп. Положим ап ~ кп —ап, Ь„= fcn + pn, где а„, (*« — положительные бесконечно малые. Имеем из (19), (20)

/ Ы - (-1)п =

так что осп^ U/i2nn). Такую же асимптотику имеет величина рп, так что ширина Дп лакуны 1п равна

Вычислим асимптотику мультипликаторов, отвечаю- щих /г-й зоне неустойчивости. Пусть р(А) — тот из мультипликаторов, модуль которого больше единицы. Из формулы (16) следует, что максимум 1р(/с)| при к^1п достигается в точке кп и из (20) находим

Так как 1р(/сп)1 мало отличается от единицы

при п > 1 ,

то это означает, что с ростом номера зоны

неустойчи-

вость проявляется все слабее.

 

Аналогично устроены зоны устойчивости и неустойчивости для уравнения (12).

Теорема . Пусть qix) непрерывная вещественная (^-периодическаяфункция, q(x) Ф const. Тогда на полуоси

к>0 имеется бесконечно много зон неустойчивости 1п =

«=[дп, fcj, а дополнительные к

ним интервалы суть

зоны устойчивости. Ширина Д„ п-й зоны неустойчивост

стремится к нулю

при п -> ©о.

 

Если функция

q(x) имеет р>1

непрерывных произ-

водных, то

 

 

Д

Доказательство теоремы и этой формулы см. в [48]. Одним из наиболее полно исследованных уравнений с периодическими коэффициентами является уравнение

§ I2J

ДЕЛЬТА-ФУНКЦИЯ И ЕЕ ПРИМЕНЕНИЯ

225

Матье:

у" + {а + Ьcos 2х)у — 0.

В последние годы были исследованы конечнозонные потенциалы q(x). Это такие периодические функции q(x), что уравнение (12) имеет конечное число протяженных зон неустойчивости (т. е. интервалов), а все остальные зоны неустойчивости вырождаются в точки. Для конечнозонных потенциалов получены явные выражения.

§ 12. Дельта-функция и ее применения

1. Дельта-функция Дирака. Эта функция — она обозначается 60r) — была введена английским физиком П. Дираком. Ее определение таково:

б (х) = 0 , х Ф 0, j б (я)dx.- 1.

(1)

Дельта-функцию можно рассматривать, например, как плотность единичной массы, сосредоточенной в точке £ = 0 (это только одна из известных физических интерпретаций дельта-функции). Действительно, обозначим эту плотность р(х); тогда р(х) = 0 при х Ф 0 (вся масса со-

средоточеиа в точке ж= 0), и J р(о:)йя~1, ибо масса

равна единице, так что р(х)

Разумеется, дельта-функция не есть функция в обыч-

зом смысле слова. Это

обобщенная функция.

Теория

)бобщенных функций была построена советским

матема-

тиком G. Л. Соболевым

и французским математиком

F1. Шварцем.

 

 

Основы теории обобщенных функций читатель может тйти в [15, 18]. Мы ограничимся тем, что пpивeдeм^ >сновные формулы, относящиеся к дельта-функции, и не >удем излагать строгую теорию обобщенных функций. 1итатель, незнакомый с обобщенными функциями, будет сходиться примерно в том же положении, в каком нахоились в начале нашего века инженеры, которые испольовали метод Хевисайда (см. гл. 1, § 11). Впрочем, это не овеем так: математического обоснования метода Хевисай- а тогда не было, а математическое обоснование привеенных ниже результатов — теория обобщенных функ- ,ий — существует.

226 ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ [ГЛ. 3

Введем понятие обобщенной функции. Пусть К — множество всех функций ф(х), которые бесконечно дифференцируемы на всей оси —°о < х < ©о и финитны. Последнее

означает, что каждая функция

ф ш ^ й

тождественно

равна

нулю вне

некоторого

отрезка: ф(я) ^,0, если

х Ф-

Ф (а,

Ь)\ числа

а, Ь-—свои

для

каждой

функции

фЫ.

Множество К есть линейное пространство; его элементы О) называются основными (или пробными) функциями.

Обобщенной функцией f (над пространством К) называется линейный функционал, определенный на простран-

стве К. Именно, каждой функции ф Ы е й

ставится в

соответствие число /(ф), причем

 

/(с^ф! + а2ф2)в cCi/Сф!)+ а2/(ф2)

(2)

для любых основных функций <pi(&)-, ф2 Ы и для любых чисел а4, <х2. От функционала / требуется также непрерывность; мы не будем вводить соответствующее определение, поскольку во всех рассматриваемых примерах это свойство выполняется.

Будем употреблять обозначения, принятые в физической литературе. Именно, значение /(ф) будем записывать в виде интеграла

/(ф)- J / Ф Ф (*)**.

(3)

Класс обобщенных функций содержит все «обычные» функции. Действительно, если fix) — непрерывная на оси

—оо<#<оо функция, то интеграл из (3) существует и обладает свойством линейности. Если fi(x), fz(x) — непрерывные на всей оси функции и

J h (*) Ф (*) dx - | /,(х)Ф (х)dx

(4)

для любой основной функции ф(я), то f%(x) s

/2(^) (гл. 6,

§ 3, основная лемма вариационного исчисления). По аналогии с этим фактом будем считать, что две обобщенные

функции /Дя), f2ix)

равны: fi(x)**fz(x),

если соотноше-

ние (4) выполняется для любой основной функции у(х).

Введем производные от обобщенных функций. Пусть

fix} — непрерывно дифференцируемая на

всей

оси функ-

ция, ф(х) — основная

функция. Интегрируя

по частям,

§ 12]

ДЕЛЬТА-ФУНКЦИЯ И ЕЕ ПРИМЕНЕНИЯ

227

получаем

 

 

 

оо

 

 

 

(5)

поскольку внеинтегральная подстановка / (х) <р(х) J!!» равна нулю — функция ф(#) финитна. Формулу (5) примем в качестве определения обобщенной функции fix). Заметим, что правая часть формулы (5) определена, так как ср'Ы — основная функция, и потому интеграл

) / (я) ф' (?)dx определен. Аналогично определяются выс-

~- оо

шие производные:

 

ОО

 

 

00

 

 

J /(п> (*) Ф(«) dar - (—l ) n

J /(яг) Ф'п ) (ос) их.

(6)

 

—оо

 

—оо

 

Приведем основные формулы

для дельта-функции

и

поясним их.

 

 

 

1°.

J

b(x — a) y{z)dx = ф(а).

 

Здесь

фЫ — непрерывная на

всей оси функция. Эта

формула есть определение дельта-функции.

 

Дельта-функцию можно представить

 

как

«предел» обычных

функций.

 

Пусть

6Ах) — ступенчатая

функция

 

(рис. 18):

бе (4

так что

i

Рис, 18.

Функцию 6«Ы можно интерпретировать как плотность единичной массы, «размазанной» на интервал (—е, е). Покажем, что

Н т б , (*)-«(*),

(7)

15*

228

ЛИНЕЙНЫЕ УРАВНЕНИЯ Й СИСТЕМЫ

{ГЛ. 3

где предел понимается так:

 

 

lim J

6z{x)<({x)dx** f

S(z)(?{z)dx**(¥(O)

(8)

для любой непрерывной

на всей оси функции фЫ. Дей-

ствительно, по теореме о среднем

 

 

 

 

е

 

 

 

 

J

6,(*)q> (*

 

 

 

 

откуда следует (8).

 

 

 

 

Обозначим

символом

0(я) ступенчатую

функцию

Хе-

еисайда:

 

 

 

 

 

 

 

•w-fc

x<o°:

 

»

20

£W \

А / \

 

 

 

 

 

• ' ' ' и \Х) === Оlu/j»

 

 

 

 

Действительно, пусть фЫ — основная

функция. Со-

гласно определению (5),

 

 

 

 

с©

 

оо

 

 

 

 

J 6'(х) ф(х) dx - — j 6 (х) ф'(z) dx

 

 

о

(так как фЫ — финитная функция)

Поскольку первый и последний интегралы равны для любой основной функции фСг), то Q'(x) ~б(а:), согласно определению равенства обобщенных функций (см. (4)).

3°. Пусть функция fix) непрерывно дифференцируема на полуосях х < а и х > а^-ъ в точке а может иметь разрыв. Тогда

/Чх) «{/'(*)>+ Л6(х - а) .

Здесь А =/( а + 0) —/(я — 0) (т. е. величина скачка функции /Ы), а {/Чх)> — обычная производная.

12]

ДЕЛЬТА-ФУНКЦИЯ И ЕЕ ПРИМЕНЕНИЯ

229

Действительно, по определению (5) имеем

 

 

 

а

 

 

- -

J /(s)q/(*)d*

 

 

 

 

—оо

а

(интегрируем по частям)

 

 

 

 

 

= J {f'(z)}4>(z)dz+A<?{a).

 

 

 

 

—'ОО

 

 

 

00

 

Последнее

слагаемое равно J

48 (#—•а) ф{х) dx, и фор-

 

 

 

— ОО

мула 3°

доказана.

 

 

ОО

 

 

 

4°.

j

б( п)

(х - а) ф (*) < t o -v

„•1)п ф( я ) (а).

— ОО

Здесь (fix) — любая основная* функция; если п фиксировано, то достаточно, ч^обы функция (fix) была п раз непрерывно дифференцируема на оси х.

 

5°. 6(ах) = \а\~1Ш)

(а*•()).

 

 

Пусть а > 0 (для

определенности),

<рСя)— непрерыв-

ная

на оси д: функция. Тогда

 

ОО

 

ОО

 

J 6(а*)Ф(*)<**=-! j6(f) ? ( - f)^

=4-<P(°) =

— ОО

—ОО

 

(мы

сделали замену ах *» ^)

 

 

 

 

ОО

 

 

=

J а~1б (я:) ф (д:) dx.

 

 

 

— ОО

Так

как первый и^роследний интегралы равны при любой

непрерывной функции ф(#), то 5° доказано.

ОО

6°.

— ОО

Здесь мы ограничимся совсем уже формальным выводом. Воспользуемся формулой обращения для преобразования' Фурье:

230

ЛИНЕЙНЫЕ УРАВНЕНИЯ И СИСТЕМЫ

[ГЛ.

где /(а) — преобразование Фурье функции / Ы :

со

/( a )- J e-iaxf(x)dx.

—оо

Пусть fix) = 8Ы, тогда f (a) ™ 1; подставляя это выраж ние в формулу обращения, получаем 6°.

7°. Пусть fix) — гладкая функция, имеющая конечнс число нулей xiy ..., х», и все нули — простые (/'Up ^ О Тогда

Рассмотрим вначале случай, когда функция /Ы им ет ровно один нуль х0 и f(xo)>O. Так как 6 ( ( ) ) при х Ф х0, то

где е > 0

может быть выбрано сколь угодно малым. Сд

лаем

в

интеграле замену переменной

fix) f{xQ) «•

тогда

существует обратная функция x^giy)

и

Действительно, числа е4, е2 положительны, f'{x)g'(y) Если /'(яо) < 0, то получим

•f

так что в обоих случаях