CC BY
0
РАССАДИН А. Э.
ОБ ОДНОМ НЕЛИНЕЙНОМ ИНТЕГРО-ДИФФЕРЕНЦИАЛЬНОМ УРАВНЕНИИ Аннотация. В работе дана иллюстрация взаимного влияния эффектов нелинейности и нелокальности, а именно, найдено точное решение задачи Коши для нелинейного интегро-дифференциального уравнения, обладающего следующим свойством: изменение знака у ограниченного начального условия за конечное время приводит к неограниченному возрастанию по модулю соответствующего ему решения. Получено общее решение задачи Коши для рассматриваемого уравнения, а также продемонстрирован метод расширения таблиц преобразования Лапласа по двум переменным с помощью частных решений этого уравнения.
Ключевые слова: оригинал, изображение, свёртка, функция Бесселя первого рода, модифицированная функция Бесселя, комплексная плоскость.
Abstract. The paper illustrates the mutual influence of the effects of nonlinearity and nonlocality, namely, an exact solution of the Cauchy problem for a nonlinear integro-differential equation is found, which has the following property: a change in the sign of a bounded initial condition in finite time leads to an unbounded increase in the modulus of the corresponding solution. A general solution of the Cauchy problem for the equation under consideration is obtained, and a method for expanding the Laplace transform tables by two variables using partial solutions of this equation is demonstrated.
Keywords: pre-image, transform table, convolution, Bessel functions of the first kind, modified Bessel functions, complex plane.
Одним из способов повышения точности описания процессов различной природы является учёт при протекании этих последних нелокальности взаимодействий [1; 2]. Однако комбинация нелинейности математической модели с простейшей нелокальностью в виде пространственного сдвига уже приводит к резкому возрастанию математических трудностей (см., например, [3]), поэтому поиск точно решаемых модельных примеров, обладающих как нелинейностью, так и нелокальностью, приобретает важное значение.
Рассмотрим следующее интегро-дифференциальное уравнение:
RASSADIN A. E. ON A NONLINEAR INTEGRO-DIFFERENTIAL EQUATION
в котором неизвестная функция и(х, у,г) определена в 1-м квадранте: х > 0, у > 0 .
Снабдим уравнение (1) начальным условием:
и(х, у,0) = и0 (х, у) , х > 0, у > 0, (2)
тогда для решения задачи Коши (1)-(2) справедлива следующая теорема.
Теорема. Пусть существуют такие константы М0 > 0, \ > 0 и к0 > 0, что в 1-м квадранте функция (2) удовлетворяет неравенству:
I у)|< М0 • ехр(^х + ^0у), (3)
тогда точное решение задачи Коши (1)-(2) имеет вид:
аЬ
и(х,у,0 = Г Г ^^ • ехр(р • х + д • у) , (4)
( у ) а! ь-1 1 +1 • и0(р, д) д у) (2 • ¿)2 ' ( )
где
и0 (р, д) = Г Ги0 (х, у) • ехр(-р • х - д • у) • йхйу, (5)
0 0
а в формуле (4) прямые Яер = а и Яе д = Ь выбраны так, чтобы все особенности функции ио (р, д)/(1 +г • ио (р, д)) в комплексных плоскостях р ид соответственно оставались слева от них.
Доказательство. Нелинейный член в уравнении (1) представляет собой двойную свёртку лапласовского типа по обеим пространственным переменным, поэтому будем искать его решение с помощью двойного преобразования Лапласа — и по х, и по у .
В рамках этого формализма переход от оригинала и( х, у, г) к изображению и (р, д, г) имеет вид [4]:
и(р, д, г) = Г Ги(х, у, г) • ехр(-р • х - д • у) • йхйу, (6)
0 0
а исходное интегро-дифференциальное уравнение (1) для оригинала трансформируется в обыкновенное дифференциальное уравнение для изображения (6):
ди(р'д*) + и2(р,д,г) = 0, и(р,д,0) = ио(р,д) . (7)
дг
Начальным условием для уравнения (7) является функция (5), причём условие (3) является достаточным условием её существования в области Яе р > й0 и Яе д > к0 [4]. Задача Коши (7) имеет следующее точное решение:
и(р, д, *) = 7-777гЧ , (8)
1 + г • 7,(р, д)
а формула (4) представляет собой обращение двумерного преобразования Лапласа (6).
Однако формула (4), выражающая общее решение задачи Коши (1)-(2), неудобна для практического применения, поскольку для её использования надо проводить интегрирование в комплексном пространстве С2. Гораздо проще находить точные решения задачи Коши (1)-(2) в явном виде следующим образом: вычислить по формуле (5) изображение начального условия (2), а затем проверить по таблице двумерных преобразований Лапласа, приведённой в книге [4], соответствует ли изображение, найденное по формуле (8), какой-либо функции из этой таблицы.
Продемонстрируем применение этого приёма на примере.
Пусть начальное условие для уравнения (1) имеет вид:
щ(х,у) = А ■ J(}(2^a~xy), х>0, у >0, а0 > 0, (9)
где ^ (г) — это функция Бесселя, тогда интеграл (5) от функции (9) легко вычисляется с помощью известного разложения функции ^ (г) в степенной ряд:
ио(р,Ч) =—^ . (10)
РЧ + ао
Подставляя выражение (10) в формулу (8), получим для изображения решения:
и(р, ч, г) =-. (11)
рч+а+Аг
Функция (11) получается из функции (10) заменой а0 ^ % + Аг , поэтому вследствие единственности обратного преобразования Лапласа по переменным р и ч [4] оригинал для изображения (11) получается из функции (9) той же заменой:
и(х,у,г) = Ао ■ Jо(2л|(aо + Аг)ху). (12)
При А > о функция (12) обладает колебательным характером. На рисунке 1 приведён график функции (9) при А = 1 и а0 = 1. С течением времени характерная «длина волны» уменьшается — для этого достаточно сравнить рис. 1 с рис. 2, на котором представлен график функции (12) при тех же параметрах А = 1 и а0 = 1, но в момент времени г = 3 . Если же А < о , то при гс = а0/| А I в решении (12) происходит переход от колебательного режима к режиму с неограниченно возрастающей по модулю функцией и(х, у,г) . С формальной точки зрения этот переход выражается в замене функции Бесселя ^ (г) на модифицированную функцию Бесселя /0 (г) при г > гс. Рисунки 3 и 4 иллюстрируют этот эффект. В этом случае при г ^ гс - о характерная «длина волны» увеличивается, что
приводит к уплощению графика функции (12) до тех пор, пока при г = гс он не становится плоскостью и(х,у,гс) = — А\.
Таким образом, при переходе в решении (12) параметром А0 нулевого значения из ограниченного начального условия (9) за конечное время развивается неограниченное решение, то есть система (1) в своём поведении демонстрирует некоторое сходство с фазовым переходом второго рода [5] или с бифуркацией Андронова-Хопфа [6].
В заключение необходимо отметить, что если точное решение и( х, у, г) задачи Коши (1)-(2) из каких-либо соображений известно, то двойная свёртка лапласовского типа начального условия (2) самого с собой может быть вычислена по формуле:
х у
Иио(ху-Л) ио(^л) ■ =
ди( х, у, г)
о о
дг
(13)
г=о
а при известной функции (5) применение двойного преобразования Лапласа по переменным х и у к левой части формулы (13) может привести к пополнению таблиц преобразований Лапласа по двум переменным.
В частности, подставляя в соотношение (13) функции (9) и (12), найдём:
х у
^С^Оху). (14)
Л Ло(2л/ а0(х -%)(у -г()) 0(27■ й&ц
О 0
Наконец, с помощью выражения (10), взятого при А0 = 1, из формулы (14) без вычислений можно получить, что:
0 0 1
— ■ ^(2^аоху ) ■ ехр(-р ■ х - д ■ у) ■ йхйу = --1—-
ао (рд + ао)
СПИСОК ЛИТЕРАТУРЫ
1. Ефимов Г. В. Нелокальные взаимодействия квантованных полей. - М.: Наука, 1977. - 367 с.
2. Учайкин В. В. Метод дробных производных. - Ульяновск: Артишок, 2008. - 512 с.
3. Алёшин С. В., Глызин С. Д., Кащенко С. А. Особенности динамики уравнения Колмогорова-Петровского-Пискунова с отклонением по пространственной переменной // Моделирование и анализ информационных систем. - 2015. - Т. 22, № 5. - С. 609-628.
4. Диткин В. А., Прудников А. П. Операционное исчисление по двум переменным и его приложения. - М.: ГИФМЛ, 1958. - 179 с.
5. Ландау Л. Д., Лифшиц Е. М. Статистическая физика. Часть 1. - М.: Наука, 1976. -584 с.
6. Шильников Л. П., Шильников А. Л., Тураев Д. В., Чуа Л. Методы качественной теории в нелинейной динамике. Часть 2 / пер. с англ. В. А. Осотовой. — М.; Ижевск: НИЦ «РХД», ИКИ, 2009. - 548 с.
