УДК 519.62
БОТ: 10.14529/шшр140310
МЕТОДЫ ТИПА АДАМСА ДЛЯ РЕШЕНИЯ ВЫРОЖДЕННЫХ ИНТЕГРО-ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
М.В. Булатов, До Тиен Тхань
В статье рассмотрены линейные интегро-дифференциальные системы уравнений первого порядка с тождественно вырожденной матрицей перед производной. Для данных систем задано начальное условие, которое предполагается согласованным с правой частью. Рассматриваемые в статье постановки задач возникают при математическом моделировании сложных электрических цепей. Используя аппарат матричных полиномов, выделен класс задач, имеющих единственное решение. Обсуждаются трудности численного решения таких задач, в частности неустойчивость многих неявных методов. Для численного решения такого класса задач предложены многошаговые методы, которые основаны на явной квадратурной формуле Адамса для интегрального слагаемого и на экстраполяционных формулах. Сформулированы достаточные условия сходимости таких алгоритмов к точному решению.
Приведены результаты численных расчетов, которые хорошо согласуются с теоретическими выкладками.
Ключевые слова: многошаговые методы; интегро-дифференциальные уравнения; матричные полиномы.
1. Постановка задачи и необходимые сведения
При математическом моделировании сложных электрических цепей возникают взаимосвязанные интегро-дифференциальные и интегральные уравнения Вольтерра второго и первого родов [1, 2]. Если данные уравнения объединить, то получим систему интегро-дифференциальных уравнений с тождественно вырожденной матрицей перед производной.
Настоящая работа посвящена численному решению системы
I
А(Ь)ж (Ь) + В(£)ж(£) + У К(Ь, з)ж(з)^5 = /(Ь), Ь € [0,1], (1)
о
с начальным условием
ж(0) = жо, (2)
где А(Ь),В(Ь),К(Ь, з) - (п х п)-матрицы, /(Ь),ж(Ь) - п-мерные известная и искомая вектор-функции. Элементы А(Ь), В(Ь), К(Ь, 8), /(Ь) предлагаются достаточно гладкие.
В статье рассмотрен случай, когда матрица А(Ь) тождественно вырожденная, т.е.
det А(Ь) = 0 УЬ € [0,1]. (3)
Под решением системы (1) мы будем понимать любую вектор ж(Ь) € С1(Т), которая обращает исходную систему в тождество и удовлетворяется начальным условием (2).
Систему вида (1) с условием (3) будем называть системой интегро-дифференциальных уравнений с тождественно вырожденной матрицей. Задачи (1) - (2) при выполнении условия (3) принципиально отличаются от систем интегро-дифференциальных уравнений, разрешенных относительно производной. В частности, для существования решения системы (1) необходима разрешимость системы линейных алгебраических уравнений вида
А(Ьо)ж'(Ьо) = /(Ьо) - В(Ьо)жо
относительно вектора х (£о). Для этого необходимо и достаточно выполнения условия:
гапк(А(£о)) = гапк(А(^)|/(£о) - В(£о)хо).
Даже если решение задачи (1) - (2) существует, то оно может быть неединственным. Ниже мы приведем достаточные условия существования единственного достаточно гладкого решения рассматриваемой задачи и обоснуем многошаговые методы ее численного решения.
Приведем класиффикацию рассматриваемых задач:
• Если К(£,т) тождественно нулевая матрица, то такие системы принято называть дифференциально-алгебраическими уравнениями (ДАУ);
• Если А(£) тождественно нулевая матрица, detB(t) = 0 В(£) = 0, то мы будем иметь интегро-алгебраическое уравнение (ИАУ);
• Если А(Ь) и В(Ь) тождественно нулевые матрицы, то мы будем иметь систему интегральных уравнений Вольтерра первого рода (СИУВ).
Степень сложности систем вида (1) определяется с помощью характеристики, называемой индексом. По аналогии с [3], будем пользоваться следующим определением индекса:
Определение 1. Минимальное число г, при котором существует квазилинейный дифференциальный оператор
rd fir = £ Wj (t)( )j,
r
j=o
где Wj (t) — (n x n) - матрицы из C(t), со свойствами
t
fir о (A(t)x (t) + B(t)x(t) + J K(t, s)x(s)ds) =
A(t)x (t) + B(t)x(t) + J K(t, s)x(s)ds = fir о /(t),
o
причем
A(t) = О, Vt Є [0,1],
назовем индексом, системы (І).
Приведем необходимые определения и вспомогательные результаты для дальнейших исследований:
Определение 2. [ІІ] Матричный полином ЛA(t) + ^B(t) + C(t), t Є [0,1], где A(t),B(t),C(t) — (n x n) матрицы, Л, ^ - скалярные параметры, имеет простую структуру, если выполнены условия:
1) rankA(t) = const = r, Vt Є [О, І];
2) rank[A(t)|B(t)] = r + l = const Vt Є [0,1];
3) det^A(t) + ^B(t) + C(t)] = ao(t^r^l + ..., причем ao(t) = О Vt Є [0,1].
Теорема 1. [3] Пусть для задачи (1), (2) выполнены условия:
1) Элементы A(t),B(t),/(t),K(t, s) принадлежат классу Cp+j;
2) Матричный полином ЛA(t) + ^B(t) + K(t, t) имеет простую структуру;
o
t
3) rankA(0) = rank(A(0)|/(0) — B(0)x(0)).
Тогда задача (І), (2) имеет на отрезке [0,1] единственное решение x(t) из класса C^ 1]. Лемма 1. [б] Пусть r-мерные векторы Qi вычисляются по правилу
г-1
Qi = CQi-1 + h Di,iQi + Fi,i = 2,3,..., N, h = 1/N, (4)
1=1
где векторы Qo,Ql, ...Qk-l - заданы, ||Di , i У < Li < то, || Fi У < L2 < то и собственные значения матрицы C удовлетворяют условию ^j | < 1, и на границе единичного круга нет кратных собственных чисел.
Тогда справедлива оценка
||Qi|| < L3||^i||, i = 2, 3,...,N,L3 < то,
где ||Фг|| = max{|Fi|, ||Ql||}.
Лемма 2. [ІІ] Если матричный полином ЛA(t) + ^B(t) + C(t) имеет простую структуру на всем отрезке [0,1], и элементы A(t),B(t),C(t) принадлежат классу функций C^ ф то существуют невырожденные матрицы P (t) и Q(t) с элементами из C^ 1] такие, что
/Er 0 0\ /Bl(t) 0 B3(t)\
P(t)A(t)Q(t) = 0 0 0 ,P(t)B(t)Q(t) = I 0 Ei 0
\ 0 0 0/ V 0 0 0 J
/Cl(t) C2 (t) 0 \
P (t)CQ(t) = C4 (t) C5 (t) 0 ,
V 0 0 En-r-J
где Em - единичная матрица размера m.
Из леммы 2 следует
Лемма 3. Если матричный полином ЛA(t) + ^B(t) + C(t) имеет простую структуру на всем отрезке [0,1], и матрицы A(t), B(t) и C(t) имеют блочный вид
/Al(t)\ /Al(t) A2(t) Aa(t)\ /Bl(t)\ /Bl(t) B2(t) B3(t)\
A(t) = 0 = 0 0 0 , B(t) = B2(t) = B4(t) B5(t) Be(tn ,
V 0 / V 0 0 0 ) V 0 / V 0 0 0 )
/Cl(t)\ /Cl(t) C2(t) C3(t)\
C (t) = C 2(t) = C4 (t) C5 (t) Ce(t) ,
VC3(t)/ \C7(t) Cs(t) Cg(t)/
где Al(t), Bl(t), Cl(t) - матрица размер (r x n), B2(t), C2(t) - матрица размер (l x n), C3(t) - матрица размер (n—r—l)xn) и rank(Al(t)) = r = const Vt Є [0,1], rank((A(t)|B(t)) = r + l = const Vt Є [0,1], то
/ Al(t) \ det qB2(t) =0, Vt Є [0,1],
\sC 3(t)/
где q = 0, s = 0.
Вначале приведем описание численных методов решения интегральных уравнений Вольтера I рода, основанных на явных квадратурных формулах типа Адамса.
Зададим на отрезке [0,1] равномерную сетку ^ = гЛ, г = 1, 2, Н = 1/Ж, и пусть
известны значения у ~ х(^). Тогда
к г к г
= Н X)в+ Е Н Е 7гУ-г = Н Е ^г+1 >(5)
г=о ^=к+1 г=о г=о
где Ьк+1(уг-к, Уг-к+1, •••, Уг, £) - интерполяционный полином степени к, проходящей через точки (уг-к, ^г—к), (Уг-к+1, ^-к+1), •••, (УгЛ). Выпишем коэффициенты 7г (см., например [6]) в таблице 1.
Таблица 1
Коэффициенты
к 7о 71 72 7э 74 75 Общий множитель
1 1 - - - - - 1
2 3 -1 - - - - 1/2
3 23 -16 5 - - - 1/12
4 55 -59 37 -9 - - 1/24
5 1901 -2774 2616 -1274 251 - 1/720
6 4277 -7923 9982 -7298 -2877 -475 1/1440
Приведем коэффициенты ^г+1 ,[ для к = 0,1, 2, 3. (см. например [6-8])
^г+1 ,1
1 1 1 1 1 ...
1
11
11
11
\
^г+1 ,1
.
^г+1 ,1 =
12
/ 9 9 9 9 9 ...
0
5
5
5
5
/4
3
3
3
3
...
27
11 23
16 7
16 12 16 12
23
7
12
23
7
23
\
1
3
1
шг+1,1
1
24
/64 2 СО - 64
55 5 5 55
55 -4 42 -4 55
55 -4 33 33 -4 55
55 -4 33 24 24 -4
55 -4 33 24 24 33
55
V
■)
Отметим, что для нечетных к коэффициент Шг+1,0 = 0, поэтому нам не потребуется начальное значение Ж0 (см. [6-8]) и в этих таблицах при нечетных к первый нулевой столбец опущен.
Отметим, что из самих формул приближенного вычисления интеграла (5) следует рекуррентное соотношение:
шг+1,7 = шг,7,.] = 0 — к — 1,
шг+1,г—к = шг,г—к + Тк) шг+1,г—к— 1 = шг,г—к— 1 + Ок— 1) •••) шг+1,г = ш0^
(6)
В работах [6, 8] были построены многошаговые методы, которые основаны на формуле (5) для численного решения интегрального уравнения
I
J К(£, в)ж(в)^в = /(£), £ € [0,1], 0 < в < £,
(7)
с достаточно гладкими ядром К(£, в) и /(£) с условиями
К(М) =0 V* € [0,1], /(0) = 0-Данные к-шаговые методы имеют вид
г
Л^шг+МКг+1,гжг = /г+1,г = к, к + 1, •••,Ж - 1, 1=0
(8)
где Кг+1,1 — К (^г+Ъ ^1), /г+1 — / (^г+1), Ж1 ж(^), Л, = 1/Ж, ш называется весами квадратурной формулы.
В статье [10] были выделены задачи вида:
где
А(£)ж (£) + В(^)ж(^) ^ ^ К(£,в,ж(в))^в = /(£), £ € [0,1],
0
^еМ(£) = 0,
(9)
с начальным условием
ж(0) = ж0,
имеющие единственное решение. Для численного решения таких задач были предложены и обоснованы многошаговые методы, основанные на формулах дифференцирования назад для первых двух слагаемых и на явных квадратурных формулах Адамса для интегрального слагаемого. Эти методы имеет вид:
кг
Аг+1 у: а Жг+1—^ + ЛВг+1Жг+1 + Л2 у: Шг+1,1 Кг+1,гжг = Л/г+1, г = к, к + 1, •••, Ж, (10)
7=0
.7=0
0
где а.% - коэффициенты формулы дифференцирования назад, а шг+1,г - веса квадратурной формулы явного метода Адамса. Для задач, рассмотренных в данной статье, методы (10) принципиально неприменимы в силу того, что матрица одАг+1 + ЬВг+1 может быть тождественно вырожденной (детали приведены в заключительном разделе).
2. Численные методы решения
В данном разделе приведены и обоснованы моногошаговые методы численного решения задачи (1), (2), которые удовлетворяют условиями теоремы 1. Приведем численный метод решения задачи (1), (2) с условием (3). Зададим на отрезке [0, 1] равномерную сетку £г = гЛ, г = 1, 2,..., Ж, Ь = 1/Ж. Обозначим Аг = А(£г),Кг7 = К (£г,£7-),/г = / (£г),жг ^ ж(£г). Тогда общие многошаговые методы имеют вид:
к к г+1
Аг+1 а7жг_7+1 + ЬВг+1 ^ ' в7жг+1_7 + Ь ^ ' шг+11 Кг+11ж1 — Ь/г+1 • (11)
7=0 7=0 г=о
Предполагается, что при реализации данных алгоритмов стартовые значения Ж1,Ж2, ...,Жк-1 вычислены достаточно точно и Хо = ж(0).
Формулы (11) могут быть:
1) явными при а0 = 0, в0 = 0, ш^+1,г+1 = 0;
2) полуявными при а0 = 0, в0 = 0, ш^+1,г+1 = 0, или при а0 = 0, в0 = 0, ш^+1,г+1 = 0;
3) неявными при а0 = 0, в = 0, ш^+1,г+1 = 0.
В силу вырожденности матрицы А(£) явные методы для рассматриваемых задач применять нельзя. Полуявные методы, также неприменимые для задачи (1), (2), удовлетворяющих условиями теоремы 1 в силу того, что матрица («0^+1 + Ьв0^г+1) или (одАг+1 + Ь2Шг+1;г+1Кг+1;г+1), могут быть вырожденными. Напомним (см. теорему 1 и лемму 2), что рассматриваемый класс задач включает в себя и интегральные уравнения Вольтера первого рода с ядром на диагонале не равным нулю. А для таких уравнений многие неявные многошаговые методы порождают неустойчивый процесс.
Для численного решения рассматриваемых задач мы предлагаем модифицированные методы, основанные на явных формулах Адамса. А именно, для вычисления слагаемого в уравнении (1) будем использовать к—1-шаговый явный метод Адамса, описанный в предыдущем параграфе. Вырождение А(£)ж (£) + В(£)ж(£) в точке £ = £г+1 будем находить по экстраполяционным формулам. Будем вычислять жг+1 как значение интерполяционного полинома степени к — 1, проходящего через (жг,£г), (жг_1,£г_1),..., (ж^-к+1, £^—к+1) в точке £ = £г+1, т.е.
к-1
жг+1 — ^к_1 жг_1, ..., жг_к+1,4) — ^ ' в7жг_7.
7=0
Аналогично, жг+1 ~ ж' (£)|*=*4+1 - значение производной интерполяционного полинома степени к, проходящего через точки (жг,£г), (жг_1,£г_1), —, (жг_к, £г_к) в точке £ = £г+1, т.е.
' ' 1 к
жг+1 = ^к(жг, жг_ 1, ..., жг_к) = Ь ^ у а7жг_7.
7=0
Таким образом предложенные многошаговые методы имеют вид:
к к_1 г
Аг+1 а7жг_7 + ЬВг+1 ^ ' в7жг_7 + Ь ^ ' шг+11Кг+11ж1 — Ь/г+1. (12)
7=0 7=0 г=0
Приведем коэффициенты в7' для различных к = 0,1, 2, 3, 4, 5 в табл. 2.
Таблица 2
Коэффициенты в7'
к в0 в1 в2 вз в4 в5
0 1 - - - - -
1 2 -1 - - - -
2 3 -3 1 - - -
3 4 -6 4 -1 - -
4 5 -10 10 -5 1 -
5 6 -15 20 -15 6 -1
Приведем коэффициенты а7- для различных к = 1, 2, 3, 4, 5, 6 в табл. 3.
Таблица 3
Коэффициенты а7-
к «0 «1 «2 аз а4 а5 «6 Общий множитель
1 1 -1 - - - - - 1
2 5 -8 3 - - - - 1/2
3 26 -57 42 -11 - - - 1/6
4 127 -414 534 -322 75 - - 1/12
5 522 -1755 2540 -1980 810 -137 - 1/60
6 669 -2637 4745 -4920 3015 -1019 147 1/60
Непосредственные вычисления показывают, что корни характеристических уравнений
к
5>7 Рк"7 =0 (13)
7=0
при к < 6 лежат в единичном круге, и на границе круга нет кратных корней, а при к > 7 есть хотя бы один корень, по модулю больший единицы.
Аналогично, для характеристических уравнений [4]
к
5>/"7 =0 (14)
7=0
и для уравнений (6) - (8)
к
5>7 рк_7 =0 (15)
7=0
при к < 5 корни лежат в единичном круге, и на границе круга нет кратных корней.
Теорема 2. Пусть для задачи (1), (3) выполнены условия теоремы 1, тогда справедлива
оценка ||ж7 — ж(£7-)|| = 0(Ьк) У? = 1, N где ж7- определены из системы (12), ||жг — ж(£)|| <
КгЬк, Кг < то, 1 = 0,1,..., к — 1 .
Доказательство. В силу леммы 2 и второго условия теоремы, существует невырожденная матрица P(t) с элементами из Cтакая, что
/A+iN
Pi+iAi+i = о
V 0
где rankAj+i = r = const. Здесь и в дальнейшем принято обозначение Pi+i = P(ti+i)-Умножая (12) на матрицу Pi+i, получим
^Ai+i\ i k /Bi+i\ k i /Ki+i,A /gi+i\
0 I h Eaxi-j + I Bi+i I S exi-j + ^I]wi+i,i | K3+v) xi = I g+i I , (16)
0 / j=o V 0 / j=° i=o \KiVu/ \giVi/
где
/a1 A2 A3\ /Ai+i
Pi+iAi+i = 0 0 0=0
0 0 0 0
/B1 B2 B 3\ /B3+
Pi+iBi+i = B4 B5 B6 = B2+i
0 0 0 0
/Ki K2 K 3\ \ (g+i-
Pi+iKi+i,i = к4 к5 к6 = Ki+M , Pi+ifi+i = I g2+i
\K7 k8 kV \Kf+M/ 3
Положим £i = Xi — x(ti), тогда уравнения ошибки имеет вид:
к кг
А+1 X] а7£г_7 + ЬВг+1 ^ в7£г_7 + Ь2 ^ Шг+1,гКг+1,г= Мг+Ъ (17)
7=0 7=0 г=0
кг
ЬВг+1 ^ в7£г_7 + шг+1,гКг2+1,г= Ь^г+Ъ (18)
7=0 г=0
г
ь2 ^ шг+1>г Кг+1,г £г = ЬРг+1. (19)
г=0
Вычитая из (г + 1)-го уравнения (19) г-е, умножая полученное на Ь_2, обозначая
АКг+1,г = Кг+1,г — Арг+1 = рг+1 — рг
и учитывая формулу (8), получим
г_к_1
шг+1,г АКг3+1;г£г + ((шг,г_к + 7к )Кг3+1,г_к — шг,г_к Кг3г_к )ег_к +
г=0
((шг,г_к+1 + 7к_1)Кг3+1,г_к+1 — шг,г_ к+1 Кг3,г_ к+1) ^г_ к+1 + ... + 70Кг3+1,гег = Ь 1АРг+Ъ
г = к, к + 1,..., N — 1.
Эти равенства перепишем в виде
ХЛ?К3+М_,-Єг-? + ^ ^г+1,1 ДКІ+11|ЄІ = Л 1Арг+1, І = к, к + 1, ..., N - 1.
(20)
І=0
1=0
По первому условию теоремы (достаточная гладкость входных данных) справедливо представление
АКг+1,г = ||мт,гI < м3 < г = к,к + 1,...,Ж — 1.
АКг+1,г = Ь^М+1,г, ||мг2+1,г! < м2 < г = к,к + 1,...,Ж — 1.
Используя стандартные оценки для приближенного вычисления определенного интеграла [5] и учитывая первое условие теоремы, получим:
||Арг+1|| < Ь2Ьк+1, г = к, к + 1,...,Ж — 1,
|Ь_1Арг+11 < Ь2Ьк, Ь2 < то, г = к, к + 1,..., N — 1.
или
По аналогии имеем:
РтІІ < £3Лк, £з < то, І = Л,Л + 1,...,N — 1. Учитывая (20), перепишем уравнения (17),(18),(19) в виде:
аоА1+1 + Ьв0Вг1+1
«1^1+1 + ^1^+1
вов2
7оК3
+1
Єг +
+1
в1В?
71К3
+1
Єг-1 + ••• +
+1
ак Аі+1 + ^к
+
вк В 2+1
+1
7к К3
+1
Єг-к + Л X] +1,1 1=0
К1+1,1
М+и
мг+1,і
Єї =
(21)
(22)
(23)
где ^ Ь ^Р^+Д г = k,к + 1..., Ж — 1.
В силу второго условия теоремы и леммы 3, при достаточно малых Ь, матрицы
^ао А1+1 + ЛвоВ /+1^
вгВг2+1 | , І = Л, Л + 1,.., N — 1
7оК3
+1,
являются невырожденными, а по первому условию теоремы (достаточная гладкость входных данных) справедливо представление
^аоА1+1 + ЛвоВ1
где
-1
9 Аі+1 + в1
7о*3
г+1
7д К3+1, г -
9
Ег + М1 «0 Г 1
0
0
Е9 + Л^г+1,г-д,
0
0
0 £ Еп-г-1+Мзу
к
0
а ||М+1;г_д|| < ^ < то, г = к, к + 1,..., N — 1.
Умножая уравнения (23) на матрицы
^0^+1 + Ьв0 в1+^ в0Вг+1
70К3+1,г
_1
получим
к
ег + Егег_1 + Ь Мг+1,7е7 =
г=1 7=1
(24)
где
Мг+1,7 — шг+1,7
^0А1+1 + Ьв0Вг1+1\ 1 / К1
в0В2+1 70К?
г+1,г
г+1 >7
Мг2+1,7
^0А!+1 + Ьв0В ^ Зг = | в0Вг2+1
70Кг3+1,г
-1
рг, г = к, к + 1,..., N.
Обозначая 0г = (еТ, е^,..., е^к)Т в формулах (24) стандартным образом, совершим переход от многошаговых методов к одношаговым. Будем иметь
г _ 1
(25)
г=1
где
(—Е —Е2 ■ ■ ■ — Ек^
Е 0 .. .0
С= 0 Е .. .0
0 0 .. .Е
нормы матриц ^7 равномерно ограничены, а векторы Фг определены по правилу
Фг = (ЗТ , 0, 0,..., 0)т .
Собственные числа матрицы С совпадают с корнями характеристических уравнений (13), (14), (15), которые, как было отмечено выше, при к < 6 по модулю меньше единицы. Таким образом, в силу условия теоремы
||00|| < р4Ьк, Ь4 < то,
а в силу оценок (21) и (22), имеем
||Фг|| < р5Ьк, Ь5 < то, г = к, к + 1,..., N.
Вспоминая, что ^ = (еТ, £?_!, ..., £?_к)Т из приведенных выше оценок и леммы 1 следует,
что
|ег|| < Ь6Ьк, Ь6 < то, г = к, к + 1,..., N — 1.
□
3. Численные эксперименты
Приведем очень простой модельный пример:
'1 0 0\ /1 0 1\ ‘ / е
0 0 0 ж' (*) + 0 1 0 ж(*) + /
0 0 0 0 0 0
000
,*+« 0
0 е*-5
0 0
0
0
е*+2«
'е-2* + *е* ж(т)^т = | (1 + *)е* *е*
где ж-трехмерная искомая вектор-функция с начальным условием ж(0) = (1,1,1)т. Данный пример удовлетворяет всем условиям теоремы 1, и точное решение ж(*) = (е-*, е*, е-2*)т. Легко заметить, что методы (10) принципиально нельзя применить для данного примера, т.к. ^е£(а:А + ЛВ) = 0.
Усложним данный пример, а именно умножим исходную систему на матрицу
Р (*) =
1 0 0>
е* 10
2* е* 1,
е
и произведем замену переменной ж(*) = ф(*)у(*), где
/1 2* ^ £(*) = 0 1 3*
0 0 1
В итоге получим
*
А(%'(*) + В(*)у(*) + ^ К(*, в)^в = / (*), * € [0,1],
(26)
где А(*) =
К(*,8) =
1 2*
е* 2£е*
^е2* 2е2*
/ ' е*+5
= е*е*+5
\ е2*е*+«
е**2
е2**2
1 2* + 2 *2 + 1 + 2*
, В(*) = | е* 2*е* + 1 + 2е* е**2 + 3* + е* + 2*е*
2* 2е2** + е* + 2е2* е2**2 + 3*е* + е2* + 2е2*
2е*+^ е*+5*2
+ е*-5 е*е*+^2 + 3е*-^
е 2* + *е*
/ (*)=| е*(е-2* + *е*) + (1 + *)е*
че2* (е-2* + *е*) + (1 + *)(е*)2 + *е*у
Точное решение у(*) = ^-1 (*)ж(*), начальное условие у(0) = (1,1,1)т.
Численные расчеты данной задачи при различных значениях к приведены в табл. 4.
Таблица 4
Численные расчеты задачи при различных значениях к
о
о
2
с
егг к=1 к=2 к=3
N=5 N=10 N=20 N=40 N=80 1,309600415814891 0,7497289570481798 0,3988507964835724 0,2051764163549656 0,1039752161311108 0,6015407275019990 0,1844243516458794 0,0503707677718254 0,0129986398315527 0,0032742356352037 0,21171281782986052430 0,04761740960151257878 0,00732509005266374868 0,00097017989140169301 0,00012382133627371258
ЗДесь err = max д/(y1 - y1(ti))2 + (У2 - y2(ti))2 + (У3 - У3(ti))2-
В качестве стартовых значений для у 1, У2, Ук— 1 были взяты точные значения
у(t 1), у(^2), •••, y(tfc-i)- Из приведенной табл. 4 видно, что расчеты хорошо согласуются с теоритическими выкладками.
Исследования были поддержаны грантом РФФИ 13-01-93002Вьет-а.
Литература
1. Ушаков, Е.И. Статическая устойчивость электрических систем / Е.И. Ушаков. - Новосибирск: Наука. сиб. отд-ние, 1988.
2. Сенди, К. Современные методы анализа электрических цепей / К. Сенди. - М.: Энергия, 1971. - 360 с.
3. Булатов, М.В., Об одном семействе вырожденных интегро-дифференциальных уравнений / М.В. Булатов, Е.В. Чистякова // Журнал вычислительной математики и математической физики. - 2011. - Т. 51, №9. - С. 1665-1673.
4. Булатов, М.В. Численное решение интегро-алгебраических уравнений многошаговыми методами / М.В. Булатов, О.С. Будникова // Журнал вычислительной математики и математической физики. - 2012. - Т. 52, № 5. - С. 829-839.
5. Бахвалов, Н.С. Численные методы / Н.С. Бахвалов. - М.: Наука, 1975.
6. Тен Мен Ян Приближенное решение линейных интегральных уравнений Вольтерра I рода: дис.... канд. физ. мат. наук / Тен Мен Ян. - Иркутск, 1985.
7. Brunner, H. The Numerical Solution of Volterra Equations / H. Brunner, P.J. van der Houwen. - Amsterdam: North-Holland, CWI Monographs 3, 1986.
8. Linz, P. Analytical and Numerical Methods for Volterra Equations / P. Linz. - SIAM, Philadelphia, 1985.
9. Булатов, М.В. Об интегродифференциальных системах с вырожденной матрицей перед производной / М.В. Булатов. - Дифференциальные уравнения. - 2002. - Т. 38, №5. -С. 692-695.
10. Булатов, М.В. Численное решение интегро-дифференциальных систем с вырожденной матрицей перед производной многошаговыми методами / М.В. Булатов, Е.В. Чистякова // Дифференциальные уравнения. - 2006. - Т. 42, № 9. - С. 1248-1255.
11. Булатов, М.В. Применение матричных полиномов к исследованию линейных дифференциально алгебраических уравнений высокого порядка / М.В. Булатов, Ли М.Г. // Дифференциальные уравнения. - 2008. - Т. 44, № 10. - С. 1299—1305.
Михаил Валерьянович Булатов, доктор физико-математических наук, профессор, Институт динамики систем и теории управления СО РАН (г. Иркутск, Российская Федерация), mvbul@icc.ru.
До Тиен Тхань, аспирант, Национальный исследовательский Иркутский государственный технический университет (г. Иркутск, Российская Федерация), thanhdotien278@yahoo.com.
Поступила в редакцию 16 мая 2014 г.
Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming & Computer Software",
2014, vol. 7, no. 3, pp. 93—106.
MSC 65R20 DOI: 10.14529/mmp140310
Multistep Method for Solving Degenerate Integral-Differential Equations
M. V. Bulatov, Institute for System Dynamics and Control Theory of Siberian Branch
of Russian Academy of Sciences, Irkutsk State Technical University, Irkutsk, Russian
Federation, mvbul@icc.ru,
Thanh Do Tien, National Research Irkutsk State Technical University, Irkutsk, Russian
Federation, thanhdotien278@yahoo.com
In this work we consider the linear integral-differential equations of the fist order with the identically singular matrix at the derivative. For these systems, the initial conditions is given and assumed consistent with the right part. Considered tasking in this paper arise in the mathematical modeling of complex electric circuits. By using the apparatus of matrix polynomials a class of problems, which having a unique solution, is marked out. The difficulties of the numerical solution of such problems, in particular the instability of many implicit methods is considered. For numerical solution of this class of problems we have suggested multistep methods, which are based on an explicit quadrature formula for the integral term Adams and extrapolation formulas. Sufficient conditions for the convergence of such algorithms to the exact solution is formulated.
Keywords: integral-differential equations; multistep methods; matrix polynomials.
References
1. Usakov E.I. Staticheskaya ustoychivost’ elektricheskikh sistem [Static Stability of Electrical Systems]. Novosibirsk, Nauka, 1988.
2. Scendy K. Sovremennyye metody analiza elektricheskikh tsepey [Modern Methods of Analysis of Electric Circuits]. Moscow, Energiya, 1971. 360 p.
3. Bulatov M.V., Chistyakova E.V. On a Family of Singylar Integro-Differential Equations. Computational Mathematics and Mathematical Physics, 2011, vol. 51, issue 9, pp. 1558-1566. DOI: 10.1134/S0965542511090065
4. Budnikova O.S., Bulatov M.V. Numerical Solution of Integral-algebraic Equations for Multistep Methods. Computational Mathematics and Mathematical Physics, 2012, vol. 52, issue 5, pp. 691-701. DOI: 10.1134/S0965542512050041
5. Bakhvalov N.S. Chislennyye metody [Numerical Methods]. Moscow, Nauka, 1975.
6. Ten Men Yan Priblizhennoye resheniye lineynykh integral’nykh uravneniy Vol’terra I roda [Approximate Solution of Linear Volterra Integral Equations of the First Kind. Candidate’s Dissertation in Mathematics and Physics]. Irkutsk, 1985.
7. Brunner H., van der Houwen P.J. The Numerical Solution of Volterra Equations. Amsterdam: North-Holland, CWI Monographs 3, 1986.
8. Linz P. Analytical and Numerical Methods for Volterra Equations. SIAM, Philadelphia, 1985. DOI: 10.1137/1.9781611970852
9. Bulatov M.V. Regularization of Singular Systems of Volterra Integral Equations. Computational Mathematics and Mathematical Physics, 2002, vol. 42, issue 3, pp. 315-320.
10. Bulatov M.V., Chistyakova E.V. Numerical Solution of Integro-Differential Systems with a Degenerate Matrix Multiplying the Derivative by Multistep Methods. Differential Equations, 2006, vol. 42, issue 9, pp. 1317-1325. DOI: 10.1134/S0012266106090102
11. Bulatov M.V., Lee M.-G. Applications of Matrix Polynomials to the Analysis of Linear Differential-Algegraic Equations of Higher Order. Differential Equations, 2008, vol. 44, issue
10, pp. 1353-1360. DOI: 10.1134/S0012266108100017
Received May 16, 2014