- •Таганрогский государственный радиотехнический университет
- •Введение
- •1 Численное интегрирование
- •1.1 Основные методы численного интегрирования
- •1.2 Пример выполнения лабораторной работы №1
- •1.3 Варианты заданий к лабораторной работе №1
- •2 Методы решения нелинейных уравнений
- •2.1 Метод половинного деления
- •2.2 Метод хорд (метод линейной интерполяции)
- •2.3 Метод секущих
- •2.4 Метод Ньютона
- •2.5 Пример выполнения лабораторной работы №2
- •2.6 Варианты заданий к лабораторной работе №2
- •3 Прямые методы решения систем линейных алгебраических уравнений
- •3.1 Метод Гаусса для решения систем линейных алгебраических уравнений.
- •3.2. Алгоритм lu-разложения.
- •3.3. Метод прогонки.
- •3.4 Пример выполнения лабораторной работы №3
- •Вводим функцию, реализующую алгоритм прогонки
- •3.5 Варианты заданий к лабораторной работе №3
- •4 Итерационные методы решения систем линейных алгебраических уравнений
- •4.1. Итерационные методы Якоби и Зейделя.
- •4.2. Каноническая форма итерационных методов.
- •4.3. Вариационно-итерационные методы решения слау.
- •4.4 Пример выполнения лабораторной работы №4
- •Задаем матрицу коэффициентов и столбец свободных членов
- •Вводим функцию, реализующую алгоритм метода Зейделя
- •4.5 Варианты заданий к лабораторной работе №4
- •5 Методы решения задачи Коши
- •5.1. Метод Эйлера.
- •5.2. Метод Рунге-Кутта.
- •Пример выполнения лабораторной работы №5
- •5.4 Варианты заданий к лабораторной работе №5
- •6 Методы приближения функций
- •6.1. Интерполяционный полином Лагранжа и Ньютона.
- •6.2 Интерполяционный кубический сплайн.
- •6.3 Понятие о методе наименьших квадратов.
- •6.4 Интерполяционный тригонометрический полином
- •6.5 Пример выполнения лабораторной работы №6
- •Интерполяционный полином Лагранжа.
- •6.6 Варианты заданий к лабораторной работе №6
5.2. Метод Рунге-Кутта.
Повышение порядка точности осуществляется путем усложнения разностной схемы. На практике широко распространенными являются разностные схемы Рунге-Кутта второго и четвертого порядка точности.
Метод Рунге-Кутта второго порядка точности
Вычисления по этому методу осуществляются в два этапа. На первом этапе по схеме Эйлера находится промежуточное значение
. |
(5.5) |
На втором этапе находится значение yi+1 по схеме
, |
(5.6) |
где >0, >0 - параметры. Подставляя из (5.5) в (5.6), имеем
. |
(5.7) |
Нетрудно проверить (разложение по формуле Тейлора), что схема (5.7) имеет второй порядок аппроксимации при условии =1/2. Частные случаи разностной схемы (5.7)
. |
(5.8) |
Эта разностная схема носит название предиктор-корректор, или счет-пересчет. Первая схема из (5.8) - схема Эйлера с шагом (предиктор), вторая - схема со значением на полушаге (корректор)
Метод Рунге-Кутта четвертого порядка точности
Используется схема
. |
(5.9) |
где k1, k2, k3, k4 - поправки, вычисляемые по формулам
. |
(5.10) |
При определении yi+1 по заданному yi необходимо четыре раза вычислять правую часть (5.9) в следующей последовательности: k1, k2, k3, k4. Если предположить достаточную гладкость u(t) (непрерывную дифференцируемость вплоть до производных 4-го порядка) и разложить ui+1, k1, k2, k3, k4 в окрестности t=ti, нетрудно показать, что невязка = 0(4), т.е. разностная схема (5.9) имеет 4-й порядок аппроксимации.
-
Пример выполнения лабораторной работы №5
Требуется решить задачу Коши при помощи численных методов
Вводим функцию f(t,y)
Задаем шаг, количество шагов по времени и начальное условие
Выводим длину расчетного временного интервала
Задаем точное значение решения (считается аналитически)
Метод Эйлера для решения задачи Коши.
Вводим функцию, реализующую алгоритм метода Эйлера
Вызываем данную функцию
Строим графики точного решения и приближенного рассчитанного методом Эйлера
Выводим значение погрешности
Метод Рунге–Кутта второго порядка точности для решения задачи Коши.
Вводим функцию, реализующую алгоритм метода Рунге-Кутта второго порядка точности
Вызываем данную функцию
Строим графики точного решения и приближенного рассчитанного методом Рунге-Кутта второго порядка точности
Выводим значение погрешности
Метод Рунге–Кутта четвертого порядка точности для решения задачи Коши.
Вводим функцию, реализующую алгоритм метода Рунге-Кутта четвертого порядка точности
Вызываем данную функцию
Строим графики точного решения и приближенного рассчитанного методом Рунге-Кутта четвертого порядка точности
Выводим значение погрешности
5.4 Варианты заданий к лабораторной работе №5
Решить задачу Коши, используя методы Эйлера, Рунге – Кутта второго и четвертого порядков точности.
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
Содержание отчета
Отчет должен содержать:
-
титульный лист;
-
постановку задачи (согласно варианту);
-
точное решение задачи;
-
краткое описание методов решения задачи Коши;
-
программную реализацию данных методов;
-
выводы о проделанной работе.