УДК 517.968.78
ОБ ОДНОЙ КРАЕВОЙ ЗАДАЧЕ ДЛЯ СИСТЕМЫ ГИПЕРБОЛИЧЕСКИХ УРАВНЕНИЙ С ВОЛНОВЫМ ОПЕРАТОРОМ И СИНГУЛЯРНЫМ КОЭФФИЦИЕНТОМ ПРИ МЛАДШЕЙ ПРОИЗВОДНОЙ
Р. Р. Раянова
Самарский государственный технический университет,
Россия, 443100, Самара, ул. Молодогвардейская, 244.
E-mail: г ay anova. г ina@gmail. com
Рассмотрена краевая задача в характеристическом квадрате с данными на параллельных характеристиках для. системы гиперболических уравнений с волновым оператором и сингулярным матричным коэффициентом при младшей производной. В характеристических координатах указанная система дифференциальных уравнений редуцируется к системе уравнений Эйлера—Пуассона—Дар-бу. С использованием известного решения задачи Коши с данными на линии сингулярности матричного коэффициента задача редуцируется к системе интегральных уравнений Карлемана. На основе проведённых ранее автором данной статьи исследований по разрешимости систем обобщённых интегральных уравнений Абеля в 'работе найдено в явном виде решение указанной краевой задачи.
Ключевые слова: дробное исчисление Римана—Лиувилля, функции от матриц, интегро-дифференциальные операторы матричного порядка, система обобщенных интегральных уравнений Абеля, интегральное уравнение Карлемана.
Обозначим через Мп — множество постоянных матриц порядка п. Пусть А(Сг)—спектр матрицы О Є Мга; Лі, г = 1,2— собственные значения матрицы О.
В области И = {(ж, у) : 0 < х + у < 1, 0 < х — у < 1} рассмотрим систему уравнений
20
Ихх ^уу ~ У'у — 0) (1)
где и(х,у) = («і(ж, у), «г(ж, у),..., ип(х, у))Т — вектор искомых функций,
Л(С) Є (0,1/2).
Пусть Д) = -О П {у > 0} = {(х, у) : 0 < х + у < 1,0 < х — у < 1,у > 0}, £>і = Б П {у < 0} = {(ж,у) : 0 < х + у < 1,0 < х — у < 1,у < 0}.
Обозначим (ж/2,ж/2) и $і((1 + ж)/2, (ж — 1)/2) —точки пересечения характеристик х — у = 0 и х — у = 1 с характеристикой другого семейства, выходящей из точки (ж, 0).
Задача 1. Найти вектор-функцию и(х, у), удовлетворяющую следующим свойствам:
1) и(х, у) Є С (В) П С2(Д) и £>і),
2) и(х,у) удовлетворяет системе (1) в области Д и Д,
3) и(во) = ір(ж), ж Є [0,1],
4) и(ві) = ф(х), х Є [0,1],
5) Дш (-у)Жди=^у2Сди}
Рина Ринатовна Раянова, ассистент, каф. прикладной математики и информатики.
где ір(х) = (ірі(х),ір2(х),...,ірп(х)), ф{х) = (фі(х),ф2(х),...,фп(х)), у2С -показательная функция матричного аргумента [1].
Нетрудно показать, что в характеристических координатах £ = х — у и г] = х + у область _Д преобразуется в область Н = {(£, ?у) :0 <{ < г] < 1}, а система дифференциальных уравнений (1) редуцируется к системе дифференциальных уравнений Эйлера—Пуссона—Дарбу следующего вида:
щ - = 0. (2)
Г]~(
В характеристических координатах £ = х+у и г] = х—у область Д преобразуется в ту же самую область Н, а система дифференциальных уравнений (1) — в систему дифференциальных уравнений Эйлера—Пуссона—Дарбу (2).
Для решения поставленной задачи воспользуемся решением задачи Коши с данными
Ит
57-?->0+
и(£,0=т(0, С Є [0,1], (з)
(^)2%, - Щ) = КО, ( є (0,1), (4)
где т, V — заданные п-мерные вектор-функции.
Решение уравнения Эйлера—Пуссона—Дарбу в скалярном случае [2] нетрудно обобщить на матричный. Справедливо следующее утверждение.
Утверждение 1. Пусть спектр матрицы О принадлежит интервалу (0,1/2). Тогда регулярное в характеристической области Н решение задачи Коши (3), (4) может быть записано в виде
«(£, г?) = *4 (со (г? - оЕ~2С !\п -т- о\-{Е-с)тт+
+ к2(С)22С~Е ^ [(?7 — £)(£ — £)]_сг/(£)сЙ, (5)
гдекг(С) = Г(2С)Г“2(С), к2{С) = Т{Е-2а)Т~2{Е-О), т(£),К0 еС2(0,1),
Е — единичная матрица, причём функция г/({) на концах интервала (0,1) может иметь интегрируемые особенности.
Располагая решением задачи Коши (5) для систем уравнений (2), можно показать, что решение системы уравнений (1) в области Д с данными 3) и
1 • 2 С / \
пт у — = 1У+(х), у^о+у ду ь
и в области Д — с данными 4) и
1 • 2 С / \
1шп[ У 7Г = 1У-(Х)’
у^о- ду
соответственно, имеют вид
и+(х,у)=Т 1(С) [ [в(1-в)] {Е с)(х-у + 2уз)Е °х X 0^х_у+2уз(х -у + 2У8)2°-Е(р(х -у + 2уз№+
+ Т~2(Е - С)Т(Е - 2 С)у2С~Е [ [в(1 - 8)}~си+(х -у + 2уз)й8-
Jo
- Т-1(С)Т(Е - 20)Т~1(Е - О)22С~ех
X [ № ~ *)}~{Е~С) 1о,х-у+2у8''+(х-У + 2у8)<18', (6)
и-(х,у) = -Т~1(0) [ [в(1 -в)]_(в_с)(1 -х + у-2у8)Е-°х ■)о Х Вх-у+2у3,1(1 -Х + у- 2у8)2С-Еф{х -у + 2уз)<18--Т~2(Е -С)Т(Е -2С){-у)е~2С [ [8{1-8)]-си_{х-у + 2у8)й8+ ■) о + Т-1(С)Т(Е - 20)Т~1(Е - О)22С~ех
х [ ^-{х-у+ 2У8)(18. (7)
Здесь (формулы (6) И (7)) И далее предполагается, ЧТО функции (Рг(х), фг(х), щ(х), г = 1,2,..., гг, удовлетворяют следующим условиям:
щ{х) € С[0,1], щ{х) € С1(0,1), г/'(ж) € ЦО, 1), щ(0) = щ(1) = 0;
^(х)еС[ 0,1], <£г(х)еС2( 0,1), ^'(ж)еД0,1), (^г(О) = <а(1) = 0;
фг(х)еС[0,1], фг(х)еС2(0,1), ф"(х) € ЦО, 1), фг(0) = фг(1) = 0.
Найдём значения функций и+(х,у) и ь,-(х,у) на границе у = 0:
и+(х, 0) = Т(С)Т~1(2С)хе~°0$+х2С~е<£(х) —
- Т~\Е - С)Т(Е - 2С)Т(С)Т-1(2С)22С-е 1Е~2Си+(х), (8)
и-(х, 0) = -Г(С)Г“1(2С)(1 -ж)е“с1^_(1 -х)2а~Еф{х)+
+ Т~\Е - С)Т(Е - 2С)Т(С)Т-1(2С)22С~Е 1Е-2Си_(х). (9)
По условию задачи предполагается, что искомая вектор-функция и(х, у) € 0(0). Это означает, что и+(х, 0) = и-(х,0) = и(х, 0). Приравнивая правые части равенств (8) и (9), приходим к системе интегральных уравнений типа системы обобщённых уравнений Абеля:
11(С)1Е~2Си(х) + 1\(С)^~2С:и(х) = = 12(С)хЕ-°в£+х2С~Еф) - 12(С)( 1 - х)Е~сВ?_(1 - х)2С~Еф(х),
где 1г(С) = Т-\Е - С)Т(Е - 2С)Г(С)Г“1 (2С)22С~Е и 12{С) = Г(С)Г“1(2С).
Полученная система интегральных уравнений сводится к системе интегральных уравнений Карлемана. Действительно, записывая интегральные операторы 1Е~2С и 1Е_Г2 по определениям левостроннего и правостороннего матричного интегрального оператора Римана—Лиувиля [3,4], получим
/(х) = 12(С)хЕ СД^_Ж2С е ір(х) — 12(С)(1 — х)е с£>(р_(1—ж)2С Еф(х).
д(х) =Т(Е — С)2е~2С (хе~°0%+х2С~е<р(х) - (1 - х)Е~с(1 - х)2С~Еф(х)).
Используя результаты работы [5], решение системы интегральных уравнений (10) можно найти как частный случай решения системы обобщённых интегральных уравнений Абеля
при А = В = Е.
Пусть О = [а, 6], где — оо < а < Ъ < +оо. Будем говорить, что вектор-функция /(ж) = (/і(ж), /2(ж),..., /п(х)) с областью определения Df = О, удовлетворяет условию Гёльдера порядка Л = (Аі, Аг,..., Ага), если
\/жі,ж2єП ІЛ(жі) - Д(ж2)| ^ Ак\хі - х2\Хк, к = 1,2,... ,п,
где Ак — постоянные, а А& — показатель Гёльдера.
Определение 1. Через Нх = НХ(П) обозначим класс вектор-функций /(ж) = (/і(ж), /2(ж),..., /п(ж)) с областью определения Df = О, каждая компонента которых удовлетворяет на О условию Гёльдера со своими фиксированными значениями Ак и А^.
Определение 2. Через Н* = Н*(а,Ь) обозначим класс вектор-функций /(ж) = (/і(ж), /2(ж),..., /п(ж)), Для которых существуют такие мультииндексы А = (Аі, А2,..., Ага), є = (єі, є2, • • •, єп) и ё = (6і,62, ■■■,6П), что
где /*(ж) Є Нх([а, Ь]) и каждой компоненте вектора /(ж) соответствуют свои фиксированные значения мультииндексов А, є и 5.
Определение 3. Пусть а = (скі, ск2,..., ага) — мультииндекс. Через Я* обозначим класс вектор-функций /(ж) = (/і(ж), /2(ж),..., /п(ж)) таких, что
іі(С)Г 1(_Б — 2С) Г / (ж — і) 2Сг/(і)сЙ + ! (і — ж) 2Си(і)сІі \ = /(ж)
где
В итоге имеем
(10)
где
л/йг2^ + д/^-2^ = ^)
(ж-а^Ь-ж)1-5’
/(ж)
/*0*0
(ж — а)1_а_е(6 — ж)1_а_<5 ’
где 0 < Єк < 1 — (У-к-, 0 < 5к < 1 — (У-к (к = 1, 2,..., п), а каждая компонента вектора /*(ж) принадлежит своєму классу Hak, который определён в [6].
Справедлива следующая теорема.
Теорема. Пусть матрица G Є Мп, спектр Л(С) С (0,1/2). Пусть матрица Т Є Мп является матрицей преобразования G к жордановому виду Ас = TGT~l. Пусть вектор Тд(ж) Є Н*, где а = (а\, а2, ■ ■ ■, ап), а,к = = 1 — 2 А к, А к Є Л (G). Тогда единственное решение системы интегральных уравнений (10) в классе вектор-функций, таких, чтоТи(х) Є Н*, имеет вид
Подставляя найденное в формуле (11) выражение для вектора v (ж) в формулы (6), (7), найдём окончательно выражения искомых вектор-функций и+(х,у) и v,-(x,y), определяющих решение задачи 1 в областях Do и D\ соответственно, а следовательно, и решение и(х,у) в области D = Do U D\.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Ф. Р. Гантмахер, Теория матриц. М.: Наука, 1988. 549 с. [F. R. Gantmakher, Theory of matrices. Moscow: Nauka, 1988. 549 pp.]
2. А. А. Андреев, “О методе Римана для одной системы уравнений гиперболического типа с кратными характеристиками” / В сб.: Корректные краевые задачи для неклассических уравнений математической физики. Новосибирск: ИМ СОАН СССР, 1981. С. 13-16. [A A. Andreev, “On the Riemann method for one system hyperbolic equations with multiple characteristics” / In: Correct Boundary Value Problems for Nonclassical Equations of Mathematical Physics. Novosibirsk: IM SOAN SSSR, 1981. Pp. 13-16].
3. А. А. Андреев, “Нелокальные краевые задачи для одной модельной вырождающейся системы гиперболического типах” / В сб.: Краевые задачи для уравнений математической физики. Куйбышев: Куйбыш. гос. пед. ин-т, 1990. С. 3-6. [A A. Andreev, “Nonlocal boundary value problems for a model degenerate system of hyperbolic type” / In: Boundary value problems for equations in mathematical physics. Kuybyshev: Kuybyshev. Cos. Ped. Inst., 1990. Pp. 3-6].
4. А. А. Андреев, “Об одном обобщении операторов дробного интегро-дифференцирования и его приложениях”/ В сб.: Интегральные уравнения и краевые задачи математической физики: Матер. Всессоюзной конф.. Владивосток, 1990. С. 91. [A. A. Andreev, “On one generalization of fractional integro-differentiation operators and its applications” / In: Integral Equations and Boundary Value Problems of Mathematical Physics. Vladivostok, 1990. Pp. 91].
5. P. P. Исмагилова, “Решение полного матричного аналога обобщённого уравнения Абеля с постоянными коэффициентами”// Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2011. №1(22). С. 93-98. [Л. R. Ismagiliva, “The solution of the full matrix analogue of the generalized Abel equation with constant coefficients” // Vestn. Samar. Cos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki, 2011. no. 1(22). Pp. 93-98].
ф) = (ЩЕ - 2G))" V+-2G(£ - Z(x)Dlp^Z~l(x))g(x), (11)
где
6. С. Г. Самко, А. А. Килбас, О. И. Маричев, Интегралы и производные дробного порядка и некоторые их приложения. Минск: Наука и техника, 1987. 688 с. [5. С. Samko, A. A. Kilbas, О. I. Marichev, Integrals and derivatives of fractional order and some of their applications. Minsk: Nauka i Tekhnika, 1987. 688 pp.]
Поступила в редакцию 29/X/2012; в окончательном варианте — 27/1/2013.
MSC: 35L80
ON A BOUNDARY VALUE PROBLEM FOR A SYSTEM OF HYPERBOLIC EQUATIONS WITH THE WAVE OPERATOR AND A SINGULAR COEFFICIENT AT LOWER DERIVATIVE
R. R. Rayanova
Samara State Technical University,
244, Molodogvardeyskaya St., Samara, 443100, Russia.
E-mail: rayanova. rinaSgmail. com
The boundary value problem with data given on the parallel characteristics for the system of hyperbolic equations with the wave operator and the singular matrix coefficient at the lower derivative is considered in the characteristic square. This system of differential equations in the characteristic coordinates can be reduced to the system of Euler-Poisson-Darboux equations. Using the known solution of Cauchy problem with data given on the singularity line of matrix coefficient, we reduce the problem to the Carleman system of integral equations. The explicit solution of the considered boundary value problem is constructed using the results of previous research on the solvability of the systems of generalized Abel integral equations, made by the author.
Key words: Riemann-Liouville fractional calculation, matrix functions, integral-differentional operators of matrix order, system of generalized Abel integral equations, Carleman integral equation.
Original article submitted 29/X/2012; revision submitted 27/1/2013.
Rina R. Rayanova, Assistent, Dept, of Applied Mathematics & Computer Science.