
- •1. Погрешности и их классификация
- •2. Приближенные числа и действия с ними.
- •Правила записи приближенных чисел.
- •Округление чисел
- •3. Уменьшение погрешностей
- •4. Устойчивость. Корректность. Сходимость.
- •5. Вычисление значений функций. Вычисление значений многочлена. Схема Горнера.
- •Вычисление значений рациональных дробей
- •Приближенное нахождение сумм числовых рядов.
- •Вычисление значений аналитической функции
- •Вычисление значений показательной функции
- •Вычисление значений логарифмической функции
- •Вычисление значений синуса и косинуса.
- •Вычисление значений гиперболического синуса
- •Вычисление значений гиперболического косинуса.
- •Применение метода итерации для приближенного вычисления значений функции
- •Вычисление квадратного корня
- •6. Приближение функций При решении задач постоянно встречается необходимость замены одной функции некоторой другой функцией .
- •6.1.2. Равномерное приближение.
- •6.2. Многочлены Тейлора.
- •6.3. Интерполяционный многочлен Лагранжа.
- •6.4. Линейная интерполяция.
- •7. Минимизация погрешности интерполяции. Многочлены Чебышева.
- •8. Интерполяция с равноотстоящими узлами.
- •9. Конечные и разделенные разности.
- •10.Интерполяционный многочлен Ньютона.
- •11. Интерполирование сплайнами.
- •12. Численное дифференцирование.
- •Погрешность численного дифференцирования.
- •13. Численное интегрирование.
- •Для удобства в вычислениях вместо дробных индексов вводят целочисленные, т.Е.
- •И формула Симпсона принимает вид
- •Еще одной формулой численного интегрирования является следующая, которая называется формулой Ньютона:
- •Остаточный член имеет вид
- •В последней формуле число узлов обязательно равно .
- •Особые случаи численного интегрирования.
- •Применение метода Монте-Карло к вычислению определенных и кратных интегралов.)
- •Приближенное решение алгебраических и трансцендентных уравнений. Отделение корней.
- •Теорема об оценке погрешности приближенного корня.
- •Графическое решение уравнений.
- •Метод половинного деления (дихотомия).
- •Метод хорд (метод пропорциональных частей)
- •Метод Ньютона (метод касательных).
- •Видоизмененный (упрощенный) метод Ньютона.
- •Комбинированный метод (хорд и касательных).
- •Метод итераций
- •Оценка погрешности приближения
- •15. Метод итерации для системы двух уравнений
- •Метод Гаусса.
- •Применение метода Гаусса для вычисления определителей.
- •Метод итерации для решения систем уравнений
- •Решение обыкновенных дифференциальных уравнений.
- •Так, например, уравнение
- •Является уравнением 1-го порядка.
- •В ходе поиска общего решения часто приходят к соотношению вида
- •Запишем это разложение в виде
- •Метод Эйлера с пересчетом.
- •Решение это приводим:
- •Видно, что при большом числе узлов метод Эйлера может привести к заметным погрешностям, в таких случаях предпочитают пользоваться численными методами высших порядков.
- •Метод Адамса.
2. Приближенные числа и действия с ними.
Пусть а – точное числовое значение некоторой величины; а - приближенное числовое значение этой величины. Величина
(2.1)
называется абсолютной погрешностью приближенного числа а, а величина
(2.2)
называется его относительной погрешностью.
Любое
число
,
о котором известно, что
(2.3)
называется предельной абсолютной погрешностью приближенного числа а.
Любое
число
,
о котором известно, что
(2.3')
называется предельной относительной погрешностью приближенного числа а.
По
определению принято, что
и
связаны отношением
.
Учет погрешностей в арифметических операциях.
Принимается за очевидное, что если с = а + b и с = а + b или с = а - b из с = а – b, то
(2.4)
и, следовательно, в качестве предельной абсолютной погрешности естественно взять
(2.5).
Таким образом, при сложении и вычитании двух приближенных чисел их предельные абсолютные погрешности складываются.
Рассмотрим пример.
Пусть а* = 100 5; b* = 50 5;
;
.
Здесь брались граничные случаи. Понятно, что если будем брать значения а* или b* внутри интервалов 100 ± 5 и 50 ± 5 соответственно, сумма или разность будут попадать в пределы максимального интервала с*.
Подобные
рассуждения можно проводить и для
произведения двух приближенных чисел
а
b,
их частного
и, в общем случае, для любой дифференцируемой
функции от этих приближенных чисел.
Получающиеся при этом формулы очень
напоминают уже известные нам формулы
математического анализа. Вспомним
формулу для полного дифференциала
функции z
двух переменных x
и y
.
(2.6).
В действиях с приближенными числами пользуются такой же формулой, только частные производные берут по модулю.
Выведем формулу (2.5), воспользовавшись формулой, аналогичной (2.6).
Пусть х = а*; у = b* и c* = a* b*. Тогда
Аналогично выведем формулы для погрешностей произведения и частного приближенных чисел а и b, используя формулу 2.6. Пусть u* = (a* b*). Тогда
(2.7).
Пусть
теперь
.
Тогда
и
(2.8).
Итак:
т.е. при сложении и вычитании приближенных чисел складывают предельные абсолютные погрешности, а при умножении и делении приближенных чисел складывают их предельные относительные погрешности.
По аналогии с приведенными примерами могут быть вычислены погрешности приближенных величин, являющихся функциями произвольного количества приближенных чисел. При этом пользуются формулой
(2.9),
где
.
Обратная
задача. Часто
приходится решать такую задачу: с какой
точностью надо задать значения аргументов
функции z
= z
(a1,
,
an),
чтобы погрешность
z
(
),
не превосходила заданную величину
?
Пусть
точки (a1,
,
an)
и (),
соответствующие истинным и приближенным
значениям параметров aj
(j
= 1,
2,
,
n),
принадлежат некоторой выпуклой области
G
и
.
Тогда
(см.
2.9)
Определение. Область G элементов а1, а2, , ап называется выпуклой, если прямая, соединяющая любые две точки этой области, нигде не пересекает границ этой области.
-
максимальное значение модуля частной
производной по аргументу aj
в области определения G
функции z
(a1,
a2,
,
an).
Очевидно, что любая совокупность
абсолютных погрешностей, удовлетворяющих
неравенству
(2.10)
будет обеспечивать требуемую точность.
Если
функция z
зависит только от одного аргумента (п
= 1),
то имеем неравенство
и для достижения требуемой точности
достаточно взять
.
В
случае п
> 1
иногда рекомендуют отвести для погрешности
каждого аргумента равную долю, то есть
выбрать
из условия
,
т.е.
.
В других случаях предлагают взять
погрешности равными и максимально
возможными, т.е. положить
,
где
.
Но это возможно в простейших случаях. Более сложные случаи мы пока рассматривать не будем.