2 курс / Дифференциальные уравнения / Федорюк М., Обыкновенные дифф.уравнения, 1985
.pdf§ 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 + j£ 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
так что в обоих случаях