Механическое удерживание земляных масс: Механическое удерживание земляных масс на склоне обеспечивают контрфорсными сооружениями различных конструкций...

Семя – орган полового размножения и расселения растений: наружи у семян имеется плотный покров – кожура...

Встроенные функции для решения системы ОДУ

2017-05-23 270
Встроенные функции для решения системы ОДУ 0.00 из 5.00 0 оценок
Заказать работу

В MathCAD имеются три встроенные функции, которые позволяют решать поставленную в форме (2-3) задачу Коши различными численными методами.

1. rkfixed(y0,t0,t1,M,D) –метод Рунге-Кутты с фиксированным шагом.

2. Rkadapt(y0,t0,t1,M,D) - метод Рунге-Кутты с переменным шагом.

3. Bulstoer(y0,t0,t1,M,D) – метод Булирша-Штера.

В этих функциях: y0 – вектор начальных значений в точке t0 размера N×1; t0 – начальная точка расчета; t1 – конечная точка расчета; M – число шагов, на которых численный метод находит решение; D – векторная функция размера N×1 двух аргументов – скалярного t и векторного y. При этом y – искомая векторная функция аргумента t того же размера N×1.

Каждая из приведенных функций выдает решение в виде матрицы размера (M+1)×(N+1). В ее левом столбце находятся значения аргумента t, а в остальных N столбцах - значения искомых функций y0(t), y1(t),…, yN-1(t), рассчитанные для этих значений аргумента. Всего в этой матрице M+1 строка.

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

Задание 4. Реализуйте следующий пример, содержащий решение системы ОДУ модели осциллятора с затуханием. Проанализируйте результаты решения.

 

 

Функция D(t,y), определяющая систему ОДУ, имеет ряд особенностей:

1. Функция D должна быть обязательно функцией двух аргументов.

2. Второй ее аргумент должен быть вектором того же размера, что и сама функция D.

Точно такой же размер должен быть и у вектора начальных значений y0.

 

Решение системы осуществлено на промежутке (0, 50). Матрица решений имеет размер (M+1)×(N+1), т.е. 101×3. Просмотреть все компоненты матрицы u можно с помощью вертикальной полосы прокрутки, которая появляется при щелчке мышкой по таблице решений. Так расчетное значение первой искомой переменной y0 на 12 шаге u12,1=0.07. Для вывода элементов решения в последней точке интервала используется выражение типа uM,1=7.523·10-3.

Чтобы построить график решения, надо отложить соответствующие компоненты решения по координатным осям: значение аргумента u<0> - вдоль оси X, а u<1> и u<2> - вдоль оси Y (Чтобы разместить на рисунке два графика параметры u<1> и u<2> для оси Y надо вводить через запятую).

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

Задание 5. По условиям задание 4 постройте фазовый портрет. Постройте второй фазовый портрет при M= 200. Проанализируйте результаты решения.

В общим случае, если система состоит из N ОДУ, то фазовое пространство является N -мерным. При N>3 наглядность теряется, и для визуализации фазового портрета приходится строить его различные проекции.

Все численные методы решения ОДУ основаны на аппроксимации дифференциальных уравнений разностными аналогами. В зависимости от конкретной формы аппроксимации, получаются алгоритмы различной точности и быстродействия. В MathCAD использован наиболее популярный алгоритм Рунге-Кутты четвертого порядка, описанный в большинстве книг по методам вычислений. Он обеспечивает малую погрешность для широкого класса систем ОДУ за исключением жестких систем. Поэтому в большинстве случаев применяют функцию rkfixed. Если по различным причинам время расчетов становится критичным или точность неудовлетворительна, стоит попробовать применить другие встроенные функции.

Функция Rkadapt может быть полезна в случае, когда известно, что решение на рассматриваемом интервале меняется слабо, либо существуют участки медленных и быстрых изменений. Метод Рунге-Кутты с переменным шагом разбивает интервал не на равномерные шаги, а более оптимальным способом. Там, где решение меняется слабо, шаги выбираются более редкими, а в областях его сильных изменений – частыми. В этом случае для достижения одинаковой точности требуется меньшее число шагов. Метод Булирша-Штера (Bulstoer) часто оказывается более эффективным для поиска гладких решений.


Поделиться с друзьями:

Эмиссия газов от очистных сооружений канализации: В последние годы внимание мирового сообщества сосредоточено на экологических проблемах...

Опора деревянной одностоечной и способы укрепление угловых опор: Опоры ВЛ - конструкции, предназначен­ные для поддерживания проводов на необходимой высоте над землей, водой...

Организация стока поверхностных вод: Наибольшее количество влаги на земном шаре испаряется с поверхности морей и океанов (88‰)...

Индивидуальные очистные сооружения: К классу индивидуальных очистных сооружений относят сооружения, пропускная способность которых...



© cyberpedia.su 2017-2024 - Не является автором материалов. Исключительное право сохранено за автором текста.
Если вы не хотите, чтобы данный материал был у нас на сайте, перейдите по ссылке: Нарушение авторских прав. Мы поможем в написании вашей работы!

0.007 с.