Метод конечных разностей

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

Так, заменив производную в обыкновенном дифференциальном уравнении

$~u^\prime(x) = 3u(x) + 2$

на конечную разность

$~\frac {u(x+h)-u(x)} {h} \approx u^\prime(x)$,

получаем аппроксимированную форму (конечно-разностную схему)

$~u(x+h) = u(x) + h\cdot(3u(x)+2)$.

Последнее выражение носит название конечно-разностного уравнения, а его решение соответствует приближённому решению первоначального дифференциального уравнения.

Применение для интерполяции/экстраполяции

p(0)=2,0
	2,0−1,72=0,28
p(0,1)=1,72		0,28−0,24=0,04
	1,72−1,48=0,24
p(0,2)=1,48		0,24−0,20=0,04
	1,48−1,28=0,20
p(0,3)=1,28		0,20−0.16=0.04
	1,28−1,12=0,16
p(0,4)=1,12

Заметим, что значения в третьем столбце одинаковы. Это не случайность. Фактически, если мы начинаем таблицу для любого полинома степени n, колонка с номером n + 1 будет всегда содержать константу. Этот важный факт делает работоспособным указанный метод.

Мы составляли таблицу слева направо, но точно также её можно заполнить и справа налево, вычислив недостающие значения полинома.

Для вычисления p(0,5) мы воспользуемся значениями самой нижней диагонали. Начнём с самого нижнего значения в последней колонке, 0,04. Затем продолжим вторую колонку, вычтя 0,04 из 0,16 и получив значение 0,12. Таким же образом мы заполним первую колонку, вычитая из её нижнего значения 1,12 полученное нами на предыдущем шаге число 0,12 из второй колонки. p(0,5) будет равным 1,12-0,12 = 1,0. Для вычисления p(0,6) используется тот же самый алгоритм: берётся 0,04 из третьей колонки, вычитается из нижнего значения (теперь уже 0,12) во второй колонке, получившееся число 0,08 прописывается в нижнюю часть второй колонки и затем вычитается из нижнего значения в первой колонке (как мы помним — 1,0). Результат — 0,92 является значением p(0,6).

См. также

• Конечные разности
• Разностная схема
• Метод конечных разностей во временной области
• Интерполяция/Экстраполяция
• Метод конечных элементов
• Метод подвижных клеточных автоматов