Многоточечные аппроксимации производной
До сих пор для расчета разностных производных применялись
двухточечные конечные разности. Однако, порядок
аппроксимации можно существенно улучшить, если задействовать большее количество точек, в которых
рассчитывается f(x
0). Объем компьютерных вычислений при этом, конечно, возрастет, но и точность
численного дифференцирования улучшится.
Чтобы получить многоточечные разностные аппроксимации производной,
следует действовать точно так же, как и раньше (см. разд.
о погрешностях).
Вдобавок к трем точкам: x=x
0-
D,
x=x
0 и x=x
0+
D выберем
дополнительные точки аргумента x, например: x=x
0-2
D и x=x
0+2
D . Подставим эти значения x в разложение Тейлора
с тем, чтобы можно выразить через f
(x0-2D), f (x0+2D
) и слагаемые, содержащие множитель D2 и D3.
Полученная после приведения подобных слагаемых формула будет выглядеть следующим
образом:
Это - пятиточечная
аппроксимация 1-й производной. То, что в ней присутствуют значения f(x) лишь в
4-х точках не должно смущать читателя. Просто оказалось, что f(x0) входит в нее с множителем
0.
Аналогично можно получить и семиточечную аппроксимацию 1-й производной:
Для краткости, в данной формуле мы обозначили
символом f
±D значение функции f(x) в точках, отстоящих влево или
вправо от x0 на расстояние, кратное
D.
Аналогичным образом можно получить формулы и для
производных высших порядков. Следует записать
разложение Тейлора с точностью до слагаемых, содержащих нужную производную, а затем подставить в него
соответствующие разностные формулы для низших производных. Однако надо учесть, что следует обеспечить
надлежащий порядок точности, подставляя в ряд Тейлора соответствующие многоточечные аппроксимации низших
производных.