Простое решение интерполяционной задачи Эрмита Текст научной статьи по специальности «Математика»

DOI: https://doi.org/10.23670/IRJ.2018.70.028

Шустов В.В.

ORCID: 0000-0002-2465-7475, Кандидат технических наук, ФГУП Государственный научно-исследовательский институт авиационных систем, Москва ПРОСТОЕ РЕШЕНИЕ ИНТЕРПОЛЯЦИОННОЙ ЗАДАЧИ ЭРМИТА

Аннотация

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

Ключевые слова: задача Эрмита, многочлен Тейлора, формула для интерполяционного многочлена Эрмита, интерпретация многочлена Эрмита.

Shustov V.V.

ORCID: 0000-0002-2465-7475, PhD in Engineering, FSUE State Research Institute of Aviation Systems, Moscow SIMPLE SOLUTION OF HERMIT INTERPOLATION PROBLEM

Abstract

The Hermite interpolation problem on the construction of a polynomial taking given values of a function and its derivatives at the nodal points is considered. A new way of solving this problem is proposed. Using the Taylor formula, we simplified the derivation of the Hermite interpolation polynomial formula. An interpretation of this polynomial in the representations of the Taylor polynomial is proposed. An algorithm is given for obtaining finite formulas for the Hermite polynomial in cases when the maximal orders of the derivatives given at the node points are the same.

Keywords: Hermite problem, Taylor polynomial, the formula for the Hermite interpolation polynomial, interpretation of the Hermite polynomial.

Введение

Интерполяционный многочлен Эрмита известен давно [1], [2, С. 163], [3, С. 71], [4, С. 154], [5, С. 143], [6] и используется в теории и практике для различных целей. В частности, он применяется для интерполяции функций, заданных в виде таблиц, в которых заданы значения функции и значения ее производных, и может использоваться для аппроксимации функций, принадлежащих определенному классу гладкости. В этом случае для повышения точности приближения используются данные о значениях функции и ее производных в нескольких точках отрезка.

Данный подход в случае выбора аппроксимирующей функции в классе многочленов приводит к задаче интерполирования Эрмита [2, С. 163-165] и рассмотренной в [1, С. 70-79]. Для решения задачи интерполяции строится многочлен, называемый интерполяционным многочленом Эрмита.

Вывод представления для интерполяционного многочлена Эрмита в общем случае для произвольного числа узловых точек приводится, например, в [2, С. 169-172], [3, С. 71-73], [4, С. 154-155] и является сложным и достаточно объемным. Вывод представления связан с предельными переходами и при некоторых подходах требует использования методов теории функций комплексного переменного.

Сложность вывода общей формулы для многочлена Эрмита побуждает ряд авторов тем или иным путем упростить его. Например, в работе [7] упрощение ищется в рамках линейной алгебры путем обращения матрицы.

Целью настоящей работы является с одной стороны дальнейшее упрощение вывода представления для интерполяционного многочлена Эрмита, с другой стороны разработка способов представления этого многочлена в конечном виде для частных случаев.

Постановка и решение задачи

Пусть в n+1 различных точках отрезка [x0,xn], удовлетворяющих условию

Х0 < Xj < ... < xn ,

заданы значения функции fx) и ее производных до порядка а—1 включительно:

fa)(x,) = f,u\ j = 0,1,...,«, -1, i = 0,1,..., n. (1)

Вслед за [3, С. 71] будем считать, что производная нулевого порядка есть сама функция, т.е. f0/)(x)=fx).

Необходимо построить интерполирующую функцию в классе многочленов, которая удовлетворяла бы условиям (1). Для решения задачи введем следующие обозначения.

Поставим в соответствие каждой узловой точке xi функцию ю,(х) определяемую формулой:

n

(x) = (x - Xo)«0...(x - xt-i)ai-1( x - xM)a'+1...(x - Xn Г =П (x - x Г . (2)

I=0,1 *t

Функцию ю,(х), которая также является многочленом, можно записать в виде:

©.( X) = Q(x) , где (3)

'( ) (x - X )"•

n

