Том 1 (1160083), страница 29

Файл №1160083 Том 1 (И.С. Березин, Н.П. Жидков - Методы вычислений (1962)) 29 страницаТом 1 (1160083) страница 292019-09-19СтудИзба
Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

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

=( ) у 2т 2 . Г 5!22х 2 ) 'рг2 2 $ 11. Общая задача интерполирования алгебраическими многочлеиамн 1. Интерполяцнонный многочлен Эрмнта. Рассмотрим теперь более общую задачу, чем та, которую мы решали. Пусть нам задана базисная система интерполяцнонных функций !рв(х), р,(х), в„(х), ... на [а, Ь[. Требуется найти такую их линейную комбина- цию ср(х) = ч~р с,рв(х), !то <р(хв) =у, <р'(х,) =у,', ..., ф" -и (хв) = у!"-'1, <р(х ) =у„, ср'(х„) =у'„, <р('» ')(х ) — у('» '), ~ (2) где у!р! — заданные числа, а хв~[а, Ы (1=О, 1, 2, и; х! ~ ХР при 1 ~ )). Так как число условий, которые мы накладываем Проверкой убеждаемся, что этот многочлен удовлетворяет поставленным условиям.

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

сов х, сов 2х, ., сов ах и чя х, 5!Н2х, ... 51пих первая иа [О, л), а вторая на (О, л) образуют системы Чебышева. 164 теОРия интеРпОлиРОВАния и некотовые ее пРилОжения (гл. 2 на у(х), равно ао+.а,+ ... +а„, то для того, чтобы наша задача всегда имела единственное решение, требуется, чтобы яо = ао + а, + ... + а„— 1 (3) и 9в (Хо) ~р '(х,) То (Хо) в', (х,) ~о, (хо) о,' (хо) ~р~," -П (х ) жоао-О (х ) то (хо) 'Рг (хд <р1 о-О (Х ) т (.'о) чь О.

(4) ~7( о г)(х). ~( о ~)(х ) <7( и ~)(х ) Н (х) — уя(х) = в„(х) Н „(х), в„(х) =(х — х,)(х — х) ... (х — х„). где (6) При любом мцогочлене Нв „(х) функция Нв(х) = Ь„(х)+в„(х) Нв „(х) (7) принимает в узлах интерполирования значения уо Подберем теперь Нв „(х) так, очтобы были выполнены и остальные условия. )1ифференцируя обе части равенства (7), получим: П (х) = б„(х) + в„(х) Н „(х) -~- в„(х) Н „(х). Полагая здесь х=хо будем иметь: Н (хо) = 7.„(хо)+в„(хо) Н „(хг). (9) Так как в'„(х,) чь О, то в каждой точке, в которой задано Н' (х,), мы найдем Н „(ха. 11ифференцируя еще раз, получим: Нв (Х) = 7 о (х) +.

ооя (Х) Нв „(Х) + 2во (х) Нв-и (Х) + +- в„(х) Н „(х). (10) Мы не будем здесь входить в подробности общего случая, а ограничимся лишь алгебраическими многочленами, т. е. положим ~ро(х) = х*. Итак, нам требуется построить алгебраический многочлен степени не выше т, удовлетворяющий поставленным выше условиям. Предположим, что такой многочлен существует, и обозначим его через Н (х). Наряду с Н (х) рассмотрим интерполяционный многочлен Лагранжа уо(х), принимающий в точках х,, х,, ..., х„значения ув у,, ..., у„. Разность Н„(х) — 7.„(х) должна быть много- членом степени не выше т, обращающимся в нуль в точках хв хп ха, ...

х„. Следовательно. й 11[ задача интвгполигования алгввгаичвскими многочлвнлми 165 Полагая снова х= хи найдем: Нио(хо) = ьи(хг)+ми(хо) Н и(х))+2о)и(хг) Нио и(х)). (11) . ни,-а) (а,-а) Нт-и (хо) = ао ° ° ° г)ио-и (ло) = ао l ,.И,-а) )и,-а) Н и(х,)=зп ..., и и (х,)=л)' Н и (хо ) = го; Н и(х,)=л,; (12) (" а) ( и а) Н „(х„)=ли; Н,„-и(х)=ли,..., Н "и (Х )=гии где а)й — известные числа.

К Н и(х) применим точно такой же прием. Получим некоторые условия, наложенные на Н и(х). В конце концов, нам потребуется построить интерполяционный многочлен Лагранжа по данным в некоторых точках х). Посмотрим, какова же будет степень полученного таким образом многочлена Нт(х). Эта степень будет равна числу узлов, в которых заданы уо плюс число узлов, в которых заданы у',, плюс число узлов, в которых заданы у,.

