birmaga.ru
добавить свой файл

1
ИНТЕРПОЛЯЦИЯ МНОГОЧЛЕНАМИ

Введение
Весьма часто ученым и инженерам приходиться сталкиваться с таблично заданными функциями. Такого типа функции возникают при работе с экспериментальными, статистическими данными. Однако нередко работа с табличными данными является слишком трудоемким занятием, особенно при больших объемах информации. Поэтому весьма полезным в таких случаях является приближение табличной функции аналитической. Это является полезным для компактного хранения информации, а также для дальнейшей практической и теоретической работы с полученной информацией.

Довольно часто табличные функции приближаются с помощью многочленов. Одной из разновидностей приближения табличной функции является задача интерполяции. Сформулируем одну из самых распространенных задач интерполяции функции одной переменной: Пусть дано n точек на координатной плоскости (x,y). Требуется построить многочлен степени не превосходящей n-1, проходящий через указанные точки. Нередко возникают задачи интерполяции, когда помимо значения функции в узловых точках задается еще значение первой производной (наклон функции). Существуют и другие формулировки задачи интерполяции. Мы, в основном, ограничимся указанными двумя задачами.
§ 1. Решение задачи интерполяции через решение СЛАУ
Многочлен степени n

(1)

Имеет n+1 коэффициент. Естественно полагать, что n+1 условие, наложенное на многочлен в общем виде, позволит однозначно определить коэффициенты. В частности, можно потребовать, чтобы многочлен проходил через п+1 точку (xi , yi) , i = 1, 2,..., n+1 с xi xj. To, что многочлен проходит через точки (xi , yi) означает выполнение условий


, i = 1, 2, …, n+1. (2)

Таким образом, мы получили систему n+1 линейных алгебраических уравнений (СЛАУ) для нахождения неизвестных коэффициентов a0 , a1, …, an . Определитель этой системы относительно не­известных ak -есть определитель Вандермонда

который не равен нулю, если xi xj для i j . Действительно, определитель есть функция от x1, x2, …, xn+1. Если считать его сначала функцией от xn+1, то он есть многочлен степени п и обращается в нуль всякий раз, когда xn+1 = xj (для j = 1, 2,n). Таким образом, f (x1, x2, ..., xn+1) содержит множители



Далее, рассматривая определитель как функцию от xn, мы видим точно так же, что существуют множители


Проведя аналогичные рассуждения для оставшихся узлов xn-1, xn-2, …, x1 , нетрудно заметить, что в нашем определителе встречаются вообще все множители




Произведение всех этих множителей есть многочлен степени



Но определитель также есть многочлен той же самой степени. Следовательно,

(3)

где С – некоторая константа которую еще надо найти. Чтобы определить С, рассмотрим выражение, которое получается при умножении главной диагонали



Раскрывая произведение, мы находим точно такой же член, и, сле­довательно, С = 1. Таким образом, определитель Вандермонда не равен нулю, если xi xj для i j

Возвращаясь к нашей главной задаче о нахождении многочлена по п+1 точке (xi , yi) , мы видим, что ее всегда можно решить и найти коэффициенты ak , решив систему (2) каким-либо известным способом. Подставив эти значения в общий вид многочлена (1), мы можем представить результат в виде
(4)

Впрочем, этот результат можно проверить, рассуждая следующим образом: во-первых, выражение (4) должно быть многочленом степени п по х, так как определитель можно разложить по элемен­там верхней строчки, и, во-вторых, он проходит через п+1 точку (xi , yi) так как подстановка этих значений в верхнюю строку делает две строки определителя одинаковыми.


Все изложенное можно подытожить, сказав, что если есть п+1 узловая точка функции, то можно найти многочлен степени п, кото­рый совпадает с функцией в узло­вых точках. И таким образом, мы можем использовать интерполяционный многочлен вместо функции в дальнейших аналитических процессах: интегрировании, дифференцировании, отыскании нулей и т. д.

Вместо требования, чтобы многочлен проходил через некоторые данные точки, можно потребовать, чтобы он в некоторых заданных точках имел данный наклон. Так, уравнение

определяет многочлен третьей степени по х, проходящий через (x1, y1) с наклоном y1 и через (x2, y2) с наклоном у2 (так как для дифференцирования этого определителя по х достаточно продиффе­ренцировать почленно его верхнюю строку).

Если в точке дано значение y то вовсе не обязательно в ней должно быть указано значение у, то же относится и к более высо­ким производным. Ограничением на выбор условий является требова­ние, чтобы минор члена у обязательно не был равен нулю. В противном случае многочлена может не существовать. Это означает, что должно быть дано по крайней мере одно значение функции yi, по крайней мере два условия на yi и уi’, три на yi, yi’, yi и т. д. вплоть до производной самого высокого порядка.

Только что сформулированное условие не является достаточным. Рассмотрим, например, три равноотстоящие точки, которые мы обо­значим -1, 0, 1. В каждой точке нам даны функция уi и вторая производная.


Минор при у равен

Таким образом, эти шесть условий не определяют многочлена пятой степени.

Упражнения



