Серия «Математика»
2012. Т. 5, № 1. С. 80-95
Онлайн-доступ к журналу: http://isu.ru/izvestia
УДК 517.968
О разрешимости уравнений Вольтерра I рода с кусочно-непрерывными ядрами в классе обобщенных функций *
Д. Н. Сидоров
Институт систем энергетики им. Л. А. Мелентьева СО РАН Иркутский государственный университет
Аннотация. Получены достаточные условия существования и единственности решений уравнений Вольтерра I рода с кусочно-непрерывным ядром в классе обобщенных функций с точечным носителем. Построено асимптотическое приближение параметрического семейства обобщенных решений и предложен способ уточнения его регулярной части последовательными приближениями.
Ключевые слова: уравнение Вольтерра I рода; разрывное ядро; обобщенное решение; асимптотика; последовательные приближения.
Введем в плоскости в,Ь треугольную область Б = {в,Ь;0 < в < Ь < Т} и зададим непрерывные функции в = аг(Ь),г = 1,и, имеющие непрерывные производные при Ь е (0,Т). Предполагается, что аг(0) =0, 0 < а1(Ь) < ■■■ < ап-1(Ь) < Ь при Ь е (0,Т), 0 < а1(0) < ••• < а'п-1 (0) < 1, причем кривые в = аг(Ь), г = 0,и, где а0(Ь) = 0, ап(Ь) = Ь, разбивают область Б на непересекающиеся секторы Б1 = {в,Ь : 0 < в < а\(Ь)},
Бг = {в, Ь : аг-1(Ь) < в < аг(Ь), г = 2, и}, Б = и Бг. Введем непрерывные
* Работа частично поддержана РФФИ, проекты № 12-01-00722, 11-08-00109, при поддержке Федеральной целевой программы «Научные и научно-педагогические кадры инновационной России» Минобрнауки П696 от 30.05.2010, а также при поддержке гранта Минобрнауки, номер госрегистрации НИР: 01200804682.
70-летию профессора А. С. Апарцина
посвящается
1. Введение
п
1
функции Кг(Ь, в), определенные и дифференцируемые по Ь при Ь,8 € г = 1,и.
Введем интегральный оператор
п а1(£)
J К(Ь,8)и(8)й8 = ^ У Кг(Ь, 8)х(8)й8 (1)
с кусочно-непрерывным ядром
К\(Ь, в), Ь,8 € Б\,
к (1,8) = { ......................... (2)
Kn(І, 8^ t, 8 € Бп.
Рассмотрим уравнение
£
У К (Ь, 8)и(8) й8 = / (Ь), 0 <Ь<Т <<х>, (3)
о
где функция /(Ь) имеет непрерывную производную при Ь € (0,Т),
/(0) = 0. Уравнение (3) назовем уравнением Вольтерра I рода с кусочно-
непрерывным ядром. Требуется построить в классе обобщенных функций [3] решение уравнения (3). Отметим, что в силу условия /(0) = 0 уравнение (3) не имеет классических решений. Дифференцирование уравнения (3) приводит к интегро-функциональным уравнениям и его решение в общем случае не единственно [15]. Поэтому построение решений уравнения (3) в общем случае не может быть проведено только классическими методами теории вольтерровых уравнений [1,6,15,14]. В данной работе, продолжающей цикл работ [8,9,11,17,19] по уравнениям Вольтерра, уравнение (3) рассматривается с использованием элементарных результатов интегральных, разностных уравнений, функционального анализа, распределений Соболева-Шварца и теории уравнений с функционально возмущенным аргументом нейтрального типа [7]. Статья организована следующим оборазом. В п. 2 вычисляется сингулярная компонента решения и строится интегро-функциональное уравнение для определения регулярной компоненты. В п. 3 получены достаточные условия существования и единственности обобщенного решения уравнения (3) вида и(Ь) = а5(Ь) +х(Ь), где 5{Ь) — функция Дирака, х(Ь) — регулярная непрерывная функция. Такие решения, удовлетворяющие уравнению (3) в смысле теории распределений Соболева - Шварца [3], при исследовании уравнений (3) ранее не расматривались. В п. 3 искомая регулярная часть х(Ь) единственного обобщенного уравнения (3) строится сочетанием известного в теории функциональных уравнений «метода шагов» [13] с методом последовательных приближений.
£
В п. 4 и 5 рассмотрен теоретически наиболее интересный случай, когда уравнение (3) имеет семейство обобщенных решений, зависящих от свободных параметров. Предложен метод построения асимтотиче-ских приближений параметрических решений и способ их уточнения последовательными приближениями. Важную роль при этом играет известный метод А. О. Гельфонда (см. [4, с. 338]) решения разностных уравнений.
2. Определение регулярной компоненты решения
Продолжив f(t) на отрицательную полуось нулем и продифференцировав уравнение (3), получим эквивалентное функционально-интегральное уравнение
n—1 s
F(u) d=f Kn(t,t)u(t) + ^2 ai(t)< Ki(t,ai(t))— (4)
i=1 ^
. n ai(t)
—Ki+l(t,ai(t))\u(ai(t)) + ^2 J K(1(t, s)u(s) ds = f (1\t) + f (0)5(t),
i=1 ai-i(t)
где a0 = 0, an(t) = t. Далее везде будем предполагать, что K1 (0, 0) =
0, Kn(t,t) = 0, при t Є [0,Т]. Введем функциональный оператор
n-1
Au = ^2 K—1(t,t)ai(t){Ki(t,ai(t)) — Ki+l(t,ai(t))}u(ai(t))
i=1
def n ai(t) (1)
и интегральный оператор Ku = f K—1(t,t)K() (t, s)u(s)ds.
i=1 ai-i (t)
С учетом этих обозначений уравнение (4) приводится к виду
u(t) + Au + Ku = K— 1(t, t)f (1\t) + K—1(0,0)f (0)S(t). (5)
Будем искать решение вида u(t) = aS(t)+x(t), где a—const, x(t) Є C(0,t). Легко проверить справедливость тождеств:
ai(t)
f dKl(t, s) f л dKl(t, 0)
J -щ-s(s)ds = —Щ—'
0
ai (t)
I Kpsms = 0
ai— i (t)
при г = 2,п. Действительно, первое тождество выполняется, т.к. а\(Ь) >
а-\_(1)
0, дКт’^^(я) = ЭК1дь’0д(в), I д(в)Лв = в(а^Ь)) = 1 при Ь > 0, где
о
в — функция Хевисайда. Второе тождество тоже становится очевидец ( £')
ным, если учесть, что при г = 2, п, вирр5(в) П Di = 0, / 5(в)д,в =
аг—1( ^)
в(а^Ь)) — в(а— 1(Ь)) = 0, т.к. 0 < а\(Ь) < а2(Ь) < ■■■ < ап(Ь) = Ь.
Напомним еще тождество 5(а^Ь)) = |0(0)| (см., например, [3], стр.
34). В силу отмеченных тождеств замена и = ад(Ь) + х(Ь) приводит уравнение (5) к виду К-1(0, 0)К1 (0, 0)а5(Ь) + К-1 (Ь,Ь) дК1д1’° а + х(Ь) + Ах+Кх = К- 1(Ь,Ь)/(1)(Ь)+К- 1(0, 0)/(0)5(Ь). Приравнивая в последнем равенстве коэффициенты при 5(Ь), получим а = к{(00). Регулярную часть остается определить из уравнения
х(Ь) + Ах + Кх = / (х), (6)
где / (Ь) = К- ЧМ) | / (1)(Ь) — дКт'0° кЦщ}. Отметим, что в силу операторного тождества
Кп(Ь, Ь)(1 + А + К)х = Г(х) уравнение (6) можно переписать в виде
Г(х)= /'(Ь) — дК1(Ь, 0 /(0) (7)
( ) ! () дЬ К1(0,0). ()
3. Достаточные условия существования единственного обобщенного решения уравнения (3)
Так как K1(0,0) = 0, то однородное уравнение (4) имеет только тривиальное решение среди сингулярных функций
m
Using d= ^ Ci5(%\t)
0
с точечным носителем в нуле. Поэтому существование и единственность обобщенного решения уравнения (4) вида
u(t) = Using + x(t),
где x(t) € C(o,t) эквивалентно доказательству существование единственного решения в классе С(о,т) уравнения (6). Введем функцию
n -1 Г 'I
A(t) a<i1)(t)K-1(t,t)l Ki(t,ai(t)) - Ki+i(t,ai(t))\. (*)
i=1 ^ ^
Пусть выполнены условия
(A) \A(0)\ < 1, sup \Knl(t,t)K (t, s)\< с < те.
0<s<t<T
Условие (A) выполняется, если производные а(1)(0) достаточно малы.
П
Здесь и далее ядро K(t,s) в области |J Di определено формулой (2).
1
П
Его производная по t в обычном смысле при t,s Е |J Di определяется
1
формулой
K^\t, s), t,s Е D1,
K W(t,s) = { ............
K(1(t,s), t,s Е Dn.
Теорема 1. (Достаточные условия существования и единственности обобщенного решения). Пусть выполнены условия (A), все ядра Ki(t,s) в представлении (2) непрерывны, а по t имеют и непрерывные производные, функция f (t) имеет непрерывную производную, f (0) = 0. Пусть K1(0, 0) = 0. Тогда уравнение (3) имеет единственное обобщенное решение
u(t) = Kf^0) 6(t)+ x(t), где x(t) Е C(0,t) При этом x(t) можно найти методом шагов, сочетая его с методом последовательных приближений.
Доказательство.
Так как сингулярная часть решения уже определена, то рассмотрим уравнение (6), которому удовлетворяет регулярная составляющая x(t). Зафиксировав q < 1 выберем h1 > 0 так, чтобы sup \A(t)\ =
0<t<hx
q < 1. В силу условия (A) такое h1 > 0 найдется. Положим 0 < h < min{h1, }, где постоянная с определена в условии (A). Разобьем
интервал [0,Т] на промежутки
[0, h], [h,h + eh], [h + eh,h + 2eh],.... (8)
Обозначим через xo(t) сужение искомого решения x(t) на интервал
[0, h], а через xm(t) — его сужения на интервалы
Im = [(1 + (m - 1)e)h, (1 + me)h], m = 1,2,---
Выберем e из промежутка (0,1] так, чтобы при t Е Im «возмущен-
т—1 _______ (!)
ные» аргументы ai (t) Е 1J Ik, i = 1,n- 1. Если 0 < a( (t) <
k=1
при Ь £ [0,Т), г = 1,и — 1, то указанное выше включение выполняется на промежутке [0,Т). Это включение дает возможность при построении решения х(Ь) применить известный в теории функциональнодифференциальных уравнений метод шагов (см., например, [13], стр. 199).
Для вычисления элемента Хо(Ь) £ С[о,н] построим последовательность {хП(Ь)} :
Хо№ = —Ах%-1 — Кхп-1 + 1(Ь), х0(Ь) = 1 (Ь), Ь £ [0,Ъ].
В силу выб°ра Ъ имеем оценку ||А + К||^(с(0,^)^с(0,^)) < 1
Поэтому при Ь £ [0,Ъ] существует единственное решение х0(Ь) уравнения (6). Последовательность х'П(Ь) равномерно сходится к нему. Продолжим процесс построения искомого решения при Ь > Ъ, т.е. на промежутках 1п, и = 1, 2,.... Для определенности пусть далее в (8) е =
1.
Вычислив элемент х0(Ь) £ С[о^] будем искать элемент х1(Ь) в пространстве С(Н,2Н) непрерывных вектор-функций. Найдем х1(Ь) из интегрального уравнения Вольтерра второго рода
хЦ) + ! К-1(Ь,Ь)К^,в)х(в) йв = 1 (Ь) — Ах0 — J К-1Ы)Щ,а)хо{8) йв
Н 0
последовательными приближениями. При этом хо(Ъ) = х^Ъ).
Введем непрерывную функцию
ошк {9)
являющуюся сужением искомого непрерывного решения х(Ь) на интервал [0, 2Ъ]. Тогда элемент х2(Ь) £ С(2н,зн) можно будет вычислить последовательными приближениями из интегрального уравнения Вольтерра
II рода
Ь
1
х(Ь) + J К-1М,в)х(в) йв =
2Н
2Н
= 1 (ь) — Ах1 — ! К-1(г,г)К1(г, в)х1(в) йв.
0
Продолжая этот процесс, за N шагов (Ж > Н) построим искомое решение х(Ь) £ С(0,т) уравнения (3).
Н
ь
4. Построение асимптотического приближения Х(г) регулярной части параметрических семейств обобщенного решения уравнения (3)
Рассмотрим уравнение (7), которому удовлетворяет регулярная часть обобщенного решения. Пусть выполнено условие
N ____
^) Существуют полиномы Тг = ^ в^, г = 1,п,
и+ц=1
N N ________
/N(г) = Е и^, а^1 (г) = ^ аги^,г = 1,п - 1, где 0 < ап < а2\ <
и=1 и=1
■ ■■ < ап-1}1 < 1, такие, что при г +0, в +0 справедливы оценки \Кг(г,в) -Тг(г,в)\ = о ((г + в^+1), г = Т~П, \/(г) - /N (г)\ = 0(^+1), \аг(г) - а^(г)\ = 0(^+1),г = 1,п - 1.
Разложения по степеням г, в, представленные в условии (В), далее будем называть «полиномами Тейлора» соответствующих функций. Введем функцию
-1
В(3) = Кп(0,0) + ^2(аШ1+] (Кг(0,0) - Кг+1(0,0)),
г=1
зависящую от целочисленного аргумента ], 3 € N и 0. Функцию В(3), отвечающую главной «функциональной» части уравнения (7), назовем характеристической функцией уравнения (7). Рассмотрим построение асимптотического решения уравнения (7), то есть асимптотики регулярной части искомого решения.
В отличии от п. 3, в п. 4 и в п. 5 мы не предполагаем, что однородное уравнение, отвечающее (3), имеет только тривиальное решение. Поэтому теперь решение интегро-функционального уравнения (7) может быть не единственным. Следуя работе [16] будем искать асимптотическое приближение частного решения неоднородного уравнения (7) в виде полинома
N
х(г) = ^ хз 0П г)г3. (10)
3=0
Покажем, что коэффициенты хз в общем нерегулярном случае зависят от 1п г и свободных параметров. Это согласуется с возможностью существования нетривиальных решений у однородного уравнения.
При вычислении коэффициентов Хз возможны регулярный и нерегулярный случаи.
Определение 1. Число 3 * — регулярная точка характеристической функции В(]), если В(3*) = 0 и нерегулярная точка в противном случае.
4.1. Регулярый случай: характеристическая функция В (і) = 0 при і є (0,где N достаточно велико
В этом случае коэффициенты хз будут постоянными, то есть не зависят от 1п і. Действительно, подставляя разложение (10) в уравнение (7), методом неопределенных коэффициентов с учетом условия (В), приходим к рекуррентной последовательности линейных алгебраических уравнений относительно Х3 :
В(0)хо = !(0) _ 1 (0) 1 (0) дк«' 0)
Кі(0,0) Кі(0,0) ді
(11)
і=0
В(з)хз = М3 (хо,.. .,Хз-1), 3 = 1,...,^. (12)
Правая часть М3 выражается определенным образом через решения Хо,... ,Хз-1 предыдущих уравнений и коэффициенты «полиномов Тейлора» из условия (В).
Так как в регулярном случае В(3) = 0, то коэффициенты х0,..., ХN определяются единственным образом и асимптотика (10) будет построена.
4.2. Нерегулярый случай: характеристическая функция В(]) в массиве (0, 1,..., N имеет нули
Покажем, что в нерегулярном случае коэффициенты хз будут полиномами по степеням 1п г и зависят от произвольных постоянных. Порядок полиномов и число произвольных постоянных связаны с кратностями целочисленных решений уравнения В(3) = 0.
Действительно, т.к. коэффициент хо в нерегулярном случае может зависеть от 1п г, то на основании метода неопределенных коэффициентов хо следует искать как решение разностного уравнения
-1
Кп(0,0)хо(г) + ^2 а'г(0)(Кг(0,0) - Кг+1(0, 0))хо(г + аг) = (13)
г=1
/(0) _ /(0) Кі(і, 0)
Кі(0,0) ді
і=0
где аг = 1па'(0), г = 1пг. Здесь возможны три случая:
Случай 1 (В(0) = 0).
В этом случае коэффициент хо от г не зависит и определится единственным образом из уравнения (11).
Случай 2 (В(0) = 0).
Пусть 3 = 0 — простой нуль функции В(з), то есть В(0) = 0, В'(0) = 0.
Тогда коэффициент хо(г) будем искать из разностного уравнения (13) в виде линейной функции
хо (г) = хо1г + хо2. (14)
Подставляя (14) в (13), получим для определения коэффициентов хо1, хо2 два уравнения:
В(0)хо1 = 0, (15)
В^в-»= /'(0) - /0)д
где В(0) = 0, В(1)(0) = 0. Поэтому коэффициент хо(г) линеен относительно г и зависит от произвольной постоянной. Итак, в случае
2
1
-г + с,
(16)
XQ(Z) = (f (i)(0) _f^ dKi(0,0) V
Xo(Z) У (0) Ki(o, 0) dt )
,0) dt ) B(1)(0)' где c — const.
Случай 3. Пусть j = 0 — корень уравнения B(j) =0 кратности к + 1, то есть B(0) = B'(0) = ...B(k')(0) = 0, B(k+1)(0) = 0, к > 1. Решение Xq(z) разностного уравнения (12) будем искать в виде полинома
X0(z) = X0lZk+1 + X02Zk +-------+ X0k+lZ + x0k+2. (17)
Подставляя полином (17) в уравнение (13), учитывая тождество
dk n-1
—кB(j) = ^2(а'т1+1 ak(Ki(0,0) - K+i(0,0)),
j i=1
где ai = ln ai (0) и приравнивая коэффициенты при степенях
Zk+1 zk z Z
z ,z ,...,z,z
нулю, получим рекуррентную последовательность линейных алгебраических уравнений относительно коэффициентов Xo1,Xo2,..., X0k+2 :
' B(0)x01 = 0,
B(0)X002 + B(1)(0) (^к + 1 ) X01 = 0, (18)
B(0)XOl+1 +B{1) (0)^ k+++i1_^X01 + B(l-1)(0)^ к + к1 _ ^X02 +...
• • • + B(1)(0)( к + + __ + ^Xoi = 0,l = 1,...,к,
B(0)X0k+2 + B(k+1)(0)XQ1 + B(k)(0)X02 + ... B(1)(0)X0k+1 = (19)
f (0) dK1(t, 0)
f'(0) _
K1(0,0) dt
t=o
В рассматриваемом случае В(0) = В'(0) = • • • = В(к\0) = 0, В(к+1\0) = 0. Поэтому в полиноме (17) следует положить
г 1 !(°) КЬ 0) \
хт В(к+1)(0)У' К1(0,0) Ы ).
Уравнения системы (18) превращаются в тождества В(0)х0^ = 0, 2 =
1,к + 1, т.к. В(0) = 0. Поэтому коэффициенты х02,..., х0к+2 полинома (17) остаются произвольными постоянными.
Далее применяя метод неопределенных коэффициентов с учетом тождества
_ , к(к - 1) . ..(к - (5 — ^ ] ^к-Б
в=0
I # 1пк ш = 11+1у (-1)Б к(к -1) ...(к~(*-1)) 1пк-Б г, .1 ^ 2 + 1)Б+1 ,
построим разностные уравнения для определения коэффициента Х]_(г) (г = 1п Ь) и последующих коэффициентов асимптотического приближения (10). Действительно,
def
п—1
Ь(х)
х=хо(г)+х\(г)'Ь
Кп (0, 0)х1(г) + У(аг(0))2(Кг (0, 0)- (20)
г=1
-Кг+1(0,0))х1(г + аг) + Р1(хо(г))
Ь + г(Ь), г(Ь) = о(Ь).
Здесь Р1(х0(г)) — определенный полином от г, степень которого по доказанному равна кратности решения 2 = 0 уравнения В(2) = 0. Из соотношения (20) в силу оценки г(Ь) = о(Ь) при Ь ^ 0 следует, что коэффициент х]_(г) должен удовлетворять разностному уравнению
п1
Кп(0,0)х1(г) + ^(а'(0))2(Кг(0,0) - Кг+1(0, 0))х1(г + аг)+ (21)
г=1
+Р1(хо(г)) = 0.
Если В(1) = 0, то уравнение (21) имеет решение х1(г) в виде полинома того же порядка, что и кратность решения 2 = 0 уравнения В(2) = 0. Если 2 = 1 — тоже является решением уравнения В(2) = 0, то решение х]_(г) строится в виде полинома степени к0 + к1, где к0 и к1 — кратности решений 2 = 0 и 2 = 1 уравнения В(2) = 0 соответственно. Коэффициент х1(г) будет зависеть от ко+к1 произвольных постоянных. Введем условие
(^ Пусть уравнение В(2) =0 в массиве (0,1,...,К) имеет решения 21,. ..,2 V кратностей кг, г = 1,и.
Тогда аналогичным образом можно вычислить остальные коэффициенты х2(г),...,хм(г) асимптотического приближения х(Ь) решения уравнения (7) из последовательности разностных уравнений вида
П— 1
Кп(0,0)Х](г) + ''^(а'(0))1+3 (Кі(0,0) - Кш(0, 0^х^(г + аі)+
і=1
+Гз (хо(г),..., х-1(г))) = 0, 2 = 2, N. Из изложенного вытекает
Лемма 1. Пусть выполнены условия (В), (С). Тогда существует
N
функция х(Ь) = ^ хг(1пЬ)Ьг, такая, что при Ь ^ +0 невязка решения
г=о
уравнения (7) удовлетворяет оценке
= о(1*).
При этом коэффициенты хг(1п Ь) являются полиномами от 1п Ь возрастающих степеней, не превосходящих суммы кратностей ^ кз реше-
з
ний уравнения В(2) = 0 из массива (0,1,...,г). Коэффициенты хг(1пЬ) зависят от^3=0 кз произвольных постоянных.
г
Замечание 1. Если В(2) = 0, то в сумме ^ кз соответствующие кз
з=о
полагаем равными нулю.
5. Теорема существования непрерывных параметрических семейств обобщенных решений
Так как 0 < а'г(0) < 1, аг(0) = 0, г = 1,п - 1, то для любого 0 < е < 1 найдется Т' € (0, Т] такое, что
тах \а'г(Ь)\ < е
і=1,п-1,іЄ[0,Т']
і=і,п-і,ге(о,т']
Введем условие
аі(1) ,
8Ир —і— < в.
и
ф) Пусть функция Кп(Ь,Ь) = 0 при Ь € [0, Т'] и N * выбрано настолько большим, чтобы выполнялось неравенство
где функция А(Ь) определена в п.3 формулой (*).
Лемма 2. Пусть выполнено условие (D). Пусть в классе С(о,т') функций, непрерывных при Ь € (0,Т'] и имеющих предел (возможно бесконечный) при Ь ^ +0 существует элемент х(Ь) такой, что при Ь ^ +0 невязка решения уравнения (7) удовлетворяет оценке
причем N > N*. Тогда уравнение (7) в классе С(о,т') имеет решение
где у(Ь) определяется единственным образом последовательными приближениями.
Доказательство.
Подставляя (22) в уравнение (7) получим для определения функции у(Ь) интегро-функциональное уравнение
вир є**\А(1)\<д< 1,
ІЄ(0,Т')
Г(х(1)) - Ґ(1) + К1(1,0) к^0 0) | =о(1*),
х(1) = Х(1) +1* ь(1),
(22)
аі — 1 (і)
Введем линейные операторы
Тогда уравнение (23) перепишется в компактной форме
u + (M + K )u = Y(t),
где j(t) — правая часть уравнения (23), являющаяся непрерывной функцией в силу условия леммы 2. Введем банахово пространство X непрерывных по t функций v(t) с нормой
\\v\\i = max e-lt\v(t)\, l > 0.
0<t<T'
Тогда в силу неравенств sup < е < 1 и условия (D) при У1 > 0
te(0,T']
норма линейного функционального оператора M удовлетворяет оценке
\\M\\c(X^X) < q < 1.
Кроме того, для интегрального оператора K при достаточно большом
l справедлива оценка
\\K\\c{x^x) < qi < 1 - q.
Следовательно, при достаточно большом l > 0
\\M + K\\l(X^X) < 1)
т.е. линейный оператор M + K является сжимающим в пространстве X. Поэтому последовательность {vn}, где vn = — (M + K)vn-i + Y(t), v0 = j(t), сходится.
Теорема 2. Пусть выполнены условия (B), (C) и (D), f (0) = 0,
Ki(0, 0) = 0. Тогда уравнение (3) при 0 < t < T' < T имеет решение
x(t) = Kl(o\)5(t) + x(t) +tN *v(t)l
V
зависящее от hi произвольных постоянных, где числа hi определя-
i=i
ются в условии (С). Более того, функция X строится в виде лога-рифмо-степенной суммы (10), затем v(t) вычисляется единственным образом последовательными приближениями и имеет место асимптотическая оценка '-----------
x(t) - Kl(0fi)5(t) - X(t)
= O(tN ) при t ^ +0.
Доказательство. На основании леммы 1 в силу условий теоремы возможно построение асимптотического приближения регулярной части х(Ь) искомого решения в виде логарифмо-степенного полинома
N
Xi,(\n t)ti.
i=0
При этом по построению коэффициенты Xi(\n t) будут зависеть от указанного числа произвольных постоянных. В силу леммы 2, применяя подстановку x(t) = X(t) + tN u(t), непрерывную функцию u(t) можно построить методом последовательных приближений. Теорема доказана.
Как и в теореме 1 построенное на интервале [0, T1] параметрическое семейство решений можно продолжить на весь интервал [0,T], используя метод шагов [13, с.199].
В простых случаях, решая эквивалентное уравнение (4), решение интегрального уравнения (3) можно построить в замкнутом виде.
Пример 1.
t/2 t
J x(s)ds + 2 J x(s)ds = 2 + t, t > 0.
0 t/2
Здесь эквивалентное уравнение (4) имеет вид 2x( 1) + 2x(t) = 25 (t) + 1. Искомое решение имеет вид x(t) = 25 (t) + 2/3.
Пример 2.
t/2 t
ix(s)ds 'Ix(s)ds = 1 + t't> 0
0 t/2
Здесь эквивалентное уравнение (4) имеет вид x(2,) — x(t) = 5(t) + 1. Оно имеет с-параметрическое семейство обобщенных решений x(t) = 5(t) + с — n2, с — const.
6. Заключение
Если в условиях теорем 1, 2 /(0) = 0, то уравнение (3) будет иметь непрерывное решение и придем к результатам работы [16].
Автор благодарит участников семинара (рук. А. С. Апарцин, ИСЭМ СО РАН) лаборатории неустойчивых задач вычислительной математики института систем энергетики им. Л. А. Мелентьева СО РАН за обсуждение и поддержку тематики данной статьи.
Список литературы
1. Апарцин А. С. Неклассические уравнения Вольтерра I рода: теория и численные методы / А. С. Апарцин. - Новосибирск : Наука, 1999.
2. Вайнберг М. М. Теория ветвления решений нелинейных уравнений / М. М. Вайнберг, В. А. Треногин. - М : Наука, 1969.
3. Владимиров В. С. Обобщенные функции в математической физике / В. С. Владимиров. - М : Наука, Физматлит, 1976.
4. Гельфонд А. О. Исчисление конечных разностей / А. О. Гельфонд. - М. : Физматлит, 1959.
5. Магницкий Н. А. Асимптотика решений интегрального уравнения Вольтерры первого рода / Н. А. Магницкий // ДАН СССР. - 1983. - Т. 269, № 1. — C. 29-32.
6. Маркова Е. В. О моделях развивающихся систем типа Глушкова и их приложениях в электроэнергетике / Е. В. Маркова, И. В. Сидлер, В. В. Труфанов // Автоматика и телемеханика. - 2011. - № 7. - C. 20-28.
7. Сидоров Н. А. Нелинейные операторные уравнения с функционально возмущенным аргументом нейтрального типа / Н. А. Сидоров, А. В. Труфанов // Дифференц. уравнения. - 2009. - Т. 45, № 12. - C. 1804-1808.
8. Сидоров Н. А. О малых решениях нелинейных дифференциальных уравнений в окрестности точек ветвления / Н. А. Сидоров, Д. Н. Сидоров // Изв. вузов. Математика. - 2011. - № 5. - C. 53-61.
9. Сидоров Н. А. Существование и построение обобщенных решений интегральных уравнений Вольтерры первого рода / Н. А. Сидоров, Д. Н. Сидоров // Дифференц. Уравнения. - 2006. Т. 42, № 9. - C. 1243-1242.
10. Сидоров Н. А. О решении операторно-интегральных уравнений Вольтерры в нерегулярном случае методом последовательных приближений / Н. А. Сидоров, Д. Н. Сидоров, А. В. Красник // Дифференц. уравнения. - 2010. - Т. 40, № 6. - C. 874-882.
11. Сидоров Д. H. Обобщенные решения в задаче моделирования нелинейных динамических систем полиномами Вольтерра / Д. H. Сидоров, Н. А. Сидоров // Автоматика и телемеханика. - 2011. - № 6. - C. 127-132.
12. Треногин В.А. Функциональный анализ / В. А. Треногин. - М : Наука, 1993.
- 439 с.
13. Эльсгольц Л. А. Качественные методы в математическом анализе / Л. А. Эльсгольц. - М : ГИТТЛ, 1955.
14. Яценко Ю. П. Интегральные модели систем с управляемой памятью / Ю. П. Яценко. - Киев : Наукова думка, 1991.
15. Denisov A. M. On a special Volterra integral equation of the first kind / A. M.
Denisov, A. Lorenzi // Boll. Un. Mat. Ital. B. Vol. - 1995. - Vol. 7, N 9. - P. 443-
457.
16. Sidorov D. Volterra Equations of the First kind with Discontinuous Kernels in the Theory of Evolving Systems Control / D. Sidorov // Studia Informatica Universalis. - 2011. - Vol. 9, N 3. - P. 135-146.
17. Sidorov D. N. Convex majorants method in the theory of nonlinear Volterra
equations / D. N. Sidorov, N. A. Sidorov // Banach J. Math. Anal. - 2012. -
Vol. 6, N 1. - Р. 1-10.
18. Lyapunov-Schmidt Methods in Nonlinear Analysis and Applications / N. Sidorov, B. Loginov, A. Sinitsyn, M. Falaleev. - Dordrecht-Boston-London : Kluwer Academic Publisher, 2002. - 568 p.
19. Sidorov D. On impulsive control of nonlinear dynamical systems based on the Volterra series / D. Sidorov // 10th IEEE International Conference on Environment and Electrical Engineering (EEEIC), 8-11 May 2011. Rome, Italy, 2011. - P. 1-6.
D. N. Sidorov
Solution to the Volterra equations of the 1st kind with discontinuous kernels in the class of generalized functions
Abstract. Sufficient conditions for existence and uniqueness of solutions of the equations of Volterra equation of the 1st kind with a piecewise continuous kernels in the class of generalized functions with point support are derived. An asymptotic approximation of a parametric family of generalized solutions is constructed. A method of the regular part of the solution’s improvement employs the method of successive approximations.
Keywords: Volterra integral equations; discontinuous kernel; generalized solution; successive approximations.
Сидоров Денис Николаевич, с.н.с., кандидат физико-математических наук, доцент, Институт систем энергетики СО РАН, 664033, Иркутск, ул. Лермонтова, 130, тел. (3952)429440
Институт математики, экономики и информатики, Иркутский государственный университет, 664000, Иркутск, ул. К. Маркса, 1 тел.: (3952)242210 ([email protected])
Sidorov Denis, PhD, Associate Professor, Senior Research Fellow 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])