Вычисление значений функции с помощью интерполяционного полинома Лагранжа

Страницы работы

Содержание работы

Интерполяция нелиненых функций

Лабораторная работа № 2.8. Вычисление значений функции с помощью интерполяционного полинома Лагранжа

Цель: Изучить интерполяционный полином Лагранжа, область его применения и освоить решение прикладных задач с его помощью.

Оборудование и программное обеспечение: Персональный компьютер под управлением операционной системы семейства Windows, пакет Microsoft Office, интегрированная среда программирования Turbo C.

Введение. При решении прикладных задач иногда возникает необходимость вычисления значения таблично заданной функции f(x) при некотором значении аргумента x0, которое не совпадает ни с одним из узлов xi (i=0, 1, ¼ , n), в которых известны её значения f(xi). Поэтому возникает проблема доопределения функции f(x), используя известные её значения f(xi), для всех xÎ [x0, xn]. Один из способов доопределения функции заключается в построении интерполяционного полинома Лагранжа Ln(x), значения которого в узлах xi совпадают с  f(xi). Такой полином имеет вид

 .

(4.8.1)

Формы представления интерполяционного полинома Лагранжа. После того как построен полином Лагранжа (4.8.1), то в качестве приближенного значения функции f(x0) берется его значение в точке x0, то есть

 .

(4.8.2)

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

Таблица 4.8.1

Представление разностей полинома Лагранжа

(x-x0)

(x0-x1)

(x0-x2)

¼

(x0-xn)

(x1-x0)

(x-x1)

(x1-x2)

¼

(x1-xn)

(x2-x0)

(x2-x1)

(x-x2)

¼

(x2-xn)

¼

¼

¼

¼

¼

(xn-x0)

(xn-x1)

(xn-x2)

¼

(x-xn)

Обозначив произведение элементов i – ой строки таблицы 4.8.1 через Di(x) (i = 0,1, ¼ , n), а произведение элементов главной диагонали – через Õn+1(x)

 ,

(4.8.2)

 ,

(4.8.3)

выражение для полинома Лагранжа (4.8.1) можно переписать в виде

 .

(4.8.4)

Формулы (4.8.2) –(4.8.4), как правило, используются для вычисления значения функции f(x), которая задана таблично в неравноотстоящих узлах, при значении аргумента x0, не совпадающего ни с одним из узлов.

На практике не редки случаи, когда функция f(x) задается таблично в узлах xi, соседи которых удалены на одинаковое расстояние. В таком случае говорят, что функция задана таблично в равноотстоящих узлах xi (i=0, 1, ¼ , n). Рассмотрим, какой вид примет интерполяционный полином Лагранжа (4.8.4) при таком выборе узлов.

Показано, что если ввести следующие обозначения

 ,

(4.8.5)

 ,

(4.8.6)

то полином Лагранжа (4.8.4) можно записать в виде следующего соотношения

 .

(4.8.7)

Полученные выражения (4.8.5) – (4.8.7) следует применять для вычисления значения функции f(x) таблично заданной в равноотстоящих узлах.

Таким образом, в зависимости от того, каким образом заданы узлы xi (i=0, 1, ¼ , n), в которых определены значения функции f(xi), имеем две формы представления полинома Лагранжа, а именно (4.8.2) – (4.8.4) и (4.8.5) – (4.8.7).

Пример выполнения задания. Предположим, что требуется вычислить значение функции f(x) = y(x) при x0 = 0,263, которая определяется табличными значениями

Таблица 4.8.2

Значения функции f(x)

x

y

0,05

0,050042

0,10

0,100335

0,17

0,171657

0,25

0,255342

0,30

0,309336

0,36

0,376403

Так как узлы xi расположены не равномерно, воспользуемся формулами (4.8.2) – (4.8.4). Результаты вычислений приведены в таблице 4.8.3.

Таблица 4.8.3

Данные расчетов

i

Значения разностей

Di

f(xi)/Di

0

0,213

-0,05

-0,12

-0,20

-0,25

-0,31

-0,19809×10-4

-2526,2

1

0,05

0,163

-0,07

-0,15

-0,20

-0,26

0,44499×10-5

25547,7

2

0,12

0,07

0,093

-0,08

-0,13

-0,19

-0,154365×10-5

-111202,0

3

0,20

0,15

0,08

0,013

-0,05

-0,11

0,1716×10-6

1488007,0

4

0,25

0,20

0,13

0.05

-0,037

-0,06

0,7215×10-6

428740,0

5

0,31

0,26

0,19

0,11

0,06

-0,097

-0,980402×10-6

-38392,7

Похожие материалы

Информация о работе

Тип:
Задания на лабораторные работы
Размер файла:
472 Kb
Скачали:
0