Информационная система "Конференции"


Международная конференция молодых ученых по математическому моделированию и информационным технологиям

29-31 октября 2002 г., Новосибирск, Академгородок

Гайдомак С.В.
Институт динамики систем и теории управления СО РАН, Иркутск

О заменах переменных в системах уравнений не типа Коши-Ковалевской

  • Введение
  • Рассматривается система дифференциальных уравнений

    где -постоянные матрицы, искомая и заданная вектор-функции соответственно, . Предполагается, что

    и вектор-функция достаточно гладкая в области . Систему вида (1), удовлетворяющую условию (2), будем называть системой не типа Коши-Ковалевской. Еще И.Г. Петровский указывал на необходимость изучения общих систем дифференциальных уравнений, неразрешенных относительно старшей производной по времени. Такие системы возникают в различных прикладных задачах: в гидродинамике, газовой динамике, теории малых колебаний жидкости, в теории тепломассообмена [1], [2], [3].

    Для системы (1) с нулевой матрицей ранее была получена теорема существования и единственности начально-краевой задачи и сформулирован ряд важных утверждений [4]. Поскольку применение разностных схем, разработанных и обоснованных для систем в нормальной форме, к системам не типа Коши-Ковалевской не всегда возможно, то возникает задача нормализации: приведения системы (1) к виду, разрешенному относительно производных . Этот вопрос решён для случая регулярного пучка и нулевой матрицы [5].

    Поэтому при исследовании системы (1) с ненулевой матрицей естественен вопрос о ее приведении к системе у которой, используя соответствующую замену переменной. Показать возможность такого перехода, а также его применение, снабдив слушателя предварительными фактами о системах не типа Коши-Ковалевской, когда , является целью доклада.

  • Теорема существования для системы с нулевой матрицей С.
  • Будем предполагать, что в системе (1) пучок матриц регулярен, где под символами и понимаются некоторые числовые параметры или операторы.

    Определение 1.

    Пучок матриц регулярен, если ,

    где числовые параметры из .

    Лемма 1[6].

    Пусть пучок матриц регулярен. Тогда существуют невырожденные -матрицы и с постоянными элементами такие, что

    где нильпотентные матрицы индекса соответственно: , некоторая -матрица, единичные матрицы размерности равной индексу,

    Определение 2[6].

    Первые два слагаемые из правой части равенства (3) называется канонической (кронекеровой) структурой пучка матриц . Параметры называются индексами пар матриц и соответственно.

    Лемма 2[6].

    Пучок матриц регулярен тогда и только тогда, когда существует матричный многочленсо свойством где-постоянные матрицы.

    Лемма 3.

    Общее решение системы уравнений

    где искомая вектор-функция, заданная вектор-функция, имеет вид

    где произвольная вектор-функция, .

    Доказательство существования решения в виде (5) можно провести либо непосредственным вычислением, либо используя лемму 2, предварительно проинтегрировав систему (4). Единственность доказывается методом от противного.

    Следствие 1.

    Общее решение уравнения имеет вид:

    Замечание 1. Вместо (4) рассмотрим систему

    где скалярные операторы с постоянными коэффициентами. Тогда можно построить аналог формулы (5) и операторы в ней будут иметь вид где оператор является левым обратным к оператору Во втором слагаемом оператор действует не просто на произвольную вектор-функцию , а на ядро оператора

    Перейдем к анализу разрешимости системы (1) в предположении, что Для этого воспользуемся леммой 1. Умножая систему (1) на матрицу и, производя замену, мы расщепим ее на три подсистемы:

    где .

    Очевидно, что задавая условия для второго и третьего уравнений системы (6), мы согласно лемме 3 и следствию к ней однозначно определим их решения в области :

    Следует отметить, что при задании начальных условий на других границах области приводит к невозможности произвольного выбора задаваемых функций. Рассмотрим для простоты случай однородных систем. С необходимостью должны выполняться соотношения

    где произвольные вектор-функции. Из (9) следует, что при выполнении условий совместности решение задачи Коши системы (4) существует, но не единственно.

    Лемма 4.

    Пусть: 1) многочлен ненулевой и его ненулевые корни положительные и простые;

    2)

    где ; 3) ; ; 4)

    Тогда в области система (1) имеет единственное решение.

    Доказательство. Решения двух последних уравнений (6) можно получить подставляя начальные условия и правые части в формулы (7), (8). Первая подсистема имеет вид:. Пусть число ненулевых корней многочлена равно . Тогда матрицы в лемме 1 можно выбрать так, что и первая подсистема системы (6) будет гиперболической [7]. При наших предположениях характеристики этой системы: где произвольная постоянная, образуют острый угол с положительным направлением оси по , начальные и краевые условия совместны, следовательно краевая задача разрешима и имеет единственное решение [7]. Лемма доказана.

    Замечание 2. Зачастую на практике начальные и краевые условия задаются в виде тогда для того, чтобы краевая задача была совместна необходимо и достаточно выполнение условий

    В лемме 4 для упрощения выкладок требования на гладкость свободного члена и краевых условий существенно загрублены. В принципе можно учесть в требованиях на гладкость размеры нильпотентных блоков.

  • Теорема существования для системы с ненулевой матрицей С.
  • Перейдем к утверждению, которое в некоторых случаях, позволяет осуществить переход от системы с ненулевой матрицей к системе у

    которой

    Лемма 5.

    Если система матричных уравнений

    разрешима относительно матриц и , то заменой переменных она сводится к системе

    Доказательство. После замены переменной система (1) примет вид где В силу предположения о разрешимости системы (11) относительно матриц и матричная экспонента перестановочна с и , Лемма доказана.

    Следствие 2. Если найдется матрица и то в лемме 5 можно принять Если то можно принять

    Теорема 1.

    Пусть: 1) выполнены условия 1), 2), 3) леммы 4;

    2) в произведении из равенства (4) блоки

    3) содержит не более двух нулевых блоков вне диагонали и расположенных в одной строке или одном столбце, причем блоки и (не совпадающие с выделенными нулевыми) связаны условиями: или

    4) Тогда в области система (1) имеет единственное решение.

    Доказательство. Пусть После умножения системы (1) на матрицу и замены получим

    Рассмотрим последнюю подсистему, в предположении, что Воспользуемся следствием леммы 5 и выполним замену переменной где получим равенство

    Откуда Подставляем во второе уравнение получим систему

    в которой искомой является только вектор-функция определяемая по формуле (7) после соответствующей замены. Здесь учтено, что в силу свойств блоков и и перестановочности матриц

    После этого мы можем найти Так как компоненты и вычислены, их можно подставить в первую подсистему, которая является гиперболической и удовлетворяет условиям известной теоремы [7]. Другие случаи доказываются аналогично. Теорема доказана.

    Безусловно, проверка условий 2), 3) теоремы 1 связана с большими трудностями, поэтому желательно иметь критерии в терминах входных данных. Приведем один такой частный критерий.

    Лемма 6.

    Для того, чтобы матрица имела структуру со свойствами из теоремы 1 необходимо, чтобы где параметры из равенства (4), а некоторые

    коэффициенты.

  • Литература
    1. Campbell S.L., Marzalek W. The Index of Infinite Dimensional Implicit System // Mathematical and Computer Modelling of System, 1999, vol.5, № 1, p.18-42.
    2. Демиденко Г.В., Успенский С.В. Уравнения и системы, неразрешенные относительно старшей производной. Новосибирск: Сибирская издательская фирма РАН "Наука", 1998.
    3. Таиров Э.А., Запов В.В. Интегральная модель нелинейной динамики парогенирирующего канала на основе аналитических решений // ВАНТ. Сер. Физика ядерных реакторов, 1991, вып. 3, с.14-20.
    4. Гайдомак С.В. О существовании решений системы дифференциальных уравнений не типа Коши-Ковалевской (II Всесибирский конгресс женщин математиков), тезисы докладов, 15-20 января 2002г., Красноярск).
    5. Чистяков В.Ф., Гайдомак С.В. О существовании нормализатора для системы дифференциальных уравнений не типа Коши-Ковалевской, тезисы докладов, ИДСТУ СО РАН, Иркутск, 2001.
    6. Бояринцев Ю.Е., В.Ф.Чистяков. Алгебро – дифференциальные системы. Методы численного решения и исследования.- Новосибирск: Сибирская издательская фирма РАН "Наука", 1998.
    7. Годунов С.К. Уравнения математической физики. М.: Наука, 1971. - 416 с.