- •А.П. Бырдин н.В. Заварзин а.А. Сидоренко л.П. Цуканова основы численных методов Учебное пособие
- •А.П. Бырдин н.В. Заварзин а.А. Сидоренко л.П. Цуканова
- •Введение
- •1. Действия над приближенными числами
- •1.1. Основные источники погрешностей
- •1.2. Приближенные числа. Абсолютная и относительная погрешности
- •1.3. Правила записи приближенных чисел
- •Решение. В нашем случае и . Следовательно,
- •Решение. Имеем .
- •2. Интерполирование функции
- •2.1. Постановка задачи интерполирования
- •2.2. Вычисление значений многочлена по схеме Горнера
- •2.3. Интерполяционный многочлен Лагранжа
- •Интерполяционный многочлен Ньютона
- •Сопоставление интерполяционных формул Лагранжа и Ньютона. Погрешность интерполяции
- •2.6. Интерполирование функции кубическими сплайнами
- •Методы численного решения систем
- •3.1. Метод Гаусса
- •3.2. Метод итерации
- •3.3. Метод Зейделя
- •4. Методы численного решения
- •4.1. Отделение корней
- •Метод половинного деления
- •4.3. Метод хорд
- •4.4. Метод Ньютона
- •Комбинированный метод
- •4.6. Метод итерации
- •4.7. Метод Ньютона для системы двух уравнений
- •Метод итерации для системы двух уравнений
- •5. Численное дифференцирование
- •5.1. Постановка вопроса
- •5.2. Формулы приближенного дифференцирования, основанные на первой интерполяционной формуле Ньютона
- •5.3. Конечно-разностные аппроксимации производных
- •6. Среднеквадратичное приближение функций
- •6.1. Метод наименьших квадратов
- •7. Численное интегрирование
- •Численное решение обыкновенных дифференциальных уравнений
- •8.1. Понятие о численном решении задачи Коши
- •8.2. Метод Эйлера
- •8.3. Методы Рунге-Кутта
- •8.4. Численные решение систем дифференциальных уравнений первого порядка
- •9. Численное решение дифференциальных уравнений в частных производных
- •Решение задачи Дирихле для уравнения Лапласа методом сеток
- •Решение смешанной задачи для уравнения гиперболического типа методом сеток
- •9.3. Решение смешанной задачи для уравнения теплопроводности типа методом
- •Библиографический список
- •Оглавление
- •1. Действия над приближенными числами ..……………...4
- •2. Интерполирование функций ……………………….…... 9
- •3. Методы численного решения систем линейных
- •4. Методы численного решения нелинейных уравнений
- •10. Библиографический список …….…….……..……… 110
- •Составители: Бырдин Аркадий Петрович
- •394026 Воронеж, Московский просп., 14
2.2. Вычисление значений многочлена по схеме Горнера
Поскольку большинство методов аппроксимации функций сводится к замене этих функций многочленами, то одной из важнейших задач является отыскание удобного способа вычисления значения многочлена в произвольной точке. Одним из критериев удобства является уменьшение числа сложений и умножений, необходимых для вычисления многочлена.
Рассмотрим многочлен ой степени
. (2.2)
Если непосредственно вычислить значение этого многочлена в произвольной точке , то потребуется арифметических действий. При этом возникает потеря точности за счет погрешности вычисления.
Вычисление многочлена удобнее производить следующим образом. Представим его в виде
. (2.3)
Отсюда, последовательно вычисляя числа
(2.4)
находим .
Описанный метод вычисления значения многочлена в заданной точке носит название схемы Горнера. Она позволяет сократить число сложений и умножений до операций, при этом повышается точность вычисления.
Пример. Вычислить значение многочлена
при .
Решение. Имеем Последовательное применение формулы (2.4) дает
Следовательно,
Практически схему Горнера удобно реализовывать в виде следующей таблицы
Тогда рассмотренный пример будет иметь вид
2.3. Интерполяционный многочлен Лагранжа
Будем искать многочлен степени не выше , имеющий в заданных узлах интерполяции те же значения, что и функция . Запишем этот многочлен в виде (2.2). Неизвестные коэффициенты , , …, можно найти из условия . Это условие приводит к системе линейных уравнений
(2.5)
Определитель этой системы
называется определителем Ван-дер-Монда. Он отличен от нуля, так как все узлы интерполирования различны. Следовательно, система (2.5) имеет единственное решение и интерполяционный многочлен определяется единственным образом. Однако, для большого числа узлов интерполяции, решать эту систему сложно. Многочлен можно построить другим способом: строятся
многочленов таких, что
Так как искомый многочлен обращается в нуль в точках , то он имеет вид
, (2.6)
где - постоянный коэффициент. Полагая в формуле (2.6) и учитывая, что , получим
Отсюда
.
Подставив это значение в формулу (2.6), будем иметь
. (2.7)
Тогда, многочлен , удовлетворяющий условиям , имеет вид
. (2.8)
Подставив в формулу (2.8) значение из (2.7), получим
. (2.9)
Многочлен , определяемый по формуле (2.9) , называется интерполяционным многочленом Лагранжа. Очевидно, что степень многочлена Лагранжа не превышает числа . Часто многочлен Лагранжа записывают виде
.
При имеем две точки, и формула Лагранжа представляет в этом случае уравнение прямой , проходящей через две заданные точки
,
где - абсциссы этих точек.
При получим уравнение параболы , проходящей через три точки
.
Пример. Для функции построить интерполяционный многочлен Лагранжа, выбрав узлы
Решение. Многочлен Лагранжа для трех узлов интерполирования запишется так
Вычисляем соответствующие значения функции
Применяя формулу Лагранжа, получим
или