Материал из MachineLearning.
![L_{2,i}'(x)=\frac {1}{\bar(h_i)}[(x-x_{i-\frac{1}{2}}) \frac{u_{i+1}-u_i}{h_{i+1}} + (x_{i+\frac{1}{2}}-x) \frac{u_i-u_{i-1}}{h_i}]](../mimetex/? L_{2,i}'(x)=\frac {1}{\bar(h_i)}[(x-x_{i-\frac{1}{2}}) \frac{u_{i+1}-u_i}{h_{i+1}} + (x_{i+\frac{1}{2}}-x) \frac{u_i-u_{i-1}}{h_i}])
)
=\frac{1}{2}(\frac{h_i}{\bar{h_i}}\frac{u_{i+1}-u_i}{h_{i+1}}+\frac{h_{i+1}}{\bar{h_i}}\frac{u_i-u_{i-1}}{h_i}))
=u_{\dot{x},i})
)
≈=\frac{1}{\bar{h_i}}(\frac{u_{i+1}-u_i}{h_{i+1}}- \frac{u_i-u_{i-1}}{h_i}))
Введение
Постановка математической задачи
Численное дифференцирование применяется, если функцию )
трудно или невозможно продифференцировать аналитически - например, если она задана таблицей. Оно нужно также при решении дифференциальных уравнений при помощи разностных методов.
Изложение метода
При численном дифференцировании функцию аппроксимируют легко вычисляемой функцией )
и приближенно полагают
=\varphi'(x))
. При этом можно использовать различные способы аппроксимации.
Интерполирование полиномами Лагранжа
Рассмотрим неравномерную сетку 
и обозначим за 
, 
шаги этой сетки. В качества примера получим формулы численного дифференцирования, основанные на использовании многочлена Лагранжа )
, построенного для функции )
по трем точкам 
.
Многочлен имеет вид
=\frac {(x-x_i)(x-x_{i+1})}{h_i(h_i+h_{i+1})}u_{i-1}-\frac {(x-x_{i-1})(x-x_{i+1})}{h_ih_{i+1}}u_i+\frac {(x-x_{i-1})(x-x_i)}{h_{i+1}(h_i+h_{i+1})}u_{i+1})
Отсюда получим
=\frac{(2x-x_i-x_{i+1})}{h_i(h_i+h_{i+1})}u_{i-1}-\frac {(2x-x_{i-1}-x_{i+1})}{h_ih_{i+1}}u_i+\frac {(2x-x_{i-1}-x_i)}{h_{i+1}(h_i+h_{i+1})}u_{i+1})
Это выражение можно принять за приближенное значение )
в любой точке 
∈ ![[x_{i-1},x_{i+1}]](../mimetex/?[x_{i-1},x_{i+1}])
.
Его удобнее записать в виде
![L_{2,i}'(x)=\frac {1}{\bar{h_i}}[(x-<tex>x_{i-\frac{1}{2}}) \frac{u_{i+1}-u_i}{h_{i+1}} + (x_{i+\frac{1}{2}}-x) \frac{u_i-u_{i-1}}{h_i}]](../mimetex/? L_{2,i}'(x)=\frac {1}{\bar{h_i}}[(x-<tex>x_{i-\frac{1}{2}}) \frac{u_{i+1}-u_i}{h_{i+1}} + (x_{i+\frac{1}{2}}-x) \frac{u_i-u_{i-1}}{h_i}])
, где
, .
В частности, при 
получим
,
И если сетка равномерна, 
, то приходим к центральной разностной производной, .
При использовании интерполяционного многочлена первой степени точно таким образом можно получить односторонние разностные производные 
и 
.
Далее вычисляя вторую производную многочлена , получим приближенное выражение для при ∈:
≈
На равномерной сетке это выражение совпадает со второй разностной производной 
. Ясно, что для приближенного вычисления дальнейших производных уже недостаточно многочлена , надо привлекать многочлены более высокого порядка и тем самым увеличивать число узлов, участвующих в аппроксимации.
