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

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

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

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

Рассмотрим интеграл з ) /(х) лх ь и выберем в качестве узлов точки хо=О, хз ОД хз !. Квадратуриая формула иитерполяциоииого типа, построенная по заданным узлам, совпадает с формулой Симпсона 1 ! ( (х) Нх — — (( (О) + 4( (05) + ( (1)). б е 173 где х,е=(а, Ь] н с,— числа, у=О, 1,..., п. В отличие от предыдущего параграфа, не будем разбивать отрезок (а, Ь] на частичные отрезки, а получим квадратурные формулы путем замены ((х) интерполяционным миогочленом сразу на всем отрезке (а, Ь]. Полученные таким образом формулы называются квадратурнглми формулами интерполяционного типа.

Как правило, точность этих формул возрастает с увеличением числа узлов интерполирования. Рассмотренные в ~ 1 формулы прямоугольников, трапеций и Симпсона являются частными случаями квадратурных формул интерполяциониого типа, когда п=О, 1, 2, р(х) — = 1. Получим выражения для коэффициентов квадратуриых формул интерполяциоиного типа. Пусть на отрезке (а, Ь] заданы узлы интерполирования хь„я=О, !, ..., п. Предполагается, что среди этих узлов нет совпадающих, в остальном они могут быть расположены как угодно на (а, Ь]. Заменяя в интеграле (1) функцию ((х) интерполяционным многочленом Лагранжа (7» !(х) = й„(х) +и„+,(х), где Е„(х) — интерполяциопный многочлен для ((х), построенный по узлам х„х„..., х„и г„ь,(х) — погрешность интерполирования. Тогда получим ь ь ь 1 р (х) / (х) ссх .= ( р (х) еа (х) ссх + ( р (х) г„(х) ссх = а ь = ~~ сь! (хь) + ( р (х) г„ , [х)с)х Таким образом, погрешность ф„+, квадратурной формулы (2), (4) ранна ь ф„„= ) р (х) г„„(х) с(х, (8) а где г„+, (х) — погрешность интерполирования.

Вспоминая выражение для погрешности интерполирования (см. (3) из $2гл. 3) г„.х (х) = сп (х), с а и ( ь ( х ) ) (а+ 1)1 получаем ф„„= ) р(х)ш(х)г'" м(п(х))с(х. (9) (л+ 1)1,) а Отсюда приходим к следующей оценке погрешности квадратурной формулы интерполяционного типа: ь !ф„„)( "" ( р(х))со(х))с(х, (и+!)1,) (10) а где М„„= шах 11~""н(х)). Из формулы (10) вядно, что справедхе[а,Ь] лино следующее утверждение.

174 Заменим теперь в (5) функцию !(х) миогочлеиом ф(х) наилучшего среднеквадратичиого приближения первой степени. Согласно (18) иэ й 5 гл. 3 этот миогочлеи имеет вид 1 ср (х) = (! (! ) — ! (О)) (х — 0,5) + — (! (0) + ) (0,5) + ) ( Ц]. 3 Отсюда приходим к квадратуриой формуле 1 =з 1 Г (х) с!х — — !/ (О) + ! (0,5) + ! (1)), о ие совпадающей с (6). 2. Оценка погрешности. Получим выражение для погрешности квадратурной формулы интерполяционпого типа.

Представим функцию ((х) в виде Квадратурная формула интерполяционного типа, построенная по и+ 1 узлу х„х„..., х, является точной для любого многочлена степени и, т, е. если 1(х) — многочлен степени и и с,— коэффициенты, вычисленньсе согласно (4), то имеет место точное равенство ь л ) р (х) 1 (х) дх = '~~' с»1 (х»). » »=» Справедливо и обратное утверждение.

Т е о р е м а 1, Если квадратурная формула ь » ) р(х) ~(х)дх=~ д»1(х») (12) имеющие степень п, и вычислим интегралы 1» = ( р(х) сг»(х) »1х. а По условию теоремы справедливы точные равенства л )» =;~,' дгр» (х1). »=о Поскольку /гФ1, а=1, 10, ф»(хд = ~ получаем 1„=д„, я=О, 1,..., п. С другой стороны, согласно (4) имеем ь !» = ~р(х) ь» (х) =ем (х — х»)»ь' (х») » Таким образом, д„=с„, )»=-0, 1,..., п, что и требовалось.

3. Симметричные формулы. Для некоторых квадратурных формул оценка погрешности (10) является грубой, так как оиа не учитывает симметрии формул. Рассмотрим, например, формулу Симпсона 1 ~ (х) дх .= — Д ( — 1) + 4~ (0) + )' (1)) 3 -1 17з точна для любого многочлена степени и, то она является квадратурной формулой интерполяционного типа. Доказательство. Достаточно показать, что д,=с„, где с„ определены согласно (4), я=О, 1, ..., п. Рассмотрим многочлены <р»(х) =,, А=О, 1, „и, (х — х»)»»' (х») для функции ((х) =х'. В данном случае имеем п=2, ['"+о(х) =6, поэтому согласно (9) погрешность можно представить в виде 1 фа=~ ы(х)ах, -1 где от(х) =х(х' — 1).

Благодаря симметричному расположению узлов имеем ф,=0. В то же время правая часть неравенства (10), равная 1 ( ы (х) [йх = —, -1 отлична от нуля. Таким образом, оценка (!О) не является точной для формулы Симпсона. Квадратурная формула (2) назь|вается симметричной, если 1) и четно; 2) узлы расположены симметрично относительно середины отрезка [а, Ь1, т, е.

а+Ь о+Ь вЂ” — ха=х„а — —, (с=О, 1, ..., п(2; 2 2 (Пй является нечетной функнней относительна точки х„п= (а+6)(2. Имеем л(а-1 м (х) = (х — х„) !~~ (х — ха) (х — х„а), а=о откуда, учитывая условия (14), получим ч(3-1 ге(х) = (х — х„,) ~ [ 1(х — «„,)' — (хл — ха~а)ч!. л=е Пв) 176 3) с„=с„„, (с=0, 1, ..., п(2. Свойство (16) коэффициентов квадратурной формулы определяется не только симметричным расположением узлов, но и симметрией весовой функции р(х). Говорят, что р(х) — четная функция относительно середины отрезка [а, Ь[, если р( — +х) =р( — — х) (16) для всех хек[0, (Ь вЂ” а)/2).

Л е м м а 1. Если с„определены согласно (4) и и четно, то соотношения (16) являются следствием условий (14), (16). Доказательство. Покажем сначала, что ч о>(х) = Ц (х — ха) Следовательно, при любом / Л/2-1 12(хл/ +/) = / [ [ [Р— (хь — хл )'] = — ы(х„/2 — 2), 2 —.-2 т. е. функция ы(х) нечетяа относительно точни х,ы. Из формулы (18) следует также, что ы'(х) — четная функция относительно ТОЧКИ»п/2 ПОЭТОМУ ы'(х „2) = ы'(хл) .

(19) Рассмотрим тецерь разность ь с2 — сл ь —— ~ Р (х) ы (х) Р (х) 2(х, О (20) где ь ( р (х) (х — х„/,)л" 2(х = О л вследствие нечетности подынтегральной функции, необходимо доказать, что л 1л= 'Я с4(хь) =О. 2=2 Представим 1„в виде суммы 10/-]- 1,"1, где л/2-1 1'„'= '~ с*(ха — х„/2)"", 2 2 л 1," ,= '~~~ с» (ха — х,/,)" ". а=л/221 1ут 1 1 и (х)— (» — Ха)Ы'(Хь) (Х вЂ” Хл 2)Ы' (Хл 2) Учитывая (19) видим, что 1 хп хл-ь 1 Хь — Х р (х) ы' (»2) (» — ха)(х — хл е) ы'(хь) [(х — хл„)2 в (хп/2 — ха)'] откуда следует четность м(х) относительно точки х„/2.

Таким образом,иодынтергаль~ая функция в (20) нечетна относительно середины отрезка [а, ь], н, следовательно, интеграл равен нулю. Лемма 1 доказана. Покажем теперь, что наличие симметрии повышает точность квадратурных формул. Справедлива Т е о р е м а 2. Пусть р(х) — четная функция относительно точки (а+Ь)12 и пусть выполнены условия (14), где и — четное число. Тогда, если квадратурная фор/нула интерполяционного типа (2) точна длл любого многочлена степени и, то она точна и для любого /нногочлена степени и+ ]. Доказательство. Достаточно показать, что формула точна для многочлена ](х) = (х — х.„)"+', х„„=0,5(а+Ь).

Поскольку Из условий (14) получим л 1„'') = 'Я с«(к«н — к -«)"" «=«Мы или л(«-) «(«-« («) ««1 «ы 1« = ~~' с«((х«м х() = — ~~ с«-«(х« — х«(«) (=а «=-« Последнее равенство справедливо в силу того, что и — четное число. Таким образом, получаем л('«-« 1„= 1(') + 1(') = ~~ (с« — с„«) (х« — х«м) ". «=а Согласно лемме 1 имеем с,=с„,, Ь=О, 1, ..., п/2 — 1, т.

е. 1.=0, что н завершает доказательство теоремы 2. 4. Формулы Ньютона — Котеса. Численная устойчивость ква- дратурных формул. Формулами Ньютона — Котгса называются квадратурные формулы интерполяционного типа, построенные на равномерной сетке, когда х,— х«,=й, (г=1, 2,..., и. Различают два типа формул Ньютона — Котеса: формулы замк- нутого типа и формулы открытого типа, В формулах замкнутого типа х„=а и х„=Ь, а в форл(улах открытого типа хотя бы один из узлов х, или х„не совпадает с соответствующей граничной точкой отрезка (а, Ь). Для простоты изложения рассмотрим лишь случай формул замкнутого типа, когда х„=а+И(), 1=0„1,..., п, йп=Ь вЂ” а.

В случае равномерной сетки можно упростить выражения для коэффициентов квадратурных формул. В формуле (4) сделаем за- мену к=а+(1), 0(1«=п, Простые выкладки, которые мы опускаем, показывают, что в результате замены формула (4) примет вид с« = (Ь вЂ” а) Ь,'"), где л ((и) ( — )] / ) ( 1()) ( (( — )) ° ° (( и) Ц( (21) И (и — ((1( ), и,) ( — () й Отметим, что формулы Ньютона — Котеса с п~!0 редко используются из-за нх численной неустойчивости, приводящей к резкому возрастанию вычислительной погрешности.

Причиной такой неустойчивости является то, что коэффициенты формул Ньютона— Котеса при больших и имеют различные знаки, а именно при и) )10, р(х) =1 существуют как положительные, так и отрицательные коэффициенты. Остановимся подробнее на значении знакопостоянства коэффициентов для устойчивости вычислений. Рассмотрим квадратурную сумму « 1„= 'Я с«1 (х(,), «=о 178 где 67, = ч~~ с»6». »=а (23) Поскольку квадратурная формула (2), (4) точна для )(х) = — 1, имеем л ь э' с» = ') р (х)»(х, »=ь л т.

е. при р(х) )О сумма л 'Я с» =,И »=ь (24) ограничена числом М)0, не зависящим от н. Предположим теперь, что все коэффициенты с„неотрицательны. Тогда из (23), (24) получим оценку л л /67„/('~~ /с»/ /6»/='Я с»/6»/ =(шах /6»/) М, »=о »=-О ь~»~л которая означает, что при больших и погрешность в вычислении квадратурной суммы (22) имеет тот же порядок, что и погрешность в вычислении функции. В этом случае говорят, что сумма (22) вычисляется устойчиво. Если коэффициенты с„ имеют различные знаки, то может оказаться, что сумма л '~~~ /с»/ не является равномерно ограниченной по и и, следовательно, погрешность в вычислении 1„ неограниченно возрастает с ростом п. В этом случае вычисления по формуле (22) будут неустойчивы и пользоваться такой формулой при больших п нельзя.

Таким образом, если необходимо сосчитать интеграл (1) более точно, то имеются две возможности. Во-первых, можно разбить весь отрезок /а, 6) на несколько частичных отрезков и на каждом нз частичных Отрезков применить формулу Ньютона — Котеса с небольшим числом узлов. Полученные таким образом формулы называются составными квадритурными формулами. Они часто применяьотся на практике, хотя и не являются достаточно экономнчнымн, поскольку требуют многократного вычисления значений функции г79 Предположим, что значения функции 1(х) вычисляются с некоторой погрешностью, т.

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

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

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

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