Метод исключения Гаусса с выбором главного элемента по столбцу

100779
знаков
18
таблиц
23
изображения

3.3 Метод исключения Гаусса с выбором главного элемента по столбцу

Хотя метод Гаусса является точным методом, ошибки округления могут привести к существенным погрешностям результата. Кроме того исключение по формулам (3.7) нельзя проводить, если элемент главной диагонали a равен нулю. Если элемент a мал, то велики ошибки округления при делении на этот элемент. Для уменьшения ошибок округления применяют метод исключения Гаусса с выбором главного элемента по столбцу. Прямой ход так же, как и для схемы единственного деления, состоит из n – 1 шагов. На первом шаге прежде, чем исключать переменную x1, уравнения переставляются так, чтобы в левом верхнем углу был наибольший по модулю коэффициент ai1, i = 1, 2, …, n. В дальнейшем, на k-м шаге, прежде, чем исключать переменную xk, уравнения переставляются так, чтобы в левом верхнем углу был наибольший по модулю коэффициент aik, i = k, k + 1, …, n. После этой перестановки исключение переменной xk производят, как в схеме единственного деления.

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

Пример 3.2.

Применим метод исключения Гаусса с выбором главного элемента по столбцу для решения системы уравнений (3.10) из примера 3.1. Прямой ход. 1-ый шаг. Так как коэффициент a11 = 2.0 наибольший из коэффициентов первого столбца, перестановки строк не требуется и 1-ый шаг полностью совпадает с 1-ым шагом примера 3.1. Из второго, третьего и четвертого уравнений исключается переменная x1 и система приводится к виду (3.11).

2-ой шаг. Наибольший по модулю коэффициент при x2 в системе (3.11) a = –1.15. Поэтому переставим уравнения следующим образом:

2.0x1 + 1.0x2 – 0.1x3 + 1.0x4 = 2.7

–1.15x2 + 1.015x3 + 5.05x4 = – 4.305 (3.14)

0.3x2 + 4.02x3 – 8.70x4 = 21.36

– 0.30x2 + 2.55x3  – 1.50x4 = 8.55

Вычислим множители:

 

m =  =  = –0.26087 m =  =  = 0.26087.

Вычитая из третьего и четвертого уравнений системы (3.14) второе уравнение, умноженное соответственно на m и m, приходим к системе:

2.0x1 + 1.0x2 – 0.1x3 + 1.0x4 = 2.7

–1.15x2 + 1.015x3 + 5.05x4 = – 4.305 (3.15)

4.28478x3 – 7.38261x4 = 20.23696

2.28522x3  – 2.81739x4 = 9.67305

3-ий шаг. Вычислим множитель:

 

m =  =  = 0.53333.

Вычитая из четвертого уравнения системы (3.15) третье, умноженное на m, приведем систему к треугольному виду:

2.0x1 + 1.0x2 – 0.1x3 + 1.0x4 = 2.7

–1.15x2 + 1.015x3 + 5.05x4 = – 4.305 (3.16)

4.28478x3 – 7.38261x4 = 20.23696

1.11998x4 = –1.11998

 

Обратный ход. Обратный ход полностью совпадает с обратным ходом примера 3.1. Решение системы имеет вид:

 

x1 = 1.000, x2 = 2.000, x3 = 3.000, x4 = – 1.000.

 


Информация о работе «Вычислительная математика»
Раздел: Математика
Количество знаков с пробелами: 100779
Количество таблиц: 18
Количество изображений: 23

Похожие работы

Скачать
10356
0
0

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

Скачать
7721
1
1

... если  - предельная абсолютная погрешность приближённого числа , то  (1.2) отсюда следует, что  (1.3) Значение предельной абсолютной погрешности, обычно, подбирается интуитивно по смыслу задачи. Пример 2: Определить предельную абсолютную погрешность числа , заменяющего число , точное значение которого нам неизвестно. Так как мы знаем, что , ...

Скачать
37333
0
0

... удивили меня…, хоть речь идёт обо мне самой. Они действительно написаны прекрасным стилем, который превосходит стиль самого очерка" /2/. 2.3. Рождение первенца и критическое перенапряжение Августа Ада Лавлейс работает с большим напряжением. В письмах к Бэббиджу она неоднократно жалуется на утомление, болезни, плохое самочувствие. Наконец, 6 августа Бэббидж отсылает Аде свои последние замечания ...

Скачать
54819
0
0

... в Украине, бывшем Советском Союзе и за рубежом научная школа теоретического программирования. В 2001-м году ее не стало... Но не только в научном плане велика роль женщин в развитии вычислительной техники. Со временем образуется огромное количество различных фирм по разработке и продаже программного и аппаратного обеспечения. Следовательно, разыгрываются человеческие трагедии капиталистического ...

0 комментариев


Наверх