простейшая одночленная формула для первой производной является формулой первого порядка точности.
Рассмотрим теперь формулу второй производной (10). представим первое слагаемое правой части в виде разложения по формулу Тейлора в координате одного из узлов. Рассмотрим производную в окрестности точки х0, левой точки промежутка
[х0, x2]. Для этого нужно разложить f1; и f2 B окрестности точки х0
Б12.
1. Формулы прямой и обратной прогонки
МЕТОД ПРОГОНКИ
и эти равенства дают 2N соотношений для определения неизвестных констант аi, и bi
и в точности совпадает с интерполяционным многочленом Лагранжа в том случае, когда интерполяция проводится по двум точкам. Но в отличие от линейной лагранжевой интерполяции, с помощью сплайновой технологии можно аппроксимировать и неоднозначные кривые.
Интерполяция кубическими сплайнами
Если обозначить через у(х) отклонение от положения равновесия на каждом промежутке между узлами, то условием равновесия линейки является соотношение
Отсюда получаем, что кривая между соседними узлами имеет следующий вид: у(х) = А + Вх + Сх2 +Dx3.y”(x0)=0,y”(N)=0
Как и в случае формулы (9) мы предполагаем, что исходная функция обладает необходимым количеством производных для представления ее по формуле Тейлора. Тогда имеем
одночленная формула второй производной тоже является формулой первого порядка точности. Определим теперь погрешность в среднем узле промежутка х1. Для этого нужно разложить f0 и f2в окрестности этой точки:
в средней точке промежутка погрешность оказалась пропорциональной не первой степени h, а квадрату шага h2-второй порядок точности формулы (10) в средней точке интервала.
Двучленные формулы численного дифференцирования для получения формул численного дифференцирования используем первые два слагаемых.
ЕДИНСТВЕННОСТЬ РЕШЕНИЯ КУБИЧЕСКОЙ ИНТЕРПОЛЯЦИИ
КВАДРАТНЫЕ СПЛАЙНЫ
Оценим погрешность формулы (11), используя формулу Тейлора (по-прежнему предполагаем, что исходная функция обладает необходимым количеством производных):
Подставляем эти выражения в правую часть (11). Тогда имеем:
Таким образом, согласно формуле (11) первая производная по трехточечной схеме имеет вид:
Произведем вычисления по формуле (12)
Подстановка в правую часть (12) дает
и сравнение с правой частью приводит к следующему конечно-разностному выражению.
Из формул (13), (14) можем сделать несколько выводов. Во-первых, двучленные формулы первой производной (трехточечные схемы) имеют второй порядок точности. Справедлив более обший результат -нее двучленные формулы имеют второй порядок точности. Во-вторых, из сравнения (13) и (14) находим, что формула для центрально-разностной схемы вдвое точнее.
Чаще всего точная величина а неизвестна, и тогда вводят предельную абсолютную погрешность Δа, которую определяют как всякое число, не меньшее абсолютной погрешности
Погрешность результата арифметической операции с приближенными числами.
Т1. Абсолютная погрешность алгебраической суммы нескольких приближенных чисел не превышает суммы абсолютных погрешностей этих чисел.
Доказательство. Пусть - приближенные числа. Рассмотрим алгебраическую сумму y = ±x1±x2±…±xn (5)
Формулу (5) будем рассматривать как функцию n переменных. Очевидно, что. Используя (4), получим:
Δy = ±Δx1 ± Δх2 ±… ± Δхn à |Δy| ≤ |Δx1 |+ |Δх2 |+… +| Δхn|
Что и требовалось доказать.
Теорема Относительная погрешность произведения нескольких приближенных чисел, отличных от нуля, не превышает суммы относительных погрешностей этих чисел.
Б13.
1. Порядок аппроксимации производных на сетке
Одночленные формулы дифференцирования. Согласно одночленным формулам (7) и (8), первые две производные равны:
Вычисление первой производной проводится на отрезке [х0, x1], второй - на [х0, x2]. Рассмотрим оценку точности вычислений основанную, на использовании ряда Тейлора (в общем случае формулы Тейлора) для представления исходной функции.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.