4.11. Решение варианта задания

Пример 1. Методом квадратного корня решить систему уравнений:
(4.47)

Решение. Все вычисления помещаем в табл. 4.1 (см. [9, с. 68]).

Прямой ход

  1. В первый раздел таблицы записываем коэффициенты системы и свободные члены.
  2. Суммируем элементы по каждой строке и получаем величины

которые записываем в последнем столбце табл. 4.1 (раздел I).

       3. Находим  с помощью формул:
(4.48)

Последовательно получаем:

4. По формулам, аналогичным (4.48), вычисляем:

Обратный ход

5. По формулам (4.45) находим :

Полученные  помещаем в табл. 4.1 в раздел II, где записаны  и .

6. Проводим контроль вычислений. Для этого находим суммы:

Сравниваем их с  Видим, что результаты совпадают до

Переходим к последнему этапу.

7. По формулам (4.46) находим :

Полученные значения  записываем в раздел III табл. 4.1.

8. Найденные значения  подставляем в исходную систему. Получаем:

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

Таблица 4.1

I

3,2

1,0

1,0

4,0

9,2

1,0

3,7

1,0

4,5

10,2

1,0

1,0

4,2

4,0

10,2

II

1,78885

0,55902

0,55902

2,23607

5,14297

 

 

1,84052

0,37353

1,76580

3,97984

 

 

 

 

1,93597

1,07974

3,01574

III

 

 

0,81127

0,84621

0,55772

 

 

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

,

где

 – симметрическая матрица.

Решение.

Будем использовать табл. 4.2:

Таблица 4.2

1.

–0,41433

0,8144

5,1518

2,1413

0,5400

8,6475

 

2.

–0,28359

5,1518

0,8114

1,4656

0,2210

7,6498

 

3.

 

2,1413

1,4556

5,1681

0,1210

8,8960

 

4.

 

–0,07280

4,54456

0

0,48987

4,96163

4,96162

5.

–0,08709

4,54455

0,39577

0

0,18669

5,12701

5,12698

6.

 

4,55089

0

0

0,14403

4,69492

4,69487

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

Решая эту систему “снизу вверх”, последовательно определим:

Для проверки правильности обратного хода схемы вычислим соответствующие :

Сравнивая  и , видим, что:

т.е.  отличаются от  на 1, а  и  на 0,99999.

Это вполне допустимая точность вычислений.

Замечание. Пояснения к табл. 4.2, приведенной выше, можно посмотреть в [9, с. 48], а смысл обозначений, используемых в табл. 4.2, – в разделе 4 данной главы.