All lectures pdfs / Диф.уравнения / Лекция_15
.pdf§15. Линейное уравнение 2го порядка.
В этом параграфе мы рассмотрим линейное уравнение 2го порядка
a0(x)y00 + a1(x)y0 + a2(x)y = f(x): |
(1) |
Будем рассматривать уравнение (1) на отрезке [a; b], причем будем полагать, что a0(x) =6 0 на [a; b]. Как известно, Задача Коши для уравнения (1) ставится так:
a0(x)y00(x) + ::: = f(x); |
x |
[a; b]; |
(5) |
|
(З.К.) (y(x0) = y0; y0(x0) = y1; |
x20 |
2 |
(a; b): |
|
|
|
|
|
Задача Коши разрешима и решение единственно, если коэффициенты a0(x); a1(x); a2(x) и правая часть f(x) непрерывны на отрезке [a; b]. Задача Коши с нулевыми начальными данными имеет только тривиальное решение для однородного уравнения (т.е. при f(x) ´ 0 на [a; b]).
Заметим, теперь, что вместо уравнения (1) мы можем не нарушая общности рассматривать либо уравнение (мы полагаем, что
f(x) ´ 0 в (1)) |
y00(x) + Q(x)y(x) = 0; |
(10) |
||
|
||||
либо уравнение в так называемом самосопряженном виде |
|
|||
f |
P (x)y0(x) |
g |
0 + q(x)y(x) = 0: |
(100) |
|
|
|
Действительно, при a0(x) =6 0 и f(x) ´ 0 уравнение (1) после замены
y(x) = exp |
8 |
1 |
x |
a1 |
(») |
d»9z(x) |
|
Z |
a0 |
|
|||
|
< |
2 |
(») |
= |
||
|
|
x0 |
|
|
||
|
: |
|
|
|
|
; |
1
Лекция №15, НГУ, ММФ, 2010 |
2 |
приводится к уравнению (10) для функции z(x). Умножая уравнение (1) на функцию
8
1 <Zx
¹(x) = a0(x) exp :
x0
9
a1(»)d»=z(x) a0(») ;
(опять при условии, что a0(x) =6 0; f(x) ´ 0), приведем его к виду
(100).
Пример. Простейшими уравнениями 2го порядка являются следующие уравнения:
y00 + a2y = 0;
y00 ¡ a2y = 0; a = const:
Общее решение этих уравнений имеет вид:
y(x) = C1 sin(ax) + C2 cos(ax); y(x) = C1eax + C2e¡ax:
Из этих формул следует, что у уравнений 2го порядка есть решения, имеющие бесконечно много нулей, а есть решения не обращающиеся в ноль.
Определение. Решение y(x) уравнения (1) (при f(x) ´ 0) назовем колеблющимся на отрезке [a; b], если оно более 1го раза обращается в нуль на этом отрезке. В противном случае - неколеблющимся.
Без доказательства приведем несколько теорем, связанных с нулями решений уравнения (10):
y00(x) + Q(x)y(x) = 0: |
(10) |
Теорема 1. Колеблющееся решение не может иметь точки накопления нулей. Иными словами, все нули решения уравнения (10) суть изолированные точки.
Теорема 2. Пусть для уравнения (10): Q(x) · 0; x 2 [a; b]. Тогда, все решения его не колеблющиеся на [a,b].
Лекция №15, НГУ, ММФ, 2010 |
3 |
Теорема 3. Между двумя последовательными нулями одного решения уравнения 2го порядка лежит точно один нуль другого решения того же уравнения, линейно-независимого с первым. Иными словами, нули линейно-независимых решений разделяют друг друга.
Теорема 4 (Теорема сравнения). Пусть заданы два уравнения вида (10):
y00(x) + Q(x)y(x) = 0;
z00(x) + q(x)z(x) = 0; x 2 [a; b]
и Q(x) ¸ q(x) на [a; b]. Тогда между двумя последовательными нулями решения второго уравнения найдется, по крайней мере, один нуль решения 1го уравнения. Иными словами, решение 1го уравнения колеблется быстрее.
Рассмотрим краевую задачу следующего вида:
8
>y00(x) + Q(x; ¸)y(x) = f(x); x 2 (a; b);
<
где - |
> |
y(a)l1 + y0(a)l2 = 0; |
(2) |
||
|
: |
y(b)r1 |
+ y0(b)r2 |
= 0; |
|
|
|
|
¸ некоторый, вообще говоря комплексный параметр. Поскольку, вместо уравнения 2го порядка мы можем рассмотреть систему
Ãy0(x)! |
0 |
= A |
Ãy0(x)! |
+ |
Ãf(x)! |
; A = µ |
¡Q |
0¶ |
; |
y(x) |
|
y(x) |
|
0 |
|
0 |
1 |
|
то краевая задача (2) есть частный случай краевой задачи из §13
(N = 2; M = 1), при этом: L = (l1; l2); R = (r1; r2). Пусть в (2): f(x) ´ 0; Q(x; ¸) = ¸¡q(x), q(x) - непрерывная функция на [a; b],
величины l1; l2; r1; r2 не зависят от ¸ и l12 + l22 =6 0, r12 + r22 6= 0. Пусть, далее
|
|
|
|
l1 |
|
|
l2 |
|
|||
|
|
p |
|
= cos ®; |
p |
|
|
= sin ®; |
|||
|
|
r1 |
r2 |
||||||||
|
|
|
|
l2 |
+ l2 |
l2 |
+ l2 |
||||
1 |
2 |
|
1 |
2 |
|
|
|||||
|
|
|
= cos ¯; |
|
= sin ¯: |
||||||
p |
|
p |
|
||||||||
|
r12 + r22 |
r12 + r22 |
Лекция №15, НГУ, ММФ, 2010 |
4 |
Тогда задача (2) перепишется в виде так называемой краевой задачи Штурма-Лиувилля (З.Ш.Л.):
|
> |
y00(x) + [¸ ¡ q(x)]y(x) = 0; x 2 (a; b); |
|
|
<y(b) cos ¯ + y0(b) sin ¯ = 0: |
|
|
|
8y(a) cos ® + y0(a) sin ® = 0; |
(3) |
|
Число |
¸ = ¸> |
|
|
:0 называется собственным значением З.Ш.Л., |
если она при ¸ = ¸0 имеет нетривиальное решение y = y(x; ¸0), которое называется собственной функцией.
Пример. Пусть в (3): q(x) ´ 0; ® = ¯ = 0; a = 0. Пусть ¸ · 0. Тогда все решения уравнения y00 + ¸y = 0 неколеблющиеся и следовательно граничным условиям y(0) = y(b) = 0 удовлетворить не могут. Иными словами, при ¸ · 0 нет собственных значений. Пусть ¸ > 0. Тогда, общее решение уравнения будет (положим
¸ = a2)
y(x) = C1 sin(ax) + C2 cos(ax):
Отсюда:
y(0) = C2 = 0; y(b) = C1 sin(ab) = 0:
Поскольку C1 6= 0, то sin(ab) = 0, т.е. ak = |
k¼ |
; k = 1; 2; ::: . Сле- |
||||
|
||||||
b |
||||||
2¼2 |
|
|
k¼ |
|
||
довательно, ¸k = |
k |
- собственные значения, yk(x) = sin µ |
|
x¶ |
||
b2 |
6 |
- собственные функции.
Обсудим, теперь, некоторые свойства собственных функций и собственных значений.
1) Все собственные функции, соответствующие собственному значению ¸0, пропорциональны, т.е. представимы в виде: C1y(x; ¸0). В самом деле, если y(x; ¸0) - собственная функция, то любая функция вида C1y(x; ¸0) - тоже будет удовлетворять уравнению и краевым условиям. Предположим, что нашлась собственная функция !(x; ¸0), линейно-независимая с y(x; ¸0). Тогда, любое решение
уравнения
y00(x) + [¸0 ¡ q(x)]y(x) = 0
Лекция №15, НГУ, ММФ, 2010 |
5 |
представимо в виде:
C1y(x; ¸0) + C2!(x; ¸0);
т.е. любое решение этого уравнения удовлетворяет дополнительным граничным условиям, что, естественно, невозможно.
2) Собственные функции y(x; ¸1); y(x; ¸2), соответствующие разным собственным значениям ¸1; ¸2: ¸1 6= ¸2, удовлетворяют усло-
вию:
Zb
y(x¸1)y(x; ¸2)dx = 0:
a
Это свойство собственных функций называется свойством ортогональности собственных функций. Докажем этот факт.
Поскольку: |
|
y00 |
(x; ¸1) + [¸1 ¡ q(x)]y(x; ¸1) = 0; |
|
> |
||
|
8y(a; ¸1) cos ® + y0(a; ¸1) sin ® = 0; |
||
|
: |
|
¡ |
|
<y(b; ¸1) cos ¯ + y0(b; ¸1) sin ¯ = 0; |
||
|
>y00 |
(x; ¸2) + [¸2 q(x)]y(x; ¸2) = 0; |
|
|
> |
|
|
|
8y(a; ¸2) cos ® + y0(a; ¸2) sin ® = 0; |
||
|
: |
|
|
|
<y(b; ¸2) cos ¯ + y0(b; ¸2) sin ¯ = 0; |
||
то: |
> |
|
|
y(x; ¸2)y00(x; ¸1) ¡ y(x; ¸1)y00(x; ¸2) =
=[y(x; ¸2)y0(x; ¸1) ¡ y(x; ¸1)y0(x; ¸2)]0 =
=¡(¸1 ¡ ¸2)y(x; ¸1)y(x; ¸2);
т.е.
Zb
(¸1 ¡ ¸2) y(x; ¸1)y(x; ¸2)dx = 0:
a
Так как ¸1 =6 ¸2, то
Zb
y(x; ¸1)y(x; ¸2)dx = 0;
a
Лекция №15, НГУ, ММФ, 2010 |
6 |
что и требовалось доказать.
3) Все собственные значения З.Ш.Л вещественные. Допустим противное: ¸0 - комплексное собственное значение, а y(x; ¸0) - собственная функция:
|
|
|
|
|
> |
y00(x; ¸0) + [¸0 ¡ q(x)]y(x; ¸0) = 0; |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
8y(a; ¸0) cos ® + y0(a; ¸0) sin ® = 0; |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
> |
|
|
|
|
|
|
00 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
<y(b; ¸0) cos ¯ + y0(b; ¸0) sin ¯ = 0; |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
8 ³ |
|
|
|
|
´ + [¸0 ¡ q(x)] |
|
|
|
= 0; |
|
|
||||||||||||||
|
|
|
|
|
y(x; ¸0) |
y(x; ¸0) |
|
|||||||||||||||||||||||
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
³ |
|
´ |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
< |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
> |
|
|
|
cos ® + |
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
y(a; ¸0) |
³ |
y0(a; ¸0) sin ® = 0; |
|
||||||||||||||||||||||
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
´ |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
> |
y(b; ¸0) |
cos ¯ + |
|
y0(b; ¸0) |
|
|
|
sin ¯ = 0; |
|
|
|
|
|||||||||||||
т.е. |
|
|
|
- тоже |
|
собственное> |
значение, а |
|
|
|
|
|
|
|
) - собствен- |
|||||||||||||||
¸ |
0 |
y(x; ¸ |
) = y(x; ¸ |
|||||||||||||||||||||||||||
|
|
: |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
0 |
|
|||||||||
ная функция. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Тогда: |
|
|
|
|
|
|
y00(x; ¸0) ¡ y(x; ¸0) ³ |
|
|
|
|
´00 |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
y(x; ¸0) |
y(x; ¸0) |
= |
|
|
= ¡(¸0 ¡ ¸0)y(x; ¸0)y(x; ¸0):
Так как ¸0 =6 ¸0, то
Zb Zb
y(x; ¸0)y(x; ¸0)dx = jy(x; ¸0)j2dx = 0;
a a
т.е.
y(x; ¸0) ´ 0:
В §13 было показано, что все собственные значения З.Ш.Л. суть нули целой аналитической функции
|
det |
0 |
¡ ¢ ¡ ¢ ¡1 |
|
|
|
|
LY (a; ¸) |
|
¢(¸) = |
|
@ RY (b; ¸)A |
: |
|
|
|
|||
|
|
det Y (a; ¸) |
Лекция №15, НГУ, ММФ, 2010 |
|
|
|
|
|
|
|
7 |
|
|
||
В нашем случае она принимает следующий вид: |
|
|
¯ |
|
||||||||
¯ |
1( ¡ )¡ ¡ ¡ ¡10 |
¡ ¡ ¡ jj |
¡2 |
¡ ¡ ¡ ¡ ¡20 |
¡ ¡ ¡ |
|
||||||
¯ |
y¡(b; ¸) cos ¯ + y (b; ¸) sin ¯ |
|
y (b; ¸) cos ¯ + y (b; ¸) sin ¯ |
¯ |
|
|||||||
|
y a; ¸ |
cos ® + y (a; ¸) sin ® |
|
y (a; ¸) cos ® + y (a; ¸) sin ® |
|
|
||||||
¯ |
|
|
¯ |
y1(a; ¸) y2(a; ¸) |
|
|
|
¯ |
|
|||
¯ |
1 |
10 |
y10 (a; ¸) y20 (a; ¸)¯ |
20 |
|
¯ |
|
|||||
¯ |
|
|
j |
2 |
|
|
|
¯ |
|
|||
¢(¸) = ¯ |
|
|
¯ |
|
|
|
¯ |
|
|
¯; |
||
|
|
|
¯ |
|
|
|
|
¯ |
|
|
|
|
|
|
|
¯ |
|
|
|
|
¯ |
|
|
|
|
где y1(x; ¸); y2(x; ¸) - два линейно-независимых решения уравне-
ния
y00 + [¸ ¡ q(x)]y = 0:
Там же было указано, что если собственные значения (или что то же самое, что нули функции ¢(¸)) существуют, то их не более чем счетное число и они изолированы.
Для задачи (3) существование собственных значений можно доказать и своим путем. Приведем, также, без доказательства уточненную теорему сравнения (она используется при доказательстве следующей теоремы 5).
Теорема 40 (Уточненная теорема сравнения.) Пусть заданы два уравнения вида (10) на [a; b]:
y00 + Q(x)y = 0; z00 + q(x)z = 0
с дополнительными начальными условиями Коши:
y(a) = z(a) = sin ®; y0(a) = z0(a) = ¡ cos ®
и пусть Q(x) > q(x) на [a; b]. Тогда решение Задачи Коши y(x) имеет на [a; b] не меньше нулей, чем z(x), причем если нули z(x); y(x) перенумеровать при продвижении от a к b, то kый нуль y(x) расположен левее kго нуля z(x). Иными словами, если мы увеличим коэффициент Q(x) в (10), то все нули решения сдвинутся влево.
Теорема 5 (Теорема осцилляции). У Задачи ШтурмаЛиувилля существует бесконечно много (счетное число) собствен-
ных значений ¸0 < ¸1 < ¸2 < ::: < ¸n < ::: ( lim ¸n = 1), причем
n!1
Лекция №15, НГУ, ММФ, 2010 |
8 |
собственные функции yn(x) = y(x; ¸n) имеют ровно n нулей на
[a; b].
Вместо З.Ш.Л. (3) будем далее рассматривать З.Ш.Л. более
общего вида: |
|
|
|
> |
[p(x)y0]0 + [¸ ¡ q(x)]y = 0; x 2 (a; b); |
|
|
8y(a; ¸) cos ® + y0(a; ¸) sin ® = 0; |
(30) |
||
<y(b; ¸) cos ¯ + y0(b; ¸) sin ¯ = 0: |
|
||
полагать,> |
что З.Ш.Л. не имеет собственного значе- |
||
Будем далее : |
|
|
|
ния ¸ = 0.
Определение. Функцией Грина З.Ш.Л. (30) будем называть функцию g(x; s), удовлетворяющую условиям:
1)g(x; s) непрерывна при a · x; s · b;
2)при x =6 s:
[p(x)gx0 (x; s)]x ¡ q(x)g(x; s) = 0;
3) Функция Грина g(x; s) удовлетворяет краевым условиям:
g(a; s) cos ® + gx0 (a; s) sin ® = 0; g(b; s) cos ¯ + gx0 (b; s) sin ¯ = 0;
4) при x = s gx0 (x; s) испытывает разрыв 1го рода: gx0 (s + 0; s) ¡ gx0 (s ¡ 0; s) = p(1s):
Докажем, теперь, существование и единственность функции Грина. Рассмотрим такую задачу Коши:
(З.К.) |
(p(x)y0)0 ¡ q(x)y = 0; x 2 (a; b]; |
||
(y(a) = sin ®; y0(a) = |
¡ |
cos ® |
|
|
|
|
и найдем решение ее y1(x). Тогда очевидно (!?), что все решения уравнения (p(x)y0)0 ¡ q(x)y = 0, удовлетворяющие первому (левому) краевому условию y(a) cos ® + y0(a) sin ® = 0 выражаются в
Лекция №15, НГУ, ММФ, 2010 |
9 |
виде: C1y1(x). Аналогично через y2(x) обозначим решение задачи
Коши: |
(p(x)y0)0 ¡ q(x)y = 0; x 2 [a; b); |
|||
(З.К.) |
||||
(y(b) = sin ¯; y0(b) = |
¡ |
cos ¯; |
||
|
|
|
при этом все решения уравнения (p(x)y0)0¡q(x)y = 0, удовлетворяющие второму (правому) краевому условию y(b) cos ¯+y0(b) sin ¯ = 0 выражаются в виде: C2y2(x). Кроме того, эти решения (т.е. y1(x) и y2(x)) линейно-независимы, ибо в противном случае они являлись бы собственными функциями З.Ш.Л. (30), отвечающими собственному значению ¸ = 0„ а это противоречит нашему предполо-
жению. Определим функцию Грина g(x; s) так:
(
g(x; s) = |
® y1(x)y2(s); x < s; |
|||||||
b |
|
s |
y |
|
|
x |
; x > s; |
|
|
|
|
|
|||||
|
¯ y |
1( |
2 |
( |
||||
|
b |
) |
|
) |
|
b
где ®;b ¯ - пока произвольные постоянные. При этом ясно, что усло-
b
вия 2), 3) выполнены автоматически. Если ®b = ¯, то выполнено и условие 1). Потребуем, также, выполнения условия 4):
®b fy20 (s)y1(s) ¡ y2(s)y10 (s)g = ®Wb (s) = p(1s);
т.е.
1
®b = W (s)p(s):
Здесь W (x) - определитель Вронского, причем по формуле Лиувилля:
W (x) = W (x0)e¡xR0 |
p0(») d» = W (x0)pp((xx0)): |
|
|||
|
x |
p (») |
|
||
Итак: |
|
|
|
|
|
1 |
y1(x)y2(s); x < s; |
|
|||
g(x; s) = |
|
(y1(s)y2(x); x > s; |
(4) |
||
W (s)p(s) |
чем и заканчивается доказательство существования. Заметим, что окончательный вид функции Грина не зависит от выбора решений
Лекция №15, НГУ, ММФ, 2010 |
10 |
y1(x), y2(x). Из этого замечания следует единственность функции
Грина.
Рассмотрим,8теперь, краевую задачу вида:
>(p(x)y0)0 ¡ q(x)y = f(x); x 2 (a; b); y(a) cos ® + y0(a) sin ® = 0;
|
<y(b) cos ¯ + y0(b) sin ¯ = 0: |
|
|
|
|
|||||||
Перепишем (5)>в виде краевой задачи для системы (см. §13): |
||||||||||||
|
: |
|
|
|
|
|
|
|
|
|
|
|
|
|
Z0(x) = A(x)Z(x) + F (x); x |
2 |
(a; b); |
|
|||||||
|
8LZ(a) = 0; |
|
|
|
|
|
|
|
|
|||
|
> |
|
|
|
|
|
|
|
|
|
|
|
Здесь: |
<RZ(b) = 0: |
|
|
|
|
|
|
|
|
|
||
> |
|
|
|
|
|
|
|
|
|
|
|
|
|
: |
|
|
|
|
|
|
|
|
0 |
1; |
|
|
|
|
|
|
0 |
1 |
|
|
||||
Z(x) = |
y(x) |
; A(x) = |
|
q |
|
p |
|
; F (x) = 0 f |
||||
|
Ãy0(x)! |
|
à |
|
|
|
! |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
p |
¡p0 |
|
@ p A |
|||||||
|
|
|
|
|
|
|
L = (cos ®; sin ®); R = (cos ¯; sin ¯):
Матрица Грина для краевой задачи (50) определяется так:
1)G 0x = AG для всех x =6 s;
2)G(s + 0; s) ¡ G(s ¡ 0; s) = I2;
3)LG(a; s) = 0; RG(b; s) = 0.
Выберем матрицу Грина так: |
|
|
|
||||
|
|
|
8 |
G(x; s) = (Gij(x; s)) = |
|
||
|
1 |
|
y1(x)y20 (s) y1(x)y2(s) |
; x < s; |
|||
= |
|
¡y10 (x)y20 |
(s) y10 |
(x)y2(s) |
|||
W (s) |
> |
µ ¡ |
|
|
¶ |
||
|
> |
y2(x)y10 (s) y2(x)y1(s) |
|
||||
|
|
|
> |
|
|
|
|
|
|
|
< |
µ ¡y20 (x)y10 |
(s) y20 |
(x)y1(s)¶; x > s; |
|
|
|
|
> |
¡ |
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
> |
|
|
|
|
|
|
|
: |
|
|
|
|
(5)
(50)
(6)
где y1(x); y2(x) - вышеупомянутые линейно-независимые решения. Легко убедиться, что свойства 1), 2), 3) выполнены. Сравнивая (6) и (4), видим, что
G12(x; s) = p(s)g(x; s):