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

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

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

(з — ав). Подынтегральная функция аналитична на С и внутри С, за исключением точек хв, вн..., з„, Следовательно, интеграл будет равен сумме вычетов относительно каждой из этих точек Но ч (1) — и (х) ю (х) йщ . ("„„... /(~)(С вЂ”.,) =/(х,) „, „„,,„,-. Отсюда Х ~ ( в) ч' (лл) (в — хл) ь-о т. е. Р(х) является интерполяционным многочленом Лагранжа. Далее, представим Р(х) в виде разности Р (х) = / — с(ч — -з — /— 1 1 У(С) ч(З) Г У(~) 2яГ ./ яг ./ ч (С) (à — х) Ж(. С с 196 твовия интввполивовлния и нвкотовыв вв пвиложвния (гл.

2 Первый член в силу интегральной формулы Коши равен /(г). Следовательно, ч (л) (' у(С) У(з) =Р(Я)+ 2вГ,/ -(()(à — ) с Итак, остаточный член интерполяционной формулы Лагранжа в нашем случае может быть представлен в виде ч (л) )' У(С) д 2вг / ч (с) (1 — л) о ф 14. Применение интерполирования для составлении таблиц Теория интерполирования имеет большие приложения при составлении таблиц функций. Получив зздание на составление таблиц тех или иных функций, математик должен решить перед началом вычислений ряд вопросов. Должна быть выбрана формула, по которой будут производиться вычисления.

Эта формула может изменятся от участка к участку. Обычно формулы для вычисления значений функции, использующие способ задания функции, бывают громоздкими и поэтому их используют для получения некоторых опорных значений и затем путем субтабулирования сгущают таблицу. Формула, дающая опорные значения функции, должна обеспечивать нужную точность таблиц с учетом последующего субтабулирования. Если предполагается составить таблицы с постоянным шагом, то должен быть определен шаг таблицы. Шаг таблицы связан с двумя факторами: объемом таблиц и интерполяционной формулой, по которой будут вычисляться промежуточные значения уже в готовой таблице. Чем больше будет шаг, тем больше членов интерполяционной формулы придется использовать при пользовании этой таблицей на практике.

Это создает некоторые неудобства при использовании таблицы. С другой стороны, чем меньше шаг, тем больше объем таблиц, что также не очень удобно. Математик должен как-то согласовать действие этих противоположных факторов с учетом средств вычислений, имеющихся в распоряжении потребителя. Если таблица должна быть введена в быстродействующую машину, то особенно важно уменьшить ее объем, При этом можно отказаться от постоянства шага и использовать, например, узлы Чебышева на отдельных участках, для которых, как мы видели, получается наилучшая оценка остаточного члена интерполяционной формулы. При определении шага таблицы будут иметь значение и такие факторы, как наличие вычислительных средств и время, отведенное на вычисления.

Мы не можем здесь входить в детали каждого из поставленных вопросов и остановимся лишь на выборе шага и субтабулировании, ф 14] пвимкнкник интквполигования для состлвлкния тавлиц 197 Чаще всего таблицы функций составляются так, чтобы была допустима линейная интерполяция (т. е.

интерполяция с использованием первых двух членов формулы). В этом случае остаточный член будет иметь вид гс (х) = ~ () дзг(г — 1). Здесь $ принадлежит интервалу между двумя соседними табличными значениями аргумента, в котором лежит х. а 1 заключено между 0 и 1, Произведение 1(1 — 1) принимает наибольшее по модулю 1 1 значение при г= —. Это значение равно —. Следовательно, 4 8 Чтобы ошибка интерполирования не превышала по абсолютной величине а, необходимо выбирать Ь, удовлетворяющим условию Нужно помнить, что наряду с втой ошибкой — ошибкой метода, при практическом вычислении промежуточных значений будут воз- никать еще неустранимая погрешность и погрешность округлений.

Как мы видели ранее, неустранимая погрешность при линейной интерполяции будет равна погрешности табулированных значений функций. Погрешность округления будет зависеть от вычислитель- ных средств и от программы вычислений. Поэтому здесь мы ее касаться не будем. Совершенно аналогично можно исследовать квадратичную интер- поляцию и интерполяцию более высоких порядков.

Если, например, используется интерполяционная формула Эверетта, то остаточный член будет иметь вид ущ (1) Яз(х)= ~ „) Д~(1' — 1)(1 — 2), и в этом случае наибольшее значение для 1г(Р— 1)(1 — 2)! на [0,1] 1 й будет достигаться при г= —. Это значение равно 1 .

Таким 2 16 ' образом, з < 12К л'М ° н для того чтобы ошибка квадратичной интерполяции не превышала а, нужно, чтобы шаг и удовлетворял условию д<у Ф. И здесь, кроме этой ошибки метода, возникают неустранимая ошибка и ошибка округления. Неустранимая ошибка будет такова же, квк 198 тяовия интявполивовлния и нвкотовыв ве пвиложвния (гл. 2 И-') . ~„=~(.+ — ",)=ИМ+ — , 'д~( 9+ ' ' дЧ( )+ — ( — — 1)( — — 2) й ( — — 1)( — — 2)( — — 3) — ( — — 1)( — — 2) ( — — 3)( — — 4) +, Д'|(х,). (З~ й Обозначим через Д, разность с шагом —. Тогда й' д г(Ха) ду(х)+ дзг(»)+ (! )(! й)дзг(ха) (1 — й) (1 — 2й) (1 — Зй) 4! й4 (1 — й) (1 — 2й) (1 — Зй)(1 — 4й) да у( 5! йз Таким образом, операторы д и д, связаны следующим соотношением: (1 — й) (1 — 2й) (! — Зй) (1 — й) (1 — 2й) (1 - - Зй) (1 — 4й) 4(йа + 5!йз (5; и для формулы Гаусса, взятой до третьих разностей. Как видно из приведенной ранее таблицы, неустранимая погрешность не может больше чем в 1,4 раза превысить абсолютные погрешности табулированных значений.

Аналогично можно исследовать и другие формулы. Рассмотрим теперь вопрос о сублгабулировамии. Как применяется формула Эверетта для субтабулирования, мы уже знаем. Приведем здесь еще один способ субтабулирования. Пусть ув, уо ... — данные последовательные значения функции, соответствующие шагу аргумента, равному й, и д, д', ... — их разности. Предположим, что нам нужно сгустить таблицу в й раз, й т, е. новый шаг будет —. Обозначим новые значения функции й ' соответственно черезУоо Уш Уев ° ° ° го,к-ь Ло Л» - ° Л,й-ь ° ° ° здесь (!в= уо а )ы — последовательные значения функции между х! и х;+и Будем для определенности считать, что разности пятого порядка исходных значений функции постоянны, Найдем выражения для разностей значений функции с новым шагом через разности с прежним шагом.

Для этого используем операторное исчисление, примененное нами ранее для вывода интерполяционных формул. По формуле Ньютона для интерполирования вперед будем иметь: ф 14! пгименение интевполиеовлния для составления теялиц 199 Следовательно, степени их будут связаны таким соотношением: дч 11д 1 — Ддз ! (1 — а)(1 — 2а) „г <1 —. а) (1 — 2Гг) (1 — Згг) (! — а) (1 — 2д) (! — 3<г) (1 — 4<г) ! и Отсюда последовательно получаем: дз ! дг ! а дз (! а)(7 11а) д» (! а) (! 2а) дь аь аз 12а» 6дз дз ! дз 3 (1 — <г) д» (1 — а) (5 — 7а) ! = аз За» ' 4аз ! 2(1 — а) д! = — — д'+ —,— д', аТ аь ь 1 д,= — ьд.

аь 9с 0,156434 0,207912 0,258819 0,3090! 7 0,358368 51478 12' — 571 — 709 — 847 50907 15' 50198 !За 49351 21а В нашем случае 7»=6. Последовательно получаем: сбп 9' = з1п 9', з!п 9'30' = згп 9 + У! Уз ! 6 гд 72 ! 1296 'Ь' з!п 1О" = э<п 9' + —,7' — — уя + — ~~, 3 15 9 ! 81 "й' После того как получены разности, нетрудно, используя постоянство разностей пятого порядка, произвести субтабулирование.

Сначала заполняем столбец пятых разностей, затем четвертых и т. д., пока не придем к значениям функции. Эти формулы, связывающие д и д,, можно получить и без операторного исчисления. Рассмотрим, как это делается на примере. Возьмем те же значения з!п х, которые мы использовали ранее в 9 7, и получим таблицу с шагом в 30'. Таблица исходных значений функции и их разностей выглядит следующим образом: 200 твовия интвеполивования и нвкотовыв ив пеиложвния !гл. 2 з!и 10' 30' = з1п 9'+ — ут — — Г" -+ — ~а 1 1 1 2 Ч» 6» 16 '/»' 2 1 4 яп 11' = а1п 9'+ — 7/ — — 7Я+ — Я . 3 Ч» 9 г 81 'Ь' Будем обозначать разности с новым шагом чертой сверху.

Первые разности будут выражаться так через предыдущие разности: 1 5 55 1 г 1 25 1 г 1 1 1 1 17 Х', = — У' -+ — 7' — — 7.' '/» 6»/» 72 г 1226 '/»' Вторые разности будут иметь вид: ! 5 — 7' — — У' 36 г 216 'ь' 1 4 — 7' — — 7",, 36 г 216 Н»' 1 3 — Р— — /-; 36 ! 216 "/' Третьи разности будут равны: 1 /*/. 216/ч»' Подставляя сюда числовые значения. получим: 0.051478+ — 0,000571 — , . 0,000138 =0,00861346, ~д —— — 6 ° 0,000571+ — ° 0,000138 = — 0,00001267, 7з . 0 000138 0 00000064 '/ 216 Далее вычисления проводим так же, как в $ б, где мы продолжали таблицу многочлена. Сначала заполняем столбец третьих разностей, затем вторых, первых и, наконец, столбец значений функции. В узловых точках записываем данные нам значения.

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

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

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

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