Донецкий национальный технический университет
Опубликован: 15.03.2007 | Доступ: свободный | Студентов: 5956 / 2116 | Оценка: 4.11 / 3.78 | Длительность: 12:32:00
Специальности: Математик
Лекция 12:

Компьютерное моделирование и решение нелинейных уравнений

< Лекция 11 || Лекция 12: 12345678910

Решение дифференциальных уравнений высоких порядков

Методы Рунге-Кутта можно использовать не только для решения дифференциальных уравнений первого порядка

y'=f(x,y)
но и для решения дифференциальных уравнений более высоких порядков

y''=f(x,y,y'),\\ y'''=f(x,y,y',y''),\\ \ldots\\ y^{(m)}=f(x,y,y',y''',\ldots,y^{(m-1)}).

Любое дифференциальное уравнение m-го порядка

y^{(m)}=f(x,y,y',y",\ldots,y^{(m-1)}) ( 12.8)

можно свести к системе, состоящей из m уравнений первого порядка при помощи замен.

Заменим:

y_1=y',\\ y_2=y''=y'_1,\\ y_3=y'''=y'_2,\\ \ldots\\ y_m=y^{(m)}=y'_{(m-1)}.

В результате дифференциальное уравнение m -го порядка (12.8) сводится к системе, состоящей из m дифференциальных уравнений первого порядка:

\left\{ \begin{array}{l}
y'=y_1,\\ y'_1=y_2,\\ y'_2=y_3,\\ \ldots\\ y'_{m-1}=f(x,y,y_1,y_2,\ldots,y_{m-1}). \end{array} \right. ( 12.9)

Решением системы (12.2), а значит и дифференциального уравнения m -го порядка (12.1) является m табличных функций y, y_1=y', y_2=y_1^{''}, \ldots, y_m=y_{(m-1)}

Решение дифференциальных уравнений второго порядка

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

Общий вид дифференциальных уравнений второго порядка:

F(x,y,y',y'')=0 ( 12.10)

Нормальная форма дифференциальных уравнений второго порядка:

y''=f(x,y,y') ( 12.11)

Пример

Уравнение в общем виде

x^2y'' - xy' + (x^2 - 1)y = 0.

Его нормальная форма

y''=\frac{1}{x}y' - \frac{x^2-1}{x^2}y. ( 12.12)

Дифференциальное уравнение второго порядка (12.11) можно свести к системе, состоящей из двух дифференциальных уравнений первого порядка при помощи замен.

Заменим y1=y',

Тогда y'1=y".

В результате уравнение (12.11) сводится к системе, состоящей из двух дифференциальных уравнений первого порядка:

\left\{ \begin{array}{l} y_1=y'\\ y'_1=f(x,y,y_1). \end{array} \right.

Для примера (12.12) эта система имеет вид:

\left\{ \begin{array}{l} y'=y_1\\ y_1^{'}=\frac{1}{x}y_1-\frac{x^2-1}{x^2}y \end{array} \right. ( 12.13)

Решением этой системы являются две функции y(x) и y1(x),

где

y_1(x)=y'(x).

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

Дана система

\left\{ \begin{array}{l} y'=f_1(x,y,y_1),\\ y_1^{'}=f_2(x,y,y_1). \end{array} \right. ( 12.14)

Даны два начальных условия:

y(x_0) = y_0,\\
y_1(x_0)=(y_1)_0.

Необходимо проинтегрировать систему на участке [a, b] с шагом h.

В численных методах задача Коши для системы (12.14) ставится следующим образом:

Найти табличные функции

y_i(x_i) и (y_1)_i(x_i),  i=\overline{1,n}, т.е. найти таблицу

i x y y1
0 x0 y0 (y1)0
1 x1 y1 (y1)1
2 x2 y2 (y1)2
3 x3 y3 (y1)3
 ...  ...  ...  ...
n xn yn (y1)n

Здесь

h - шаг интегрирования дифференциального уравнения,

a=x0 - начало участка интегрирования уравнения,

b=xn - конец участка,

n=(b-a)/h - число шагов интегрирования уравнения.

На графике решением задачи Коши для системы, состоящей из двух дифференциальных уравнений первого порядка, является совокупность узловых точек (рис. 12.15).

При этом на каждом шаге, т.е. для каждого значения xi решением являются две узловые точки с координатами (xi, yi), (xi, (y1)i).


Рис. 12.15.

Для решения системы дифференциальных уравнений используем те же методы, что и для решения одного дифференциального уравнения первого порядка. При этом необходимо соблюдать условие: на каждом шаге интегрирования, т. е. в точках с координатами х1 , х2 , х3 , : , хn все уравнения системы надо решать параллельно.

< Лекция 11 || Лекция 12: 12345678910
Равиль Султанов
Равиль Султанов

В уравнениях движения кривошипно-шатунного механизма вместо обозначения радиуса кривошипа "r" ошибочно записан символ "γ" (гамма).

P.S. Может быть это слишком очевидно, но не упомянуто, что угол поворота кривошипа φ считается малым.

Александр Никитин
Александр Никитин

Добрый день.

В расчете параметра Т4 xi суммируется с величиной h/2 ?