Page 21 - 4744
P. 21

x
                                                 a N ,N  1 N 1   a N  ,N  x N   b N  ,
                                                                                                 (1.21)
                                                  x     M  x   L  ,
                                                  N 1    N  N   N
               інші значення  x  одержуємо за формулою (1.14).
                                     k
                      Якщо матрицю  A можна розкласти в добуток:
                                                          A   LU ,                               (1.22)
               де  L–  нижня  трикутна  матриця  з  елементами  рівними  1  на

               головній  діагоналі,  U –  верхня  трикутна  матриця,  то  розв’язок
               системи (1.8) зводиться до розв’язку двох систем з діагональними
               матрицями LY  – відносно Y , та UX  – відносно  X .
                                    B
                                                                   Y
                      Фактично  елементи  матриці  U   співпадають  з  елементами
               матриці  A, приведеної до трикутного виду, елементами матриці  L
               є  коефіцієнти,  на  які  треба  множити  відповідні  рівняння  для
               виключення змінних при реалізації метода Гауса.




                      1.1.4  Ітераційні  методи  розв’язку  систем  алгебраїчних

               рівнянь.
                      Розглянемо методи Якобі та Гауса-Зейделя.

                      Ітераційний  метод  Якобі  –  метод  розв’язання  системи
               рівнянь  Ax    за  умови  a        ii    0  –  передбачає  розв’язання  кожного
                                b
               рівнянняґ  системи  окремо  відносно  тільки  однієї  змінної  в

               припущенні, що всі інші змінні фіксовані. Ітераційна процедура
               методу Якобі має такий вигляд:
                                                                n
                                                         b        a ij  x  k
                                                          i
                                                                        j
                                                  k 1      j  i , 1   j
                                                 a i                     ,  i 1  2 ,   n.      (1.23)
                                                                a ii
                      Запишемо  метод  Якобі  в  матричній  формі.  Розділимо
               матрицю  A  на  матриці  –  діагональну  D  та  нижню  L  і  верхню

               U трикутні:
                                                        A   D (L   E   U ) ,                   (1.24)


















                                                                                                       21
   16   17   18   19   20   21   22   23   24   25   26