Метод обратной матрицы

Автор работы: Пользователь скрыл имя, 03 Декабря 2012 в 13:21, контрольная работа

Краткое описание

Процесс вычисления обратной матрицы достаточно трудоемкий. Решение системы линейных уравнений методом обратной матрицы не является самым удобным, но такой способ есть... Данный метод применим, если определитель, составленный из коэффициентов при переменных, не равен нулю.

Содержимое работы - 1 файл

Метод обратной матрицы.docx

— 107.29 Кб (Скачать файл)

Метод обратной матрицы. (Решение систем линейных уравнений)


 
Процесс вычисления обратной матрицы  достаточно трудоемкий. Решение системы  линейных уравнений методом обратной матрицы не является самым удобным, но такой способ есть... Данный метод применим, если определитель, составленный из коэффициентов при переменных, не равен нулю.


Введите исходные данные  
целые числа и ( или ) десятичные дроби ( например -0.125   -5   2.12   10 )


* x1 +

* x2 +

* x3 =

 

* x1 +

* x2 +

* x3 =

 

* x1 +

* x2 +

* x3 =

 

Обратную матрицу  найти   

     
     
     
     

Решение вашей  системы методом Гаусса ...


Решение вашей  системы методом Крамера ...


Решим систему уравнений


7 x1

+ 5 x2

-20 x3

=

-2

2 x1

+ 2 x2

-3 x3

=

17

x1

+ x2

- x3

=

11


Запишем систему  уравнений в матричной форме


A * X = B


 

   

7

 

   

5

 

 

-

20

 

*

   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


 

   

x1

 

=

   

x2

 

   

x3

 


 

 

-

2

 

 
   

17

 

   

11

 



Найдем матицу A-1, обратную к матрице А, методом алгебраических дополнений


Будем обозначать элементы матрицы A маленькими буквами аij. Первый индекс i обозначает номер строки , а второй j - номер столбца, где находится элемент матрицы аij.


A =

   

a11

 

   

a12

 

   

a13

 

 
   

a21

 

   

a22

 

   

a23

 

   

a31

 

   

a32

 

   

a33

 


Обратную матрицу A-1, будем искать в следующем виде:


A -1 = 1 / det A *

   

A11

 

   

A21

 

   

A31

 

 
   

A12

 

   

A22

 

   

A32

 

   

A13

 

   

A23

 

   

A33

 


где Aij = ( -1 ) i+j * M ij


Обратите внимание на множитель 1 / det A - стоящей перед матрицей. Очевидно, если определитель матрицы А равен нулю , то обратной матрицы не существует.


М ij это минор элемента а ij, т.е. определитель , полученный вычеркиванием из матрицы А строки с номером i и столбца с номером j. А ij - это алгебраическое дополнение элемента а ij, или, проще говоря, минор взятый с определенным знаком. Если сумма номера строки и номера столбца элемента аij четная , то алгебраическое дополнение это минор. Если сумма номера строки и номера столбца элемента аij нечетная , то алгебраическое дополнение это минор, взятый со знаком минус. Математически это выражается выражением ( -1 )i+j. Не забудьте обратить внимание на индексы алгебраических дополнений в обратной матрице.


·  Найдем определитель матрицы А.


det A =

 

7

5

-20

 

=

2

2

-3

1

1

-1


Из элементов  столбца 1 вычитаем соответствующие  элементы столбца 2 .


=

 

2

5

-20

 

=

0

2

-3

0

1

-1


Разлагаем определитель по элементам первого столбца.


= ( - 1 )1+1 * 2*

 

2

-3

 

+

1

-1


( - 1 )2+1 * 0*

 

5

-20

 

+

1

-1


( - 1 )3+1 * 0*

 

5

-20

 

=

2

-3



= 2*

 

2

-3

 

=

1

-1



= 2* ( 2 * ( -1) - ( -3) * 1 ) =


= 2 * 1

= 2


более подробное  вычисление det А


Определитель матрицы А отличен от нуля, следовательно обратная матрица A-1 существует.


·  Найдем алгебраическое дополнение A11 элемента a11 . В матрице А вычеркиваем строку 1 и столбец 1.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M11 ) элемента a11.