Q(x) = (x-x0 Г (xx Г...(x-^ Г = П(x - x Г . (4)

I=0

Функции ю,(х), как это следует из формулы (2), в точках xl обращаются в нуль вместе со своими производными до порядка а,—1 включительно и не обращаются в нуль в собственных точках xi, т.е.

($P(Xj)=0 при ijj , /=0,1,.. ,,а—1 и

■|( j) Jx=x,

(7)

. ...... (5)

юДх^фО при /=/. (6)

Через Tjа' I(х)) обозначим многочлен Тейлора степени ai-1 для функции_Дх), построенный для точки х,, по значениям производных, определенных условием (1):

Г-1(/ (X)) = ? [/ (х)]

7=0 7!

В некоторых случаях верхний индекс в обозначении Tiа' 1 (I(х)) будем опускать.

Представим Н(х) в виде суммы многочленов Н,(х), где Н,(х) строится для каждой узловой точке х,, т.е.

п

Н (х) = 3 Я,. (х), (8)

1=0

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

I ( х)

f (х) = Щ ( X)

л(х)

т.е. выделим из функции fx) первый сомножитель ra,(x), который в соответствии с (2) обращается в нуль во всех узловых точках, кроме xi, и второй сомножитель f(x)/rai(x).

Многочлен Hi(x) для каждой узловой точке xi построим в виде произведения двух сомножителей. В качестве первого сомножителя возьмем функцию ra,(x), которая в соответствии с (5) и (6) обладает теми свойствами, что она обращается в нуль сама вместе со своими производными до порядка а,-1 включительно в каждой узловой точке x/ , отличной от точки xi. Второй сомножитель f(x)/ra,(x) , так как rai(xi)^0, представим в виде многочлена Тейлора для точки xi степени а—1 включительно с учетом условия (1), что в точке xi заданы производные до порядка а,—1 включительно. В соответствии с рассуждениями, изложенными выше, можно записать, что

H (x) = л( x)T

f f (x) ^ x)

или, записывая многочлен Тейлора в явном виде согласно (7), как:

Иг (х) = ®,( х)]Г( Х " х)'

j=0

j!

f ( X) х)

(j)

(9)

(10)

Подставляя Н,(х) из (9) в (8) получим, что суммарный интерполяционный многочлен Н(х) представляется в следующем простом виде:

H ( x) = £®,.( x)T

f f (x) ^

(11)

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

И ")(х ) = И °)(х,) .

Это объясняется тем, что в соответствии с (5) и (6) значения Н,(х/) и их производных в точках ху при уф равны нулю.

Многочлен Н(х) с учетом представления многочлена Тейлора формулой (7) может быть записан в виде:

н (x)=I

i=0 j=0 j !

f ( x)

x)

(j)

(12)

Степень интерполяционного многочлена

Степень 5 многочлена Н(х) определяется степенями • многочленов Н,(х), которые в свою очередь определяются в

соответствии с формулой (9) суммой • степеней сомножителей ю,(х) и Tiа' 1 (I(х)), т.е.

= (а0 + а1 +...+ а,-1+ а,+1...+ ап)+(а, - 1), в которой скобками выделены степени сомножителей. Эту формулу для степени • многочлена Н,(х) можно записать в виде:

( п Л ^ =\3аг - 1 ,

V 1=0 )

из которого видно, что степень • многочлена Н,(х) для всех значений ,, определяемых условиями (1), одна и та же.

Из этой формулы и из того, что степень суммы многочленов не превосходит максимума из степеней многочленов суммы, следует, что максимальная степень • многочлена Н(х), определяемого (8), выражается формулой:

5 =

I

Л

а.

ч i = 0 /

-1.

Проверка заданных условий для многочлена

Покажем, что построенный многочлен Н(х), представленный формулой (12), удовлетворяет и условиям (1).

х=х

i=0

Используя формулу для производной /-го порядка произведения двух функций (см. например, [8, С. 149])

(/(х) я (х))( Л =£ с]/(^)(х) я(к)(х)

(14)

к=0

и свойство многочлена Тейлора, выраженное формулой

Т(-)(/(х))| = /(-)(х,),

1х=х,

найдем производную/-го порядка многочлена Н(х), определенного формулой (11):

н(ЛЧх) = ££с>г (Л-к Т

(к)

г=0 к=0

Л

/ ( х)

чк(х) У

Далее, вычислим значение этой производной в точке х=х,- , заменив обозначение переменной внешнего суммирования / на I:

п - г

н(-)(х,)=ц < ®(^)(х)

Т

(к)

г

1=0 к=0

/ ( х) к( х)

к ]-к)(х)

=1с

к=0

/ ( х) к(х)

(к)

(15)

