Том 1 (1160083), страница 34
Текст из файла (страница 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 Далее вычисления проводим так же, как в $ б, где мы продолжали таблицу многочлена. Сначала заполняем столбец третьих разностей, затем вторых, первых и, наконец, столбец значений функции. В узловых точках записываем данные нам значения.