Вычислительная математика. Ч. 1. Асламова В.С - 40 стр.

UptoLike

79
Начало
ввод n,
матрицы А
u := 0
i: = 1, n
j
:= i
;
k:=0
2 3
j
n
да
нет
61],[ > EjiA
нет
j := j+1
k: = j
j: = n+1
да
k< >i
j
:= i, n
C := A[k, j]
A[i, j] := A[k, j]
A[k, j] := C
u: = u+1
1
да
нет
Приведение матрицы
к треугольному виду
k=0
нет
Печатьв
столбце’,i,’все
элементы=0
Exit
Перестановка
строк i и k
да
80
1
k
:= i+1, n
2
C := -A
[
k
,
i
]
/A
[
i
,
i
]
j
:= i
,
n
A[k, j] := A[k, j]+C*A[i, j]
Получение нулей ниже
главной диагонали
3
det := 1
i:=1, n
det := det*A
[
i
,
i
]
u mod 2<>0
det := -det
да
нет
Печать
оп
е
елитель=’
det
конец
Рис. 41. Блок-схема вычисления определителя матрицы по методу
Гаусса с выбором главного элемента в столбце.
                   Начало
                                                                                                             1
                                                                                                                                       Получение нулей ниже
                   ввод n,                                                                                                             главной диагонали
                 матрицы А
                                                                                                       k := i+1, n                 2

                    u := 0
                                                     Приведение матрицы                            C := -A[k, i]/A[i, i]
                                                     к треугольному виду
     2             i: = 1, n               3
                                                                                                         j := i, n
                 j: = i; k:=0
                                                                                              A[k, j] := A[k, j]+C*A[i, j]
                                    нет
                     j≤ n
                                                                                                            3
                           да
                                          нет
               A[i , j ] > 1E − 6
                                                                                                        det := 1
                           да                  j := j+1
                     k: = j                                                                             i := 1, n
                   j: = n+1                                         да
                                                          k=0
                                                            нет                                    det := det*A[i, i]


                                    нет
                     k< >i                                        Печать ‘в                                                  нет
Перестановка                                                      столбце’,i,’все                    u mod 2<>0
 строк i и k                да                                    элементы=0’
                                                                                                                да
                   j: = i, n
                                                                         Exit                          det := -det
                 C := A[k, j]
               A[i, j] := A[k, j]                                                                      Печать
                 A[k, j] := C
                                                                                                 ‘определитель=’, det


                   u: = u+1                                                                              конец

                                                                                    Рис. 41. Блок-схема вычисления определителя матрицы по методу
                       1                                                            Гаусса с выбором главного элемента в столбце.



                                 79                                                                                  80