Заметим, что при суммировании по переменной I после подстановки х=хг, в силу свойства (5) функции ю,(х) остаются только слагаемые, отвечающие условию 1=1, т.е. знак суммирования по I исчезает. Далее, так как ю,(х,)^0, то при дифференцировании тождества

/ (х) = к (х)-/ (х)

к(х)

/ раз, использовании формулы (14) и подстановке х=х„ получим что

/(-)(х1) = £ек к

к=0

к ^ (1-к)

(х)

/ (х) к (х)

(к)

(16)

Так как правые части формул (15) и (16) равны между собой, то по свойству числовых равенств и левые части равны друг другу, т.е. для всех/ = 0,1,.. .а-1 и для всех / = 0,1,.. .и имеет место:

Н\*)(х.) = /(Л)(х.),

а это и означает, что построенный многочлен Н(х) удовлетворяет условиям (1).

Выделение производных в явном виде

Для того, чтобы выделить значения производных в узловых точках в явном виде, выполним некоторые преобразования. Используя формулу (14) для определения производной /-го порядка от произведения двух функций, формулу (12) для Н(х) можно записать в виде:

Н (х) = ( х)£

1=0 -=0

(х - х,) - >

- к=0

I ск / (х)1

( к)

1

или с учетом обозначений (1) как

п т, ( у -

Н (х) = 1 к (х)1 (к)

к, (х) 1

к(х)

(--к)

(--к)

_..............(17)

1=0 -=0 — к=0

Здесь и, далее, буквами тг , /=0,1 обозначены порядки наивысших производных, используемых для построения двухточечного многочлена Н(х), т.е.

тг = а-1, откуда аг= тг +1, (18)

и используется формула для биноминальных коэффициентов

т!

- . (19)

к !(т - к)!

Для того, чтобы вынести обозначения производных из-под знака внутреннего суммирования, в формуле (17) для Н(х) поменяем порядок суммирования по переменным/ и к и, соответственно, пределы суммирования:

ск =-

т

щ т

Н (х) = Лк (х)1 I с-/Г

1=0

1

к(х)

(--к)

к=0 -=к —

Перейдем к новой переменной суммирования I, определенной формулой /=/-к, из которой следует, что /=/+к. После перехода к переменной I изменятся пределы внутреннего суммирования: при/=к имеет место 1=0, при /=тг имеет место 1= т—к. Тогда формула для Н(х) примет вид:

п т, т,-- (у. у \к+'

н (х) = 1 к, (I с1,/;к)

1=0 к=0 I=0 (к + 1)!

1

к(х)

(I)

Еще раз поменяем обозначения для внутренних переменных суммирования: к=/ и 1=к. Тогда Н(х) можно записать как

х=х

х=х

х=х

х=х

т

х=х

х=х

п т т-] _у-у+к

Н(х) = 3^(х)! ъ^^хкг)

1=0 ]=0 к=0 (] + к)!

щ( х)

(к)

После подстановки в эту формулу выражения для биноминальных коэффициентов (19) и сокращения числителя и знаменателя дроби на (¡+к)! получим, что

п т, щ- ] / \]+к

И(х)=?щ(х)? '

1=0 1=0 к=0 ГК!

1

Щ( х)

(к)

(20)

Если внести сомножитель ю,(х) под знак суммирования по переменной к, то формулу (20) можно записать как:

п Щ Щ-7 (V _ V \]+к

И(х) = ? ? ?!I(

1=0 7=0 к=0 у к!

1

щ(х)

(к)

(х) .

Эта формула для многочлена Н(х) с точностью до обозначений (см. формулы (3) и (4)) совпадает с формулой для этого многочлена, приведенной, например в [2, С. 172] и полученной другим путем, что дает дополнительное подтверждение верности представления многочлена в форме (11).

Если в формуле (20) вынести члены, не зависящие от индекса суммирования по к в внешний цикл, то формула для многочлена Н(х) примет вид:

п т ^ (з) т-] 1

Н (х) = ? щ (х)? (х - х )3 ? -

1=0 3=0 3 ! к=0 к!

1

Щ( х)

( к )

(х - х,. )к.

(21)

В соответствии с полученными результатами имеет место следующая теорема.

Теорема

Пусть функция /(х) определена и достаточное число раз дифференцируема на отрезке [х$,хп], а также удовлетворяет условиям (1). Тогда существует многочлен степени не выше •, определенной формулой (13), который представляется в любой из форм (11), (12) или (21).

Интерпретация интерполяционного многочлена

