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

5. Решение системы линейных уравнений методом Ритца

 

Если  - симметричная и положительно-определённая матрица, то задача решения линейной системы уравнений:

(36)

эквивалентна задаче нахождения точки минимума функции многих переменных:

(37)

где скалярные произведения понимаются в смысле , т.е.

(38)

Иначе говоря, решение системы линейных уравнений (36) доставляет минимум функции многих переменных:

(39)

И наоборот, точка минимума функции (39) является решением системы линейных уравнений (36).

Таким образом, метод Ритца позволяет решение линейной системы уравнений с симметричной и положительно-определённой матрицей свести к задаче нахождения точки минимума функций многих переменных. А эту задачу мы уже умеем решать.

6. Решение системы линейных уравнений с трехдиаганальной матрицей методом прогонки Томаса

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

(40)

для решения этой линейной системы уравнений, конечно, можно применять метод Гаусса, но тогда пришлось бы делать много необязательных операций с нулями. Чтобы сэкономить время вычислений и не работать лишний раз с нулями, Томас (1949г.) разработал специальный алгоритм расчета. Рассчитывая по алгоритму Томаса элементы получаемой треугольной матрицы, мы следуем методу Гаусса, с уточнением, что с нулями никаких действий не производим; алгоритм Томаса называют – методом прогонки.

Для решения системы (40) методом прогонки – Томаса действуем следующим образом:

а) прямой ход:

(41)

Замечание: после проведения прямого хода предполагается, что все , и - неизменны (что очевидно).

б) обратный ход:

(42)

Таким образом, для системы линейных уравнений с трехдиаганальной матрицей наиболее экономным является алгоритм прогонки – Томаса, который является «отфильтрованным» методом Гаусса.

Метод минимизации невязки для решения линейной системы уравнений (метод наименьших квадратов).

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

Число наблюдаемых величин больше числа неизвестных . Пусть известно, что величины  связаны между собой линейной зависимостью:

, , . (43)

Коэффициенты  - считаются известными и неотягощенными случайными ошибками. Система (43) называется системой условных уравнений.

Если бы все числа  были точными, то неизвестные ,  могли бы быть определены из любых  - уравнений системы . Но так, как  - определены с ошибками, то система условных уравнений несовместна (переопределена, т.к. ), существуют «невязки»:

,   (44)

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

В методе наименьших квадратов, в качестве нормы рассматривают дискретную норму Гаусса:

(45)

Очевидно, что эта норма минимальна тогда, когда минимально подкоренное выражение, т.е. сумма квадратов невязок .

(46)

Условия существования минимума для функций специального вида  имеют вид:

,, (47)

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

Для примера рассмотрим уравнений с тремя неизвестными, система условных уравнений имеет вид:

(48)

Тогда система соответствующих нормальных уравнений имеет вид:

(49)

Решение системы (49) дает решение задачи (48) наилучшим приближением, в смысле дискретной нормы Гаусса.

Замечания:

1) классический метод Гаусса, метод Гаусса с выбором главного элемента, метод Якоби и метод минимизации невязки являются общими методами и применяются для определения решения невырожденных систем линейных уравнений, когда ведущие (большие по модулю) элементы матрицы системы расположены в окрестности главной диагонали (система хорошо обусловлена), если же система плохо обусловлена, тогда нужно менять соответствующую модель, чтобы она приводила к приемлемой системе уравнений;

2) для ускорения сходимости методов разработаны специальные методы – метод Гаусса-Зейделя, методы релаксации и др., которые применимы лишь для узкого класса систем – с симметрической, положительно-определенной матрицей; с ненулевыми диагональными элементами;

3) для нужд разностных уравнений разработаны специальные алгоритмы прогонки Томаса, которые являются «экономными» методами Гаусса для трехдиагональных матриц системы линейных уравнений.

 

Литература

1.   Т. Шуп. Решение интегральных задач на ЭВМ. Мир., М.,2002

2.   Л. Коллатц, Ю. Альберхт. Задачи по прикладной математике. Мир., М.,1998.

3.   Т.А. Обгадзе. Элементы математического моделирования. Учебное пособие. Грузинский Политехнический Институт им. В.И. Ленина, Тбилисси, 1999.