Вычисление квадратного корня
В качестве одного из многочисленных примеров рассмотрим нахождение точного значения квадратного корня.
Пусть .
Положим .
Тогда .
Применяя формулу (5.23),
,
имеем
,
или
,
, (5.25)
(процесс Герона).
Рис.5.3. К вычислению квадратного корня.
Последовательные приближения получаются по методу Ньютона, примененному к параболе
. Если за у0 принять табличное значение, дающее
с относительной погрешностью
, то у1, определенное по формуле (5.25), дает новое значение
приблизительно с относительной погрешностью
.
Действительно, полагая
и пренебрегая степенями d, выше третьей, будем иметь:
.
Отсюда получаем вывод: при применении процесса Герона число верных цифр примерно удваивается на каждом этапе по сравнению с первоначальным количеством.
Пример 1. Для приближенно имеем:
.
Уточняя это значение, получаем
.
Еще раз повторяя процесс, будем иметь:
, причем восемь или семь десятичных знаков являются верными. Действительно,
Приближение функций
При решении задач постоянно встречается необходимость замены одной функции некоторой другой функцией
.
Условнотакого рода задачи можно разделить на два типа:
1).Вид связи между параметрами и
не известен, но эта связь задана в виде таблицы
, т.е. дискретному множеству значений аргумента
поставлено в соответствие множество значений функции
.
Но при выполнении расчетов требуются и другие значения .
Эта цель достигается решением задачи о приближении (аппроксимации) функций: данную функцию требуется приближенно заменить (аппроксимировать) некоторой функцией
так, чтобы отклонение
от
в заданной области было наименьшим. Функция
при этом называется аппроксимирующей.
Для практики весьма важен случай аппроксимации ф-ций многочленами вида
(6.1)
В дальнейшем для аппроксимации будут рассматриваться лишь такие функции.
2). Вид связи известен. Например,
. Очевидно, что при ручном счете могут быть использованы таблицы, где с определенной погрешностью приведены значения
. Но при машинном счете ввод таблиц требует больших затрат памяти. Поэтому для вычисления значений функций на ЭВМ используются разложения этих функций в степенные ряды. Например, функция
вычисляется с помощью ряда
(6.2).
Если приближение строится на заданном дискретном множестве точек , то аппроксимация называется точечной.
При построении приближения на непрерывном множестве точек аппроксимация называется непрерывной.
6.1 Точечная аппроксимация.
Одним из основных типов точечной аппроксимации является интерполирование. Оно состоит в следующем: для данной функции строим многочлен (6.1), принимающий в заданных точках
те же значения
, что и функция
, т. е.
(6.3)
При этом предполагается, что среди значений нет одинаковых, т.е. xi ≠ xk при
. Точки
называются узлами интерполяции, а многочлен
-интерполяционным многочленом. Максимальная степень интерполяционного многочлена равна
.
В этом случае мы имеем дело с глобальной интерполяцией, поскольку один многочлен
используется для интерполяции функции на всем интервале изменения аргумента
. Коэффициенты
находятся из системы уравнений 6.3.
Если интерполяционные многочлены построить отдельно для разных частей рассматриваемого интервала изменения , то получим кусочную (или локальную) интерполяцию.
Если интерполяционные многочлены используются для приближенного вычисления функции вне рассматриваемого отрезка , то такое приближение называют экстраполяцией.
Кроме интерполирования, где требуется выполнение условий возможны и другие виды аппроксимации. Например, в случае глобальной интерполяции при большом количестве узлов интерполяции получается высокая степень многочлена (6.1). В этом случае можно пойти другим путем: выбирается многочлен меньшей степени, график которого проходит близко от данных точек (штриховая линия на рис.6.1).
![]() |
Рис. 6.1. К вопросу об экстраполяции.
6.1.1. Одним из таких видов является среднеквадратичное приближение функций с помощью многочлена. При этом ; случай
соответствует интерполяции.
На практике стараются подобрать многочлен как можно меньшей степени (как правило,
.
Мерой отклонения многочлена от заданной функции
на множестве точек
при среднеквадратичном приближении является величина
, равная сумме квадратов разностей между значениями многочлена и функции в данных точках:
(6.4)
Для построения аппроксимирующего многочлена нужно подобрать коэффициенты так, чтобы величина
была наименьшей. В этом состоит метод наименьших квадратов.