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

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

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

Но тогда, если занумеровать их в по- 2 рядке возрастания х,СХЯ( ... хв, то получим: (х» — а)(хь — Ь)а>' (хь) = (хв » 1 — а)(х » — Ь)а>' (х ь. 1). Следовательно, (17) с<в> = с<в> Ь В-»4-1 т. е. коэффициенты при >(хь) и у(х„»41) будут совпадать. Полученные здесь формулы численного интегрирования были впервые найдены Гауссом. Поэтому мы будем называть их формулами Гаусса. Было бы невыгодно каждый раз, как нам нужно использовать формулу численного интегрирования Гаусса, заново находить еа„(х), вычислять корни уравнения е>„(х) = 0 и подсчитывать коэффициенты с[у>. )[ля случая р(х)= — 1 и отрезка интегрирования [ — 1, 1[ такие вычисления были произведены для различных и. Произвольный отрезок [а, Ь[ может быть приведен к отрезку [ — 1, 1[ простой заменой переменной интегрирования: х= — + Ь+а Ь вЂ” а 2 2 Приведем некоторые значения коэффициентов с)з> и абсцисс для формул численного интегрирования -а-1 аа 2~Фа(л!)а ! > ( !)+[(2л)![з(2п [ !) а (6) — 1 ! 1 х,= О.

—,са — — 1, >та= —,7'(5); 1 11> 1 2 ' ' 1 3 =х,=0,577 350 269 189 6258, — с!1~>= — с",>= —, 2 2 2' 135 — х, п=3: — х, =х,=0,774 596 669 241 4834, ха=О, — са = — сз = —, — са = —, Йз= . У (1); 1 !з» и> 5 1 аз> 4 1 и> 2 13' 2 9' з 15750 л=4: — х, = х, = 0,861 136 311 594 0492, — ха = ха —— 0,339 981 043 584 8646, — с'," = — са ' = 0,173 927 422 568 7284, 2 — с11> =---сана =0,326 072 577 431 2716, 2 2 1 <в> 3 472 375 262 численное диФФеРенциРовАние и интегРНРовАние [гл.

3 9 5! ФОРМУЛЫ ЧИСЛЕННОГО ИНТЕГРИРОВАНИЯ ГАУССА 20'3 а=5: — х, = х, = 0,906 179 845 938 6640 — ха = х = 0,538 469 310 105 6830, — с(,! = †, с(а' = 0,1!8 463 442 528 0945 — с(а" = — с~~! — — 0,239 314 335 249 6832 — са(а( = —,=0.284 444 444 444 4444 2 '225 1 Оо! 1 237 732 550 х,=Π— х, = ха = 0,932 469 514 203 1520 — х, х,=0.661 209 386 466 2644 — х,=х,=0,238 6!9 186 083 1970 — с(," = —,са(м =0,085 662 246 189 5852 — с1" = — с(Г" = О,!80 380 786 524 0693 2 2 — са(" = — сз"' = 0,233 956 967 286 3455 ! ((а( 643 Ч34 486 !50 ~ — х( = ха = 0,949 ха=ха 0 741 107 9!2 342 7596 531 185 599 3944 845 !51 377 3970, — ха —— х,= 0,405 — с, = —, са = 0.064 1 ги ! Го 2 2 ха=О 742 483 084 4348 †, са(и = †, с!(! = 0.139 852 695 744 6384 — саи! — са('! = 0,190 915 025 252 5595 — са(И = — =0,208 979 591 836 7347 2 1225 ((а( 470050192 !11 500 ~ 1(ак видно из этой таблицы, коэффициенты сДч и абсциссы х» очень громоздки.

Поэтому формулы Гаусса следует применять в тех случаях, когда требуется большая точность и значения функции при большом числе аргументов получить затруднительно. 4. Формула численного интегрирования Эрмита, Если взять ч качестве р(х) функцию р(х) = 1 (18) и в качестве отрезка интегрирования отрезок 1 — 1, +Ц то получим формулу численного интегрирования е! и Ых = ~ с~л!/(х!) +Й(Г). 1 ! -1 (19) Если х; являются корнями многочлена я„(х), ортогонального с весом ! произвольному многочлену а(х) степени ~а — 1, то Й(Т! )' 1 — хе обращается в нуль, когда г'(х) являетсч произвольным многочленон степени (2л — 1.

Как будет показано в глаге. б, в качестве жл(х) можно взять 1 1 вл(х) = —, Т„(х)= —,соз(л агс сов х), (20) т. е. многочлены Чебышева, о которых говорилось з предыдущей главе. Следовательно, х; = соз , (! = 1, 2, ..., л) (21 — 1) и 2л (21) ьз ел (х) гГх Тл (х) (х — х) н (х!) У 1 — хз Г (х — хг) Т„(хб У! — ~~ (22) Для с!"! можно найти числовые значения. Для этого произведем под знаком интеграла замену переменного, положив х= соя В.

Г!олучим: с!"! = соз лВ лв (соз 6 — соз ВН Т„(соз Вб Но Т„(соз 8;) = 3~ 1 — созе в! з!и В! а соз лв соя 6 — со 6 оо+ В!соя В+ Взсоз28+... +Вл !соя(л — 1) В. 2б4 численное диееегенцпеованив и интвггигование (гл. 3 5 51 ФОРмулы численнОГО интеГРиРОЕАния ГАуссА 265 6" и 6 бб=кво. сов 6 — сов 6, о Чтобы найти В, в предпоследнем равенстве, положим там 9=9|, 0,.+-", 9,-1-2'--', ..., 9|+(п — 1)— и воспользуемся соотношениями сов п (6|+ й — ) сов пб ! сов' пб ! и в|и пб| сов 0 — сов О» !о о сов 0 |о о в|и О» сов 1пб|+ 2ав) 2ит сов (01+ й — ~ — сов О| 2в1 соо(0|+ Д вЂ” ) — сов 6» П Получим: В +-В,совб;+В сов20|+...

+В„,сов(п — 1)0|= В +В,сов(9, -+ — ")+В,сов2(9|+ — )+ ... 2и! ... + В„, сов(п — 1)(8! + — ) = б, Во + В, сов ~0! + (п — 1) — ~ + В, сов 2 ~ 9, + ( и — 1) — „~ + ... 2 ! ... + В„, сов (и — 1) ~0, + (и — 1) — ~ = О. Складывая эти равенства, будем иметь: п в|и пб| в|и ив| о — |и б так как и-1 ~) сов!Г(0|+у — )= ~' 2сов!Г(8|+/ — )в!п — = Р-О 2 в|и — в-о п и-1 "~6',~в!НА~9|+('~+" ~ — в!Нй~б|+ "У ') Д=.

2 в|и — ' о п ~ в!НА~9! + 1 — в!п !в [0! — — 1~= О 2 в|и— и (и = 1, 2, ..., и — 1). так как левая часть последнесо равенства является четной функцией 0 и числитель ее обязан делиться на знаменатель нацело, ибо Тп(х) делится на х — хл Отсюда 266 численное ДНФФеРенци! свАние и интеГРИРОЕАние [гл.

3 Отс!ода с1,! "В~ Т„(х!) (23) Получили формулу численного интегрирования Ф! а .!. 1 — у( !)+, ! ' У(х) сРх а У У1за) (1) [' Т„(х) ах У1 — хь п аа! ' ' (2п)1,/ 2!а-з )' 1 — хе — г(х!)+ — ! !.=1 -1 2! — 1 х. = соя —;-а,!. 2л (24) Остаточный член мо!кет быть упрощен, если воспользоваться результатами главы 5: (25) 5.

формулы численного интегрирования Маркова. Если пользоваться классификацией, приведенной в предыдущем параграфе, то формулы численного интегрирования Гаусса следует отнести к формулам открытого типа. так как концы отрезка интегрирования не принадлежат к числу узлов. А. А. Марков рассмотрел формулы численного интегрирования, лля которых )с()) обращается в нуль, когда Т(х) является произвольным многочленом степени (2и — 2 при дополнительном требовании. что либо х, = и, либо х„= Ь, и формулы численного интегрирования, для которых Й(Т) обращается в нуль, когда Т'(х) является произвольным многочленом степени < 2п — 3 при дополнительном требовании. что х, = а и х„= д.

Рассмотрим сначала первый случай. Обозначим через 8„(х) много- член степени л со старшим козффициентом 1, корни которого равны искомым узлам. Рассуждениями, как и при выводе формул Гаусса, покажем, что необходимым и достаточным условием для того. чтобы В„(х) удовлетворял поставленным условиям, будет ь ~ р(х) 9„(х) у(х) ь(х= О, (26) а где п(х) — произвольный многочлен степени <и — 2. Отсюда ь ь ~ р(х) [9„(х) — а!а(х)[ д(х) !1х= ~ р(х) 8„(х) п(х) а!х— а а ь — / р (х) а1„(х) и (х) Ух = О. а Эта формула является частным случаем формул Гаусса и носи! на- звание 99ормулы Эрмитп. 9 51 ФОРмулы численнОГО интеГРНРОЕАния ГАуссА 267 Поэтому 8„(х) — о1а(х) может лишь постоянным множителем отличаться от о1„1(х): 8„(х) — „(х) = а„„, (х), Для получения постоянной а„положим в этом равенстве х=а.

Получим: еи (а) -1 (')" Итак, 8„(х) = о1„(х) + а„о1„1(х). (27) И в этом случае 8„(х) имеет а действительнык корней, один из которых равен а, а остальные расположены между а и Ь. Доказательство полностью совпадает с тем, которое проводится в главе 5 ала о1а(х). Для отыскания коэффициентов А)~ ~ формулы численного интегрирования возьмем в качестве 7(х) функцию х — х где х',.— один иа корней уравнения 8„(х)=0. 7(х) — многочлен степени а — 1. Поэтому ь а )" Р (х) У(х) ~х =,Уа МА'7 (ХА) = и',гв 8„'(х;) а А 1 и о 7)1а1 / р (х) З (х) ах / (х — х,) 6„(х,) Остаточный член будет иметь вид: ь Г11 -Н(П /' 81,(х) (29) а Аналогично можно построить формулу численного интегрирования ь а о ума-Н (Р) ° ф'„(х) р (х) у (х) г(х = ~~~ Е',аЧ(ха)+ —,, ~ р (х) " ах.

(30) а 1 1 а где х» — корни уравнения ф„(х) =0; фа(х) = о1а(х)+ Ра(аа 1(х), Р„= — ", (31) ь Е1„> / р (х) фн (х) (32) (х — х,) ф (х,) 268 численное дивьвгвнциговлнив и интвггиговлние [гл. 3 Наконец, в последнем случае узлы х',« будут корнями уравнения (я(х) = О, где 44 (Х) 44« (а) аа (Ь) ",„(Х) = 41„,(Х) ш„1(а) «„1(Ь) 1 ! . (33) ~„-1 (а) ~ч-1 (й 44ч 1 (а) шч 2 (Ь) 44„2 (х) 44„, (а) ш„г (Ь) Коэффициенты Г(сч формулы / р(х) Г(х) 4(х= ~ Р~"~~(х; )+Ге„()) (34) определятся из равенств (35) и остаточный член будет иметь видь ь )41 (4') =, ' 1 " 4(х. Г42«-г! (() /' р (х) (~ (к) (2л — 2)1 ./ (х — а) (х — Ь) (35) Приведем значения х',." и Р)") для некоторых значений и при р(х) = 1 и а = — 1, Ь = 1.

хг = О, х,= — 1, Хг= 1 4 ~2 3' 1 3' 1 Рг= — ' х, =- — 0,2, Х2=02 х,= — 1, 1 ГЧ1 = —, б' б Рг= б 1 Г = —; б' 3 х = — —, 2 7 ' а=5: х,=О, Хь 49 64 ~З ' ~4 99' гч1 = 0,1, В заключение приведем вычислительный пример на мул Гаусса и Маркова. Прим е р. Вычислить по формулам Маркова и взяв п = 5.

(Ы ( ( (х) р (х) 41х 1 1 (1 «4~„4~ 1 4(Х 1+хг ' ь 3 — х =1, 7 ' ь 49 Рв — 0,1. применение форГаусса интеграл 9 б) ФОРмулы численнОГО интеГРНРОИАния чевышеВА 269 Преобразуем интеграл к промежутку [ — 1, -) — 11. Получим: с'х 7=2 / + —,. — 1 В формуле Гаусса: /(х,) = 0,24945107, /(хД = 0.23735995, /(х) = 0.2, /(х ) = 0,15706261. /(х )=0.13!00114. Вычисления лают /ж0,78539816. Все знаки верны. В формуле Маркова: !(х1) = 0,25, /(ха) = 0,24276181, /(х ) = 0,2, /(х1) = 0,14841465, /(х ) = 0,125, Вычисления дают /ж 0,78539214.

ф 6. Формулы численного интегрмрования Чебышева 1, Построение формул. В предыдушсм параграфе мы получили как частный случай формул Гаусса формулы численного интегрирования Эрмита. Они характеризуются тем. что все коэффициенты при /(х;) равны. Это оказывается существенным, когда значения /(х;) подвержены случайным ошибкам (напримсо, получены из эксперимента). Тогда выражение с, /(хг) + с,/ (хя) + ...

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

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

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

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