- •Лекция №1
- •§ 1. Задача вычисления.
- •§ 2. Абсолютная и относительная погрешности
- •§ 3. Неустранимая погрешность значения функций для приближенных значений аргументов. Погрешности арифметических операций.
- •Лекция №2 Численные методы линейной алгебры
- •Формальное решение. Устойчивость.
- •Обусловленность матрицы. Погрешности.
- •Лекция №3
- •1. Схема единственного деления
- •3. Расчетные формулы
- •Лекция № 4 Метод Гаусса с выбором главного элемента (оптимальный метод).
- •Применения метода Гаусса к вычислению определителей и обратных матриц.
- •Лекция № 5 Итерационные методы решения систем линейных алгебраических уравнений.
- •Лекция № 6 Метод Зейделя (модификация метода итераций).
- •Тогда условие окончания итерационного процесса Зейделя будет иметь вид:
- •Лекция № 7 Методы решения нелинейных уравнений и систем нелинейных уравнений.
- •1. Метод деления отрезка пополам (метод бисекций или дихотомия).
- •Метод хорд (метод линейной интерполяции).
- •3. Метод Ньютона (метод касательных или метод линеаризации).
- •4. Метод итераций (задача о неподвижной точке).
- •Оценка погрешности приближений:
- •Лекция № 8
- •1. Метод итераций для системы двух уравнений.
- •2. Метод Ньютона для системы двух уравнений.
- •Лекция №9 Алгебраическая проблема собственных значений.
- •Лекция № 10 Приближение функций и их производных.
- •Постановка задачи приближения функций.
- •2. Оценка погрешности полиномиальной интерполяции.
- •Лекция № 11 Интерполяционный многочлен Ньютона с конечными разностями.
- •Лекция № 12 Метод наименьших квадратов и наилучшие среднеквадратические приближения.
- •О нормальной системе мнк при полиномиальной аппроксимации.
- •Лекция №13 Сплайн интерполяция
- •Лекция № 15
- •Метод Эйлера – разные подходы к построению.
- •Методы Рунге – Кутта.
- •Лекция № 16
- •Лекция № 17 Разностные схемы для уравнений параболического типа.
- •Лекция №18
- •Лекция № 19 Разностные схемы для уравнений эллиптического типа.
Лекция № 10 Приближение функций и их производных.
Постановка задачи приближения функций.
Простейшая задача приближения функций заключается в следующем.
Пусть известны значения некоторой функции f(x) при заданных x0 <x1 < <xn на некотором отрезке [x0 , xn]. Требуется получить значения f(x) для такого значения аргумента x, которое не входит в отрезок [x0 , xn ]., нo и не совпадает ни с одним из значений xi, i=0, ,n
x |
x0 |
x1 |
|
xn |
f(x) |
f0 |
f1 |
|
fn |
При этом основная цель интерполяции получить быстрый и экономичный алгоритм вычисления значений функции F(x), для значений x не содержащихся в исходной таблице, т.е. xÎ[a, b] и x¹xi.
f(x) @ F(x, a0, a1,…, an)
Параметры a0, a1,…, an - определяются из условий совпадения f(x) и приближающей функции в точках x0 , x1 , , xn - узлах интерполяции
F(xi, a0, a1,…, an) = f(xi), i = 0, , n
Такой способ называется интерполированием.
Далее будем рассматривать задачу интерполирования многочленами, но это не единственный способ. Иногда удобнее приближать функцию тригонометрическими функциями или ln(f(x)). Интерполяционный многочлен Лагранжа.
На практике применяют аналитический способ нахождения коэффициентов aj , применяя другой способ записи многочлена.
Определим - символ Кронекера.
Задача интерполирования будет решена, если мы построим такие многочлены Fi(x) степени не выше n, такие, что .
Тогда многочлен будет искомым интерполяционным многочленом. Действительно,
Кроме того, Fn - многочлен степени n .
Поскольку обращается в ноль в точкахx0 , x1 , ,xn в n точках, то Fi(x) делится на (x-xj)
Из условия находим, что
и тогда ,
тогда искомый многочлен имеет вид
Называется интерполяционным многочленом Лагранжа.
Существуют другие формы записи интерполяционного многочлена, Формула Ньютона и ее варианты. Этим вариантам записи соответствует различие в величине вычислительной погрешности, а также различное количество арифметических операций.
Покажем единственность полинома Лагранжа.
Предположим обратное, пусть - полином степени не вышеn и
Тогда полином - обращается в нуль в(n+1) точках x0 , x1 , ,xn и степени не выше n, т.е.
По следствию из основной теоремы алгебры многочленов – многочлен n –ой степени не может иметь более n корней.
Формуле Лагранжа можно придать более сжатый вид:
Дифференцируя по х это произведение, получим:
При x=xi
Отсюда, получаем
Замечание. Нетрудно оценить число арифметических действий, в главном порядке по n это есть величина O(n2).
Рассмотрим 2 частных случая интерполяционного полинома Лагранжа.
При n =1 имеем 2 узла и формула Ln представляет собой уравнение прямой, проходящей через 2 заданные точки:
При n =2 получим уравнение параболы L2(x), проходящей через 3 точки:
L1(x) и L2(x), называются соответственно формулами линейной и квадратичной интерполяции.