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

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

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

3, 4,..., будем последовательно получать все В„(~). При этом непосредственно видно, что все Ва(!) будут многочленами. Многочлены Бернулли обладают двумя характеристическими свойствами. Рассмотрим Р(1+1, х) — Ч (Г, х). С одной стороны, зта разность равна хо(!+1! а хоаг А Ааг хо!о еа — 1 еа — 1 И 2! — =хост = х+- — + — + ..., с другой, Ва «+ 1) — В„(!) Х" л! а-о Приравнивая коэффициенты при х", имеем: ЛГ = Ва«+1) — Ва(Г), (4) Это и есть одно из характеристических свойств В„(1).

Продифференцируем теперь (1) по 1. Получим: О> / или Приравнивая опять коэффициенты при х", будем иметь: В„' (!) = пВ„г (1). (5) Это — второе характеристическое свойство многочленов Бернулли. Свойства (4) и (5) з свою очередь определяют В„(Г). В самом деле, на основании формулы Тейлора п1" = В„(Е + 1) — Вз (1) = ~ —, В~„~ (1).

Используя (5), получим: ВГяа)Я п(п 1) (и А +1)В а(!) Следовательно, ~аСвВ г (1) а-г Снова получили (3), однозначно определяюшее В„(!). Рассмотрим еше ряды х 1 — +- — х = ~) Ь„х", ел — 1 2 в О х = ~ Р„(1) х". а'г — 1 %т ;1 г~ ° (6) (7) Сравнивая (6), (7) с (1) и (2). видим. что Ь„= —,, если п ~ 1, в„ в„(~) — в„ Изучим свойства Ь„и Р„(!). Заменим в (6) х на — х. Получим: слева х 1 хз~ 1 — — х= — — х= е л — 1 2 лт 1 2 =х+ " — —,х= +- — х х 1 2 286 численное диееегенцигование и интеггиговлние [гл. 3 287 % 81 ФОРМУЛА ЭйЛВРА и справа СО 171( и и-а Следовательно, все Ь„с нечетными индексами равны нулю: (8) оы-1 — — О Найдем теперь Р'„Я: в'„ (г) В„ ,(г) Р'„(1) = "„, Таким обрааом, для всех п > 2 Р'„(1) = Р, (1) -+ 11 При п= 2 получим Р„'(Г) =В,(Г).

