
- •Конспект лекций
- •Пермь 1995
- •Лекция №1 Введение. Понятие о численных методах. История развития численных методов.
- •Приближенные числа. Абсолютная и относительная погрешности.
- •Погрешность числа (а).
- •Абсолютная погрешность().
- •Предельная абсолютная погрешность.
- •Относительная погрешность ().
- •Предельная относительная погрешность (a).
- •Основные источники погрешности.
- •Значащие и верные цифры.
- •Верные цифры.
- •Округление чисел.
- •Связь относительной погрешности приближённого числа с количеством верных цифр этого числа.
- •Лекция № 2 Погрешность суммы.
- •Погрешность разности.
- •Погрешность произведения.
- •Число верных знаков.
- •Погрешность частного.
- •Общая формула для погрешности.
- •Обратная задача теории погрешности.
- •Лекция №3 и 4 Интерполяция функций.
- •Постановка задачи.
- •Конечные разности различных порядков.
- •Первая интерполяционная формула Ньютона.
- •Вторая интерполяционная формула Ньютона.
- •Общая характеристика интерполяционных формул с постоянным шагом.
- •Интерполяционная формула Лагранжа.
- •Частные случаи.
- •Лекции №5 Фрмула Ньютона для неравностоящих узлов Разделённые разности
- •Интерполяционная формула Ньютона для неравностоящих значений аргумента
- •Погрешность формулы Ньютона для неравностоящих узлов
- •Интерполяция сплайками
- •Многочлены Чебышева
- •Выбор узлов интерполирования
- •Обратное интерполирование для равноотстоящих узлов
- •Обратное интерполирование для неравноотстоящих точек
- •Общие выводы по задаче интерполяции
- •Оценка погрешности интерполяционной формулы Лагранжа
- •Лекция n 6 и 7 приближенное решение алгебраических и трансцендентных уравнений.
- •Определение корней Графический способ.
- •Аналитический способ.
- •Графическое решение уравнений.
- •Метод половинного деления.
- •Алгоритм метода
- •Метод пропорциональных частей (метод хорд)
- •Алгоритм метода
- •Метод ньютона (метод касательных)
- •Сходимость метода итераций
- •Сходимость метода
- •Лекция №8 и 9 Решение системы линейных уравнений Общая характеристика методов решения систем линейных уравнений
- •Метод Гаусса
- •Трудоёмкость метода Гаусса
- •Достоинства метода
- •Метод итераций
- •Сходимость метода итераций для решения системы алгебраических уравнений
- •Достоинства метода итераций
- •Метод Зейделя
- •Лекция №10 Численное решение систем линейных уравнений
- •Метод Ньютона
- •Сходимость метода Ньютона
- •Теорема о существовании корней и сходимости процесса Ньютона
- •Метод скорейшего спуска (градиентный метод)
- •Градиент функции u
- •Сходимость градиентного метода
- •Лекции №11 и 12 Приближённое дифференцирование.
- •Графический способ дифференцирования.
- •Алгоритм построения графика производной.
- •Лекция n 13
- •Особенности метода Эйлера.
- •Метод Руте-Кутта.
- •Метод Адамса.
- •Лекция №14 Краевая задача. Методы её решения.
- •Метод конечных разностей для линейных дифференциальных уравнений.
- •Метод прогонки.
- •Метод конечных разностей для нелинейных уравнений второго порядка.
- •Лекция №15 методы обработки эксперементальных данных. Постановка задачи
- •Узловые точки
- •Класс функций
- •Критерий согласия
- •Среднеквадратический критерий
- •Минимальный критерий или критерий чебышева
- •Вероятностно-зональный критерий
- •Точность
- •Метод наименьших квадратов постановка задачи
Конечные разности различных порядков.
Пусть: y =f(x) - заданная функция
-
фиксированная величина приращения
аргумента
Тогда
- называется первой конечной разностью
функции y,
или конечной
разностью первого порядка.
Вторая конечная разность, или конечная разность второго порядка.
Третья конечная разность, или конечная разность третьего порядка.
Т.о., в общем виде:
Конечная разность n-го порядка.
Пример:
Конечные разности различных порядков удобно располагать в форме таблиц двух видов: горизонтальной и диагональной таблиц разностей
x |
y |
|
|
|
x0 |
y0 |
|
|
|
x1 |
y1 |
|
|
|
... |
... |
... |
... |
... |
Диагональная таблица разностей.
x |
y |
|
|
|
X0 |
Y0 |
|
|
|
|
|
|
|
|
X1 |
Y1 |
|
|
|
|
|
|
|
|
X2 |
Y2 |
|
|
|
|
|
|
|
|
X3 |
Y3 |
|
|
|
Пример: горизонтальная таблица функции y = f(x) = x2 при
, x0 = 0 начальное значение, x6=5 конечное значение
x |
y |
|
|
|
0 |
0 |
1 |
2 |
0 |
1 |
1 |
3 |
2 |
0 |
2 |
4 |
5 |
2 |
0 |
3 |
9 |
7 |
2 |
|
4 |
16 |
9 |
|
|
5 |
25 |
|
|
|
Диагональная таблица
x |
y |
|
|
|
0 |
0 |
|
|
|
|
|
1 |
|
|
1 |
1 |
|
2 |
|
|
|
3 |
|
0 |
2 |
4 |
|
2 |
|
|
|
5 |
|
|
3 |
9 |
|
|
|
При составлении
таблиц возможные ошибки вычисляются и
диагональная таблица наглядно показывает
нам, как отразится ошибка
в значении yn.
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Следует заметить,
что максимальная ошибка
– в той же горизонтальной строке, где
и табличная величина yn.
Пример: исправить ошибку в таблице (неверные цифры взяты в скобки).
|
|
|
|
Ошибка |
15 |
13,260 |
|
0 |
|
|
|
884 |
|
|
16 |
14,144 |
|
0 |
|
|
|
884 |
|
|
17 |
15,912 |
|
0 |
|
|
|
884 |
|
|
18 |
15,028 |
|
(-4)0 |
|
|
|
88(0)4 |
|
|
19 |
16,79(2)6 |
|
(8)0 |
|
|
|
88(8)4 |
|
|
20 |
17,680 |
|
(-4)0 |
|
|
|
884 |
|
|
21 |
18,564 |
|
0 |
|
|
|
884 |
|
|
22 |
19,448 |
|
0 |
|
|
|
884 |
|
|
23 |
20,332 |
|
0 |
|
Как видно из таблицы, ход вторых разностей нарушается при x=19. Ошибка распространяется на 3 строки. Находим среднее арифметическое значение второй разности для средней из 3 точек:
=
,
=
Внося исправление в табличное значение y для x=19, получим верное значение функции:
n=(yn+ )- =16.792-(-0.004)=16.796.