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

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

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

Представляем 7(х) в виде )'(х) = р(х)+й(х), (3) где р(х) — интерпеляционный многочлен, а Й(х) — остаточный член. 'Тогда ь ь ь ) р(х)7(х)йх= ~ р(х)~(х)с(х+~ р(х)Й(х)дх. (й) Первый член справа будет давать формулу численного интегрирования, а второй — остаточный член этой формуты. Интерполяционный многочлен >о(х) можно предатавить в виде у (х) = У (хо) ФО (х) + У (х1) Ф1(х) + ... + У (х„) Ф„(х). (5) Будем предполагать, что интегралы ь ~ р (х) Ф> (х) с(х = с; мы умеем вычислять точно. Они не зависят от функции>" (х). Поэтому их можно вычислить раз и навсегда и использовать лля вычисления интегралов ь ~ р (х) у (х) 1>'х а при произвольных >'(х).

Сама формула численного интегрирования будет иметь вид ь ~ р(х)у(х)ах=соу(хо)+с>1(х>)+ ... + О„У(х„). (7) а При численном интегрировании (а также и при численном лифференцировании) можно использовать интерполяционные формулы с кратными узлами. Тогда формула численного интегрирования примет вид ь ~ р (х) у (х) с(х со > (хо) + со>'(х>) + ... + с„> (х„) + + с),'> г" ' (х ) + с<,'> >' (х) + ... + с~>у' (х„) + +С>"-1>утч-1>(Х )+С>а-1>~>а-1>(Х )+ ... + С> -1>ут"а-1>(Х ). (8) о О 1 ! ''' а" а ' Специальным выбором узлов хь иногда удается добиться того. что часть коэффициентов с(у> обратится в нуль. Мы не будем пока входить в подробности этого случая и ограничимся формулами вида ~ р(х)> (х)>>х= ) сс> (хг)+ЙЩ.

(9) >-0 Остаточный член этой формулы обращается в нуль, если в качестве /(х) взять любую нз функций Ро(х), тд(х), ..., Ра(к). Учитывая последнее заме1ание, ны можем встретиться со случаем когда он обратится в нуль н для ОекотоРых дРУгих фУнкцнй Уа+1(х), Ра+з (х), ..., т„, (х), таких, что й>38 численное ЛНФФеРенциРОВАние и иптеггиговьние [гл. 3 239 % 3] злдлчл численного интегвивовлния %'[те.

Ть, ..., ТВ] ~ 0 из ]а, Ь]. Тогда, согласно второй главе (см, (3) е 5 гл. 2), запишем /(х) в зиле /(х) =~~~ ауру(х)+ ~ К(х, з)й +,[/(з)]лье (т) и). (10) У-О а Умножим обе части равенства на р(х) и проинтегрируем в пределах от а до Ь. Получим; Ь ьн Ь 1 Р(х)/(х) ь(х = ) а ° ~ Р(х) у ° (х) ь(х+ а О а [ л + ] Р(х) ~ К(х, з) йьиьь[/(л)]аьл ~аьх. а а Но / р(х)уу(х)пхаи ~ сьчу(хь) (/=О, 1, ...', т), !=О а Первый член справа даст приближенное значение интеграла, а остальные— остаточный член. Остаточный член может также быть записан в виде: ь [ ь )([/] = — ~ р(х)] / К(х, з) А „[/(з))льл ) Нх+ а ]В В Ь + У сь / К(хь з) й +ь [/(л)] ь(з.

(13) т О ВЬ Полусумма двух представлений остаточных членов ласт иам ь Й[/]=~О(з)й [/(з)]ил. а Позтому ЬВ т ] В (и ВЬ Х ау ~ р(х) ьгу(х) ь(х = ~ ~«у ~ ~~~ стьг ° (хь) = 1) сь ) а уу (хт) = О а 1-О 1-О ь-О 1-О В В ст/(хь) — ~ ст ~ К(хь, л) й .„]/(з)] ь(з. 1-О ь О а Итак, Ь В ь ] Р(х)/(л) тих= ~~~й сь/(хг)+ ~ Р (х)( ~ К(х, 3) й +ь[/(з)]аьз )ь(х а ь О а [а и сь ~ К(хь, з) 1 +ь [/(з)]ь(л. (12) 240 численнОе диФФеРенциРОеАние и интРГРиРОЕАние [гл. 3 где Ь ч 20(х) = ~ Р (х) К(х, г) в1ип (х — г) из — ~ с,К(хь з) з(яп (хг — з), (15) г 0 На этом мы закончим изложение общих методов численного интегрирования и перейдем к более подробному изучению формул, получающихся при использовании интерполирования алгебраическими многочленами.

5 4. Формулы Ньютона — Котеса 1. Вывод формул. В этом параграфе мы рассмотрим формулы для приближенного вычисления интегралов ~ г'(х) с(х. которые получаются путем замены подынтегрального выражения интерполяционным многочленом Лагранжа с узлами, разбивающими промежуток интегрирования на равные части. Эти формулы носят название формул Ньютона — Котеса. Пусть узлы интерполироаания расположены так: хг=а+1)г (1=1, 2, ..., л). (2) с = а+ (1 — й) (г, Н = а+ (и + )г) Й. (3) Для формул открытого типа й=1, а для формул замкнутого типа н=б.

Обозначим Р(у)= у'(а+(гу). Тогда ~ ) (х)с(х=)г [ Р(у)г(у, (4) г — А Здесь а либо совпадает с с, и тогда будем предполагать, что г( = а+(и+ 1) й, либо а+А = с, и тогда предполагаем, что а+п)г=Н. В первом случае узлы интерполирования не содержат точек с и г(, а промежуток интегрирования разбивается этими узлами на и+1 равных частей. Во втором случае концы промежутка интегрирования являются узлами интерполирования и промежуток интегрирования разбивается узлами на и — 1 равных частей. Формулы численного интегрирования, которые получатся в первом случае, будем называть формулами открытого типа, а во втором случае — формулами замкнутого типа.

Чтобы не проводить рассуждения дважды, положим 241 з 41 ФОРМУЛЫ НЬЮТОНА — КОТЕСА Интерполяционный многочлен Лагранжа, построенный для функции у(х) по узлам х), при такой замене независимого переменного перейдет в (у — 1) (у — 2) ... (у — ) + 1) (у — Š— 1) ... (у — и) ). (у) = ). (Š— 1) () — 2) ... (1) ( — 1)... (à — ) Г ( )' (5) Таким образом, ~ ~(х) е)х= И ~ Р'(у) с)у = И ~ ).(у)с)у-+ в 1-Ь 1-Ь .+ И ~ (у — 1)(у — 2)... (у — п)гт(у) 1; 2;...; п)пу= 1-Ь в в -1-Ь И ~т Р,)) ~' ( 1)в- (.у — 1)(у — 2)" (у — п),2у+ () — 1)1 (и — ))1 (у — 1) 1-1 1-а + И ~ (у — 1)(у — 2)...

