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

UptoLike

ВУЗ: 

АГТА | Ангарск

Составители: 

Рубрика: 

63
Начало
ввод N,
ε
, расши-
ренная матрица А
i := 1, N
2
3
j := i
max := 0
j
≤
N
да
k := j; max:=abs(A[j,i])
max < abs(A[j,i])
нет
да
нет
j := j+1
Прямой ход
Поиск максимального по
модулю элемента в i –ом столбце
j
:= i, N+1
j
<> k
нет
да
C := A[k, j]
A[k, j] := A[i, j]
A[i, j] := C
1
Меняем местами
строки i и k
max <
ε
нет
печать
“Бесчисленное
множество
решений”
Exit
да
Ведущий элемент не на
главной диагонали?
64
C := A[i, i]
A[i, j] := A[i, j]/C
1
j
:= i, N+1
k := i+1, N
2
C := -A[k, i]
A[k, j] := A[k, j]+C*A[i, j]
j
:= i, N+1
Деление строки
i на A[i,i ]
Получение единицы
на главной 
                            Начало


                      ввод N, ε, расши-
                      ренная матрица А
                                                                 Прямой ход
            2              i := 1, N                   3
                                                                                                            1
                                                                                                                               Деление строки
                                                                                                                               i на A[i,i ]
                            j := i
                           max := 0                                                                   C := A[i, i]
                                                           Поиск максимального по
                                                           модулю элемента в i –ом столбце
                                                                              да
                                                                                                       j := i, N+1
                                            нет
                             j≤N                                 max < ε

                                 да                                  нет                           A[i, j] := A[i, j]/C                  Получение единицы
                                                                               печать
                                                           нет             “Бесчисленное                                                 на главной
                       max < abs(A[j,i])                                     множество
                                                                             решений”
                                 да
                   k := j; max:=abs(A[j,i])
                                                                               Exit                   k := i+1, N                   2

                           j := j+1
                                                                                                      C := -A[k, i]

                                                                       Ведущий элемент не на
                                                                       главной диагонали?              j := i, N+1
                                            нет
                             j <> k
Меняем местами                                                                                 A[k, j] := A[k, j]+C*A[i, j]
                                 да
    строки i и k
                          j := i, N+1                        1


                       C := A[k, j]
                       A[k, j] := A[i, j]
                       A[i, j] := C




                                                  63                                                                      64

Страницы