- •Конспект лекций
- •Пермь 1995
- •Лекция №1 Введение. Понятие о численных методах. История развития численных методов.
- •Приближенные числа. Абсолютная и относительная погрешности.
- •Погрешность числа (а).
- •Абсолютная погрешность().
- •Предельная абсолютная погрешность.
- •Относительная погрешность ().
- •Предельная относительная погрешность (a).
- •Основные источники погрешности.
- •Значащие и верные цифры.
- •Верные цифры.
- •Округление чисел.
- •Связь относительной погрешности приближённого числа с количеством верных цифр этого числа.
- •Лекция № 2 Погрешность суммы.
- •Погрешность разности.
- •Погрешность произведения.
- •Число верных знаков.
- •Погрешность частного.
- •Общая формула для погрешности.
- •Обратная задача теории погрешности.
- •Лекция №3 и 4 Интерполяция функций.
- •Постановка задачи.
- •Конечные разности различных порядков.
- •Первая интерполяционная формула Ньютона.
- •Вторая интерполяционная формула Ньютона.
- •Общая характеристика интерполяционных формул с постоянным шагом.
- •Интерполяционная формула Лагранжа.
- •Частные случаи.
- •Лекции №5 Фрмула Ньютона для неравностоящих узлов Разделённые разности
- •Интерполяционная формула Ньютона для неравностоящих значений аргумента
- •Погрешность формулы Ньютона для неравностоящих узлов
- •Интерполяция сплайками
- •Многочлены Чебышева
- •Выбор узлов интерполирования
- •Обратное интерполирование для равноотстоящих узлов
- •Обратное интерполирование для неравноотстоящих точек
- •Общие выводы по задаче интерполяции
- •Оценка погрешности интерполяционной формулы Лагранжа
- •Лекция n 6 и 7 приближенное решение алгебраических и трансцендентных уравнений.
- •Определение корней Графический способ.
- •Аналитический способ.
- •Графическое решение уравнений.
- •Метод половинного деления.
- •Алгоритм метода
- •Метод пропорциональных частей (метод хорд)
- •Алгоритм метода
- •Метод ньютона (метод касательных)
- •Сходимость метода итераций
- •Сходимость метода
- •Лекция №8 и 9 Решение системы линейных уравнений Общая характеристика методов решения систем линейных уравнений
- •Метод Гаусса
- •Трудоёмкость метода Гаусса
- •Достоинства метода
- •Метод итераций
- •Сходимость метода итераций для решения системы алгебраических уравнений
- •Достоинства метода итераций
- •Метод Зейделя
- •Лекция №10 Численное решение систем линейных уравнений
- •Метод Ньютона
- •Сходимость метода Ньютона
- •Теорема о существовании корней и сходимости процесса Ньютона
- •Метод скорейшего спуска (градиентный метод)
- •Градиент функции u
- •Сходимость градиентного метода
- •Лекции №11 и 12 Приближённое дифференцирование.
- •Графический способ дифференцирования.
- •Алгоритм построения графика производной.
- •Лекция n 13
- •Особенности метода Эйлера.
- •Метод Руте-Кутта.
- •Метод Адамса.
- •Лекция №14 Краевая задача. Методы её решения.
- •Метод конечных разностей для линейных дифференциальных уравнений.
- •Метод прогонки.
- •Метод конечных разностей для нелинейных уравнений второго порядка.
- •Лекция №15 методы обработки эксперементальных данных. Постановка задачи
- •Узловые точки
- •Класс функций
- •Критерий согласия
- •Среднеквадратический критерий
- •Минимальный критерий или критерий чебышева
- •Вероятностно-зональный критерий
- •Точность
- •Метод наименьших квадратов постановка задачи
Трудоёмкость метода Гаусса
Прямой ход
Обратный ход
Общее число выполняемых арифметических действий
то есть для
Предложенный метод Гаусса ориентирован на то, чтобы ведущие элементы не равнялись 0. А если на каком-то шаге возникает ситуация, что ведущий элемент равен 0, то тогда схема “формально” непригодна, хотя заданная система может иметь единственное решение.
Тогда применяют разновидность метод Гаусса. Схема с выбором главного элемента:
Выбираем элемент
- наибольший по модулю и неявляющийся
свободным членом.Вычисляем коэффициенты
,
для всех
-тая
строка называется главной строкой.
Из каждой неглавной строки вычитаем главную строку, умноженную на
.
В результате получим матрицу, у которой
в
-ом
столбце все коэффициенты нулевые.Преобразуем матрицу следующим образом: отбрасываем - (главную) строку и -й столбец. Получим матрицу
.Делаем подобные преобразования над матрицей до тех пор, пока не получим одну строку из двух столбцов, которая является главной.
Для определения . Объединим все главные строки, начиная с последней. После надлежащего изменения неизвестных получается система с треугольной матрицой.
При работе на ЭВМ при вывод главного элемента может оказаться достаточно трудоёмкой задачей. Поэтому практически в качестве главной строки берут первую строку, а в качестве главного элемента - наибольший по модулю элемент этой строки.
Пример:
Схема Гаусса с выбором главного элемента
|
|
|
|
|
|
|
|
|
|
1 |
-0,6 |
3 |
1 |
-1 |
2 |
6 |
11 |
I |
2 |
|
5 |
1 |
3 |
-4 |
-12 |
-17 |
|
3 |
-0,4 |
2 |
0 |
1 |
-1 |
1 |
3 |
|
4 |
-0,2 |
1 |
-5 |
3 |
-3 |
3 |
-1 |
|
1 |
-0,333 |
|
1,6 |
0,8 |
-0,4 |
-1,2 |
0,8 |
II |
2 |
-0,083 |
|
0,4 |
2,2 |
-2,6 |
-3,8 |
-3,8 |
|
3 |
|
|
-4,8 |
3,6 |
-3,8 |
0,6 |
-4,4 |
III |
1 |
0,571 |
|
|
2,0 |
-1,665 |
-1,0 |
-0,665 |
|
2 |
|
|
|
2,5 |
-2,915 |
-3,75 |
-4,165 |
IV |
1 |
|
|
|
0,572 |
|
1,141 |
1,713 |
V |
|
|
|
|
1 |
|
2,0 |
3 |
VI |
|
|
|
|
|
1 |
3,0 |
4 |
VII |
|
|
|
1 |
|
|
-1,0 |
0 |
VIII |
|
|
1 |
|
|
|
1,0 |
2 |
