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

Волков Функция Грина самосопряженной краевой задачи второго 2012

.pdf
Скачиваний:
5
Добавлен:
12.11.2022
Размер:
251.11 Кб
Скачать

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ

НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ ЯДЕРНЫЙ УНИВЕРСИТЕТ «МИФИ»

В.Е. Волков, Е.В. Сумин

Функция Грина самосопряженной краевой задачи второго порядка

Учебно-методическое пособие

Москва 2012

УДК 517.968(076) ББК 22.161.6я7 Ф 94

Волков В.Е., Сумин Е.В. Функция Грина самосопряженной краевой задачи второго порядка: учебно-методическое пособие.

М.: НИЯУ МИФИ, 2012. – 20 с.

Данная работа представляет собой учебно-методическое пособие к практическим занятиям по интегральным уравнениям.

В работе дано определение функции Грина и приведен метод ее построения для однородной краевой задачи второго порядка в общем случае. Введено понятие самосопряженной краевой задачи и описывается способ построения функции Грина в этом случае. Приведены примеры построения функции Грина, а также задачи для самостоятельного решения.

Предназначено для студентов НИЯУ МИФИ, изучающих курс интегральных уравнений. Будет также полезно преподавателям, ведущим практические занятия.

Рецензенты: канд. физ.-мат. наук, доц. НИЯУ МИФИ В.Б. Шерстюков;

канд. физ.-мат. наук, доц. МПГУ О.Н. Быкова

Рекомендовано к изданию редсоветом НИЯУ МИФИ

ISBN 978-5-7262-1706-2

©Национальный исследовательский ядерный университет «МИФИ», 2012

Редактор М.В. Макарова

Подписано в печать 20.06.2012. Формат 60х84 1/16

Уч.-изд. л. 1,25. Печ. л. 1,25. Тираж 430 экз.

Изд. № 002-1. Заказ № 110

Национальный исследовательский ядерный университет «МИФИ». Типография НИЯУ МИФИ.

115409, Москва, Каширское ш., 31

§ 1. ПОСТРОЕНИЕ ФУНКЦИИ ГРИНА ДЛЯ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА

С ОДНОРОДНЫМИ КРАЕВЫМИ УСЛОВИЯМИ

Пусть задано линейное однородное обыкновенное дифференциальное уравнение второго порядка с переменными коэффициентами

′′

(x) y(x) = 0

(1)

L[ y] p0 (x) y (x) + p1

(x) y (x) + p2

где p0 (x) , p1 (x) , p2 (x) – известные непрерывные на отрезке [a, b] вещественнозначные функции. При этом предполагается, что

p0 (x) 0 , x [a, b] .

Запишем для уравнения (1) однородные краевые условия в общем виде:

 

 

 

(1)

 

 

(1)

 

 

V1

( y) ≡α1 y(a) 1

y (a) 1 y(b) 1

 

y (b) = 0;

(2)

 

 

 

y(a)

(1)

y (a) y(b)

(1)

y (b) = 0.

V ( y) ≡α

 

 

2

 

 

 

 

 

 

 

2

 

2

 

 

2

 

2

 

 

 

Здесь линейные формы V1 ( y) , V2 ( y) являются линейно независимыми, т.е.

от y(a) , y(b) , y(a) , y(b)

ранг матрицы

 

α

α(1)

β

β(1)

 

 

1

1

1

1

 

 

α2

α(1)2

β2

β(1)2

 

равен 2; коэффициенты линейных форм α1 , α1(1) , …, β2 , β(1)2 – ве-

щественные числа.

Определение. Функцией Грина краевой задачи (1), (2) называется функция G(x, ξ) , построенная для любой точки ξ, a <ξ<b , и обладающая следующими четырьмя свойствами:

1)G(x, ξ) непрерывна при a x b ;

2)первая производная по x в точке x = ξ имеет разрыв первого

рода, причем скачок равен

 

1

 

, т.е.

 

 

 

 

 

p (ξ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

G(x, ξ)

 

 

 

G(x, ξ)

 

 

=

1

;

 

 

 

 

 

x

 

x=ξ+0

 

x

 

x=ξ−0

p (ξ)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

3

3) в каждом из полуинтервалов [a, ξ) и (ξ, b] функция G(x, ξ) , как функция от x, является решением уравнения (1):

