Правило Рунге практической оценки погрешности

5.5 Правило Рунге практической оценки погрешности

Оценки погрешности по формулам (5.4), (5.8) и (5.12) являются априорными. Они зависят от длины элементарного отрезка h, и при достаточно малом h справедливо приближенное равенство:

 I – I^h » Ch^k, (5.15)

где I^h приближенное значение интеграла, вычисленное по одной из формул (5.3), (5.5), (5.9), C ¹ 0 и k > 0 – величины, не зависящие от h.

Если уменьшить шаг h в два раза, то, в соответствии с (5.15) получим:

 I – I^h/2 » Ch^k » ( I – I^h). (5.16)

Непосредственное использование оценок погрешности (5.4), (5.8) и (5.12) неудобно, так как при этом требуется вычисление производных функции f(x). В вычислительной практике используются другие оценки.

Вычтем из равенства (5.15) равенство (5.16):

I^h/2 – I^h » Ch^k(2^k – 1). (5.17)

Учитывая приближенное равенство (5.16), получим следующее приближенное равенство:

I – I^h/2 »  . (5.18)

Приближенное равенство (5.18) дает апостериорную оценку погрешности. Вычисление этой оценки называется правилом Рунге. Правило Рунге – это эмпирический способ оценки погрешности, основанный на сравнении результатов вычислений , проводимых с разными шагами h.

Для формул прямоугольников и трапеций k = 2, а для формулы Симпсона k = 4. Поэтому для этих формул приближенное равенство (5.18) принимает вид:

I – I_пр» , (5.19)

I – I_тр» , (5.20)

I – I_С » . (5.21)

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

Пример 5.4.

Найдем значение интеграла  с точностью e = 10^-4, используя формулу трапеций и применяя вышеизложенную процедуру дробления шага. В примере 5.2 было получено значение I при h₁ = 0.1, I^h =0.74621079. Уменьшим шаг вдвое: h₂ = 0.05 и вычислим I= 0.74667084, e₂  = ( I- I) = (0.74667084 – 0.74621079) » 1.5×10^-4. Так как |e₂| > e, то снова дробим шаг: h₃ = 0.025, вычисляем I= 0.74678581, e₂  = ( I- I) = (0.74678581 – 0.74667084) » 4×10^-5. Поскольку |e₃| < e, требуемая точность достигнута и I » 0.7468 ± 0.0001.

Тема 6. Численное решение дифференциальных уравнений

6.1 Постановка задачи Коши

Известно, что обыкновенное дифференциальное уравнение первого порядка имеет вид:

y' (t) = f(t, y(t)). (6.1)

Решением уравнения (6.1) является дифференцируемая функция y(t), которая при подстановке в уравнение (6.1) обращает его в тождество. На рис. 6.1 приведен график решения дифференциального уравнения (6.1). График решения дифференциального уравнения называется интегральной кривой.

Рис. 6.1

Производную y'(t) в каждой точке (t, y) можно геометрически интерпретировать как тангенс угла a наклона касательной к графику решения, проходящего через эту точку, т е.: k = tga = f(t, y).

Уравнение (6.1) определяет целое семейство решений. Чтобы выделить одно решение, задают начальное условие:

y(t₀ ) = y₀, (6.2)

где t₀ – некоторое заданное значение аргумента t, а y₀ – начальное значение функции.

Задача Коши заключается в отыскании функции y = y(t), удовлетворяющей уравнению (6.1) и начальному условию (6.2). Обычно определяют решение задачи Коши на отрезке, расположенном справа от начального значения t₀, т. е. для t Î [t₀, T].

Разрешимость задачи Коши определяет следующая теорема.

Теорема 6.1. Пусть функция f(t, y) определена и непрерывна при t₀ £ t £ T, -¥ < y < ¥ и удовлетворяет условию Липшица:

| f(t, y₁) – f(t, y₂)| £ L| y₁ – y₂|,

где L некоторая постоянная, а y₁ , y₂  – произвольные значения.

Тогда для каждого начального значения y₀ существует единственное решение y(t) задачи Коши для t Î [t₀, T].

Даже для простых дифференциальных уравнений первого порядка не всегда удается получить аналитическое решение. Поэтому большое значение имеют численные методы решения. Численные методы позволяют определить приближенные значения искомого решения y(t) на некоторой выбранной сетке значений аргумента t_i, (i = 0, 1, …). Точки t_i называются узлами сетки, а величина h_i = t_i+1 – t_i – шагом сетки. Часто рассматривают равномерные сетки, для которых шаг h_i постоянен, h_i = h = . При этом решение получается в виде таблицы, в которой каждому узлу сетки t_i соответствуют приближенные значения функции y(t) в узлах сетки y_i » y(t_i).

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

Сходимость численных методов решения задачи Коши. Пусть y(t) – решение задачи Коши. Назовем глобальной погрешностью (или просто погрешностью) численного метода функцию e_i = y(t_i) – y_i , заданную в узлах сетки t_i. В качестве абсолютной погрешности примем величину R = | y(t_i) – y_i|

Численный метод решения задачи Коши называется сходящимся, если для него R ® 0 при h ® 0. Говорят, что метод имеет p-ый порядок точности, если для погрешности справедлива оценка R £ Ch^p, p > 0, C – константа, C ¹ 0.

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

Лекции по вычислительной математике

Скачать

10356

0

0

... . Рассмотрение метода ветвей и границ для решения задачи о коммивояжере удобнее всего проводить на фоне конкретного примера. Пользуясь введенными здесь обозначениями, мы проводим это описание в следующей лекции. Введем некоторые термины. Пусть имеется некоторая чис- ловая матрица. Привести строку этой матрицы означает выде-лить в строке минимальный элемент (его называют константой приведения) ...

Основные задачи вычислительной математики

Скачать

7721

1

1

... если  - предельная абсолютная погрешность приближённого числа , то  (1.2) отсюда следует, что  (1.3) Значение предельной абсолютной погрешности, обычно, подбирается интуитивно по смыслу задачи. Пример 2: Определить предельную абсолютную погрешность числа , заменяющего число , точное значение которого нам неизвестно. Так как мы знаем, что , ...

Выдающиеся личности в истории вычислительной техники. Августа Ада Лавлейс

Скачать

37333

0

0

... удивили меня…, хоть речь идёт обо мне самой. Они действительно написаны прекрасным стилем, который превосходит стиль самого очерка" /2/. 2.3. Рождение первенца и критическое перенапряжение Августа Ада Лавлейс работает с большим напряжением. В письмах к Бэббиджу она неоднократно жалуется на утомление, болезни, плохое самочувствие. Наконец, 6 августа Бэббидж отсылает Аде свои последние замечания ...

Роль женщин в развитии вычислительной техники

Скачать

54819

0

0

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