
- •Численные методы (Вычислительная математика) Санкт-Петербург
- •Оглавление
- •Глава 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. Метод штрафных функций
- •Литература
Глава 1. Элементарная теория погрешностей
1.1. Приближенные числа и их погрешности
В самых разнообразных теоретических и прикладных исследованиях широко используются методы математического моделирования, которые сводят решение задач данной области исследования к решению адекватных (или приближенно адекватных) им математических задач. Необходимо довести решение этих задач до получения числового результата (вычисления различного рода величин, решения различных типов уравнений и т.п.). Целью вычислительной математики является разработка алгоритмов численного решения обширного круга математических задач. Методы должны быть разработаны так, чтобы их можно было эффективно реализовать с помощью современной вычислительной техники. Как правило, рассматриваемые задачи не допускают точного решения, поэтому речь идет о разработке алгоритмов, дающих приближенное решение. Для возможности замены неизвестного точного решения задачи приближенным необходимо, чтобы последнее было достаточно близко к точному. В связи с этим возникает необходимость оценки близости приближенного решения к точному и разработки приближенных методов построения приближенных решений, сколько угодно близких к точным.
Схематически вычислительный процесс заключается в том, чтобы для данной величины x (числовой, векторный и т.д.) вычислить значение некоторой функции A(x). Разность между точным и приближенным значениями величины называют погрешностью. Точное вычисление значения A(x) обычно невозможно, и вынуждает заменить функцию (операцию) A ее приближенным представлением Ã, которое можно вычислить: вычисление величины A(x), заменяется вычислением- Ã(x). Возникающую при этом погрешность A(x) - Ã(x) называют погрешностью метода. Способ оценки этой погрешности должен быть разработан вместе с разработкой метода вычисления величины Ã(x). Из возможных методов построения приближения следует использовать тот, который при имеющихся средствах и возможностях дает наименьшую погрешность.
Значение величины
x,
то есть исходных данных, в реальных
задачах получается или непосредственно
из измерений, или в результате предыдущего
этапа вычислений. В этих случаях
определяется лишь приближенное значение
x
o
величины
x.
Поэтому вместо величины Ã(x)
может быть вычислено лишь приближенное
ее значение Ã(xo).
Возникающую при этом погрешность A(x)
- Ã(xo)
называют неустранимой.
В результате неизбежных в ходе вычислений
округлений, вместо величины Ã(xo)
вычисляется ее «округленное» значение
,
что приводит к появлению погрешности
округления
Ã(xo)
-
.
Полная погрешность вычислиниия
оказывается равной A(x)
-
.
Представим полную погрешность в виде
A(x)
-
= [A(x)
-
]
+
[
- Ã(xo)]
+
+ [Ã(xo) - ] (1)
Последнее равенство показывает, что полная погрешность вычисления равна сумме погрешности метода, неустранимой погрешности и погрешности округления. Первые две составляющие погрешности можно оценить до начала вычислений. Погрешность округления оценивается лишь в ходе вычислений.
Рассмотрим следующие задачи:
а) характеристика точности приближенных чисел
б) оценка точности результата при известной точности исходных данных (оценка неустранимой погрешности)
в) определение необходимой точности исходных данных для обеспечения заданной точности результата
г) согласование точности исходных данных и вычислений с возможностями имеющихся вычислительных средств.