L[G] = 0 ;

4) G(x, ξ) , как функция x, удовлетворяет граничным условиям

(2):

V1 (G(x, ξ)) = 0;V2 (G(x, ξ)) = 0.

Утверждение 1. Функция Грина однородной краевой задачи (1),

(2) существует и единственна тогда и только тогда, когда краевая задача (1), (2) имеет лишь тривиальное решение ( y(x) 0 ,

x [a, b]) .

Предположим, что задача (1), (2) имеет только тривиальное решение. Построим функцию Грина для (1), (2) в этом случае.

Пусть функции y1 (x) и y2 (x) образуют фундаментальную систему решений уравнения (1). Тогда согласно свойству 3 функцию

Грина G(x, ξ)

можно представить в виде:

 

 

 

 

 

 

 

G(x, ξ) =

a y (x) +a y (x), a x ≤ξ;

 

 

 

 

1 1

 

2 2

 

 

 

 

 

 

 

 

 

 

b1 y1 (x) +b2 y2 (x), ξ≤ x b,

 

 

где a1 , a2 , b1 ,

b2

– функции от ξ.

 

 

 

 

 

 

 

Из свойств 1 и 2 функции Грина получим соотношения:

 

[b y (ξ) +b y

2

(ξ)] [a y (ξ) +a y

2

(ξ)] =0;

 

 

 

1

1

2

1

1

2

 

1

 

(3)

 

 

 

 

 

 

 

[b1 y1 (ξ) +b2 y2

(ξ)] [a1 y1

(ξ) +a2 y2

(ξ)] =

 

.

 

p0 (ξ)

 

 

 

 

 

 

 

 

 

 

 

 

 

Введем обозначения:

c1 (ξ) =b1 (ξ) a1 (ξ); c2 (ξ) =b2 (ξ) a2 (ξ).

Тогда из (3) получим систему линейных алгебраических уравнений относительно c1 и c2 :

4

c y (ξ) +c y

2

(ξ) =0;

 

 

 

1 1

2

 

1

 

 

(4)

 

 

 

 

 

 

c1 y1(ξ) +c2 y2

(ξ) =

 

 

.

 

p0 (ξ)

 

 

 

 

 

 

 

Заметим, что определитель

 

 

системы (4) представляет собой

определитель Вронского W функций y1 (x)

и y2 (x) при x = ξ :

 

=y1 (ξ) y2 (ξ) =W (ξ) . y1(ξ) y2 (ξ)

Как известно, определитель Вронского для фундаментальной системы решений дифференциального уравнения отличен от нуля во всех точках области определения. Следовательно, определитель системы (4) также отличен от нуля. Отсюда по теореме Крамера существует единственное решение c1 (ξ) , c2 (ξ) системы (4).

Для того чтобы найти вид функций a1 , a2 , b1 , b2 , необходимо воспользоваться краевыми условиями (2). Ради сокращения объема дальнейших записей представим линейные формы V1 и V2 в виде:

V1 ( y) = A1 ( y) + B1 ( y);V2 ( y) = A2 ( y) + B2 ( y),

где

 

y (a) ;

Ak ( y) k y(a) k

 

(1)

Bk ( y) k y(b) k

y (b) , k = 1, 2.

(1)

 

Тогда можно записать:

 

 

Vk (G) = a1 Ak ( y1 ) +a2 Ak ( y2 ) +b1Bk ( y1 ) +b2 Bk ( y2 ) = 0 , k = 1, 2.

Учитывая, что ak =bk ck , k = 1, 2, а значения ck уже найдены,

из свойства 4 функции Грина получим систему для неизвестных b1 и b2 :

(b1 c1 ) A1 ( y1 ) +(b2 c2 ) A1 ( y2 ) +b1B1 ( y1 ) +b2 B1 ( y2 ) = 0;(b1 c1 ) A2 ( y1 ) +(b2 c2 ) A2 ( y2 ) +b1B2 ( y1 ) +b2 B2 ( y2 ) =0.

Преобразуем эту систему к виду:

5

bV ( y ) +b V ( y

) =c A ( y ) +c A ( y

);

1

1

1

2

