Численное решение задачи Коши
Контрольная работа - Математика и статистика
Другие контрольные работы по предмету Математика и статистика
Оглавление
Аннотация.
Оглавление.
. Задача № 1 (1.4)
.1Постановка задачи.
.2Исходные данные.
.3Решение поставленной задачи.
.Задача № 2 (2.2)
.1Постановка задачи.
.2Исходные данные.
.3Решение поставленной задачи.
.Задача № 3 (6.2)
.1Постановка задачи.
.2Исходные данные.
.3Решение поставленной задачи.
Заключение.
Список литературы.
1. Задача № 1 (1.4)
1.1Постановка задачи
Найти приближенное решение задачи Коши для обыкновенного дифференциального уравнения (ОДУ) 1 порядка
(1)
и оценить погрешность решения задачи.
Порядок решения задачи:
1. Задать исходные данные: функцию f правой части, начальное значение .
. Используя функцию eyler (см. ПРИЛОЖЕНИЕ B), найти приближенное решение задачи Коши с шагом h=0.1 по явному методу Эйлера.
. Используя встроенную функцию rkfixed пакета MATHCAD, найти приближенное решение задачи Коши с шагом h=0.1 по методу Рунге-Кутты 4 порядка точности (см. ПРИЛОЖЕНИЕ B).
. Найти решение задачи Коши аналитически.
. Построить таблицы значений приближенных и точного решений. На одном чертеже построить графики приближенных и точного решений.
. Оценить погрешность приближенных решений двумя способами:) по формуле ; здесьи - значения точного и приближенного решений в узлах сетки , i=1,..N;) по правилу Рунге (по правилу двойного пересчета) (см. ПРИЛОЖЕНИЕ C).
7. Выяснить, при каком значении шага h=h* решение, полученное по методу Эйлера, будет иметь такую же погрешность (см. п. 6а), как решение, полученное с помощью метода Рунге-Кутты с шагом h=0.1.
УКАЗАНИЕ. В п. 7 рекомендуется провести серию вычислений решения по методу Эйлера, дробя шаг h пополам.
1.2Исходные данные
Nf(t,y)t0Ty01.4011
.3Решение поставленной задачи
1)Задача Коши: y(t)=, t0=0, T=1, y0=1.
Исходные данные:
Начальное значение:
Концы отрезка:
Шаг сетки:
Число узлов сетки:
2)Функция, реализующая явный метод Эйлера, возвращает вектор решения:
Входные параметры:
f - функция правой части;
y0 - начальное значение;- начальная точка отрезка;- шаг сетки;- число узлов сетки.
3)Приближенное решение задачи Коши с шагом h=0.1 по методу Рунге-Кутты 4 порядка точности с помощью встроенной функции rkfixed пакета MATHCAD.
Функция rkfixed возвращает матрицу, первый столбец которой содержит узлы сетки, а второй - приближенное решение в этих узлах.
4)Аналитическое решение задачи:
,
,
,
= ,
,
По методу вариации произвольной постоянной заменим постоянную С на функцию C(t) и решим неоднородное уравнение:
Подставляем в исходное уравнение:
=,
,
Решение в MathCad:
5)Решения, полученные различными способами:
Метод Эйлера:
Метод Рунге-Кутты:
Точное решение:
Графики приближенных и точного решений:
6)Рассчитаем погрешность полученных приближенных решений:
Погрешность метода Эйлера:
Вычисление погрешности по правилу Рунге:
Вычисление приближенных решений с шагом h/2:
Вычисление погрешностей:
Значение погрешностей:
7)Проведём серию вычислений решения по методу Эйлера, дробя шаг h пополам.
Первая итерация:
Вторая итерация:
Третья итерация:
И т.д.
Девятая итерация:
При значении шага hd=h/1024=0.1/1024=0,0000977 решение, полученное по методу Эйлера, будет иметь примерно такую же погрешность, как и решение, полученное с помощью метода Рунге-Кутты с шагом h=0.1.
Погрешность решения по методу Эйлера с шагом hd=0,0000977 :
Погрешность решения по методу Рунге-Кутты с шагом h=0.1 :
Вывод: явный метод Эйлера - это численный метод 1-го порядка точности. Метод Рунге-Кутты - это метод 4-го порядка точности. Это означает, что при одном и том же значении шага, метод Рунге-Кутты даёт более точное значение. Поэтому погрешности методов сильно (на несколько порядков) отличаются. В рассмотренном выше примере с помощью метода Рунге-Кутты было получено решение, которое совпадает с решением, полученным аналитическим путём.
2.Задача № 2 (2.2)
2.1Постановка задачи
Задача Коши для ОДУ 2 порядка
,
описывает движение груза массы m, подвешенного к концу пружины. Здесь x(t) - смещение груза от положения равновесия, H - константа, характеризующая силу сопротивления среды, k -коэффициент упругости пружины, f(t) - внешняя сила. Начальные условия: - смещение груза в начальный момент времени t=0, - скорость груза в начальный момент времени. Промоделировать движение груза на временном отрезке [0,T] при заданных в индивидуальном варианте трех наборах (I, II, III) значений параметров задачи. Для каждого набора по найденной таблице (или графику) решения задачи определить максимальное и минимальное значения функции x(t) и моменты вр