1. Напишите многочлен, проходящий через (x1, y1) с тангенсом угла наклона y1’ и второй производной у1” и имеющий тангенс угла наклона y2’ в x2.

2. Показать, что если использовать условия на у, у', у",..., у(n-1) в одной точке x1, то многочлен будет представлять оборванный ряд Тейлора.

3. Найти нетривиальный многочлен у(х) такой, что

y(-1) = y(0) = y(1) = y(-1) = y(0) = y(1) = 0.

4. Построить многочлен 2-й степени проходящий через узловые точки

а)(0,1), (1,3), (2,2); б)(0,-1), (1,2), (2,1); в)(0,1), (1,-1), (2,3).

5. Построить многочлен 3-й степени проходящий через узловые точки

а)(0,2), (1,3), (2,4), (3,1); б)(0,-1), (1,2), (2,1) (3,4); в)(0,2), (1,-1), (2,3), (3,5).
  1. Построить многочлен 3-й степени для следующей табличной функции:





а)

x

1

2

б)

x

1

2




y

0

2




y

1

0




y’

-1

0




y’

1

2


в)

x

1

2

г)

X

1

2




y

0

1




Y

1

2



y’


1

0




Y’

0

1


§ 2. Интерполяционные формулы Лагранжа и Эрмита.
Другой подход к задаче интерполяции - метод Лагранжа. Основ­ная идея этого метода состоит в том, чтобы прежде всего найти многочлен, который принимает значение 1 в одной узловой точке и 0 во всех других. Легко видеть, что функция


(где штрих у знака произведения означает отсутствие j-го множителя) является требуемым многочленом степени n; он равен 1, если x = xj, и 0, когда x=xi, i j.

Многочлен Lj (x) yj принимает значение yi в i-й узловой точке равен нулю во всех других узлах. Из этого следует, что

(5)

есть многочлен степени n, проходящий через n+1 точку (xi, yi).

Возникает вопрос: совпадает ли многочлен Лагранжа с многочле­ном (1), построенным в предыдущем, параграфе. Их разность - многочлен степени не выше п и она обращается в нуль в n+1 узловой точке x=xi, следовательно, она тождественно равна нулю. Таким образом, можно сделать следующее важное замечание: если дана п+1 узловая точка, то соответствующий многочлен степени п, про­ходящий через эти точки, однозначно (в пределах ошибок округле­ния) определен, независимо от того, как он строится и какие обо­значения использованы.


Интерполяционная формула Лагранжа неудобна для практического использования. Преобразования позволяют привести ее к несколько более удобному для вычисления виду. Рассмотрим частный случай, когда все yi = 1. Тогда у(х)=1 для всех х, так что тождественно



Мы можем теперь разделить правую часть формулы Лагранжа на это выражение, положив



после деления числителя и знаменателя на получим

. (6)
Формула Лагранжа вида (6) часто является более удобной для применения , чем (5).

Рассмотрим теперь задачу нахождения интерполяционного много­члена, который в каждой точке xi принимает значение yi и имеет производные yi’ , i=1, ... , п+1. Такой многочлен называется мно­гочленом Эрмита. Ясно, что он имеет степень 2n+1 (так как должны быть удовлетворены 2n+ 2 условия).

Для решения указанной задачи воспользуемся методом, являющимся естественным продолжением метода Лагранжа. Будем искать такие многочлены Hi(x) и hi(x) степени не выше 2п+1, что при i j

Hj(xi)=0, hj(xi)=0,

Hj’(xi)=0 , hj’(xi)=0,

Hj(xj)=1, hj(xj)=0,

Hj(xj)=0, hj’(xj)=1.


После того как эти многочлены будут найдены, решение поставлен­ной задачи можно будет записать в виде

.

Прежде всего построим hj(х) . Если hji) = 0 и hj’(хi) = 0, то hj) должен содержать множитель

(x - xi)2 (для i j).

Если hjj) = 0, но hj’(хj) 0, то существует простой множитель (x - xj). Чтобы получить hj’(хj) = 1, следует взять

Нетрудно видеть, что построенный многочлен hj(х) обладает требуемыми свойствами.

Построение многочлена Hi(x) почти такое же. Искомый многочлен должен содержать множители

(x - xi)2 (для i j),

но в то же время требуется, чтобы

Hj(xj)=1, Hj’(xj)=0.

Будем искать Hj(x) в виде

Условие Hj(xj)=1 требует, чтобы

axj + b =1.

Взяв производную от Hj и положив в ней x = xj получим


Таким образом,


и требуемая формула построена.

Упражнения


1. Указать способ получения интерполяционной фор­мулы, которая имеет значения yi, yiи yiв каждой точке хi.

2. С помощью интерполяционной формулы Лагранжа построить многочлен 2-й степени, проходящий через узловые точки:

а)(0,1), (1,3), (2,2); б)(0,-1), (1,2), (2,1); в)(0,1), (1,-1), (2,3).

3. С помощью интерполяционной формулы Лагранжа построить многочлен 3-й степени проходящий через узловые точки:

а)(0,2), (1,3), (2,4), (3,1); б)(0,-1), (1,2), (2,1) (3,4); в)(0,2), (1,-1), (2,3), (3,5).

