 
        
        - •А.П. Бырдин н.В. Заварзин а.А. Сидоренко л.П. Цуканова основы численных методов Учебное пособие
- •А.П. Бырдин н.В. Заварзин а.А. Сидоренко л.П. Цуканова
- •Введение
- •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.2)
Если  непосредственно вычислить  значение
этого многочлена в произвольной точке
 ,
то потребуется
,
то потребуется 
 арифметических действий. При этом
возникает потеря точности за счет
погрешности вычисления.
арифметических действий. При этом
возникает потеря точности за счет
погрешности вычисления.
Вычисление многочлена 
 удобнее производить следующим образом.
Представим его в виде
удобнее производить следующим образом.
Представим его в виде
 
 .
  (2.3)
.
  (2.3)
Отсюда, последовательно вычисляя числа
 (2.4)
                                     (2.4)
находим 
 .
.
Описанный метод вычисления значения
многочлена в заданной точке носит
название схемы Горнера. Она позволяет
сократить число сложений и умножений 
до 
 операций, при этом повышается точность
вычисления.
операций, при этом повышается точность
вычисления.
Пример. Вычислить значение многочлена
 при
  при   
 .
.
Решение. Имеем  
 
 
 
 
 
 
 Последовательное применение формулы
(2.4)  дает
Последовательное применение формулы
(2.4)  дает 
 
Следовательно,  
 
Практически схему Горнера удобно реализовывать в виде следующей таблицы
 
Тогда рассмотренный пример будет иметь вид
 
2.3. Интерполяционный многочлен Лагранжа
Будем искать многочлен 
 степени не выше 
,
имеющий в заданных узлах интерполяции
те же значения, что и функция 
.
Запишем этот  многочлен в виде  (2.2).
Неизвестные коэффициенты
степени не выше 
,
имеющий в заданных узлах интерполяции
те же значения, что и функция 
.
Запишем этот  многочлен в виде  (2.2).
Неизвестные коэффициенты 
 ,
,
  
 ,
 …,
,
 …, 
 можно найти из условия
 можно найти из условия 
 .
Это условие  приводит к системе линейных
уравнений
.
Это условие  приводит к системе линейных
уравнений
                     
            
 (2.5)
                   (2.5)
Определитель этой системы
 
называется
определителем Ван-дер-Монда. Он отличен
от нуля, так как все узлы интерполирования
различны. Следовательно,  система (2.5) 
имеет единственное  решение и 
интерполяционный многочлен  
определяется единственным образом.
Однако, для большого  числа узлов
интерполяции,  решать эту систему 
сложно. Многочлен 
можно   построить 
    другим     способом:
  строятся      
 
многочленов
 
 таких,  что
таких,  что 
 
Так как искомый
многочлен обращается  в нуль в 
точках 
 ,
 то он имеет вид
,
 то он имеет вид
 ,
   (2.6)
,
   (2.6)
где 
 - постоянный коэффициент. Полагая
- постоянный коэффициент. Полагая  
 в формуле  (2.6)  и учитывая, что
в формуле  (2.6)  и учитывая, что 
 ,
получим
,
получим
 
Отсюда
 .
.
Подставив это значение в формулу (2.6), будем иметь
 .
  (2.7)
.
  (2.7)
Тогда, многочлен 
,
удовлетворяющий условиям 
 ,
 имеет вид
,
 имеет вид
 .				(2.8)
.				(2.8)
Подставив в
формулу (2.8) значение 
 из (2.7), получим
из (2.7), получим
 .
   (2.9)
.
   (2.9)
Многочлен , определяемый по формуле (2.9) , называется интерполяционным многочленом Лагранжа. Очевидно, что степень многочлена Лагранжа не превышает числа . Часто многочлен Лагранжа записывают виде
 .
.
При  
 имеем две точки, и формула Лагранжа
представляет в этом случае уравнение
прямой
имеем две точки, и формула Лагранжа
представляет в этом случае уравнение
прямой 
 ,
проходящей через две заданные точки
,
проходящей через две заданные точки 
 ,
,
где 
 - абсциссы  этих точек.
 - абсциссы  этих точек.
При  
получим уравнение параболы 
 ,
проходящей через три точки
,
проходящей через три точки
 .
.
Пример. Для функции 
 построить интерполяционный многочлен
Лагранжа, выбрав узлы
построить интерполяционный многочлен
Лагранжа, выбрав узлы
 
Решение. Многочлен Лагранжа для трех узлов интерполирования запишется так
 Вычисляем
соответствующие значения функции
Вычисляем
соответствующие значения функции
 
Применяя формулу Лагранжа, получим
 
  или   
 
