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

§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

(»)

9z(x)

 

Z

a0

 

 

<

2

(»)

=

 

 

x0

 

 

 

:

 

 

 

 

;

1

Лекция №15, НГУ, ММФ, 2010

2

приводится к уравнению (10) для функции z(x). Умножая уравнение (1) на функцию

8

1 <Zx

¹(x) = a0(x) exp :

x0

9

a1(»)=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 = 1; 2; ::: . Сле-

 

b

2¼2

 

 

 

довательно, ¸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(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(») = 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):

Соседние файлы в папке Диф.уравнения