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

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

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

За приближенное значение второй производной при численном диф- ференцировании будет приниматься 222 численнОе диФФеренциРОВАние и интегрнРОЕАние [гл. 3 Таким образом, остаточный член в етом случае примет следующий ВИД1 + 2 "„ ~ (Х; Х; Х11 Х,; ...;Х„) +2аз„ (Х)у (Х; Х; Х; Хр; Х,; ...; Х„) Имл (х) (10) или Ф (х) У~ + ~(е) + И"' (х) У~ '61) 1 2,„(х)У~ (са) ахз (л+ 1)! Их (л+ 2)! " (л+ 3)! Если х принимает одно из значений хе, х,, ..., хан то последка[и член справа обратится в нуль и остаточный член упростится.

ааналогичные рассуждения можно провести и для любого й (л. В общем случае получим: уЧЮ(х) = й[[1'(хр; х,; ...; хь)+(ае-[-аа-+... -+а )у(хе; х,; ...; ха+1)+ -[-(арпа+враз-+ ... -[-аааа,„,) )С' ХП .;,;...; „Д+ ... -[-(;,.... а+... ° ° ° +аа+ааааа ° ° а~-1)У(хе: ха: ° ° ° ' х~)[ -[- ~а +-„— „[аз„(х))'(х; хе; х,;...; х„)[. (12) Для упрощения остаточных членов нам понадобятся выражения И за — у(х; х„х„...; х„) (т (п).

Покажем, что — /(х; хз; х,; ...; х„) = т! Г" (х; х; ...; х; х,; х,;...; х„). (13) тР1 раз Как зто следует из предыдущего, при ла = 1 и 2 зта формула справедлива. Предположим, что она справедлива при ла=г, и докажем ее справедливость при ла = г[-1, В силу нашего предположения Иа „вЂ” „г" (х; х; ...; х„) =г[Г(х; х(...; х; хе; х,; ...; х„) гз! раз — „—, г(х; х„; х,;...; х„) = г! — у (х; х;...; х; хр; х,;...; х„). а+1 раа 223 $ 2! ФОРМУЛЪ| ЧИСЛЕННОГО ДИФФЕРЕНЦИРОВАНИЯ Воспользуемся опять определением производной ар"+'У(Х; Х,; Х11 ..., 'Х„) Лх"" у(х', х'...,; х'! х,; х,;...; х„) — у(х; х; ...; х; х„; х,;...; х„) Н-1 раз г+! Раз = Г! !Нп Ф' ФЖ х' — х Выражение в числителе последней дроби можно записать в виде гт1 ~~.', (1(х'! х';...; х'; х; х;...; х; х;...; х„)— » раз ае1-» раз — й*' " " " . *: ' : .' " ' " : *.1) = +1 = (х' — х)~ у(х'; х'! ..,; х'; х; х;...; х; х; х,;...; х„).

» раз «+2-» раз Таким образом, и формула (!3) доказана. В силу доказанной формулы остаточный член при численном отыскании производной порядка л может быть представлен в виде я = ~)„С» — „«,у(х; хе, ...,. х„) „„,", 2-0 у (.; ...; '; ..; „ ..

' ..) »Ф1 раа (14) или » Й = ~~)~~ У<" з."1! (с ) в!"-1! (х), (15) »! (» — 1)! (л+ Г+ 1)! 1 л а 0 Где 21 — некоторые точки, заключенные в интервале между наиболь- шим и наименьшим из чисел х, хе, х,, ..., х„. Еа " /(Х; Хз, Ха,...', Хл) лх'+1 «+! = Г1 !1ш р, У(х'!...; х'! х;...; х; хр; х,; ...; х„) = » раз з-~-2-» раз =(У+1)! у(х; х;...; х; хе; х,;...; х„). г+2 раз '224 числвннов дияявгвнциговлнив и интвггиговлнив (гл.

3 Если точка х находится вне отрезка, солержащего точки ха, х, , х„, то остаточный член может быть представлен более простым выражением. Для этого рассмотрим многочлен <,»(х) = Е„(х) +Се„(х) (С = сопл<). Он совпадает с функцией 7"(х) в точках ха, хн ..,. х . Подберем постоянную С так, чтобы в точке х', для которой производится оценка, имело место равенство <)<"» (х') = 7.я<" (х') +- Со4" ,(х') = у<"» (х'». Это возможно, так как все корни уравнения и<а»(х) = О лежат в наименьшем отрезке, солержащем ха, х„..., х„. Рассмотрим вспомогательную функцию т (х) = 7 (х) — Ь„(х) — См„(х). Эта функция обращается в нуль в точках ха, х„..., х„. Следовательно, первая производная ее обращается на наименьшем отрезке, солержащем точки ха, х,...., х„, в нуль по крайней мере и раз.

Проводя те же рассужления лальше. получим, что производная порядка я обратится на этом отрезке в нуль по крайней мере и + 1 — А раз. Но в силу выбора С она обратится в нуль и в точке х', лежащей вне этого отрезка. Таким образом, она обращается в нуль по крайней мере в и +2 — х точках. Снова будем последовательно применять теорему Рояля.

В конце концов, придем к выводу, что производная порядка л + ! обращается в нуль по крайней мере в одной точке Е Но фя+'»(х) = у<" г'» (х) — См<"э'» (х) = у<"+'»(х) — С(л-(- !)< Отсюда и <яэ!) ! (х') — 7.„» (х') = ( ) м<л» (х'). (л+ 1)< я (! 6) Получили более простое выражение остаточного члена. Рассмотрим пример на применение формул численного лифференцирования.

П р и м е р. По таблице 16' 20' 0,275637 0,342020 225 ФОРМУЛЫ ЧИСЛЕННОГО ДИФФЕРЕНЦИРОВАНИЯ используя формулы численного дифференцирования, найти соя !5' н сйп[5'. Составляем таблицу разделенных разностей: у(хн х,е,) У(Х4 Х44Ь Х442) У(Х4 Х444 Х442, Х442) х, 7(х) 10' 0,173648 144 0 241922 16' 0,275837 17068,50 16857,50 16595,75 — 35,17 — 43,62 20' 0,342020 Отсюда получаем,) учитывая, что в нашем случае ае = 5, а, = 1, П2— сов[5'= [7 (ха; х4)+(аа+аг) 7'(хе; х,; х,)+ 180 +(пап, +свая+агав) 7 (ха; х,; хя, ха)[ — = = [0,0!70685 — 0,0002!1 +0,00000084[ ° 57,295779 = 0,965912.

180 Множитель — справа появился за счет' того, что у нас х взято в градусном измерении. Точное значение с шестью верными знаками соз 15'= 0,965926. Используя формулу для второй производной, получим: /180)2 з)п 15' = — 2 [) (ха; х,; хв) +(ае+а, + аД ) (хе; х,; ха; х,) [1 — ~ = = 2[0.00003517 +5 0,00000084[3282,8063 = 0,257027. Прн этом Фа (х) = (х — 10') (х — 14') (х — ! 6') (х — 20~ [ — 1, )130) ' Фа (х) = [ (х — 14') (х — 16') (х — 20') + + (х — 10') (х — ! 6') (х — 20') + (х — 10') (х — 14') (х — 20') + +(х — 10')(х — 14')(х — 16') [ ! — ' [ . 1180г Точное значение з!и !5' с шестью знаками равно 0,258819. Расхождения получились довольно значительными.

Это и естественно, так как функции могут быть и очень близки друг к другу, но иметь сильно различающиеся производные. Произведем оценку погрешности. В первом случае остаточный член будет иметь следующий вид: Ум) (80 4 У(ь) (4 ) 41 2 51 22б численное диеевгвнциговлнив и интвггиговлнив [гл, 3 При х = 15' получим: вз(!5 ) = 25.0 000000092, в,'(!5 ) = О.

Таким образом, [ Я[ С вЂ” 0,000000092 — 0,0000000!9 . !(= ') в" (х)+2 ') в'(х)+2 (') в (х) 4! вз х 3! з 61 При этом вз(15)=(!80) [ 52[ — — 0 0003 52= — 0 00!56. Таким образом, [!э[ < '66136 + '26 ОООО(ЯЮ692 = О,ооообб. 24 720 И в этом случае вычислительная погрешность очень велика. 2. Формулы численного дифференцирования для равноотстоязцих узлов. Если узлы интерполирования расположены через равные промежутки, то удобнее использовать соответствуюшие интерполяционные формулы.

Так, например, взяв интерполяционную формулу Ньютона для интерполирования вперед у(.) = у (.. +гй) =-.у.+г.у', +'",, "у',-+ '" ! (! — !) (! — 2) (! — 3) 4! з в результате последовательного дифференцирования получим: ЛУЛ! 1 [ у~ +2! — ! +3!я — 6!+2 .з + ШЙх л~ 1 2! з 3! з 4Р— 18!я+ 2И вЂ” б 41 ~в+ (13) ( х ) 1 .а+Я вЂ” б .з+ 12!я — 361+22 + 3! 4! з У',+' 4!" У;+.. 1. ув(~)= р Эта величина значительно меньше фактически полученной погрешности.

В данном случае вычислительная погрешность значительно перекрывает погрешность метода. Во втором случае 227 9 2] ФОРМУЛЫ ЧИСЛЗННОГО ДИФФВРВНЦИРОВАНИЯ В частности, при х = х будем иметь: 1 «1 з 1 (! 9) (20) Здесь предполагается, что формальное разложение !и (! +Ь) = ая д« 2 =Ь вЂ” — + — — ..., доведенное до постоянных разностей, формально возводится в степень как многочлен. Дадим операторный л' вывод этой формулы. Если оператор — обозначить буквой А), то лх формула Тейлора 1(хе+И) — Г(хд+Ч (хэ)+ Г (хе)+ может быть записана так: 1(хо+И) =~!+ И() + — 0'+- ...) Г(хе) (1+д)Г(хо) = е" .Г(хо).

«и или Отсюда Беря логарифмы от обеих частей равенства, получим: ()= 1 1п(!+Ь) или О 1 (!и(!+Ь)) Получили как раз то выражение, которое было дано выше. Проверим наши формулы на примере многочлена, для которого они должны давать точные значения производных. 11ри мер. Найти методом численного дифференцирования производные первых трех порядков для многочлена х' — 2х — 5 в точке х = !. Если использовать значок Ь для разностей, то последние формулы будут иметь следующий, легко запоминающийся операторно-символический вид: 228 численное диФФВРенциРОВАние и интеГРВРОВАние )гл.

3 Составляем таблицу разностеи: уз 17 1б 18 35 51 24 59 5 По нашим формулам получаем: х М7 — — + —, )з=1, У(1)=б 6+2 — 1 ду' азу Азу Дх 2 3 Последовательные производные будут иметь вид йЧа(~) =7„+11,+'"',, 'У:+" . 7з'у" (х) = 70+ Г7ч+ ... (22) В частности, при х= ха О) .70 31 А) + 51 з (21)з )ззгаа(х ) агз 31(1 + 2 ) гз+ 0 0 5! 0 (23) Если использовать другие формулы интерполирования, то можно получить другие формулы численного дифференцирования. Возьмем, например, формулу Стирлинга ~(х) =У(хо+7з() =.70+~У,'+ 21 У,'+ еогмклы численного диеевввнцивовлния Если взять формулу Бесселя «(!в У(~)=У, +(Š— —,') У', +'" — ') У', + 3 а а Е (Еа — 1) (Š— 2) 4! (24) то получится: ЗЕз — ЗЕ +— 1 2 2Š— 1 а + —,— у, з уз + 4! (25) Еьау'"(х) =г', + 41 а а и при к=ха 1 — — У'+" 12 а У ( з) =Уь — 2 г'ь 3 а лау" (х,) = г', — — г, з а азуьл (х ) = уз ! у4 ь з Мы уже получили выражение оператора дифференцирования Е) через операторы Е1, Е!а, ЕГз, ...

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

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

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

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