Пусть имеется квадратная матрица n-го порядка
Матрица А -1 называется обратной матрицей по отношению к матрице А, если А*А -1 = Е, где Е — единичная матрица n-го порядка.
Единичная матрица — такая квадратная матрица, у которой все элементы по главной диагонали, проходящей от левого верхнего угла к правому нижнему углу, — единицы, а остальные — нули, например:
Обратная матрица может существовать только для квадратных матриц т.е. для тех матриц, у которых число строк и столбцов совпадают.
Теорема условия существования обратной матрицы
Для того чтобы матрица имела обратную матрицу необходимо и достаточно, чтобы она была невырожденной.
Матрица А = (А1, А2,...А n) называется невырожденной , если векторы-столбцы являются линейно независимыми. Число линейно независимых векторов-столбцов матрицы называется рангом матрицы . Поэтому можно сказать, что для того, чтобы существовала обратная матрица, необходимо и достаточно, чтобы ранг матрицы равнялся ее размерности, т.е. r = n.
Алгоритм нахождения обратной матрицы
- Записать в таблицу для решения систем уравнений методом Гаусса матрицу А и справа (на место правых частей уравнений) приписать к ней матрицу Е.
- Используя преобразования Жордана, привести матрицу А к матрице, состоящей из единичных столбцов; при этом необходимо одновременно преобразовать матрицу Е.
- Если необходимо, то переставить строки (уравнения) последней таблицы так, чтобы под матрицей А исходной таблицы получилась единичная матрица Е.
- Записать обратную матрицу А -1 , которая находится в последней таблице под матрицей Е исходной таблицы.
Для матрицы А найти обратную матрицу А -1
Решение: Записываем матрицу А и справа приписываем единичную матрицу Е. Используя преобразования Жордана, приводим матрицу А к единичной матрице Е. Вычисления приведены в таблице 31.1.
Проверим правильность вычислений умножением исходной матрицы А и обратной матрицы А -1 .
В результате умножения матриц получилась единичная матрица. Следовательно, вычисления произведены правильно.
Ответ:
Решение матричных уравнений
Матричные уравнения могут иметь вид:
АХ = В, ХА = В, АХВ = С,
где А,В,С — задаваемые матрицы, Х- искомая матрица.
Матричные уравнения решаются с помощью умножения уравнения на обратные матрицы.
Например, чтобы найти матрицу из уравнения , необходимо умножить это уравнение на слева.
Следовательно, чтобы найти решение уравнения , нужно найти обратную матрицу и умножить ее на матрицу , стоящие в правой части уравнения.
Аналогично решаются другие уравнения.
Пример 2Решить уравнение АХ = В, если
Решение : Так как обратная матрица равняется (см. пример 1)
Матричный метод в экономическом анализе
Наряду с другими в находят применение также матричные методы . Эти методы базируются на линейной и векторно-матричной алгебре. Такие методы применяются для целей анализа сложных и многомерных экономических явлений. Чаще всего эти методы используются при необходимости сравнительной оценки функционирования организаций и их структурных подразделений.
В процессе применения матричных методов анализа можно выделить несколько этапов.
На первом этапе осуществляется формирование системы экономических показателей и на ее основе составляется матрица исходных данных , которая представляет собой таблицу, в которой по ее отдельным строкам показываются номера систем (i = 1,2,....,n) , а по вертикальным графам — номера показателей (j = 1,2,....,m) .
На втором этапе по каждой вертикальной графе выявляется наибольшее из имеющихся значений показателей, которое и принимается за единицу.
После этого все суммы, отраженные в данной графе делят на наибольшее значение и формируется матрица стандартизированных коэффициентов .
На третьем этапе все составные части матрицы возводят в квадрат. Если они имеют различную значимость, то каждому показателю матрицы присваивается определенный весовой коэффициент k . Величина последнего определяется экспертным путем.
На последнем, четвертом этапе найденные величины рейтинговых оценок R j группируются в порядке их увеличения или уменьшения.
Изложенные матричные методы следует использовать, например, при сравнительном анализе различных инвестиционных проектов, а также при оценке других экономических показателей деятельности организаций.
Уравнения вообще, линейные алгебраические уравнения и их системы, а также методы их решения занимают в математике, как теоретической, так и прикладной, особое место.
Это связано с тем обстоятельством, что подавляющее большинство физических, экономических, технических и даже педагогических задач могут быть описаны и решены с помощью разнообразных уравнений и их систем. В последнее время особую популярность среди исследователей, ученых и практиков приобрело математическое моделирование практически во всех предметных областях, что объясняется очевидными его преимуществами перед другими известными и апробированными методами исследования объектов различной природы, в частности, так называемых, сложных систем. Существует великое многообразие различных определений математической модели, данных учеными в разные времена, но на наш взгляд, самое удачное, это следующее утверждение. Математическая модель – это идея, выраженная уравнением. Таким образом, умение составлять и решать уравнения и их системы – неотъемлемая характеристика современного специалиста.
Для решения систем линейных алгебраических уравнений наиболее часто используются методы: Крамера, Жордана-Гаусса и матричный метод.
Матричный метод решения - метод решения с помощью обратной матрицы систем линейных алгебраических уравнений с ненулевым определителем.
Если выписать коэффициенты при неизвестных величинах xi в матрицу A, неизвестные величины собрать в вектор столбец X, а свободные члены в вектор столбец B, то систему линейных алгебраических уравнений можно записать в виде следующего матричного уравнения A · X = B, которое имеет единственное решение только тогда, когда определитель матрицы A не будет равен нулю. При этом решение системы уравнений можно найти следующим способом X = A -1 · B , где A -1 - обратная матрица.
Матричный метод решения состоит в следующем.
Пусть дана система линейных уравнений с n неизвестными:
Её можно переписать в матричной форме: AX = B , где A - основная матрица системы, B и X - столбцы свободных членов и решений системы соответственно:
Умножим это матричное уравнение слева на A -1 - матрицу, обратную к матрице A : A -1 (AX ) = A -1 B
Так как A -1 A = E , получаем X = A -1 B . Правая часть этого уравнения даст столбец решений исходной системы. Условием применимости данного метода (как и вообще существования решения неоднородной системы линейных уравнений с числом уравнений, равным числу неизвестных) является невырожденность матрицы A . Необходимым и достаточным условием этого является неравенство нулю определителя матрицы A : detA ≠ 0.
Для однородной системы линейных уравнений, то есть когда вектор B = 0 , действительно обратное правило: система AX = 0 имеет нетривиальное (то есть не нулевое) решение только если detA = 0. Такая связь между решениями однородных и неоднородных систем линейных уравнений носит название альтернативы Фредгольма.
Пример решения неоднородной системы линейных алгебраических уравнений .
Убедимся в том, что определитель матрицы, составленный из коэффициентов при неизвестных системы линейных алгебраических уравнений не равен нулю.
Следующим шагом будет вычисление алгебраических дополнений для элементов матрицы, состоящей из коэффициентов при неизвестных. Они понадобятся для нахождения обратной матрицы.
Тема 2. СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ.
Основные понятия.
Определение 1 . Системой m линейных уравнений с n неизвестными называется система вида:
где и - числа.
Определение 2 . Решением системы (I) называется такой набор неизвестных , при котором каждое уравнение этой системы обращается в тождество.
Определение 3 . Система (I) называется совместной , если она имеет хотя бы одно решение и несовместной , если она не имеет решений. Совместная система называется определенной , если она имеет единственное решение, и неопределенной в противном случае.
Определение 4 . Уравнение вида
называется нулевым , а уравнение вида
называется несовместным . Очевидно, что система уравнений, содержащая несовместное уравнение, является несовместной.
Определение 5 . Две системы линейных уравнений называются равносильными , если каждое решение одной системы служит решением другой и, наоборот, всякое решение второй системы является решением первой.
Матричная запись системы линейных уравнений.
Рассмотрим систему (I) (см. §1).
Обозначим:
Матрица коэффициентов при неизвестных
,
Матрица – столбец свободных членов
Матрица – столбец неизвестных
.
Определение 1. Матрица называется основной матрицей системы (I), а матрица - расширенной матрицей системы (I).
По определению равенства матриц системе (I) соответствует матричное равенство:
.
Правую часть этого равенства по определению произведения матриц (см. определение 3 § 5 главы 1 ) можно разложить на множители:
, т.е.
Равенство (2) называется матричной записью системы (I) .
Решение системы линейных уравнений методом Крамера.
Пусть в системе (I) (см. §1) m=n , т.е. число уравнений равно числу неизвестных, и основная матрица системы невырожденная, т.е. . Тогда система (I) из §1 имеет единственное решение
где Δ = det A называется главным определителем системы (I), Δ i получается из определителя Δ заменой i -го столбца на столбец из свободных членов системы (I).
Пример.Решить систему методом Крамера:
.
По формулам (3) .
Вычисляем определители системы:
,
,
,
.
Чтобы получить определитель , мы заменили в определителе первый столбец на столбец из свободных членов; заменяя в определителе 2-ой столбец на столбец из свободных членов, получаем ; аналогичным образом, заменяя в определителе 3-ий столбец на столбец из свободных членов, получаем . Решение системы:
Решение систем линейных уравнений с помощью обратной матрицы.
Пусть в системе(I) (см. §1) m=n и основная матрица системы невырожденная . Запишем систему (I) в матричном виде (см. §2 ):
т.к. матрица A невырожденная, то она имеет обратную матрицу (см. теорему 1 §6 главы 1 ). Умножим обе части равенства (2) на матрицу , тогда
. (3)
По определению обратной матрицы . Из равенства (3) имеем
Решить систему с помощью обратной матрицы
.
Обозначим
; ; .
В примере (§ 3)мы вычислили определитель , следовательно, матрица A имеет обратную матрицу . Тогда в силу (4) , т.е.
. (5)
Найдем матрицу (см. §6 главы 1 )
, , ,
, , ,
, , ,
,
.
Метод Гаусса.
Пусть задана система линейных уравнений:
. (I)
Требуется найти все решения системы (I) или убедиться в том, что система несовместна.
Определение 1. Назовем элементарным преобразованием системы (I) любое из трёх действий:
1) вычёркивание нулевого уравнения;
2) прибавление к обеим частям уравнения соответствующих частей другого уравнения, умноженных на число l;
3) перемена местами слагаемых в уравнениях системы так, чтобы неизвестные с одинаковыми номерами во всех уравнениях занимали одинаковые места, т.е. если, например, в 1-ом уравнении мы поменяли 2-ое и 3-е слагаемые, тогда то же самое необходимо сделать во всех уравнениях системы.
Метод Гаусса состоит в том, что система (I) с помощью элементарных преобразований приводится к равносильной системе, решение которой находится непосредственно или устанавливается её неразрешимость.
Как было описано в §2 система (I) однозначно определяется своей расширенной матрицей и любое элементарное преобразование системы (I) соответствует элементарному преобразованию расширенной матрицы:
.
Преобразование 1) соответствует вычёркиванию нулевой строки в матрице , преобразование 2) равносильно прибавлению к соответствующей строке матрицы другой её строки, умноженной на число l, преобразование 3) эквивалентно перестановке столбцов в матрице .
Легко видеть, что, наоборот, каждому элементарному преобразованию матрицы соответствует элементарное преобразование системы (I). В силу сказанного, вместо операций с системой (I) мы будем работать с расширенной матрицей этой системы.
В матрице 1-ый столбец состоит из коэффициентов при х 1 , 2-ой столбец - из коэффициентов при х 2 и т.д. В случае перестановки столбцов следует учитывать, что это условие нарушается. Например, если мы поменяем 1-ый и 2-ой столбцы местами, то теперь в 1-ом столбце будут коэффициенты при х 2 , а во 2-ом столбце - коэффициенты при х 1 .
Будем решать систему (I) методом Гаусса.
1. Вычеркнем в матрице все нулевые строки, если такие имеются (т.е. вычеркнем в системе (I) все нулевые уравнения).
2. Проверим, есть ли среди строк матрицы строка, в которой все элементы, кроме последнего, равны нулю (назовём такую строку несовместной). Очевидно, что такой строке соответствует несовместное уравнение в системе (I) , следовательно, система (I) решений не имеет и на этом процесс заканчивается.
3. Пусть матрица не содержит несовместных строк (система (I) не содержит несовместных уравнений). Если a 11 =0 , то находим в 1-ой строке какой-нибудь элемент (кроме последнего) отличный от нуля и переставляем столбцы так, чтобы в 1-ой строке на 1-ом месте не было нуля. Будем теперь считать, что (т.е. поменяем местами соответствующие слагаемые в уравнениях системы (I)).
4. Умножим 1-ую строку на и сложим результат со 2-ой строкой, затем умножим 1-ую строку на и сложим результат с 3-ей строкой и т.д. Очевидно, что этот процесс эквивалентен исключению неизвестного x 1 из всех уравнений системы (I), кроме 1-ого. В новой матрице получаем нули в 1-ом столбце под элементом a 11 :
.
5. Вычеркнем в матрице все нулевые строки, если они есть, проверим, нет ли несовместной строки (если она имеется, то система несовместна и на этом решение заканчивается). Проверим, будет ли a 22 / =0 , если да, то находим во 2-ой строке элемент, отличный от нуля и переставляем столбцы так, чтобы . Далее умножаем элементы 2-ой строки на и складываем с соответствующими элементами 3-ей строки, затем - элементы 2-ой строки на и складываем с соответствующими элементами 4-ой строки и т.д., пока не получим нули под a 22 /
.
Произведенные действия эквивалентны исключению неизвестного х 2 из всех уравнений системы (I), кроме 1-ого и 2-ого. Так как число строк конечно, поэтому через конечное число шагов мы получим, что либо система несовместна, либо мы придём к ступенчатой матрице (см. определение 2 §7 главы 1 ) :
,
Выпишем систему уравнений, соответствующую матрице . Эта система равносильна системе (I)
.
Из последнего уравнения выражаем ; подставляем в предыдущее уравнение, находим и т.д., пока не получим .
Замечание 1. Таким образом, при решении системы (I) методом Гаусса мы приходим к одному из следующих случаев.
1. Система (I) несовместна.
2. Система (I) имеет единственное решение, если в матрице число строк равно числу неизвестных ().
3. Система (I) имеет бесчисленное множество решений, если число строк в матрице меньше числа неизвестных ().
Отсюда имеет место следующая теорема.
Теорема. Система линейных уравнений либо несовместна, либо имеет единственное решение, либо – бесконечное множество решений.
Примеры. Решить систему уравнений методом Гаусса или доказать ее несовместность:
а) ;
б) ;
в) .
а) Перепишем заданную систему в виде:
.
Мы поменяли местами 1-ое и 2-ое уравнение исходной системы, чтобы упростить вычисления (вместо дробей мы с помощью такой перестановки будем оперировать только целыми числами).
Составляем расширенную матрицу:
.
Нулевых строк нет; несовместных строк нет, ; исключим 1-ое неизвестное из всех уравнений системы, кроме 1-го. Для этого умножим элементы 1-ой строки матрицы на «-2» и сложим с соответствующими элементами 2-ой строки, что равносильно умножению 1-го уравнения на «-2» и сложению со 2-ым уравнением. Затем умножим элементы 1-ой строки на «-3» и сложим с соответствующими элементами третьей строки, т.е. умножим 2-ое уравнение заданной системы на «-3» и сложим с 3-им уравнением. Получим
.
Матрице соответствует система уравнений
Пусть дана система линейных уравнений снеизвестными:
Будем
предполагать, что основная матрица
невырожденная.
Тогда, по теореме 3.1,
существует обратная матрица
Помножив матричное уравнение
на матрицу
слева, воспользовавшись определением
3.2, а также утверждением 8) теоремы 1.1,
получим формулу, на которой основан
матричный метод решения систем линейных
уравнений:
Замечание. Отметим, что матричный метод решения систем линейных уравнений в отличие от метода Гаусса имеет ограниченное применение: этим методом могут быть решены только такие системы линейных уравнений, у которых, во-первых, число неизвестных равно числу уравнений, а во-вторых, основная матрица невырожденная.
Пример. Решить систему линейных уравнений матричным методом.
Задана
система трёх линейных уравнений с тремя
неизвестными
где
Основная матрица системы уравнений невырожденная, поскольку её определитель отличен от нуля:
Обратную
матрицу
составим одним из методов, описанных в
пункте 3.
По формуле матричного метода решения систем линейных уравнений получим
5.3. Метод Крамера
Данный метод так же, как и матричный, применим только для систем линейных уравнений, у которых число неизвестных совпадает с числом уравнений. Метод Крамера основан на одноимённой теореме:
Теорема 5.2. Система линейных уравнений снеизвестными
основная матрица которой невырожденная, имеет единственное решение, которое может быть получено по формулам
где
определитель
матрицы, полученной из основной матрицысистемы уравнений заменой её
го
столбца столбцом свободных членов.
Пример.
Найдём решение системы линейных
уравнений, рассмотренной в предыдущем
примере, методом Крамера. Основная
матрица системы уравнений невырожденная,
поскольку
Вычислим определители
По формулам, представленным в теореме 5.2, вычислим значения неизвестных:
6. Исследование систем линейных уравнений.
Базисное решение
Исследовать систему линейных уравнений – означает определить, какой является эта система – совместной или несовместной, и в случае её совместности выяснить, определённая эта система или неопределённая.
Условие совместности системы линейных уравнений даёт следующая теорема
Теорема 6.1 (Кронекера–Капелли).
Система линейных уравнений совместна тогда и только тогда, когда ранг основной матрицы системы равен рангу её расширенной матрицы:
Для совместной системы линейных уравнений вопрос о её определённости или неопределённости решается с применением следующих теорем.
Теорема 6.2. Если ранг основной матрицы совместной системы равен числу неизвестных, то система является определённой
Теорема 6.3. Если ранг основной матрицы совместной системы меньше числа неизвестных, то система является неопределённой.
Таким
образом, из сформулированных теорем
вытекает способ исследования систем
линейных алгебраических уравнений.
Пусть n
– количество неизвестных,
Тогда:
Определение 6.1. Базисным решением неопределённой системы линейных уравнений называют такое её решение, в котором все свободные неизвестные равны нулю.
Пример. Исследовать систему линейных уравнений. В случае неопределённости системы найти её базисное решение.
Вычислим ранги основной и расширенной матрицданной системы уравнений, для чего приведём расширенную (а вместе с тем и основную) матрицу системы к ступенчатому виду:
Вторую
строку матрицы сложим с её первой
строкой, умноженной на
третью строку – с первой строкой,
умноженной на
а четвёртую строку – с первой, умноженной
наполучим матрицу
К
третьей строке этой матрицы прибавим
вторую строку, умноженную на
а к четвёртой строке – первую, умноженную
на
В результате получим матрицу
удаляя из которой третью и четвёртую строки получим ступенчатую матрицу
Таким
образом,
Следовательно, данная система линейных
уравнений совместна, а поскольку величина
ранга меньше числа неизвестных, система
является неопределённой.Полученной
в результате элементарных преобразований
ступенчатой матрице соответствует
система уравнений
Неизвестные
иявляются главными, а неизвестныеи
свободными. Придавая свободным неизвестным
нулевые значения, получим базисное
решение данной системы линейных
уравнений.
Решение систем линейных алгебраических уравнений матричным методом (с помощью обратной матрицы).
Пусть система линейных алгебраических уравнений задана в матричной форме , где матрицаA имеет размерностьn наn и ее определитель отличен от нуля.
Так как , то матрицаА – обратима, то есть, существует обратная матрица. Если умножить обе части равенстванаслева, то получим формулу для нахождения матрицы-столбца неизвестных переменных. Так мы получили решение системы линейных алгебраических уравнений матричным методом.
матричным методом.
Перепишем систему уравнений в матричной форме:
Так как то СЛАУ можно решать матричным методом. С помощью обратной матрицы решение этой системы может быть найдено как.
Построим
обратную матрицу
с
помощью матрицы из алгебраических
дополнений элементов матрицыА
(при
необходимости смотрите статьюметоды
нахождения обратной матрицы):
Осталось вычислить - матрицу неизвестных переменных, умножив обратную матрицуна матрицу-столбец свободных членов(при необходимости смотрите статьюоперации над матрицами):
или в другой записи x 1 = 4, x 2 = 0, x 3 = -1 .
Основная проблема при нахождении решения систем линейных алгебраических уравнений матричным методом заключается в трудоемкости нахождения обратной матрицы, особенно для квадратных матриц порядка выше третьего.
Более подробное описание теории и дополнительные примеры смотрите в статье матричный метод решения систем линейных уравнений.
К началу страницы
Решение систем линейных уравнений методом Гаусса.
Пусть нам требуется найти решение системы из n линейных уравнений сn неизвестными переменнымиопределитель основной матрицы которой отличен от нуля.
Суть метода Гаусса состоит в последовательном исключении неизвестных переменных: сначала исключаетсяx 1 из всех уравнений системы, начиная со второго, далее исключаетсяx 2 из всех уравнений, начиная с третьего, и так далее, пока в последнем уравнении останется только неизвестная переменнаяx n . Такой процесс преобразования уравнений системы для последовательного исключения неизвестных переменных называетсяпрямым ходом метода Гаусса . После завершения прямого хода метода Гаусса из последнего уравнения находитсяx n , с помощью этого значения из предпоследнего уравнения вычисляетсяx n-1 , и так далее, из первого уравнения находитсяx 1 . Процесс вычисления неизвестных переменных при движении от последнего уравнения системы к первому называетсяобратным ходом метода Гаусса .
Кратко опишем алгоритм исключения неизвестных переменных.
Будем считать, что , так как мы всегда можем этого добиться перестановкой местами уравнений системы. Исключим неизвестную переменнуюx 1 из всех уравнений системы, начиная со второго. Для этого ко второму уравнению системы прибавим первое, умноженное на, к третьему уравнению прибавим первое, умноженное на, и так далее, кn-ому уравнению прибавим первое, умноженное на. Система уравнений после таких преобразований примет видгде, а.
К такому же результату мы бы пришли, если бы выразили x 1 через другие неизвестные переменные в первом уравнении системы и полученное выражение подставили во все остальные уравнения. Таким образом, переменнаяx 1 исключена из всех уравнений, начиная со второго.
Далее действуем аналогично, но лишь с частью полученной системы, которая отмечена на рисунке
Для этого к третьему уравнению системы прибавим второе, умноженное на , к четвертому уравнению прибавим второе, умноженное на, и так далее, кn-ому уравнению прибавим второе, умноженное на. Система уравнений после таких преобразований примет видгде, а. Таким образом, переменнаяx 2 исключена из всех уравнений, начиная с третьего.
Далее приступаем к исключению неизвестной x 3 , при этом действуем аналогично с отмеченной на рисунке частью системы
Так продолжаем прямой ход метода Гаусса пока система не примет вид
С этого момента начинаем обратный ход метода Гаусса: вычисляем x n из последнего уравнения как, с помощью полученного значенияx n находимx n-1 из предпоследнего уравнения, и так далее, находимx 1 из первого уравнения.
Решите систему линейных уравнений методом Гаусса.
Исключим
неизвестную переменную x
1
из второго и третьего уравнения системы.
Для этого к обеим частям второго и
третьего уравнений прибавим соответствующие
части первого уравнения, умноженные наи
насоответственно:
Теперь
из третьего уравнения исключим x
2
,
прибавив к его левой и правой частям
левую и правую части второго уравнения,
умноженные на:
На этом прямой ход метода Гаусса закончен, начинаем обратный ход.
Из последнего уравнения полученной системы уравнений находим x 3 :
Из второго уравнения получаем .
Из первого уравнения находим оставшуюся неизвестную переменную и этим завершаем обратный ход метода Гаусса .
x 1 = 4, x 2 = 0, x 3 = -1 .
Более детальную информацию и дополнительные примеры смотрите в разделе решение элементарных систем линейных алгебраических уравнений методом Гаусса.
К началу страницы
" |