
- •Смагин в.И., Решетникова г.Н.
- •Содержание
- •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.9. Многочлены наилучшего равномерного приближения
Введем абсолютное отклонение аппроксимирующего многочлена от аппроксимируемой на интервале непрерывной функции
.
(2.42)
Если многочлен выбран отрезке так, что отклонение минимально, то называется многочленом наилучшего равномерного приближения. Точки, в которых реализуется максимальное отклонение многочлена от называются точками чебышевского альтернанса.
Отметим, что существование и единственность многочлена наилучшего равномерного приближения на ограниченном отрезке может быть строго доказана. Однако, ни общий вид многочленов наилучшего равномерного приближения, ни способы их построения, неизвестны. Имеются лишь методики построения многочленов наилучшего равномерного приближения для достаточно узкого класса функций и некоторые способы построения многочленов близких к многочленам наилучшего равномерного приближения.
Пример 2.5. Пусть
для непрерывной функции
на интервале
требуется построить многочлен наилучшего
равномерного приближения нулевой
степени
.
Для того чтобы найти значение
воспользуемся свойством непрерывной
на ограниченном замкнутом интервале
функции, согласно которому на этом
интервале найдутся по крайней мере две
точки, в которых функция принимает
максимальное и минимальное значение.
Введем обозначения:
,
.
Тогда
,
при этом абсолютное отклонение определится по формуле
.
Пример 2.6. Пусть
для непрерывной функции
на интервале
требуется построить многочлен наилучшего
равномерного приближения степени
.
В силу определения наилучшего равномерного
приближения многочлен
должен быть наименее отклоняющимся от
нуля на интервале
.
Но такой многочлен степени
нам известен это многочлен Чебышева
.
Тогда из равенства
=
легко определить многочлен :
.
В
частном случае при
имеем:
.
Таким
образом получаем, что наилучшим
равномерным приближением для
на интервале
является многочлен третьей степени.
Сравним точность
построенного многочлена наилучшего
равномерного приближения
с многочленом Лагранжа 3-ей степени,
построенного по чебышевским узлам
(2.36) (см. п. 2.9). Для нашего случая этот
многочлен будет иметь вид:
.
На рис 2.2. приведен график ошибок аппроксимации
для многочлена наилучшего равномерного приближения (сплошная линия) и график ошибки аппроксимации многочленом Лагранжа
(пунктирная
линия). Как видно из графиков точность
аппроксимации на интервале
с использованием многочлена наилучшего
равномерного приближения выше и величина
максимальной ошибки
в 2 раза меньше чем при аппроксимации
многочленом Лагранжа (максимальная
ошибка достигается на концах интервала
и равна 0,125, см. рис. 2.2).
Рис. 2.2. Графики ошибок аппроксимации
2.10. Экономизация степенных рядов
В некоторых задачах достаточно просто получить аппроксимацию функции с помощью ряда Тейлора вида:
,
сходящегося
при
.
Тогда может быть применен следующий
метод построения аппроксимирующей
функции:
1) подбирается значение , такое, чтобы многочлен
,
аппроксимировал
функцию
с погрешностью, не превышающую величину
;
2) степень
многочлена
понижается на единицу посредством
замены
на наилучшее равномерное приближение
степени
.
Выполнив эти два этапа, получим аппроксимирующий многочлен степени :
.
Погрешность аппроксимации функции определится на интервале оценкой
Такой подход позволяет осуществить экономизацию степенного ряда, не снижая значительно точность аппроксимации. Если полученная оценка позволяет дальнейшую экономизацию степенного ряда, можно попытаться понизить степень аппроксимирующего многочлена еще на единицу.
Экономизацию
степенного ряда можно осуществлять
также с помощью замены аппроксимирующего
многочлена
эквивалентным разложением по многочленам
Чебышева
.
Для этого можно воспользоваться
выражениями степеней
через многочлены Чебышева:
,
,
,
,
,
и т.д.
Отметим, для достижения одной и той же точности в отрезке разложения по многочленам Чебышева можно брать, как правило, меньшее число членов, чем в степенной аппроксимации.