
- •Глава 1. Интерполирование
- •Постановка задачи интерполирования
- •1.2. Единственность интерполяционного многочлена
- •. Интерполяционный многочлен Лагранжа
- •1.4. Погрешность интерполяционной формулы Лагранжа
- •. Конечные разности и их свойства
- •1.6. Первая интерполяционная формула Ньютона
- •1.7. Вторая интерполяционная формула Ньютона
- •1.8. Интерполяционная формула Гаусса
- •1.10. Численное дифференцирование
- •Глава 2. Численное интегрирование
- •2.1. Общие замечания
- •Глава 3. Численные методы решения алгебраических и трансцендентных уравнений
- •3.1. Вводные замечания
- •3.2. Отделение корней
- •3.3. Метод половинного деления
- •3.5. Метод Ньютона (метод касательных)
- •3.6. Комбинированный метод
- •3.7. Метод итерации (метод последовательных приближений)
- •Глава 4. Решение систем линейных уравнений
- •4.1. Метод Гаусса (схема единственного деления)
- •4.2. Вычисление определителей
- •4.3. Вычисление обратной матрицы
- •4.4. Некоторые сведения из линейной алгебры
- •Рассмотрим систему линейных алгебраических уравнений
- •Вычитая второе равенство из первого, получаем
- •4.6. Приведение линейной системы к виду,
- •4.7. Метод Зейделя Пусть система линейных алгебраических уравнений приведена к виду
- •Глава 5. Обработка результатов наблюдений
- •5.1. Постановка задачи
- •5.2. Приближение функции, заданной таблично,
- •Глава 6. Численные методы решения
- •6.1. Методы решения задачи Коши. Вводные замечания
- •6.2. Решение дифференциальных уравнений
- •6.3. Метод Эйлера Пусть дано дифференциальное уравнение
- •6.4. Уточненный метод Эйлера Рассмотрим дифференциальное уравнение
- •6.6. Методы Рунге – Кутта Пусть дано дифференциальное уравнение
- •6.7. Метод Милна
- •6.8. Приближенное решение систем дифференциальных
- •6.9.Краевые задачи для обыкновенных дифференциальных
- •6.10. Решение краевых задач
- •6.11. Метод прогонки Рассмотрим линейную систему уравнений (6.54). Преобразуя первые (n-1) уравнения:
- •Глава 7. Численное решение дифференциальных уравнений в частных производных
- •7.1. Классификация дифференциальных уравнений
- •7.2. Метод сеток решения краевых задач
- •7.3. Погрешность аппроксимации дифференциальных уравнений
- •7.4. Аппроксимация граничных условий
- •7.6. Метод сеток для уравнений
. Конечные разности и их свойства
Пусть известны
значения некоторой функции y=f(x)
для
равноотстоящих значений аргумента
:
,
,
…,
.
Постоянная величина h называется шагом таблицы.
Конечными
разностями первого порядка
называются величины
,
,
…,
… . Аналогично вводятся конечные разности
второго порядка
,
,
…,
,
… и т.д. Конечные
разности
(m+1) порядка определяются через конечные разности m-го порядка следующим образом:
;
;
…;
.
Табл. 1.1 конечных разностей для функции y=f(x) обычно записывают в следующем виде,
Таблица 1.1.
х |
y |
|
|
|
|
x0 |
y0 |
|
|
|
|
|
|||||
x1 |
y1 |
|
|||
|
|
||||
x2 |
y2 |
|
|
||
|
|
||||
x3 |
y3 |
|
|
||
|
|
||||
x4 |
y4 |
|
|||
|
Приведенная таблица называется диагональной таблицей конечных разностей. Иногда используется горизонтальная таблица разностей (х0, y0, , , и т.д. располагаются в одной верхней строке).
Если табличные значения функции заданы с одинаковым числом десятичных знаков, то при оформлении таблицы конечных разностей разности записываются в единицах последенго разряда табличных значений функции без нулей впереди.
Так,
например, табл. 1.2 конечных разностей
для функции
будет
выглядеть следующим образом .
Таблица 1.2
х |
y |
|
|
|
|
2.0 |
0.6931 |
|
|
|
|
488 |
|||||
2.1 |
0.7419 |
-22 |
|||
466 |
0 |
||||
2.2 |
0.7885 |
-22 |
4 |
||
444 |
4 |
||||
2.3 |
0.8329 |
-18 |
|
||
426 |
|
||||
2.4 |
0.8755 |
|
|||
|
Свойства конечных разностей
Конечная разность суммы двух функций
равна сумме конечных разностей слагаемых
и .
Если
, где с – постоянный множитель, то
Свойства 1 и 2 верны для конечных разностей любого порядка и следуют из определения конечных разностей.
Конечная разность порядка n от многочлена степени n равна постоянной величине, и, следовательно, конечные разности более высокого порядка равны нулю.
Рассмотрим произвольный многочлен степени n:
.
Покажем, что конечная разность первого порядка представляет собой многочлен степени (n-1).
Действительно, имеем
.
Аналогично получаем
и т.д.:
,
.
Разность порядка k можно выразить непосредственно через табличные значения функции.
Покажем это. По определению конечных разностей имеем:
;
;
.
При помощи метода математической индукции можно показать, что
.
(1.14)
Запишем выражение (1.14) в компактной форме. Для этого введем оператор Е, применение которого к функции f(x) увеличивает ее аргумент на величину h:
;
;
;
.
Формулу (1.14) можно записать в следующем виде:
.
Это символическое
равенство следует понимать так. Выражение
записывается по формуле бинома Ньютона
и полученный многочлен умножается на
ym.
Табличные значения функции могут быть выражены через конечные разности различных порядков.
По определению
,
отсюда
.
Аналогично
,
отсюда
.
При помощи метода математической индукции можно показать, что
.
(1.15)
Эту формулу можно записать компактно в символической форме:
.