CC BY
81
УДК 517.956.3
ЗАДАЧА ГУРСА ДЛЯ ОДНОЙ СИСТЕМЫ ГИПЕРБОЛИЧЕСКИХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ТРЕТЬЕГО ПОРЯДКА С ДВУМЯ НЕЗАВИСИМЫМИ ПЕРЕМЕННЫМИ
А. А. Андреев, Ю. О. Яковлева
Самарский государственный технический университет,
443100, Самара, ул. Молодогвардейская, 244.
E-mail: julia.yakovleva@mail.ru
На основе метода Римана получено решение задачи Гурса для системы дифференциальных уравнений третьего порядка. Получена матрица Римана, выраженная через гипергеометрические функции матричного аргумента, с помощью которой найдено решение задачи Гурса для системы линейных гиперболических уравнений третьего порядка.
Ключевые слова: система гиперболических дифференциальных уравнений третьего порядка, сопряжённое уравнение, задача Гурса, метод Римана, матрица Римана.
Введение. Известно [1], что основной вклад в начало современной теории гиперболических уравнений в частных производных был внесён Б. Риманом, получившим в работе [2] интегральное представление задачи Коши в форме, аналогичной представлениям решений краевых задач для эллиптических уравнений второго порядка с помощью функций Грина.
Метод, применяемый Б. Риманом, предполагал существование вспомогательной функции, теперь называемой функции Римана, обладающей известными свойствами [1]. Функция Римана играет фундаментальную роль в теории линейных уравнений гиперболического типа и с её помощью удаётся, как правило, записать решение задач Коши и Гурса в явном виде (в квадратурах).
Однако и до сегодняшнего дня известен небольшой класс уравнений, для которых функция Римана может быть записана в замкнутом виде.
Метод Римана решения краевых задач для уравнений гиперболического типа второго порядка с двумя независимыми переменными обобщался на другие уравнения и системы. В монографиях А. В. Бицадзе [3] и И.Н. Векуа [4] приведено обобщение метода Римана на один класс гиперболических систем второго порядка с двумя независимыми переменными и кратными характеристиками, при этом показано, что вопрос о нахождении матрицы Римана редуцируется к решению системы интегральных уравнений Вольтерра второго рода, которая всегда имеет единственное решение.
Дальнейшему развитию метода Римана для уравнения гиперболического типа с двумя независимыми переменными были посвящены работы А. П. Солдатова, М.Х Шханукова [5], О. М. Джохадзе [6], О. С. Зикирова [7], где дано определение функции Римана, которая вводится как решение некоторой специальной задачи Гурса. Такой способ обеспечивает существование функ-
Александр Анатольевич Андреев (к.ф.-м.н., доц.), доцент, каф. прикладной математики и информатики. Юлия Олеговна Яковлева, аспирант, каф. прикладной математики и информатики.
ции Римана, но вопрос о ее явном построении остается открытым. В работах
В. И. Жегалова, А.Н. Миронова и Е. А. Уткиной [8,9] предлагается способ нахождения функции Римана как решения интегрального уравнения.
Известно, что с помощью интерполяционного многочлена Лагранжа-Сильвестра [10] можно определить значение аналитической функции на множестве постоянных квадратных матриц. Если ограничиться множеством матриц, являющихся значениями некоторых аналитических функций от одной матрицы, то определение легко обобщается на случай аналитических функций многих комплексных переменных, что позволяет, в свою очередь, доопределять целый ряд специальных функций на матричные значения входящих в них параметров. Некоторые примеры приведены в работах [11,12], где также дано определение матрицы Римана, которое мы приведём ниже.
1. Постановка задачи. В односвязной области И = {(х,у) : 0 < х < 1, 0 < у < 1} независимых переменных (х,у) рассмотрим систему линейных уравнений третьего порядка в частных производных, не содержащую производные порядка меньше третьего:
где и(х,у) —искомая п-мерная вектор-функция, Л — постоянная квадратная матрица порядка п.
Регулярным в области И решением системы уравнений (1) называется вектор-функция и(х,у), имеющая в И все непрерывные частные производные, входящие в уравнение (1), и удовлетворяющая ей в обычном смысле.
Для системы уравнений (1) поставим следующую граничную задачу (задачу Гурса): найти регулярное в области И решение II (х, у) системы уравнений (1), удовлетворяющее условиям
где А(у), В (у), С(х) — заданные вектор-функции такие, что А(0) = С(0), С'{ 0) = 5(0).
где Р = У11ху + Ухуи, <5 = —УХ17Х1 на области Д) = {(х,у) : 0 < х < хо, 0 < у < уо} придём к следующему соотношению:
ми = ихху + Аи = 0,
(1)
и(х,у)\х=0 = А(у), 0^у^1 их(х,у)\х=0 = В(у), 0^у^1 и(х,у)\у=о = С(х), 0^ж^1
(2)
(3)
где Г = АДА)-36
Полагая, что II(ж, у)— решение системы уравнений (1), а V = V(хо,Уо',х, у) удовлетворяет матричному уравнению М*У = 0, после некоторых преобразований с криволинейными интегралами второго рода получим
/хо
Ух(хо,уо]х, у)их(х, у) |У=0<1х-
ГУО
- / {У(хо,Уо]х,у)иХу(х,у) + УХу(х0,уо]х,у)и(х,у))\ с!у+
ГУО
+ / {У(хо,Уо]х,у)иХу(х,у)+ УХу(х0,уо]х,у)и(х,у))\х=хос1у-
ГХ о
ГХ о
Ух(х0,у0]х,у)их(х,у)\у=уоёх = 0. (5)
2. Определение матрицы Римана. Матрицей Римана для уравнения (1) называется решение V = У(хо,уо',х,у) задачи
М*У = О
при выполнении условии
У(х0,уо]х,у)\х=хо = О,
У(хо,уо]х,у)\у=у0 = (х-хо)Е, (6)
Ух(Х0> 2/0) X, у) |ж=жо —
где (жо, уо) — произвольная точка области И, Е — единичная матрица порядка п.
Интегрируя по частям полученное равенство (5) и учитывая условия (2) и (6), получим
[■X о
/ Ух{х0,у0]х^)с1{х)йх-
- У(хо,уо]0,уо)их (0, уо) + У(хо,уо]0,0)их (0, 0) -гуо гуо
- / Уу(хо,уо]0,у)их(0,у)с1у - / УХу(хо,уо]0,у)А(у)с1у+
Jo J о
ГХ о
+ / их(х,у0)йх = 0,
Jo
ГХ о
Ух(хо,уо]Хо, О)С(жо) - Ух(хо,уо] 0,0)С(0) - / Ухх(х0,уо]х, 0)С(х)йх-
J о
ГУО
- У(хо,уо]0,уо)В(уо) + У(хо,уо]0,0)В(0) + / Уу(хо,уо]0,у)В(у)с1у-
ГУО
- / УХу(хо,уо] О, у)А(у)с1у - и(ж0, Уо) + А(у0) = 0. Вектор-функция, являющаяся решением уравнения (1), имеет вид
ГУО
и(х0,уо) = Ух(хо,уо] О, О)А(О) + / Ух(х0,уо] 0, у)А'(у)с1у-
•УО
гуо гхо
- / У(х0,у0]0,у)В,(у)йу + / Ух(хо,уо-,х,0)С'(х)йх. (7)
.)о .)о
3. Построение матрицы Римана. Введём новые переменные I = х — хо, в =
= у — уо- Матрицей Римана У{Ь,в) является решение матричного уравнения
Уш -УА = 0, (8)
удовлетворяющее условиям
У(0,в) = 0, У(1,0)=1Е, У(0,з)=Е. (9)
Решение задачи (8), (9) будем искать в виде мультипликативного пре-
образования V = £И^(<7), где <7 = Матричное уравнение (8) при этом преобразуется к виду
А(т2\У"'((т) + 1А(т\У" + 6И^ - ШК = 0 (10)
при \¥(0) = Е
а ___
(1(7
Пусть 5 = а 4- [13], тогда
г12
2 а — г2 г а г- = 6 — й,
(1<т2
<т3-^ = £3-3£2+2£. (11)
аа6
Подставляя (11) в (10), получим
(4(£3 - 3£2 + 25) + 14(£2 - 5) + 65) \У = а\УЛ,
(^462 (5 + = 0. (12)
Ищем решение матричного уравнения (12) в виде (о) = ^т=о АпСг"г, где Ап — постоянные квадратные матрицы порядка п. Так как
4т2 + 0 Дп = Лп_ь Л = Я,
получаем следующую рекуррентную формулу:
Ат =-----------------Ат.
(1)т(§)тт!22т
Следовательно,
ОО
(7\т
“'п (1)™(|)™.т1 VI
т=0
Пользуясь определением обобщенной гипергеометрической функции [13,14], получим матрицу Римана:
У(х0,уо]х,у) = (ж-ж0) 0^2 ^1; (0 Л) • (13)
При подстановке (13) в (7) после некоторых преобразований решение задачи Гурса уравнения (1) принимает вид
и(х0,уо) = ^0^2(1; ^;т0л) + ^00^2(2; А(о)+
+ ^ ^^(іфтіЛ) + |п0^2(2; 7^тіЛ)) А'(у)г1у+
+ іо Х° (°І?2(1ФТіЛ)) В'(№у+
+ ^ ^0^2(і; ^;т2Л) + ^7-2 0^2^2; ^;т2л)^ С{х)с1х, (14)
где т0 = —0,25жоУо, ті = 0,25ж^(у - у0), т2 = -0,25(ж - х0)2уо-
С помощью интерполяционного многочлена Лагранжа-Сильвестра [10,14] можно определить значение функции матричного аргумента на множестве квадратных матриц следующим образом:
/ 3 \ 8 тк
акз(А - \кЕу~1ф\к),
к=1j=1
где фк(Х) = (лзх^г, Які = (^["У(ЛУ7]Л"=Л\ З = 1,2,...,тк, к = = 1,2,...,8.
Непосредственной проверкой убеждаемся, что (14) удовлетворяет задаче (1), (2). Таким образом, резюмируя вышеизложенное, основные результаты можно сформулировать в виде теоремы. Теорема. Если А(у), В (у), С{ ж) Є С1 (/), I = (0,1), то существует единственное решение задачи Гурса (1), (2), которое имеет вид (14).
fb(l; f ;Л)
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Courant R., Hilbert D. Methods of mathematical physics. Vol. II: Partial differential equations. New York - London: Interscience Publishers, 1962. 830 pp.; русск. пер.: Курант P. Уравнения с частными производными. М.: Мир, 1964. 831 с.
2. Riemann В. Ueber die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite // Abh. d. Gott. Ges. d. Wiss., 1860. Vol. 8. Pp. 43-68.
3. Бицадзе А. В. Некоторые классы уравнений в частных производных. М.: Наука, 1981. 448 с. [Bitsadze А. V. Some classes of partial differential equations. Moscow: Nauka, 1981. 448 pp.]
4. Векуа И. H. Новые методы решения эллиптических уравнений. M.-JL: ОГИЗ, 1948. 296 с. [Vekua I. N. New Methods for Solving Elliptic Equations. Moscow-Leningrad: OGIZ, 1948. 296 pp.]
5. Солдатов А. П., Шхануков М. X. Краевые задачи с общим нелокальным условием А. А. Самарского для псевдопараболических уравнений высокого порядка // ДАН СССР, 1987. Т. 297, №3. С. 547-552; англ. пер.: Soldatov А. P., Shkhanukov М. Kh. Boundary value problems with A. A. Samarskii’s general nonlocal condition for higher-order pseudoparabolic equations// Soviet Math. Dokl., 1988. Vol. 36, no. 3. Pp. 507-511.
6. Джохадзе О. М. Функция Римана для гиперболических уравнений и систем высокого порядка с доминированными младшими членами // Дифференц. уравн., 2003. Т. 39, №10. С. 1366-1378; англ. пер.: Dzhokhadze О. М. The Riemann function for higher-order hyperbolic equations and systems with dominated lower terms // Differ. Equ., 2003. Vol. 39, no. 10. Pp. 1440-1453.
7. Зикиров О. С. Локальные и нелокальные краевые задачи для гиперболических уравнений третьего порядка// Соврем, мат. прилож., 2011. Т. 68. С. 101-120; англ. пер.: Zikirov О. S. Local and nonlocal boundary-value problems for third-order hyperbolic equations// J. Math. Sci.. Vol. 175, no. 1. Pp. 104-123.
8. Жегалов В. И., Миронов А. Н. О задачах Коши для двух уравнений в частных производных// Изв. вузов. Матем., 2002. №5. С. 23-30; англ. пер.: Zhegalov V.I., Mironov А. N. Cauchy problems for two partial differential equations // Russian Math. (Iz. VUZ), 2002. Vol. 46, no. 5. Pp. 21-28.
9. Жегалов В. И., Уткина Е. А. Об одном псевдопараболическом уравнении третьего порядка // Изв. вузов. Матем., 1999. № 10. С. 73-76; англ. пер.: Zhegalov V. Utkina Е. А. Pseudoparabolic equation of the third order// Russian Math. (Iz. VUZ), 1999. Vol. 43, no. 10. Pp. 70-73.
10. Гантмахер Ф. P. Теория матриц. М.: Наука, 1988. 549 с. [Cantmakher F. R. Theory of matrices. Moscow: Nauka. 549 pp.]
11. Андреев А. А. Построение элементарных решений и решение задачи Коши для уравнений и систем уравнений гиперболического типа: Дис. ... канд. физ.-мат. наук. Куйбышев, 1981. 100 с. [Andreev A. A. Construction of elementary solutions and solution of Cauchy problem for equations and hyperbolic systems of equations: Ph. D. Thesis (Phys. & Math.). Kuibyshev, 1981. 100 pp.]
12. Андреев А. А. О некоторых приложениях ассоциированных гипергеометрических функций/ В сб.: Дифференциальные уравнения (математическая физика): Тез. докл. участников Куйбышевского облает, межвуз. науч. совещания-семинара. Куйбышев, 1984. С. 8-9. [Andreev A. A. On some applications of associated hypergeometric functions / In: Differential Equations (Mathematical Physics). Kuibyshev, 1984. Pp. 8-9].
13. Erdelyi A., Magnus W., Oberhettinger F., Tricomi F. C. Higher transcendental functions. Vol. I / ed. H. Bateman. New York - Toronto - London: McGraw-Hill Book Co, Inc., 1953. 302 pp.; русск. пер.: Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. В 3-х т. Т. 1: Гипергеометрическая функция. Функция Лежандра. М.: Наука, 1973. 296 с.
14. Андреев А. А, Огородников Е. Н. Матричные интегро-дифференциальные операторы и их применение// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 1999. №7. С. 27-37. [Andreev A. A., Ogorodnikov Е. N. Matrix integro-differential operators and their application// Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 1999. no. 7. Pp. 27-37].
15. Кошляков H. С., Глинер Э. Б., Смирнов М. М. Дифференциальные уравнения математической физики. М.: Государств, изд-во физ.-мат. лит., 1962. 767 с. [Koshlyakov N.S., Cliner Е. В., Smirnov М. М. The differential equations of mathematical physics. Moscow: Gosudarstv. Izd-vo. Fiz.-Mat. Lit., 1962. 767 pp.]
Поступила в редакцию ll/V/2011; в окончательном варианте — 11 /VIII/2011.
MSC: 35L35
THE GOURSAT PROBLEM FOR ONE HYPERBOLIC SYSTEM OF THE THIRD ORDER DIFFERENTIAL EQUATIONS WITH TWO INDEPENDENT VARIABLES
A. A. Andreev, J. O. Yakovleva
Samara State Technical University,
244, Molodogvardeyskaya St., Samara, 443100, Russia.
E-mail: julia.yakovleva@mail.ru
The solution of the Goursat problem for the system of the differential equations of the third order is received on the basis of Riemann method. The Riemann matrix expressed in terms of hypergeometric functions with matrix argument is received, using it the solution of Goursat problem for the system of the linear hyperbolic equations of the third order is constructed.
Key words: hyperbolic system of the third order differential equations, adjoint equation, Goursat problem, Riemann method, Riemann matrix.
Original article submitted ll/V/2011; revision submitted ll/VIII/2011.
Aleksander A. Andreev (Ph.D. (Phys. & Math.)), Associate Professor, Dept, of Applied Mathematics & Computer Science. Julia O. Yakovleva, Postgraduate Student, Dept, of Applied Mathematics & Computer Science.
