КОМПЬЮТЕРНАЯ АЛГЕБРА -
УДК 517.929
СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ С ВЫДЕЛЕННОЙ ЧАСТЬЮ НЕИЗВЕСТНЫХ
© 2015 г. A.A. Панферов
Московский государственный университет имени М.В. Ломоносова 119991 Москва, ГСП-1, Ленинские горы, д. 1 E-mail: ast.a_s@mail.ru Поступила в редакцию 25.09.2014
Пусть дифференциальная система A\y' + Aoy = 0 с возможно вырожденной матрицей Ai имеет полный ранг. Пусть выделена часть неизвестных, т. е. часть компонент вектора у. Предлагается алгоритм, который в этой ситуации позволяет получить для некоторой части у компонент вектора у систему вида у' = By, при этом те выделенные компоненты у, которые не вошли в у, линейно выражаются через вошедшие в у выделенные компоненты.
1. ВВЕДЕНИЕ
Пусть К — дифференциальное поле характеристики 0 с производной д = '. Рассматривается система дифференциальных уравнений
Лгу' + Аоу = 0, (1)
где Аг, Ао е Ктхт, Аг = 0, а у = (уг,..., ут)Т -вектор неизвестных системы. Матрицу А\ будем
Ао
говой матрицей системы.
Пусть система (1) имеет полный ранг, т.е. ее уравнения независимы (последнее будет уточнено в разделе 2). Предположим, часть неизвест-
у
болыний интерес, чем остальные. Назовем эти компоненты выделенными. В этом предположении можно рассматривать ряд задач: поиск только выделенных компонент решений; проверка существования решений, выделенные компоненты которых принадлежат заданным классам; частичная устойчивость решений по выделенным компонентам ([1]) и др.
Если ведущая матрица А\ системы (1) обратима, то можно перейти к нормальной системе
У' = Ау, (2)
где А = —А-гА0.
Для систем вида (2) С.А. Абрамовым и М. Бронштейном в работе [2] был предложен алгоритм, позволяющий для выделенных компонент вектора неизвестных перейти к нормальной системе
г' = Кг,
где компоненты г суть выделенные компоненты у
позволяет решать указанные выше задачи.
Алгоритм Абрамова и Бронштейна был реализован в среде Марье как часть встроенного пакета 0геТоо1Б в виде процедуры 11ес1исес13у^ет.
В случае, когда ведущая матрица А\ исходной системы (1) вырождена, мы имеем дело с дифференциально-алгебраической системой и непосредственно алгоритм Абрамова и Бронштейна применим быть не может.
Ниже предлагается алгоритм, который позволяет в этой ситуации получить для некоторой части у компонент вектора у систему вида
у' = ву, (3)
у
не вошли в у, линейно выражаются через вошедшие в у выделенные компоненты.
К системе (3) далее может быть применен алгоритм Абрамова и Бронштейна, что дает систему, в которую уже не будут входить "лишние" неизвестные.
Таким образом мы получаем возможность решать упомянутые выше задачи не только для систем с невырожденной ведущей матрицей, но и для любых систем вида (1), имеющих полный ранг.
Заметим, что дифференциальные системы более высоких порядков
Аг у(г) + Лг-гу(г-1) + ... + А1У + Лоу = 0 (4) с помощью перехода к вектору неизвестных
Y = (У1,... ,ym,yi,... ,... ,y(r-1),... ,ym-1))T L(y) = 0, где L e K[d]
Пусть Ь е К[д]тхт. Через Ьг,* обозначим г-ю
строку матрицы Ь. Пусть 7 С {1,..., т}, тогда строки Ьг,*, где г е называются К[д] -линейно зависимым,и, если существуют дифференциальные операторы С К[д], не все одновременно равные нулю, такие, что YlгeJ = 0; в противном случае они называются К[д]-линейно независимыми.
Ь
как максимальное число ее К[д]-линейно независимых строк. Дифференциальная система
сводятся к системе первого порядка AiY ' + AOY = 0,
в которой
Ai =
Ao =
0 -I„
Ao Ai
.Ar
0 -I„
... Ar_
r-1
где I„
единичная матрица размера m x m.
2. ПРЕДВАРИТЕЛЬНЫЕ СВЕДЕНИЯ
Обозначим через К [д] кольцо дифференци-
К
Система (1) может быть записана в операторной форме Ь(у) = 0 где Ь = Л1д + Л0 — операторная матрица, т. е. матрица размера т х т с элементами, принадлежащими кольцу К[д].
Произвольная ненулевая операторная матрица Ь е К[д]тхт может быть записана в виде
Ь = Лг дг + Лг-1дг-1 + ■ ■ ■ + Ло,
где г е N Л е Ктхт (г = 0,1,..., г), Лг = 0. Матрица Лг называется ведущей матрицей операторной матрицы Ь. Число г называется порядком Ь и обозначается огё Ь. Для случая Ь = 0 будем считать огё Ь = —те.
Матрица и е К[д]тхт называется унимо-дулярной, если для нее найдется матрица и е
K [d ]m
такая, что UU = UU = Im
, имеет полный ранг тогда и только тогда, когда ранг ¿равен m.
Положим ord Li,* = max1<j<m ord Lj. Вектор-строка S = (¿1,..., Sm), где Si = max(0, ord Li,*),
L
тим, что для операторной матрицы, соответствующей дифференциальной системе (1), Si e {0,1}.
L
матрица -£crow(L) того же размера, что и L, элемент (i,j) которой равен коэффициенту при (i, j)-oro элемента L для всех 1 < < m. Опе-L
строкам, если ненулевые строки ее фронтальной матрицы (L) линейно независимы над K.
Согласно [3], для любой матрицы L e K [d]mxn существует унимодулярная матрица U e K[d]mxm такая, что матрица L = UL является приведенной по строкам. Процесс, послуживший основой доказательства теоремы 2.2 в работе [3], определяет алгоритм, получивший название Row-Reduction. Его формальное описание можно также найти в [4, разд. 2.1].
Входом рассматриваемого ниже алгоритма является дифференциальная система (1) первого порядка, соответствующая операторная матрица которой является приведенной по строкам.
3. АЛГОРИТМ
Основная цель предлагаемого алгоритма — преобразовать исходную систему (1) к виду, пригодному для применения алгоритма Абрамова и Бронштейна, т. е. получить систему вида (2) с выделенной частью неизвестных. Как уже было сказано во введении, для случая невырожденной ведущей матрицы A1 исходной системы (1), это легко проделать.
0
m
m
0
0
0
Основную сложность представляют системы с вырожденной ведущей матрицей, поэтому в дальнейшем будем рассматривать именно этот случай.
На вход алгоритма поступает дифференциальная система с приведенной по строкам операторной матрицей. При этом, так как ведущая матрица вырождена, вектор порядков строк § будет содержать нулевые компоненты, а система разделится на две части: дифференциальную — состоящую из уравнений, соответствующих единичным элементам §, и алгебраическую — соответствующую нулевым элементам §.
Суть предлагаемого алгоритма состоит в том, чтобы из дифференциальной части исключить некоторые неизвестные, используя уравнения из алгебраической части, и тем самым получить дифференциальную систему с числом уравнений, равным числу неизвестных, и имеющую обратимую ведущую матрицу. При исключении необходимо учитывать, какие неизвестные были изначально выделены. В связи с этим, исключения проводятся в два этапа: на первом этапе из дифференциальной части системы исключаются невыделенные неизвестные, а на втором — выделенные.
Алгоритм состоит из трех этапов:
этап 1: исключение (из дифференциальной части системы) невыделенных неизвестных;
этап 2: исключение выделенных неизвестных;
этап 3: выражение исключенных выделенных неизвестных через оставшиеся в дифференциальной системе выделенные неизвестные.
Первый этап выполняется, если среди алгебраических уравнений системы есть такое, в которое входит какая-либо невыделенная неизвестная с ненулевым коэффициентом. Этот этап состоит из следующих шагов:
1. Из множества алгебраических уравнений выбрать уравнение, в которое какая-нибудь невыделенная неизвестная входит с ненулевым коэффициентом.
2. Используя выбранное на предыдущем шаге уравнение, произвести исключение этой невыделенной неизвестной из всех остальных уравнений системы как дифференциальных, так и алгебраических.
3. Использованное алгебраическое уравнение удалить из системы.
4. Если в системе остались алгебраические уравнения, в которые входят невыделенные неизвестные с ненулевыми коэффициентами, то вернуться к шагу 1.
В терминах операторных матриц исключение, выполняемое на шаге 2 первого этапа, можно описать как домножение операторной матрицы, соответствующей дифференциальной системе, слева на унимодулярную матрицу. Представим операторную матрицу, соответствующую системе, в виде
L11 . . . L1m
L
L
m1
Lm
Будем считать, что используемое алгебраическое уравнение соответствует ш-й строке этой матрицы (огёЬту = 0 для з = 1,...,ш), Ьтт = 0 ш
гда можно добиться перестановкой уравнений и перенумерацией неизвестных). Соответствующая унимодулярная матрица будет иметь вид:
U =
1 0 0 '•.
0 ...
_Т т-1 L1mLmm
1 —Т Т-1
1 Lm 1,mL
0
Jm-1,mT 1
(5)
Если после первого этапа в системе не осталось алгебраических уравнений, то алгоритм заканчивает свою работу. Если же еще остались алгебраические уравнения, то в них с ненулевыми коэффициентами будут присутствовать только выделенные неизвестные. Исключение из дифференциальной части некоторых из выделенных неизвестных осуществляется на втором этапе работы алгоритма.
Второй этап состоит из следующих шагов:
1. Из множества алгебраических уравнений выбрать произвольное уравнение.
2. Используя выбранное уравнение, произвести исключение какой-либо выделенной неизвестной, входящей в это уравнение с ненулевым коэффициентом, из всех остальных уравнений системы также, как и на шаге 2 первого этапа.
3. Использованное уравнение удалить из системы. Впоследствии (на третьем этапе алгоритма) оно найдет применение для выражения выделенных неизвестных, не вошедших в результирующую дифференциальную систему, через оставшиеся в ней выделенные неизвестные.
4. Если в системе остались алгебраические уравнения, то вернуться к шагу 1.
На заключительном третьем этапе алгоритма строится матрица Т, посредством которой выделенные неизвестные, не вошедшие в результирующую дифференциальную систему, выражаются через выделенные неизвестные, вошедшие в Т
I — числу выделенных неизвестных, не вошедших в результирующую дифференциальную систему (что соответствует числу исключений на втором этапе), а число столбцов равно числу выделенных неизвестных, попавших туда.
Третий этап состоит из следующих шагов:
1. Если на втором этапе были удалены все выделенные неизвестные, то
Т=
0
0
алгоритма), поэтому они исключаются. Тем самым осуществляется переход к матрице II е К1хк, где к — количество выделенных неизвестных в исходной системе.
4. Оставшиеся столбцы упор
Для дальнейшего прочтения статьи необходимо приобрести полный текст. Статьи высылаются в формате PDF на указанную при оплате почту. Время доставки составляет менее 10 минут. Стоимость одной статьи — 150 рублей.