4. С помощью интерполяционной формулы Эрмита построить многочлен 3-й степени для следующей табличной функции:



а)

x

1

2

б)

x


1

2




y

0

2




y

1

0




y’

-1

0




y’

1

2



в)

x

1

2

г)

x

1

2




y

0

1




y

1

2




y’

1

0




y’

0

1


§ 3. Интерполяционная формула Ньютона

Рассмотренные методы нахождения интерполяционного многочлена по n+1 точке подразумевали, что множество используемых узлов известно. Часто известным является лишь требуемая точность, а мно­жество узлов, которые следует использовать, определяется информа­цией о том, как вычисляется функция. Интерполяционная формула Ньютона, которая будет изложена, представляет собой просто другой способ написания интерполяционного многочлена. Она полезна, потому что число используемых узлов может быть легко увеличено или уменьшено без повторения всего вычисления.


Обозначим многочлен, проходящий через n+1 точку i, yi), i=1, 2, ..., n+1, через Pn = Рn (x). Можно написать



где Pn-1(x) – некоторый многочлен степени n-1. Ясно, что Pn(x1)=y1 и следует заботиться лишь об n точках (i=2, 3, …, n+1). Предыдущее уравнение может быть записано в виде


и следовательно, требуется, чтобы


Тем самым нужно, чтобы многочлен Pn-1(x) проходил через точки



Величины



называются разделенными разностями и обычно записываются в квад­ратных скобках.

На следующем шаге запишем

и потребуем, чтобы многочлен Pn-2(x) принимал значения разделенных разностей от разделенных разностей


Нетрудно видеть, что разделенные разности первого порядка не зависят от порядка аргументов в квадратных скобках. Покажем теперь, что это верно и для разделенных разностей второго порядка. Если начать с трех точек


(x1, y1), (x2, y2), (x3, y3),
то получим единственный многочлен второй степени, проходящий через три точки. Его можно записать так:



Если теперь взять точки в таком порядке: (xa, ya), (xb, yb), (xc, yc), то получим


Так как оба этих уравнения определяют один и тот же квадратный трехчлен, то коэффициенты при x2 должны быть одинаковы, так что два символа
[ x3, x2, x1 ] = [ xc, xb, xa ]

выражают одно и то же.

Вообще определяем


и в точности таким же способом показываем, что разности не зависят от порядка xi. Заметим, что знаменатель есть разность неповторяющихся значений x, взятых в том же порядке.

Один из способов построения таблицы необходимых значений следующий:

Таблица 1.

x1

y1*

















[x2, x1]*










x2

y2




[x3, x2, x1]*













[x3, x1]




[x4, x3, x2, x1]*




x3

y3




[x4, x2, x1]

….










[x4, x1]









x4

y4

























Звездочки означают числа, использующиеся как опорные значения при вычислении значений многочлена.

Из этой таблицы можно написать
y(x)=y1+( x – x1 )([x2, x1]+(x - x2){[x3, x2, x1]+(x – x3)[…]}).
В качестве примера рассмотрим таблицу логарифмов


x

lgx

[,]

[,,]

[,,,]

1

0,0000*

0,30100*







2

0,3010

-0,06245*

+0,01230*

0,23855

3

0,4771

-0,05015

0,20070

4

0,6021








Звездочки обозначают опорные значения. Отсюда получаем
y(x)=0+(x – 1){0,3010+(x - 2)[(-0,06245)+(x - 3)(0,01230)]}.
В частности,



Верное значение

lg2,5 = 0,3979
Таблица разделенных разностей, которая лежит в основе интерполяционной формулы Ньютона, может быть написана в другой, иногда более полезной форме. Заметим, что


В частности,


Таким образом, можно написать


Таблица 2.




[,]

[,,]

[,,,]










x1

y1























[x2, x1]













x2

y2




[x3, x2, x1]
















[x3, x2]




[x4, x3, x2, x1]







x3

y3




[x4, x3, x2]
















[x4, x3]












x4


y4











Несмотря на то, что числа в этой таблице отличаются от чисел в Таблице 1, числа в верхнем ряду в обеих таблицах (которые только и используются в формуле Ньютона), те же самые.


Упражнения


  1. Используя таблицу




x

2

0

3

1

y

8

0

27

1


вычислить интерполяционный многочлен Ньютона. Показать, что y = x3.

2. Вычислить таблицу разностей обоими способами для y=sinx с шагом в 30о:


x

0o

30o

60o

90o

sin x

0,0000

0,5000

0,8660

1,0000


Найти аппроксимационный многочлен.

3. С помощью интерполяционной формулы Ньютона построить многочлен 2-й степени, проходящий через узловые точки:

а)(0,1), (1,4), (2,3); б)(0,-2), (1,3), (2,1); в)(0,2), (1,-1), (2,3).

4. С помощью интерполяционной формулы Ньютона построить многочлен 3-й степени проходящий через узловые точки:

а)(0,2), (1,4), (2,4), (3,3); б)(0,0), (1,2), (2,0) (3,4); в)(0,1), (1,-2), (2,4), (3,5).