2.3 Метод Рунге – Кутты
Методом Рунге-Кутты четвертого порядка точности называют одношаговый метод, относящийся к широкому классу методов Рунге-Кутты. В этом методе величины yi+1 вычисляются по следующим формулам:
yi+1 = yi + h (k1 + 2k2 + 2k3 + k4)/6 , i = 0, 1, ...
k1 = f(xi , yi),
k2 = f(xi+h/2, yi+hk1/2),
k3 = f(xi+h/2, yi+hk2/2),
k4 = f(xi+h, yi+hk3).
Найдем на [0, 1] приближенное решение задачи Коши y' = sin x – cos y, y(0)=1методом Рунге-Кутты 4-го порядка с шагом h=0.2 и методом Эйлера с тем же шагом.Изобразим оба приближенные решения графически
Для решения задачи методом Рунге-Кутты воспользуемся функцией rkfixed
Определим начальное условие - решение в начальной точке
Для того чтобы ввести нижний индекс переменной, щелкните по соответствующей позиции в панели Matrix или в панели Calculator
Определим правую часть уравнения
Знак присваивания можно ввести щелчком по соответствующей позиции в панели Evaluation.
Вычислим приближенное решение на отрезке [0,1], выполнив n=1/h=5 одинаковых шагов, методом Рунге-Кутты 4-го порядка; обозначим приближенное решение Y
Выведем в рабочий документ вычисленное приближенное решение
Для того чтобы вывести значение переменной в рабочий документ, введите имя переменной, знак равенства и щелкните по рабочему документу вне выделяюшей рамки
В первом столбце приведены значения x, во втором столбце - соответствующие значения приближенного решения
Решим ту же задачу методом Эйлера
Выведем в рабочий документ вычисленное приближенное решение, и, для сравнения, решение, вычисленное методом Рунге-Кутты
Построим графики приближенных решений
Для того чтобы построить график приближенного решения, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции возле оси абсцисс имя первого столбца матрицы Y, содержащего значения x в узлах сетки, а в позиции возле оси ординат - имя второго столбца, содержащего значения приближенного решения в узлах сетки; затем щелкните по свободному месту в рабочем документе вне поля графиков.
Для того чтобы одновременно построитьграфики нескольких функций от разных аргументов, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции у оси абсцисс имя первого аргумента, запятую, имя второго аргумента, и т.д., разделяя имена аргументов запятой.
Аналогично, в позиции возле оси ординат введите имя функции первого аргумента, запятую, имя функции второго аргумента и т.д. разделяя имена функций запятой.
Когда функции определены, щелкните по рабочему документу вне поля графиков.
Для того чтобы ввести номер столбца, щелкните по соответствующему символу в панели Matrix
Для того чтобы изменить стиль изображения, щелкните дважды по полю графиков и установите в окне соответствующие параметры
Заключение
Миллионы людей занимаются математическими расчетами, иногда в силу влечения к таинствам математики и ее внутренней красоте, а чаще в силу профессиональной или иной необходимости, не говоря уже об учебе. Ни одна серьезная разработка в любой отрасли науки и производства не обходится без трудоемких математических расчетов. Система Mathcad пользуется огромной популярностью во всем мире, позволяя готовить вполне профессиональные документы, имеющие вид статей и книг по математике.
Программа Microsoft Excel входит в офисный пакет Microsoft Office и предназначена для подготовки и обработки электронных таблиц под управлением операционной системой Windows. Microsoft Excel – это многофункциональный, мощный редактор электронных таблиц. Он представляет возможность производить различные расчеты, составлять списки, сметы и что немаловажно, строить наглядные графики и диаграммы.
MathCAD – это мощная и в то же время простая универсальная среда для решения задач в различных отраслях науки и техники, финансов и экономики, физики и астрономии, строительства и архитектуры, математики и статистики, организации производства и управления… Она располагает широким набором инструментальных, информационных и графических средств. Сегодня MathCAD – одна из самых популярных математических систем. Она пользуется большим спросом у студентов, инженеров, экономистов, менеджеров, научных работников и всех тех, чья деятельность связана с количественными методами расчета.
Microsoft Excel ‑ средство для работы с электронными таблицами, намного превышающее по своим возможностям существующие редакторы таблиц, первая версия данного продукта была разработана фирмой Microsoft в 1985 году. Microsoft Excel ‑ это простое и удобное средство, позволяющее проанализировать данные и, при необходимости, проинформировать о результате заинтересованную аудиторию, используя Internet. Microsoft® Excel разработан фирмой Microsoft, и является на сегодняшний день самым популярным табличным редактором в мире. Кроме стандартных возможностей его отличает следующие возможности, он выводит на поверхность центральные функции электронных таблиц и делает их более доступными для всех пользователей. Для облегчения работы пользователя упрощены основные функции, создание формул, форматирование, печать и построение графиков.
Данная курсовая работа позволила мне более близко познакомится с пакетом прикладных программ MathCAD и Microsoft Excel. Мной было рассмотрено несколько способов решения дифференциальных уравнений.
Всё это позволило в полном объеме усвоить лекционный материал и понять перспективы использования вычислительной техники при решении различных задач практического характера.
Список литературы
1. Индейкин В. В. Табличный редактор Microsoft Excel. Учебное пособие. – Казань, 1999. – 75с.
2. Кудрявцев Е. М. MathCAD 2000 Pro. – М.: ДМК Пресс, 2001. – 571с.
3. Руководство пользователя Mathcad 6.0 и Mathcad PLUS 6.0. – Компания SoftLine. http://www.exponenta.ru/soft/Mathcad/UsersGuide/0.asp
... в точке хк проводят вычисления ук с шагом h, затем с шагом 2h и берут 1/3 разницы этих значений: | ук*-у(хк)|=1/3(yk*-yk), (2.5.9) где у(х)-точное решение дифференциального уравнения. Таким образом, методом Эйлера можно решать уравнения любых порядков. Например, чтобы решить уравнение второго порядка y//=f(y/,y,x) c начальными условиями y/(x0)=y/0, y(x0)=y0, ...
... . , т.е. таких уравнений, у которых правая часть не является ненпрерывной по x функций рассмотрены в статье [5]. Теория систем автоматического управления, описываемых дифференциальными уравнениями с разрывными правыми частями рассматривается в книгах [13, 14, 15]. В работе С.В. Емельянова [13] излагается один из разделов теории автоматичесеого управления – теория систем с переменной структурой, ...
... пакетах. Заключение Результатом выполнения курсового проекта является готовый программный продукт, позволяющий решать задачу Коши для системы дифференциальных уравнений при помощи неявной схемы Адамса 3-го порядка, демонстрирующий возможности численного решения поставленной задачи с заданной степенью точности. Готовый программный продукт может найти широкое применение при решении многих ...
... силы взаимодействия между рыбой и рабочими органами машин, поскольку изменяется площадь контакта, обусловливающая силы трения. До настоящего времени структурно-механические характеристики в основном оцениваются органолептическим методом. Рыбу сдавливают пальцами и оценивают ее консистенцию. Не достатком такого метода оценки структурно-механических характеристик мышечной ткани является его ...
0 комментариев