Метод Рунге-Кутта второго порядка точности

Вычисления по этому методу осуществляются в два этапа. На первом этапе по схеме Эйлера находится промежуточное значение

(5.5)

На втором этапе находится значение y_i+1 по схеме

(5.6)

где a> 0, s> 0 - параметры. Подставляя из (5.5) в (5.6), имеем

(5.7)

Нетрудно проверить (разложение по формуле Тейлора), что схема (5.7) имеет второй порядок аппроксимации при условии sa=1/2. Частные случаи разностной схемы (5.7)

(5.8)

Эта разностная схема носит название предиктор-корректор, или счет-пересчет. Первая схема из (5.8) - схема Эйлера с шагом (предиктор), вторая - схема со значением на полушаге (корректор)

Метод Рунге-Кутта четвертого порядка точности

Используется схема

(5.9)

где k₁, k₂, k₃, k₄ - поправки, вычисляемые по формулам

(5.10)

При определении y_i+1 по заданному y_i необходимо четыре раза вычислять правую часть (5.9) в следующей последовательности: k₁, k₂, k₃, k₄. Если предположить достаточную гладкость u(t) (непрерывную дифференцируемость вплоть до производных 4-го порядка) и разложить u_i+1, k₁, k₂, k₃, k₄ в окрестности t=t_i, нетрудно показать, что невязка y = 0(t⁴), т.е. разностная схема (5.9) имеет 4-й порядок аппроксимации.

Пример выполнения лабораторной работы №5

Требуется решить задачу Коши при помощи численных методов

Вводим функцию f(t, y)

Задаем шаг, количество шагов по времени и начальное условие

Выводим длину расчетного временного интервала

Задаем точное значение решения (считается аналитически)

Метод Эйлера для решения задачи Коши.

Вводим функцию, реализующую алгоритм метода Эйлера

Вызываем данную функцию

Строим графики точного решения и приближенного рассчитанного методом Эйлера

Выводим значение погрешности

Метод Рунге–Кутта второго порядка точности для решения задачи Коши.

Вводим функцию, реализующую алгоритм метода Рунге-Кутта второго порядка точности

Вызываем данную функцию

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

Выводим значение погрешности

Метод Рунге–Кутта четвертого порядка точности для решения задачи Коши.

Вводим функцию, реализующую алгоритм метода Рунге-Кутта четвертого порядка точности

Вызываем данную функцию

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

Выводим значение погрешности

Варианты заданий к лабораторной работе №5

Решить задачу Коши, используя методы Эйлера, Рунге – Кутта второго и четвертого порядков точности.

10.

Содержание отчета

Отчет должен содержать:

1. титульный лист;

2. постановку задачи (согласно варианту);

3. точное решение задачи;

4. краткое описание методов решения задачи Коши;

5. программную реализацию данных методов;

6. выводы о проделанной работе.

Методы приближения функций

Интерполяция в общем понимании - это нахождение промежуточных значений величины по некоторым известным ее значениям.

Интерполяция применяется во многих задачах, связанных с вычислениями, например:

1) восстановить функцию y(x) для всех значений x на отрезке [a, b], если известны ее значения в некотором конечном числе точек этого отрезка;

2) заменить функцию близким выражением, удобным для проведения вычислений, если исходная функция представлена выражением, трудным для вычислений;

3) выполнить субтабулирование, т.е. сгущение таблицы значений;

4) решить задачу обратного интерполирования, т.е. по заданной таблице y_i=y(x_i) значений функции найти x как функцию от y.

В результате возникает следующая математическая задача: пусть на отрезке [a, b] задана сетка , в ее узлах заданы значения функций y(x), равные y(x₀) = y₀, y(x₁) = y₁, …, y(x_i) = y_i, …, y(x_n) = y_n. Требуется построить интерполянту, т.е. функцию j(x), совпадающую с функцией y(x) в узлах сетки

j(x_i) = y_i, i = 0, 1, …, n.

(6.1)

Интерполянта строится, как правило, в виде линейной комбинации некоторых линейно независимых на отрезке [a, b] элементарных функций

{j_k(x)}, k = 0, 1, …, n;

(6.2)

где c_k -неизвестные коэффициенты. Из условия (6.1) для определения этих коэффициентов получимсистему (n+1) уравнений

(6.3)

Система функций {j_k(x)}, k = 0, 1, …, n такова, что при любом выборе сетки узлов отличен от нуля определитель системы (6.3)

Поэтому из системы (6.3) однозначно определяются коэффициенты с_к.

В качестве функции j(x) в дальнейшем будем брать алгебраический многочлен степени n, т.е. j(x) = P_n(x). В этом случае интерполяция называется алгебраической. Алгебраическая интерполяция является наиболее распространенной, т.к. многочлены легко вычисляются, их легко дифференцировать и интегрировать. В этом случае в качестве системы линейно независимых функций {j_k(x)} обычно выбирают следующую систему: j_k(x) = x^k, k = 0, 1, …, n.

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

1 2 3 456