
- •Численные методы (Вычислительная математика) Санкт-Петербург
- •Оглавление
- •Глава 1. Элементарная теория погрешностей........................10
- •Глава 2.Численные методы решения уравнений……………24
- •Глава 3. Приближенное представление функций.................50
- •Предисловие
- •Введение
- •Глава 1. Элементарная теория погрешностей
- •1.1. Приближенные числа и их погрешности
- •1.2. Правила приближенных вычислений и оценка погрешностей при вычислениях
- •1.3. Неустранимая погрешность вычисления значения функции
- •1.4. Обеспечение заданной точности результата
- •Глава 2.Численные методы решения уравнений
- •2.1. Итерационные методы решения нелинейных уравнений
- •2.2. Метод половинного деления
- •2.3. Метод хорд
- •2.4. Метод касательных (метод Ньютона).
- •2.5. Комбинированный метод хорд и касательных
- •2.6.Метод итераций
- •2.7. Решение систем линейных уравнений методом итераций.
- •2.9. Метод итераций для нелинейных систем уравнений
- •Глава 3. Приближенное представление функций (аппроксимация, интерполяция)
- •3.1. Интерполирование функции алгебраическими многочленами
- •3.2. Единственность интерполяционного многочлена
- •3.3. Интерполяционная формула Лагранжа
- •3.4. Погрешность интерполяционной формулы Лагранжа
- •3.5. Погрешность линейной интерполяции
- •3.6.Конечные разности и их свойства
- •3.7.Интерполяционная формула Ньютона (для равностоящих узлов)
- •3.10. Эмпирические формулы и метод наименьших квадратов
- •Глава 4. Численное дифференцирование и интегрирование
- •4.1. Дифференцирование на основе интерполяционной формулы Ньютона
- •4.2. Дифференцирование на основе формулы Стирлинга
- •4.3. Дифференцирование на основе интерполяционной формулы Лагранжа
- •Глава 5. Численное интегрирование функций
- •5.1. Квадратурные формулы Ньютона-Котеса
- •5.2. Формула трапеций
- •5.3. Формула Симпсона (формула парабол)
- •5.4. Формулы Ньютона – Котеса высших порядков
- •Глава 6. Численное интегрирование обыкновенных дифференциальных уравнений
- •6.1. Метод Эйлера
- •6.2. Метод Рунге – Кутта
- •6.3. Метод Адамса
- •Глава 7. Численные методы оптимизации
- •7.1. Градиентные методы
- •7.2. Метод Франка – Вулфа
- •7.3. Метод наискорейшего спуска
- •7.4. Метод штрафных функций
- •Литература
3.5. Погрешность линейной интерполяции
Вычислим погрешность интерполяционной формулы для наиболее часто употребляемого случая линейной интерполяции (при n=1)
В этом случае
Далее,
Для
функция
принимает
своё наибольшее значение в точке, где
.
Так как
,то
Положим
.Тогда
.
Вычислим неустранимую погрешность:
=
=
Неустранимая погрешность совпадает с погрешностью вычисления значений функции (табличной погрешностью).
Окончательно,
(5)
При использовании
таблиц интерполяция допустимо, чтобы
сумма
не превосходила одной единицы последнего
табличного разряда. Если используется
таблица с к
верными знаками после запятой, то должно
выполняться условие:
,
так как
,то
приходим к условию:
.
В случае линейной интерполяции это
условие принимает вид:
или
Вычисляя sin
,
при использовании линейной интерполяции
получим:
,то
есть
Следовательно,
Использованы
пятизначные таблицы, поэтому к=5 и нужное
условие
оказалось выполненным.
Полная погрешность равна:
3.6.Конечные разности и их свойства
Конечная разность первого порядка функции f(x) при шаге h>0 определяется равенством
(6)
Разности высших порядков определяются индуктивно: разность (к+1) – порядка есть разность первого порядка от разности к-го порядка:
Например,
Методом индукции можно установить общую формулу.
,
(7)
Значение разности
полностью определяется значениями
функции f(x)
в точках x,x+h,…,x+nh.
Вычисление разностей производится по таблице значений функции. При этом конечную разность в точке x записывают между значениями функции в точках x и x+h. Получение разностей показано в таблице:
x |
f(x) |
|
|
|
x |
|
|
|
|
x + h |
|
|||
x +2h |
|
|||
x +3h |
|
Свойства конечных разностей:
Конечная разность от постоянной величины равна нулю.
Операция вычисления конечной разности является линейной, т.е. для ,
R:
Доказательство:
Cвойство верно для разностей любого порядка.
Конечной разностью от многочлена степени n является многочлен степени (n-1), т.е. при вычислении разности, степень многочлена понижается на единицу.
С учётом предыдущего
свойства достаточно доказать это
свойство для
.
Имеем:
Конечная разность n-го порядка от многочлена n -ой степени постоянна, все разности более высокого порядка равны нулю.
Доказательство свойства следует из предыдущего.
5. Для многочлена
справедлива
формула:
Применяя последовательно формулу из доказательства свойства 3, получим
Следовательно,
для многочлена
получим
Поскольку по свойству 4:
Без доказательства сформулируем следующее свойство.
Если функция f(x) имеет непрерывную в точке x производную к-ого порядка, то
Следовательно,
при малых справедливо приближенное
равенство
Приведённая формула
дает грубое приближение вследствие
малого числа
в знаменателе. Например, если вычисляется
производная
и
погрешность вычисления значения функции
равна
,
то возникающая из-за этого неустранимая
погрешность равна
.
При h=0,1;0,01;0,001
получаются
ошибки
.