(у — п)тч(у; 1; 2; ...; п)пу= 1-Л = (11 — с) ~ 1(:ЛьУ(п+ )И)+ + И / (У вЂ” 1)(У вЂ” 2)... (У вЂ” п)Р(У; 1; 21...; п)г)У. (6) 1-Ь олесь через 1в), обозначены выражения (в) веь )ь а — ° 1 йо ( — 1)в ' 1 (у — 1) (у — 2) ... (у — в) (и — 1 + 2Л) (Š— 1)1(п — Е)1,/ У вЂ” 1' )у. (7) 1-Ь Они не зависят от промежутка интегрирования и могут быть вычислены раз и навсегда, Кроме того, вычисления облегчаются благодаря тому, что )1,Ь=) 1+1,Ю (в) )в) (8) т.

е. разноотстояц)не от копцов коэффициенты формулы Ньютона— Котеса равны. В самом деле, в+а ) (в) ( — 1)'-' (у — 1) (у — 2) ... (у — и) " '+Н Л (и — 1-)- 2И) (п — 1)1() — 1)!,) (у — и -)-1 — 1) 1'А 242 числвннов диеезгннцигозлниг и интвгниновлниз [гл. 3 Заменяя под знаком интеграла у на и — г+ 1, получим: 7)в) ( — 1) 1 (в — г)(п — г — 1)... ( — г-).1) (н — 1+ 2а) (и — ))1(т — 1)!,/ и+л (н — 1+ 2а) (л — ))1 (1 — 1)! „)' г — 1 а'г=у,,л, ( — 1)" ' / (г — 1) (г — 2)... (г — и) (в) 1'л что и требовалось доказать.

