
- •Смагин в.И., Решетникова г.Н.
- •Содержание
- •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
4.6. Квадратурные формулы с равными коэффициентами
Квадратурное правило, все коэффициенты которого одинаковы,
(4.30)
особенно удобно при графических расчетах, так как сумму ординат можно достаточно легко снять с чертежа.
Будем считать, что весовая функция не эквивалентна нулю, т.е.
.
Формула (4.30) содержит
параметр:
Эти параметры будем выбирать таким
образом, чтобы правило (4.30) выполнялось
точно для многочленов степени
,
что равносильно выполнению равенств
(4.31)
Правила интегрирования, обладающие этим свойством, называются квадратурными правилами Чебышева.
Теорема 4.9. Если , то квадратурное правило (4.30), точное для всех алгебраических многочленов степени не выше , с действительными или комплексными узлами всегда может быть построено и при этом единственным способом.
Доказательство.
Условие точного выполнения равенства
(4.30) для многочлена нулевой степени даст
уравнение для нахождения коэффициента
и
. (4.32)
Если равенство
(4.31) записать для
,
то для нахождения значений
получится система уравнений
(4.33)
Полученная система
является нелинейной и процесс нахождения
ее решения достаточно сложный. Поэтому
будем искать не узлы
,
а строить многочлен, для которого эти
узлы будут корнями
.
Коэффициенты
,
являются симметрическими функциями
корней
,
(соотношениями Виетта) и
также являются симметрическими функциями
корней, причем справа в (4.33) указаны их
значения. Из алгебры известны соотношения,
связывающие
и
,
(4.34)
Соотношения (4.34)
позволяют по известным значениям
найти последовательно и при этом
единственным образом значения
.
Затем по коэффициентам
,
можно построить многочлен
и,
решив уравнение
,
найти узлы
,
квадратурного правила (4.30).
Замечание
4.3. Бернштейном
доказано, что при
и для всех
,
среди
есть комплексные.
Рассмотрим частный
случай квадратурного правила Чебышева,
когда весовая функция
и интеграл имеет вид
.
Так как всякий конечный отрезок интегрирования линейной заменой переменных
преобразуется в отрезок , то квадратурное правило записывается в виде:
.
Будем считать, что отрезок интегрирования приведен к интервалу и квадратурное правило имеет вид
.
Так как
,
то
и
.
Значения
,
согласно (4.34), определяются следующим
образом:
.
Тогда уравнения (4.34) запишутся в виде
,
,
,
,
Значения с нечетными индексами получились равными и в многочлене сохраняются только либо одни четные, либо одни нечетные степени
.
В таблице 4.5.
приведены значения узлов
правила Чебышева для
.
Таблица 4.5.
|
|
|
|
|
|
|
|
|
|
|
= 0,0000000 |
|
|
|
|
Замечание 4.4. Можно уточнить значение интеграла путем деления интервала интегрирования на равных частей. В этом случае используется только табличных узлов, а значение интеграла вычисляется по узлам.
Квадратурная формула Чебышева для вычисления интеграла путем деления интервала интегрирования на равных частей имеет вид
.
Пример 4.14. Вычислим интеграл
.
С помощью замены переменной
интеграл и соответствующая квадратурная формула преобразуются к виду
Если
вычислять интеграл при
,
используя данные таблицы 4.5., то получим
.
Вычисляя интеграл при тех же значениях
путем деления интервала
на две части
,
получим
,
а при
−
.