Главная страница
Навигация по странице:

  • Второй способ нахождения обратной матрицы.

  • РТС. Обращение матриц Первый способ нахождения обратной матрицы


    Скачать 138.28 Kb.
    НазваниеОбращение матриц Первый способ нахождения обратной матрицы
    Дата12.03.2019
    Размер138.28 Kb.
    Формат файлаdocx
    Имя файлаРТС.docx
    ТипДокументы
    #70111

    Обращение матриц

    Первый способ нахождения обратной матрицы. Пусть дана матрица a = \begin{pmatrix} 0 & 1 & 3 \\ 2 & 3 & 5 \\ 3 & 5 & 7 \end{pmatrix}. Обратную матрицу можно вычислить по формуле a^{-1}=(\det a)^{-1} \cdot a^{t}, где a^{t} — транспонированная матрица алгебраических дополнений. Найдем определитель этой матрицы по правилу треугольника. \det a=0 \cdot 3 \cdot 7+1 \cdot 5 \cdot 3+2 \cdot 5 \cdot 3-3 \cdot 3 \cdot 3-5 \cdot 5 \cdot 0-2 \cdot 1 \cdot 7=4. Если бы определитель был равен нулю, то обратная матрица не существует. Дальше найдем алгебраическое дополнение матрицы. Чтобы найти алгебраическое дополнение каждого элемента матрицы, нужно вычеркнуть строку и столбец содержащий этот элемент, найти определитель минора каждого элемента и умножить на -1 в степени суммы номера строки и столбца в которых располагается элемент.
    a_{11}=(-1)^{1+1} \begin{pmatrix} 3 & 5 \\ 5 & 7 \\ \end{pmatrix}=-4
    a_{12}=(-1)^{1+2} \begin{pmatrix} 2 & 5 \\ 3 & 7 \\ \end{pmatrix}=1
    a_{13}=(-1)^{1+3} \begin{pmatrix} 2 & 3 \\ 3 & 5 \\ \end{pmatrix}=1
    a_{21}=(-1)^{2+1} \begin{pmatrix} 1 & 3 \\ 5 & 7 \\ \end{pmatrix}=8
    a_{22}=(-1)^{2+2} \begin{pmatrix} 0 & 3 \\ 3 & 7 \\ \end{pmatrix}=-9
    a_{23}=(-1)^{2+3} \begin{pmatrix} 0 & 1 \\ 3 & 5 \\ \end{pmatrix}=3
    a_{31}=(-1)^{3+1} \begin{pmatrix} 1 & 3 \\ 3 & 5 \\ \end{pmatrix}=-4
    a_{32}=(-1)^{3+2} \begin{pmatrix} 0 & 3 \\ 2 & 5 \\ \end{pmatrix}=6
    a_{33}=(-1)^{3+3} \begin{pmatrix} 0 & 1 \\ 2 & 3 \\ \end{pmatrix}=-2
    Матрица алгебраических дополнений a = \begin{pmatrix} -4 & 1 & 1 \\ 8 & -9 & 3 \\ -4 & 6 & -2 \end{pmatrix}. Транспонируем Матрицу алгебраических дополненийa^{t} = \begin{pmatrix} -4 & 8 & -4 \\ 1 & -9 & 6 \\ 1 & -3 & -2 \end{pmatrix}. Теперь найдем обратную матрицу a^{-1}=\frac{1}{4} \begin{pmatrix} -4 & 8 & -4 \\ 1 & -9 & 6 \\ 1 & -3 & -2 \end{pmatrix}= \begin{pmatrix} -1 & 2 & -1 \\ 1/4 & -9/4 & 3/2 \\ 1/4 & -3/4 & -1/2 \end{pmatrix}. Если обратная матрица найдена правильно, то при умножение обратной матрицы на исходную получим матрицу, у которой на главной диагонали единицы, а все остальные элементы равны нулю. \begin{pmatrix} 0 & 1 & 3 \\ 2 & 3 & 5 \\ 3 & 5 & 7 \end{pmatrix} \begin{pmatrix} -1 & 2 & -1 \\ 1/4 & -9/4 & 3/2 \\ 1/4 & -3/4 & -1/2 \end{pmatrix}= \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}. Так как получили единичную матрицу, то обратная матрица найдена верно.
    Второй способ нахождения обратной матрицы. Запишем рядом с исходной матрицей единичную \begin{pmatrix} 0 & 1 & 3 \\ 2 & 3 & 5 \\ 3 & 5 & 7 \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}. Любую матрицу можно привести к единичной, это мы и сделаем с нашей матрицей a, выполняя действия по привидению матрицы a к единичному виду, будем выполнять такие же с единичной матрицей.
    \begin{pmatrix} 0 & 1 & 3 \\ 2 & 3 & 5 \\ 3 & 5 & 7 \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix}
    Умножим вторую строку на -1 и прибавим к третьей.
    \begin{pmatrix} 0 & 1 & 3 \\ 2 & 3 & 5 \\ 1 & 2 & 2 \end{pmatrix} \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & -1 & 1 \end{pmatrix}
    Поменяем первую и третью строки местами.
    \begin{pmatrix} 1 & 2 & 2 \\ 2 & 3 & 5 \\ 0 & 1 & 3 \end{pmatrix} \begin{pmatrix} 0 & -1 & 1 \\ 0 & 1 & 0 \\ 1 & 0 & 0 \end{pmatrix}
    Первую строку умножим на -2 и прибавим ко второй.
    \begin{pmatrix} 1 & 2 & 2 \\ 0 & -1 & 1 \\ 0 & 1 & 3 \end{pmatrix} \begin{pmatrix} 0 & -1 & 1 \\ 0 & 3 & -2 \\ 1 & 0 & 0 \end{pmatrix}
    Вторую строку прибавим к третьей.
    \begin{pmatrix} 1 & 2 & 2 \\ 0 & -1 & 1 \\ 0 & 0 & 4 \end{pmatrix} \begin{pmatrix} 0 & -1 & 1 \\ 0 & 3 & -2 \\ 1 & 3 & -2 \end{pmatrix}
    Поделим третью строку на четыре.
    \begin{pmatrix} 1 & 2 & 2 \\ 0 & -1 & 1 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} 0 & -1 & 1 \\ 0 & 3 & -2 \\ 1/4 & 3/4 & -1/2 \end{pmatrix}
    Умножим вторую строку на -2 и прибавим к первой.
    \begin{pmatrix} 1 & 4 & 0 \\ 0 & -1 & 1 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} 0 & -7 & 5 \\ 0 & 3 & -2 \\ 1/4 & 3/4 & -1/2 \end{pmatrix}
    Умножим третью строку на -1 и прибавим ко второй.
    \begin{pmatrix} 1 & 4 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} 0 & -7 & 5 \\ -1/4 & 9/4 & -3/2 \\ 1/4 & 3/4 & -1/2 \end{pmatrix}
    Умножим вторую строку на -1.
    \begin{pmatrix} 1 & 4 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} 0 & -7 & 5 \\ 1/4 & -9/4 & 3/2 \\ 1/4 & 3/4 & -1/2 \end{pmatrix}
    Вторую строку умножим на -4 и прибавим к первой.
    \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} \begin{pmatrix} -1 & 2 & -1 \\ 1/4 & -9/4 & 3/2 \\ 1/4 & 3/4 & -1/2 \end{pmatrix}
    Полученная матрица является обратной.

    Пусть имеется система линейных уравнений

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image713.gif (33).

    Если уравнение (33) умножить слева и справа на обратную матрицу C–1

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image715.gif ,

    то, учитывая, что

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image717.gif ,

    где E – единичная матрица, получим формулу для решения системы методом обращения матриц:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image719.gif . (34)

    Сложность этого метода заключается в нахождении обратной матрицы С-1, которая рассчитывается следующим образом.

    Находится транспонированная матрица СТ.

    Если https://konspekta.net/studopediainfo/baza1/868369093313.files/image721.gif , то https://konspekta.net/studopediainfo/baza1/868369093313.files/image723.gif .

    Затем рассчитывается матрица алгебраических дополнений:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image725.gif ,

    где Сi,j – алгебраические дополнения элементов Сi,j ( https://konspekta.net/studopediainfo/baza1/868369093313.files/image727.gif ), которые находятся следующим образом: из транспонированной матрицы вычеркивается i-я строка и j-й столбец, определитель оставшейся части записывается в элемент матрицы алгебраических дополнений С*i,j. Знак «–» ставится перед определителем в том случае, если сумма индексов определителя является нечетным числом.

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image729.gif , https://konspekta.net/studopediainfo/baza1/868369093313.files/image731.gif , https://konspekta.net/studopediainfo/baza1/868369093313.files/image733.gif ,

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image735.gif , https://konspekta.net/studopediainfo/baza1/868369093313.files/image737.gif , https://konspekta.net/studopediainfo/baza1/868369093313.files/image739.gif ,

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image741.gif , https://konspekta.net/studopediainfo/baza1/868369093313.files/image743.gif , https://konspekta.net/studopediainfo/baza1/868369093313.files/image745.gif .

     

    Элементы обратной матрицы ищутся из элементов матрицы алгебраических дополнений по формуле:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image747.gif ,

    где det C – определитель матрицы С.

    В Mathcad существует встроенная функция для расчета обратной матрицы. Она вызывается нажатием кнопки Inverse (Инверсия) на панели Matrix (Матрицы) (рис. 41).

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image749.jpg https://konspekta.net/studopediainfo/baza1/868369093313.files/image751.jpg

    Рис. 41. Вызов вычисления обратной матрицы

    Так каксогласно (34)https://konspekta.net/studopediainfo/baza1/868369093313.files/image753.gif ,

    имеем

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image755.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image757.gif ,

    где zij – элементы обратной матрицы С-1.

    Проведя умножение матрицы на столбец, получим выражения для каждого коэффициента:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image759.gif ,

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image761.gif ,

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image763.gif .

    Пример решения системы линейных алгебраических уравнений методом обращения матриц.

    Имеем систему (32):

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image638.gif

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

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image237.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image648.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image650.gif

    2. Транспонируем матрицу:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image765.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image767.gif

    3. Найдем матрицу алгебраических дополнений:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image769.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image771.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image773.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image775.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image777.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image779.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image781.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image783.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image785.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image787.gif

    4. Найдем обратную матрицу и осуществим проверку с помощью встроенной функции Mathcad:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image789.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image791.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image793.gif

    5. Найдем решение системы:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image795.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image797.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image799.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image657.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image659.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image661.gif

    6. Осуществим проверку решения:

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image801.gif

    https://konspekta.net/studopediainfo/baza1/868369093313.files/image663.gif https://konspekta.net/studopediainfo/baza1/868369093313.files/image665.gif

    Результаты совпали, следовательно, решение верно.


    написать администратору сайта