- •Содержание
- •Глава 1. Численные методы в теории приближений
- •Глава 2. Численное интегрирование.
- •Глава 3. Численные методы алгебры.
- •Глава 4. Численные методы решения обыкновенных дифференциальных уравнений.
- •Предисловие
- •Глава 1. Численные методы в теории приближений.
- •1.1. Структура погрешности в численном анализе.
- •1.3 Понятие близости в метрическом пространстве.
- •1.2 Распространение ошибок округления в арифметических операциях.
- •1.4. Задача интерполяции как простейшая задача приближения функций.
- •1.5 Конечные разности и их свойства.
- •1.6 Интерполяционный полином Ньютона
- •1.7 Многочлены Чебышева, их свойства и применение в задаче интерполяции.
- •1.7.1. Основные определения.
- •1.7.2. Простейшие свойства многочленов Чебышева.
- •1.7.3. Применение многочленов Чебышева в задаче интерполяции.
- •1.8 Среднеквадратичное приближение функций.
- •1.8.1. Общая постановка задачи и ее разрешимость.
- •1.8.2. Среднеквадратичное приближение функций алгебраическими многочленами.
- •1.8.3. Среднеквадратичная ошибка аппроксимации полиномами Лежандра.
- •Глава 2. Численное интегрирование.
- •2.1. Использование функциональных рядов.
- •2.2. Квадратурные формулы на основе интерполяции.
- •2.3. Квадратурные формулы Ньютона-Котеса.
- •2.4. Некоторые общие свойства ортогональных с весом полиномов.
- •2.5.Квадратурные формулы Гаусса-Кристоффеля.
- •Глава 3. Численные методы алгебры.
- •3.1. Принцип сжатых отображений.
- •3.2. Метод простых итераций для функциональных уравнений.
- •3.3. Метод Ньютона.
- •3.4. Численные методы решения систем лау.
- •3.4.1. Прямые методы решения систем лау.
- •3.4.2. Нормы векторов и матриц.
- •3.4.3. Обусловленность матриц и систем уравнений.
- •3.4.4. Итерационные методы решения систем лау.
- •3.4.5. Стационарные итерационные процедуры.
- •Глава 4. Численные методы решения обыкновенных дифференциальных уравнений (оду).
- •4.1. Численное дифференцирование на основе интерполяции.
- •4.2. Численное дифференцирование на равномерной сетке.
- •4.3. Задача Коши для оду.
- •4.3.1. Постановка задачи.
- •4.3.2. Метод Эйлера и его модификации.
- •4.4. Численные методы решения краевых задач для оду.
- •4.4.1. Постановка задачи для диффернциального уравнения 2-го порядка.
- •4.4.2. Метод конечных разностей (метод сеток).
- •4.4.3. Аппроксимация, устойчивость и сходимость разностных схем.
- •, Где с не зависит от h (т.Е. От n). Без доказательства (см.[1,2]).
- •Литература
1.7 Многочлены Чебышева, их свойства и применение в задаче интерполяции.
1.7.1. Основные определения.
Определение 1.
Говорят, что
функция
,
если
.
При этом
называется весовой
функцией и
удовлетворяет условиям:
на
и
.
Определение 2.
Функции
и
называются ортогональными
на
с весом
,
если их скалярное произведение
.
Замечание.
Из неравенства Коши-Буняковского-Шварца
для интегралов следует, что скалярное
произведение
существует
Определим на отрезке [-1,1] следующие многочлены Чебышева:
|
(24) |
Найдем два первых многочлена Чебышева по формуле (24):
Для больших n
неудобно работать с формулой (24). Выведем
более удобную рекуррентную
формулу.
Полагая
и подставляя в формулу тригонометрии:
,
получаем:
|
(25) |
Формула (25) начинает
работать, начиная со значений
.
Последовательно получаем:
1.7.2. Простейшие свойства многочленов Чебышева.
1. Многочлены
Чебышева ортогональны на отрезке
с весом
.
Рассмотрим интеграл
в силу ортогональности
системы функций
на отрезке
.
Вычислим квадрат нормы:
.
2. При четных (нечетных) n многочлен Чебышева Tn(x) содержит только четные (нечетные) степени х, т.е. является четной (нечетной) функцией.
Доказывается по индукции с помощью рекуррентной формулы (25).
3. Коэффициент при старшей степени xn многочлена Tn(x) равен 2n-1.
Доказывается по индукции с помощью рекуррентной формулы (25).
4. Многочлен Tn(x)
имеет на отрезке
ровно n
различных действительных корней,
определяемых формулой:
Действительно:
5.
и достигается в точках экстремума:
.
Из определения
(1) следует, что
для любого
.
Очевидно, что
.
6. Многочлен
среди всех многочленов
n-ой
степени с коэффициентом при старшей
степени an=1
обладает тем свойством, что
|
(26) |
Доказывается от противного: пусть это не так и существует многочлен
такой, что выполняется противоположное:
|
(27) |
Разность (
)
– многочлен (n-1)-ой
степени, причем в силу (4)
.
Обозначим
и заметим что, в силу (27)
|
|
Продолжим рассмотрение разности:
Таким образом, при
переходе от точки
к
разность (
)
меняет знак.
Всего при переходе от точки
к
произойдет ровно n
смен знака. Отсюда следует, что разность
имеет на отрезке
ровно n
действительных корней (нулей), что
противоречит теореме Гаусса, т.к. это
многочлен (n-1)-ой
степени.
Замечание.
Благодаря
свойству 6 многочлен Чебышева
называется многочленом,
наименее отклоняющимся от нуля.
1.7.3. Применение многочленов Чебышева в задаче интерполяции.
В следующей теореме решается следующая задача: как следует выбрать узлы интерполяции, чтобы минимизировать погрешность на всем отрезке.
Теорема 1.3. Пусть
,
тогда наименьшая максимальная абсолютная
погрешность интерполяции полиномом
Лагранжа
на отрезке
достигается при выборе в качестве узлов
интерполяции нулей функции
Обозначим
- корень многочлена
.
Согласно свойству 4
|
(28) |
Пусть
– некоторая система узлов на
.
Запишем формулу максимальной абсолютной
погрешности интерполяции по Лагранжу
(формула (10) тз п.п.1.4):
,
где
,
Как следует из свойства 6:
.
Выберем в качестве
узлов точки
,
определяемые формулой (28). Тогда
.
Отметим, что
многочлен
имеет одну и ту же степень n+1,
что и
,
один и тот же коэффициент при старшей
степени
и одни и те же нули на
.
Отсюда немедленно
следует, что
,
и соответствующая оценка погрешности:
.
Данная оценка является наилучшей среди всех возможных способов выбора узлов интерполяции.
Замечание. Для оптимальной интерполяции на произвольном конечном отрезке [a;b] предварительно необходимо сделать линейное преобразование:
и преобразовать
формулу для нулей функции
к следующему виду:
,
Пример
15. Вывести
следующую формулу для максимальной
абсолютной погрешности интерполяции
по Лагранжу на отрезке
:
