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

Определим начальное условие - решение в начальной точке

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50091.JPG

Для того чтобы ввести нижний индекс переменной, щелкните по соответствующей позиции в панели Matrix или в панели Calculator

Определим правую часть уравнения

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50094.JPG

Знак присваивания можно ввести щелчком по соответствующей позиции в панели Evaluation.

Вычислим приближенное решение на отрезке [0,1], выполнив n=1/h=5 одинаковых шагов, методом Рунге-Кутты 4-го порядка; обозначим приближенное решение Y

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50097.JPG

Выведем в рабочий документ вычисленное приближенное решение

Для того чтобы вывести значение переменной в рабочий документ, введите имя переменной, знак равенства и щелкните по рабочему документу вне выделяюшей рамки

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50100.JPG

В первом столбце приведены значения x, во втором столбце - соответствующие значения приближенного решения

Решим ту же задачу методом Эйлера

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50135.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50136.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50137.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50138.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50139.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50140.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50141.JPG

Выведем в рабочий документ вычисленное приближенное решение, и, для сравнения, решение, вычисленное методом Рунге-Кутты

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50133.JPG


Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50134.JPG

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50145.JPG

Построим графики приближенных решений

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50148.JPG

Для того чтобы построить график приближенного решения, щелкните в панели Graph по пиктограмме декартова графика, введите в помеченной позиции возле оси абсцисс имя первого столбца матрицы Y, содержащего значения x в узлах сетки, а в позиции возле оси ординат - имя второго столбца, содержащего значения приближенного решения в узлах сетки; затем щелкните по свободному месту в рабочем документе вне поля графиков.

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

Аналогично, в позиции возле оси ординат введите имя функции первого аргумента, запятую, имя функции второго аргумента и т.д. разделяя имена функций запятой.

Когда функции определены, щелкните по рабочему документу вне поля графиков.

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50150.JPG

Для того чтобы ввести номер столбца, щелкните по соответствующему символу в панели Matrix

Описание: http://www.exponenta.ru/educat/class/courses/ode/theme2/mathcad/ex5/ex50151.JPG

Для того чтобы изменить стиль изображения, щелкните дважды по полю графиков и установите в окне соответствующие параметры


Заключение

Миллионы людей занимаются математическими расчетами, иногда в силу влечения к таинствам математики и ее внутренней красоте, а чаще в силу профессиональной или иной необходимости, не говоря уже об учебе. Ни одна серьезная разработка в любой отрасли науки и производства не обходится без трудоемких математических расчетов. Система 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


Информация о работе «Решение дифференциальных уравнений. Обзор»
Раздел: Информатика, программирование
Количество знаков с пробелами: 21527
Количество таблиц: 0
Количество изображений: 12

Похожие работы

Скачать
24266
4
0

... в точке хк проводят вычисления ук с шагом h, затем с шагом 2h и берут 1/3 разницы этих значений: | ук*-у(хк)|=1/3(yk*-yk), (2.5.9) где у(х)-точное решение дифференциального уравнения.  Таким образом, методом Эйлера можно решать уравнения любых порядков. Например, чтобы решить уравнение второго порядка y//=f(y/,y,x) c начальными условиями y/(x0)=y/0, y(x0)=y0, ...

Скачать
41135
2
10

... . , т.е. таких уравнений, у которых правая часть не является ненпрерывной по x функций рассмотрены в статье [5]. Теория систем автоматического управления, описываемых дифференциальными уравнениями с разрывными правыми частями рассматривается в книгах [13, 14, 15]. В работе С.В. Емельянова [13] излагается один из разделов теории автоматичесеого управления – теория систем с переменной структурой, ...

Скачать
39446
2
12

... пакетах.   Заключение   Результатом выполнения курсового проекта является готовый программный продукт, позволяющий решать задачу Коши для системы дифференциальных уравнений при помощи неявной схемы Адамса 3-го порядка, демонстрирующий возможности численного решения поставленной задачи с заданной степенью точности. Готовый программный продукт может найти широкое применение при решении многих ...

Скачать
51291
3
14

... силы взаимодействия между рыбой и рабочими органами машин, поскольку изменяется площадь контакта, обусловливающая силы трения. До настоящего времени структурно-механические характеристики в основном оцениваются органолептическим методом. Рыбу сдавливают пальцами и оценивают ее консистенцию. Не достатком такого метода оценки структурно-механических характеристик мышечной ткани является его ...

0 комментариев


Наверх