CC BY
18
ВЫЧИСЛИТЕЛЬНАЯ МАТЕМАТИКА
УДК 519.624.2
DOI: 10.18101/2304-5728-2019-3-32-41
О ВЫБОРЕ КРАЕВЫХ УСЛОВИЙ ДЛЯ ДИФФЕРЕНЦИАЛЬНО-АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА1
© Ботороева Мария Николаевна
старший преподаватель,
Иркутский государственный университет, Педагогический институт Россия, 664011, г. Иркутск, ул. Нижняя Набережная, 6
научный сотрудник,
Бурятский государственный университет имени Доржи Банзарова Россия, 670000, г. Улан-Удэ, ул. Смолина, 24а E-mail: [email protected]
© Будникова Ольга Сергеевна
кандидат физико-математических наук, доцент, Иркутский государственный университет, Педагогический институт Россия, 664011, г. Иркутск, ул. Нижняя Набережная, 6 E-mail: [email protected]
© Соловарова Любовь Степановна
кандидат физико-математических наук,
Институт динамики систем и теории управления
имени В. М. Матросова СО РАН
Россия, 664033, г. Иркутск, ул. Лермонтова, 134, а/я 292
E-mail: [email protected]
В статье рассмотрены системы линейных обыкновенных дифференциальных уравнений второго порядка на конечном отрезке интегрирования с тождественно вырожденной квадратной матрицей перед второй производной — так называемые дифференциально-алгебраические уравнения второго порядка. Такие постановки задач достаточно часто возникают в приложениях. Отмечены сложности качественного исследования и построения численных методов решения рассматриваемых уравнений. Предполагается, что для систем, исследуемых в данной работе, заданы краевые условия, число которых меньше, чем удвоенная размерность исходной системы. Выделен класс дифференциально-алгебраических уравнений второго порядка, для которого, используя известные результаты теории проекторов, предложен алгоритм выбора недостающих краевых условий. Для иллюстрации разрабатываемого авторами метода рассмотрен простой пример. В заключении статьи уделено внимание перспективам дальнейшего исследования данного алгоритма.
1 Ботороева М. Н. поддержана грантом РФФИ, № 18-51-54001_Вьет_а, и Министерством образования и науки РФ, проект 1.5049.2017/БЧ; Будникова О. С. поддержана грантом РФФИ, № 18-51-54001_Вьет_а; Соловарова Л. С. поддержана грантом РФФИ, № 18-29-10019 мк.
Ключевые слова: краевая задача; краевые условия; дифференциально-алгебраические уравнения; вырожденные обыкновенные дифференциальные уравнения второго порядка; метод матричной прогонки; начальная задача; ал-гебро-дифференциальные уравнения; численные методы; дифференциальные уравнения второго порядка; трехточечная аппроксимация.
Для цитирования:
Ботороева М. Н., Будникова О. С., Соловарова Л. С. О выборе краевых условий для дифференциально-алгебраических уравнений второго порядка // Вестник Бурятского государственного университета. Математика, информатика. 2019. № 3. С. 32-41.
Введение
Ряд природных и технических явлений можно описать взаимосвязанными обыкновенными дифференциальными уравнениями (ОДУ) первого и второго порядков и конечномерными (трансцендентными) уравнениями [8; 10; 14]. Объединяя эти уравнения в виде системы, мы получим систему ОДУ второго порядка с тождественно вырожденной матрицей перед второй производной. Такие системы принято называть дифференциально-алгебраическими уравнениями (ДАУ) второго порядка.
В настоящее время бурно развивается качественная теория и численные методы решения ДАУ первого порядка. Это связано с тем, что такие уравнения имеют принципиальные отличия от систем ОДУ и, кроме того, имеют многочисленные приложения [8-10; 12; 14].
С ДАУ порядка выше первого, как правило, поступают следующим образом: вводят новую вектор-функцию, размерность которой равна kn ( n — размерность исходной задачи, k — порядок ДАУ), и переписывают исходную задачу в виде ДАУ первого порядка [13].
Описанный подход значительно усложняет как качественное исследование, так и создание численных алгоритмов для полученной задачи. Этот факт и инициировал исследования, приведенные ниже, которые являются продолжением работ [3; 4; 7; 11].
1 Постановка задачи
Рассмотрим систему
A(t)x"(t) + B(t)x'(t) + C(t)x(t) = f (t), t î [0,1] (1)
где A(t), B(t) и C(t) — (n x n) матрицы, f (t) — заданная, а x(t) — неизвестная n -мерные вектор-функции.
В дальнейшем изложении предполагается, что
det A(t) ° 0. (2)
Системы (1) с условием (2) будем называть дифференциально-алгебраическими уравнениями второго порядка (ДАУ 2).
Предполагается, что для системы (1) заданы краевые условия
(3)
Г«о х(0)+& х'(0 )=Уо [«1 х(1) + Д х'(1) = 7 где а0, Ь0 — (т0 х и) -матрицы, Д — (т1 хи)-матрицы. В частности, начальную задачу для (1) можно переписать в виде (3) с матрицами
( Е„ Л ( 0 Л
0
Е
п
, /30 = , а1 и Ь — нулевые матрицы, а Еп — единичная
матрица, 0 — нулевая матрица размерности (и х п).
В отличие от ОДУ второго порядка для выделения единственного решения ДАУ2 начальные или краевые условия можно не задавать. Например, ДАУ2 вида
(а, (() 0Л( х1 ()Л" (Ь1 (() 0Vх1 (0 1 Л(х1 ()Л ( 0 Л
0 0
V 2
0 0
^ 2У >0
1 0
Х2 ()0
[0,1] (4)
при достаточно гладких функциях а1 (), Ь), ) имеет единственное решение: х1 ^) = ), х2 ^) = -а1 ())- Ь1 () ).
Краевые условия (3) должны также быть заданы корректно, т. е. быть согласованными с правой частью (1).
Приведем критерий корректности задания краевых условий (3). Нам понадобится известное определение.
Определение 1. (см, например, [2; 5]). Матрица называется полуобратной к матрице Л, если она удовлетворяет матричному уравнению
ЛЛ~ Л = Л.
Из этого определения следует, что
УЛ = (Е - ЛЛ~)Л = 0. (5)
Умножаем систему (1) на матрицу У в точках t = 0 и t = 1. В силу (5) получим:
(6)
|У ( 0) С ( 0) х ( 0) + У ( 0) В ( 0 ) х ' ( 0 ) = У (0 ) / ( 0 ),
IУ (1) С (1) х (1) + у (1) в (1) х'( 1) = У (1) / (1).
Вводя обозначения у0 = (хт (0), х'Т (0)) у = (хт (1), х'т (1)), У(0)С(0) = Ь0 У(0)В(0) = М0, У(0)/(0) = N0, У(1)С(1) = Ц, У(1)В(1) = М„ У(1)/(1) = N1, краевые условия (3) и условия (6) перепишем в виде:
Г(а01 Ь0 К = 70 Ц«! А Ь =7:'
Г(^0 | М0 )У0 = Nо (8)
[(/,1 М Ь = N1 . ()
Здесь и всюду в дальнейшем изложении обозначение (а|ь) означает расширенную матрицу. Таким образом, если задача (1), (3) имеет решение
(7)
(возможно, не единственное), то множество решений системы линейных алгебраических уравнений (СЛАУ) (7) является и множеством решений СЛАУ (8). Обратное неверно.
Для доказательства этого факта достаточно рассмотреть обыкновенное линейное дифференциальное уравнение
x"(t) + B(t )x'(t) + С (t ) = f (t),
где B(t) и С(() — функции. В этом случае = 1, V(о) = V(l) = L0 = = M0 = L = M1 = N0 = N1 = 0, то есть решением системы (8) является любой вектор из R2.
2 Методы решения краевой задачи
Опишем алгоритм численного решения краевой задачи (1), (3). Стандартным образом введем на отрезке интегрирования сетку
ti = ih, i = 0,1,...,N, h = 1 и обозначим Ai = A(ti), Bi = B(ti), Ci = C(ti),
f = f(ti \ x » x(ti).
Классическая трехточечная дискретизация [1; 6] системы (1) дает
Di^ + E x + F,.xi_1 = Gi, i = 1,2,...,N -1, (9)
где (n x n) -матрицы D. , E i и F. зависят от формулы, применяемой для аппроксимации (1).
Краевые условия (3) после двухточечной аппроксимации примут вид: iDo x + E о xo = Go, (10a)
[E nxN + FNXN-1 = Gn , (106)
где D0 и E0 — (m0 x n) -матрицы, EN и FN — (m1 x n)-матрицы.
Условий (10a), (106) недостаточно для того, чтобы применить метод матричной прогонки для решения СЛАУ (9), (10a), (106) , так как эта система является недоопределенной (содержит (m0 + m1 + n ■ (N - 2)) уравнений, а неизвестных n ■ N). Поэтому мы предлагаем поступить следую -щим образом. Умножим (1) на матрицу V(t) = E — A(t)A~(t) и вычислим полученное выражение в точках t = 0 и t = 1. В силу формулы (4) имеем:
V (0)С (0)x(0) + V (0)B(0)x'(0) = V (0)f (0) и
V (1)c(1)x(1) + V (1)B(1)x'(1) = V (1)f (1). Аппроксимируя эти условия, например простейшей формулой численного дифференцирования, получим:
S0 x1 + T0 x0 = H o, (11a)
SNxN + '7NxN-1 = HN , (11б)
где ^о = V (о)[в(0) + НС (о)], Т0 = -V (о)в(о), И0 = V (о)/ (о),
^ = V(1)[в(1)+нс(1)], гм = -V(1)5(1), Им = V(1)/(1).
С краевыми условиями (Юа), (Ю6) поступим следующим образом: умножим системы
N Е о;
Л
^ 0
а
(Е
NI ^
о,
N -1 0
на матрицы Р о
Еп и Е,
FN соответственно. Получим:
РоТ Ро X + РоТ Е о X
ОоГ Оо,
о 1 "о 1-оЛо
Е о РоХ1 + Е о Е охо = Е о Оо,
Е N Е NXN + Е N FNXN-1 = Е N GN , FN ЕNXN + FN FNXN-1 = FN GN .
С краевыми условиями (11а), (116) поступим аналогично: умножим системы
(12)
(13)
(14)
(15)
(?оТ
Л
= И
Но 0
Т N
N
^N-1 0
= И
N
на матрицы 8о
ТоТ и , ТN . Получим:
8о 8оХ1 + 8о Тохо = 8о Ио, То 8о Х1 + То То хо = То И о,
SN SNXN + SN Т NXN-1 = И
N■>
Т N SNXN + Т N TNXN-1 = Т N ИN ■
(16)
(17)
(18) (19)
Итогом суммирования уравнений (12), (13), (16), (17) и (14), (15), (18), (19) является
|Ко X1 + Ро Xо = Чо, (2о)
[KNXN + PNXN-1 = qN , (21)
KN и — квадратные (п х п) -матрицы, определенные по
где Kо, Ро правилу
Ко = РоТ Ро + Е оТ Ро + 8оТ 8о + ТоТ ¿о, Ро = РоТ Е о + Е оТ Е о + 8оТ То + ТоТ То,
X
и
X
N
т
т
X
и
т
К = Р т Р + Р т Р + К т К + Т Т К Р = Р ТР + Р ТР + К Тт + т Тт
1N~^-N ^N^^N N N 1N^1N 1N>
40 = °0Т с0 + Р 0Т с0 + К 0ТЯ 0 + TоTно,
qN = Р N ^ + РN ^ + КN НN + ^ НN .
Для иллюстрации предлагаемого метода выбора недостающих краевых условий рассмотрим простой пример. Пример 1. Система
Г1 0 Л( и'^ )Л (1 0 Л( и(( )Л Г 0 ^
0 0
/'(()
Краевые условия имеют вид
/ V
u
0 1
V ^ 0
t е [0,l].
/(0) + u '(0 ) = 1,
[u (l) + u '(l) = sin 1 + cos 1.
Точное решение данной системы u (t) = sin t, v(t) = et.
Применим описанный алгоритм поиска недостающих краевых усло-
Г1 0| .,
(псевдообратная матрица), имеем
вий. Полагая A =
00
= A ,
V (t) =
Г0 0|
V0 10
. Умножим исходную систему слева на V (t) и вычислим
полученное выражение в точках t = 0 и t = 1. В силу (5) получим:
f 0 0 Y u (t )| f 0 1 v(()
0 1
Из второго уравнения видно, что
V e 0
, t е[0,1].
,(0)= 1,
[у(1) = е.
Обозначим x(t) = (и^), v(t ))т и с учетом введенных ранее обозначений система (7) примет вид:
0 0 \ 0 00
ч
а системе (8) соответствует
1 0
1 0 00 10 00
Г11
У0 = V 0 0
У1 =
Г sin 1 + cos 1|' 0
Т 0 0 0 01 Г 01
0 V 1 0 0 / У0 = V 11
f 0 0 0 0 Л Г 01
0 Л 1 0 0 У У1 = V e 0
Далее применим двухточечную дискретизацию для исходных и полученных краевых условий и выполним дальнейшие указания описанного алгоритма. В результате, опустив все громоздкие выкладки, получим искомые краевые условия в виде:
(1 о^ (н -1 о^ ( и\
о
(Н +1
о
о
0 ^
1
x1 +
о 1
(1
о
о ^ о
V 10
XN -1 =
Н^т 1 + 0081)Л '
Заключение
В статье не рассматривается конкретный вид трехточечной аппроксимации (9). Отметим, что ряд известных формул численного дифференцирования нельзя применять для (1) даже с корректно заданными (согласованными с правой частью) краевыми условиями
| ^ ^ (22)
и в предположении, что решение задачи (1), (22) существует, единственно и достаточно гладкое. Это происходит, в частности, потому, что соотношения (9) и условия (22) дают вырожденную СЛАУ или процесс матричной прогонки является неустойчивым.
Также отметим, что предлагаемый подход не всегда дает СЛАУ (2о), (21) полного ранга. Например, если применить предлагаемый подход для ДАУ (4), имеющей единственное решение, не зависящее от условий (3), и полагая а1(о) Ф о, а1(1) Ф о, то краевые условия (19), (2о) имеют вид:
Г^о) = р(о)
I х(1) = ?(1).
Условия (3) для данного примера не заданы, а условий (23) недостаточно, чтобы формально реализовать метод матричной прогонки для данного примера (после трехточечной аппроксимации (9)).
В дальнейшем планируется выделить классы ДАУ (1) с краевыми условиями (3), для которых предлагаемый алгоритм выбора недостающих краевых условий и подходящая трехточечная аппроксимация порождали бы устойчивый процесс матричной прогонки.
(23)
Литература
1. Бахвалов Н. С. Численные методы (анализ, алгебра, обыкновенные дифференциальные уравнения). М.: Наука, 1975. 632 с.
2. Бояринцев Ю. Е. Регулярные и сингулярные системы линейных обыкновенных дифференциальных уравнений. Новосибирск: Наука, 198о. 222 с.
3. Булатов М. В., Рахвалов Н. П., Р1шоп§ Т. Б. Численное решение краевой задачи для линейных дифференциально-алгебраических уравнений второго порядка // Журнал Средневолжского математического общества. 2о1о. Т. 12, № 1. С. 52-58.
Xо _
XN
е
4. Булатов М. В., Рахвалов Н. П., Фыонг Та Зуй. О методе матричной прогонки для одного класса дифференциально-алгебраических уравнений второго порядка // Известия Иркутского государственного университета. Математика. 2011. Т. 4, № 4. С. 2-11.
5. Ваарман О. Обобщенные обратные отображения. Таллинн: Валгус, 1988. 120 с.
6. Самарский А. А. Теория разностных схем. М.: Наука, 1989. 616 с.
7. Соловарова Л. С. О выборе начальных условий для дифференциально-алгебраических уравнений // Вестник Бурятского государственного университета. Математика, информатика. 2017. № 1. С. 18-22. DOI: 10.18101/2304-57282017-1-18-22.
8. Хайрер Э., Ваннер Г. Решение обыкновенных дифференциальных уравнений. Жесткие и дифференциально-алгебраические задачи: пер. с англ. М.: Мир, 1999. 685 с.
9. Чистяков В. Ф. Алгебро-дифференциальные операторы с конечномерным ядром. Новосибирск: Наука, 1996. 280 с.
10. Bock H. G., Schloder J. P., Schulz V. H. Differential-algebraic equations and their connections to optimization. Heidelberg: Interdisziplinares Zentrum fur wiss. Rechnen der Univ. Heidelberg, 1996. 188 p.
11. Bulatov M. V., Lee M. G. Application of matrix polynomials to the analysis of linear differential algebraic equations of higher order // Differential Equations. 2008. V. 44. P. 1353-1360.
12. Kunkel P., Mehrmann V. Differential-Algebraic Equations. Analysis and Numerical Solution. Zurich, Switzerland: EMS Publishing House, 2006. 192 p.
13. Mehrmann V., Shi C. Transformation of high order linear differential-algebraic systems to first order // Numerical Algorithms. 2006. V. 42, Issue 3-4. P. 281-307.
14. Lamour R., März R., Tischendorf C. Differential-Algebraic Equations: A Projector Based Analysis. Berlin: Springer, 2013. 649 p.
SAMPLING BOUNDARY CONDITIONS FOR SECOND ORDER DIFFERENTIAL-ALGEBRAIC EQUATIONS
Mariya N. Botoroeva Senior Lecturer, Irkutsk State University
6 Nizhnyaya Naberezhnaya, Irkutsk 664011, Russia Researcher,
Dorzhi Banzarov Buryat State University 24a Smolina St., Ulan-Ude 670000, Russia E-mail: [email protected]
Olga S. Budnikova
Cand. Sci. (Phys. and Math.), A/Prof., Irkutsk State University
6 Nizhnyaya Naberezhnaya, Irkutsk 664011, Russia E-mail: [email protected]
Lubov S. Solovarova Cand. Sci. (Phys. and Math.),
Matrosov Institute for System Dynamics and Control Theory SB RAS 134 Lermontova St., Irkutsk 664033, Russia E-mail: [email protected]
The article considers systems of second-order linear ordinary differential equations with an identically degenerated square matrix preceding the second derivative (second order differential-algebraic equations). Such problem statements often arise in applications. We have noted the difficulties in qualitative study and development of numerical methods for solving the equations under consideration. It is assumed that boundary conditions are given and their number is less than the doubled dimension of the initial system. We have defined a class of second order differential-algebraic equations, and using the known results of projection theory proposed an algorithm for sampling the missing boundary conditions. The method developed is illustrated by a simple example. At the conclusion of the article we have paid attention to the further research of this algorithm.
Keywords: boundary value problem; boundary conditions; differential-algebraic equations; ordinary second order degenerate differential equations; tridiagonal matrix algorithm; initial boundary value problem; algebraic-differential equations; numerical methods; second order differential equations; three-point approximation.
References
1. Bakhvalov N. S. Chislennye metody (analiz, algebra, obyknovennye different-sialnye uravneniya) [Numerical Methods (analysis, algebra, ordinary differential equations)]. Moscow: Nauka Publ., 1975. 632 p.
2. Boyarintsev Yu. E. Regulyarnye i singulyarnye sistemy obyknovennykh differ-entsialnykh uravnenii [Regular and Singular Systems of Ordinary Differential Equations]. Novosibirsk: Nauka Publ., 1980. 222 p.
3. Bulatov M. V., Rakhvalov N. P., Phuong T. D. Chislennoe reshenie kraevoi zad-achi dlya lineinykh differentsialno-algebraicheskikh uravnenii vtorogo poryadka [Numerical Solution of Boundary Value Problem for Linear Second Order Differential-Algebraic Equations]. Zhurnal Srednevolzhskogo matematicheskogo obshchestva. 2010. V. 12. No. 1. Pp. 52-58.
4. Bulatov M. V., Rakhvalov N. P., Phuong T. D. O metode matrichnoi progonki dlya odnogo klassa differentsialno-algebraicheskikh uravnenii vtorogo poryadka [On the Tridiagonal Matrix Algorithm for One Class of Second Order Differential-Algebraic Equations]. Bulletin of Irkutsk State University. Series: Mathematics. 2011. V. 4. No. 4. Pp. 2-11.
5. Vaarman O. Obobschennye obratnye otobrazheniya [Generalized Inverse Mappings]. Tallin: Valgus Publ., 1988. 120 p.
6. Samarskii A. A. Teoriya raznostnykh skhem [Theory of Difference Schemes]. Moscow: Nauka Publ., 1989. 616 p.
7. Solovarova L. S. O vybore nachalnykh uslovii dlya differentsialno-algebraicheskikh uravnenii [Sampling the Initial Conditions for Differential-Algebraic Equations]. Vestnik Buryatskogo gosudarstvennogo universiteta. Seriya Matematika, informatika. 2017. No. 1. Pp. 18 - 22. DOI: 10.18101/2304-5728-2017-1-18-22.
8. Hairer E. Wanner G. Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems. Berlin: Springer-Verlag, 1996. 614 p.
9. Chistyakov V. F. Algebro-differentsialnye operatory s konechnomernym yadrom [Differential-Algebraic Operators with a Finite-Dimensional Kernel]. Novosibirsk: Nauka Publ., 1996. 280 p.
10. Bock H. G., Schloder J. P., Schulz V. H. Differential-Algebraic Equations and Their Connections to Optimization. Heidelberg: Interdisziplinares Zentrum fur wiss. Rechnen der Univ. Heidelberg, 1996. 188 p.
11. Bulatov M. V., Lee M. G. Application of Matrix Polynomials to the Analysis of Linear Differential Algebraic Equations of Higher Order // Differential Equations. 2008. V. 44. Pp. 1353-1360.
12. Kunkel P., Mehrmann V. Differential-Algebraic Equations. Analysis and Numerical Solution. Zurich, Switzerland: EMS Publishing House, 2006. 192 p.
13. Mehrmann V., Shi C. Transformation of High Order Linear Differential-Algebraic Systems to First Order. Numerical Algorithms. 2006. V. 42. Iss. 3-4. Pp. 281-307.
14. Lamour R., März R., Tischendorf C. Differential-Algebraic Equations: A Projector Based Analysis. Berlin: Springer, 2013. 649 p.
