Справочная информация
Задача Коши для нормальной системы обыкновенных дифференциальных уравнений записывается в виде
,
или в матричной форме
, ,
где
, , .
Система дифференциальных уравнений связывает независимую переменную x, искомые функции y1, y2,..., yn и их первые производные. В данном случае решение задачи Коши заключается в отыскании функций y1 = y1(x), y2 = y2(x),..., yn = yn(x), обращающих каждое уравнение системы в тождество на конечном или бесконечном интервале (a, b) и удовлетворяющих начальным условиям.
Такая форма записи задачи Коши является канонической для систем обыкновенных дифференциальных уравнений. К ней могут быть приведены как любые другие формы представления систем дифференциальных уравнений, разрешённых относительно старших производных, так и дифференциальные уравнения высших порядков. Приведение дифференциальных уравнений высших порядков к нормальной системе дифференциальных уравнений осуществляется по следующей схеме. Если дана задача Коши следующего вида
,
, , , … , ,
то замена переменных
, , , … , ,
сводит её к нормальной системе дифференциальных уравнений с начальными условиями
, ,
образующих задачу Коши.
Для решения такой задачи Коши используются те же методы, что для обыкновенных дифференциальных уравнений 1-го порядка. Это обуславливается тем, что матричная форма записи задачи Коши для нормальной системы полностью совпадает с её формулировкой для этих уравнений. Аналогична для неё и теорема о существовании единственного решения. Единственным отличием здесь является то, что вместо функций y(x) и f(x, y) используются вектор-функции y и f, состоящие из n функций y1(x), y2(x),..., yn(x) и f1(x, y1,..., yn), f2(x, y1,..., yn),..., fn(x, y1,..., yn), соответственно. При этом расчётные схемы методов и оценки их погрешностей сохраняются.
Метод Эйлера
Соотношения метода Эйлера для нормальной системы в матричной форме имеют вид
,
или в развёрнутой форме
, ,
где верхний индекс показывает номер шага по аргументу x.
Геометрическая интерпретация работы метода Эйлера решения задачи Коши для нормальной системы идентична его геометрической интерпретации для дифференциальных уравнений 1-го порядка. Однако в данном случае движение осуществляется вдоль некоторой кривой в (n + 1)-мерном пространстве переменных x, y1, y2,..., yn, которая является геометрическим представлением вектор-функции y.
Дата добавления: 2022-02-05; просмотров: 224;