M11 =

 

2

-3

 

= 2 * ( -1) - ( -3) * 1 = ( -2) - ( -3) = 1

1

-1


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a11, есть число четное ( 1 + 1 = 2 ) и выражение ( -1 )1+1 = 1, то алгебраическое дополнение элемента a11 равно минору данного элемента.


A11 = ( -1 ) 1+1 * M 11 = ( -1 ) 1+1 * 1 = 1


·  Найдем алгебраическое дополнение A12 элемента a12 . В матрице А вычеркиваем строку 1 и столбец 2.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M12 ) элемента a12.


M12 =

 

2

-3

 

= 2 * ( -1) - ( -3) * 1 = ( -2) - ( -3) = 1

1

-1


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a12, есть число нечетное ( 1 + 2 = 3 ) и выражение ( -1 )1+2 = - 1, то алгебраическое дополнение элемента a12 равно минору данного элемента взятого со знаком минус.


A12 = ( -1 ) 1+2 * M 12 = ( -1 ) 1+2 * 1 = -1


·  Найдем алгебраическое дополнение A13 элемента a13 . В матрице А вычеркиваем строку 1 и столбец 3.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M13 ) элемента a13.


M13 =

 

2

2

 

= 2 * 1 - 2 * 1 = 2 - 2 = 0

1

1


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a13, есть число четное ( 1 + 3 = 4 ) и выражение ( -1 )1+3 = 1, то алгебраическое дополнение элемента a13 равно минору данного элемента.


A13 = ( -1 ) 1+3 * M 13 = ( -1 ) 1+3 * 0 = 0


·  Найдем алгебраическое дополнение A21 элемента a21 . В матрице А вычеркиваем строку 2 и столбец 1.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M21 ) элемента a21.


M21 =

 

5

-20

 

= 5 * ( -1) - ( -20) * 1 = ( -5) - ( -20) = 15

1

-1


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a21, есть число нечетное ( 2 + 1 = 3 ) и выражение ( -1 )2+1 = - 1, то алгебраическое дополнение элемента a21 равно минору данного элемента взятого со знаком минус.


A21 = ( -1 ) 2+1 * M 21 = ( -1 ) 2+1 * 15 = -15


·  Найдем алгебраическое дополнение A22 элемента a22 . В матрице А вычеркиваем строку 2 и столбец 2.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M22 ) элемента a22.


M22 =

 

7

-20

 

= 7 * ( -1) - ( -20) * 1 = ( -7) - ( -20) = 13

1

-1


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a22, есть число четное ( 2 + 2 = 4 ) и выражение ( -1 )2+2 = 1, то алгебраическое дополнение элемента a22 равно минору данного элемента.


A22 = ( -1 ) 2+2 * M 22 = ( -1 ) 2+2 * 13 = 13


·  Найдем алгебраическое дополнение A23 элемента a23 . В матрице А вычеркиваем строку 2 и столбец 3.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M23 ) элемента a23.


M23 =

 

7

5

 

= 7 * 1 - 5 * 1 = 7 - 5 = 2

1

1


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a23, есть число нечетное ( 2 + 3 = 5 ) и выражение ( -1 )2+3 = - 1, то алгебраическое дополнение элемента a23 равно минору данного элемента взятого со знаком минус.


A23 = ( -1 ) 2+3 * M 23 = ( -1 ) 2+3 * 2 = -2


·  Найдем алгебраическое дополнение A31 элемента a31 . В матрице А вычеркиваем строку 3 и столбец 1.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M31 ) элемента a31.


M31 =

 

5

-20

 

= 5 * ( -3) - ( -20) * 2 = ( -15) - ( -40) = 25

2

-3


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a31, есть число четное ( 3 + 1 = 4 ) и выражение ( -1 )3+1 = 1, то алгебраическое дополнение элемента a31 равно минору данного элемента.


A31 = ( -1 ) 3+1 * M 31 = ( -1 ) 3+1 * 25 = 25


·  Найдем алгебраическое дополнение A32 элемента a32 . В матрице А вычеркиваем строку 3 и столбец 2.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M32 ) элемента a32.


