Системы дифференциальных уравнений более высокого порядка
Методика решения системы дифференциальных уравнений, каждое из которых может содержать производные от неизвестных функций выше первого порядка, та же, что описана выше. Заметим, что любое уравнение вида
x(n) = f( t, x(n-1), x(n-2),..., x,...)
посредством замены
x1(t) = x(t) |
x2(t) = x'(t) |
... |
xn(t) = x(n-1)(t) |
может быть приведено к совокупности уравнений
x'1(t) = x2(t) |
x'2(t) = x3(t) |
... |
x'n(t) = f( t, xn, xn-1,... x1,...) |
Заметьте, что в приведённых выше уравнениях не встречаются производные выше первого порядка. Преобразовав подобным образом каждое из уравнений, входящих в исходную систему, получим систему с большим количеством неизвестных функций, но с производными только первого порядка. Методика решения таких систем описана в предыдущем подразделе. Пример на Рисунке 6 показывает, как поступать при решении системы дифференциальных уравнений второго порядка:
u''(t) = 2v(t)
v''(t) = 4v(t) - 2u(t)
Рисунок 6: Система линейных дифференциальных уравнений второго порядка.
Функция rkfixed возвращает матрицу, в которой:
Первый столбец содержит точки, в которых должны быть найдены решения и их производные.
Остальные столбцы содержат значения решений и их производных, соответствующие точкам из первого столбца. Порядок, в котором появляются решение и его производные, повторяет порядок их расположения в функции D(x,y) и векторе начальных условий y.
Функция rkfixed, описанная выше, использует универсальный метод для решения дифференциальных уравнений. Хотя метод Рунге-Кутты, используемый этой функцией, не всегда является самым быстродействующим, он почти всегда приводит к искомому результату. Однако имеются случаи, в которых можно попробовать использовать некоторые специальные функции Mathcad, предназначенные для решения дифференциальных уравнений. Эти случаи относятся к трем широким категориям задач:
Решается не задача с начальными условиями (задача Коши), а краевая задача (sbval и bvalfit).
Cистема дифференциальных уравнений может иметь некоторые специфические свойства, используя которые, можно решить её более точно. Для решения системы в этих случаях лучше использовать специализированные функции Mathcad, чем функцию rkfixed. Например, для решения жесткой системы (Stiffb, Stiffr); для решения систем с гладкими решениями (Bulstoer) или медленно изменяющимися решениями (Rkadapt).
Требуется найти решение только в одной фиксированной точке (bulstoer, rkadapt, stiffb и stiffr).
Можно также пробовать решать дифференциальное уравнение разными методами, чтобы понять, какой из них работает лучше всего. Иногда имеются тонкие различия между дифференциальными уравнениями, которые делают один метод решения предпочтительнее другого.
Следующие разделы описывают использование различных методов решения дифференциального уравнения и те случаи, в которых они могут быть полезны.
- Буквенные индексы
- Ниже приводится полный список предопределенных переменных Mathcad и их значений по умолчанию:
- Используемые числа
- Специальные операции над комплексными числами
- Многозначные функции
- Создание вектора
- Создание матрицы
- Изменение размера матрицы
- Нижние индексы и элементы вектора
- Изменение способа отображения массивов
- Графическое представление матриц
- Ограничение входных массивов
- Ограничение отображаемых массивов
- Ограничение размеров массива
- Размеры и диапазон значений массива
- Специальные типы матриц
- Специальные характеристики матрицы
- Формирование новых матриц из существующих
- Собственные значения и собственные векторы
- Разложения
- Решение линейной системы уравнений
- Определение составного массива
- Отображение составных массивов
- Операторы и функции для составных массивов
- Определение и использование дискретного аргумента
- Многократные вычисления по дискретному аргументу
- Множественные дискретные аргументы и двойные индексы
- Рекурсивные вычисления с несколькими переменными
- Рекурсивные вычисления с вектором
- Советы по набору операторов
- Переменный верхний предел суммирования
- Оператор суммирования элементов вектора
- Производные более высокого порядка
- Переменные пределы интегрирования
- Изменение точности вычисления интегралов
- Криволинейные и двойные интегралы
- Определение пользовательского оператора
- Использование пользовательского оператора
- Запись функций как операторов
- Тригонометрические функции и обратные им.
- Гиперболические функции
- Логарифмические и показательные функции
- Функции Бесселя
- Специальные функции
- Введение в дискретное преобразование Фурье
- Функция if
- Циклы “while”
- Оператор “break”
- Циклы “for”
- Подпрограммы
- Рекурсия
- Что делать, когда функция root не сходится
- Некоторые советы по использованию функции root
- Решение уравнений с параметром
- Нахождение корней полинома
- Как использовать найденное решение
- Что делать, когда Mathcad не может найти решения
- Что делать, когда имеется слишком мало ограничений
- Многократное решение уравнений
- Решение одинаковых задач относительно разных переменных
- Приближенные решения
- Использование символьного решения уравнений
- Дифференциальные уравнения первого порядка
- Дифференциальные уравнения второго порядка
- Уравнения более высокого порядка
- Системы оду первого порядка
- Системы дифференциальных уравнений более высокого порядка
- Гладкие системы
- Медленно изменяющиеся решения
- Нахождение приближенного решения только в конечной точке
- Двухточечные краевые задачи
- Дифференциальные уравнения с частными производными