Page 61 - 363_
P. 61

62




                        ans =
                                  8     -3   -10
                                 -5    10    16
                                 -2     4      9
                   >> A^(-2)
                        ans =
                                  l.5385e-001    -7.6923е-002    3.0769e-001
                                 7.6923e-002     3.0769e-001   -4.6154e-001
                                 2.1328e-018    -1.5385e-001    3.8462e-001
                    Розглянемо більш складні операції.

                   Обчислення  визначника  матриці  здійснюється  за  допомогою  функції
            det(A). Матриця А повинна бути квадратною.

                   Функція  обернення  матриці  (inv(A))  –  обчислює  матрицю,  обернену  до

            заданої  матриці  А.  Для  застосування  цієї  функції  початкова  матриця  А
            повинна бути квадратною і її визначник не повинен дорівнювати нулю.

                   Наведемо приклад:
                    >> C = A* B'
                         C =
                                  -40    115   -94    299
                    >> det(C)
                         ans =
                                   -1150
                    >> inv(C)
                         ans =
                                    -2.6000e-001    1.0000e-001
                                    -8.1739e-002    3.4783e-002
                    >> inv(ans)
                         ans =
                                   -4.0000e+001    1.1500e+002
                                   -9.4000e+001    2.9900e+002
                   Вельми оригінальними у мові MatLAB є дві нові, невідомі у математиці
            функції ділення матриць. Причому вводяться поняття ділення матриць зліва

            направо та ділення матриць справа  наліво.  Перша операція записується за
            допомогою знака “ / ”, а друга – за допомогою знака “ \ ”, які вміщуються між

            іменами матриць, що діляться одна на одну.
                   Операція  В  /  А  є  рівносильною  послідовності  дій  В*inv(A).  Її  зручно

            використовувати для розв'язання матричного рівняння:

                                                          X * A = В.
                   Аналогічно операція  A\B є рівносильною сукупності операцій  inv(A)*В,

            що є розв'язком матричного рівняння:

                                                          А * X = В.
                   Для прикладу розглянемо задачу знаходження коренів системи лінійних
   56   57   58   59   60   61   62   63   64   65   66