Динамические системы, том 1(29), №2 (2011), 227-241 УДК 517.983
Слабонелинейные краевые задачи для операторных уравнений в критическом случае
В. Ф. Журавлев
Житомирский национальный агроэкологический университет, Житомир 10008. E-mail: vfz2008iQukr.net
Аннотация. Рассмотрены слабонелинейные краевые задачи для операторных уравнений с нете-ровым оператором в линейной части краевой задачи в критическом случае. Получены необходимые и достаточные условия существования единственного решения, построена сходящаяся итерационная процедура для его построения.
Ключевые слова: слабонелинейная краевая задача, нетеров оператор, критический случай первого порядка.
1. Вступление
Решение слабонелинейных краевых задач существенным образом зависит от возможности построения решений соответствующей линейной (порождающей) краевой задачи и предполагает использование информации о линейном операторе исходного операторного уравнения. Оно основано на построении обобщенного оператора Грина порождающей линейной полуоднородной краевой задачи. Такие краевые задачи рассматривались для систем обыкновенных дифференциальных уравнений в периодическом случае [5], для систем обыкновенных дифференциальных уравнений и систем функционально-дифференциальных уравнений в общем нетеровом случае [1], [2], для импульсных систем обыкновенных дифференциальных уравнений [12] и импульсных систем с запаздывающим аргументом [3].
Одной из отличительных особенностей этих задач является то. что исходные дифференциальные системы являются всюду разрешимыми [9]. Для не всюду разрешимых операторных уравнений подобные краевые задачи еще мало изучены. К таким задачам относятся, например, слабонелинейные краевые задачи для интегро-дифференциальных, функционально-дифференциальных уравнений и др., у которых исходное линейное операторное уравнение является нетеровым, то есть, не всюду разрешимым.
В данной статье ставятся следующие задачи. С использованием конструкций обобщенно-обратного оператора в банаховом пространстве найти критерии разрешимости и формулы для представления решений линейных краевых задач для нетеровых операторных уравнений, построить обобщенный оператор Грина для
© В. Ф. ЖУРАВЛЕВ
этих задач. Получить необходимые и достаточные условия разрешимости слабонелинейных краевых задач с нетеровым оператором в линейной части исходного операторного уравнения. Применяя методы типа Ляпунова-Шмидта [1], [2], [5], [7], [10], найти условия перехода от исходной краевой задачи к операторной системе, для решения которой применимы итерационные алгоритмы. Используя принцип неподвижной точки, построить сходящиеся итерационные алгоритмы для отыскания решений таких краевых задач.
2. Постановка задачи
Обозначим через Bi банахово пространство вектор-функций z: I ^ R", определенных на конечном промежутке I; B2 — банахово пространство вектор-функций f: I ^ R"1; I£0 = [0,е0] — промежуток значений малого параметра е; Cn[e] — пространство непрерывных по е функций c: I£0 ^ R".
Рассмотрим нелинейную краевую задачу с малым неотрицательным параметром е, е Е I£0
Lz(,e)(t) = f (t) + eZ (z(t,e),t,e), (2.1)
£z(,e) = a + eJ (z(,e),e), (2.2)
(al) L : B1 ^ B2 — нетеров оператор, (/ = dim N(L) < ж, v = dim N(L*) < ж, / = v);
(a2) Z : B1 x I x I£0 ^ B2 x I£0 — нелинейный no z ограниченный оператор, который в окрестности порождающего решения \\z — z0\\ < q имеет производную Фреше по z и непрерывен по е q — достаточно малая
КОНСТ&НТсЦ
(аЗ) Z(0,t, 0) = 0, Z'z(0, t, 0) = 0;
(а4) f Е B2; 1 ' )
(а5) £ : B1 ^ Rm — линейный ограниченный вектор-функционал; (аб) J : B1 x I£0 ^ Rm x I£0 — нелинейный no z ограниченный вектор-функционал, который в окрестности порождающего решения \\z — z01| < q имеет производную Фреше по z и непрерывный по е; (а7) J(0, 0) = 0, JZ(0, 0) = 0; (а8) a Е Rm.
Наряду с задачей (2.1), (2.2) рассмотрим линейную краевую задачу
Lzo(t) = f (t), (2.4)
£zo() = a, (2.5)
е=0
По аналогии с подобными задачами для систем обыкновенных дифференциальных уравнений (см. [1], [2], [5]), краевую задачу (2.4), (2.5), которая получается из
задачи (2.1), (2.2) при е = 0, будем называть порождающей краевой задачей для задачи (2.1), (2.2). Решения задачи (2.4), (2.5) будем называть порождающими.
Рассмотрим задачу об условиях существования и построении решений г(Ь,е) краевой задачи (2.1), (2.2), принадлежащих пространству В1 по Ь, пространству Сп[е] по е и обращающихся при е = 0 в порождающее решение краевой задачи (2.4), (2.5). В дальнейшем это пространство будем обозначать В1Сга[е].
3. Линейные краевые задачи
Сначала рассмотрим задачу об условиях разрешимости и представлении решений порождающей линейной краевой задачи (2.4), (2.5) с нетеровым оператором.
Обозначим через (Ь) = ¡1 < то и (Ь*) = V < то — размерности
нуль-пространств операторов Ь и ему сопряженного Ь*, соответственно. Пусть {¡г}= С N(Ь), ¡.г = со1(!(1), I (2\ ..., ¡( ) — базис подпространства N(Ь), а {рДОК^ С ^(Ь), р.в() = со1(р.^ (•), р."1 (•),..., р^О) — базис подпространства ^(Ь) С В*. Для элементов {¡ги функционалов {р«(^)}^=1 существуют сопряженно биортогональные система функционалов (•)}^=1 С ВЦ, (•) =
соЦт^ (•), (•),..., Т^О) и полная система эле ментов {фк }'к=1 С Уь С В2, фк = со1(фк}\ ф^фк)• Подпространство Уь С В2 изоморфно нуль-пространству
N(Ь*) сопряженного оператора Ь*.
Функционалы (•)}1^=1 и {р^ОК^ определенные на подпространствах N(Ь) С В1 и Уь С В2 (по теореме Хана-Банаха), могут быть продолжены с сохра-
В1 В2
Обозначим
^ = (¡1, ¡2, ..., М, Г (•) = Ы^),ъ(^),...,ъ()т, Ф(^) = (р1(^),Р2(^),...,Ри ()Т, & =(ф1, Ф2, ..., Фи)
матрицы размерностей (и х ц), (1 х и), (V х п1) и (п1 х V), соответственно. Отметим, что Г(X) = Ец, Ф(&) = Еи — единичные матрицы.
По аналогии с [6], операторы проектирования Рм(ь) ■ В1 ^ N(Ь) и Руь ■ В2 ^ Уь построим по формулам:
Гмщ(г) = ХГ(г), Уг е Вь (3.1)
Туь (у) = ФФ(у), У у е В2. (3.2)
Сначала рассмотрим условия существования и представление решений неоднородной задачи Коши
Ьг(Ь) = I (Ь), (3.3)
г(Ьо) = го,Ьо е1. (3.4)
Известно [2], [12], что операторное уравнение (2.4) с нетеровым оператором, вследствие его нормальной разрешимости, имеет решение для тех и только тех f (t) Е B2, которые удовлетворяют условию
(Pyl f )(t) = Ф (t)(Фf )() = 0. (3.5)
Условие (3.5), в силу линеинои независимости столбцов матрицы Ф (t), эквивалентно условию
Ф )(0 = 0.
При выполнении условия (3.5) общее решение уравнения (2.4) представимо в виде:
zo(t) = (VN (L)z)(t) + (L- f )(t), (3.6)
где z = z(t) — произвольный элемент пространства B1, L- — ограниченный
L
ленпе (3.1) проектора PN(L), для каждого z(t) Е N(L) получим
z (t) = (Pn (L)z)(t) = X (t)(rz)() = X (t)c, (3.7)
где c = co\(y1z(^),y2z(^), ... ,7^2(0) — произвольный вектор-столбец констант, c Е R^.
С учетом (3.7) общее решение нетерового уравнения (2.4) примет вид:
zo(t) = X (t)c + (L- f )(t), (3.8)
где X(t) — (/ x /)-мерная фундаментальная матрица уравнения (2.4).
Обозначим через X(t0) : RM ^ R"— (n x /)-мерную постоянную матрицу, t0 Е I, PN(X(t0)) : RM ^ N(X(t0)) — (/ x /)-мерную матрицу-ортопроектор, PN(x*(t0)) : R" ^ N(X*(t0)) — (n x п)-мерную матрицу-ортопроектор, X +(t0) — (/xn)-MepHyro матрицу, псевдообратную к матрице X(t0) [1], [2], [11].
(t) N(L)
в виде
z(t) = X (t)c.
Это соотношение справедливо для любого t Е I, в том числе и для t = t0. Следовательно, для любого z0 = z0(t0) существует элемент c Е RM такой, что выполняется равенство
z0 = z0(t0) = X (t0 )c. Существование для любого z0 = z0 (t0) элемента c Е RM означает, что уравнение
X (t0)c = z0
разрешимо при любой правой части z0 Е R". Это значит, что rankX(t0) = п, нуль-пространство N(X*(t0)) = {0} и, как следствие, ортопроектор на нуль-пространство N(X*(t0)) С R" равен нулю, PN(X*(t0)) = 0, Wt0 Е I.
При t = t0, из (3.8) получим алгебраическое уравнение
X (t0)c = z0 — (L-f )(t0),
которое, с учетом PN(X*(to)) = 0, всюду разрешимо и имеет общее решение
c = Pn (X (to))cM + X+(t0)[z0 — (L-f )(t0)], (3.9)
где cц Е RM — произвольный вектор-столбец констант.
c
дачи Коши (3.3), (3.4)
z(t) = X (t) {Pn (x (to))c, + X+(t0)[z0 — (L-f )(*,)]} + (L-f )(t) =
= X0(t)c, + X(t)X + (t0)z0 + (G0f)(t), [ ' }
где X0(t) = X(t)PN(X(to)), (G0f)(t) = (L-f)(t) — X(t)X+(t0)L-f)(t0) — линейный ограниченный оператор, который будем называть оператором Грина полуоднородной (z0 = 0) задачи Коши (3.3), (3.4).
Таким образом, ДЛ Я 3 ад с1Ч и Коши (3.3), (3.4) справедливо следующее утверждение.
Теорема 1. Пусть L : B1 ^ B2 — нетеров оператор. Задача Коши (3.3), (3.4) разрешима для любых z0 = z(t0), t0 Е I и тех и только тех f (t) Е B2; которые удовлетворяют условию (3.5), и ее общее решение представимо в виде
z(t) = X0(tc + X (t)X+(t0)z0 + (G0f )(t).
Замечание 1. Если операторное уравнение (3.3) всюду разрешимо, то PyL = 0
L
ограниченный правый обратный оператор L-1 [6], а обобщенный оператор Грина будет иметь вид:
(G0f )(t) = (L-f )(t) — X (t)X+(t0)L-1)(t0).
Далее рассмотрим условия разрешимости и представление решений линейной порождающей краевой задачи (2.4), (2.5).
Обозначим через Q = £X(•) (т x /)-мерную постоянную матрицу, PN(q) : RM ^ N(Q) — (/ x i^^^^^^^^ro матрицу-ортопроектор; PN(q*) : Rm ^ N(Q*) — (m x m)-мерную матрицу-ортопроектор; PNp(Q) — (/ x р)-мерную матрицу, составленную из полной системы р = / — rank Q линейно независимых столбцов матрицы-ортоироектора PN(q), PNd(q*) — (d x т)-мерную матрицу, составленную из полной системы d = т — rank Q линейно независимых строк матрицы-ортопроектора PN(Q*)-, Q+ — ^^^^^^^^^^таую к матрице Q (/ x т)-мерную матрицу, которую можно построить по формулам Q+ = (Q*Q + Pn(q))-1Q* = Q*(QQ* + PN(Q*))-1 [11].
Для порождающей краевой задачи (2.4), (2.5) справедливо следующее утверждение.
Теорема 2. Пусть L : B1 ^ B2 — нетеров оператор. Тогда, если rank Q < minim, ß), то соответствующая (2.4), (2.5) однородная (f (t) = 0,а = 0) краевая задача имеет р и только р линейно независимых решений
z(t) = Xp(t)cp, cp е Rp,
где Xp(t) = X(t)PNp(Q) — in x р)-мерная фундаментальная матрица однородной краевой задачи.
Неоднородная краевая задача (2.4), (2.5) разрешима для тех и только тех f (t) е B2 и а е Rm, которые удовлетворяют v + d линейно независимым условиям
PNdm{a - £L-f )(•)} = 0 р
z(t, cp) = Xp(t)cp + (Gf )(t) + X(t)Q+a, (3.12)
где G : B2 ^ ker£ С Bi,
(Gf)(t) = (L-f)(t) - X(t)Q+£(L-f)(•) (3.13)
(а = 0)
Доказательство. Подставив (3.8) в (2.5), получим алгебраическую систему
Qc + £(L-f )(•) = а, которая разрешима тогда и только тогда, когда выполнено условие [2, с. 92]
Рщ(Q-){a - £(L-f )(•)} = 0, р
c = PNp(Q)cp + Q+{a - £(L-f )(•)}. (3.14)
Подставляя (3.14) в (3.8), получим общее решение краевой задачи (2.4), (2.5)
zo(t) = z(t, cp) = X(t){PNp(Q)cp + Q+[a - £(L-f )(•)]} + (L-f )(t) =
= Xp(t)cp + (L- f )(t) - X(t)Q+£(L- f )(•) + X(t)Q+a = = Xp(t)cp + (Gf )(t)+ X (t)Q+a.
□
В качестве следствий из доказанной теоремы, рассмотрим два «крайних» случая, когда rank Q = ß и rank Q = m.
Следствие 1. Если rank Q = то / < т. В этом случае однородная краевая задача (2.4), (2.5) не имеет решений, кроме тривиального.
Неоднородная краевая задача (2.4), (2.5) с нетеровым оператором L : B1 ^ B2 разрешима для тех и только тех f (t) Е B2 и а Е Rm; которые удовлетворяют v + d, (d = т — /) линейно-независим,ым условиям
Ф )(-) = 0,
PNd(Q*){a — IL-f )(•)} = 0
и при этом имеет единственное реиление
z(t) = (Gf )(t)+ X (t)Q+a.
Доказательство. Если rank Q = то PN(q) = 0 и Xp(t) = 0. В этом случае псевдообратная матрица Q+ находится по формуле Q+ = (Q*Q)-1Q*. □
Следствие 2. Если rank Q = т, то т < /.В этом случае однородная краевая
рр
z(t) = Xp(t)cp, cp Е Rp, р = i — т.
Неоднородная краевая задача (2.4), (2.5) с нетеровым оператором L : B1 ^ B2 разрешима для тех и только тех f (t) Е которые удовлетворяют v линейно независимым условиям
W )(^) = 0
р
z(t, cp) = Xp(t)cp + (Gf )(t) + X(t)Q+a.
Доказательство. Если rank Q = т, то PN(q*) = 0 и второе условие в (3.11) будет всегда выполнено. В этом случае псевдообратная матрица Q+ находится по формуле Q+ = Q*(QQ*)-1. □
Замечание 2. Если операторное уравнение (2.4) является всюду разрешимым, то краевая задача (2.4), (2.5) будет нетеровой. Такие задачи подробно рассмотрены в случаях, когда (Lz)(t) = z(t) — A(t)z(t) — обыкновенный дифференциальный оператор [1], [12], дифференциальный оператор с сосредоточенным запаздыванием (Lz)(t) = z(t) — A(t)(Sh)z(t) [2], [3]. В этих случаях правый обратиый оператор L-1 имеет интегральное представление [8]
(L-1f )(t)= / "к (t,s)g(s)ds,
J a
где K(t,s) — матрица Коши.
Обобщенный оператор Грина имеет вид:
г-Ъ г-Ь
(Gf )(t)= к (t,s)f (s)ds — X (t)Q+i K (;s)f (s)ds.
Замечание 3. Если функционал I удовлетворяет условию
г-Ъ г-Ь
е к (,в)!(в^в = ек (в) ¿в,
¿а ¿а
то обобщенный оператор Грина имеет представление
С! )(ь) = Г С(г,в)! (в) ¿в,
а
ядро которого С(Ь, в) = К(Ь, в) — X(Ь^+еК(•, в) называется обобщенной матрицей Грина.
4. Слабонелинейные краевые задачи. Критический случай первого порядка
Рассмотрим краевую задачу (2.1), (2.2) в случае, когда порождающая краевая задача (2.4), (2.5) неоднозначно разрешима. Такой случай для аналогичных задач для всюду разрешимых обыкновенных дифференциальных систем назван [1] критическим случаем первого порядка.
Будем искать условие существования и алгоритм построения решений г(Ь,е) Е Б!Сп[е\ краевой задачи (2.1), (2.2), обращающихся при е = 0 в одно из порождающих решений уравнения краевой задачи (2.4), (2.5). Из теоремы 2 следует, что порождающая краевая задача (2.4), (2.5) имеет решение тогда и только тогда, когда !(Ь) Е Б2 и а Е И™ удовлетворяют условиям (3.11), при выполнении которых она имеет семейство решений (3.12) г(Ь,ср) = г0(Ь,ср).
4.1. Необходимое условие существования решений
Предположим, что !(Ь) Е Б2 и а Е И™ таковы, что условия (3.11) выполнены. Покажем, что имеет место следующая теорема.
Теорема 3. Пусть краевая задача (2.1), (2.2) удовлетворяет условиям (2.3) и имеет реиление г(Ь,е), непрерывное по е Е [0,е0], обращающееся при е = 0 в некоторое порождающее реиление г0(Ь,ср) вида (3.12), полученное при ср = с0. Тогда элемент с0 Е Ир удовлетворяет системе уравнений
ФЕ(го(, со) + х(,е), •,е) = 0,
{З(%)(•, со) + х(,е),е) — еь-Е(го(^со) + х(^,е), •,е)} = 0.
Доказательство. Пусть выполнены условия теоремы. Тогда при всех 0 < е < е0 справедливы то^кдества
Ьг(^е)(Ь) = ! (Ь)+ еЕ (г(Ь,е),Ь,е), (ег)(^) = а + еЗ (г(^,е),е).
Используя теорему 1 и учитывая справедливость условий (3.11), для этой краевой задачи получим, что при всех е Е Х£0, е = 0 нелинейные оператор Z(z(t,e),t,e) и функционал J(z(,e),e) удовлетворяют условиям
ФZ(z(;е), •,£) = 0,
PNd(Q*){J(z(,e),e) — £L-Z(z(,e), ;е)} = 0. ^
Тогда, поскольку нелинейные оператор Z(z, t, е) и функционал J(z(,e),e) удовле-
(е = 0)
zo(t, c0), то, переходя к пределу при е ^ 0, z(t,e) ^ z0(t,c0), получим систему равенств
F (со)
ФZ(z(-,C0), •, 0)()
PNd(Q*){J (z(;Oo), 0) - £L-Z (z (;во), •, 0)()}
0,
которые доказывают теорему. □
Если система (4.2) имеет некоторое решение с0 = с0 € И.р, то константа с0 определяет то порождающее решение г0(Ь,с}0), которому может отвечать решение г(Ь, е) исходной системы (2.1), (2.2), обращающееся в г0(Ь,с0) при е = 0.
Система операторных уравнений (4.2) аналогична известному в теории периодических нелинейных колебаний [1], [2] уравнению для порождающих амплитуд. Поэтому, в дальнейшем, будем ее называть системой уравнений для порождающих констант краевой задачи (2.1), (2.2).
Если система уравнений (4.2) не имеет решений, то краевая задача (2.1), (2.2) не обладает искомым решением. Таким образом, необходимое условие (4.2) не одно~ значно разрешимой краевой задачи может быть удовлетворено выбором элемента с0 = ср в семействе порождающих решений (3.12).
4.2. Достаточное условие существования решений
Найдем достаточные условия существования решения краевой задачи (2.1), (2.2). Выполняя в задаче (2.1), (2.2) замену переменных
г (Ь,е) = го(Ь,с0о) + х(Ь,е),
в которой константа с0 € И.р удовлетворяет системе уравнений для порождающих констант (4.2), приходим к следующей задаче: найти условия существования и алгоритм построения решения х(Ь,е) € Б^^е], обращающегося в нуль при е = 0 краевой задачи
Ьх(,е)(Ь) = £Z (го(Ь,с0) + х(Ь,е),Ь,е), , ,
£х(^,е) = еЗ (го(,с0)+ х(,е),е). 1 '
Используя свойства (2.3) нелинейных оператора Z(г,Ь,е) и функционала З(г, е), выделим у них линейные части по х и члены нулевого порядка по е. В
результате получим разложения:
Е (го(Ь,с}0) + х(Ь,е),Ь,е) = Ео(г,с}0) + Ьох( • ,е)(Ь) + Я(х(Ь,е),Ь,е), , ,
З(го(• + х(• ,е),е) = Зо(• ,с0)+ еох(• ,е)+ е^х(• ,е),е), 1 ' ;
где с0) = Е (го(Ь, с0),1, 0) : Б1 х! ^ Б2;
Зо( • ,с0) = ЗоЫ • ,с0, 0): Б1 х!^ И™;
Ь0 : Б1 х !е0 ^ Б2 х !е0 — линейный ограниченный оператор, представляющий собой производную Фреше нелинейного оператора Е (г,Ь,е) по ¿при г = г0(1,с}0)-,
е0 : Б1 х !е0 ^ И™ х !е0 — линейный ограниченный вектор-функционал, пред-ставляютции собой производную Фреше от нелинейного функционала З(г,е) по г при г = г0(Ь, с0);
Я : Б1 х!х!е0 ^ Б2 х!е0 — нелинейный оператор, удовлетворяющий условиям (а2) и (аЗ) из (2.3);
е1 : Б1 х! х !ео ^ И™ х !ео — нелинейный вектор-функционал, удовлетворяющий условиям (аб) и (а7) из (2.3).
Рассматривая нелинейности в краевой задаче (4.3) как неоднородности и применяя к ней теорему 1, получим для ее решения х(Ь,е) следующее представление
х(Ь,е) = Хр(Ь)с(е) + х{1\г,е), с(е) Е Ср[е\, с(е)
Ф{Ео( • ,4о) +
рмл^*){Зо( • ,с0) + еох( • ,е) + е1(х( • ,е),е) —
+Ьох( • ,е) + Я(х( • ,е)г ,е) = 0,
—еЬ-\Ео( • ,с0) + Ьох( • ,е) + Я(х( • ,е)г ,е)\ = 0.
Здесь Ср[е\ — пространство непрерывных по е функций со значениями в евклидовом пространстве Ир, с(•) : !е0 ^ Ир.
Неизвестная вектор-функция х(1 (Ь,е) определяется по формуле
х(1\г,е) = е\{СЕ (го( • ,с0)+ х( • ,е)г ,е))(Ь)+ X (1)Я+З (го( • ,с0)+ х( • ,е),е)\,
С
Используя разложения (4.4) и тот факт, что векторная константа со необходимо удовлетворяет системе уравнений для порождающих констант (4.2), для нахождения решения х(Ь,е) Е Б1Сп[е\ слабонелинейной краевой задачи (2.1), (2.2) приходим к эквивалентной операторной системе
х(Ь,е) = Хр(Ь)с(е)+ х{1)(г,е),
Ф{Ьох(1) (•, е)+
Во е(е) = -
PNd(Q*){tox(1)( • ,е)+ • ,е),е)-
+R(x(,e),^e)}
-iL-[Lox(l)(•, e) + R(x(, e), •,£)]}
(4.5)
x{l)(t,e) = e[(G{Z{zo{;cl) + Lo[XPt)c{e) + x{l)(•,e)]) + R(x(,e), ;e)})(t)+ +X (t)Q+{J (za(,cl)+ £0[Xp()c(e) + x{l)(,e)]) + il(x(, e), e)}],
Во
ФLn
PNd(Q*)[io + iL-Lo]}
Xp()
((и + х р)-мерная постоянная матрица. Обозначим через РМ(Во) : ^ N(В0) и Рмв^ : ^ УВо матрицы-
ортопроекторы, а через В+ — (р х (и + ^))-мерную псевдообратную матрицу к матрице В0.
4.3. Построение единственного решения
Предположим, что dimker В0 = 0, т.е. Рм(Во) = 0. Тогда второе уравнение операторной системы (4.5) разрешимо тогда и только тогда, когда выполнено условие
Р \ Ф[Ьох(1) ( • ,£) + Я(х( • ,£),•,£)}
мВоЧ Рм,(д*)[£ох(1)(• ,£)+ £г(х(• ,£),£) - £Ь-[Ьох(1)(• ,£) + Я(х(• ,£),• ,£)]}
при выполнении которого оно будет однозначно разрешимо
= 0, (4.6)
c(e) = -B,
+
Ф{Lox{l)(•, e) + R(x(• ,e), • ,e)} PNd(Q*){iox(l)( •,e) + il(x( • ,e),e) - iL-[Lox{l)(• ,e) + R(x(• ,e)r,e)]}
Пусть Pn(b*0)
Ф P
Nd Q)
0
РМ(Во) = 0 операторная система (4.5) примет вид:
х(г, £) = хр(ь)с(£) + х(1) (г,£),
c(e) = -B,
+
Ф{Lоx{l)(•, e)+ PNd{Q*){iox{l)(•, e) + il(x( • ,e),e)-
+R(x(• ,e)r, e)}
-iL-[Lox{l)( •, e) + R(x(•, e), • ,e)]}
x{l\t,e) = e[(G{Z (zo( • ,c») + L0[XP( • )c(e) + x{l)( • ,e)]) + R(x( • ,e)r ,e)})(t)+
{l)i
(4.7)
+eX (t)Q+{J (zo( •, c0) + io [Xp( • )c(e) + x{l)( •, e)]) + il(x( •, e), e)}].
Используя [1], [2], [4], [5], приведем операторную систему (4.7) к системе, для решения которой применим метод простых итераций. Введем следующие обозначения:
у(Ь,е) = оо1(х(1,е),с(е),х(1\1,е)) — вектор-столбец из банахова пространства В = В1Сп[е\ х Ср[е\ х Б^И;
(WM • ,e)(t) = -В+
ФЬо
PNdQ)[t0 - а-Ьо]
•,e)
ограниченный оператор;
W
О Xp(t) 1в1Сп[е]
0 0 Wi
ОО
0
— клеточно-матричный оператор верхнетреугольного вида, где 1В1сп[е] ственный оператор в пространстве В^^е^
(4.8)
тожде-
-В+
U (y(t,e),t,e) =
ФЯ(х( • ,e)r,e) PNd(Q*){^(x(•, e),e) - £Ь-Я(х(• ,e), • ,e)}
(i)i
е[(С{Е (го( • А) + Ьо[Хр( • )с(е) + х(1)( • ,е)\ + Я(х( • ,е)г ,е)Ш)+ [ +Х (1)Я+{З(го( • ,с0)+ £о[Хр( • )с(е) + х(1)( • ,е)\) + ¿1(х( • ,е),е)}\
векторная операторная функция, которая удовлетворяет условиям (а2), (аЗ) из (2.3).
Операторы, которые входят в клеточно-матричный оператор (4.8), лин6ины6 ограниченные по определению и как суперпозиции линейных ограниченных операторов.
Используя введенные обозначения, операторную систему (4.7) можно представить в виде
у(1,е) = (Му(• ,е))(1) + (Ы(у(• ,е), • ,е))(1). (4.9)
В силу структуры клеточно-матричного оператора М верхнетреугольного вида с нулевыми клетками на главной диагонали и ниже операторная система (4.9) преобразуется в систему
Уу(1,е)= Ы (у(1,е),1,е),
V = Ib - W
IB1C„[e] -Xp(t) -IBiC„[e]
0 0
I
BiC„[e] 0
I
-Wi
BiCn[e]
0
Оператор V — верхнетреугольный клеточный, с тождественными операторами на главной диагонали. По аналогии с [4] можно показать, что оператор V имеет ограниченный обратный. Поскольку оператор № — нелинейный и зависит от е, а оператор V-1 ограничен, то за счет выбора значений параметра е можно добиться того, чтобы суперпозиция операторов V-1№ была сжимающим оператором. Из принципа сжимающих отображений следует, что система (4.7) имеет единственное решение.
Используя метод простых итераций для нахождения решений краевой задачи
е
е = 0, получаем следующий итерационный процесс.
Приближения х)1+^1{1,е) к х(1\1,е) будем искать как частные решения краевых задач
Ьх^ • ,е)(1) = е^о(1,о1) + (Ь0[ХР( • )ск (е) + х^( • ,е)])(*) + К(хк (1,е),1,е)}, ¿х^•, е) = е{30(• ,с»)+ £о[Хр( • )ск(е) + хк\• ,е)] + ^(хк( • ,е),е)}.
(4.10)
По теореме 1 с учетом разложений (4.4) решение задачи (4.10) представимо в виде
хк+1(*,е) = е[(С^ (го( • ,с») + Ьо[Хр( • )ск (е) + х[1\• ,е)]) + Я(хк ( • ,е),-,е)})(*)+ + X (1)Я+{3 (го( • ,с0)+ £о[Хр( • )ск (е) + х^ (• ,е)])+ к\{хк ( • ,е),е)}].
Из необходимого и достаточного условия разрешимости краевой задачи (4.10), с учетом выбора с0 £ И.р из системы уравнений для порождающих констант (4.2), приходим к уравнению
BoCk(е) = -
Ф{Ьох{к1}(•, е) + Я(хк( • ,е),- ,е)} Рмл(д*){£охк1)( •,е) + к1(хк( • ,е),е) - кЬ-[Ьохк1)( • ,е) + Я(хк( • ,е)г,е)]
(4.11)
из которого находится к-е приближение ск (е) к с(е). Разрешимость систем (4.11) для каждого к обеспечивается выполнением условия
PN (B0)
Ф
PNd(Q*)
0,
а единственность — условием Рм(в0) = 0.
Тогда (к + 1)-е приближение хк+1(1,е) к x(í, е) запишется в виде!
хк+1(1,е) = Хр(г)ск (е) + х^^е), к = 0,1, 2,...,
где хо(Ь, е) = х^1 е) = 0 х1(Ь, е) = х1\1, е)
Таким образом, для краевой задачи (2.1), (2.2) справедлива следующая теорема.
Теорема 4. Пусть краевая задача (2.1), (2.2) удовлетворяет условиям (2.3), а соответствующая, порождающая краевая задача (2.4), (2.5) при выполнении условия (3.11) имеет семейство порождающих решений (3.12). Тогда для каждого элемента, с0 = с0 Е И.р; удовлетворяющего системе уравнений для порождающих констант (4-2), при выполнении условий
PN (Во) = ° PN (В*)
Ф
PNd(Q*) J
0
краевая, задача, (2.1), (2.2) имеет единственное решение г(Ь,е) непрерывное по е, обращающееся в порождающее решение г(Ь,с0о) при е = 0. Это решение можно найти с помощью сходящегося, на, [0,е*\ С 1ео итерационного процесса
гк+1(Ь,е) = го(Ь,с0) + хк+1(Ь,е), хк+1(г,е) = Хр(г)ск (е) + х^^^),
с {£) = _ в+\ Ф{рох{к)( • ,е) + Я(хк (• ,е),;е)}
0 [ рмл^){£охк\• ,е)+ ¿1 (хк(• ,е),е) — £Ь-[Ьох{1})(• ,е) + Я(хк(• ,е)г,е)\}
хк^1(г,е) = е[(С&Ы • ,с0) + Ьо[Хр( • )ск (е) + х^ • ,е)\)+
+Я(хк (• ,е)г ,е)}(1)+ Хр(г)д+{З (го( • ,с0)+
+£о[Хр( • )ск (е)+ хк1\• ,е)\) + ¿1(хк (• ,е),е)}\,
х0(г,е) = х{01)(г,е) = 0, х1 (г,е)= х[1)(г,е), к = 0,1, 2,....
Замечание 4. Подобные краевые задачи для всюду разрешимых систем обыкно-В6 н н ых дифференциальных уравнений рассмотрены в [1], [2], а для дифференциальных систем с сосредоточенным запаздыванием в [2], [3].
Список цитируемых источников
1. БойчукА. А. Конструктивные методы анализа краевых задач. — Киев: Наук, думка, 1990. - 96 с.
2. БойчукА. А., Журавлев В. <!>.. Самойленко А. М. Обобщенно-обратные операторы и нетеровы краевые задачи. — Киев: Изд-во ИМ 11Л11У. 1995. — 320 с.
3. БойчукА. А., Журавлев В. <!>.. Самойленко А. М. Линейные нетеровы краевые задачи для импульсных дифференциальных систем с запаздыванием // Дифференц. уравнения. - 1994. - 30, № 10. - С. 1677 - 1682.
4. Бойчук О. А., Панасенко 6. В. Слабконелшшш крайов1 задач1 для диференщальних р1внянь у критичному випадку у банаховому простор! // Нелшшш коливання. — 2010. - Т. 13, № 4. - С. 483 - 496.
5. Гребенников Е. А., Рябов Ю. А. Конструктивные методы анализа нелинейных систем. — М.: Наука, 1979. — 432 с.
6. Журавлев В. Ф. Критерий разрешимости и представление решений линейных n— (d—) нормальных операторных уравнений в банаховом пространстве //УМЖ. — 2010. - Т. 62, №2. - С. 167 - 182.
7. Канторович Л. В. Некоторые дальнейшие приложения принципа мажорант Ляпунова // Докл. АН СССР. - 1951. - 80, N◦ 6. - С. 848 - 851.
8. Канторович Л. В., ВулихБ. 3., Пинскер А. Г. Функциональный анализ в полуупорядоченных пространствах. — М.; Л.: Гостехиздат, 1950. — 548 с.
9. Крейн С. Г. Линейные уравнения в банаховом пространстве. — М.: Наука, 2071. — 104 с.
10. Ляпунов A. M. Общая задача об устойчивости движения. — М.: Гостехиздат, 1950. — 472 с.
11. Турбин А. Ф. Формулы для вычисления полуобратной и псевдообратной мтрицы // Журн. вычислит, математики и мат. физики. — 1974. — Т. 14, № 3. — С. 772 - 776.
12. BoichukA. A., Samoilenko А. М. Generalised inverse operators and Fredholm boundary-value problems. — Utrecht; Boston: VSP, 2004. — 317 p.
Получена 10.12.2011