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

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

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

+- с„/ (х„) (с,+с,+ ... +-с„фиксировзно) будет иметь наименьшую случайную ошибку при с,=ся=... =с„. В связи с этим П. Л. Чебышев поставил следу1ошую задачу: найти абсциссы х,, х,,..., х„и коэффициент К так, чтобы з формуле численного интегрирования ~ р (х) / (х) и'х = К 1) /(хг) + й (/) остаточный член )7(/) обрашался в нуль, когда /(х) является произвольным многочленом возможно большей степени.

Так как в нашем распоряжении находится и + 1 величин К, х,, ха, .... х„, то степень эта не меньше п. Коэффициент К находится без труда. Полагая/(х)= 1, получим: ~ р(х)п'х=КЛ. -1 (2) Как и в предыдущем параграфе, вместо того чтобы отыскивать сами абсциссы х,, хз,..., х„, найдем многочлен (3) м„(х) =(х — х,) (х — хз)... (х — х„). Возьмем в качестве 7(х) функцию У(.) =.' „. (4) где л — произвольное число, 1г~ ) 1. Тогда .~.! ~ Р(х)л» К~ ! + ! + + 1 ~ ~ )>( 1 ) (б) -1 Последнее выражение можно ааписать в виде (6) Проинтегрируем обе части равенства по щ Получим: +1 ~ Р (х) 1п 1 л — х ! !7х = К 1п — "" — + ~ й ( — ~ г(л.

-1 Здесь С вЂ” постоянная интегрирования. Потенцированием находим: 1 к .) — ! Р1х) Ы1е-х 1лх ю„(з)е~ ~ ' =Се (7) Представим показательную функцию, стоящую множителем в левой части равенства. в виде ряда по убывающим степеням щ Имеем: 1 1 1 х хт х" х" ь! х — х х а + хе + аз + ' ' ' + е" ч! +е"чз + ' ' ' г ! —— е Ряд справа при 1х! (1 и 1г~) и) 1 будет равномерно и абсолютно сходиться. Далее, (9) 270 числвннов дивввркнцировлнив и интвгрировлнив (гл, 3 Отсюда 272 численное диеевгвнцигование и интеггиговлнив [гл. 3 Таким образом, мы показали, что в„(г) в формулах численного интегрирования Чебышева определяется из равенства — р (х) щ Ы вЂ” е ~ ле г л,[ в„(г) = СЕе (17) Здесь С подбирается так, чтобы коэффициент при з" в м„(г) был равен 1. Дадим еще один способ получения ш„(г).

