Линейная интерполяция.

Итак, пусть мы имеем функцию, заданную таблично. Решая задачу интерполяции, найдем в таблице два соседних значения аргумента (обозначим их х_k и x_k+1), между которыми лежит заданное значение х (х_k <x<x_k+1), пусть y_k=f(x_k) и y_k+1=f(x_k+1) – соответствующие им значения функции. Будем считать, что в промежутке (х_k , x_k+1) данную функцию с достаточной степенью точности можно заменить линейной функцией, т.е. дугу графика функции можно заменить стягивающей ее хордой (рис.2). Такая замена называется линейной интерполяцией.

Рис. 2

Уравнение прямой, проходящей, через точки (х_k , y_k) и (х_k+1 , y_k+1), имеет следующий вид:

или в более привычной форме уравнения с угловым коэффициентом:

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

Будем предполагать также, что вторая производная функции f(x) на рассматриваемом участке непрерывна и удовлетворяет неравенству

,

где .

Используя аппарат математического анализа можно доказать [27], что для любого х из интервала (x_k, x_k+1) оценка погрешности линейной интерполяции будет иметь следующий вид:

Заметим, что вторая производная функции f(x) имеет конкретный механический смысл. Если f(x) описывает закон движения материальной точки, то вторая производная этой функции задает ускорение этой точки в момент времени х. Факт существования ограничения на ускорение (ограниченность второй производной) с физической точки зрения означает, что процесс описываемый функцией f(x) протекает относительно равномерно и функция изменяется не очень быстро. Таковой, например, будет функция, задающая изменение суточной температуры воздуха от времени. На практике именно этим критерием "плавности" скорости изменения процесса можно вполне воспользоваться для ответа на вопрос об обоснованности применения линейной интерполяции.

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

и условие применимости линейной интерполяции запишется в виде неравенства:

(2)

Шаг и точность таблицы обычно стараются согласовать так, чтобы условие (2) было выполнено.

Бывает, однако, что для выполнения этого условия требуется выбирать слишком малый шаг. В таком случае не считаются с этим условием, а для отыскания промежуточных значений функции пользуются более сложной квадратичной интерполяцией или другими приемами [17].