Главная » Просмотр файлов » Самарский А.А. Гулин А.В. - Численные методы

Самарский А.А. Гулин А.В. - Численные методы (1078412), страница 25

Файл №1078412 Самарский А.А. Гулин А.В. - Численные методы (Самарский А.А. Гулин А.В. - Численные методы) 25 страницаСамарский А.А. Гулин А.В. - Численные методы (1078412) страница 252018-01-10СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 25)

Интерполирование алгебраическими многочленами 1. Интерполяционная формула Лангранжа. Пусть на отрезке ач-'х(Ь заданы точки хь (с=О, 1,..., и (узлы интерполировании), в которых известны значения функции ((х). Задача интерполирования алгеораичеекими многочленами состоит в том, чтобы построить многочлен 1.„(х) = а, + а,х+... + а„х" (1) степени и, значения которого в заданных точках х„, й= О, 1, ..., л, совпадают со значениями функции ((х) в этих точках. Для любой непрерывной функции ("(х) сформулированная задача имеет едннственное решение.

Действительно, для отыскания коэффициентов а„а„..., а„получаем систему линейных уравнений а,+а,х;+а,х,'.+ ... +а„х',!=((х,), с=О, 1, ..., и, (2) определитель которой (определитель Вандермонда 112, с. 33]) отличен от нуля, если среди точек х„с'=О, 1,..., и, нет совпадаюшнх. Многочлен (.„(х), удовлетворяющий условиям (.„(х,) =((хс), (=О, 1, .... и, (3) называется интерлоляционным многочленом для функции ((х), построенным по узлам (хс)0.

127 значений функции 1(х) в узлах интерполирования. Найдем явное выражение для коэффициентов с„(х). Из условий интерполирования (3) получаем а 'Я сь(хД)(хь) =~(х;), 1=0, 1,, и. ь=о Зги соотношения будут выполнены, если на функции с,(х) нало- жить условия О, (~й, 1, 1 = к, 1 = О, 1, ..., а, которые означают, что каждая из функций с„(х), й=О, 1, ..., и, имеет не менее и нулей на 1а, 5). Поскольку 1.„(х) — многочлен степени и, коэффициенты с,(х) естественно искать также в виде многочленов степени и, а именно в виде с~(х) =1~(х — хр) (х — х~) ° .

