Как сделать обратную матрицу в smath studio
Продолжаем разговор о действиях с матрицами. А именно – в ходе изучения данной лекции вы научитесь находить обратную матрицу. Научитесь. Даже если с математикой туго.
Что такое обратная матрица? Здесь можно провести аналогию с обратными числами: рассмотрим, например, оптимистичное число 5 и обратное ему число . Произведение данных чисел равно единице: . С матрицами всё похоже! Произведение матрицы на обратную ей матрицу равно – единичной матрице, которая является матричным аналогом числовой единицы. Однако обо всём по порядку – сначала решим важный практический вопрос, а именно, научимся эту самую обратную матрицу находить.
Что необходимо знать и уметь для нахождения обратной матрицы? Вы должны уметь решать определители. Вы должны понимать, что такое матрица и уметь выполнять некоторые действия с ними.
Есть? Тогда поехали дальше. А хотя… ехать могут все, если что-то не знаете, я буду ставить нужную ссылку по ходу объяснений.
Существует два основных метода нахождения обратной матрицы:
с помощью алгебраических дополнений и с помощью элементарных преобразований.
Сегодня мы изучим первый, более простой способ.
Начнем с самого ужасного и непонятного. Рассмотрим квадратную матрицу . Обратную матрицу можно найти по следующей формуле:
, где – определитель матрицы , – транспонированная матрица алгебраических дополнений соответствующих элементов матрицы .
Обозначения: Как вы уже, наверное, заметили, обратная матрица обозначается надстрочным индексом
Найти обратную матрицу для матрицы
Решаем. Последовательность действий удобно разложить по пунктам.
1) Сначала находим определитель матрицы.
Если с пониманием сего действа плоховато, ознакомьтесь с материалом Как вычислить определитель?
Важно! В том случае, если определитель матрицы равен НУЛЮ – обратной матрицы НЕ СУЩЕСТВУЕТ.
В рассматриваемом примере, как выяснилось, , а значит, всё в порядке.
2) Находим матрицу миноров .
Для решения нашей задачи не обязательно знать, что такое минор, однако, желательно ознакомиться со статьей Как вычислить определитель.
Матрица миноров имеет такие же размеры, как и матрица , то есть в данном случае .
Дело за малым, осталось найти четыре числа и поставить их вместо звездочек.
Возвращаемся к нашей матрице
Сначала рассмотрим левый верхний элемент:
Как найти его минор?
А делается это так: МЫСЛЕННО вычеркиваем строку и столбец, в котором находится данный элемент:
Оставшееся число и является минором данного элемента, которое записываем в нашу матрицу миноров:
Рассматриваем следующий элемент матрицы :
Мысленно вычеркиваем строку и столбец, в котором стоит данный элемент:
То, что осталось, и есть минор данного элемента, который записываем в нашу матрицу:
Аналогично рассматриваем элементы второй строки и находим их миноры:
Готово.
– матрица миноров соответствующих элементов матрицы .
3) Находим матрицу алгебраических дополнений .
Это просто. В матрице миноров нужно ПОМЕНЯТЬ ЗНАКИ у двух чисел:
Именно у этих чисел, которые я обвел в кружок!
– матрица алгебраических дополнений соответствующих элементов матрицы .
4) Находим транспонированную матрицу алгебраических дополнений .
Что такое транспонирование матрицы, и с чем это едят, смотрите в лекции Действия с матрицами.
– транспонированная матрица алгебраических дополнений соответствующих элементов матрицы .
5) Ответ.
Вспоминаем нашу формулу
Всё найдено!
Таким образом, обратная матрица:
Ответ лучше оставить в таком виде. НЕ НУЖНО делить каждый элемент матрицы на 2, так как получатся дробные числа. Более подробно данный нюанс рассмотрен в той же статье Действия с матрицами.
Как проверить решение?
Необходимо выполнить матричное умножение либо
Получена уже упомянутая единичная матрица – это матрица с единицами на главной диагонали и нулями в остальных местах.
Таким образом, обратная матрица найдена правильно.
Если провести действие , то в результате тоже получится единичная матрица. Это один из немногих случаев, когда умножение матриц перестановочно, более подробную информацию можно найти в статье Свойства операций над матрицами. Матричные выражения. Также заметьте, что в ходе проверки константа (дробь) выносится вперёд и обрабатывается в самом конце – после матричного умножения. Это стандартный приём.
Найти обратную матрицу для матрицы
Обратную матрицу найдем по формуле: , где – транспонированная матрица алгебраических дополнений соответствующих элементов матрицы .
1) Находим определитель матрицы.
Здесь определитель раскрыт по первой строке.
Также не забываем, что , а значит, всё нормально – обратная матрица существует.
2) Находим матрицу миноров .
Я подробно рассмотрю парочку миноров:
Рассмотрим следующий элемент матрицы:
МЫСЛЕННО вычеркиваем строку и столбец, в котором находится данный элемент:
Ну и для закрепления – нахождение еще одного минора в картинках:
Остальные миноры попробуйте вычислить самостоятельно.
Окончательный результат:
– матрица миноров соответствующих элементов матрицы .
То, что все миноры получились отрицательными – чистая случайность.
3) Находим матрицу алгебраических дополнений .
В матрице миноров необходимо СМЕНИТЬ ЗНАКИ строго у следующих элементов:
В данном случае:
– матрица алгебраических дополнений соответствующих элементов матрицы .
4) Находим транспонированную матрицу алгебраических дополнений .
– транспонированная матрица алгебраических дополнений соответствующих элементов матрицы .
5) Ответ:
Таким образом, обратная матрица найдена правильно.
Как оформить решение на чистовик? Примерный образец чистового оформления задания можно найти на странице Правило Крамера. Метод обратной матрицы в параграфе, где идет речь о матричном методе решения системы линейных уравнений. По существу, основная часть упомянутой задачи – и есть поиск обратной матрицы.
В ряде учебников, методичек можно встретить несколько другой подход к нахождению обратной матрицы, однако я рекомендую пользоваться именно вышеизложенным алгоритмом решения. Почему? Потому что вероятность запутаться в вычислениях и знаках – гораздо меньше.
Иногда обратную матрицу требуется найти методом Гаусса-Жордана, но второй способ доступен для студентов с приличной техникой элементарных преобразований.
Автор: Емелин Александр
(Переход на главную страницу)
cкидкa 15% на первый зaкaз, при оформлении введите прoмoкoд: 5530-hihi5
Решение нелинейных уравнений
Пример 3.1. Найти корень уравнения x – sin x – 0,25 = 0 на отрезке [0,2] с точностью 0,0001.
MathCAD
Для решения одного нелинейного уравнения с одной неизвестной система MathCAD имеет встроенную функцию, которая в зависимости от типа задачи может иметь или два или четыре аргумента и, соответственно, работает несколько по-разному: root(f(x),x); root(f(x),x,a,b), где f(x) - скалярная функция, определяющая исходное нелинейное уравнение (4.1); х – скалярная переменная, относительно которой решается уравнение; а, b – границы интервала, внутри которопй происходит поиск корня.
Иногда удобнее задавать не начальное приближение к корню, а интервал [а, b], внутри которого заведомо находится корень. В этом случае следует использовать функцию root с четырьмя аргументами; присваивать начальное значение переменной х в этом случае не нужно. Поиск корня будет осуществлен в промежутке между а и b альтернативным численным методом (Риддера или Брента).
| Решение (распечатка MathCAD) |
| f(x) := x – sin(x) – 0.25 x := 0,0.1..2 х := 1 – начальное приближение s := root(f(x), x) s =1.171 - корень уравнения f(s)=–4.836?10 –5 - погрешность |
SMath Studio
Для решения одного нелинейного уравнения с одной неизвестной система SMath Studio имеет встроенную функцию, которая в зависимости от типа задачи может иметь или два или три аргумента: roots(f(x); x) или roots(f(x); x; a), где f(x) – скалярная функция, определяющая исходное нелинейное уравнение (4.1); х – скалярная переменная, относительно которой решается уравнение; а – начальное приближение переменной х.
Решение (распечатка SMath Studio)
Решение систем линейных алгебраических уравнений
Пример 3.2. Решить систему уравнений
MathCAD
Для решения систем уравнений в среде MathCAD можно применять вычислительный блок Given/Find. Он состоит из трех частей:
1. ключевое слово Given;
2. система уравнений, записанная с помощью логических операторов;
3. вызов встроенной функции Find(xl, х2. хп).
Перед применением блока Given/Find необходимо задать начальные значения переменным xl, x2,хп.
Однако более наглядным является решение СЛАУ в матричной форме. В этом случае используется встроенная функция lsolve(A,b), где А - матрица коэффициентов системы, b - вектор правых частей.
Решение (распечатка MathCAD)
| Решение с помощью вычислительного блока Given/Find | Решение с помощью встроенной функции lsolve |
| xl :=1 х2:=1 х3:= 1 х4:=1 Given 2·xl + 4·х2 – 3·хЗ + х4 = 2 3·xl + х2 – 2·хЗ - х4 = 0 4·xl +11·х2 + 7·х3 + 2·х4 = 3 xl – x2 + 5·x3+2·x4 = 8 y := Find(xl ,х2,хЗ,х4) |
SMath Studio
Для решения систем уравнений в среде SMath Studio применяется матричный метод решения: коэффициенты левых частей уравнений записываются в квадратную матрицу A размером 4?4, свободные коэффициенты записываются в столбец B размером 4?1. Решением будет столбец X, который находится умножением обратной матрицы на столбец B: . Решение существует, если определитель матрицы A не равен нулю. Получив решение, нужно сделать проверку, т.е. убедиться, что .
Читайте также: