Вычислительные технологии
Том 16, № 5, 2011
О разрешимости вырожденных систем квазилинейных интегро-дифференциальных уравнений общего вида*
Е. В. Чистякова, В.Ф. Чистяков Институт динамики систем, и теории управления СО РАН, Иркутск, Россия e-mail: [email protected], [email protected]
Рассматриваются вырожденные системы квазилинейных интегро-дифференциальных уравнений первого порядка общего вида. Доказаны локальная и глобальная теоремы существования и предложен численный метод решения.
Ключевые слова: интегро-дифференциальные уравнения, вырожденные задачи, нелокальная разрешимость, разностный метод первого порядка.
1. Введение и постановка задачи
Рассмотрим векторное уравнение
A(t,x, V>x,v)x + B(t,x,Vx,v) = 0, t е [0,1] = T, (1)
с начальными данными
x(0) = xo, (2)
где A(t,v,y,v) — (n х п)-матрнца, B(t,v,y,v) — n-мерная вектор-функция, определенные в области W = T х Rn х Rm х N, v — числовой параметр из отрезка
t
N = [—v0,v0] С R1, x = x(t, v) — искомая вектор-функция, V0x = J K0(t, s,x(s))ds,
o
t
Vx = / K(i, s, - Вольтсрра, K0(i, s, v), K(t, s, „) -^-морпыо вектор-
0
функции, определенные в области V = T х T х Rn, x0 — заданный вектор из Rn. Предполагается, что характер вырождения задается соотношением
det A(t,v,y,v) = 0 V (t,v,y,v) е W. (3)
В частности, в виде системы (1) можно записать систему дифференциальных и алгебраических уравнений и интегральных уравнений Вольтерра (первого и второго рода), связанных по части переменных.
Под решением системы (1) на отрезке T£ = [0, е] С T при значении па,рам,етра, v = v* будем понимать любую вектор-функцию x(t,v*) е C1(T£), которая обращает
T£
*Работа выполнена при финансовой поддержке Междисциплинарного проекта СО РАН № 107 и РФФИ (гранты № 09-08-00201 и 10-01-00571).
Интерес к уравнениям, описываемым в общем случае формулой (1), стимулируется проблемами математического моделирования в таких прикладных областях как теория электронных схем и электрических цепей, механика и химическая кинетика, гидродинамика и теплотехника [1-3]. Самым распространенным частным случаем этих уравнений являются дифференциально-алгебраические уравнения (ДАУ): системы А(Ь, х)Х + Б(Ь, х) = 0, удовлетворяющие условию (3). Для данного случая построена достаточно обширная качественная теория и разработан целый класс численных методов (см., например, [4-9]).
Если исследуемый процесс обладает последействием, то наряду с ДАУ модель процесса может содержать интегральные и интегро-дифференциальные уравнения. Ряд специальных случаев рассматривался в [10-15]. В работах [14, 15] существование решения постулируется и конструируются численные методы решения. Таким образом, вопрос о критериях разрешимости вырожденных задач вида (1), (2) достаточно актуален.
Основной задачей настоящей работы является получение достаточных условий сходимости численного метода для задачи (1), (2).
2. Вспомогательные результаты
Приведем ряд необходимых для дальнейших рассуждений сведений.
В работе используются нормы п-мерного вектора v = (^1 ... vn)т и (^хп)-матрицы V = (^, г = 1, 2,..., п, ] = 1, 2,..., вычисляемые по правилам
||VII = шах{^|, г = 1, 2,... ,п}, ¡V|| = шах |^ ^г = 1, 2,... | .
Под символами |^(т)|| и ¡V(т)|| понимаются нормы вектор-функции v(w) и матрицы V(ад), вычиненные в точке т € Б С И5, Включения v(w), V(т) € С 1(Б) означают, что все частные производные элементов вектор-функции v(w) или матрицы V(т) имеют непрерывные частные производные порядка до I включительно по всем компонентам вектора т в любой точке области Б. Непрерывности соответствуют включения v(w), V (т) € С (Б). Если v(w), V (т) € С(Б), то их нормы также непрерывны в Б.
Лемма 1 (см., например, [16]). Пусть для неотрицательных чисел иг, г = 2, 3,.. выполнено неравенство
г-1
иг < с + т^^щ, и1 < с,
3 = 1
где с, т — неотрицательные постоянные. Тогда, справедлива, оценка (разностный аналог леммы, Гронуолла — Беллмана)
иг < с(1 + т)г-1, г = 1, 2,...
Рассмотрим задачу (1), (2) в нормальной форме = Еп при некотором
фиксированном значении параметра V = V*
ж(Ь) = Б(Ь,ж(Ь)^ж), Ь € Т, (4)
х(0) = хо, (5)
где ^(¿, V, у) = — В(и, V, у, V*), и докажем следующие утверждения.
Лемма 2. Пусть для задачи (4), (5) выполнены условия гладкости
^(м,у) е С1(^), ^ = т X кга X и™, к(м^) е С1(У). (6)
Тогда, существует отрезок [0,£], £ > 0, на котором определено единственное решение задачи (4), (5).
Доказательство основано на выделении некоторой окрестности (шара) с центром в точке (0, х0, 0) в облаети переходе от задачи (4), (5) к системе интегральных уравнений
г
х(и) = х0 + У ^(в, х(в), Ух)^в о
и применении к последней метода последовательных приближений (Пикара) с буквальным повторением рассуждений из [17, с, 139] с учетом многомерности.
Лемма 3. Пусть для начальной задачи (4), (5) выполнено условие гладкости (6) и, кроме того,
1) Н^1 (М,у)|| < ко + К1|М| + К2||у|| V(t,v,y) е^,
2) ЦК(М^)|| < По + т|М| V(t,s,v) е V,
где Ki,nj, г = 0,1, 2, ] =0,1, — некоторые положительные постоянные.
Тогда, начальная, задача, (1), (2) имеет единственное решение из класса, С1(Т).
Доказательство, Мы не предполагаем выполнения глобального условия Липшица для вектор-функций ^(¿^,у) и К(¿,в^) и не можем воспользоваться методом последовательных приближений, как это делается, например, в работе [18, с, 393], Введем сетку {и = ¿К, г = 0,1, 2,..., N} С Т, К = 1/ЛГ.
Заменим в (4) производные конечными разностями х(^) ~ (х^ — xi)/К, а интегралы конечными суммами по формуле левых прямоугольников
х^1 — Xi = К^ ^ ,xi, К К (и ,xj . (7)
Сложив правые и левые части равенств (7), получим
х^1 = хо + К^^ I ¿к,хй,К^К(¿к,tj,Xj) I . (8)
к=0 \ j=о /
Учитывая условия 1 и 2 леммы, перейдем к нормам
i
||хт|| < ¿1 + К^2 ^ ||хк ||,
к=о
где ¿1 = 11х01 + к0 + к2п0, ¿2 = к1 + Отсюда по лемме 1
||хт|| < ¿1(1 + < ¿1(1 + ¿2^)^ < ¿1вЬ2.
Рассмотрим неравенство
||Хг+1|| + ||Хо|| < ¿1вЬ2 + ||жо||,
N+1 - Хо|| < |хг+1| + ||хо|| < ¿1вЬ2 + ||хо|| < 2Ь1вЬ2. (9)
Таким образом, согласно неравенству (9) все значения лежат в шаре радиуса 2Ь1еЬ2 с центром в точке жо. Соединим точки (¿¿,жг) отрезками. Получим семейство ломаных функций |ф(£)|, которые целиком содержатся в цилиндре
С = Т х (г : ||г - жо|| < 2Ь1вЬ2}, г € И",
в силу того что цилиндр является выпуклым множеством. По определению семейство функций (ф(^)} равномерно ограничено. Покажем, что оно равностепенно непрерывно, В силу своей гладкости Б(¿, ж, у) также непрерывна, а значит, достигает на Т своего наибольшего и наименьшего значений, т. е, найдется такая постоянная М > 0, что ||Б(¿, V, у) | < М У € Т, У- € И", у € И™. Следовательно, задав для аргумента £ приращение Д£, для любой функции из семейетва (ф(£)} можно записать
||ф(£ + Д£) - ф(£)| < Д£М = е.
Тогда, каким бы ни было е > 0, всегда можно принять Д£ = е/М так, что ||ф(£ + Д£) — ф(£)| < е. Таким образом, семейство функций (ф(£)} равностепенно непрерывно. Следовательно, по теореме Арцеля [19] из этого семейства можно выбрать равномер-Т
функции
ф(1\ф(2\...,ф(г\... (10)
Предельную функцию обозначим через ф(£), £ € Т. Покажем, что она является решением начальной задачи (4), (5),
Согласно (8) подпоследовательности ломаных (10) соответствуют равенства
жф+1 = жо + Ы ¿ Б ( ¿к, жф, Ы ^ К (¿к, Ъ, жфн . (11)
к=о \ ^=о /
Перейдем в (11) к пределу при I — то (Ы — 0), Получим
г / к
Ы/У Б | ¿к ,ф(4
<^(£¿+1) = жо + Ит Ы ^ Б I ¿к, ), Ы ^ К (¿к, ^, )) ™ к=о \ ^=о
В силу непрерывности вектор-функций ф(^), Б(¿,у,у) и К(¿,8, у) их суперпозиции также непрерывны. Для произвольного ^ согласно определению интеграла Римана можно записать
Ф(*) = жо +1 б у к
Сравним правую и левую части данного равенства, В правой части ф(^) стоит под знаком интеграла, следовательно, ф(^) € С1(Т), Таким образом, ф(^) удовлетворяет интегральному уравнению, эквивалентному начальной задаче (4), (5),
Докажем единственность решения Пусть через точку (t0,0(t0)), t0 £ T, проходит еще одно решение задачи (4), (5), являющееся решением начальной задачи
y(t) = F + K(t,s,y(s))dsj , y(to) = 0(to), (12)
to
где = / K(i,Начальная Wa (12) согласно лемме 2 нмеел еди„слвен„„е
o
решение на отрезке [t0, t0 + е]. Получено противоречие. Лемма доказана.
Замечание 1. Доказательство леммы, не меняется, если под символом, Vx понимать набор операторов Вольтерра Vx = (Vix, V2x, ..., V^x) с линейным условием роста, ядер.
Определение 1 (см., например, [4]), Матрица, обозначаемая, в дальнейшем как S-, называется, полуобратной к матрице S, если, она, удовлетворяет уравнению SS -S = S.
S
ли выполнено условие Кронекера — Капелли rank S = rank(S|u), то си,стем,а, уравнений Sy = u разрешима и все ее решения описываются, формулой
y = S-u + [E - S-S]C,
где E — единичная матрица подходящей размерности, C — произвольный вектор.
Определение 2. Пучок квадратных матриц AA(w) + B(w), w £ D С Rq, где A — скалярный параметр (в общем случае комплексный), удовлетворяет критерию "ранг-
D
1) max rank A(w) = r, w £ D,
2) det[AA(w) + B(w)] = Ara0(w) + где a0(w) = 0 Vw £ D.
Лемма 5. Пусть пучок постоянных матриц А (А) = А А+В регулярен det А (А) ф 0. Тогда, справедливо неравенство rank А > degdet А(А), где deg — сим,вол, степени многочлена.
Доказательство, Если пучок матриц A(A) регулярен, то существуют постоянные неособенные матрицы P, Q со свойством
pawq = a(«< + ¿d),
где Nk = 0, J — некоторый блок размерности (d х d) [21, с, 354], Отсюда
deg det А (А) = deg [det (АД* + J) det(AiV + En_d)] = deg det (АД* + J) • 1 = d. Следовательно, d < rank A, d = rank А тогда и только тогда, когда N = 0,
Следствие 1. Для, пучка матриц AA(w) + B(w), w £ D, удовлетворяющего критерию "ранг-степень", справедливо равенство rank A(w) = const = r V w £ D.
Лемма 6 [201. Матричный пучок Р(А; ¿) = А ( А1(и,) | + ( , т € Б, где блок
\ 0 ) \В2(т))
А1(т) имеет полный ранг для любого т € Б, удовлетворяет критерию "ранг-степень" тогда и только тогда, когда
(А; т) = Агао(т) + ..., = ае^ А^}) =0 У ™ € Б.
3. Теоремы о разрешимости вырожденных систем
Системы вида (1), удовлетворяющие условию (3), обладают сложной внутренней структурой, Мерой сложности является целочисленная характеристика, называемая индексом — минимально возможным порядком интегро-дифференциального оператора, действуя которым на (1) получаем систему с невырожденной в окрестности начальной точки матрицей при производной искомой вектор-функции.
Теорема 1. Пусть для задачи, (1), (2) при = А(^), V = 0, выполнены
условия:
1) ¿(¿) € С2(Т), В(*,у,у, V) € С2^), К(М,у) € С1(У),
2) rank A(0) = rank
3) ||dB(t,v,y,v)/dv
A(0)| b , b = -B(Z), z = (0,xo, 0, 0),
< L, ||dB(t,v,y,v)/dy|| < L V(t,v,y,v) G W, ^ ||B(t,v,y,v)/|| < Ko + Ki||v|| + K2|y| V(t,v,y,v) G W,
5j ||K(t,s,v)|| + ||dK(t, s, v)/dt| < no + ni||v|| V(t,s,v) G V, где L, щ,Пз, i = 0,1, 2, = 0,1 — некоторые положительные постоянные, 6) старший коэффициент многочлена
det
д
AA(i) + — B{t,v,y,u)
ao(t, v, y, v)Ar + ...,
gcte r = max {rank A(t), t G T}, удовлетворяет условию
|ao(t,v,y,v)| > Co > 0 V(t, v, y, v) G W.
Тогда,
1) индекс системы, равен 1,
2) начальная, задача, (1), (2) имеет единственное решение x*(t, 0) из класса, C1(T).
Доказательство, Согласно лемме 5 и условию 6 теоремы rank A(t) = const = r V t G T. Тогда существует матрица P (t) G C1(T) со свойствами
det P(t) = 0 Vt G T, P(t)A(t) = ^ Al0(t)
где блок A1(t) имеет размерность (r x n) [20], Формально подействуем на систему (1) оператором вида
Получим новую систему
где х = (t, х, Vx, 0),
A(t,v,y,v)= ( дВ2^у,иЛ , VlX=j dK^^a))d8i (14)
dv
0
P(t)B(t, v, y, V) = Bi(t, v, y, V), P2(t)B(t, v, y, V) = B2(t, v, y, v),
B2(t, ж, Fx, Fx, //) = ^^ + [K(t, t, x) + Fx]. (15)
dt dy
Покажем, что матрица, стоящая перед производной искомой вектор-функции в (13), обратима и
||A-1(t,v,y, v)|| < Li, V (t,v,y,v) G W, Li = const. Согласно лемме 6 и условию 6
| det A(t, v, y, v)| = | det P(t)a0(t, v, y, v)| > c01 det P(t)|
и все элементы матрицы A(t, v, y, v) ограниченны на W (см. условие 3 теоремы). Из условий 3-6 следует, что правая часть системы по своим аргументам растет не быстрее, чем линейно.
Приведенная к нормальной форме система (13) удовлетворяет условиям леммы 3, и задача (13), (2) имеет единственное решение x*(t, 0) G C1(T), Покажем, что оно удовлетворяет исходной задаче (1), (2). Подставим эту вектор-функцию в исходную задачу. Условие 2 теоремы удовлетворяется автоматически. Предположим, что
A(t)X* + B(t,x*,Vx*, °) = ^(t), t G T,
где ^(t) = 0. В силу условия 2 вектор-функция x*(t, 0) удовлетворяет уравнению (1) в точке t = 0, Следовательно, ^(0) = 0. Согласно формуле (13)
Pi{t)\,d_( 0 о у dt \ p2(t)
^ = 0, t G T,
где detP(х) = 0 Vt G T. Произведем замену ^ = P 1(t)h(t^. Получим на T систему вида
E 9h + (0 EL) h = ()- 'г<()> = а <16>
Теорема доказана. □
Для общего случая доказана только локальная теорема существования.
V=0
1) A(t,v,y,v) G C2(W), B(t,v,y,v) G C2(w), K(t,s,v), Ko(t,s,v) G C1(V), rank A(t,v,y, v) = r = const e окрестности mочки Z = (0,x0, 0, 0),
rank A(Z) = rank |~A(Z)| Ь , b = -B(Z),
4) гапк) = degdet[AA(C) + £«)], где
д
а»
, (17)
здесь с — одно из решений системы, )с = — £(£). Тогда,
1) индекс системы, равен 1,
2) существует отрезок [0,е] С Т, на котором, определено единственное решение я* (г, 0) задачи (1), (2).
Доказательство, В силу условий 1, 2 теоремы в окрестности Wo С ^ ^отки ( найдется неособенная матрица Р(г, V, у, V) € C2(W0) такая, что
Р (М,.^ )А(Мм/^)=( ^^
Формально подействуем на систему (1) оператором вида
здесь х = (г, я, Уж, 0), Получим новую систему
А<*>ж« + ( ^№*ж,*))=0- <18>
где
Л(г,ь,у,и)= I дВ2{г,у,у,у)
Р1(г,г,У ^ = В1(г,г,У p2(t,v,y,v)в(t,v,y, v) =
(Обозначения для символов VI, -В2(г, ж, Уж, У1ж, V) смотри в формулах (14), (15),) Покажем, что матрица, стоящая перед производной искомой вектор-функции в (18), невы-рожденна в начальной точке Учитывая, что Р2В = В2, имеем
/ А
\ дг дг
Докажем, что
Л(0 = Щ(), 7г= (р9В Л' д ]. (19)
Для упрощения записи указание зависимости от г, здесь и ниже опущено. Рас-
смотрим равенство
Для выполнения (19) необходимо показать, что
в начальной точке Заметим, что Р2(АЖ) = 0 У(£,г,у, V) € W0, следовательно,
д дР д дР д
— [Р2Ах] = Ах + Р27Г(Ах) = 0, —1Ах = -Р2—(Ах). дг дг дг дг дг
дР2
В начальной точке ( произведение ——Ах не меняется, если вектор ¿(0, 0) заменим
дг
произвольным решением уравнения )с + ) = 0, которое в силу условия 3 теоремы разрешимо. Тогда в точке ( справедливо равенство
Согласно условию 4 теоремы пучок из (20), имеющий в начальной точке вид АА(£) + Б(С), удовлетворяет критерию "ранг-степень", и из леммы 2 следует, что det А(С) = 0, Следовательно, по непрерывности существует окрестность точки (, в которой матрица А обратима, дифференцируема и согласно лемме 5 на некотором отрезке [0,е] С [0,1] существует решение начальной задачи (18), (2) ж*(г, 0),
Докажем, что ж* (г, 0) является решением задачи (1), (2), Подставим эту вектор-функцию в исходную задачу. Условие (2) удовлетворяется автоматически. Предположим, что
А(г,ж*, Уж*,0)ж* + в(г,ж*, Уж*,0) = ^(г), г € [0,е],
где = 0, В силу условия 3 теоремы вектор-функция ж* (г, 0) удовлетворяет уравнению (1) в точке г = 0. Следовательно, ^(0) = 0, Согласно формуле (18)
Q оц
рЛх*)\,1( о о rdt\P2(x*)
ц = 0, t G [0,е],
где х* = (t,x*,Vx*, 0), detP(t) = 0 Vt G [0,е]. Произведем замену ц = P 1(x*)h(t), Получим на отрезке [0, е] задачу вида (16), Теорема доказана, □
Замечание 2. Если условие 4 выполнено для одного с : A(Z)c+B(Z) = 0, mo оно выполнено и для, всех других с. Согласно лемме 4 в (11) можно принять c = — A-(Z)B(Z)• Частным случаем S- является псевдообратная матрица S + [21], которую можно вычислить как предел G(5) = (52E + STS)-1ST, 5 — 0, где ||G(5) — S + || = O(52) [4].
Следствие 2. Если, на, T существует решение x*(t, 0) системы (1), (2) и выполнено условие
rank А(х*) = deg Ф(А; t) = const V t G T, где х* = (t,x*,Vx*, 0), e качестее с принято X*(t, 0),
Ф(А; t) = det [АА(х*) + B(x*)] , (21)
mo существует окрестность начального данного Q = {w : ||x0 — w)|| < q} такая, что при достаточно малом q на T найдется решение системы (1) x*(t, 0) : x*(0, 0) G Q. Здесь II — равномерная норма в Rn.
Доказательство, Докажем вначале локальный вариант следствия. Из (21) и рассуждений при доказательстве теоремы следует, что определитель (18) не равен нулю для любой точки £ € Т. Следовательно, по непрерывности существует окрестность точки х*, в которой матрица А невырожденна. Таким образом, в этой окрестности систему (1) можно привести к нормальной форме и для нее применить схему доказательства теоремы о погружении [22, с, 36],
Следствие 3. Пусть на отрезке Т определено решение х*(£, 0) задачи, (1), (2) при значении параметра, V = 0. Тогда при достаточно малом значении величины, ¡V| в окрестности х*(£, 0) определено решение задачи, (1), (2).
Отметим еще одно обстоятельство. Пусть на Т существует решение х*(£, 0) системы (1), (2) и многочлен из (21) имеет вид
Ф(А; £) = ас(£)Аг + € Т, (22)
где г = шах(гапкА(х*), £ € Т}, Если уравнение а0(£) = 0 имеет изолированное решение £* : а0(£*) = 0, то £* может являться особой точкой системы (1), Поясним сказанное. Уравнение
г, V0. + ф, V)г + г, V0. + ф, V) = 0, £ € [£0, £0 + 8] с начальными данными х(0) = х(£0), где
¿о г
ф^ф<() = /V-0. = /
о ¿о
при достаточно малом 8 > 0 имеет решение, совпадающее с х* (£, 0) на отрезке [£0, £0 + 8], и удовлетворяет условиям теоремы 1, При £0 — £* видим, что £0 + 8 — £*, Следовательно, в формуле (18) А(х*)к+г^г* — 0, гДе X* = (£,х*, Ух*, 0), и
А(х*)|г* =0, А(х*) = 0, £ = £*.
Численные способы поиска особых точек обсуждаются ниже.
Класс систем, описываемый теоремами 1, 2, не исчерпывает всего класса систем, имеющих индекс 1,
Пример 1. Рассмотрим задачу
Г«.X,*,УХ)= ( 0 )( X2 ) + ( * ) =0, х(0)=0,
г г
где = / „«Л, у2х = / штора, ,шеет еОиисп.еииое Решеиие = 0.
00 Для любой вектор-функции, х € С2(Т) имеем,
. , . , /х 1 + . /0 - (х^х + хгх 1) \ / 1 -х 1^х \ (
Л1о^,х,х,Ух)=( ^ , Л1=1 0 ) + (о 1
Отметим, что широко используемый метод дифференцирования конечных связей, в частности в механике, здесь не дает результата: после первого шага процесса получим систему
ж^ж 0 \ / же 1 \ /ж! + У2ж
0 1 у \ ^ \ 0 ,0
с особенной в начальной точке матрицей при х. Как обратить эту матрицу, так и сделать новый шаг процесса мы не можем.
4. Выбор численного метода
Отметим, что применение явных схем (предпочтительных здесь в силу простоты их реализации) для вырожденных систем интегро-дифференциальных уравнений в настоящее время представляется затруднительным.
Пример 2. Рассмотрим задачу
А-+вж=(; 0)(у+(; 1)(ж0=/<«>=(®■<€т, (2з)
где (ж1(0) ж2(0) )т = (с /2(0) )т , с — произвольное число.
Явная схем,а, Эйлера, А(жг+1 — жг)/Л + Вжг = /(гг), где жг — значения сеточной функции, Л = 1/Ж, и = гЛ, г = 0,1,..., N — 1, после простого преобразования, имеет вид
1 1 \ / ж1,г+Л = / 01 Л = / жм + ж2,г — Л [ж 1,г + ж2,г — Л^П (ОД)
0 ж2,г+^ V 02,гУ V ж^ — /2 (¿г) У' ( )
Система (24) имеет множество решений, и не известно, как в нем, выбирать решение, нужное для, сходим,ости. Попытаемся выбрать какое-то обобщенное решение, например нормальное псевдорешение
-+■ = А+(£) ■ А+ = ^ , (25)
где А + — псевдообратная, матрица (см. замечание 2). В выражении (25) обе компоненты вектора жг+1 совпадают, а в си,стем,е (23) они, в общем случае различны.
Здесь можно воспользоваться, следующим обстоятельством,. Система (23) имеет явную алгебраическую связь ж2 = /2(г). Учитывая, что ж2,г = /2(гг), (ж2,г+1 — ж2,г)/Л ~ [/2(гг) — /2(гг-г)]/Л, г > 1, получим сходим,ость. В общем случае системы (1) такой подход не всегда, возможен.
г
ную систему. Для упрощения используют прием линеаризации. Известно (см., например, [23]), что для задачи в нормальной форме х + В (г, ж) = 0 неявную схему Эйлера (жг+1 — жг)/Л + В(гг+1, ж г+1) = 0 можно заменить сходящейся с первым порядком схемой (жг+1 — жг)/Л + В(гг+1, жг) + 3(£¿+1, жг)(жг+1 — ж;) = 0, 3( ж, £) = дВ(£, ж)/дж, ГД6 При г
Аналогично на отрезке [0,е] из теоремы 1 для задачи (1), (2) введем сетку гг = гЛ, Л = е/^ N — число узлов сетки. Запишем неявную схему
— ( г+1
А(У1+1)Хг+1к Хг+В(¥г+1) = 0, Уг+1 = \^г+1,хг+1,Н^2к(и,з3,х3),и
в которой проведем разложение входных данных в ряд Тейлора похс удержанием первого члена разложения. Интеграл будем считать не по формуле правых прямоугольников, а по формуле левых прямоугольников. Получим схему следующего вида:
где
А + + Аг) (хт - жг) = — ЛВ, г = 1, - 1,
А = А(^), В = В(^), Zi = ,х.),И ,
(26)
.7=0
В
ди
Л
ди
Р
{ftV,y,v)=Zi
х xi— 1 Л
г
Теорема 3. Пусть
1) выполнены условия теоремы, 2,
2) известны, начальные значения х(0) = х0 и х(Л) = х1.
Тогда, для схемы (26) равномерно по г справедлива, оценка, ||xi — х(^)|| = О(Л). Кратко изложим общую схему доказательства. Действуя на (26) оператором вида
0=
А(^) 0
А
0
ВД)
где А — разностный оператор: А^ = (д^ — получим разностную схему, с точ-
ностью О(Л) аппроксимирующую систему (18), у которой матрица перед производной обратима в условиях теоремы 1 в окрестности начальной точки. Несложные, но громоздкие вычисления показывают, что для невырожденной задачи разностный метод сходится, а следовательно, сходится и для вырожденной.
Пример 3. Рассмотрим тестовую задачу вида (1), (2)
' } ^
е *х1(^) е "(х2(£) + х^^х^в)^)
х1(*)х2(*) х2(*) + х2(£) / х1(в)х2(в)^5
V
хх (¿) —
/
\
2е24х2(*) + е х^) + 2е24х2(:£) у А^)^
0
0, Т =[0,1],
где при V = 0 х(£)
е4
24 > х0 = х(0)
1
е1 теорем,ы, 2. Имеем, V = ^2)т, у = (у1 у2)т € И2
Убедим,ся, в выполнении условий
)=( е—4 е—+ У1) 1 VlV2 VI + V2Уl
В(^У *) — ( 2VlV2 + е"/3 + 2^/3 + ^
г
г
Очевидно, что первое условие выполнено. Далее
det A(t, v, y, v) = 0 V (t, v, y, v), rank A(t, v, y, v) = 1 e окрестности точки Z = (0,x0, 0, 0). Условие 3 теоремы, 2 выполнено:
11 11
rank
1 1 3 1 1 3
rank A(Z) = rank [A(Z)| b] = rank Уравнение для вектора c из (4) имеет вид
1 11 С = ( 31
1 1) С \3 у ~ \3 - С1
где c1 — произвольное чиело из R1. Условие 4 также выполнено: 1 1\ . / -14/3 + ci 1 - ci
1.
c=
c1
deg det
Л
11
1
1 - 2d
deg [(-17/3)Л + ...] = 1 V ci.
Таким, образом,, существует окрестность нуля, в которой определено решение нашего примера. Проверим, задачу на наличие особых точек. Выпишем многочлен из (21):
Ф(Л; t) = det
1 (302 + 03 - 1)/30\ + Л|03 tf2(3tf2 + 03 - 1)/3] +
1-2в2- §03 0 203 0
ао^)Л
где
ao(t) = -
14
з V з
-1
= е.
Очевидно, что старший коэффициент a0(t) = 0 V t G T и условия следствий 2 и 3 выполнены.
Были произведены, расчеты, по схеме (26). Результаты приведены, в таблице (здесь Err = max ||x - x(tj)|| — глобальная, погрешность).
N Err (v = 0) Err (v = 0.2)
100 0.11 0.11
1000 0.011 0.077
10000 0.0011 0.075
В заключение сделаем замечания о поиске особых точек. При аномалиях счета (например, авосте) для систем небольшой размерности можно попытаться определить наличие особых точек решения на Т. Для этого следует по точкам х^ г = 0,1,...которые удалось вычислить, построить многочлен ¿(¿), подставить его в многочлен (21) и вычислить аг (¿) — приближение для аг (¿) из (22),
Пусть ¿(¿) близок к решению задачи (1), (2) в С1[0,+ е], е — некоторое число, большее К, Если нули у аг (¿) простые, то у функции аг (¿) при достаточной малости нормы ||аг(¿) — аг(¿)||с1[07Ъ+е] существуют близкие нули, Нулям функции аг(¿) соответствуют особые точки решений.
г
Список литературы
[1] Ушаков Е.И. Статическая устойчивость электрических систем. Новосибирск: Наука, 1988.
[2] Серов Е.П., Корольков Б.П. Динамика парогенераторов. М.: Энергоиздат, 1981.
[3] Muller Р.С.: Aspects of modeling dynamical systems by differential-algebraic equations // Math, and Computer Modelling of Dynamical Systems. 2001. Vol. 7. P. 133-143.
[4] Бояринцев Ю.Е. Регулярные и сингулярные системы линейных обыкновенных дифференциальных уравнений. Новосибирск: Наука, 1980.
[51 Griepentrog е., Maerz R. Differential-Algebriac Equations and Their Numerical Treatment. Leipzig: BSB B. G. Teubner Verlags gesellschaft, 1986. 220 p.
[61 Brenan K.E., Campbell S.L., Petzold L.R. Numerical Solution of Initial-Value Problems in Differential-Algebraic Equations (Classics in Applied Mathematics; 14). Philadelphia: SIAM, 1996.
[71 Rheinboldt W.C., Raiber P.T. Theoretical and Numerical Analysis of Differential-Algebaric Equations. Handbook of Numerical Analysis: Vol. VIII. Amsterdam, 2002.
[81 Хайрер Э., Bahhep Г. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи. М.: Мир, 1999.
[91 Kunkel P., Mehrmann V. Differential-Algebraic Equations: Analysis and Numerical Solution. Europ. Math. Soc., 2006. 337 p.
[101 Булатов M.B. Об интегро-дифференциальнызх системах с вырожденной матрицей перед производной // Дифференциальные уравнения. 2002. Т. 38, № 5. С. 692-697.
[Ill Булатов М.В., Чистякова Е.В. Численное решение интегро-дифференциальных систем с вырожденной матрицей перед производной многошаговыми методами // Там же. 2006. Т. 42, № 9. С. 1248-1255.
[121 Bulatov M.V., Chistyakova Е.V. Numerical solution of singualr systems of integral differential equations // Proc. Intern. Conf. on Comput. Math. Novosibirsk, 2004. Vol. 2. P. 813-817.
[131 Горбунов В.К. Метод нормальной сплайн-коллокации // Журн. вычисл. математики и матем. физики. 1989. Т. 29, № 2. С. 212-224.
[141 Gorbunov V.K., Perischev V.V., Sviridov V.Y. Development of the normal spilne method for linear intego-differential equations // Proc. ICCS-2003. Part II. Springer, 2003.
[151 Дмитриев С.С., Кузнецов Е.Б. Численное решение интегродифференциально-алгебраических уравнений с запаздывающим аргументом // Журн. вычисл. математики и матем. физики. 2008. Т. 48, № 3. С. 430-414.
[161 Апарцин А.С. Неклассические уравнения Вольтерра I рода: Теория и численные методы. Новосибирск: Наука. Сиб. издат. фирма РАН, 1999.
[171 Васильева А.Б., Тихонов Н.А. Интегральные уравнения. М.: Изд-во Московского гос. ун-та, 1989.
[181 Далецкий Ю.Л., Крейн М.Г. Устойчивость решений дифференциальных уравнений в банаховых пространствах. М.: Наука, 1972. 534 с.
[191 Петровский И.Г. Лекции по теории обыкновенных дифференциальных уравнений. М.: Наука, 1964. 272 с.
[20] Чистяков В.Ф. Алгебро-дифференциальные операторы с конечномерным ядром. Новосибирск: Наука, 1996.
[21] Гантмахер Ф.Р. Теория матриц. М.: Наука, 1967.
[22] Ахиезер Н.И. Лекции по вариационному исчислению. М.: Гос. изд-во технико-теор. литры, 1956.
[23] Калиткин H.H. Численные методы. М.: Наука, 1978. 512 с.
Поступила в редакцию 2 июля 2010 г.