В формуле / р (х)у(х)(х = К,'5', у(х) + й (,у) возьмем У(х)=па+а,х-+пахе+ ... -[-а„х", (18) где ач, ам а,, .... а„— произвольные действительные числа. Полу- чим: ав / Р(х)Их+а, ~ Р(х)хНх+ ... +а„[( Р(х)х" г1х= =ааКл-[-а,К [х,+-ха-[- ... +х„[+ааК Х )~ [х;+х, "+ ... -[-х'„[+ ... +пвК[х",-[-хя+ ...

+ х.~. (19) В силу произвольности а; должны иметь место следующие равенства: рч =Кп, р, = К [х, -1- х, + ... [- х„[, ря = К [х,+ха+ ... -+ха[, (20) так как )с(У) в этом случае должен обращаться в нуль. Введем обозначение е1 рг= ~ Р (х)х'гЕх. -1 5 61 ФОРМУЛЫ ЧИСЛЕННОГО ИНТЕГРИРОВАНИЯ ЧЕБЫШЕВА 273 Первое равенство дает: К= — = — ~ р(х)»»х, — 1 т. е.

то же, что и раньше. Деля остальные равенства на К и обозначая — = т», получим систему и уравнений для определения х». Р» К х, +х, + ... +-х„=т,, з з з хз + хз + ... +. х„= т„ (21) х1 + хз + ° .. + хп — т Опять будем разыскивать многочлен з1п(х), корнями которого являются х,, хз...., х„, удовлетворяющие этой системе. Пусть т„(х) =х" +Ь,х"-'+ ... +Ь„,х+Ь„.

(22) Нам нужно найти Ь,, Ь,, ..., Ь„. Имеем: У пп (Х) .Ап» Х вЂ” х» (23) зз„'(х) = пхи-'+(п — 1) Ь,х"-'+(и — 2)Ь,х» з+ +Ь„, (24) Далее, -,"( ) =х" +(Ь,+х,)х" +(Ь,+Ь,х»+х,')х" + ... ... +(Ьп 1+Ьп-гх»+ ... +.Ь,х, +х"; ).

(26) Отсюда и из равенства (23) получаем: ...'„(х)=пх '+(ЛЬ,+т1)х" з+(ЛЬз+Ь»т,.+тз)х" з+ ... ....+(ЛЬ„1+ Ь„зт»+- ... +Ь,т„, +т„,). (26) ЛЬ, +-т,=(. — 1) Ь,. пЬз+Ь,т, +та=(п — 2)Ьз, (27) пЬ„,+Ь„,т,+ ... -+Ь,т„, 1 т„,— Ь„, Сравнивая коэффициенты при одинаковык степенях х в равенствах (24) и (26), находим: 274 числзнноз диеезгвнциговлнив и интвгвивованнз (гл. 3 Отсюда последовательно можно найти Ь,, Ьа, ..., Ь„ ,. Для получения коэффициента Ь„ рассмотрим сумму аа„(хг)-!-аа„(ха) + ... +ы„(х„).

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

Очевидно, ы„(х), полученное вторым способом, будет совпадать с ы„(х), полученным первым способом. рассмотрим теперь частный случай формул Чебышева, когда р(х)= — 1, При этом +1 К= — 1 е(х= —, 1 / 2 и,/ и' -1 (29) а абсциссы х,. ха, ..., х„являются корнями уравнения +1 — 1ч! а-е !ее а д Ее (30) Но при г )1 1п ~ я — х ~ е(х = (г+ 1)! п (я+ ! ) — (я — 1) 1п (я — 1) — 2 = -1 = 21п г+(я-1- 1)1п(1 -(- — )+(1 — г) !п(! — — ) — 2 = =2!пя.+(я+1) !с — — — +- ! —...~ — (1 — г) Х г1 1 2 за Зга Х ~ +2 а+3 а+ ..~ — 2=2!пав Отсюда л н а щ (З) — Сд~заа а а а а а г' ьыаа (31) С одной стороны, эта сумма равна нулю, так как аа„(ха) = О.

С дру- гой стороны, складывая представления аа„(ха) в виде многочлена по возрастающим степеням х; при различных ! получим: Дадим п значения 1, 2, 3, 4, 5, 6, 7, 9. Получим следующие урав- нениЯ для пОлучения хгг х' — — = О, 1 3 1 — х О, 2 1 3 х'+- — 45= О г ха+ х — 0 б 72 1 1 хв -)- — хг — = = О. 5 105 '7 П9 г 149 х +,' М=О, б + 360 6480 3 27 в 57 в 2217 — х' -!- — хв — —. хв -)- х = 0. ~ 2 40 560 22 400 (32) хв Можно было бы получить г1„(х) и для других значений л, но, как показал академик С. Н, Бернштейн, при этом уравнение в1„(х) =0 будет иметь комплексные корни и, следовательно, соответствующая формула Чебышева не может быть использована. Абсциссы формул Чебышева при различных значениях л и 1 даются следующей таблицей: и = 2 — х1 = хг = 0 577350' Л=З вЂ” х,=х,=0,707!07, ха=О; и = 4 — х, = х, = 0,794654, — х, = х, = 0,187592; и= 5 — х = х =0,832498, — хг= х = 0,374541, И= 6 — х,=хв = 0 866247, — хг= ха= 0 422519 — х, = х, = 0,266635; п = 7 — х, = х, = 0,883862, — х, = х, = 0,529657; — х,=х =0,323912, х,=О; л=й — х,=х =0,91!589.

— х,=х =0,601019; — х,=х,=0.528762, — х,=х,=0,167906, х =О. Рассметрим пример на вычисление интеграла по формуле Чебы- шева. Пример, Вычислить по формуле Чебышева интеграл 1 ! = ~ 1 г = — = 0,78539816 ..., в хг — — О взяв и = 7. 6 6! ФОРМУЛЫ ЧИСЛЕННОГО ИНТЕГРИРОВАНИЯ ЧЕБЫШЕВА 275 276 числвннов диеевгинцнгованин и интвггигованиа !гл, 3 Как и в предыдущем параграфе, преобразуем предварительно интеграл к отрезку [ — 1, 1[. Последовательные значения ординат будут таковы: У (х,) = 0,249159, У(х,) = 0,236898, 7(х,) = 0.224361, ,7(х,) = 0,200000, У(х ) = 0,173831, У(хв) = 0,157732, У (х,) = О, 132469. Вычисления дают 7 = 0,785400. Ошибка не превышает двух единиц шестого знака. 2.

Остаточный член формул Чебышева. Получим теперь остаточные члены формул численного интегрировании Чебышева. Ограничимся случаем р(х)=1. Пусть число ординат, использованных в формуле Чебышева, равно и. Соответствующий остаточный член будем обозначать гг„(у). Как известно, )Св (У) = О, когда У(х) ЯвлЯетсЯ пРоизвольным многочленом степени и.

Если и — четное число, то )ся(у) обращается в нуль и для произвольного многочлена степени не выше и+ 1. В самом деле, если У(х) = аз+ агх+ аьхт+ ... + а„х" + а„+гх"+1, )(а(У) =)Са(ае+агх+ ... +а„х")+Ля(а„ьтха+1), то Первый член справа, очевидно, равен нулю, а +1 а 2 %т А'„(а„+1х"+1) = ~ ая+тхя+гах — — г а х"+'. И аз и+1 Отсюда, в силу наших предыдущих рассуждений, Но интеграл, стоящий в правой части, равен нулю, так как х"+1 является нечетной функцией х. Точно так же обращается в нуль и сумма в правой части, так как корни уравнения ма (х) = О симметричны относительно начала координат.

Итак, Яя(ая+1ха+1) = О и Р„(у) = О. Утверждение доказано. Пусть теперь У(х) — произвольная функция, обладающая на отрезке — 1, +1[ непрерывными производными до порядка из+ 1 включительно. десь ис равно и, если и — нечетное число, и равно и + 1, если и — четное. По формуле Тейлора имеем: У(~)=у( — 1)+(х+1)у'( — 1)+ 21 уа( — 1)+ °" (х+1)Я ... -[- — Ура) ( — 1) -[- — у У(ян-1) (з) (х — з)т Дз. (33) (Х+ 1)ы 1 1" и11 ш( 3 -1 3 6) еозмзлы чнслвнного ннтвгзнзоваиия чввып1вва 277 или )7и(У) = — с1х ~ У(и+1! (з) (х — з)" 1(з— т! 1 — 1 и — — р — ! 1 («1+1! (з) (хз — з)~ с(з. (34) л лыз т1,/ 1 1 -1 Меняя порядок интегрирования в первом члене справа, получим: +1 я (1 — з) + 2 жч )!я (У) = — ~ У( +О (з) — а1з — — т / У(а1+1!(з)(хз — з)г«1(з. т1 т+1 т(л л а -1 1-1 -1 (35) Введем обозначение (х — з)т = (х — з)Я1 при х ~ з, (36) О при х(з.

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

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

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

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