Если в формуле (21) внести сомножитель ю,(х) под знак суммирования по переменной у и сгруппировать сомножители, то формула для Н(х) примет следующий вид:

п Щ Л7У'

И (х) = ? З1 (х-х1)

1=0 у=0 у-

(

Щ-1 1

щ( х) ! 1

к=0

к-

1

Щ( х)

(к) Л

(х - х, )к

Обозначим выражение в скобке как ф,у(х), т.е.

т1-11

V ( х) = Щ( х) ? -

г-=п к-

щ( х)

(к)

(х - х, )к

(22)

(23)

(24)

и назовем ф,у(х) поправочной функцией.

С учетом обозначения (23) формулу (22) для Н(х) можно записать в виде:

п т Л3у

Н (х) = ? ? Л (х - х ) Ф, (х)

1=0 3=0 3 !

Из представления (24) видно, что многочлен, построенный по заданным значениям функции и ее производных в точках отрезка [х0,хп], есть сумма модифицированных многочленов Тейлора, построенных для каждой узловой точки х, этого отрезка. Модификация разложения заключается в том, что каждый член многочлена Тейлора в данной точке х, умножается на поправочную функцию, зависящую как от расстояния до других узловых точек, так и от порядков наивысших известных производных в этих точках.

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

Частные случаи многочлена Эрмита. Отметим частные случаи многочлена Эрмита. Предварительно путем последовательного дифференцирования выпишем следующие соотношения для производных функции, выражение для которой является обратным для функции ю(х):

' — (0) ,

(25)

1

а

а

1

а а\ а

а

2а'2-щщ" _ 1 (2Щ2-аа"') а

а

а

(26)

(27)

Случай 1

Пусть во всех узловых точках отрезка заданы функции и значения ее первой производной, что соответствует условию:

х=х

х=х

х= х

1

х=х

а

а

1

ГГ

щ = 1, 1 = 0,1,...,п .

С учетом формул для значений производных (25)-(27) формула (22) для многочлена Н(х) при обозначении его как Н1(х) примет вид:

н (х) = I к \ / (1 - — Ах,.) + /АхЛ , где Ах, = х-х,- и юй(х)= ю,(х,).

.=0 к I к ]

С точностью до обозначений эта формула соответствует формуле, приведенной в [3, С. 74] .

Случай 2

Пусть во всех узловых точках отрезка заданы функции и значения ее первой и второй производной, что соответствует условию

щ = 2, 1 = 0,1,..., п .

В этом случае формулу (22) для многочлена Н(х) при обозначении его как Н2(х) с учетом формул для значений производных (25)-(27) можно записать в виде

Н2(х) = 1 к\/\Ах, +к-кк^) + /'(1 Ж +1 /"ахА.

.=0 к- I к 2!к к 2! ]

Продолжая ряд формул (25)-(27) и используя (22), можно получить аналогичные формулы для многочлена Н(х) для последующих значений т=3,4... - максимального порядка заданных в узлах значений производных, не решая частных задач интерполяции для каждого конкретного значения т, как это делается для т=1 в [2, С. 167].

Случай 3

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

Формула (3) для ю,(х) с учетом (4) записывается в виде

к (х) = (х - х^. Щ +1, I = 0,1. (28)

С использованием (28) легко устанавливается, что производная к-го порядка от обратной величины ю,(х) представляется как

1 1 _(т^ + к)! (-1)к

( x)

( x - x^l)m-+l¥]

и, соответственно, значение этой производной при x=xi выражается соотношением

-|( k )

1 (m_t + k)! i i

— ( x)

mi-, ! (x,. - xi-, )mi-i+1 (xi-, - x,. )k

(29)

С учетом (28) и (29) формулу (21) для многочлена Н(х) можно записать в виде

^ (х-х1-1 Г-+1 Щ ^ .-Щ--К- + к)! (х-х;)к

H(x) = ZZ^(x - x,)j Z(

tG(x, - xi-,y^1 £ - ti ]!ml- ! (xi-, - x,)]

Введя обозначения для относительных переменных £

Jx-x)!, , = Gl * ( xi-,- x, )' ,

k

и используя выражение для биномиального коэффициента С , определенного (19), получим конечную формулу для

двухточечного интерполяционного многочлена Эрмита

1 m f( j ) m - j

H (x) = ZC-Г+1 ZfT- (x - x, )J Z<, , (30)

i=G j=G j! ] =G

которая соответствует представлению этого многочлена, полученному в [9, С. 1080].

