- •Аппроксимация производных
- •Погрешности численного дифференцирования
- •Варианты заданий
- •Лабораторная работа №7 Численное интегрирование
- •Одномерный случай
- •Метод Гаусса
- •Метод прямоугольников
- •Метод трапеций
- •Метод парабол (метод Симпсона)
- •Увеличение точности
- •Метод Гаусса-Кронрода
- •Интегрирование при бесконечных пределах
- •Методы Монте-Карло
- •Многомерный случай
- •Варианты заданий
- •Лабораторная работа №8
- •Матричная форма
- •Методы решения
- •Метод прогонки
- •Алгоритм
- •Приложения
- •Варианты заданий
- •Лабораторная работа №9
- •Итерационные методы
- •Варианты заданий
Лабораторная работа №8
Прямые методы решения систем линейных алгебраических уравнений. Метод Гаусса, метод прогонки.
Цель работы: освоить прямые решения систем линейных алгебраических уравнений. Также научится работать с методом Гаусса, методом прогонки и методом квадратных корней.
Система линейных алгебраических уравнений
Система
- линейных
уравнений
с
неизвестными
(или, линейная
система) в
линейной
алгебре
— это система уравнений вида
8.1
Здесь
—
неизвестные, которые надо определить.
a11,
a12,
…, amn
— коэффициенты системы — и b1,
b2,
… bm
— свободные члены — предполагаются
известными. Индексы коэффициентов (aij)
системы обозначают номера уравнения
(i)
и неизвестного (j),
при котором стоит этот коэффициент,
соответственно.
Система (8.1) называется однородной, если все её свободные члены равны нулю (b1 = b2 = … = bm = 0), иначе — неоднородной.
Система (8.1) называется квадратной, если число m уравнений равно числу n неизвестных.
Решение системы (8.1) — совокупность n чисел c1, c2, …, cn, таких что подстановка каждого ci вместо xi в систему (8.1) обращает все её уравнения в тождества.
Система (8.1) называется совместной, если она имеет хотя бы одно решение, и несовместной, если у неё нет ни одного решения.
Совместная система вида (8.1) может иметь одно или более решений.
Решения c1(1), c2(1), …, cn(1) и c1(2), c2(2), …, cn(2) совместной системы вида (8.1) называются различными, если нарушается хотя бы одно из равенств:
c
8.2 |
Совместная система вида (8.1) называется определённой, если она имеет единственное решение; если же у неё есть хотя бы два различных решения, то она называется неопределённой. Если уравнений больше, чем неизвестных, она называется переопределённой.
Матричная форма
С
8.3
или:
Ax = B.
Если к матрице А приписать справа столбец свободных членов, то получившаяся матрица называется расширенной.
Методы решения
Прямые (или точные) методы, позволяют найти решение за определённое количество шагов. Итерационные методы, основаны на использовании повторяющегося процесса и позволяют получить решение в результате последовательных приближений.
Рассмотрим метод Гаусса, метод прогонки, а также метод Халецкого:
Метод Гаусса — классический метод решения системы линейных алгебраических уравнений (СЛАУ). Это метод последовательного исключения переменных, когда с помощью элементарных преобразований система уравнений приводится к равносильной системе ступенчатого (или треугольного) вида, из которого последовательно, начиная с последних (по номеру) переменных, находятся все остальные переменные.
П
8.4
Матрица A называется основной матрицей системы, b — столбцом свободных членов.
Т
8.5
При
этом будем считать, что базисный
минор
(ненулевой минор
максимального порядка) основной матрицы
находится в верхнем левом углу, то есть
в него входят только коэффициенты при
переменных
Тогда
переменные
называются главными
переменными.
Все остальные называются свободными.
Если
хотя бы одно число
,
где i
> r,
то рассматриваемая система несовместна.
Пусть,
что
для любых i
> r.
П
8.6
Если свободным переменным системы (2) придавать все возможные значения и решать новую систему относительно главных неизвестных снизу вверх (то есть от нижнего уравнения к верхнему), то мы получим все решения этой СЛАУ. Так как эта система получена путём элементарных преобразований над исходной системой (1), то по теореме об эквивалентности при элементарных преобразованиях системы (1) и (2) эквивалентны, то есть множества их решений совпадают.
Трёхдиагональная матрица
Т
8.7
Системы линейных алгебраических уравнений с такими матрицами встречаются при решении многих задач математики и физики. Краевые условия x1 и xn, которые берутся из контекста задачи, задают первую и последнюю строки. Так краевое условие первого рода F(x = x1) = F1 определит первую строку в виде C1 = 1, B1 = 0, а условие второго рода dF / dx(x = x1) = F1 будет соответствовать значениям C1 = − 1, B1 = 1.
