Интерполяционная формула Лагранжа

Лабораторная работа №2

Приближение функций, заданных таблично

Цель лабораторной работы

Приобретение практических навыков в построении интерполяционного многочлена Лагранжа и использование его для вычисления приближенных значений функций вне узлов интерполяции.

Приобретение практических навыков построения аппроксимирующих функций (аналитических зависимостей) по совокупности дискретных экспериментальных данных об изменениях значений функции при изменениях значений аргумента.

Задание

Для функций, заданных в таблице 1:

· построить интерполяционный многочлен Лагранжа и вычислить по нему значения функции для заданных значений аргумента;

· изучить технологию расчетов интерполяционных многочленов в Excel;

Для функций, заданных в таблице 2:

· Вычислить коэффициенты аппроксимирующих многочленов 1-й и 2-й степени, записать многочлены и построить их графики, на которые нанести также заданные табличные точки (расчеты выполнить в виде таблиц);

· изучить технологию регрессионного анализа с помощью Excel.

Интерполирование функций

Постановка задачи

Для функций, заданных таблицами их значений на конечном интервале, возникает необходимость вычисления значений функций для значений аргументов, отсутствующих в таблице. Тогда строят функцию, которая в заданных точках принимает заданные значения, а в остальных точках интервала приближенно представляет заданную функцию. А затем вычисления значений функции для любых значений аргумента в области определения заданной таблично функции выполняют по построенной функции. Задача интерполирования – построение такой приближенной функции. Чаще всего интерполирующую функцию отыскивают в виде алгебраического многочлена. Геометрически задача интерполирования заключается в построении кривой , проходящей через заданную таблично систему точек.

Интерполяционная формула Лагранжа

Пусть функция в точках соответственно принимает значения .Требуется построить многочлен степени не выше n, принимающий в точках (узлах интерполирования) значения . Расстояние между узлами интерполирования может быть различным. Решение этой задачи – многочлен Лагранжа.

Интерполяционная формула Лагранжа в общем виде

где – базисные функции, числитель которых содержит все разности , а знаменатель – все разности за исключением .

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

Пример 1.1. Запишем интерполяционный многочлен для функции, заданной тремя точками:

Построить интерполяционный многочлен Лагранжа и вычислить приближенно значение функции для .

Для построения интерполяционного многочлена воспользуемся формулой при :

Проверяем значения функции для узлов интерполяции:

Вычисляем

Интерполяционную формулу в Ехсеl можно построить достаточно простым способом. С практической точки зрения главная проблема заключается в вычислении в произвольной точке значений базисных функций.

Пример 1.2. Запишем интерполяционный многочлен для функции, заданной шестью точками:

На рис. 1.1 представлены исходные данные, по которым будет выполняться интерполяция. На этом же рисунке проиллюстрирован процесс определения первой базисной функции.

Рис. 1.1

В диапазоне ячеек А6:А11представлены заданные значения аргументов функции, а в диапазоне ячеек В6:В11– значения функции для узловых точек аргумента. В ячейку А2вводится значение аргумента, для которого необходимо вычислить значение интерполяционного полинома. Значение полинома будет выводиться в ячейке В2.Важным моментом является заполнение ячеек в диапазоне С6:С11,где будут отображаться значения базисных функций в точке, указанной в ячейке А2.Именно по этим значениям и значениям ячеек из диапазона А6:А11определяется значение интерполяционного полинома (ячейка В2).

Диапазон С6:С11заполняется так: отдельно первая и последняя ячейки диапазона, а все остальные ячейки — распространением одной формулы. В частности, в ячейку С6вводится формула

=ПРОИЗВЕД(А2-А7:А11)/ ПРОИЗВЕД(А6-А7:А11),

согласно которой определяется первая базисная функция. Сразу следует отметить, что и эта формула, и все прочие формулы из диапазона С6:С11,вводятся как формулы для диапазонов, т.е. с помощью нажатия комбинации клавиш Ctrl+Shift+Enter.Причина состоит в том, что аргументами функции ПРОИЗВЕД() указываются результаты арифметическихопераций с диапазонами.

Далее заполняются ячейки из диапазона С7:С10.Для этого в ячейку С7вводится формула

=ПРОИЗВЕД($A$2-$A$6:A6;$A$2-A8:$A$11)/ПРОИЗВЕД(A7-$A$6:A6;A7-A8:$A$11).

Абсолютные и относительные ссылки в формуле подобраны так, чтобы при ее копировании в следующие ячейки, ссылки на начальную А6 и конечную А11ячейки диапазона, равно как и на ячейку А2со значением переменной, для которой вычисляется базисная функция (и весь полином), оставались неизменными. Это абсолютные ссылки. Вместе с тем, в процессе копирования формулы произведения в ней вычисляются без учета значения аргумента в той строке, где размещена формула. После ввода формулы в ячейку С7данная формула с помощью маркера заполнения копируется во все ячейки, вплоть до С10 (Рис. 1.2).

Рис. 1.2

Наконец, в ячейку С11необходимо ввести формулу

=ПРОИЗВЕД(А2-А6:А10)/ ПРОИЗВЕД(А11-А6:А10).

Рис. 1.3

Поскольку формулы из начальной С6 и конечной С11ячеек диапазона С6:С11никуда копировать не предполагается, то и ссылки там относительные. Результат можно видеть на рис. 1.3.

После этого осталось только вычислить интерполяционный полином. Для этого достаточно в ячейку В2ввести формулу

=СУММПРОИЗВ(В6:В11;С6:С11)

Эта формула вводится как обычная, то есть нужно нажать клавишу Enter.Результат представлен на рис. 1.4.

Рис.1.4

Чтобы посчитать значение полинома в какой-то точке, соответствующее значение следует ввести в ячейку А2. В ячейке В2практически сразу появится результат.

Ответ:

На рис. 1.5 проиллюстрирована ситуация, когда в качестве аргумента указано узловое значение. Как и следовало ожидать, в узловой точке значение интерполяционного полинома равно табличному значению функции в этой точке, а все базисные функции, кроме той, что соответствует указанному узлу, равны нулю. Отличная от нуля базисная функция равна единице.

Рис. 1.5