main (1160446), страница 9
Текст из файла (страница 9)
. . < 6 , ( ) = , = 0, .§15. Разделенные разности51Разделенной разностью первого порядка, построенной по несовпадающимузлам и , называется отношениеОпределение. ( , ) = ( ) − ( ), − 0 6 , 6 .(1)Обычно мы будем рассматривать разделенные разности, составленные по соседним узлам.Например, (1 ) − (0 ) (2 ) − (1 ) (0 , 1 ) =, (1 , 2 ) =.1 − 02 − 1Замечание.Отношение (1) является дискретным аналогом первой производной.Разделенной разностью второго порядка, построенной по несовпадающимузлам 0 , 1 , 2 , называется отношениеОпределение. (0 , 1 , 2 ) = (1 , 2 ) − (0 , 1 ).2 − 0(2)Пусть даны ( , . .
. , + ) и (+1 , . . . , ++1 ) — разделенные разности -го порядка по соответствующим узлам, где 0 6 , 6 . Тогда разделенной разностью(+1)-го порядка, построенной по несовпадающим узлам , +1 , . . . , ++1 , называетсяотношениеОпределение. ( , +1 , . . . , ++1 ) = (+1 , +2 , . . . , ++1 ) − ( , +1 , . . . , + ).++1 − (3)Введем следующие обозначения:() =∏︁( − ) = 0, (),=0, () =∏︁( − ), = 0, 1, .
. . , , = 0, .=Очевидно, что′0,( )=∏︁( − ),′( ),=∏︁( − ), = , + 1, . . . , .≠==0̸=Покажем, как разделенная разность произвольного порядка выражается через значенияфункции () в узлах { }=0 .Утверждение.Разделенная разность -го порядка представима в виде (0 , 1 , . . . , ) =∑︁ ( )′ ( ) .0,=0Доказательство.Пусть = 1. ТогдаВоспользуемся методом математической индукции. (0 , 1 ) = (1 ) − (0 ) (1 ) (0 )=+.1 − 01 − 0 0 − 1(4)52Таким образом утверждение выполнено при k=1.
Пусть теперь утверждение верно длянекоторого = . Докажем его для = + 1.Следующие соотношения вытекают из предположения индукции: (0 , 1 , . . . , ) =∑︁ ( )′ ( ) ,0,(5)=0 (1 , 2 , . . . , +1 ) =+1∑︁ ( ).′( )=1 1,+1 (6)Запишем разделенную разность ( + 1)-го порядка: (0 , 1 , . . . , +1 ) = (1 , 2 , . . . , +1 ) − (0 , 1 , .
. . , ).+1 − 0(7)Подставим выражения (5) и (6) в равенство (7) и вынесем общий множитель за скобку:1 (0 , 1 , . . . , +1 ) =+1 − 0(︃ +1∑︁=1∑︁ ( ) ( )−′′ ( )( )1,+10,)︃.=0Вынесем за скобку ( + 1)-е слагаемое первой суммы и нулевое слагаемое второй: (+1 ) (0 )++′′(+1 )(0 − +1 )0, (0 ) (+1 − 0 )1,+1(︃ (︃)︃)︃∑︁111+− ′. ( )′+1 − 01,+1( ) 0,( ) (0 , 1 , . . . , +1 ) =(8)=1Рассмотрим отдельно некоторые элементы этого равенства. Заметим, что:′′(0 ),(0 ) = 0,+1(0 − +1 )0,′′(+1 ) = 0,+1(+1 ),(+1 − 0 )1,+1(︃)︃111− ′=′+1 − 0 1,+1( ) 0,( )1=+1 − 0(︃ − 0 − +1− ′′1,+1( )( − 0 ) 0,( )( − +1 ))︃=1.′0,+1( )Подставив найденные выражения в равенство (8), получим: (0 , 1 , .
. . , +1 ) =+1=1=0∑︁ ( )∑︁ ( ) (0 ) (+1 )+ ′+=.′′′0,+1 (0 ) 0,+1 (+1 )0,+1 ( )0,+1( )Утверждение для = + 1 доказано, и в силу индукции справедлива формула (4).Значение функции () в произвольном узле , = 0, можно выразитьчерез значение функции в узле 0 и разделенные разности до порядка включительно.Утверждение.§15. Разделенные разностиДоказательство.53Пусть = 1. Запишем разделенную разность первого порядка: (0 , 1 ) = (1 ) (0 )+.0 − 1 1 − 0Домножим обе части равенства на (1 − 0 ) ̸= 0:(1 − 0 ) (0 , 1 ) = (1 ) − (0 ).Следовательно, (1 ) = (0 ) + (1 − 0 ) (0 , 1 ).Докажем утверждение для = 2. Аналогично предыдущему случаю запишем разделеннуюразность 2-ого порядка и домножим обе части равенства на (2 − 0 )(2 − 1 ) ̸= 0:(2 − 0 )(2 − 1 ) (0 , 1 , 2 ) = −2 − 02 − 1 (0 ) + (1 ) + (2 ).0 − 10 − 1Введем обозначения:=2 − 02 − 0 (1 ) =( (0 ) + (1 − 0 ) (0 , 1 )) =0 − 10 − 1=2 − 0 (0 ) − (2 − 0 ) (0 , 1 ),0 − 1=− (0 )(2 − 1 ).0 − 1Следовательно,(2 − 0 )(2 − 1 ) (0 , 1 , 2 ) = + + (2 ) ==2 − 0(2 − 1 ) (0 ) − (2 − 0 ) (0 , 1 ) − (0 ) + (2 ) =0 − 10 − 1= (2 ) − (0 ) − (2 − 0 ) (0 , 1 ).Выразив из последнего выражения (2 ), получим: (2 ) = (0 ) + (2 − 0 ) (0 , 1 ) + (2 − 0 )(2 − 1 ) (0 , 1 , 2 ).Переход от = к = + 1 для произвольного ∈ производится по аналогиис рассмотренным переходом от = 1 к = 2, но здесь не приводится, так как сопровождается более громоздкими выкладками.
Далее мы иногда будем пользоваться таким приемом,чтобы избегать громоздкости выкладок.Обобщив полученные результаты, запишем формулу для ( ): ( ) = (0 ) + ( − 0 ) (0 , 1 ) + ( − 0 )( − 1 ) (0 , 1 , 2 )++ . . . + ( − 0 )( − 1 ) . . . ( − −1 ) (0 , 1 , . . . , ).Замечание.Формула (9) является дискретным аналогом формулы Тейлора ( ) = (0 ) + ( − 0 ) ′ (0 ) +( − 0 )2 ′′( − 0 ) () (0 ) + . . . + (0 ) + . . . .2!(9)54§16Интерполяционная формула НьютонаРассмотрим вещественную функцию (), ∈ [, ] ⊂ R,заданную в узловых точках произвольного разбиения отрезка [, ]: 6 0 < 1 < 2 < . . . < 6 , ( ) = , = 0, .Воспользуемся результатами предыдущего параграфа и запишем формулу для ( ): ( ) = (0 ) + ( − 0 ) (0 , 1 ) + ( − 0 )( − 1 ) (0 , 1 , 2 )++ .
. . + ( − 0 )( − 1 ) . . . ( − −1 ) (0 , 1 , . . . , ).(1)Подставив в эту формулу вместо , получим полином степени от : () = (0 ) + ( − 0 ) (0 , 1 ) + ( − 0 )( − 1 ) (0 , 1 , 2 )++ . . . + ( − 0 )( − 1 ) . . . ( − −1 ) (0 , 1 , . . . , ).Обозначим полученный полином как (): () = (0 ) + ( − 0 ) (0 , 1 ) + ( − 0 )( − 1 ) (0 , 1 , 2 )++ . . . + ( − 0 )( − 1 ) . .
. ( − −1 ) (0 , 1 , . . . , ).Утверждение.(2)Полином (2) интерполирует функцию ().Доказательство.Для доказательства утверждения достаточно показать, что ( ) = ( ), = 0, .Подставив в формулу (2) вместо , получим: ( ) = (0 ) + ( − 0 ) (0 , 1 ) + ( − 0 )( − 1 ) (0 , 1 , 2 )++ . . . + ( − 0 )( − 1 ) . . . ( − −1 ) (0 , 1 , . . . , ).(3)В равенстве (3) все слагаемые, начиная с -ого, содержат множитель ( − ), равный нулю.Тогда получим ( ) = (0 ) + ( − 0 ) (0 , 1 ) + ( − 0 )( − 1 ) (0 , 1 , 2 )++ . . .
+ ( − 0 )( − 1 ) . . . ( − −1 ) (0 , 1 , . . . , ) = ( ), = 0, ,что и требовалось доказать.Интерполяционный полином, задаваемый формулой (2), называется интерполяционным полиномом Ньютона.Определение.Интерполяционный полином Ньютона тождественно совпадает с интерполяционным полиномом в форме Лагранжа.Замечание 1.Этот факт следует из доказанного в первом параграфе утверждения,что для любой функции () существует единственный интерполяционный полином,построенный по ( + 1) узлу. То есть интерполяцоннный полином Ньютона и интерполяцоннный полином в форме Лагранжа являются различными вариантами записи одного итого же полинома.Доказательство.§17.
Интерполирование с кратными узлами. Полином Эрмита55Так как интерполяционный полином Ньютона тождественно совпадаетс интерполяционным полиномом в форме Лагранжа, он имеет такую же погрешность:⃒⃒+1⃒⃒|()| , где +1 = sup ⃒ (+1) ()⃒.| ()| 6( + 1)!∈[,]Замечание 2.Аналогично случаю с интерполяционным полиномом Лагранжа, если исходная функция является полиномом степени, не превышающей , то интерполяционныйполином Ньютона приближает ее точно.Замечание 3.Выбор формы записи интерполяционного полинома функции () зависитот особенностей каждой конкретной задачи.
Например, если узлы зафиксированы и ихчисло постоянно, а искомая функция меняется, то удобно использовать интерполяционный полином в форме Лагранжа. Если же появляется необходимость в добавлении илиудалении узлов при условии сохранения функции, то удобно использовать интерполяционный полином в форме Ньютона.Замечание 4.§17Интерполирование с кратными узлами. Полином ЭрмитаРассмотрим вещественную функцию (), ∈ [, ] ⊂ R,заданную в узловых точках произвольного разбиения отрезка [, ]: 6 0 < 1 < 2 < . . .
< 6 , ( ) = , = 0, .Пусть, кроме того, в узле заданы значения всех производных функции () до порядка( − 1), = 0, . Натуральное число называется кратностью соответствующего узла .Постановка задачи.щий условию:Определение.Требуется построить полином () степени , удовлетворяю() ( ) = () ( ), = 0, ( − 1), = 0, .Полином () называется интерполяционным полиномом Эрмита.Будем искать интерполяционный полином () в виде () = ∑︁ −1∑︁, () () ( ),=0 =0где , () - полиномы степени .Сформулируем условие, при котором можно найти интерполяционный полиномом Эрмита.Утверждение.Если сумма кратностей узлов функции f(x) равна ( + 1):∑︁ = + 1,=0то существует, причем единственный, интерполяционный полином Эрмита степени для функции ().56Рассмотрение задачи построения интерполяционного полинома Эрмита в общей постановке, которую мы привели выше, выходит за рамки нашего курса.
Интересующиеся могутобратиться к [1], мы же далее будем рассматривать частный случай: построение интерполяционного полинома Эрмита для функции () по трем узлам, один из которых имееткратность 2.Построение полинома Эрмита по трем узламРассмотрим функцию (), определенную вместе со своей первой производной на отрезке[, ]. Построим для функции () интерполяционный полином Эрмита 3 () по трем узлам0 , 1 и 2 : 6 0 < 1 < 2 6 , где узел 1 — кратный.По определению интерполяционного полинома Эрмита для 3 () должны выполнятьсяследующие равенства:3 (0 ) = (0 ), 3 (1 ) = (1 ), 3 (2 ) = (2 ), 3′ (1 ) = ′ (1 ).(1)Будем искать полином Эрмита 3 () в виде3 () = 0 () (0 ) + 1 () (1 ) + 2 () (2 ) + 1 () ′ (1 ),(2)где 1 () и (), = 0, 2 — полиномы третьей степени.Равенства (1) и (2) позволяют сформулировать условия нахождения коэффициентов 1 ()и (), = 0, 2:0 (0 ) = 1, 1 (0 ) = 0, 2 (0 ) = 0, 1 (0 ) = 0,0 (1 ) = 0,1 (1 ) = 1,2 (1 ) = 0,1 (1 ) = 0,0 (2 ) = 0,1 (2 ) = 0,2 (2 ) = 1,1 (2 ) = 0,′0 (1 ) = 0,′1 (1 ) = 0,′2 (1 ) = 0,′1 (1 ) = 1.Воспользуемся этими условиями и получим коэффициенты интерполяционного полинома (2) в явном виде.Из условий 0 (1 ) = 0, 0 (2 ) = 0 и ′0 (1 ) = 0 следует, что узлы 1 и 2 являются корнями полинома 0 () двойной и единичной кратности соответственно.
Поэтому коэффициент0 () будем искать в виде0 () = ( − 1 )2 ( − 2 ),где ∈ R.Для нахождения воспользуемся условием 0 (0 ) = 1:0 (0 ) = (0 − 1 )2 (0 − 2 ) = 1.Поделим это равенство на (0 −1 )2 (0 −2 ) (мы можем это сделать, так как узлы 0 , 1 , 2различны):1=.(0 − 1 )2 (0 − 2 )В дальнейшем при делении на множители, содержащие разности узлов, мыне будем оговаривать неравенство нулю этих множителей, считая это очевидным.Замечание.Запишем представление для 0 () с учетом выражения для коэффициента :0 () =( − 1 )2 ( − 2 ).(0 − 1 )2 (0 − 2 )§17. Интерполирование с кратными узлами. Полином Эрмита57Очевидно, что коэффициент 2 () имеет аналогичную структуру с двукратным корнем 1и однократным корнем 0 :( − 1 )2 ( − 0 )2 () =.(2 − 1 )2 (2 − 0 )Рассмотрим коэффициент 1 (), для которого точки 0 , 1 , 2 являются однократнымикорнями.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
