Лабораторная работа №2. Решение систем линейных уравнений Цель работы

Изучить методы решения систем линейных уравнений. Разработать алгоритмы и программы для заданных методов.

Методы решения систем линейных уравнений

Способы решения систем линейных уравнений делятся на две группы:

точные методы, представляющие собой конечные алгоритмы для вычисления корней системы (решение систем с помощью обратной матрицы, метод Гаусса, метод Халецкого и др.),

итерационные методы, позволяющие получить решение системы с заданной точностью путем сходящихся итерационных процессов (метод итерации, метод Зейделя и др.).

Можно выделить группу методов, являющихся модификациями точных или итерационных методов и учитывающих некоторые заранее известные особенности матрицы системы. Примером является метод прогонки для решения систем с трехдиагональной матрицей.

Вследствие неизбежных округлений результаты даже точных методов являются приближенными. При использовании итерационных методов, сверх того, добавляется погрешность метода.

Рассмотрим систему линейных алгебраических уравнений относительно неизвестных :

(2.1)

Итерационные методы

Эффективное применение итерационных методов существенно зависит от удачного выбора начального приближения и быстроты сходимости процесса.

Метод итерации

Пусть дана линейная система (2.1). Предполагая, что диагональные коэффициенты не равны 0 , разрешим первое уравнение системы (2.1) относительно , второе - относительно и т. д. Тогда получим эквивалентную систему

(2.2)

где при не равно и при .

Введя матрицы

систему (2.2) можно записать в матричной форме, а любое приближение вычисляется по формуле . Напишем формулы приближений в развернутом виде:

(2.3)

Метод Зейделя

Метод Зейделя представляет собой некоторую модификацию метода итераций. Основная его идея заключается в том, что при вычислении -го приближения неизвестной учитываются уже вычисленные ранее -е приближения неизвестных

Пусть получена эквивалентная система (2.2). Выберем произвольно начальные приближения корней . Далее, предполагая, что k-ые приближения корней известны, согласно Зейделю будем строить -е приближения корней по формулам:

Достаточное условие сходимости метода

Нормой квадратной матрицы называется функционал, обозначаемый , удовлетворяющий условиям:

1. , причем тогда и только тогда, когда

2. , для любого

3. 

4. 

Процесс итерации для приведенной линейной системы (2.1) сходится к единственному ее решению, если какая-нибудь каноническая норма матрицы меньше единицы, т.е. для итерационного процесса (2.3) достаточное условие есть

Нормы матриц

1. ( -норма или неопределенная норма)

2. ( -норма или норма L1)

3. ( -норма или Евклидова норма)

Следствие. Процесс итерации сходится, если выполнены неравенства:

,

или

,

Условия окончания итерационного процесса

В качестве условия окончания итерационного процесса можно взять условие - заданная погрешность приближенного решения .

Точные методы Метод Гаусса

Метод Гаусса, его еще называют методом Гауссовых исключений, состоит в том, что систему (2.1) приводят последовательным исключением неизвестных к эквивалентной системе с треугольной матрицей:

Перейдем к вопросу как же добиться того, чтобы система стала треугольной. Существует несколько вычислительных схем метода. Представим стандартную версию метода схему единственного деления.

Из линейной алгебры известно, что если к некоторой строке системы уравнений прибавить любую линейную комбинацию любых других строк этой системы, то решение системы не изменится. Под линейной комбинацией строк понимается сумма строк, каждая из которых умножается на некоторое число.

Произвольно выберем в первом столбце ненулевой элемент, если необходимо, переставим эту строку первой, теперь нужно добиться того, чтобы в остальных строках коэффициенты при равнялись нулю. Для этого прибавим к этим строкам первую, умноженную на ,

Поступим аналогично с коэффициентами при остальных переменных. Алгоритм можно записать следующим образом:

,

,

Если какой-то ведущий элемент не равен нулю, но просто близок к нему, то в процессе вычислений может происходить сильное накопление погрешностей. Для предотвращения этой ситуации разработаны модификации стандартного метода Гаусса схемы с выбором главного элемента, которые обеспечивают, численную устойчивость процедуры решения.

В качестве ведущего элемента выбирают максимальный по модулю коэффициент при неизвестной в уравнениях с номерами . Затем уравнение, соответствующее выбранному коэффициенту с номером меняют местами с уравнением системы для того, чтобы главный элемент занял место коэффициента .

Аналогичным образом выглядит метод Гаусса с выбором главного элемента по строкам. В этом случае переставляются столбцы матриц; т.е. изменяется порядок исключения неизвестных. Описанные модификации называют схемами с частичным выбором.

Решение полученной системы, находят по рекуррентным формулам:

.

В матричной записи это означает, что сначала (прямой ход метода Гаусса) элементарными операциями над строками приводят расширенную матрицу системы к ступенчатому виду:

а затем (обратный ход метода Гаусса) эту ступенчатую матрицу преобразуют так, чтобы в первых столбцах получилась единичная матрица: