Страница
3
Численное интегрирование. Аппарат интерполирования функции лежит в основе построения многих квадратурных и кубатурных формул. Такого рода формулы строятся путем замены интегрируемой функции на всей области или на её составных частях интерполяционными многочленами того или иного вида и последующим интегрированием этих многочленов. Например квадратурные формулы наивысшей алгебраической степени точности, так называемые квадратурные формулы Гаусса:
где - знакопостоянная весовая функция, получаемая в результате замены функции
интерполяционным алгебраическим многочленом, построенным по корням
ортогонального относительно веса
многочлена степени n.
Изложенная выше схема построения формул для приближенного вычисления интегралов применима и в многомерном случае
Формулы численного дифференцирования, в основе которых лежит интерполирование, получаются в результате дифференцирования интерполяционных многочленов. Ввиду неустойчивости задачи численнго дифференцирования относительно ошибок использования значений функций в узлах шаг интерполирования должен согласоваться с погрешносьтью значений функций. Поэтому на практике нередки случаи, когда известная на густой сетке функция используется в данной задаче не во всех точках, а на более редкой сетке.
При численном решении интегральных уравнений, известная функция заменяется в интегральном уравнении каким-либо интерполяционным приближением (интерполяционным алгебраическим многочленом, интерполяционным сплайном и т.д.) с узлами интерполирования
, а приближенные значения
для
находятся из системы, полученной после подстановке вместо независимости переменной x узлов интерполирования
. В случае нелинейных интегральных уравнений приближенные значения
находятся соответственно из нелинейной системы.
Интерполяционная формула- для приближенного вычисления значений функции , основанного вычисления на замене приближаемой функции
более простой в каком- то смысле функцией
|
|
наперед заданного класса, причем параметры выбираются так чтобы значения
совпадали с известными заранее значениями
для данного множества
попаро различных значений аргумента:
такой способ приближенного представления функций называется интерполированием, а точки , для которых должны выполняться условия
, - узлами интерполяции.
В ряде случаев (например, при интерполировании алгебраическими многочленами) параметры могут быть явно выражены из системы
, и тогда
непосредственно используется для приближенного вычисления значений функции
.
Интерполяционный процесс- процесс получения последовательности интерполирующих функций при неограниченном возрастании числа n узлов интерполирования. Если интерполирующие функции
представлены в виде частных сумм некоторого функционального ряда, то последний иногда называется интерполяционным рядом.Целью построения интерполяционного полинома чаще всего является, по крайней мере в простейших первоначальных задачах интерполирования, приближение в каком- то смысле по средствам интерполирующих функций
, о которой или имеется неполная информация, или форма которой слишком сложна для непосредственного использования.
Интерполяционная формула Эверетта:
Интерполяционные формулы Грегори- Ньютона построенные по нисходящим или восходящим разностям, наиболее целесообразно применять в начале или конце таблицы. При этом для достижения высокой степени точности иногда приходится рассматривать разности, отстоящие достаточно далеко от интересующих нас значений функции или
. Поэтому на средних участках таблицы лучше результаты дают интерполяционные формулы, построенные на базе центральных разностей, то есть разностей, которые ближе всего расположены к центральной сотке, содержащей
.
К интерполяционным формулам с центральными разностями относятся формулы Гаусса, Стирлинга, Бесселя, Эверетта и многие другие; формула Эверетта получила наибольшее распространение, она была получена 1900 г.:
где ;
;
.
Формуле Эверетта так же можно придать форму, наиболее удобную для вычисления:
если для ее коэффициентов ввести обозначения