Справочное руководство по небесной механике и астродинамике -
Скачать (прямая ссылка):
того же столбца. Например,
/! =4 (/L,,+/},), n.=i(n+tl) (7.1.03)
и т. д.
Если функция достаточно гладкая и ее значения вычислены достаточно точно,
а шаг таблицы мал, то разности довольно быстро убывают вместе с
увеличением порядка. Ошибки округления значений / приводят к тому, что
значения разностей, начиная с некоторого порядка, определяются именно
этими ошибками; разности теряют плавный ход и, как говорят, "скачут", так
что пользоваться ими нет смысла.
§ l.OZj
ГЛ. I. ИНТЕРПОЛИРОВАНИЕ И ПРИБЛИЖЕНИЕ ФУНКЦИИ
637
Если m-е разности, выраженные в единицах последнего десятичного знака
значений функции, превосходят по абсолютной величине 2m_1, то в значениях
функции обязательно имеется погрешность, не меньшая половины последнего
выписанного десятичного знака.
§ 1.02. Интерполяционные формулы
Пусть функция f(t) задана табл. 80. Интерполяционные формулы позволяют
вычислять (вообще, приближенно) значения этой функции для промежуточных
значений аргумента, не совпадающих с узлами таблицы. Эти формулы строятся
с помощью так называемых интерполяционных полиномов, которые представляют
собой полиномы по степеням независимой переменной, точно совпадающие с
f(t) в узлах. Для удобства вычислений интерполяционные формулы
располагаются по соответствующим разностям табл. 81.
Обозначим через t* значение аргумента, для которого требуется вычислить
значение функции f(t), через tQ - ближайший к t* узел таблицы, и положим
я "-нг2-. <7Л-04>
так что t* = to +' qh.
На практике используются чаще всего следующие интерполяционные формулы.
Формула Ньютона для интерполяции вперед
(применяется, если ?* > ^0. 0 < q f(t0 + qh)xsf0 + qf\/i +
, <?(<?- О г2 I , д(д- 1) ... (<?-" + 1) ta п , nt-v
+ 21 /1 "г ¦ • • ~г п! /п/2- (7.1.05)
Формула Ньютона для интерполяции назад (применяется, если t* < t0, у q <
о): f(t0 + qh)ttf0 + qfl-4l +
I Я (? + 1) (2 , I ?(?+!) ¦¦¦ Ь+Я-1) (П щ 1 пс\
+-----21---/-1+ +----------п\---------'-Ш2- (7Л-06)
Формула Стирлинга:
/(*о + '7Л)~/о + '7/о +
л. f2 л. I ч (ч2 - П (<?' ~ 22) (д2-(п- 1)") f2n_,
' 21 '0 ' - ' ' (2п- 1)! '0 I
+ , ¦(, ¦_ ¦- у, ^ {7Ш)
638
Ч. VII. ЧИСЛЕННЫЕ МЕТОДЫ
[§ 1.02
Формула Бесселя:
f%+1. (7.1.08)
(2/i + D!
?(?г - I2) п' ?(?2 ~ I2) (?' ~п2) f2,
31 I [~г ¦ ¦ • т (2п+ 1)1 11
Р(Р2 - I2) 12 г I Р(Рг- I2) ¦¦¦ (рг-п;)
3, /о -Г ... -г (2я + 1)!
-Пп] +
/'"] (р=1-</). (7.1.09)
Формула Ньютона (7.1.06) используется также для экстраполяции вперед, т.
е. для вычислений f(t) при значениях аргумента t, находящихся за
пределами табл. 80 справа.
Если через t0 обозначен последний справа узел табл. 80, то таблица
разностей вида 81 будет содержать лишь разности с отрицательными нижними
индексами. Для вычисления f(t*) при t* > to, t* - ta-\-qh используют
формулу (7.1.06), но в этом случае q = {t* - t0)/h > 0.
Точно так же можно экстраполировать f(t) для значений t слева от начала
табл. 80. Тогда применяется формула (7.1.05) при q <. 0.
Ценность формул Ньютона состоит в простоте структуры коэффициентов, в
возможности экстраполяции и в возможности непосредственной интерполяции
между узлами вблизи краев табл. 80, когда соответствующие разности с
нижними индексами 0, 1 или 1/2 отсутствуют.
Формулы Стирлинга и Бесселя, вообще, более точные, чем формулы Ньютона.
Формула Бесселя особенно удобна при интерполяции в точках, близких к
середине интервала.
Формула Эверетта часто применяется при субтабулировании, т. е. для
составления новой таблицы значений функции f(t) с более мелким шагом.
Формула Лагранжа:
П
(7.1.10)
§1.03] ГЛ. 1. ИНТЕРПОЛИРОВАНИЕ И ПРИБЛИЖЕНИЕ ФУНКЦИЙ 639
где
(о" (0 = (t - t0) ... (t- tn), ft = f (*,).
Это - самая общая интерполяционная формула, выражающая значение функции в
точке t непосредственно через значения этой функции в п + 1 узлах t0,
..., tn, причем расположение этих узлов произвольное, a t принадлежит
интервалу (^о, tn). Правая часть (7.1.10)- полином по t.
Если узлы равноотстоящие, с шагом А, то формула Лагранжа приобретает
следующий вид:
f (о "(- D"q'[q ~ 0'{q ~п)- ^ (- D1 - A- fh (7.1.11)
1=0
где С'п- биномиальные коэффициенты, q = {t - t0)/h.
Полином в правой части (7.1.10) называется полиномом Лагранжа, и он дает
приближенное представление функции f{t) на всем отрезке [fo, fn]- Он
может применяться не только для вычислений промежуточных значений f(t),
но и для различных операций с этой функцией (дифференцирование,
интегрирование и др.).
Наивысший порядок разностей, сохраняемых в интерполяционной формуле,
называется порядком этой формулы. Если таблица заданных значений функции
достаточно обширна, то принципиально возможно составить разности
достаточно высо* кого порядка и использовать соответственно
интерполяционные формулы такого же высокого порядка. На практике обычно
ограничиваются интерполяционными формулами не выше четвертого или пятого
порядков. Отбрасываемые при вычислениях члены не должны, вообще говоря,
превышать по абсолютной величине погрешность, которой обладают сами
