|
Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология Образование Политология Производство Психология Стандартизация Технологии |
Лекция 3. Численные методы решения систем линейных дифференциальных уравнений
Рассмотрим задачу Коши для линейной динамической системы
с заданным для каждого момента времени вектором входа
представляющего собой систему линейных неоднородных дифференциальных уравнений. Для построения решения неоднородной системы предварительно изучим свойства решений линейной однородной системы
Для этого рассмотрим n-линейно независимых векторов
которые примем в качестве начальных условий линейной однородной системы
По этим n начальным условиям получим n решений однородной системы дифференциальных уравнений
каждое из которых удовлетворяет соотношению
Матрица
Любой столбец фундаментальной матрицы удовлетворяет линейной системе, поэтому
то есть фундаментальная матрица удовлетворяет матричному линейному дифференциальному уравнению
Основное свойство фундаментальной матрицы дает следующая теорема. Теорема: если существует В связи с тем, что векторы начальных условий линейно независимы, то Система линейных дифференциальных уравнений имеет бесконечное число фундаментальных матриц в зависимости от принятого набора линейно независимых векторов начальных условий. Так как
называемую переходной матрицей. Переходная матрица обладает следующими свойствами: 1. 2. Переходная матрица невырожденная для любых
3. Переходная матрица удовлетворяет матричному дифференциальному уравнению
Действительно:
Из свойств 3 и 1 следует, что 4. Действительно
Переходную матрицу используют для построения решения систем линейных однородных и неоднородных дифференциальных уравнений.
Решение задачи Коши систем линейных однородных дифференциальных уравнений через фундаментальную матрицу Рассмотрим следующую задачу Коцш:
Докажем, что решение этой задачи имеет вид
Доказательство. Продифференцируем по t левую и правую часть решения:
Второе равенство следует из свойства 3 переходной матрицы, третье равенство – из определения решения через переходную матрицу. Таким образом, предложенное решение удовлетворяет системе линейных однородных дифференциальных уравнений. Осталось показать, что предложенное решение удовлетворяет начальному условию, для чего вычислим
Следовательно, переходная матрица позволяет вычислить решение ЛОУ для любого момента времени t через решение в другой момент времени
Вычисление переходной матрицы Переходная матрица удовлетворяет матричному дифференциальному уравнению
с начальным условием Требуется найти аналитическое решение этого матричного уравнения. Используем следующий прием: проинтегрируем обе части уравнения в пределах
Учитывая, что в левой части содержится полный дифференциал, а матрица А не зависит от
или
Подынтегральное выражение можно представить по этой же формуле, то есть
Подставив (4.7) в подынтегральное выражение, получим:
продолжая итерационно этот процесс, получим выражение для переходной матрицы в виде сходящегося матричного ряда
Рассмотрим скалярную функцию
Производные от экспоненты вычисляют по формуле
В результате получим для экспоненты
откуда
Из сравнения полученного ряда с выражением для переходной матрицы следует, что переходная матрица представляет собой ряд для матричной экспоненты
Матричная экспонента обладает всеми свойствами, присущими скалярной экспоненциальной функции, в частности
Популярное:
|
Последнее изменение этой страницы: 2016-03-17; Просмотров: 1352; Нарушение авторского права страницы