С возрастанием н коэффициенты 1, л становится все более и (в) более громоздкими. Как было показано Р. О. Кузьминым, (7(нв)-и ~ с возрастанием т неограниченно возрастает, Так как, с другой стороны, (9) то среди 7; о должны иметься значении различных знаков. При(Ъв+ 0 ведем числовые значения 1)т"л длн различных (, А и и.

Каждый из коэфф~ци~нтов 7(т"'л является рациональной дробью, Длн сокращения таблиц мы будем брать знаменатели этих дробей при фиксированном и одинаковыми и эти общие знаменатели указывать в последнем столбце. В предшествующих столбцах будут даны только числители. л = 0 (формулы замкнутого типа): 4 ' 2 3 ~ 4 ) 5 ~ 6 Знамена- тели — 4 540 5 778 — 260 550 427 368 2 3 4 5 б 7 8 9 10 11 1 1 4 1 3 7 32 19 75 41 216 751 3577 989 5888 2 857 15 741 16 067 106 300 12 50 27 272 1 323 2 989 — 928 10 496 1 О 80 19344 — 48 525 272 400 2 б 8 90 288 840 17 280 28 350 89 600 598 752 ФОРмулы ньютонА — котесА Л =1 (формулы открытого типа): 2 3 Знамена- тели — 2 459 — 1 711 — 55 070 17 085 616 67 822 9 072 8 891 258 7 257 б 00 2.

Остаточные члены формул. Исследуем теперь остаточные члены формул Ньютона — Котеса. Как мы видели, они имеют вид: Л„л(У) =11 / (у — 1)(у — 2)... (у — п)1 (у; 1: 2: ...; п)г(у. (10) Преобразуем интеграл, стоягций в правой части. Возьмем сначала и = 2т — 1 и рассмотрим Згл-1+1 рл -1= ( (У вЂ” 1)(у — 2)... (у — 2т+.1) Х 1-Ь Х Г(у' 1; 2;...; 2т — 1) 1(у, (11) Введем вспомогательную функцию <р(х) = ) (у — 1)(у — 2)... (у — 2т+ 1)г(у, (12) 1-Л Очевидно, 18(1 — )з) =О. Точно так же и у(2т — 1+)л) =О. В са- мом деле, Лм-1+А 1р (2т — 1 + и) = / (у — 1) (у — 2)... (у — 2т + 1) 1(у, Произведем замену переменных под знаком интеграла, положив 1 2 11 11 611 460 1 787 4 045 2 752 447 — 1 1 — 14 — 453 — 954 — 2 803 — 11 690 — 6603 199 26 562 2 196 4 967 33 340 — 15 673 880 2 3 24 20 1 440 945 4 480 '244 численное дияавгенпигованив и интеггигование [гл.

3 Тогда 'а(2т — 1+й)= — ~ (2т — з — 1)(2т — г — 2)... ( — з+1)~Ы= а -1еь ат-ььь = ( — 1) + ~ (г — 1) (г — 2)...(г — 2т+ 1) 2г= а-ь = — ~р(2т — 1+ й). Отсюда и следует утверждение. Покажем далее, что е(х) нигде не обращается в нуль на интервале (О, 2т). Для этого исследуем сначала подынтегральную функцию ф(у) =(у — 1)(у — 2)... (у — 2т+ 1). (13) На отрезке (О, 2т1 она обращается в нуль в точках 1. 2, ... 2т — 1 и только в них. Она меняет знак при переходе через эти точки.

Далее, ф(2т — у) =(2т — у — 1)(2т — у — 2)... (2т — у — 2ги+1) = =( — 1)" Ф()) = — Ф(у). (14) е. график этой функции центрально симметричен относительно точки у = т. Покажем, что абсолютные величины интегралов 1+1 (е= ) ф(у) ау убывают, когда 1 возрастает от нуля до т — 1. В самом деле, ееа („,= ~ ф(у) Ьу.

!е1 Произведем замену у=а+1. Тогда 1.ь1 еч-1 У,~,= ( г(г — 1)... (г — 2т+2)йг= г л!у. уФ (у) 1 в Так как ф(у) не меняет знака на отрезке (1, 1+ Ц, то Й-1 Е /' Е '+г Š— 2т+ 1 / т (У) .! Š— 2т+ 1 где 1с. Еч, 1+ 1. Но 1- Е т Š— 2т+!( 2т — 1 — т+! = 1. 241» з 41 ФОРмулы нъютонй — котвсй Итак, ~/е~>~/1!>... >1/. 1~ ~/0~=~/2»1 1~, ~/1)=)/2т 2), так как 1р(1 — й)=22(2лй — 1+й) =О, и — Р(у; 1; 2; ° ..; 2яг — 1) = Р(у; у; 1; 2; ...; 2лг — 1). В силу знакопостояиства 1/(у) можно применить теорему о среднем. Поэтому М»-1ьй рй„, 1 — — — /'(Е; Е; 1; 2; ...; 2т — 1) ~ 1/(у) 1/у (15) 1-й или после замеиы разделенной разности производной яа-1~-й »чз»О (Е) /' (2»2)1,/ 1'й (16) Рассмотрим теперь случай четного а, я = 2лй. При этом Ъ„= ~ (у — 1)(у — 2) ...

(у — 2лй)Г(у; 1; 2; ...; 2яг)2(у. (17) 1-2 Разобьем последнии интеграл нз сумму двух интегралов: 2»»»2-1 рз„,= ~ (у — 1) ... (у — 2т)Р(у; 1; 2; ...; 2лй)О1у+ 1-й 2~»-~-й + ( (у — 1)... (у — 2лй)Р(у; 1; 2; ...; 2т)2/у=51+82. (18) 2~».12-1 Отсюда и следует, что 1р(х) + 0 при х~(1 — й, 2в2 — 1+й). Вернемся к исследованию рй О Произведем в (11) интегрирование по частям. Получим: рй„,— — 22(у)/'(у; 1; 2; ...; 2т — 1) (,~„~+ й~а-1+й — / 1/(у) — Р(у; 1; 2; ...; 2лй — 1) 2(у= 1/У 1-2 2М-1+й — </(у) Р(у; у; 1; 2; ...; 2т — 1)1/у. 1-2 246 числвнноз диеевгвнцигованив и интвггигованив (гл. 3 Заменим в Ъ1 произведение (у — 2т) Р (у; 1; ...; 2т) на равное ему, з силу определения разделенных разностей, выражение 1"-(у; 1; 2; ...; 2лг — 1) — гт(1; 2; ...; 2т).

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

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

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

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