и т. д., плюс число узлов, в которых заданы самые старшие, входящие в условия, производные и минус единица. Таким образом, эта степень равна по+ а, + ... + аи — 1 = т, что и требовалось. Построенный нами многочлен единственный, который удовлетворяет поставленным условиям. Действительно, если бы имелось два таких миогочлена Нт(х) и Йт(х), то их разность Н (х) — Й (х) представляла бы собой многочлен степени не выше т, имеющий на отрезке [а, 1)[ т+-1 корней (с учетом кратности корней), что невозможно.

Да и сам процесс построения в силу того, что все Н „,(х) определялись единственным образом, дает основания утверждать единственность построенного многочлена. Приведем примеры на построение таких многочленов, которые мы в дальнейшем будем называть интерполяйионными мноаооленами Эрмита. Из этого равенства мы сумеем найти Нио „(хо) в тех точках, в коо торых заданы И (х). Продолжим этот процесс далее. Каждый раз коэффициентом при старшей производной от Нт и(х) в точках хо будет о)'„(х,). Таким образом, мы сведем нашу задачу об отыскании Н„,(х) к задаче об отыскании Н и(х), удовлетворяющего условиям: 166 тногия интавиолиговлния и никотогыв яв пгиложиния (гл.

2 Пример. Пусть значения 7"(х) и ее производных в точках х = О, 1, 2 заданы таблицей: 129 448 1344 Найти интерполяционный многочлен Эрмита. Интерполяционный многочлен Лагранжа 5„(х) в этом случае будет равен 7 ( ) 1(х — 1)(х — 2) ( 2 х(х — 2) ! 129х(х — 1) 63ха — 1 2 Отсюда Н,(х)=63х' — 62х+! +х(х — 1)(х — 2)Н,(х).

Дифференцируя это выражение, находим: Н; (х) = 126х — 62 + (Зх' — бх + 2) Н (х) + (х' — Зх'-! — 2х) Н, '(х). Подставляя сюда значения х= О, 1, 2, получим: На(О) = 31; Н4(1) = 57; Н (2) = 129. Вторая производная от Н,(х) имеет вид Н; '(х) = 126+(6х — 6) Н, (х) + 2 (Зх' — бх+ 2) Н,'(х) +- + (х' — Зх' -1- 2 х) Н' (х), Отсюда находим значения Н'(х) в точках 0 и 2: Н' (0) = 15; Н,' (2) = 111. Итак, нам надо найти многочлен Н,(х), удовлетворяющий условиям: Н (о) 31 На(1) 57 Н (2) 129 Н,'(0) = 15; 4( Записываем его в виде Н (х)=31( )( ) -)-57 (х ) -«-129 4 2 — 1 2 + х(х — 1) (х — 2) Н, (х) = 2 Зх' + Зх+ 31 +- (ха — Зхя -)- 2х) Н, (х).

й 111 злдлчл ннтвгполиговлння ллгввглнчвскнмн многочлвнлмн 167 Дифференцируя, находим: Н, '(х) = 46»+- 3 -)-(Зха — бх+. 2) Н, (х) -+ (х' — Зх'+ 2х) Н' (х). Подставляя сюда»=0 н 2, определяем Н)(0) = 6. Н,(2) = 8. Н,(х)=х+6 Отсюда Н,(х) =(63х' — 62х+1)+-х(х — 1)(х — 2)[2зха-)-Зх+31 ~- + х (х — 1) (» — 2) (х + 6)1 = 1 -1- х'. Н (х)=А„(х) +ы„(х) Н „(х), в„(х) = (х — хе) (х — х,)... (х — х„). (13) где Дифференцируя, найдем: Н' (х) = Е„' (х) +-ы„'(х)Н „(х) + ы„(х) Н'„ „(х), (14) Отсюда у,'.

= Е„'(х;) + вч (х~) Н „(х;), (15) у — Е„(»<) Н „(х,)= ч„(»,) Итак, Уг — Е„(» ) (»,) Н„, „(х) = а (»д (» а) "' (»г) =Х' (16) Пусть ч (Х) = У.щ (х). (х — х,.) ч„(» ) (17) Тогда интерполяционный многочлен Эрмнта можно вапнсать в виде Н,„(х) = ~~ у,1.„;(х) + ~~ ы„(х) , " ~ 7.гм (х). а О я 0 (18) Рассмотрим еще один пример.

Найти интерполяционный многочлен Эрмнта, принимающий в точках х,, х,, ..., »„значения уы у,...., у„н имеющий там прона- водные, равные у', у'„..., у'. В этом случае э 1Ц злллчл ннтзгполигования алгвзванчвскнми многочлвнамн 169 П,ри х=ха получим: ь,„(х) Рг(х) = ~ [А+В(х — хг)), Определим коэффициенты А н В. Полагая в последнем равенстве х= хо получим: гз 1 = в„(х,) А.

(23) Отсюда А= 1 э„(хг) Лифференцнруя равенство (23) и полагая затем х= хм найдем: 0 = Р' (хг) = а'„(хг) в„"(хг) А + м'„(хг) В. Отс юла  —— м (Х) га в„(ха) Коэффициент Рг(х) можно представить в другом виде; 1 — " (х — хг) ~ ~„(х,.) Ф 2 ~и( $) = Е„,(х) 1(1 — —,— (х — хг) ч(х) ю (Х) Рг(х) = (х — х )тм„(х,.) Итак э Ф ю Н~(х) =,')' уг ~1 — ", ' (х — хг)~ Е.'„,(х)+ '5' у,',(х — хг) Е'„,(х).

Окончательно получаем: Н (х)= ~~) ~Уз ~1 — ", (х — хг)~+У',.(х — хг)~ Е„н(х). (24) -„(-) 2. Общий вид интерполяцнонного многочлена Эрмита. Найдем теперь общий внд ннтерполяцнонного многочлена Эрмнта. Для этого ь„г (х.) Рг (хь) = Е' г (хь) — м' (х ) ~~~, Е у (ха) = т о "(д) =Е г(х„) — У.„г(хь)=0. (22) / Таким образом, Р,(х) имеет лвукратный корень прн всех х= ха, й~д Следовательно, этот многочлен имеет в качестве множителя н~, (х) э . Так как степень Рг(х) равна 2л+. 1, то его можно запи(х — хз)э сать в виде 170 теОРия интеРпОлиРОВАния и некОтОРые ее пРилОжения (гл. 2 построим многочлены Н; (х) степени не выше л«, удовлетворяющие следующим условиям: Н«(х„)=Н,' (х„)=...

=Н,'"д "(х,)=0, Н,~=1 (1=0, 1, 2, ..., в; /=О, 1, 2, ..., ૠ— 1). (25) Н,(х) =А«««+А~ ~(х — х,)+ ... +А(,.«у )(х — х«)"«« . (27) Для определения коэффициентов А„воспользуемся тем же приемом, (ы что и в предыдущем примере. Пусть («(х) = (х — х,)"' (х — х,)"' ... (х — х,) 'и. (28) Тогда А«"1.+А««'1(х х«)+... 1 А('««)(х х ) «« (х — »«)"«« НН (х) (29) Подставляя сюда х= х«, получим: (х — хД ' НО (х) з-ьа, ««(х) (х — »«)« (30) Первое отношение непрерывно прн х=х«.

Следовательно, ~(» — х«) «~ ~(» — »«) «1 Предел второго отношения найдем по правилу Лопиталя: Н«1 (х) 1 Н$ (х) 1пп = Йп ~~~, 1 (» — х«) 3 ~~~, /1 «1 Итах, А1«1 1 1 (» — »«) «) А 1 ] (3 1) Так как НН обладает в точках х,, х,, ..., х«п х,+и.... х„нулями соответственно кратности а,, а,, .... а;,, а,ьп ..., а„, а в точке х« нулем кратности 7', то Н« =(х — х,)" (х — х,)" ... (х — х«,)"«-«(х — х«)«Х Х (х — х«««)'' «... (х — х„)" Н««(х).

