Физический смысл метода Гаусса

Метод Гаусса имеет определенный физический смысл. Посмотрим, например, чем отличаются друг от друга исходная система

и система

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

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

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

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

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

12.8. Погрешность решения

Подставив в исходные уравнения, мы получим числа, отличающиеся от правых частей этих уравнений. Это означает, что решение имеет погрешность. Основная причина – ограниченная точность расчетов. Компьютер хранит каждое число в виде 7 десятичных цифр при обычной точности (4 байта на одно число) и 15 цифр - при двойной точности (8 байт на одно число).

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

Предположим, имеется 10 миллионов одинаковых малых чисел (равных 0.1) и одно большое 1.234.567, отличающееся от них более, чем на 7 порядков. Если вначале сложить малые числа, получим миллион, тогда общая сумма составит 2.234.567. Если же вначале взять число 1.234.567, то, сколько ни прибавляй к нему по 0.1, оно останется без изменения, если результат хранится в формате с обычной точностью (восьмая цифра числа теряется).

Устойчивость решения

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

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

Понятие вырожденной матрицы поясняет рис. 12.4, где рассмотрено решение системы из двух уравнений, каждое уравнение изображается прямой линией на графике, по осям которого отложены значения неизвестных. Решение отображает точка пересечения прямых. Если матрица близка к вырожденной (ее определитель близок к нулю), то прямые пересекаются под острым углом. Это не мешает найти решение при отсутствии погрешностей в значениях коэффициентов матрицы и процедуре решения. Однако появление данных погрешностей, даже небольших, (на рисунке это показано заменой прямых на полосы, ширина которых характеризует погрешность) приводит к резкому увеличению погрешности результата решения.

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