° (х хд-~) (х х~~~)..(х х ) ° Из условия с,(х,) =1 находим Х," = (хь — х,)(хь — х,) ... (хь — хь,)(хь — хь„) ... (хз — х„). Таким образом, коэффициенты с,(х) ннтерполяционного многочлепа (4) находятся по формулам Д (х — х;1 сь(х) = П "-;) 1~э (5) Часто коэффициенты с„(х) записывают в другом виде. Введем м огочлен ы(х) степсни и+1: ы (х) = (х-- х,) (х — х,)...(х — х,,) (х — х,) (х — х„+,)...(х — х.') (6) и вычислим его производную в точке х,: м'(х,) = (х,— х,)... (х,— х,,) (х,— х„+,)... (х,— х„), Тогда получим, что сь(х) =— и (х) (х — х ) и' (х„) 1га Решение системы (2) можно записать различным образом.

Наиболее употребительна запись интерполяционного многочлена в форме Лагранжа н в форме Ньютона. Интерполяционная формула Лагранжа позволяет представить многочлен (.„(х) в виде линейной комбинации н 1.„(х) =- 'Я сь (х) ~ (хх) э=а Итак, интерполяу(ионный многочлен Лагранжа имеет вид л 2=2 или, более подробно, Т22 (х — х;) Т., (х) = 'Я и 1(хй). (8) 2=.2 ~ 1 (х — х,) (Фх 2.

Интерполяционная формула Ньютона. Эта формула позволя. ет выразить интерполяционный многочлен („(х) через значеннс 1(х) в одном из узлов и через разделенные разности функции 1(х), построенные по узлам х„х„..., х„, Она является разностным аналогом формулы Тейлора У(х) =Р(хй)+(х — хр)Г(х.)+ '" "' Г(хй) + . Сначала приведем необходимые сведения о разделенных разностях. Пусть в узлах х,~(а, Ь), й=О, 1, ..., п, известны значения функции 1(х). Предположим, что среди точек х„, (г=О, 1, ..., и, нет совпадающих. Разделенными разностями первого порядка называются отношения 1(.,) — 1(;) 1(хих()= ', (,1'=0,1, ..., и, (Ф1.

х( — х, Будем рассматривать разделенные разности, составленные по соседним узлам, т. е. выражения )(хй, х,), 1(х„хй), ..., 1(х„„х„). По этим разделенным разностям первого порядка можно построить разделенные разности второго порядка: 1(х„х,) — 1(хй, хд ('(х„х„х,) = Х2 — Хо 1(22 12) /(х1. 22) ) (х1, х„х,) = Хй — Х, /(хх 1 хх) ~(х 2 х 1) ~(хл-й, хл-д, хл)— ХП Хй 2 Аналогично определяются разделенные разности более высокого порядка. Например, если известны разности я-го порядка 1(хн Ху+1, ..., Х122), 1(х12.1, Хуег, ° ° ., Х11221), то разделенная разность (я+!)-го порядка определяется как 1(хн х;„, ..., хмй, х„й.„) = l (х' х 12 х' 2 1) 1(х> х' 1 хг22) х;,А „— х( 5 А. А.

Сйчйрсхйй, А. В. Гулни При вычислении разделенных разностей принято записывать их в виде таблицы Г (хо) хо Г(хо, хй / (хо, х„хо) ! (х,! х, Г" (х„х,! ~ (х„хо... х„) ! (х,! хо 1(х„„х„, х„) ! (х„„х„) )(х„) хо Разделенная разность л-го порядка следуюпгим образом выражается через значения функции ! (х) в узлах: тоо о( !'(х;, х;„„..., хыо) = '~" '=' П(,—,) ойдо (9) Эту формулу можио доказать методом индукции. Нам потребуется частный случай формулы (9): ! (х,! )(х„хо ..., хо) =",„' ( ! (х; — х~! С =-а.

оФо +~ (х,. — хо) (х; — х,) (х; — х; ) (х; — х; ) (х; — х„! о=а Покажем, что многочлен Р„(х) совпадает с многочленом Лагранжа (8), Рассмотрим наряду с Е„(х) миогочлены Е„(х) =~(х,), !30 Интер)голяционным многочленом Ньютона называется много- член Р„(х) = = ) (хо) + (х — х,) ) (х„ х,) + (х — х,) (х — х,) ( (х„ х„ х,) + ... ... + (х — х,) (х — хо) ... (х — х„,) ) (х„х„..., х„). ( ! ! ) Е,(х), ..., 1.„,(х) и представим 1.„(х) в виде Е„(х) = 1., (х) + '~~~ (Е! (х) — 1.!, (х)).

(! 2) Из условий интерполяции (3) получаем, что Е,,(х,) =Е.,(х„) =~(х,) при а=О, 1, ..., !' — 1 и !'=1, 2, ..., п. Следовательно, разность Ь,(х) — А,,(х) представляет собой алгебраический многочлен степени 1, который обращается в нуль в точках х„х„..., х, „т. е. Е,;(х) — Х.;,(х) =А;(х — х,)(х — х,)...(х — х,,), (13) где А, — числовой коэффициент. Этот коэффициент нзходится из условия У,(х,) †,,(х,) =А;(х, — х,)(х; — х,)...(х, — х~,), откуда, учитывая условие Ь;(х;) =~(х,), получаем А ! -» ! !( ~) — ь!»(»!) (14) (» — х»)...

!»! — х ) Подставим в (14) вместо слагаемого Ь,,(х,) его значение, вычисленное согласно (8), т. е. Е!, (х!) = (х~ — »0) (х! — »>) ... (»~ — х» ,) (х! — хы,) ... (»т — к! ,) = ~ч ', г (х») (»ь — »0) (»» — »,) ... !»» — х» )(»» — х» ,! ... (»» — »~,) х. ) ! !»»! Сопоставляя это выражение с (!0), видим, что А, совпадает с разделенной разностью 1-го порядка: А,=1(х„х„..., х;).

Отсюда и из (12), (13) приходим к интерполяционной формуле Ньютона У.,(х) = = ! (х,) + (х — х,) ! (х„х,) + (х — х„) (х — х,) ! (х„х„х,) +... ... + (х — х,)(х — х,)... (х — х,,))(х„х,, ..., х„), (15) 5* пы Тогда получим ! !»г! (х- — »0) ...(».— I с-1 р(» ! Х Х(,,„) ! =Х ь !», — ха! (хь — »~) ...(»» — »» ,! (»» — »ь ,)...(»» — х- ) (»» — »ь )!»,— х» ! ... (»» — ».

)!»,— »! Подчеркнем еще раз, что формулы (а) и (15) представляют собой различную запись одного и того же многочлена а„+ а,л + а,хз+... + а„х", удовлетворяющего условиям интерполяции (2). Интерпаляционную формулу Ньютона удобнее применять в том случае, когда интерполпруется одна и та же функция )(х), но число узлов интерполяции постепенно увеличивается. Если узлы интерполяции фиксированы и нптерполируется не одна, а несколько функций, то удобнее пользоваться формулой Лагранжа, 3 а м е ч а н и е. При выводе формулы (16) не предполагалось, что узлы хь хь ..., х„расположены и какоч-то определенном порядке. Поэтому роль точки х, а формуле (16) может играть любая из точек хм х„..., х„.

Соотаетстиуюшее множество интерполяпионных формул можно получить из (161 гере- нумераций уэлоа. Например, тот же самый многочлен Е„(х) можно предстаиить и инде 1.„(х) = 1(хп) и- (г — хп) 1(ххи х„) + — (х — хи) (х — х„,) 1(х„, хп а хп,) -1- ... ... +(х — хи)(х — х„,)... (х — хйг(х„, хп „..., л;,). (16) Если хя(х,я.хз<... <х„, то (16) называется формулой интерполирования вперед, а (!61 — формулой интерполирования назад $ 2. Погрешность интерполирования 1.

Остаточный член интерполяцианной формулы. Заменяя функцию!1(х) интерполяционным мпогочленам 1.„(х), мы допускаем погрешность г„(х) =1(х) -1.„(х), которая называется погрешностью интерполирования илн, чта та же самое, остаточным членом интерполяуионнай формулы. Ясно, что в узлах интерполирования эта погрешность равна нулю. Оценим погрешность в любой точке хе=[а, Ь).

Для этого рассмотрим вспамогательнуча функцию ь (5) =1(5) — 1.„(5) — Кь> (5), (1) где зен(а, Ь), К вЂ” постоянная и ю(5) = (5 — х,) (5 — х,)... (5 — л„). (2) Пусть требуется оценить г„(х) в заданной точке х~(а, Ь], не являющейся узлом интерполирования. Выберем постоянную К из условия д(х) =О. Для этого достаточно положить К= 1(х) — 1.и (х) ы (х) Предположим, что 1(5) имеет и+1 непрерывную производную на отрезке а (5 .Ь. Функция у(5) имеет не менее и+2 нулей на этом отрезке, а именно в точках х, х„, й=б, 1, ..., л.

Поэтому производная у'(5) имеет не менее чем и+! нулей на (а, Ь), ди(5)— 132 не менее и нулей н т. д., функция д'"оо(в) по крайней мере один раз обращается в нуль на [а, Ь). Тем самым существует точка 5~[а, Ез!, в которой й'"+о(й) =О. Поскольку ошоп (з) — [(оо-и (з) (п 1 1) 1Е( получаем ы (к) Таким образом доказано, что погрешность интерполирования можно представить в виде Еш'и ($) Е (х) — Е.„(х) = в (х), (3) (и+ 1)! еде фа=[а, Ь[ и в(х) — многочлен, определенный согласно (2), Отсюда следует оценка [Е (х) — Е.„(х) [~~ "" [в (х) [, (4) где М„., = зпр [Е"'+"(х) [. В частности, если Е(х) — алгебрапче.=Го.ь1 ский многочлен степени п, то интерполирование, проведенное по любым точкам х„х„..., х„, осуществляется точно, т.

с. Е.„(х) =— = — Е(х). 3 а меч ание. Наряду с иитерполировзязем применяют и зкстрплолиро. ванпс, т. е. вычисление значений функции Е(к) в тачках хои[а, Ь[ по приближен- ной формуле Е(к) жь (к), где ь„(к) — иятерполяцяаяяый мзагачлен. Одиака погрешность экстрзполираззяия обычно акззызается сушостяензо большей, чем погрешность иятсрполяроязяяя. К этому выводу можио прийти, рассматривая поведение мнагочлеяа в(к) внутри и зие отрезка [а, Ь), Поскольку мпогочлены Лагранжа и Ньютона отличаются толь- ко формой записи, представление погрешности в виде (3) справед- ливо как для формулы Лагранжа, так и для формулы Ньютона.

Однако погрешность интерполирования можно представить в в дру- гом виде. Для этого рассмотрим разделенную разность Е(х, х„х„..., х„)— Е (к) и + (к — ко) (к — «о) ° ° (к — «„) + Е(ко) [ [ Е( л) («о — к) (ко — кз) ... (ко — ко) («„— К) (к„— «о) ... (к„— клоп) имеющую порядок и+1. Отсюда найдем Е (х) Е (хо) (« — к,) (к — ко) ... (к — «„) " +".

Характеристики

Тип файла
DJVU-файл
Размер
6,53 Mb
Тип материала
Высшее учебное заведение

Список файлов книги

Свежие статьи
Популярно сейчас
Как Вы думаете, сколько людей до Вас делали точно такое же задание? 99% студентов выполняют точно такие же задания, как и их предшественники год назад. Найдите нужный учебный материал на СтудИзбе!
Ответы на популярные вопросы
Да! Наши авторы собирают и выкладывают те работы, которые сдаются в Вашем учебном заведении ежегодно и уже проверены преподавателями.
Да! У нас любой человек может выложить любую учебную работу и зарабатывать на её продажах! Но каждый учебный материал публикуется только после тщательной проверки администрацией.
Вернём деньги! А если быть более точными, то автору даётся немного времени на исправление, а если не исправит или выйдет время, то вернём деньги в полном объёме!
Да! На равне с готовыми студенческими работами у нас продаются услуги. Цены на услуги видны сразу, то есть Вам нужно только указать параметры и сразу можно оплачивать.
Отзывы студентов
Ставлю 10/10
Все нравится, очень удобный сайт, помогает в учебе. Кроме этого, можно заработать самому, выставляя готовые учебные материалы на продажу здесь. Рейтинги и отзывы на преподавателей очень помогают сориентироваться в начале нового семестра. Спасибо за такую функцию. Ставлю максимальную оценку.
Лучшая платформа для успешной сдачи сессии
Познакомился со СтудИзбой благодаря своему другу, очень нравится интерфейс, количество доступных файлов, цена, в общем, все прекрасно. Даже сам продаю какие-то свои работы.
Студизба ван лав ❤
Очень офигенный сайт для студентов. Много полезных учебных материалов. Пользуюсь студизбой с октября 2021 года. Серьёзных нареканий нет. Хотелось бы, что бы ввели подписочную модель и сделали материалы дешевле 300 рублей в рамках подписки бесплатными.
Отличный сайт
Лично меня всё устраивает - и покупка, и продажа; и цены, и возможность предпросмотра куска файла, и обилие бесплатных файлов (в подборках по авторам, читай, ВУЗам и факультетам). Есть определённые баги, но всё решаемо, да и администраторы реагируют в течение суток.
Маленький отзыв о большом помощнике!
Студизба спасает в те моменты, когда сроки горят, а работ накопилось достаточно. Довольно удобный сайт с простой навигацией и огромным количеством материалов.
Студ. Изба как крупнейший сборник работ для студентов
Тут дофига бывает всего полезного. Печально, что бывают предметы по которым даже одного бесплатного решения нет, но это скорее вопрос к студентам. В остальном всё здорово.
Спасательный островок
Если уже не успеваешь разобраться или застрял на каком-то задание поможет тебе быстро и недорого решить твою проблему.
Всё и так отлично
Всё очень удобно. Особенно круто, что есть система бонусов и можно выводить остатки денег. Очень много качественных бесплатных файлов.
Отзыв о системе "Студизба"
Отличная платформа для распространения работ, востребованных студентами. Хорошо налаженная и качественная работа сайта, огромная база заданий и аудитория.
Отличный помощник
Отличный сайт с кучей полезных файлов, позволяющий найти много методичек / учебников / отзывов о вузах и преподователях.
Отлично помогает студентам в любой момент для решения трудных и незамедлительных задач
Хотелось бы больше конкретной информации о преподавателях. А так в принципе хороший сайт, всегда им пользуюсь и ни разу не было желания прекратить. Хороший сайт для помощи студентам, удобный и приятный интерфейс. Из недостатков можно выделить только отсутствия небольшого количества файлов.
Спасибо за шикарный сайт
Великолепный сайт на котором студент за не большие деньги может найти помощь с дз, проектами курсовыми, лабораторными, а также узнать отзывы на преподавателей и бесплатно скачать пособия.
Популярные преподаватели
Добавляйте материалы
и зарабатывайте!
Продажи идут автоматически
6390
Авторов
на СтудИзбе
307
Средний доход
с одного платного файла
Обучение Подробнее