(26) где Й«(х) — многочлен степени ૠ— / — 1, не обращающийся в нуль прн х=х«. Представим его в виде в 11) злдлчл ннтввполигозлння ллгвзвличвскнми многочлвнлмн 17! Аналогично находим коэффициенты А(ы.. Ц, 1ь> 1 ль ~ (х — хг) а Н;д (х) А;, = — Иа ье Мх" 1 Я(х) (х — х;)У (32) Применим правило Лейбница для дифференцирования произведения лхл м (х) (х — х))( ~Г аа я (х) (х — хг)а Производные ~(х — хг)» ~ непрерывны в точке х=хо Поэтому ~ (х — хг) г ~ ~ (х — хг) а ~ 4 ~4 Для определения Нгу(х) ' (ь в) Иш аье (х — хг)а воспользуемся таким же приемом, как н для отыскания А, . Много- (В) член Ны(х) имеет степень не выше гв.

Он делится на (х — х,)т. Следовательно, его можно записать в виде Н„(х) = В(',-'(х — хг)'+ В(п(х — х;)"'+ ... -(- В(',"-'з(х х,.)™ нли . =В(',)+В~/(х — х,)+ ... +Вф "(х — хг)"-'. (х — хг)у Отсюда Нм (х) 1'~ '" = (Д вЂ” Р)1 Во-, (л я) (х — ху)а) г В нашем случае /+(з — р </+й <1+аг — ~' — 1=а — 1. Но В)", я), будучи коэффициентами разложения НН(х) по степеням (х — хг), записываются в виде езс1+ ь — в) В(ь-э) н~я (хВ (У+ Д вЂ” р)1 172 тзогня ннтагполнговлния н нзкотогыа ее пгиложения (гл. 2 в," ,= —.' .

(о) ! у! Итак ,„) ! Ль (х — х!)' Нг, (х) А!)) = — 1)ш Л! т т лх" ы(х) (х — хг)г 1 (х — хг) ро л)у! а( ) ) (ЗЗ) а -)-) 1 ц(х) кт 1 Г(х — хг) гт Фг т(й) Н;~(х) =-;, — (х — хг) . (34) ь (х — хг)» ь е ' (х) е т< Используя построенные нами функции Н)у(х), нетрудно написать выражение для Н,(х). Легко вндеть, что и ч Н„(х) = ч; ~ у()) Н, (х) г-о у-о Ц илн ч Г' «г-у-г , „)ю г-о у-я ь о (х) м а (х — ~д 3.

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

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

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

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