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

↑

⇐ ПредыдущаяСтр 7 из 9Следующая ⇒

При переходе от решения одиночных дифференциальных уравнений к реше­нию систем дифференциальных уравнений сложность решения быстро нараста­ет. Уже при решении системы из двух дифференциальных уравнений значительно усложняются рекуррентные формулы, определяющие коэффициенты в формулах Рунге—Кутта. . При этом добавление каждый раз очередного уравнения увеличивает число уравнений в их векторной записи. Естественно, это увеличивает сложность решения, ведет к увеличению .числа переменных в функциях, задающих коэффициенты Ki, но принципиально не меняет реализации алгоритма вычислений.

При решении систем дифференциаль­ных уравнений «в лоб» этот путь оказывается тупиковым из-за быстрого нара­стания сложности уравнений и их частного характера.

Встроенные в «Mathcad» функции для решения систем дифференциальных урав­нений обходят эти трудности, поскольку не требуют составления формул для решения систем дифференциальных уравнений и задают такое решение в общем виде. Надо, однако, учитывать, что реальные алгоритмы решения диффе­ренциальных уравнений оказываются заметно сложнее рассмотренных. Так, широкое распространение получили адаптивные алгоритмы, в которых на каж­дом шаге выполняется контроль погрешности решения (функция Rkadapt в «Mathcad»). Если погрешность падает ниже заданной, шаг h уменьшается. И так до тех пор, пока погрешность не станет меньше заданной. Это заметно снижает неустойчивость решения, но одновременно увеличивает время вычислений как на каждом шаге, так и в области возможной неустойчивости решения.

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

rkfixed(y,x1,x2,N,D) – матрица решений методом Рунге – Кутта системы обыкновенных дифференциальных уравнений с начальными условиями в векторе y  на интервале от x1 до x2 с фиксированным числом шагов N; правые части уравнений записаны в векторе D.

Rkadapt(y,x1,x2,N,D) – матрица решений методом Рунге – Кутта (с переменным шагом) системы обыкновенных дифференциальных уравнений с начальными условиями в векторе y  на интервале от x1 до x2 с фиксированным числом шагов N; правые части уравнений записаны в векторе D.

В качестве примера рассмотрим динамическую модель электрической схемы, изображенной на рисунке 2.

Электрическое равновесие в предложенной схеме определяется системой дифференциальных уравнений

Познавательные статьи:



Как правильно слушать собеседника

Типичные ошибки при выполнении бросков в баскетболе

Принятие христианства на Руси и его значение

Средства массовой информации США