Аппроксимация функции

Содержание

Слайд 2

Постановка задачи Простейшая задача: в дискретные моменты времени наблюдаются значения функции

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

Простейшая задача: в дискретные моменты времени наблюдаются значения функции ;

требуется восстановить её значения при других . .

26.09.2014

Слайд 3

2. Для функции известно её аналитическое представление, но вычисление каждого значения

2. Для функции известно её аналитическое представление, но вычисление каждого значения

сопряжено с большим объёмом вычислений. Например,

26.09.2014

Слайд 4

3. Функция задаётся своими значениями в узлах , из интервала .

3. Функция задаётся своими значениями в узлах , из интервала .


В вычислительном процессе используется эта таблица.

26.09.2014

Слайд 5

4. Задача численного решения определённого интеграла или дифференциального уравнения. 26.09.2014

4. Задача численного решения определённого интеграла или дифференциального уравнения.

26.09.2014

Слайд 6

Области использования аппроксимации : моделирование; планирование и статистическая обработка данных; определение

Области использования аппроксимации :
моделирование;
планирование и статистическая обработка данных;
определение значений функции

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

26.09.2014

Слайд 7

Аппроксимация Точная в узлах — интерполяция Приближённая в узлах Кусочно —


Аппроксимация

Точная в узлах — интерполяция

Приближённая в узлах

Кусочно —
линейная

Многочленами

Нелинейными

функциями

Гаусса

Ньютона

Лагранжа

Линейными функциями

Нелинейными функциями

26.09.2014

Экстраполяция

Слайд 8

Аппроксимация 26.09.2014

Аппроксимация

26.09.2014

Слайд 9

Если аппроксимация функции происходит в промежуточных узлах, т.е. , причём тогда

Если аппроксимация функции происходит в промежуточных узлах, т.е. , причём тогда

говорят о точной интерполяции.
Интер – между.

26.09.2014

Слайд 10

Интерполирование 26.09.2014

Интерполирование

26.09.2014

Слайд 11

Если аппроксимация функции происходит вне рассматриваемого отрезка , тогда говорят об экстраполяции. экстра – вне. 26.09.2014

Если аппроксимация функции происходит вне рассматриваемого отрезка , тогда говорят об

экстраполяции.
экстра – вне.

26.09.2014

Слайд 12

ИНТЕРПОЛИРОВАНИЕ Основные понятия. Задано множество точек , принадлежащих отрезку и значения

ИНТЕРПОЛИРОВАНИЕ

Основные понятия.
Задано множество точек , принадлежащих отрезку и значения функции
в

этих точках.
Требуется найти функцию , определённую на этом отрезке такую, что в указанных точках

26.09.2014

Слайд 13

ОБОБЩЁННЫЙ ПОЛИНОМ – числовые коэффициенты, – множество функций, определённых на рассматриваемом

ОБОБЩЁННЫЙ ПОЛИНОМ
– числовые коэффициенты,
– множество функций, определённых на рассматриваемом отрезке

и линейно независимых на нём.
Такие функции называются БАЗИСНЫМИ.

26.09.2014

Слайд 14

1. Последовательность степеней . 2. Последовательность тригонометрических функций . 3. Последовательность

1. Последовательность степеней .
2. Последовательность тригонометрических функций .
3. Последовательность экспонент

.
4. Последовательность дробно-рациональных функций .

26.09.2014

Слайд 15

Многочлен Тейлора m-ой степени: Для можно использовать известные разложения функций. 26.09.2014

Многочлен Тейлора m-ой степени:
Для можно использовать известные разложения функций.

26.09.2014

Слайд 16

Погрешность метода — остаточный член формулы Тейлора в форме Лагранжа: где

Погрешность метода — остаточный член формулы Тейлора в форме Лагранжа:
где c

лежит между x и x0.
Если c не принадлежит интервалу сходимости ряда Тейлора, то погрешность не уменьшается.

26.09.2014

Слайд 17

Существование и единственность интерполяционного многочлена 26.09.2014 Теорема. Для таблично заданной функции

Существование и единственность интерполяционного многочлена

26.09.2014

Теорема.
Для таблично заданной функции на множестве узлов

существует единственный интерполяционный многочлен степени , удовлетворяющий условию , .
Доказательство:
Слайд 18

Будем искать многочлен Нужно найти коэффициенты a0, a1, ..., an. Они

Будем искать многочлен
Нужно найти коэффициенты a0, a1, ..., an.
Они должны удовлетворять

СЛУ
В системе (m+1) неизвестных, (n+1) уравнение.

