
- •2202 “Автоматизированные системы обработки
- •Приближенные числа.
- •Действия над приближенными числами.
- •Общая схема метода Гаусса для систем, имеющих единственное решение
- •Обратный ход:
- •Метод Гаусса с выбором максимального элемента по столбцу
- •Достаточное условие сходимости процесса итераций
- •Основные понятия алгебры матриц и теории линейных векторных пространств.
- •Обратная матрица
- •Правило нахождения ранга матрицы:
- •Нахождение обратной матрицы методом Гаусса
- •Применение метода итераций для умножения элементов обратной матрицы
- •Пример:
- •Приближенное решение алгебраических и трансцендентных уравнений
- •Пример:
- •Блок – схема решения уравнения f(X) методом половинного деления
- •Метод хорд
- •Метод Ньютона
- •Пример:
- •Приближенное решение алгебраических и трансцендентных уравнений (продолжение) Метод итераций
- •Теорема
- •Оценка приближения.
- •Интерполирование функций.
- •1. Конечные разности различных порядков.
- •2) Таблица разностей.
- •3) Постановка задачи интерполирования.
- •4) Первая интерполяционная формула Ньютона.
- •Блок-схема построения кубического сплайна
- •Численное дифференцирование
- •2.Конечно-разностные аппроксимации производных.
- •Блок-схема вычисления производной.
- •Численное интегрирование.
- •Блок-схема вычисления определенного интеграла методом двойного пересчета по формуле Симпсона.
- •Квадратурные формулы Гаусса
- •Приближенное вычисление несобственных интегралов с бесконечными пределами.
- •Среднеквадратичное приближение функций
- •Точечное квадратичное аппроксимирование функций
- •Условия нахождения экстремума функции
- •IV Интегральное квадратичное аппроксимирование функций на отрезке.
- •Ортогональные на промежутке системы функций
- •Основные понятия гармонического анализа.
- •4.2.1. Аппроксимация тригонометрическими полиномами кусочно-гладких функций
- •4.3. Среднеквадратичное приближение функций полиномами Лежандра.
- •4.3. Полиномы, наименее уклоняющиеся от нуля.
- •Эмпирические формулы
- •Численное решение обыкновенных дифференциальных уравнений.
- •Численное решение дифференциальных уравнений и систем дифференциальных уравнений высших порядков
- •Численные методы безусловной оптимизации.
Правило нахождения ранга матрицы:
-
Начиная с миноров первого порядка (элементов матрицы), переходить к минорам больших порядков.
-
Пусть найден минор D r-го порядка, отличный от нуля, тогда нужно вычислить лишь миноры (r+1)-го порядка, окаймляющие минор D. Если все эти миноры равны нулю, то ранг матрицы равен r. Если же хотя бы один отличен от нуля, то эту операцию нужно применять к нему, увеличить ранг матрицы А на 1.
Пример:
Найти ранг матрицы А (4х5)
В матрице содержатся миноры второго порядка, отличные от нуля, например:
Окаймляющий его минор третьего порядка:
Оба минора четвертого порядка, окаймляющие
минор
,
равны нулю.
Таким образом, r=3, дефект равен 1: m-r=1.
-
Клеточные матрицы
Разобьем исходную матрицу на блоки или клетки, или подматрицы
Клетки:
Тогда
Частные случаи клеточных матриц
-
Квазидиагональные матрицы
-
Окаймленные матрицы
Операции над клеточными матрицами выполняются по аналогии с обыкновенными матрицами (вместо элементов – клетки). При умножении клеточных матриц размеры клеток должны быть согласованны.
Обращение матрицы при разбиении на клетки
Разобьем матрицу А на четыре клетки
Будем искать обратную матрицу А-1также в виде четырехклеточной матрицы:
Так как
,
то перемножая эти матрицы, получим
четыре уравнения:
(1)
Es, Er - единичные матрицы (rxr) (sxs)
Решая уравнения (1), получим:
Вначале находим β12 и β22, а затем β11 и β21
В этом случае приходится обращать матрицы меньшей размерности:
что дает существенный выигрыш в памяти ЭВМ.
-
Нахождение обратной матрицы методом Гаусса
Пусть дана неособенная матрица А и обратная – А-1
Для вычисления элементов обратной
матрицы xij используем
соотношение
Умножая матрицы А и А-1 и сравнивая каждый элемент произведения соответствующему элементу матрицы Е, получим систему уравнений с n2 неизвестными xij(i,j (1,…n))
Умноженная почленно все строки матрицы А на первый столбец матрицы А-1 получим первые n уравнений:
(1)
Умножая почленно все строки матрицы А на второй столбец матрицы А-1 получим вторые n уравнений:
(2)
И так далее.
В общем виде система n2
уравнений
может быть записана как
Все n уравнений (n систем вида (1), (2)) имеют одну и ту же матрицу А и различные свободные члены. Решаем систему методом Гаусса.
-
Применение метода итераций для умножения элементов обратной матрицы
Найдем обратную матрицу методом Гаусса.
из-за погрешностей вычислений.
Будем считать полученную матрицу
- первым приближением к обратной матрице.
Для уточнения элементов матрицы строим
итерационный процесс (Демидович и Марон,
и др.)
(1)
(2)
Если
,
то итерационный процесс сходится.
Процесс (1), (2) продолжают до тех пор, пока
элементы матрицы F по
модулю не станут меньше заданного числа
Е. Тогда полагают
(можно рассматривать какую-либо норму
матрицы F)
Пример:
Уточнить элементы матрицы А-1. Итерации продолжать до тех пор, пока элементы матрицы Fk по модулю не станут ≤5*10-5
Решение:
-
Находим F0 по формуле (1)
-
Находим D0F0
-
Находим
ЛЕКЦИЯ6