Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
154
Добавлен:
18.03.2015
Размер:
1.22 Mб
Скачать

83

 

 

 

é1

 

0,466 0ù

 

é1,45

- 0,25

0

ù

 

H = Q

(1 )

1( )

ê

 

1

ú

×

ê

0

0,537

0

ú

=

H

=

ê0

 

0ú

ê

ú

 

 

 

ê0

 

0

1ú

 

ê

0

0

0,255ú

 

 

 

 

ë

 

 

û

 

ë

 

 

 

û

 

 

 

 

 

é1,45

0

 

 

0

ù

 

 

 

 

 

 

=

ê

0

0,537

 

0

ú

 

 

 

 

 

 

ê

 

ú.

 

 

 

 

 

 

 

ê

0

0

 

 

0,255ú

 

 

 

 

 

 

 

ë

 

 

 

 

 

û

 

 

 

Вычислим результирующую верхнюю треугольную матрицу:

 

 

 

é1

0,466

0ù

 

é1

0

0,78ù

 

é1

0,466

1,45ù

Q = Q

(1 )

2( )

ê

1

ú

 

ê

1

1,43

ú

 

ê

1

 

ú

Q

=

ê0

0ú

×

ê0

ú

=

ê0

1,43ú .

 

 

 

ê0

0

1ú

 

ê0 0

1

ú

 

ê0

0

1

ú

 

 

 

ë

 

û

 

ë

 

 

û

 

ë

 

 

û

Обратную матрицу вычислим в соответствии с (4.14):

A-1 = H-1 × Q ×L

 

 

 

é0,69

0

0

ù

 

é1

0,466

1,45ù

 

é0,69

0,32

1

ù

 

H

-1

Q =

ê

0

1,862

0

ú

 

ê

1

 

ú

 

ê

0

1,86

 

ú

 

 

ê

ú

×

ê0

1,43ú

=

ê

2,66ú

,

 

 

 

ê

0

0

3,92ú

 

ê0

0

1

ú

 

ê

0

0

3,92

ú

 

 

 

 

ë

 

 

 

û

 

ë

 

 

û

 

ë

 

 

 

û

 

 

 

 

 

 

é0,69

0,32

1

ù

 

é

1

0

0ù

 

é1

1

1

ù

А

-1

= Н

-1

QL =

ê

0

1,86

 

ú

 

ê

0,172

1

ú

 

ê

3,67

 

ú

 

 

ê

2,66ú

×

ê

0ú

=

ê1

2,66ú .

 

 

 

 

 

ê

0

0

3,92

ú

 

ê0,255

0,68

1ú

 

ê1

2,66

3,92

ú

 

 

 

 

 

ë

 

 

 

û

 

ë

 

 

û

 

ë

 

 

û

4.5. Итерационные методы решения систем линейных уравнений

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

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

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