
- •И. Н. Двойцова вычислительная математика
- •Введение
- •1. Элементарная теория погрешностей; вычислительные задачи, методы и алгоритмы
- •1.1. Источники и классификация погрешностей результата численного эксперимента
- •1.2. Погрешности чисел
- •1.3. Погрешности арифметических операций
- •1.4. Погрешности функций
- •1.5. Особенности машинной арифметики
- •1.6. Лабораторная работа № 1. Определение абсолютной и относительной погрешностей приближенных чисел. Оценка погрешностей результата
- •1.7. Корректность вычислительной задачи
- •1.8. Обусловленность вычислительной задачи
- •1.9. Вычислительные методы, их классификация
- •2. Приближение функций
- •2.1. Задача приближения функций
- •2.2. Интерполяция обобщенными многочленами
- •2.3. Полиномиальная интерполяция. Многочлен Лагранжа
- •2.4. Погрешность интерполяции
- •2.5. Конечные разности и их свойства
- •Доказательство
- •2.6. Разделенные разности и их свойства
- •2.7. Интерполяционный многочлен Ньютона
- •2.8. Вычислительная схема Эйткена
- •2.9. Лабораторная работа № 2. Интерполирование и экстраполирование данных. Интерполяционный многочлен Лагранжа
- •2.10. Интерполяционный многочлен Ньютона с конечными разностями
- •2.11. Лабораторная работа № 3. Интерполирование и экстраполирование данных. Интерполяционный многочлен Ньютона
- •3. Метод наименьших квадратов и специальные интерполяционные многочлены
- •3.1. Постановка задачи и вывод формул метода наименьших квадратов
- •3.2. Лабораторная работа № 4. Аппроксимация функции по методу наименьших квадратов
- •3.3. Глобальная полиномиальная интерполяция
- •3.4. Чувствительность интерполяционного многочлена к погрешностям входных данных
- •4. Численное дифференцирование и численное интегрирование
- •4.1. Простейшие формулы численного дифференцирования для первой производной
- •4.2. Формулы численного дифференцирования для второй производной
- •4.3. Формулы численного дифференцирования, основанные на интерполяции алгебраическими многочленами
- •4.4. Обусловленность формул численного дифференцирования
- •4.5. Простейшие квадратурные методы численного интегрирования
- •4.6. Оценка погрешностей простейших квадратурных формул
- •4.7. Квадратурные формулы интерполяционного типа
- •4.8. Квадратурные формулы Гаусса
- •4.9. Лабораторная работа № 8. Численное дифференцирование и численное интегрирование функций
- •Ирина Николаевна Двойцова Вычислительная математика
- •660049, Красноярск, пр. Мира , 82
1.3. Погрешности арифметических операций
Теорема
1.4. Абсолютная погрешность
алгебраической суммы не превосходит
суммы абсолютных погрешностей слагаемых,
то есть
Доказательство
Запишем формулу для абсолютной погрешности алгебраической суммы двух величин по определению и воспользуемся свойствами модуля. Получим сразу необходимый результат
.
(1.3.1)
Теорема
1.5. Пусть
- ненулевые числа одного знака. Тогда
(1.3.2)
Доказательство
Поступим так же, как и в предыдущей теореме. Выразим абсолютную погрешность через относительную:
(1.3.3)
Из
последних двух равенств видно, что при
вычитании двух близких чисел может
произойти катастрофическая потеря
точности, так как при
в последней формуле
Теорема 1.6. Для относительных погрешностей произведения и частного приближенных чисел верны следующие оценки:
(1.3.4)
Доказательство
Первая требуемая формула получается традиционным путем:
Для
вывода второй оценки предварительно
получим, используя свойства модуля,
следующую формулу:
Тогда
Следствие.
Если
,
то для оценки границ относительных
погрешностей можно использовать
следующие приближенные равенства
Чаще всего на практике делают именно
так.
1.4. Погрешности функций
Пусть
- функция
-
переменных, дифференцируемая в
рассматриваемой области (например, на
отрезке
).
Теорема
1.7. Для абсолютной погрешности
значения
справедлива следующая формула:
. (1.4.1)
Доказательство
Вспомним
сначала формулу Тейлора
для функции нескольких переменных. Для
функции одного переменного разложение
в окрестности точки
будет иметь вид
Для функции переменных форма записи формулы Тейлора остается точно такой же, если вместо производных записать дифференциалы соответствующих порядков:
где,
.
Например, для функции двух переменных
Отбрасывая все члены второго порядка и выше, получим
Таким
образом, искомая формула сразу вытекает
из формулы Лагранжа1.
Если
достаточно мало, то для предельных
значений погрешностей можно положить
Для относительных погрешностей тогда имеем следующие формулы:
(1.4.2)
I
частный случай. Функция
- функция одного переменного. Здесь
следует положить
,
тогда
.
Для относительных погрешностей все
аналогично:
. (1.4.3)
II
частный случай. Функция
- неявная. Этот случай отличается от
исходного только формулой для нахождения
частных производных:
Пример.
Пусть корни квадратного уравнения
вычисляются при коэффициентах
Каково влияние погрешностей задания
коэффициентов на точность вычисляемых
значений?
Решение
Договоримся
сразу, чтобы избежать путаницы, об
обозначениях. Здесь роль двух функций
(
в формулах 1.4.1 и 1.4.2) играют
-
корни квадратного уравнения. Роль
переменных
играют
- коэффициенты квадратного уравнения.
Их всего два. В этих обозначениях
.
Однако, этих явных уравнений нет, они
лишь предполагаются, задана же неявная
функция
.
Найдем
выражения для частных производных
.
Отсюда
.
Тогда по формуле (1.4.2) получим
,
где
Здесь первый индекс у коэффициента
- номер переменного, а второй индекс
соответствует номеру функции, то есть
номеру корня квадратного уравнения.
Решим исходное квадратное уравнение.
Тогда
.
Аналогично,
Итак,