
- •Смагин в.И., Решетникова г.Н.
- •Содержание
- •1. Основы теории погрешностей 9
- •2. Аппроксимация функций 19
- •2.16. Контрольные вопросы 90
- •2.17. Задания к главе 2 92
- •3. Численное дифференцирование 95
- •3.5. Контрольные вопросы 110
- •4. Численное интегрирование 112
- •4.10. Контрольные вопросы 189
- •4.11. Задания к главе 4 191
- •5. Приложение. Варианты к заданиям 195
- •1. Основы теории погрешностей
- •1.1. Математические оценки точности приближенного числа
- •1.2. Запись чисел на эвм
- •1.3. Верные знаки приближенного числа
- •1.4. Классификация погрешностей
- •1.5. Погрешность вычисления функции многих переменных
- •1.6. Обратная задача теории погрешностей
- •1.7. Погрешности простейших функций
- •1.8. Контрольные вопросы
- •1.9. Задания к главе 1
- •2. Аппроксимация функций
- •2.1. Постановка задачи
- •2.2. Интерполяционный многочлен Лагранжа
- •2.3. Интерполяционная схема Эйткена
- •2.4. Остаточный член многочлена Лагранжа
- •2.5. Разделенные разности и их свойства
- •2.6. Интерполяционная формула Ньютона при неравноотстоящих узлах
- •2.7. Многочлены Чебышева и их свойства
- •2.8. Минимизация погрешности метода при аппроксимации многочленом Лагранжа
- •2.9. Многочлены наилучшего равномерного приближения
- •2.10. Экономизация степенных рядов
- •2.11. Интерполирование с кратными узлами
- •2.12. Интерполирование при равноотстоящих узлах
- •2.12.1. Конечные разности
- •2.12.2. Интерполирование в начале и конце таблицы
- •2.12.3. Формулы Гаусса
- •2.12.4. Формулы Стирлинга и Бесселя
- •2.12.5. Оценки погрешности метода и неустранимой погрешности
- •2.13. Сплайн-функции
- •2.13.1. Линейный сплайн
- •2.13.2. Параболический сплайн
- •2.13.3. Кубический сплайн
- •2.13.4. В-сплайны
- •2.13.5. Эрмитовы сплайны
- •2.14. Аппроксимация данных методом наименьших квадратов (мнк)
- •2.14.1. Аппроксимация алгебраическими полиномами
- •2.14.2. Аппроксимация ортогональными полиномами
- •2.14.3. Аппроксимация ортогональными полиномами дискретной переменной
- •2.15. Аппроксимация функций многих переменных
- •2.15.1. Построение интерполяционных многочленов
- •2.15.2.Метод последовательного интерполирования
- •2.15.3. Применение метода наименьших квадратов
- •2.16. Контрольные вопросы
- •2.17. Задания к главе 2
- •3. Численное дифференцирование
- •3.1. Численное дифференцирование при неравноотстоящих узлах
- •3.2. Численное дифференцирование при равноотстоящих узлах
- •3.3. Оценка приближений численного дифференцирования по правилу Рунге
- •3.4. Метод квадратурных формул
- •3.5. Контрольные вопросы
- •3.6. Задание к главе 3
- •4. Численное интегрирование
- •4.1. Общая интерполяционная квадратура
- •4.2. Квадратурные формулы Ньютона-Котеса
- •4.2.1. Квадратурные формулы прямоугольников
- •4.2.2. Квадратурная формула трапеций
- •4.2.3. Квадратурная формула Симпсона (парабол)
- •4.2.4. Квадратурная формула “трех восьмых” (формула Ньютона)
- •4.3. Метод Рунге оценки погрешности
- •4.4. Квадратурные формулы наивысшей алгебраической степени точности
- •4.5. Частные случаи квадратурного правила наивысшей алгебраической степени точности
- •4.6. Квадратурные формулы с равными коэффициентами
- •4.7. Приближенное вычисление несобственных интегралов
- •4.8. Приближенное вычисление неопределенных интегралов
- •4.9. Методы Монте-Карло
- •4.9.1. Простейший метод Монте-Карло
- •4.9.2. Геометрический метод Монте-Карло
- •4.10. Контрольные вопросы
- •4.11. Задания к главе 4
- •5. Приложение. Варианты к заданиям
- •5.1. Варианты к заданиям 1.1
- •5.2. Варианты к заданиям 2.1-2.5
- •5.3. Варианты к заданиям 2.6
- •5.4. Варианты к заданиям 2.7-2.8, 3.1
- •5.5. Варианты к заданиям 4.1-4.6
- •5.6. Варианты к заданиям 4.7
2.7. Многочлены Чебышева и их свойства
Многочлены Чебышева степени вычисляются по формуле
.
(2.30)
Они
определены на отрезке
.
Применяя формулы тригонометрии, вычислим
многочлены Чебышева. Для этого введем
обозначение
,
тогда имеем:
,
,
,
,…,
и т.д.
Используя тригонометрическое тождество:
,
получим рекуррентное соотношение для вычисления многочленов Чебышева в виде
,
(2.31)
со следующими начальными условиями
.
Многочлены Чебышева обладают следующими свойствами:
1) коэффициент при
старшей степени
равен
(свойство следует из соотношения (2.31));
2) корни многочленов Чебышева вычисляются по формуле
,
,
(2.32)
(определяются
из уравнения
);
3) многочлены Чебышева имеют точку экстремума
,
,
(2.33)
(определяются
из уравнения
);
4)
;
5) многочлены
Чебышева ортогональны с весом
:
;
(2.34)
6) многочлены
Чебышева являются наименее отклоняющимися
от нуля на отрезке
.
Обоснование последнего свойства следует из теоремы.
Теорема 2.2.
Для любого многочлена
степени
с единичным коэффициентом при старшей
степени справедливо неравенство
Доказательство. Допустим противное
Тогда
разность
является многочленом степени n-1.
В n+1 точках многочлен
принимает попеременно значения
,
а так как по предположению
то разность попеременно принимает, то
положительное, то отрицательное значение
в
точке. Таким образом, получается, что
разность
пересекает на интервале
ось абсцисс
раз, то есть имеет
корней. Но это противоречие, так как
разность является многочленом степени
и должна иметь
корень. Теорема доказана.
Итак, мы доказали, что из всех многочленов степени n с единичным коэффициентом при старшей степени, точная верхняя грань на интервале будет наименьшей у многочлена Чебышева . Поэтому, как указывалось в 6-ом свойстве, многочлены Чебышева называются многочленами, наименее отклоняющими от нуля.
Графики многочленов Чебышева 1, 2, 3 и 10 порядка приведены на рис.2.1.
Рис. 2.1. Графики многочленов Чебышева
Как видно из графиков корни многочлена Чебышева располагаются на интервале чаще вблизи концов интервала.
2.8. Минимизация погрешности метода при аппроксимации многочленом Лагранжа
Рассмотренные выше способы вычисления многочлена Лагранжа – формула Лагранжа, схема Эйткена и формула Ньютона, дают одинаковую погрешность метода, которую можно оценить, зная выражение для остаточного члена формулы Лагранжа.
Пусть функция , подлежащая интерполированию, на интервале имеет непрерывные производные до порядка включительно. Тогда, как показано в п. 2.4, остаточный член многочлена Лагранжа имеет вид:
(2.35)
где – некоторое число, принадлежащее интервалу . Анализируя остаточный член (2.35) видно, что при неудачном расположении узлов интерполирования значение
может оказаться очень большим. Поэтому возникает задача выбора узлов интерполирования (при заданном n) таких, чтобы величина
была
минимальной или иначе многочлен
был бы наименее отклоняющимся от нуля
на отрезке [
].
Как отмечалось в п. 2.8, наименее
отклоняющимися от нуля являются
многочлены Чебышева. Многочлен
имеет единичный коэффициент при старшей
степени и для этого многочлена справедливо
Отсюда вытекает,
что если требуется минимизировать
погрешность метода при интерполировании
на интервале
с помощью многочлена Лагранжа, то в силу
(2.35), необходимо выбрать многочлен
так, чтобы
была наименьшей. А это можно достигнуть,
если в качестве
взять многочлен Чебышева
,
который является наименее отклоняющимся
от нуля. При этом в качестве узлов
интерполирования необходимо взять
корни многочлена Чебышева степени
(2.36)
Погрешность метода при интерполировании на интервале будет минимальной и определится по формуле:
(2.37)
В случае интерполирования на произвольном интервале вводится линейное преобразование
,
(2.38)
где
.
Тогда узлы интерполирования на
будут следующие:
,
(2.39)
Для оценки
погрешности метода при интерполировании
на интервале [
],
необходимо найти верхнюю грань
.
Оценим эту величину, учитывая свойства
многочленов Чебышева:
.
(2.40)
Тогда в силу (2.40), погрешность метода при интерполировании на интервале [ ] определится по формуле
(2.41)