
- •2202 “Автоматизированные системы обработки
- •Приближенные числа.
- •Действия над приближенными числами.
- •Общая схема метода Гаусса для систем, имеющих единственное решение
- •Обратный ход:
- •Метод Гаусса с выбором максимального элемента по столбцу
- •Достаточное условие сходимости процесса итераций
- •Основные понятия алгебры матриц и теории линейных векторных пространств.
- •Обратная матрица
- •Правило нахождения ранга матрицы:
- •Нахождение обратной матрицы методом Гаусса
- •Применение метода итераций для умножения элементов обратной матрицы
- •Пример:
- •Приближенное решение алгебраических и трансцендентных уравнений
- •Пример:
- •Блок – схема решения уравнения f(X) методом половинного деления
- •Метод хорд
- •Метод Ньютона
- •Пример:
- •Приближенное решение алгебраических и трансцендентных уравнений (продолжение) Метод итераций
- •Теорема
- •Оценка приближения.
- •Интерполирование функций.
- •1. Конечные разности различных порядков.
- •2) Таблица разностей.
- •3) Постановка задачи интерполирования.
- •4) Первая интерполяционная формула Ньютона.
- •Блок-схема построения кубического сплайна
- •Численное дифференцирование
- •2.Конечно-разностные аппроксимации производных.
- •Блок-схема вычисления производной.
- •Численное интегрирование.
- •Блок-схема вычисления определенного интеграла методом двойного пересчета по формуле Симпсона.
- •Квадратурные формулы Гаусса
- •Приближенное вычисление несобственных интегралов с бесконечными пределами.
- •Среднеквадратичное приближение функций
- •Точечное квадратичное аппроксимирование функций
- •Условия нахождения экстремума функции
- •IV Интегральное квадратичное аппроксимирование функций на отрезке.
- •Ортогональные на промежутке системы функций
- •Основные понятия гармонического анализа.
- •4.2.1. Аппроксимация тригонометрическими полиномами кусочно-гладких функций
- •4.3. Среднеквадратичное приближение функций полиномами Лежандра.
- •4.3. Полиномы, наименее уклоняющиеся от нуля.
- •Эмпирические формулы
- •Численное решение обыкновенных дифференциальных уравнений.
- •Численное решение дифференциальных уравнений и систем дифференциальных уравнений высших порядков
- •Численные методы безусловной оптимизации.
IV Интегральное квадратичное аппроксимирование функций на отрезке.
До сих пор мы рассматривали аппроксимацию точечных множеств. Теперь рассмотрим аппроксимацию ф-ий на отрезке (непрерывные множества).
За меру отклонения ф-ии f(x)
от полинома
на отрезке [a,b]
принимается величина
Решение ищется из условия минимизации
т.е.
приравниваются нулю все частные
производные
.
В результате получаем систему уравнений:
(1)
где
(k=0,1,2,…,m)
Решение
единственное
решение, обеспечивающее минимум.
Пример: Найти наилучшую квадратичную
аппроксимацию посредством полинома
второй степени ф-ии
на отрезке [0,1]. Решение:
Система
(1) примет вид:
С=[6/35; 48/35; -4/7];
;
Q(0)=6/35; Q(1)=6/35+48/35-4/7=34/35
Если
-
обобщенный полином, то
Система (1) в этом случае примет вид:
C0 (φ0, φ0) + C1 (φ1, φ0) +…+ Cn (φn, φ0) = (f, φ0)
C0 (φ0, φ1) + C1 (φ1, φ1) +…+ Cn (φn, φ1) = (f, φ1) (2)
C0 (φ0, φn) + C1 (φ1, φn) +…+ Cn (φn, φn) = (f, φn)
Где
Неудобство интегральной квадратичной аппроксимации - необходимость вычисления определенных интегралов.
В этом случае способ точечной квадратичной аппроксимации предпочтительнее.
-
Ортогональные на промежутке системы функций
Если
,
при m<>n,
То
-
система ф-ий, ортогональных (1) на [a,b]
-
норма ф-ий
на [a,b]
Если нормы всех ф-ий системы (1) равны единицы, то это система называется ортонормированной.
Для ортонормированных систем
Чтобы любую систему, не содержащую функцийй с нулевой нормой, пронормировать, необходимо каждую функцию разделить на ее норму:
Если система ф-ий
ортогональна на отрезке [a,b],
то коэффициенты
обобщенного полинома
,
аппроксимирующего непрерывную ф-ию
f(x)на [a,b]
Имеют вид:
(2)
-коэфф.
Фурье ф-ии f(x)
относительно заданной ортогональной
системы
.
Обобщенный полином с коэффициентом Фурье данной ф-ии обладает наименьшим квадратичным отклонением от этой ф-ии по сравнению со всеми другими обобщенными полиномами того же порядка m.
(3)
После преобразований (раскрывая скобки,
меняя местами
и
и
приводя подобные члены) получим(без
вывода):
(4)
(5)
Из (3) следует, что
,
потому из (5) получаем:
(6) –неравенство Бесселя
При
m ∞
(7)
Если система
-ортонормированная,
то
(8)
Если
,
то система
наз. ПОЛНОЙ.
Для полной ортонормированной системы имеет место равенство Парсеваля
(9)
Свойства обобщенного полинома
с коэффициентами Фурье:
1. при увеличении числа слагаемых m
младшие коэффициенты
остаются неизменными, т.е. при
добавлении новых членов прежние
коэффициенты не пересчитываются (это
следует (2)).
2. при увеличении m
квадратичная погрешность
монотонно
убывает в широком смысле, т.е.
Т.о.
присоединение новых слагаемых увеличивает
точность аппроксимации.
-
Основные понятия гармонического анализа.
Тригонометрическая система функций:
1, sin x, cos x, sin 2x, cos 2x,…, sin nx, cos nx (1)
ортогональна на любом отрезке длины 2π (напрмер, на [-π, π]).
Нормы функций системы (1)
n=1,2…. (2)
Пусть дана непрерывная периодическая функция с периодом 2π (введением новой переменной можно область определения функции [a, b] перевести в интервал [-π,π])
Составим тригонометрический полином
(3)
Слагаемые
k=1,2…, называются
гармонитами.
Чтобы минимизировать
min
Коэффициенты
,
,
должны быть коэффициентами Фурье функции
f(х) относительно системы
(1)
т.е.
Т.о. получаем:
(4)
(k=0,1,2,…m)
Полином (3) – тригонометрический полином Фурье;
,
- тригонометрические коэффициенты Фурье
функции f(х).
Если f(х) четная, то
(k=1,2,…,m)
k=0,1,2…m
(5)
Если f(х) нечетная, то
(k=0,1,2,…,m)
(k=1,2,…,m)
(6)
Для четной функции
-
При m→∞ получаем тригонометрический ряд Фурье
Представленные
функции тригонометрическим полиномом
Фурье или тригонометрическим рядом
Фурье называется гармоническим анализом.
В простейших случаях коэффициенты тригонометрического полинома Фурье вычисляются непосредственно по формулам (4)
Среднеквадратическое отклонение
определено как
,
в общем случае функция
задана на интервале
.