Методы интерполяции и аппроксимации
Скачать 93.46 Kb.
|
4.1.4.2. Вторая интерполяционная формула НьютонаДля нахождения значений функции в конце интервала интерполирования интерполяционный полином запишется в виде Pn(x) a0 a1(x xn) a2 (x xn)(x xn1) ... ... an(x xn)(x xn1)...(x x1). Коэффициенты а0, а1, ..., аnнаходятся из условия Pn(xi) = yi. Подставляя в (7) x= xn, найдем Pn(xn) yn a0 . (7) Для x=xn-1: P(x )=y =y+a(x – x), a yn yn1 yn1 . n Для x=xn-2: n-1 n-1 n 1 n-1 n 1 h h Pn(xn 2 ) yn 2 yn yn1 (xn 2 xn) a2 (xn 2 xn)(xn 2 xn1) h yn yn1 (2h) a2 2h2 h yn 2yn1 a2 2h2 ; a2 2 yn2 . 2!h2 Формула для нахождения всех коэффициентов запишется как: iy aiп1 . i!hi Подставив выражения для определения коэффициентов aiв формулу (7), получим вторую интерполяционную формулу Ньютона: Pn(x) yn yn1 h (x xn) 2 yn2 2!h2 (x xn)(x xn1 ) 3 yn3 3!h3 (x xn)(x xn1 )(x xn2 ) ... (8) |