Аппроксимация функций методом наименьших квадратов
Аппроксимация функций методом наименьших квадратов - раздел Образование, Погрешность. Определение погрешности
В Прикладных Задачах Бывает Необходимо Аппроксимировать Функц...
В прикладных задачах бывает необходимо аппроксимировать функцию единым аналитическим выражением на большом интервале изменения независимой переменной, что невозможно методами интерполяции. К тому же, если функция определена со значительной случайной погрешностью, то интерполировать ее вообще не имеет смысла. Такая ситуация возникает, например, при аппроксимации результатов экспериментальных исследований механических свойств материалов.
В гильбертовом пространстве L2(r) функций, интегрируемых с квадратом с весом r(x)>0 на [a, b] норму вводим следующим образом:
(34)
где скалярное произведение определено так:
(35)
Для таблично заданной функции скалярное произведение имеет вид
(36)
Тогда условие (33) принимает вид условия наилучшего среднеквадратичного приближения
(37)
здесь dj — среднеквадратичное уклонение.
Выберем линейную аппроксимацию
(38)
с числом членов n<N. Тогда
(39)
и необходимые условия минимума (37) запишутся в виде системы линейных уравнений относительно искомых констант аппроксимации ak:
(40)
или
(41)
Можно показать, что эти условия минимума являются и достаточными.
Весовая функция r(xi) должна давать больший вес в тех узлах, в которых нужна большая точность. Обычно где ei — абсолютная погрешность. Если погрешность по всем узлам одинакова, то ri=1.
Определенную по всем узлам погрешность интерполяция оценивает по среднеквадратичному уклонению согласно формуле (37):
Определение погрешности Погрешность ошибка некоторой величины это разность... Способы оценки погрешности... Различают априорную и апостериорную оценки погрешности...
Определение погрешности
Физические величины почти всегда известны приближенно, с погрешностью, если только речь не идет о целочисленных величинах или физических эталонах.
Погрешность (
Источники погрешности
В теоретических исследованиях выделяют следующие источники погрешности:
- погрешность исходных данных;
- погрешность математической модели
Задачи аппроксимации и интерполяции
Часто функции либо имеют очень громоздкое аналитическое выражение, либо заданы таблично. В этом случае имеет смысл на некотором интервале заменить заданную функцию y(x) приближенным аналитическим,
Интерполяционный многочлен Ньютона.
Наиболее проста и универсальна интерполяция степенными функциями
. (2.7)
Можно показать, что условие (6) при этом всегда в
Погрешность и трудоемкость интерполяции
Известные в литературе априорные оценки погрешности интерполяции на практике неприменимы, т.к. они требуют знания производных от функции y(x).
Для апостериорной оценки погрешнос
Нелинейная интерполяция
Полиномиальная интерполяция не всегда сходится. Ряд (17) расходится для быстро изменяющихся функций (для сетки с большим шагом). В этом случае следует применить метод выравнивания
Эрмитова интерполяция
Постановка задачи эрмитовой интерполяции: таблично задана функция y(x), а также ее производные:
Интерполяция сплайнами
Сплайн - непрерывная кусочно - полиномиальная функция. Ее степень в задачах интерполяции не зависит от числа узлов и поэтому сплайны эффективны при многоузловой
Двумерная интерполяция
Пусть функция z=z(x, y) задана таблично на прямоугольной сетке:
(43)
Необ
Полиномиальные формулы
Численное дифференцирование применяется, если функцию y(x) трудно или невозможно продифференцировать аналитически, например, если она задана таблично. Оно используется также при реш
Конечноразностные формулы
Часто требуется найти производные от функций, заданных на равномерной сетке и не в произвольной точке, а в узле сетки. Тогда можно получить формулы более простые, чем в общей постан
Метод Рунге - Ромберга
Этот метод позволяет, применяя формулы р-го порядка точности, получить результат (р+1)-го порядка точности.
Пусть приближенное значение некоторой величины на сетке с шагом
Вычисление интегралов
4.1. Квадратурные формулы Ньютона-Котеса
Даже если функция f(x) задана аналитически, вычислить интеграл
Формула средних
Самая простая квадратурная формула — одноузловая. Наиболее рационально выбирать этот узел посредине интервала, в точке `х=(a+b)/2. Для вывода формулы разложим подынтегральную функци
Формула трапеций
Если аппроксимировать подынтегральную функцию в (1) многочленом Лагранжа:
(10)
Формула Симпсона
Еще одна возможность получения квадратурных формул - использование метода Рунге-Ромберга (подраздел 3.3). Для вычисления интеграла
Нестандартные случаи интегрирования
Вычисление интеграла от функции, разрывной в точке х=с (а<с<в), возможно с высокой точностью, если интеграл предварите
Корректность задачи
Задача поставлена корректно, если:
1) решение задачи существует;
2) оно единственно,
3) решение непрерывно зависит от входных данных.
Вхо
Метод прогонки
Если матрица СЛАУ ленточная трехдиагональная, то метод Гаусса принимает более компактную форму и называется методом прогонки. СЛАУ при этом имеет следующий вид:
Метод LU-разложения
Матрица A СЛАУ
, (5.21)
если все главные миноры матрицы A отличны от ну
Метод квадратного корня
Если СЛАУ имеет симметричную матрицу, то для последней возможно представление
A = STDS,(5.29)
где S- верхняя треугольная матрица,
Итерационные методы решения СЛАУ
Итерационные методы решения СЛАУ заключаются в построении последовательности векторов (k=0,1,2,…), сходящейся к вектору
Итерационный метод вращений (Якоби)
Метод применим к симметричным матрицам и состоит в приведении заданной матрицы .к диагональному виду с помощью беско
Новости и инфо для студентов