Использование программы Mathcad для решения дифференциальных уравнений

Вначале покажем решение систем ОДУ первого порядка с использованием вычислительного блока Given/Odesolve [4].

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

Решение системы ОДУ показано на графике слева. Как известно, решения ОДУ часто удобнее изображать не в таком виде, а в фазовом пространстве, по каждой из осей которого откладываются значения каждой из найденных функций (как показано на рисунке справа). При этом аргумент входит в них лишь параметрически. В рассматриваемом случае двух ОДУ такой график – фазовый портрет системы – является кривой на фазовой плоскости. В общем случае, если система состоит из N ОДУ, то фазовое пространство является N – мерным. При N > 3 наглядность теряется, и для визуализации фазового портрета приходится строить его различные проекции.

Рассмотрим решение этой же системы ОДУ первого порядка с использованием встроенной функции rkfixed [5].

Полученное решение полностью соответствует вышеприведенному решению с использованием вычислительного блока Given/Odesolve. Следует отметить, что начальные условия здесь задаются в виде вектора y, а функциям x(t) и y(t) соответствуют элементы этого вектора y1 и y2. Вектор начальных условий y и вектор правых частей F имеют размер равный двум, т.к. система состоит из двух уравнений первого порядка. Для системы ОДУ, состоящей из двух уравнений второго порядка, размер этих векторов будет равен четырем [6].

Литература

1. Дифференциальные уравнения [Электронный ресурс]. Режим доступа: http://ru.wikipedia.org/wiki/Mathcad
2. Решение дифференциальных уравнений в среде Mathcad [Электронный ресурс]. Режим доступа: 

http://nice-ne.ru/publ/mathcad/reshenie_ differencialnykh_uravnenij_i_sistem_v_mathcad/18-1-0-923

3. Обыкновенные дифференциальные уравнения [Электронный ресурс]. Режим доступа: www.exponenta.ru