Аппроксимация функциональной зависимостью вида

y = a₀ . ( 2.7 )

Прологарифмировав (2.7) и вводя обозначение имеем:

Y = a₁ x + A₀ ,

где Y = ln y ; A₀ = ln a₀ ; X = ln x .

Имея результаты опыта в виде n пар

, x1	x2	…	xn
y1	y2	…	yn

перейдем к парам

. X1	X2	…	Xn
y1	y2	…	yn

Здесь Y = ln y , X = ln x и по (2.3) находим a₁ и A₀ .

В заключение рассмотрим в качестве эмпирической функции многочлен

j(x) = a₀ + a₁ x + … + a_m x^m .

Тогда формула для определения суммы квадратов отклонения будет иметь вид

F = ² ® min.

Для составления системы уравнений найдем частные производные функции F по (a₀ , a₁, …, a_m ):

(a ₀ + a₁ x _i + … + a_m- y_i );

(a ₀ + a₁ x _i + … + a_m- y_i ) x_i ;

(a ₀ + a₁ x _i + … + a_m- y_i ).

Приравнивая эти выражения нулю в соответствии с уравнением (2.4) и группируя коэффициенты при неизвестных a₀ , a₁, …, a_m , получим следующую систему уравнений :

(2.7)

…

Решая эту систему линейных уравнений, находим коэффициенты a₀ , a₁, …, a_m , которые являются искомыми параметрами эмпирической формулы.

Систему (2.7) можно записать в более компактной форме

b₀₀ a₀ + b₀₁ a₁ + … + b_0m a_m = c₀ ;

b₁₀ a₀ + b₁₁ a₁ + … + b_1m a_m = c₁ ;

…

b_m0 a₀ + b_m1 a₁ + … + b_mm a_m = c_m ,

где b_ki = ; c_k = ; k, i = 0, m .

Относительные погрешности аппроксимации в заданных точках можно оценить по формуле

.

 

Интерполирование

В результате экспериментальных исследований часто получают таблицу значений некоторой функции f(x) при фиксированных значениях аргумента x_i , т.е. f(x_i ), i = 0, n. Аналитическая зависимость между x_i и f(x_i ) неизвестна, что позволяет вычислить значение функции f(x_i ) в промежуточных точках, отличающихся от экспериментальных точек x_i , i = 0, n. Для отыскания этих значений строят аппроксимирующую (приближенную) функцию j(x), расчеты по которой совпадают либо в некотором смысле приближаются к экспериментально наблюдаемым значениям. Построение функции j(x) называется интерполированием. К интерполированию прибегают и в случае, когда аналитический вид функции f(x) известен, но для получения ее значений необходимо провести большой объем вычислений. Замена функции f(x) приближенной формулой j(x) позволяет упростить вычисления.

Пусть y = f(x) существует для любой точки отрезка [a, b ], тогда ее значения известны только в отдельных точках x₀ , x₁ , … , x_n этого отрезка. Пусть x^* некоторая точка из [a, b], и нужно найти неизвестное значение y^* = f(x^*) по известным значениям y₀ = f(x₀) , y₁ = f(x₁) , …, y_n = f(x_n). Такая задача называется задачей интерполирования функции y = f(x).

Для решения этой задачи используют алгебраический многочлен n-й степени P_n( x ), принимающий в точках x₀ , x₁ , … , x_n те же значения, что и функция f(x), т.е.

 

Данный многочлен называется интерполяционной формулой Лагранжа.

.

Его так же можно записать в следующем виде:

Это алгебраический многочлен n-й степени, удовлетворяющий всем табличным значениям функции, т.е. если в эту формулу подставить значение x = a , то получим значение j (a) = f (a), алгоритм нахождения интерполяционного многочлена представлен на рис. 2.4.

Рис. 2.4