- •Министерство образования российской федерации
- •Московский международный институт эконометрики, информатики, финансов и права
- •И.Н. Мастяева о.Н. Семенихина
- •Численные методы
- •Учебное пособие
- •Москва 2004
- •Содержание:
- •1. Погрешность результата численного решения задачи
- •1.1. Источники и классификация погрешностей.
- •1.2. Точные и приближенные числа. Правила округления чисел
- •1.3. Математические характеристики точности приближенных чисел
- •1.4. Число верных знаков приближенного числа. Связь абсолютной и относительной погрешности с числом верных знаков. Правила подсчета числа верных знаков
- •5423,47 6 Значащих цифр,
- •0,0000605 3 Значащие цифры,
- •0,060500 5 Значащих цифр.
- •1.5. Общая формула теории погрешностей (погрешность вычисления значения функции)
- •1.6. Погрешность арифметических действий
- •1.7. Обратная задача теории погрешностей
- •2. Численные методы решения нелинейных уравнений
- •2.1. Отделение корней
- •2.2. Метод половинного деления
- •2.3. Метод хорд (секущих)
- •2.4. Метод касательных (метод Ньютона)
- •2.5. Метод итераций
- •3. Численные методы линейной алгебры
- •3.1. Метод Гаусса
- •З.2. Метод прогонки
- •3.3. Норма вектора и норма матрицы
- •3.4. Метод простой итерации
- •3.5. Частичная проблема собственных значений
- •Интерполирование.
- •4.1. Интерполяционный полином, его существование и единственность. Остаточный член.
- •4.2. Интерполяционный полином Лагранжа.
- •4.3. Разделенные разности и их свойства.
- •4.4. Интерполяционный полином Ньютона с разделенными разностями
- •4.5. Конечные разности и их свойства
- •4.6. Интерполяционные формулы Ньютона
- •4.7. Интерполяционные полиномы с центральными разностями
- •4.8.Обратное интерполирование
- •4.9. Численное дифференцирование
- •5. Интерполирование с кратными узлами и сплайны
- •5.1. Разделенные разности с повторяющимися (кратными) узлами
- •5.2. Интерполяционный полином Эрмита
- •5.3. Интерполирование сплайнами
- •6. Численное интегрирование
- •6.1. Формула прямоугольников
- •6.2. Формула трапеций
- •6.3. Формула Симпсона
- •6.4. Правило Рунге практической оценки погрешности квадратурных формул. Уточнение приближенного значения интеграла по Ричардсону
- •7. Численные методы решения дифференциальных уравнений
- •7.1. Метод Рунге-Кутта
- •7.2. Разностный метод решения краевой задачи
- •Список литературы
1.7. Обратная задача теории погрешностей
Основная задача теории погрешностей заключалась в том, что по известным погрешностям аргументов находилась погрешность функции. На практике очень важное значение имеет и обратная задача: каковы должны быть погрешности аргументов, чтобы абсолютная погрешность функции не превышала заданной величины?
На основании общей формулы теории погрешностей имеем
. (1.32)
Задача отыскания допустимых значений абсолютной погрешности аргументов по известной абсолютной погрешности функции является математически неопределенной, так как в общем случае для определения п неизвестных мы имеем одно уравнение.
Обратная задача теории погрешностей имеет однозначное решение только для функции одного аргумента. Действительно, на основании общей формулы теории погрешностей
,
следовательно,
. (1.33)
-
Принцип равных влиянии. Согласно этому принципу предполагается, что все выражения
,
одинаково влияют на образование общей абсолютной погрешности, т.е.
.
Пусть нам задана абсолютная погрешность . На основании общей формулы теории погрешностей можно написать
.
Тогда
,
откуда
.
-
Принцип равных абсолютных погрешностей. Согласно этому принципу предполагается, что
x1=x2 =…=xn .
Тогда из общей формулы теории погрешностей будем иметь
или
.
-
Принцип равных относительных погрешностей. Согласно этому принципу предполагается, что
.
По определению , тогда .
Подставляя это выражение в общую формулу, получим
,
откуда
.
2. Численные методы решения нелинейных уравнений
2.1. Отделение корней
Рассмотрим некоторую функцию f(x).
Определение. Всякое число обращающее функцию в нуль, т.е. такое, что , называется корнем (нулем) функции или корнем уравнения
f(x)=0. (2.1)
Приближенное вычисление корня, как правило, распадается на две задачи:
1 отделение корней, т.е. определение интервалов, в каждом из которых содержится только один корень уравнения;
2 уточнение корня, т.е. вычисление его с заданной степенью точности.
При отделении корней уравнения общего вида (2.1) часто используется известная из курса математического анализа теорема Больцано - Коши:
пусть функция f(x) непрерывна на отрезке и на концах отрезка принимает значения разных знаков, т.е.
. Тогда существует такая точка , принадлежащая интервалу , в которой функция обращается в нуль. Заметим, что корень будет единственным, если (или ) существует и сохраняет знак на рассматриваемом отрезке.
Остановимся более подробно на алгебраических уравнениях
. (2.2)
Верхнюю границу модулей корней уравнения (2.2) дает следующая теорема.
Пусть . Тогда любой корень уравнения (2.2) удовлетворяет условию
. (2.3)
Допустим, что существует корень уравнения (2.2), не удовлетворяющий условию (2.3), т.е.
. (2.4)
Из (2.4) следует, что
.
Тогда
Согласно (2.4)
и ,
что противоречит предположению о том, что - корень уравнения (2.2).
2.2. Метод половинного деления
Пусть функция f(x) непрерывна на отрезке [a,b] и имеет на его концах разные по знаку значения. Задача состоит в том, чтобы вычислить корень уравнения (2.1.), принадлежащий отрезку [a,b] с заданной степенью точности , т.е. найти такое приближенное значение корня Xn, ( п - номер итерации), что
. (2.5)
В методе половинного деления за приближенное значение корня принимается середина отрезка
.
При этом очевидно, что
.
Затем определяется знак и для дальнейшего деления пополам выбирается тот из двух отрезков , на концах которого функция f(x) имеет разные по знаку значения. Расчет продолжается до тех пор, пока не выполнится условие (2.5) либо условие
. (2.6)