M32 =

 

7

-20

 

= 7 * ( -3) - ( -20) * 2 = ( -21) - ( -40) = 19

2

-3


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a32, есть число нечетное ( 3 + 2 = 5 ) и выражение ( -1 )3+2 = - 1, то алгебраическое дополнение элемента a32 равно минору данного элемента взятого со знаком минус.


A32 = ( -1 ) 3+2 * M 32 = ( -1 ) 3+2 * 19 = -19


·  Найдем алгебраическое дополнение A33 элемента a33 . В матрице А вычеркиваем строку 3 и столбец 3.


A =

   

7

 

   

5

 

 

-

20

 

 
   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


Определитель состоящий  из оставшихся элементов матрицы  А, называется минором ( M33 ) элемента a33.


M33 =

 

7

5

 

= 7 * 2 - 5 * 2 = 14 - 10 = 4

2

2


Так как сумма  номера строки и номера столбца ,на пересечении которых находится элемент a33, есть число четное ( 3 + 3 = 6 ) и выражение ( -1 )3+3 = 1, то алгебраическое дополнение элемента a33 равно минору данного элемента.


A33 = ( -1 ) 3+3 * M 33 = ( -1 ) 3+3 * 4 = 4


Осталось, только записать обратную матрицу.


A -1 = 1 / 2 *

   

1

 

 

-

15

 

   

25

 

 
 

-

1

 

   

13

 

 

-

19

 

   

0

 

 

-

2

 

   

4

 


A -1 =

   

1

 
 

2


 

-

15

 
 

2


   

25

 
 

2


 
 

-

1

 
 

2


   

13

 
 

2


 

-

19

 
 

2


   

0

 

 

-

1

 

   

2

 


Вернемся  к нашему уравнению, которое мы записали в матричной форме.


A * X = B


Умножим левую и  правую часть нашего матричного уравнения  на A-1


A-1 * A * X = A-1 * B


 

   

1

 
 

2


 

-

15

 
 

2


   

25

 
 

2


*

 

-

1

 
 

2


   

13

 
 

2


 

-

19

 
 

2


   

0

 

 

-

1

 

   

2

 


 

   

7

 

   

5

 

 

-

20

 

*

   

2

 

   

2

 

 

-

3

 

   

1

 

   

1

 

 

-

1

 


 

   

x1

 

=

   

x2

 

   

x3

 


 

   

1

 
 

2


 

-

15

 
 

2


   

25

 
 

2


*

 

-

1

 
 

2


   

13

 
 

2


 

-

19

 
 

2


   

0

 

 

-

1

 

   

2

 


 

 

-

2

 

 
   

17

 

   

11

 



Произведение обратной матрицы на исходную есть единичная  матрица, т.е. A-1 * A = Е, следовательно


X = A-1 * B


 

   

x1

 

=

   

x2

 

   

x3

 


 

   

1

 
 

2


 

-

15

 
 

2


   

25

 
 

2


*

 

-

1

 
 

2


   

13

 
 

2


 

-

19

 
 

2


   

0

 

 

-

1

 

   

2

 


 

 

-

2

 

 
   

17

 

   

11

 



x1 =

   

1

 
 

2


*

(

-

2

)


+

(

-

15

)

 

2


*

   

17

 

+

   

25

 
 

2


*

   

11

 

=

(

-

1

)


+

(

-

255

)

 

2


+

   

275

 
 

2


=

   

9

 


x2 =

(

-

1

)

 

2


*

(

-

2

)


+

   

13

 
 

2


*

   

17

 

+

(

-

19

)

 

2


*

   

11

 

=

   

1

 

+

   

221

 
 

2


+

(

-

209

)

 

2


=

   

7

 


x3 =

   

0

 

*

(

-

2

)


+

(

-

1

)


*

   

17

 

+

   

2

 

*

   

11

 

=

   

0

 

+

(

-

17

)


+

   

22

 

=

   

5

 


Ответ :


x1 = 9


x2 = 7


x3 = 5



 

 


Информация о работе Метод обратной матрицы