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

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

(5.18)

Здесь – искомые функции, значения которых подлежат определению при . В момент времени х = х₀ задаются начальные условия

, (5.19)

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

Введем следующие две вектор-функции (т.е. значение которых – массив)

,

и вектор (массив) . Тогда задачу Коши (5.18), (5.19) можно записать в компактной форме

,

.

Описанные выше численные методы решения задачи Коши для одного уравнения можно использовать и для систем уравнений первого порядка, причем форма их записи претерпевает минимальные изменения. В расчетных формулах следует лишь заменить числа y_i на векторы (массивы) и все функции – на соответствующие вектор-функции (процедуры). Например, расчетная формула метода Эйлера применительно к решению системы (5.18) принимает вид

.

Покоординатная запись этого соотношения выглядит так:

Аналогично, формулы метода Рунге-Кутты 4-го порядка точности для систем дифференциальных уравнений первого порядка имеют вид:

,

; ;

; .

Решение уравнений высших порядков. Задача Коши для дифференциального уравнения n-го порядка состоит в нахождении функции y(x), удовлетворяющей при x ³ x₀ дифференциальному уравнению

(5.20)

а при x = x₀ – начальным условиям

(5.21)

Универсальным приемом решения этой задачи является ее сведение с помощью замены к задаче Коши для системы дифференциальных уравнений 1-го порядка

(5.22)

(5.23)

Для решения задачи Коши (5.20), (5.21), приведенной к виду (5.22), (5.23), можно воспользоваться известными методами.

Приложение

Дата добавления: 2014-01-07