1

2

 

1

1

1

2

1

2

 

bV ( y ) +b V ( y

 

) =c A ( y ) +c A ( y ).

1

2

1

2

2

2

1

2

1

 

2

2

2

По-прежнему система уравнений остается линейной относительно величин b1 и b2 . Определитель этой системы

=V1 ( y1 ) V1 ( y2 ) V2 ( y1 ) V2 ( y2 )

отличен от нуля в силу предположения о линейной независимости форм V1 и V2 . Поэтому существует единственное решение b1(ξ) ,

b2 (ξ) , а по формулам ak (ξ) =bk (ξ) ck (ξ) , k = 1, 2, однозначно определяется и a1(ξ) , a2 (ξ) . Таким образом, функция Грина построе-

на.

Заметим, что указанный метод построения функции Грина можно рассматривать как доказательство теоремы о существовании и единственности функции Грина.

Пример 1. Построить функцию Грина для дифференциального уравнения

y′′+k2 y =0 ( k > 0 ,

k 2πn , n ` )

(5)

с периодическими граничными условиями

 

y(0) = y(1);

(6)

y (0) = y (1).

 

Решение. Прежде всего, необходимо убедиться в том, что краевая задача (5), (6) имеет лишь тривиальное решение. Действительно, так как фундаментальной системой решений для уравнения (5) являются функции:

y1 (x) =sin kx , y2 (x) = cos kx ,

то общее решение (5) есть функция вида

y(x) = Asin kx + Bcos kx ,

где A и B – произвольные постоянные.

6

Краевые условия (6) дают следующие соотношения для определения A и B:

Asin k + B(cos k 1) = 0;

(7)

 

cos k) + Bk sin k = 0.

Ak(1

 

Система (7) обладает единственным решением A = B = 0 , так

как ее определитель

= 2k(1cos k) 0 при k >0, k 2πn ,

n `.

Итак,

задача (5),

(6) имеет только тривиальное

решение

y(x) 0,

а значит, для нее можно построить единственную функ-

цию Грина G(x, ξ) . Проведем это построение. Для этого предста-

вим искомую функцию G(x, ξ)

в виде:

 

 

 

a sin kx +a cos kx

при

0 x ≤ ξ;

(8)

G(x, ξ) = 1

2

при

ξ≤ x 1,

b1 sin kx +b2 cos kx

 

где a1 , a2 , b1 , b2 – неизвестные функции от ξ, подлежащие определению. Положим

c1 (ξ) =b1 (ξ) a1

(ξ);

(9)

c2 (ξ) =b2 (ξ) a2 (ξ)

 

и выпишем систему линейных уравнений для нахождения c2 (ξ) (см. систему (4)):

c1 sin kξ+c2 cos kξ =0;c1k cos kξ−c2k sin kξ =1.

Решая эту систему, получим

c

= cos kξ

,

с

= −sin kξ .

1

k

 

2

k

 

 

 

 

 

c1 (ξ) и

(10)

Далее воспользуемся свойством 4 функции Грина, согласно которому она должна удовлетворять краевым условиям задачи, т.е.

G(0, ξ) =G(1, ξ);Gx(0, ξ) =Gx(1, ξ).

7

В нашем случае эти соотношения принимают вид:

a2 =b1 sin k +b2 cos k;

ka1 = kb1 cos k kb2 sin k.

С помощью (9), (10) из последних равенств найдем все неизвестные коэффициенты в (8):

a

=

 

sin kξsin k

 

cos kξ ,

1

 

 

2k(1cos k)

2k

 

 

 

a2 =

 

cos kξsin k

+ sin kξ ,

 

2k(1cos k)

 

 

 

2k

b =

 

sin kξsin k

 

+ cos kξ ,

1

 

2k(1cos k)

2k

 

 

b

=

 

cos kξsin k

sin kξ .

2

 

 

2k(1cos k)

2k

 

 

 

Подставляя найденные значения a1 ,

a2 ,

b1 , b2 в (8) и производя

простейшие тригонометрические преобразования, найдем

 

 

 

 

 

1

 

 

cos k x −ξ+

2

 

 

 

 

 

 

 

 

, 0 x ≤ξ;

 

 

k

 

 

 

 

2k sin

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

G(x, ξ) =

 

 

 

 

1

 

 

 

ξ− x +

 

cos k

 

 

 

 

 

 

 

 

2

 

,

ξ≤ x 1.

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

2k sin 2

 

 

 

 

 

 

 

 

 

8

§ 2. ПОНЯТИЕ САМОСОПРЯЖЕННОЙ КРАЕВОЙ ЗАДАЧИ

В примере 1 получилась симметричная функция Грина, т.е.

G(x, ξ) =G(ξ, x)

(при замене x на ξ, а ξ на x выражение для функции Грина не меняется). Если бы мы заранее знали, что функция Грина будет симметричной, мы могли бы несколько сократить объем вычислений. В частности, не понадобилось бы находить величины b1 и b2 , а также

проводить преобразования для G(x, ξ) на промежутке ξ≤ x 1.

Оказывается, еще до решения задачи можно узнать, будет ли построенная функция Грина симметрична или нет. Сформулируем соответствующее утверждение.

Утверждение 2. Функция Грина краевой задачи (1), (2) симметрична тогда и только тогда, когда краевая задача (1), (2) является самосопряженной, т.е. являются самосопряженными уравнение (1) и отвечающие ему краевые условия (2).

Уравнение (1) является самосопряженным тогда и только тогда, когда p0(x) p1 (x) для x [a, b] , т.е. оператор L[ y] можно записать в виде L[ y] = ( p0 y(x))′+ p2 y(x) .

Краевые условия (2), отвечающие уравнению (1), являются самосопряженными тогда и только тогда, когда

p (b)

α

α(1)

= p (a)

β

β(1)

.

α

1

1

β

1

1

0

2

α(1)

0

2

β(1)

 

 

 

2

 

 

2

 

Возвращаясь к примеру 1, убедимся в том, что поставленная задача является самосопряженной. В самом деле, так как в уравнении

(5) p0 (x) =1 ,

p1 (x) = 0 , то получим

p0(x) p1 (x) , x [0, 1] , т.е.

уравнение (5)

является самосопряженным.

Для краевых условий (6) имеем:

 

 

y(0) y(1) = 0;

 

= 0,

 

y (0)

y (1)

9

откуда

с учетом

(1),

(2),

и

(5)

(здесь

 

p0 (0) = p0 (1) =1, α1 =1,

β = −1,

α(1)

(1)

=0 ,

α(1)

=1 ,

β(1)

=−1, α

2

2

= 0 ) получаем, что

1

1

1

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

1

 

1

 

0

 

 

=1

 

1 0

 

 

=1 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

1

 

 

 

 

0 1

 

 

 

 

 

 

Таким образом, периодические граничные условия (6) для самосопряженного уравнения (5) являются самосопряженными, а следовательно, задача (5), (6) также самосопряженная.

Рассмотрим важный частный случай – распадающиеся краевые условия (2).

Пусть в уравнении (1) p0 (x) 0 , x [a, b] . Вообще говоря, последнее требование можно ослабить:

p0 (a) = 0 ; p0 (x) 0 , x (a, b] .

Предположим, что краевые условия (2) являются распадающимися: одно задается на левом конце (в точке a), другое – на правом

(в точке b) (в случае

 

p0 (a) = 0

 

обычное краевое условие в точке

 

x = a заменяется

на естественное условие

ограниченности:

 

y(a)

 

<∞ ):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

(11)

 

 

 

α1 y(a) 1

 

y (a) = 0;

 

 

 

 

 

 

y(b)

(1)

 

(12)

 

 

 

β

 

 

 

 

 

2

2

y (b) = 0.

 

 

 

 

 

 

 

 

 

Заметим, что краевые условия вида (11), (12) для любого урав-

нения (1) и для любых значений α ,

α(1)

, β

2

, β(1)

являются самосо-

пряженными, так как

 

 

 

1

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

α

α(1)

 

=

 

0

0

 

=0

,

 

 

 

 

 

 

 

1

1

 

 

β2

β2(1)

 

 

 

0

0

 

 

 

 

 

 

 

 

а условие ограниченности y(a) <∞ относится к самосопряжен-

ным.

Построим функцию Грина для краевой задачи (1), (11), (12).

10

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]