CC BY
12
Серия «Математика»
2011. Т. 4, № 1. С. 109-117
Онлайн-доступ к журналу: http://isu.ru/izvestia
ИЗВЕСТИЯ
Иркутского
государственного
университета
УДК 517.968, 517.965
Об одном классе нелинейных уравнений I рода с однородными интегральными операторами *
Д. Н. Сидоров
Институт систем энергетики им. Л.А.Мелентьева СО РАН Иркутский государственный университет
Аннотация. Строятся обобщенные решения нелинейных уравнений первого рода, содержащих однородные интегральные операторы Вольтерра. Решения строятся в виде суммы сингулярной и регулярной компоненты. Сингулярная компонента имеет носитель в нуле и состоит из функционала Дирака и его производных. Коэффициент при старшей производной функционала Дирака удовлетворяет определенному полиному. Остальные коэффициенты сингулярной части вычисляются рекуррентным образом из линейных алгебраических уравнений. Регулярная компонента решения является непрерывной функцией и строится по вычисленной сингулярной части последовательными приближениями.
Ключевые слова: нелинейные интегральные уравнения первого рода, последовательные приближения, функционал Дирака, ряд Вольтерра, blow-up пределы.
1. Введение
Введем непрерывные функции Кп(И) где 1 = (¿0, ¿ь . . . , ¿п) £ Мп+1, п = 1, Ж, имеющие непрерывные производные по ¿о в окрестности нуля
_ _ го п
и представимые при |1| < р в виде рядов Кп(1) = П П Кп^(¿о, и). Рас-
.7 = 1 ¿=1
смотрим при 0 < и < р нелинейное интегральное уравнение Вольтерра
N го п /•
ЕЕП/ Кп%з (¿,в)ж(в)^в = f (и). (1)
п=17=1¿=10
Левая часть уравнения (1) разложена по п-однородным интегральным операторам и является частичной суммой известного ряда Вольтерра, играющего важную роль в теории динамических систем. Оператор в
Работа частично поддержана РФФИ, грант № 09-01-00377.
*
левой части уравнения (1) можно рассматривать и как специальный случай интегро-степенного оператора Ляпунова-Шмидта, использованного А. М. Ляпуновым и позднее Л. Лихтенштейном (L. Lichtenstein) при исследовании классических задач механики.
При построении модели (1) в теории динамических систем надо уметь находить ядра Knj. Для решения это задачи можно использовать метод А. С. Апарцина [2], основанный на подаче специальных тестовых сигналов x(t) и измерения соответствующих откликов f (t). В настоящее время этот подход получил дальнейшее развитие в теоретических и прикладных исследованиях. Вопросы построения непрерывных решений уравнения (1) изучались в частных случаях. Исследование разрешимости уравнения (1) в пространстве распределений [4] впервые строго рассмотрено в [7], но лишь для класса “степенных” уравнений
N
В этой работе, продолжающей исследование [7], [8], решение уравнения общего вида строится в классе обобщенных функций 0(-р, р)'. Показано, что число решений зависит от количества корней определенного алгебраического уравнения. Предложен метод построения искомых решений. Отмечается, что в случае /(0) = 0 при определенных условиях существует и последовательными приближениями строится непрерывное решение уравнения (1).
Статья организована следующим образом. В п. 2 рассмотрены решения с сингулярностью нулевого порядка. В п. 3 строятся решения с сингулярностью т-го порядка, при т > 1. В п. 4 содержит заключение.
2. Решения с сингулярностью нулевого порядка
Будем искать решение уравнения (1) в виде
х(1) = с5(1) + ь(Ь), (2)
где 5(г) - функционал Дирака, у(г) - непрерывная функция. Решение (2) рассматривается как элемент пространства ^|_р р) [6], определенного на множестве бесконечно-дифференцируемых финитных функций ф(Ь) из р) с носителями на (-р,р). Функция /(¿) в правой части уравнения (1) считается продолженной нулем в левую полуплоскость точки 1 = 0. В пространстве 0'_р р) имеет место тождество
£ £
J Кц7(г,в)(с5(в)+ у(в))йв = Кп7(¿, 0)с + У Кп7(¿,8)у(8)ёв.
00
Поэтому уравнение (1) подстановкой (2) сводится к уравнению относительно непрерывной функции у(г)
N го п £,
Е Е П(К™?(и’ 0)с + Кпг](г, .е)у(,е)(1е) = /(г) (3)
п=17=1г=1 0
с неизвестной постоянной . Полагая в (3) 1 = 0 получим алгебраическое уравнение
N го п
Ь(с) = Е сп ЕП Кпу (0, 0) = /(0). (4)
п=1 7=1 г=1
Уравнение (4) назовем характеристическим.
Введем условие Л) с* - простой вещественный корень уравнения (4).
Тогда, полагая в (3) с = с* и дифференцируя полученное равенство по ¿, получим уравнение
N го п , £ дК (1 е) \
ЕЕЕ (К>пи (¿, 0)с* + КпЦ (и,ф(и)+ ----п^, , X
1 • 1 / IV дг
п=17=11=1 4 о
х П (кп^(г,0)с* + К,пц(¿,8)у(в)сЛ = /'(г). (5)
г=1 , г=Л 0 '
Полагая в (5) 1 = 0 получим алгебраическое линейное уравнение относительно у(0) :
N го п / \ п
ЕЕЕ кп г, (0,0)с* + Кпц (0, 0)у(0)) П Ктз (0,0)с* = /' (0). (6)
п=1 7=1 г=1 ' ' г=1 , г=г
Коэффициент при у(0) в уравнении (6) имеет вид
N го п
ЕЕЕ Кпг7 (0, 0) П Кп%з (0, 0)с .
п=1 7=1 1=1 г=1 , г=г
Он преобразуется к виду
N го п
¿'(с*) = Е пс*п-1 ЕП Кпг^ (0, 0)
п=1 7=1 г=1
и поэтому в силу условия А) не равен нулю. Следовательно, у(0) из уравнения (6) определяется единственным образом.
Уравнение (5) преобразуем к виду
а(г)у = г (у,г) + ь(г), (7)
где
N го п п
а(Г) = Ес*”-1 ЕЕ П Кп%3 (г, 0)Кп^ (г,г),
п=1 7=1 г=1 г=1,г=г
Г (у, 0) = 0. Таким образом, вопрос свелся к построению непрерывной функции у(г) ^ 0 при г ^ 0, удовлетворяющей уравнению (7).
Т.к. а(0) = Ь'(с*) = 0, у(0) = , то существуют положительные
числа р> 0,г > 0, такие, что уравнение (7) в пространстве С[_р,р] при у € 5(у(0),г) будет удовлетворять всем условиям принципа сжимающих отображений ([9], с.381). При этом приближения
уп(г) = щ(Г(уп_1,1) + ь(t)), уо =
при |1| < р равномерно сходятся к регулярной части у(1) обобщенного решения (2).
Если выполнено условие
го
В) /(0) = 0, ^ К117(0,0) = 0, то с* будет простым решением урав-7=1
нения (4) и уравнение (1) имеет классическое решение. Конечно при выполнении условия В) уравнение (4) кроме нулевого может иметь и другие простые вещественные решения. Поэтому уравнение (1) в общем случае наряду с классическим может иметь несколько обобщенных решений. Если /(0) = 0, то непрерывных решений у уравнения (1) нет. Из изложенного вытекает
Теорема 1. Пусть выполнено условие А). Тогда уравнение (1) имеет обобщенное решение вида (2), где с* = 0 - простой корень уравнения (4), а непрерывная функция у(1) определяется единственным образом последовательными приближениями. Если выполнено условие В) , то уравнение (1) в классе непрерывных функций имеет единственное решение и его можно построить последовательными приближениями в окрестности точки 1 = 0.
Пример 1.
£ 2 £
(¡фуь) + / ФУ1, +1 = м> 0
оо
В этом примере равенство с2 + с = 0 является характеристическим уравнением. Оно имеет два простых решений с1 = 0, с2 = —1. Поэтому
на основании теоремы 1 данное интегральное уравнение имеет в пространстве D ровно два решения x\(t) = —S(t) — yj—t, x2(t) = yj—t•
Точка to = j является точкой blow-up предела этих решений, т.е. при приближении к этой точке решение принимает сколь угодно большие значения по абсолютной величине.
На основании результатов работы [4] и в общем случае решение v(t) уравнения (7) за конечное время t может уйти в то (или разветвиться). Оценку снизу интервала, на котором регулярная часть решения (2) существует и единственна можно получить, используя схему Л.В.Канторовича (см. [5], гл. 12, п.2) и работу [12] при построении и исследовании мажорантных интегральных и алгебраических уравнений для уравнения (7).
Таким образом, в конкретных примерах, сочетая результаты этого параграфа с методами работы [12], можно вычислять нижнюю границу t*, при прохождении которой у решений интегрального уравнения (1), могут быть точки ветвления или явления взрыва (blow-up пределы в зарубежной литературе). Такие явления в последнее время привлекали внимание многих математиков (см. [9], c.15, c.247). В случае нелинейных интегральных уравнений Вольтерра это явление остается неизученным, проявляет себя при численных расчетах. Некоторые результаты в этом направлении см. в [1].
3. Решения уравнения с сингулярностью порядка т, т > 1.
В этом параграфе рассматривается нелинейное уравнение в свертках вида
N го п
ЕЕП Кпи *х =/ ^ 0 (8)
п=17=1¿=о
£
где КпЦ * х = I Кпу (1 — 8)х(8)с18. о
Пусть выполнено условие
С) функции Кпу (1), /(1)- т-раз дифференцируемы, причем первые т—
1 производных функций Кпу (1) в точке 1 = 0 равны нулю, причем в
окрестности нуля Кпу (1) ~ при 1 ^ 0.
Будем строить решение уравнения (1) вида
х(1) = ст6(т)(1) + ст_1б(т_1)(1) + ••• + со6(1) + у(1), (9)
где сг - искомые постоянные, у(Ь) € С[о,р). Подстановка функции (9) в уравнение (8) преобразует его в пространство Д [6] к виду
N го п / \
Е Е П (^Й^т + ^ПГ^СФт-! + ... КпЦ (Ь)с0 + К.пЦ * И = I (0).
п=17=1г=1 '
(10)
Полагая в (10) Ь = 0 получим алгебраическое уравнение для определения ст :
N го п
¿(с) = Е ст Е П т!апц = / (ь). (и)
п=1 7=1 г=1
Пусть
о) ст - простой вещественный корень уравнения (11).
Положим в (10) ст = с*т. Далее, последовательно дифференцируя уравнение (10) по Ь и полагая затем Ь = 0, придем к рекуррентной последовательности линейных алгебраических уравнений вида
Ь'(с*т)ст-г = 1 (г) (0) + bi(c*m,c*m-1, . . . ,c*m-(г-1)), * = 1 2,...,т.
Проводя (т + 1)-ое дифференцирование по Ь, с уже вычисленными с*, г = т,т — 1,.., 0, получим эквивалентное нелинейное интегральное уравнение Вольтерра второго рода вида (7), где теперь
N го п п
а(£) = Е стп-1 Е Е кпт)(£) П т]апч.
п=1 7=1г=1 г=1г=г
Отметим, что
N го п
а(0) = Е Пспп-1 Е П т! ашз ■
п=1 7=1г=1
Поэтому в силу условия о) а(Ь) = 0 в некоторой окрестности точки Ь = 0. Следовательно, соответствующее интегральное уравнение относительно функции у(Ь) в пространстве С(0,р в шаре Б(у(0),т) € С(0,р при достаточно малых р > 0,г > 0 удовлетворит условиям принципа сжимающих отображений [9], с. 381.
Из изложенного вытекает
Теорема 2. Пусть выполнены условия С) и о). Тогда уравнение (1) в классе Д(_р имеет решение (9), где ст - простой веществен-
ный корень уравнения (11), ст-1,... ,с0 вычисляются однозначно из последовательности линейных алгебраических уравнений. Непрерывная функция у(Ь) строится однозначно последовательными приближениями.
4. Заключение
В п. 2 рассмотрены обобщенные решения с сингулярностью порядка т > 1 нелинейных уравнений с ядрами, зависящих от разностей. Покажем, что обобщенные решения можно строить и в случае ядер общего вида. Действительно, вместо условия С) введем условие С’) Пусть в уравнении (1) функции Кпу-(Ь,в) имеют вид
т
(Ь, 8) = Е к,т_" г-^ + <9((|г| + И)т+>)
V=0
и алгебраическое уравнение
N го п
ад = Е с” Е П(-1)тК» = / (0) (12)
п=1 7=1г=1
имеет простое вещественное решение с .
Тогда решение уравнения (1) в классе Дт(_р р) можно искать в виде
х(Ь) = с *5(т')(Ь) + ст_15(т-1') (Ь) +-+ с05(Ь) + у(Ь). (13)
Подставляя (13) в (1) и полагая Ь = 0 с учетом условия С’) придем к уравнению (12) относительно с *. Аналогично доказательству Теоремы 2 для определения коэффициентов ст_ 1,... ,с0 получим линейные алгебраические уравнения. Для вычисления непрерывной функции ь(Ь) после необходимого числа дифференцирований придем к нелинейному интегральному уравнению вида (7), в котором теперь а(Ь) ~ аЬр, а = 0, 0 < р < т. Напомним, что в п.1 и в п. 2 в соответствующем уравнении (7) выполнялось условие а(0) = 0, то есть оно содержало интегральный оператор Вольтерра второго рода. При выполнении условия С’) может оказаться, что р > 1. Тогда а(0) =0 и соответствующее интегральное уравнение вида (7) будет интегральным уравнением Вольтерра третьего рода. Решение интегральных уравнений Вольтерра третьего рода в окрестности особой точки Ь = 0 можно строить, сочетая метод неопределенных коэффициентов для определения начального приближения с дальнейшим его уточнением последовательными приближениями (см. детали, например, в доказательстве теоремы 1 в [4]).
Список литературы
1. Апарцин А. С. Об эквивалентных нормах в теории полиномиальных интегральных уравнений Вольтерра I рода / А. С.Апарцин // Изв. Иркут. гос. Ун-та. Сер.: Математика. - 2010. - Т. 3, № 1. - С. 19-29.
2. Апарцин А. С. Неклассические уравнения Вольтерра первого рода. Теория и классические методы / А. С. Апарцин. - Новосибирск : Наука, 1999.
3. Вайнберг М. М. Теория ветвления решений нелинейных уравнений / М. М. Вайнберг, В. А. Треногин. - М. : Наука, 1969.
4. Владимиров В. С. Уравнения математической физики / В. С. Владимиров. -М : Наука, 1988. - 512 с.
5. Канторович Л. В. Функциональный анализ в полуупорядоченных пространствах / Л. В. Канторович, Б. З. Вулих, А. Г. Пинскер. - М. ; Л., 1950.
6. Сидоров Н. А. Существование и построение обобщенных решений нелинейных интегральных уравнений Вольтерры первого рода / Н. А. Сидоров, Д. Н. Сидоров // Дифференц. уравнения. - 2006. - Т. 42, № 9. - С. 1243-1247.
7. Сидоров Д. Н. Обобщенные решения полиномиальных интегральных уравнений первого рода в одной модели нелинейной динамики / Д. Н. Сидоров, Н. А. Сидоров // Автоматика и телемеханика. - 2011 (в печати).
8. Сидоров Д. Н. Моделирование нелинейных нестационарных динамических систем рядами Вольтерра: идентификация и приложения / Д. Н. Сидоров // Сиб. журн. индустр. математики. - 2000. - Т. 3, № 1. - С. 182-194.
9. Треногин В. А. Обыкновенные дифференциальные уравнения / В. А. Трено-гин. - М. : Физматлит, 2009. — 311 с.
10. Belbas S. A. Numerical solution of multiple nonlinear Volterra integral équations / S. A. Belbas, Y. Bulka // Applied Mathematics and Computation. - 2011. - Vol. 217, Issue 9. - P. 4791-4804.
11. Lichtenstein L. Vorlesungen über einige Klassen nichtlinearer Integralgleichungen und Integro-differentialgleichungen nebst Anwendungen / L. Lichtenstein. - Berlin : Julius Springer, 1931.
12. Sidorov D. N. Convex Majorants Method in the Theory of Nonlinear Volterra Equations / D. N.Sidorov, N. A.Sidorov. - arXiv:1101.3963v2 [math.DS].
D. N. Sidorov, N. A. Sidorov
Method of monotone majorants of the theory of nonlinear Volterra equations
Abstract. The authors have constructed the main solutions of nonlinear operatorintegral equations of Volterra in sense of Kantorovich. Convergence of the successive approximations is established through studies of majorants of integral and algebraic equations. Estimates are given for the solutions and for the intervals on which right margin the solution has the blow-up limit or start branching.
Keywords: majorants, Volterra operator-integral equations, blow-up limit, successive approximations.
Сидоров Денис Николаевич, кандидат физико-математических наук, доцент, Институт систем энергетики СО РАН, 664033, Иркутск, ул. Лермнотова, 130, тел. (3952)429440; Институт математики, экономики и информатики, Иркутский государственный университет, 664000, Иркутск, ул. К. Маркса, 1 тел.: (3952)242210 ([email protected])
Сидоров Николай Александрович, доктор физико-математических наук, профессор, Институт математики, экономики и информатики,
Иркутский государственный университет, 664000, Иркутск, ул. К. Маркса, 1 тел.: (3952)242210 ([email protected])
Sidorov Denis, Energy Systems Institute SB RAS, 130 Lermontov Str., Irkutsk, 664003 Irkutsk State University, 1, K. Marks St., Irkutsk, 664003 professor, Phone: (3952)428440 ([email protected])
Sidorov Nikolay, Irkutsk State University, 1, K. Marks St., Irkutsk, 664003 professor, Phone: (3952)242210 ([email protected])