26.09.2014

Слайд 19

Если система крамеровская, то решение СУЩЕСТВУЕТ и ЕДИНСТВЕННОЕ. Пусть m+1 =

Если система крамеровская, то решение
СУЩЕСТВУЕТ и ЕДИНСТВЕННОЕ.
Пусть m+1 =

n+1, т.е. m = n. Главный определитель системы
— определитель Вандермонда.
Если все узлы различны, то .
Теорема доказана.

26.09.2014

Слайд 20

Метод нахождения коэффициентов, используемый при доказательстве теоремы, называется методом неопределённых коэффициентов. 26.09.2014

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

26.09.2014

Слайд 21

Если m Например, m = 1, n = 2, найти линейную

Если m < n, СЛУ может быть несовместна.
Например, m = 1,

n = 2,
найти линейную
функцию (прямую),
проходящую
через три точки:

26.09.2014

Слайд 22

Если m > n, СЛУ имеет бесконечно много решений. Например, m

Если m > n, СЛУ имеет бесконечно много решений.
Например, m

= 2, n = 1,
найти квадратичную
функцию (параболу),
проходящую
через две точки:

26.09.2014

Слайд 23

Интерполяционный многочлен Лагранжа Выразим многочлен Ln(x) как линейную комбинацию значений f0, f1, ..., fn: 26.09.2014

Интерполяционный многочлен Лагранжа

Выразим многочлен Ln(x) как линейную комбинацию значений f0, f1,

..., fn:

26.09.2014

Слайд 24

Рассмотрим в качестве подсказки частные случаи. 1) n = 1: x0,

Рассмотрим в качестве подсказки частные случаи.
1) n = 1: x0, x1

— узлы, f0, f1 — значения в узлах
Найти

при x0

при x1


и

26.09.2014

Слайд 25

2) n=2: x0, x1, x2 — узлы, f0, f1 , f2—

2) n=2: x0, x1, x2 — узлы, f0, f1 , f2—

значения в узлах
Найти

при x0

при x1

при x3

26.09.2014

Слайд 26

Опр. Интерполяционным многочленом Лагранжа называется полином . Опр. Лагранжевы коэффициенты —

Опр. Интерполяционным многочленом Лагранжа называется полином

.
Опр. Лагранжевы коэффициенты —

Замечание:
Лагранжевы

коэффициенты удовлетворяют тождеству

для каждого i = 0, ..., n.

, т.к. обладают свойством:

26.09.2014

Слайд 27

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

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

Так как задано

4 узла интерполяции, то степень полинома не выше третьей.

26.09.2014

Слайд 28

26.09.2014 Можно вычислить приближённое значение в точке .

26.09.2014

Можно вычислить приближённое значение в точке .

Слайд 29

26.09.2014 Погрешность интерполяции

26.09.2014

Погрешность интерполяции

Слайд 30

Погрешность интерполяционного многочлена Лагранжа Предполагаем, что . Ln(x) — многочлен Лагранжа:

Погрешность интерполяционного многочлена Лагранжа
Предполагаем, что .
Ln(x) — многочлен Лагранжа: для

всех i=0,...,n
[a,b] — отрезок, содержащий узлы x0, x1, ..., xn.
Найдем оценку отличия значения f (x) от значения Ln(x) в точке , не совпадающей ни с одним из узлов, иначе, величину остаточного члена

26.09.2014

Слайд 31

Запишем равенство где – многочлен определённый через узлы x0, x1, ...,

Запишем равенство
где – многочлен определённый через узлы x0, x1, ..., xn


С – некоторая постоянная (параметр).
Подберём параметр С так, чтобы обращалась в нуль в точке , для которой делаем оценку, т.е. и
.
К функции на каждом из отрезков применима теорема Ролля

26.09.2014

Слайд 32

Введём в рассмотрение функцию 26.09.2014

Введём в рассмотрение функцию

26.09.2014

Слайд 33

Итак, существует : . Тогда и , т. к. получаем Отсюда 26.09.2014

Итак, существует : .
Тогда
и , т. к.
получаем
Отсюда

26.09.2014

Слайд 34

Для остаточного члена получаем выражение 26.09.2014

Для остаточного члена получаем выражение

26.09.2014

Слайд 35

Точное представление f (x) через её интерполяционный многочлен Лагранжа : где

Точное представление f (x) через её интерполяционный многочлен Лагранжа :
где и

зависит от x.
Можно оценить предельную абсолютную погрешность интерполирования на отрезке с помощью формулы

26.09.2014