CC BY
33
УДК 517.929.7
ИССЛЕДОВАНИЕ РАЗРЕШИМОСТИ ВТОРОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ В ЧАСТНЫХ ПРОИЗВОДНЫХ С ИНВОЛЮТИВНЫМ ОТКЛОНЕНИЕМ В МЛАДШИХ ЧЛЕНАХ
Бжеумихова Оксана Игоревна Кабардино-Балкарский государственный университет им. Х.М. Бербекова, Нальчик, Россия
В работе исследован вопрос разрешимости второй краевой задачи для модельного уравнения в частных производных с инволютивным отклонением в младших членах. Исследование проведено на основе метода разделения переменных
UDC 517.929.7
THE INVESTIGATION SOLVABILITY OF THE SECOND BOUNDARY VALUE PROBLEM FOR PARTIAL DIFFERENTIAL EQUATION WITH INVOLUTORY DEVIATIONS IN THE LOWEST TERMS
Bzheumikhova Oksana Igorevna Kabardino-Balkarian State University named after H.M. Berbekova, Nalchik, Russia
In this article we consider the problem of solvability oа second boundary value problem for the model equation in partial derivatives with involutive deviation in the lowest terms. The investigation is based on a variable separation method
Ключевые слова: КРАЕВАЯ ЗАДАЧА, ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ В ЧАСТНЫХ ПРОИЗВОДНЫХ, ОТКЛОНЯЮЩИЙСЯ АРГУМЕНТ, ЗАДАЧА ШТУРМА-ЛИУВИЛЛЯ
Keywords: BOUNDARY VALUE PROBLEM, PARTIAL DIFFERENTIAL EQUATION, DEVIATIONS ARGUMENT, STURM-LIOUVILLE PROBLEM
Введение
Многие математические модели, применяемые при исследовании процессов, в таких важных областях как математическая биоэкология, механика, автоматизированные системы управления, теория климатических моделей, иммунология и т. д. базируются на дифференциальных уравнениях с отклоняющимся аргументом (например, [1-6]). Широкие возможности применения уравнений с отклоняющимся аргументом в качестве математических моделей способствуют росту интереса к исследованию новых задач для уравнений с частными производными [7-10], которые по сравнению с обыкновенными дифференциальными уравнениями описывают процессы еще в большей степени приближенные к процессам, протекающим на практике [11, 12].
В настоящей работе, методом разделения переменных, установлена разрешимость классической краевой задачи для уравнения в частных
производных с инволютивным отклонением аргумента в прямоугольной области.
1. Постановка задачи
Пусть О = {(х, t): -х0 < х < х0, 0 < t < t0} - односвязная область
евклидовой плоскости R2 точек (х, t).
В области О рассмотрим уравнение
Lu °uхх(х,t) + k^) utt(х,t)+ u(-х,t) = 0, (1)
где k^) - достаточно гладкая, причем k^) Ф 0.
Для уравнения (1) исследована следующая
Задача 1. Найти регулярное в области О решение u(х, t) уравнения (1) из класса C1 (о) I С^2, удовлетворяющее условиям
uх (- х0, t )=Р( \ (х0, t )=Р ( ), (2)
Щ (х,0)=^з (х), Щ (х, t0 )=^4 (х),
где р ( = 1,4) - заданные, достаточно гладкие функций.
2. Доказательство существования и единственности задачи
Для задачи 1 справедлива следующая Теорема 1. Пусть
1) р1 (),р2()е С3[0,t0], р3(х),р4(х)е С3{-х0,0)и (0,х0]}, где х0 е
А _ Л 0,
to(0,р), ^)е С1 [0,to];
х0 х0
2) Ь (х) * = \р4 (х) * = °’
0 0
тогда задача (1), (2) разрешима в требуемом классе функций.
Действительно, разобьем задачу (1), (2) в области W на две вспомогательные:
Lu1 ° 0, (З)
ulx (- x0, t )= 0, ulx (x0, t ) = 0, (4)
ult (x,0)= Рз (x), ult (x, t0 )=P4 (x) , (5)
Lu2 ° 0, (б)
u2t (x,0)= 0, u2t t0 )= 0, (7)
U2x (- X0, t )=Pl (t), U2x (Xbt )=P2 (t) , (8)
где u = u1 + u2.
Решение уравнения (З) удовлетворяющее однородным граничным
условиям (4) будем искать в виде [1З]:
u, = Х1 (x )• Tl (t). (9)
Подставляя (9) в (З) и опуская нижние индексы, получим
Х *(x) + Х (- x) = - ҐЦ) = -Я Х (x) = T (t) = ,
где Я = const.
Отсюда, с учетом (4) будем иметь
Х "(x) + Х (- x ) + ЯХ (x ) = 0, (10)
Х (-xo) = Х (x0 ) = 0, (11)
k (t )T "(t)-IT (t ) = 0. (12)
Исследуем задачу о собственных значениях (10), (11).
Дважды дифференцируя (10), приходим к соотношению:
Х1]/ (x) + Х'(- x) + ЯХ' (x )= 0. (1З)
С другой стороны из (10) имеем:
Х'(- x ) = - Х (x)-ЯХ (- x), (14)
Х (- x ) = - Х' (x)- ЯХ (x). ( )
На основании (1З) и принимая во внимание (14), получим
Х1]/ (х) + 21Х"(х )+(12 - 1)х(х ) = 0. (15)
Характеристическое уравнение соответствующее (15), будет иметь
вид:
г4 + 21 г2 + 12 -1 = 0.
Разрешая биквадратное уравнение, находим:
гх = л/1 -1, г2 = -л/1 -1,
г3 = л/-1 -1, г4 = -V-1 -1.
Таким образом, общее решение уравнения (15) может быть записано в виде:
X(х) = С]е^~1х + С2е+ С3^^~1~1х + С4в~у[~1~1х. (16)
Следуя [14, 15], получим из (16) представления решения (10) для различных 1.
Случай 1: 1 < -1. В этом случае общее решение (10) имеет вид:
X (х) = С1 sh (хл/Т-1)+ С3 ch (хл/ -1 -1).
Используя условия (11), получим
^л/Г-1С1 сЪ[х0у1-Х)-4-\-1С3 sh(х0л1 -1 -1 )= 0, л(х0л/1-1)+л/-1 -1С3 sh (х0 V -1 -1 )= 0.
Определитель этой системы
А = 2л1 Л2 -1 ch(х0^/^^л)sh(х0л/-1 -1 )= 0 только при 1 = ±1, что противоречит рассматриваемому случаю 1 < -1. Следовательно, С = С3 = 0. Откуда заключаем, что X (х) ° 0.
Случай 2: 1 = -1. При таком значении 1 решение (10) имеет вид:
X (х) = С1 sh(л/2х)+С3.
Удовлетворяя (11), имеем
72^(72х0 )= 0,
■ДС^Лх0 )= 0.
Откуда заключаем, что С1 = 0 и X (х) ° С3.
Этому собственному значению 1 = -1 соответствует
и1 (х, t ) = С3Т ^), (17)
где Т ^) - решение уравнения (12).
Требуя выполнения граничных условий (5) , получаем систему для определения постоянных входящих в (17):
Рз(х ) = ии (х,0) = С3Т (0) ,
Р4 (х) = (х, t0 )= С3Т'(t0 ) .
Таким образом, решение задачи (3)-(5) при 1= -1 определяется соотношением (17).
Случай 3: Для -1 <1< 1, удовлетворяя общее решение (10)
X (х) = С1 sh (ха/ 1 -1)+ С3 cos(хл/1 +1) условиям (11), находим
л/1 - 1С1 ch(х0^Jг—л)+^Jг+Лcз sin(х0^/l + 1)= 0, л/1 - 1С1 ch(х0^|Г—л)-y|Г+Лcз sin(х0л/ТЛ)= 0.
Определитель системы
А = -2лА -12 л(х0л/1 - 1^т(х0л/1 +1)
обращается в нуль либо при 1 = ±1, либо при 1 Следовательно, С1 = С3 = 0 и X(х) ° 0, т.к. 1 £ (-1, 1).
Случай 4: 1 = 1. При указанном значении 1 для всех х0 е (10) принимает вид
X (х ) = С2 х + С3 cos (л/2х).
Удовлетворяя (18) граничным условиям (11) получим
с2 +V2c3 c2 -V2c3
2
pn
V х0 У
1.
V
У
В силу того, что
А = -2л/2х0 sin (л/2х0 )^ 0,
имеем С2 = С3 = 0 и X (х) ° 0.
Случай 5. При 1> 1 общее решение (10) принимает вид:
X(х) = С2 sin(хл/1-1)+С3 cos(W 1 +1).
Удовлетворяя полученное выражение для X (х) граничным условиям (11), будем иметь:
лИС2 cos (х0л/Х-Т )+л/1 + 1С3 sin ^лИ )= 0,
лИс2 cos (х0л/Х-1 ^лИС3 sin ^лИ )= 0.
Равенство
А = -2л/1 -1 cos^ л/ 1 -1 )sin (х^л1+1) = 0
справедливо при 11n
2
т + 2т
2 х,
2
pn
0
V х0 у
1.
+ 1, либо при 12г
Таким образом, задача (10), (11) имеет собственные значения , 1
и соответствующие им собственные функции X1n (х) = Cn cos
2 х,
х
0 У
X2n (х) = Cn sin
Г \
pn
—х
V х0 у
, (n е N), где Cn
произвольные постоянные,
нуждающиеся в определении.
Собственным значениям 11n уравнения (12) равные
f Л 2
т + 2т
V 2 х 0 у
1 соответствуют решения
T (t )= T„ (t).
Возвращаясь к решению задачи (3)-(5), видим, что функция
и
11(^ t) = Z CnTn (t)cos
n=1
2х
(19)
0у
является решением уравнения (3) при 1
/ \ 2 т + 2т
1п
2х
1.
0 У
(19).
Условия (5) позволяют определить значение постоянных входящих в
С учетом условия 1) теоремы 1, функции р3 (х) и р4 (х), - х0 < х < х
разлагаются в ряд Фурье, который содержит только косинусы, а именно:
&
где
Р3 (х + Е5п C0s
2
п=1
гж + 2лп Л
----------х
2х
V 0 у
Р4 (х ):
Р0
2 п=1
+ X Рпcos
гж + 2т Л х
2х
0 У
2 х0 2 х0
&0 = — | р3 (х )йх, р0 = — | р4 (х )йх
0 0
2 х0
&п = — р3 (х ^
00
х
Рп = — р4 (х ^ х
00
х0
2х
V 0 у
йх.
00
-------х
2х
V 0 у
йх.
причем ряды Е & и XI 8п\ сходятся.
п=1
п=1
Учитывая граничные условия (5), получаем:
Е СТ (0)cos
п=1
¥
X СпТ'п (t)c0s
п=1
т + 2т
---------х
V 2 х0 у
т + 2т
---------х
2х
V 0 у
¥
&+Х&cos
2 п=1
т + 2т Л
х
=Р+Ер
2
cos
V 2 х0 у
гт + 2т Л
х
п=1
2х
0у
(20)
Сопоставляя соответствующие коэффициенты в полученных соотношениях, а так же учитывая условие 2) теоремы 1 определяем постоянные входящие в (19). Следовательно, ряд (19) с коэффициентами
0
определяемыми по формулам (20), удовлетворяет всем условиям задачи
(3)-(5).
Переходя к рассмотрению случая собственных значений
С \2
т
V х0 у
1 будем иметь
и
12
(х, t) = Е СпТп ^^П
п=1
А Л
т
—х
V х0 у
(21)
Используя условия (5) позволяют определим значение постоянных входящих в (21).
С учетом условия 1) теоремы 1, функции р3 (х) и р4 (х), - х0 < х < х0 разлагаются в ряд Фурье, который содержит только синусы, а именно:
Р3(х )=Е тsin
п=1
г \ т
—х
V х0 у
Р4 (х )= ЕПп SІn
п=1
Г \
т
—х
V х0 у
где
2 х0 Г т ^ 2 х0
тп =— [р3 (х)sin — х йх, пп =— [р4(х)sin
у •/ -у» у •/
0 0 V 0 у Л0 0
г \ т —х
V х0 у
йх,
а ряды Е |т| и Е1 Пп| сходятся.
п =1
п=1
Учитывая граничные условия (5), получаем:
Е cлт;(0)sin
п=1
¥
Е СпТ'п ^)sІn
п=1
Г \
т
—х
V х0 у
г \ т —х
V х0 у
= Ет
Sin
п=1
Еп,
п=1
Sin
г \ т —х
V х0 у
г \ т —х
V х0 у
(22)
Сравнивая соответствующие коэффициенты в полученных соотношениях, а так же учитывая условие 2) теоремы 1 определяем постоянные входящие в (21). Представленные выше рассуждения остаются справедливыми и для случая задачи (6)-(8). Причем функция и2(х,t) аналогично функции и1 (х, t) для различных собственных значении
находится в виде сходящихся тригонометрических рядов. Таким образом, решение задачи 1 определяется из соотношения и(х, t ) = и1 (х, t)+ и2 (х, t).
Заключение
На основе метода разделения переменных было доказано существование решения второй краевой задачи для модельного уравнения в частных производных с инволютивным отклонением в младших членах. Несмотря на то, что результаты работы носят теоретический характер, они могут иметь широкое применение, как и в дальнейших исследованиях уравнений с отклоняющимся аргументом, так и в прикладных задачах.
Список литературы
1. Wu J. Theory and Applications of Partial Functional Differential Equations, Springer-Verlag, New York, 1996.
2. Wu J. Traveling wave fronts of reaction-diffusion systems with delay/ J.Wu, X. Zou // J. Dynamics and Differential Equations, 2001. - Vol. 13, No. 3. - P. 651-687.
3. Huang J. Traveling wave fronts in diffusive and cooperative Lotka-Volterra system with delays / J. Huang, X. Zou // J. Math. Anal. Appl, 2002. - Vol. 271. - P. 455-466.
4. Faria T. Nonmonotone travelling waves in a single species reaction-diffusion equation with delay / T. Faria, S. Trofimchuk // J. Differential Equations, 2006. - Vol. 228. - P. 357376.
5. Trofimchuk E. Slowly oscillating wave solutions of a single species reaction-diffusion equation with delay / E. Trofimchuk, V. Tkachenko, S. Trofimchuk // J. Differential Equations, 2008. - Vol. 245. - P. 2307-2332.
6. Meleshko S.V. On the complete group classification of the reaction-diffusion equation with a delay / S. V. Meleshko, S. Moyo // J. Math. Anal. Appl., 2008. - Vol. 338. - P. 448466.
7. Hernandez E. A note on partial functional differential equations with state-dependent delay / E. Hernandez, A. Prokopczyk, L. Ladeira // Nonlinear Analysis, R.W.A., 2006. - No.
4. - P. 510-519.
8. Rezounenko A.V. Stability of positive solutions of local partial differential equations with a nonlinear integral delay term / A.V. Rezounenko // Electronic Journal of Qualitative Theory of Differential Equations. Proc. 8th Coll. QTDE, 2008. - No. 17. - P. 1-7.
9. Bzheumikhova O.I. Application of Fourier method to investigation of the Dirichlet problem for partial differential equations with deviating arguments / O.I. Bzheumikhova, V.N. Lesev // International Journal of Differential Equations and Applications, 2013. - Vol.
12, No. 2. - P. 103-120.
10. Лесев В.Н. Об однозначной разрешимости задачи Неймана для эллиптического уравнения с отклоняющимся аргументом / В.Н. Лесев, О.И. Бжеумихова // Экологический вестник научных центров ЧЭС, 2012. - №3. - С. 41-46.
11. Wang L. Global exponential robust stability of reaction-diffusion interval neural networks with time-varying delays / L. Wang, Y. Gao // Physics Letters A, 2006. - Vol. 350.
- P. 342-348.
12. Lu J. G. Global exponential stability and periodicity of reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions / J. G. Lu. // Chaos, Solitons and Fractals, 2008. - Vol. 35. - P. 116-125.
13. Тихонов А.Н., Самарский А.А. Уравнения математической физики. - М.: Изд-во Наука, 1977. - 735 с.
14. Лесев В.Н. Применение метода Фурье к исследованию задачи Дирихле для уравнения с отклоняющимся аргументом и оператором Лапласа в главной части / В.Н. Лесев, О.И. Бжеумихова // Научный журнал КубГАУ [Электронный ресурс]. -Краснодар: КубГАУ, 2012. - №07(81). - С. 1-10.
15. Бжеумихова О.И. Краевые задачи для модельных уравнений смешанного типа второго порядка с отклоняющимся аргументом / О.И. Бжеумихова, В.Н. Лесев // Обозрение прикладной и промышленной математики, 2011. - Т. 18, вып. 5. - С. 744745.
References
1. Wu J. Theory and Applications of Partial Functional Differential Equations, Springer-Verlag, New York, 1996.
2. Wu J. Traveling wave fronts of reaction-diffusion systems with delay/ J.Wu, X. Zou // J. Dynamics and Differential Equations, 2001. - Vol. 13, No. 3. - P. 651-687.
3. Huang J. Traveling wave fronts in diffusive and cooperative Lotka-Volterra system with delays / J. Huang, X. Zou // J. Math. Anal. Appl, 2002. - Vol. 271. - P. 455-466.
4. Faria T. Nonmonotone travelling waves in a single species reaction-diffusion equation with delay / T. Faria, S. Trofimchuk // J. Differential Equations, 2006. - Vol. 228. - P. 357376.
5. Trofimchuk E. Slowly oscillating wave solutions of a single species reaction-diffusion equation with delay / E. Trofimchuk, V. Tkachenko, S. Trofimchuk // J. Differential Equations, 2008. - Vol. 245. - P. 2307-2332.
6. Meleshko S.V. On the complete group classification of the reaction-diffusion equation with a delay / S. V. Meleshko, S. Moyo // J. Math. Anal. Appl., 2008. - Vol. 338. - P. 448466.
7. Hernandez E. A note on partial functional differential equations with state-dependent delay / E. Hernandez, A. Prokopczyk, L. Ladeira // Nonlinear Analysis, R.W.A., 2006. - No.
4. - P. 510-519.
8. Rezounenko A.V. Stability of positive solutions of local partial differential equations with a nonlinear integral delay term / A.V. Rezounenko // Electronic Journal of Qualitative Theory of Differential Equations. Proc. 8th Coll. QTDE, 2008. - No. 17. - P. 1-7.
9. Bzheumikhova O.I. Application of Fourier method to investigation of the Dirichlet problem for partial differential equations with deviating arguments / O.I. Bzheumikhova, V.N. Lesev // International Journal of Differential Equations and Applications, 2013. - Vol.
12, No. 2. - P. 103-120.
10. Lesev V.N. Ob odnoznachnoj razreshimosti zadachi Nejmana dlja jellipticheskogo uravnenija s otklonjajushhimsja argumentom / V.N. Lesev, O.I. Bzheumihova // Jekologicheskij vestnik nauchnyh centrov ChJeS, 2012. - №3. - S. 41-46.
11. Wang L. Global exponential robust stability of reaction-diffusion interval neural networks with time-varying delays / L. Wang, Y. Gao // Physics Letters A, 2006. - Vol. 350.
- P. 342-348.
12. Lu J. G. Global exponential stability and periodicity of reaction-diffusion delayed recurrent neural networks with Dirichlet boundary conditions / J. G. Lu. // Chaos, Solitons and Fractals, 2008. - Vol. 35. - P. 116-125.
13. Tihonov A.N., Samarskij A.A. Uravnenija matematicheskoj fiziki. - M.: Izd-vo Nauka, 1977. - 735 s.
14. Lesev V.N. Primenenie metoda Fur'e k issledovaniju zadachi Dirihle dlja uravnenija s otklonjajushhimsja argumentom i operatorom Laplasa v glavnoj chasti / V.N. Lesev, O.I. Bzheumihova // Nauchnyj zhurnal KubGAU [Jelektronnyj resurs]. - Krasnodar: KubGAU, 2012. - №07(81). - S. 1-10.
15. Bzheumihova O.I. Kraevye zadachi dlja model'nyh uravnenij smeshannogo tipa vtorogo porjadka s otklonjajushhimsja argumentom / O.I. Bzheumihova, V.N. Lesev // Obozrenie prikladnoj i promyshlennoj matematiki, 2011. - T. 18, vyp. 5. - S. 744-745.
