Механическое удерживание земляных масс: Механическое удерживание земляных масс на склоне обеспечивают контрфорсными сооружениями различных конструкций...
Семя – орган полового размножения и расселения растений: наружи у семян имеется плотный покров – кожура...
Топ:
Марксистская теория происхождения государства: По мнению Маркса и Энгельса, в основе развития общества, происходящих в нем изменений лежит...
Отражение на счетах бухгалтерского учета процесса приобретения: Процесс заготовления представляет систему экономических событий, включающих приобретение организацией у поставщиков сырья...
Интересное:
Аура как энергетическое поле: многослойную ауру человека можно представить себе подобным...
Средства для ингаляционного наркоза: Наркоз наступает в результате вдыхания (ингаляции) средств, которое осуществляют или с помощью маски...
Как мы говорим и как мы слушаем: общение можно сравнить с огромным зонтиком, под которым скрыто все...
Дисциплины:
2017-05-23 | 270 |
5.00
из
|
Заказать работу |
|
|
В 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 - Не является автором материалов. Исключительное право сохранено за автором текста.
Если вы не хотите, чтобы данный материал был у нас на сайте, перейдите по ссылке: Нарушение авторских прав. Мы поможем в написании вашей работы!