Подробное исследование этого случая и применение двухточечных интерполяционных многочленов Эрмита для аппроксимации функций при симметричном и несимметричном распределении производных на концах отрезка и его оптимизация проведено в [9] и [10], соответственно.

Заключение

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

Список литературы / References

1. Hermite Sh. Sur la formule d'interpolation de Lagrange I Sh. Hermite II Journal fur die reine und angewandte Mathematik. - 1S7S. - Vol. S4. - P. 70-79.

2. Березин И. С. Методы вычислений. Т. 1 I И. С. Березин, Н. П. Жидков - М.: Физматлит, 1962. - 464 с.

3. Гончаров В.И. Теория интерполирования и приближения функций. I В. И. Гончаров. - М.: Гостехтеориздат, 1934. - 316 с.

ml-,!

4. Прасолов В. В. Многочлены. / В. В. Прасолов - М.: МЦНМО. - 1999. - 336с.

5. Микеладзе Ш. Е. Численные методы математического анализа. / Ш. Е. Микеладзе - М.: Гостехтеориздат, 1953. -528 с.

6. Spitzbart A. A generalization of Hermite's interpolation formula / A. Spitzbart // Amer. Math. Monthly, 1960. - Vol. 67. - P. 42-46.

7. Stanica D. A short proof of the Hermite's formula for polynomial interpolation using operator methods / D. Stanica // Annals of the Alexandra loan Cuza University - Mathematics. - 2015. - Vol 61. - P. 209-215. doi: 10.2478/aicu-2013-0040

8. Кудрявцев Л. Д. Математический анализ. Т. 1. / Л. Д. Кудрявцев - М.: Высшая школа, 1970. - 590 с.

9. Шустов В. В. О приближении функций двухточечными интерполяционными многочленами Эрмита / В. В. Шустов // Журнал вычислительной математики и математической физики. - 2015. - Т. 55. - № 7. - С. 1091-1108.

10. Шустов В. В. Аппроксимация функций несимметричными двухточечными многочленами Эрмита и ее оптимизация / В. В. Шустов // Журнал вычислительной математики и математической физики. - 2015. - Т. 55. - № 12.

- С. 1999-2014.

Список литературы на английском языке / References in English

1. Hermite Sh. Sur la formule d'interpolation de Lagrange / Sh. Hermite // Journal fur die reine und angewandte Mathematik. - 1878. - Vol. 84. - P. 70-79.

2. Berezin I. S. Computing Methods. / I. S. Berezin, N. P. Zhidkov. Vol. 1 - Pergamon: Oxford, 1965. - 464 P.

3. Goncharov V. I. Teoriya interpolirovaniya i priblizheniya funktsij [The theory of interpolation and approximation of functions] / V. I. Goncharov. - M.: Gostekhteorizdat, 1934. - 316 P. [in Russian]

4. Prasolov V. V. Mnogochleny [Polynomials] / V. V. Prasolov. - М.: М^МШ, 1999. - 336 P. [in Russian]

5. Mikeladze Sh. E. Chislennye metody matematicheskogo analiza [Numerical methods of mathematical analysis] / Sh. E. Mikeladze. - M.: Gostekhteorizdat, 1953. - 528 P. [in Russian]

6. Spitzbart A. A generalization of Hermite's interpolation formula / A. Spitzbart // Amer. Math. Monthly, 1960. - Vol. 67.

- P. 42-46.

7. Stanica D. A short proof of the Hermite's formula for polynomial interpolation using operator methods / D. Stanica // Annals of the Alexandra Ioan Cuza University - Mathematics. - 2015. - Vol 61. - P. 209-215. doi: 10.2478/aicu-2013-0040

8. Kudryavtsev L. D. Matematicheskij analiz [Mathematical analysis] T. 1 / L. D. Kudryavtsev - M: Vysshaya Shkola, 1970. - 590 P. [in Russian]

9. Shustov V. V. Approximation of functions by two-point Hermite interpolating polynomials / V. V. Shustov // Computational Mathematics and Mathematical Physics. - 2015. - Vol. 55. - No 7. - P. 1077-1093. doi: 10/1134/S0965542515040156

10. Shustov V. V. Approximation of functions by asymmetric two-point Hermite polynomials and its optimization / V. V. Shustov // Computational Mathematics and Mathematical Physics. - 2015. - Vol. 55. - No 12. - P. 1960-1974. doi: 10/1134/S0965542515120155