Назвем оператор
спряженим з оператором Lu.
Якщо оператор L співпадає з спряженим йому оператором M, то такий оператор називають самоспряженим.
Розглянемо різницю
.При отриманні цього виразу ми додали суму
,але вона дорівнює нулю, так що значення виразу не змінилося.
Одже, вираз vLu – uMv являє собою суму частинних похідних по xi від деяких виразів Pi, тобто
,де
.Розглянемо тепер деякий n-мірний об’єм W, який обмежений кусочно-гладкою поверхнею S.
Користуючись формулою Остроградського-Гауса (3.2), будемо мати
, (5.1)де cos(nx1), cos(nx2),… - направляючі косінуси внутрешньої нормалі до S.
Формула (5.1) носить назву формули Гріна.
Розглянемо рівняння (1.1). Оператори Lu, Mv, а також функції P1 та P2 будуть мати вигляд:
При цьому формула Гріна дає (нормаль внутрішня)
(5.2)Будувати розв’язок будемо методом Рімана, який полягає на використовуванні формули Гріна та дає рішення задачі (1.1) через граничні умови (1.2).
Нехай нам потрібно знайти значення функції u у деякій точці М області (x > x0, t > t0 ) з координатами (x1, t1).
Проведемо через точку М (рис. 2) з координатами (x1, t1) дві прямі, які паралельні координатним осям. Нехай точка P(x0, t1) – це точка перети-ну прямих x = x0 та t = t1, а точка Q(x1, t0) – точка перетину прямих
x = x1 та t = t0. Прямі х = х0, х = х1, t = t0, t = t1 як було показано раніше, є характеристиками рівняння (1.1). Область W буде являти собою прямокутник MPRQ. У цій області ми можемо застосувати метод Рімана для знаходження розв’язку.
Якщо враховувати, що обіг області W відбувається проти годинни-кової стрілки, так що обігаєма площа завжди залишається зліва, формулу (5.2) можна записати у вигляді
(5.2’)З рисунку 2 бачимо, що при цьому
dx = cos(nt)dS,
dt = - cos(nx)dS.
За умови u(x0, t) = j(t) отримуємо:
= 0; = j’(t).За умови u(x, t0) = y(x), отримуємо:
= 0; = y’(x).Рис. 2
Якщо застосувати формулу (5.2’) до прямокутника MPRQ, враховуючи, що на характеристиках QM та PR змінюється лише t, а на характерис-тиках MP та RQ змінюється лише x ,будемо мати:
(6.1)Перетворимо кожен з інтегралів, який стоїть у правій частині (6.1):
(6.2.1) (6.2.2) (6.2.3) (6.2.4)Нехай тепер v(x, t, x1, t1) – деяка функція, яка задовільнює умовам:
Mv = 0, (6.4)
, .При цьому
v(x1, t1, x1, t1) = 1,
(6.5)Розв’язок v(x, t, x1, t1) однорідного спряженого рівняння (6.4), який задовільнює умовам (6.5), називається функцією Рімана. Ця функція не залежить від початкових даних (1.2), та для неї точка (x, t) грає роль аргументу, а точка (x1, t1) – роль параметру. Існування та єдиність такої функції v було доказано методом послідовних наближень.
Оскільки на прямій MP t = t1, а на прямій QM x = x1, то останні члени у формулах (6.2.1) та (6.2.2) обертаються в нуль, і ми отримаємо:
.Формулу (6.1) тепер можна записати у вигляді:
Приводячи подібні, та враховуючи, що v(x1, t1, x1, t1) = 1, u(x0,t) = j(t), u(x, t0) = y(x) та ;
= y’(x), маємо:Звідки знаходимо розв’язок нашої задачі
(6.6)Як ми бачимо, формула (6.6) дозволяє у явному вигляді написати розв’язок данної задачі, оскільки точку М(x1, t1) ми вибрали довільно.
Приклад 1.
Знайдемо функцію Рімана для рівняння
. (7.1)Зробивши заміну змінних
рівняння (7.1) приводиться до канонічного вигляду
при цьому будемо мати a = 0, b = -.
Звернемося тепер до відшукання фунції Рімана v(x, h, x1, h1). Згідно загальної теорії, вона повинна задовольняти спряженому рівнянню
(7.2)та умовам на характеристиках, які проходять через точку (x1, h1):
(7.3)неважко вконатися, що функція
задовільнює як рівнянню (7.2), так і умовам (7.3), слід, це і є шукана функція Рімана.
Приклад 2.
Знайдемо функцію Рімана для рівняння
(x > 0) (7.4)приведемо рівняння (7.4) до канонічного вигляду, для чого складемо рівняння характерстик
xdt2 – dx2 = 0
це рівняння має два різних інтеграла
+ = C1, - = C1,
слід, треба ввести нові змінні x та h за формулами
x = + , h = - (x >0)
приєднаємо до цих рівностей ще одну залежність
тоді рівняння (7.4) перетвориться до канонічного вигляду:
при цьому будемо мати a = 0, b = 0.
Для відшукання функії Рімана нам потрібно знайти частинний розв’язок спряженого рівняння
(7.5)який задовольняв би слідуючим умовам на характеристиках, проведених через точку (x1, h1)
(7.6)Будемо шукати розв’язок рівняння (7.1) у вигляді v = G(s), де
s =.
Тоді для G(s) ми отримаємо слідуюче рівняння:
s(1-s)G’’(s) + (1-2s)G’(s) - G(s) = 0
Це рівняння частинним випадком гіпер геометрічного рівняння Гаусса
s(1-s)y’’ + [g - (1 + a + b)s]y’ - aby = 0
при a = b = , g = 1.
Рівняння Гаусса припускає частинний розв’язок у вигляді гіпергеометрічного ряду
який збігається абсолютно при |s| < 1.
Звідки ясно, що взявши
v = G(s) = F
= 1 +ми задовільним рівнянню (7.5) та усмовам (7.6). Слід, функція
і є функцією Рімана.
Приклад 3.
Знайдемо функцію Рімана для телеграфного рівняння