Иа (3) следует 21= 2В,(г)-+Ве(г), В,(г) г 2 Итак, (1 0) 2' Положим в (7) 1= 0. Получим: О = ~ Р„(0) х". А=О Следовательно, Р„(0) = О. (11) Полагая в (7) 1=1, будем иметь: х= ~~~Р„(1) х . Следовательно, (12) Р„(1)=0 при п + 1. Положим еше в (7) г= —, получим: 1 2' Заменим здесь х на — хг 1 — х = ~~~~( — 1)" Р„(2).х". 288 численное диефввенциеование и интеггиеовлние [гл. 3 Вычитая друг из друга последние два равенства, будем иметь: 1 — е з + 1 1 — ж — Ж е — 1 е — 1 з 3 х, +х е' 1 — х е х — ! ее — 1 3 + х = х = 2 ) ' Р,в, ( — ) хзь '.

е-г Итак, Р,( — )= —, Реь,( — )=0 при й~~2. (13) Покажем теперь, что при й) 1 Рзь~,(1) нигде не обращается в нуль 1 на отрезке [О, 1[ кроме точек 1= О, —. 1, а Рвам нигде не обращается в нуль, кроме 1= 0 и 1. Из (10) следует а из (11) С= О. Следовательно, для Ре(1) наше утверждение справедливо. РзЯ будет многочленом третьей степени и в силу (11), (12), (13) обращается в нуль при 1=0, —, 1. Других нулей Р,(Г) 1 2' иметь не может и для него утверждение также справедливо. Допустим теперь, что утверждение справедливо для Рн,,(Г). Положим, 1 для определенности, что Ргл,(Г) ) 0 при 0 <1< — и Рц,,(г) < 0 при —, < ~ < 1.

Тогда в силу (8), (9), (11) и (12): 1 2 Р' (Е))0 при 0<г< —, Р' (Г)<0 при 2 <Г<! ° 1 Рн (0) = Рзь(1) = О. Следовательно, Р л(1) имеет максимум при Г =— 1 2 и этот максимум единственный на [О, 1[. Ргл(Г) не обращается в нуль на [О,!) и имеет знак Рза 1(с), где ! некоторая точка от- 11 резка ~0, — з1.

Далее, 2з' эвз-1 ( ) за ( ) зв' 289 ФОРМУЛА ЭЙЛЕРА Так как Р„ '(1) на интервале (О, 1) обращается в нуль только в одной точке, то Рвь(1)+-Ь,в может обращаться в нуль на отрезке 10, 11 только в авух точках, а Р,ь„,(1) — только в трех точках. Слеао- 1 вательно, с=О, —,, 1 являются единственными нулями Ры,(т). 2' +г Так как знак Р (1) не меняется на отрезке 10, 11, а знак Р',,(() меняется, то й,л имеет знак, противоположный знаку Рег (г), и по абсолютной величине меньше, чем шах ! Ргл (1) ).

Поэтому 1 нуль Ре„.гг(т) при 1= — простой. Далее, Ры(0) =О, След 1вательно, Р„',,(с), а поэтому и Р ьг,(1) при $, близких к нулю, имеет знак (ггв, противоположный знаку Рге (ч), а следовательно и знаку Рн, г(1). Таким образом, Реве,(1) будут иметь чередующиеся по л знаки. Это же будет справедливо и для дев. 2.

Формула Эйлера и примеры ее применения. После этих прсаварительных рассуждений перейдем к вывоау формулы Эйлера. Рассмотрим а+А г'(х)гх, а где Г" (х) — некоторая, достаточное число раз дифференцируемая функция. Произвеаем замену переменных, положив х = а+ггг. Тогаа наш интеграл будет равен .( (х) г(х = й ~ / (а + И) Ж = Й )Г аг (г) Н, где ~(() =Г'(а+гГг). К послеанему интегралу применим правило интегрирования по частям: рг Ф(т)г(Г=Гр(П1 — а) ГР'(Г)ит=-Ф(П вЂ” 1 1~1 — —,) (Г)дт а о 1 1 1 1 1 /'г — -;~(1)+ —:Р(О) = — Р(О)+ — Р(П вЂ” 1 (1 — —,),р (1),11. 2 2 2 г,l ', а Первые ава члена здесь лают формулу численного интегрирования трапеций, послеаний член — поправку к ней. Воспользуемся теперь формулами (9), (10) н свойствами Р„(1) и (г„ алч преобразования 290 численное диееененцигование и интвгоигононие [гл.

3 последнего члена. Снова, интегрируя по частям, получим: л~ (( — —,, ) р'(() о[Г = / Р' 00 р'(Г) о[Г = [Рв (Г) р' (ГИ 1 1 — ~Рв(г) в" (1) (г= — ~[Р',(г) — Ь,,~ ро(1) а(= о О 1 = Ьо [и' (1) — у' (О) [ — / Р' (() сро (() оУ = Ь, [о' (1) — оо' (О) !— — [Р. (() ро'(1)[„'+ ~ Ро(1) р'о (Г) (Г = о г = Ь.

[со'(1) — ю'(0)[+- ~ Ро (() оо" (() оИ. о В последнем интеграле можно заменить Рз(Г) на Р;(Г) и еше раз повторить интегрирование по частям. При этом добавится член Ь, [~'н (1) — ~" (0)[, а вместо последнего интеграла будем иметь: ,(1) рбч(г) и. о Повторив наши операции г раз, получим: ! ,/ о(г) (г= —.' ,[г(0)+-о(П! — Ьв[л(1) — о'(он— — Ь, [ро(1) — ао'(0)[ — ... — Ьв, [р~ — ~(1) — ~ — ~(0)[+Я„. где Я,„можно записать в двух видах: 1 Ло, = — $Ре,, (Г) во'он (Г) Ю = / Р,„„о (Е) ф~'~н (Г) ЖГ.

о о Возвратимся к старым переменным и выразим Ь„через числа Бернулли. Будем иметь: оео д В,Ьо 1(х) (л = — — [1(а)+1(а+ й)! — — '; [1'(а+ Ь) — 1 (а)!— о — — [1о(а+А) — 1" ( )! — ...— 291 $8! ФОРМУЛА ЗЙЛВРА где 1 К„= — И +'~Р „,(1)у!1+П(а+ей)г(1= о ! =й Ф ) Р2„12(1)у'~ ~ю(а-+гй)Ж. о Если применить последнюю формулу к отрезкам (а, а-+И1, 1а+й, а+2И1, ..., (а+(и — 1)й, а+пй1 и сложить полученные выражения, то и получим формулу Эйлера: аеоо у(х)1(х=-й~,~ ~(а)+/(а+й)+/(а+2И)+ ., + о +/(а+(и — 1)И)+- — ((а+-пй)~ — — ', (Г'(а+ай) — у'(а)1— — 1)'п(а+пй) — ('п(а)1 — ...— — [/!" 0(а+пй) — Р (а))+йо„, (14) где й' „= — И ""' ~~~~ ~ Р2„~, (Г) 1 0 (а + И (Г + й)) о(1 = А-О О о — 1 1 = й"+' „5„~ Ре„ов (1) ~' +" (а + й (1 + й)1 211.

(15) А=-О О Остаточный член можно записать в другой форме, если воспользоваться тем, что Р,„„(1) не меняет знака нв отрезке 10, 11. Тогда 1 1122 пй 1 ! ) ~ Р2! о2(~) г(~' о где 1 — некоторая точка промежутка (а, а+ай). Но по (9) 1 1 1 „ () ~ . () Следовательно, 2.-.2() Во о 292 численное диееегенциговлние и интеггиеовлние !гл. 3 Поэтому остаточный член можно записать в форме ив жив Васев у!вг-~в!(() (2г -1- 2) ! (16) Б некоторых случаях об остаточном члене можно судить по самим вычислениям по формуле Эйлера. Так, пусть все производные нечетного порядка функнии 7'(х) иве!от одинаковые постоянные знаки на рассматриваемом отрезке и монотонны там, например монотонно возрастающие.

Тогда на основании второй теоремы о среднем нз первого выражения (16) получим, что знак Вз„ будет определяться знаком 1 Рв„в,(!) п1, Е Вв (х) = 1; В, (х) = х — —,; В, (х) = х' — х + —; 1 — 2 б' В, (х) = х' — — х'+ —, х; В, (х) = х' — 2хв+ х' — —; 5 4 5, 1 В (х) =ха — — 4+-, В (х) =хв — Зхз+ — х' — — хв+-г ° 5 12 1. в 2 2 42' 7 в 7 в 7 в ! В (х)=х — — х + — хв — — х'+ —.

х; 2 2 6 6 В (х)=х — 4х + —,х — —, х + —, хв — —, 14 в 7 4 2 в 1 з 3 3 3 30' = хз — — хз+ бхв — — ха+ 2хв — — х; 9 21 в 3 2 5 10 = х" — бхз+ — хз — 7х'+ бх' — -- х'+ —,, ° 2 2 66' в 42 в ! 5 691 7, — — Вю= —..! В!в= — ! Ввв= — ! 31!!7 5!О ' Вз (х) Вьз (х) В,= Вы= которые будут чередоваться вместе с г. Таким образом, остаточные члены будут иметь череду!ошнеся знаки и истинное значение интеграла будет заключено между суммой г и г+-1 членов формулы Эйлера. З,ля удобства пользования формулой Эйлера приведем значения чисел Бернулли и выражения многочленов Бернулли для некоторых значений л; 294 численное диееиввнциговлние и интаггиговлнив [гл.

3 Здесь как раз было применимо замечание, сделанное относительно оценки остаточного члена, когда производные нечетного порядка функции 7(х) знакопостоянны и монотонны. Пример. Вычислить сумму ~'., )з~, где р — целое положительное число. А 1 Положим в формуле Эйлера а= — О, и=1, 7(х)=хв. Тогда хгг(х= —, 0+ 1 +2г+ ... +(п — !)т+ — п-— ' ° 2 2 о — — рпг — — р(р — 1) (р — 2) пг Вт, В -а 21 41 Здесь ряд обязательно оборвется на каком-то шаге, так как производные с некоторого порядка будут равны нулю. Отсюда '+ лр ляе' лв рпя ' р(р — 1) (р — 2) р и р+1 + 2 12' 720 л-г р(р — 1) (р — 2) (р — 3) (р — 4) +- 30420 + Так, для р= 2, 3, 4, 5 получим следующие выражения: л в 1 л4+ в-1 л л (п+ 1) (2п+ 1) Зп~ пе (л+ 1)з 12 4 Вычисления дают: 0,050 000 000 0 0,090 909 090 9 0.083 ЗЗЗ 333 3 0,076 923 076 9 0,07142857! 4 О, 066 666 666 7 0,062 500 000 0 0,058 823 5294 0,055 555 555 6 0,052 631 578 9 0,025 000 000 0 0,693 771 403 1 — —, ° (0,01 — 0,0025) = — 0,000 625 1 12 — ° (0,000 1 — 0,000 006 25) = ! = 0,000 000 833 — 0,000 000 052 1 —,' .О,ООООО1(! — —.')= 252 ' ( 641 = — 0,000 000 003 9 — О,ООО ООО 01 (1 — — ) = =0,0000000000 ...

— 0,000 624 222 7 !п 2 = 0,693 147 180 4 Точное значение !п 2 = 0,693 147 180 5 ... 295 ФОРМУЛА ЗйЛЕРА П Х '=-'— и' аз 4пз 24п 1 Ь' = — + — +- — — — = —, п (а + 1) (2а + 1) (Зп'+ Зп — 1). 5 2 12 720 30 Ь 1 Х '=.' и" па 5п' 60пя п ь + + ' — (2иь + бп'-+ биз а) 2 12 720 12 В 1 Пример. Вычислить сумму ряда 1 1 1 1 НП + 1Оз + 105Я '+ !О7, + В формуле Эйлера положим Ь= 2, Г'(х) = —,, Тогда а=101, п= Оо. Э 1 11х 1 1 1 1 х + + ' + хз 2 1011 1031 1051 1О1 Вычисления дают В = 0„004 999 ЗЗЗ Зб.

СОз х пх = — сОз —, +. —, соз ( — — ).+ —, ~з!п ( — ) — з)п ( — — )~ Ь Ь Ь Ь Вяля . Ь . Ь1 — — (з1п( — ) — юп( — — )~+ + ... — ~" (( — 1) з1пЯ) — ( — 1)'з1п( — ~)~+ г . Отсюда Ь Ь йч ЬявВвь, Ь 2гйп — =Ьсоз — — 2 ~( — 1)в яп — +й . 2 2 а~а (2Ь)! 2 11 А-1 ь ь ь Деля на 2з!п —, и перенося — с1и — в л 2 2 2 1 — сги —,= Р ( — 1)1— ь ь ъч вь 2 2 л я' (2Ь)! л-о евую сторону, получим: Ь21 ~'11. 2 а!п —, и 2 П р и м е р. Используя формулу Эйлера, разложить в ряд по стех х пеням х функцию — с1и —, В формуле Эйлера полагаем 7"(х) =созх, а= — —, п= 1. ь 2 ' Тогда 296 численное диееееенциговлние и интеггиеовлние Ггл. 3 Это и есть искомое разложение.

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

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

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

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