- •Конспект лекций
- •Пермь 1995
- •Лекция №1 Введение. Понятие о численных методах. История развития численных методов.
- •Приближенные числа. Абсолютная и относительная погрешности.
- •Погрешность числа (а).
- •Абсолютная погрешность().
- •Предельная абсолютная погрешность.
- •Относительная погрешность ().
- •Предельная относительная погрешность (a).
- •Основные источники погрешности.
- •Значащие и верные цифры.
- •Верные цифры.
- •Округление чисел.
- •Связь относительной погрешности приближённого числа с количеством верных цифр этого числа.
- •Лекция № 2 Погрешность суммы.
- •Погрешность разности.
- •Погрешность произведения.
- •Число верных знаков.
- •Погрешность частного.
- •Общая формула для погрешности.
- •Обратная задача теории погрешности.
- •Лекция №3 и 4 Интерполяция функций.
- •Постановка задачи.
- •Конечные разности различных порядков.
- •Первая интерполяционная формула Ньютона.
- •Вторая интерполяционная формула Ньютона.
- •Общая характеристика интерполяционных формул с постоянным шагом.
- •Интерполяционная формула Лагранжа.
- •Частные случаи.
- •Лекции №5 Фрмула Ньютона для неравностоящих узлов Разделённые разности
- •Интерполяционная формула Ньютона для неравностоящих значений аргумента
- •Погрешность формулы Ньютона для неравностоящих узлов
- •Интерполяция сплайками
- •Многочлены Чебышева
- •Выбор узлов интерполирования
- •Обратное интерполирование для равноотстоящих узлов
- •Обратное интерполирование для неравноотстоящих точек
- •Общие выводы по задаче интерполяции
- •Оценка погрешности интерполяционной формулы Лагранжа
- •Лекция n 6 и 7 приближенное решение алгебраических и трансцендентных уравнений.
- •Определение корней Графический способ.
- •Аналитический способ.
- •Графическое решение уравнений.
- •Метод половинного деления.
- •Алгоритм метода
- •Метод пропорциональных частей (метод хорд)
- •Алгоритм метода
- •Метод ньютона (метод касательных)
- •Сходимость метода итераций
- •Сходимость метода
- •Лекция №8 и 9 Решение системы линейных уравнений Общая характеристика методов решения систем линейных уравнений
- •Метод Гаусса
- •Трудоёмкость метода Гаусса
- •Достоинства метода
- •Метод итераций
- •Сходимость метода итераций для решения системы алгебраических уравнений
- •Достоинства метода итераций
- •Метод Зейделя
- •Лекция №10 Численное решение систем линейных уравнений
- •Метод Ньютона
- •Сходимость метода Ньютона
- •Теорема о существовании корней и сходимости процесса Ньютона
- •Метод скорейшего спуска (градиентный метод)
- •Градиент функции u
- •Сходимость градиентного метода
- •Лекции №11 и 12 Приближённое дифференцирование.
- •Графический способ дифференцирования.
- •Алгоритм построения графика производной.
- •Лекция n 13
- •Особенности метода Эйлера.
- •Метод Руте-Кутта.
- •Метод Адамса.
- •Лекция №14 Краевая задача. Методы её решения.
- •Метод конечных разностей для линейных дифференциальных уравнений.
- •Метод прогонки.
- •Метод конечных разностей для нелинейных уравнений второго порядка.
- •Лекция №15 методы обработки эксперементальных данных. Постановка задачи
- •Узловые точки
- •Класс функций
- •Критерий согласия
- •Среднеквадратический критерий
- •Минимальный критерий или критерий чебышева
- •Вероятностно-зональный критерий
- •Точность
- •Метод наименьших квадратов постановка задачи
Достоинства метода
Если матрица вырождения, то перед исключением неизвестной главный элемент считается равным нулю =>
С помощью метода Гаусса можно вычислить определитель
треугольной матрицы.
При большом числе неизвестных схема метода Гаусса, дающая точное решение, становится весьма сложной.
В этих случаях для нахождения корней системы лучше пользоваться приближёнными численными методами.
Метод итераций
Дана система уравнений
Можно привести
систему к такому виду, чтобы диагональные
элементы были отличные от нуля, то есть
,
тогда разрешая
-тое
уравнение относительно
,
получаем
,
(2)
где
,
,
или
при
,
.
Тогда систему уравнений 2 можно записать в виде:
- итерационная
формула.
Таким образом,
выбрав начальные значения
и так далее.
Итерации
останавливаются, когда
,
Сходимость метода итераций для решения системы алгебраических уравнений
Теорема:
Система уравнений имеет единственное
решение и сходится при любом начальном
значении
тогда и только тогда, когда все собственные
значения матрицы
по модулю меньше 1.
Если для системы уравнений
выполнено хотя бы одно из условий:
,
,
то процесс итерации сходится, независимо от выбора начального условия.
Однако этой теоремой в общем случае очень тяжело воспользоваться, поэтому на практике пользуются другим правилом менее жёстким.
Если эти условия выполняются, то в принципе логично выбрать для начальных значений. На практике в качестве начального приближения используют вектор свободных членов.
Приведение линейной системы к виду, удобному для итерациию.
Теорема сходимости накладывает жёсткие условия к коэффициентам данной линейной системы.
Однако, если
,
то эту систему всегда можно привести к
такому виду:
,
чтобы удовлетворить условиям 1
Первый способ.
Дано:
Домножим это
уравнение на матрицу
,
где
,
где
Второй способ.
Каждое
-ое
уравнение делится на
Тогда
,
,
.
Тогда уравнение сходимости имеет вид
,
,
Эти неравенствабудут выполняться, если диагональные элементы будут удовлетворять условиям:
,
,
то есть если модули диагональных коэффициентов для которого уравнения системы больше суммы модулей всех остальных коэффициентов.
Достоинства метода итераций
Если итерации сходятся быстро, то есть для сходимости требуется менее итераций, то выигрыш во времени по сравнению с методом Гаусса:
,
-
число итераций
Погрешности округления в методе итераций сказывается значительно меньше, чем в методе Гаусса. Кроме того, метод итерации является самоисправляющимся, то есть отдельная ошибка запрещается в вычислениях, не отражаясь на конечном результате, то есть ошибочное приближение можно рассматривать как новый начальный вектор.
Метод итераций становится особенно выгодным при решении систем, у которых значительное число коэффициентов равно нулю.
Метод итераций легко программируется.
Метод Зейделя
Является модификацией
метода итераций. Основная идея заключается
в том, что при вычислении
-го
приближения
-го
корня используются уже вычисленные
приближённые корни
.
Дано:
,
Выбираем начальное приближение:
На
-том
шаге, согласно Зейделю строим приближение
по следующим формулам:
Метод Зейделя даёт полную сходимость по сравнению с методом итерации, но приводящий к громоздким вычислениям.
Теорема: Для существования единственного решения системы сходимости метода Зейделя достаточно выполнение хотя бы одного из двух условий:
,
матрица - симметричная положительно определённая (все её соответственно значения положительны)
