浅谈几种求逆矩阵的方法
Discussion on Several Methods of Finding Inverse Matrix
摘要: 矩阵求逆作为线性代数的一个重要研究内容,在通信工程、计算机、物理学等领域有着广泛的应用。矩阵是大学数学中很重要的一个内容,在《高等代数》中我们学习了矩阵的一些基本知识及应用,而矩阵求逆的方法是矩阵中一个很重要的部分,那么如何判断一个矩阵是否可逆,怎样快速的去求解矩阵的逆,前人也总结了一些非常实用的方法。本文结合自身所掌握的知识,结合一些有代表性的例子进行说明,为了更便捷地解决求矩阵的逆,本文根据不同矩阵的不同特点简单介绍了几种求逆矩阵的方法:定义法、伴随矩阵法、初等变换法、分块矩阵法,同时也找出了一类特殊矩阵求逆的方法。
Abstract: As an important research content of linear algebra, matrix inverse is widely used in communication engineering, computer, physics and other fields. Matrix is an important content in university mathematics, in the “advanced algebra”, we learned some basic knowledge and application of matrix, the matrix inversion method is a very important part of the matrix, so how to determine whether a matrix is reversible, how to quickly to solve the inverse of the matrix, predecessors also summarize some very practical method. The text is explained in combination with its own knowledge and some representative examples. In order to solve the inverse matrix more conveniently, this paper simply introduces several methods of inverse matrix according to the different characteristics of different matrices: definition method, adjoint matrix method, elementary transformation method, partition matrix method. At the same time, the inverse method of a special kind of matrix is found.
文章引用:徐晓静. 浅谈几种求逆矩阵的方法[J]. 应用数学进展, 2021, 10(9): 3217-3224. https://doi.org/10.12677/AAM.2021.109336

1. 逆矩阵的概念 [1]

n阶方阵A是可逆的,如果存在数域P上的n阶方阵B,使得 A B = B A = E ,这里的E是n阶单位阵。当A可逆时,逆矩阵由A唯一确定,记为 A 1

2. 利用初等变换求矩阵的逆 [1]

初等变换法:

1. ( A | E ) ( E | A 1 )

对于阶数 ( n 3 ) 的矩阵一般采用初等行变换,但需要注意用上述方法求矩阵的逆只可以用初等行变换。

2. 可以利用 ( A E ) ( E A 1 )

3. 特别的当A是可逆的,可以用 ( A | B ) ( E | A 1 B ) ( A C ) ( E C A 1 )

快速求出 A 1 B C A 1

例1:

求矩阵 A = ( 2 2 3 1 1 0 1 2 1 ) 的逆矩阵 A 1

解:利用 ( A | E ) ( E | A 1 )

( A | E ) = ( 2 2 3 1 1 0 1 2 1 | 1 0 0 0 1 0 0 0 1 ) r 2 r 1 ( 1 1 0 2 2 3 1 2 1 | 0 1 0 1 0 0 0 0 1 ) r 3 + r 1 r 2 2 r 1 ( 1 1 0 0 4 3 0 1 1 | 0 1 0 1 2 0 0 1 1 ) r 3 r 2 ( 1 1 0 0 1 1 0 4 3 | 0 1 0 0 1 1 1 2 0 ) r 3 4 r 2 r 1 + r 2 ( 1 0 1 0 1 1 0 0 1 | 0 2 1 0 1 1 1 6 4 ) r 3 × ( 1 ) ( 1 0 1 0 1 1 0 0 1 | 0 2 1 0 1 1 1 6 4 ) r 1 r 3 r 2 r 3 ( 1 0 0 0 1 0 0 0 1 | 1 4 3 1 5 3 1 6 4 ) = ( E | A 1 )

A 1 = ( 1 4 3 1 5 3 1 6 4 )

例2 求矩阵 A = ( 1 1 1 2 1 0 1 1 0 ) 的逆矩阵 A 1

解利用 ( A E ) ( E A 1 )

( 1 1 1 2 1 0 1 1 0 1 0 0 0 1 0 0 0 1 ) r 3 + r 1 ( 1 1 0 2 1 2 1 1 1 1 0 1 0 1 0 0 0 1 ) r 2 r 1 ( 1 0 0 2 1 2 1 2 1 1 1 1 0 1 0 0 0 1 ) r 1 r 3 ( 1 0 0 0 1 2 0 2 1 0 1 1 0 1 0 1 0 1 ) r 3 + 2 r 2 ( 1 0 0 0 1 0 0 2 3 0 1 1 0 1 2 1 0 1 ) r 2 + ( 2 3 ) r 3 ( 1 0 0 0 1 0 0 0 3 0 1 3 1 0 1 3 2 1 2 3 1 ) r 2 × ( 1 ) r 3 × ( 1 3 ) ( 1 0 0 0 1 0 0 0 1 0 1 3 1 3 0 1 3 2 3 1 2 3 1 3 ) = ( E A 1 )

A 1 = ( 0 1 3 1 3 0 1 3 2 3 1 2 3 1 3 )

3. 利用伴随矩阵求矩阵的逆 [2]

伴随矩阵法: A 1 = 1 | A | A *

对于阶数比较低( n 3 )的矩阵或元素的代数余子式易于计算的矩阵可用此方法求逆矩阵。特别的对

于2阶方阵 A = ( a 11 a 12 a 21 a 22 ) 的伴随矩阵为 A * = ( a 22 a 12 a 21 a 11 ) ,即可以概括为主对角元素互换,次对角元素

变号。

例3 求矩阵 A = ( 1 2 3 4 5 8 3 4 6 ) 的逆矩阵 A 1

解: | A | = | 1 2 3 4 5 8 3 4 6 | = 1

各个元素的代数余子式分别为: A 11 = 2 A 12 = 0 A 13 = 1 A 21 = 0 A 22 = 3 A 23 = 2 A 31 = 1

A 32 = 4 A 33 = 3 ,故有伴随矩阵 A * = ( 2 0 1 0 3 4 1 2 3 ) ,有公式 A 1 = 1 | A | A * ,故

A 1 = ( 2 0 1 0 3 4 1 2 3 )

4. 利用哈密顿–凯莱定理求逆矩阵 [3]

定理:设A是数域P上的n阶方阵,则 f ( λ ) = | λ E A | 是A的特征多项式, f ( A ) = A n ( a 11 + a 22 + + a n n ) A n 1 + + ( 1 ) n A E = 0

f ( A ) = A n + a 1 A n 1 + + a n 1 A + a n E = 0

可以得到 A 1 = 1 a n ( A n 1 + a 1 A n 2 + + a n 1 E )

注:此种方法适用于阶数不多的矩阵求逆。

例4 求矩阵 A = ( 1 1 1 1 1 0 2 2 1 ) 的逆矩阵 A 1

解:A的特征多项式为 | λ E A | = ( λ 1 1 1 1 λ 1 0 2 2 λ 1 ) = λ 3 3 λ 2 + 2 λ 2 ,由哈密顿–凯莱定理有

f ( A ) = A 3 3 A 2 + 2 A 2 E = 0 ,所以得到

A 1 = A 2 3 A + 2 E 2 = 1 2 ( 1 3 1 1 1 1 0 4 2 )

5. 利用分块矩阵求矩阵的逆

分块矩阵多适用于高阶矩阵求逆。

对于分块的对角矩阵或次对角交矩阵可以有公式一、公式二成立 [4];

公式一 ( A 1 A n ) 1 = ( A 1 1 A n 1 ) ( A 1 A n ) 1 = ( A n 1 A 1 1 ) ,其中每个 A i ( i = 1 , 2 , )

均为可逆矩阵。

例5 已知 A = ( 0 0 5 2 0 0 2 1 1 2 0 0 1 1 0 0 ) ,求 A 1

将A分块如下 = ( 0 0 5 2 0 0 2 1 1 2 0 0 1 1 0 0 ) = ( 0 A 1 A 2 0 ) ,其中, A 1 = ( 5 2 2 1 ) A 2 = ( 1 2 1 1 ) ,对于二阶矩阵可求得, A 1 1 = 1 | A 1 | A 1 * = ( 1 2 2 5 ) A 2 1 = 1 | A 2 | A 2 * = 1 3 ( 1 2 1 1 ) ,所以有 A 1 = ( 0 A 2 1 A 1 1 0 ) = ( 0 0 1 3 2 3 0 0 1 3 1 3 1 2 0 0 2 5 0 0 )

可见对于矩阵分块可以简化高阶矩阵求逆,将其转化为容易求逆的低阶矩阵,进而在求逆。

设矩阵 X = ( A B C D ) ,其中ABCD均为方阵。

方法一1

则可以利用初等变换

( A B C D | E E ) r 1 + B D 1 r 2 ( A B D 1 C 0 C D | E E B D 1 E ) r 2 + C ( A B D 1 C ) 1 E r 1 ( A B D 1 C 0 0 D | E B D 1 C ( A B D 1 C ) 1 E E + C ( A B D 1 C ) 1 B D 1 ) r 1 × ( A B D 1 C ) 1 r 2 × D 1 ( E 0 0 E | ( A B D 1 C ) 1 E ( A B D 1 C ) 1 B D 1 D 1 C ( A B D 1 C ) 1 E D 1 + D 1 C ( A B D 1 C ) 1 B D 1 )

方法二:利用 ( A B C D ) 进行初等行变换或者初等列变换:

1. 进行行初等变换得到 P 1 P 2 P n ( A B C D ) = ( E E ) ( P 1 P 2 P n ) 为初等矩阵

2. 进行列初等变换得到 ( A B C D ) Q 1 Q 2 Q n = ( E E ) ( Q 1 Q 2 Q n ) 为初等矩阵

( A B C D ) 1 = P 1 P 2 P n = Q 1 Q 2 Q n

3. 或者同时进行行列的初等变换

P 1 P 2 P s ( A B C D ) Q 1 Q 2 Q t = ( A 1 A n ) ,其中每个 A i ( i = 1 , 2 , ) 均为可逆的对角矩阵, ( A B C D ) = P s 1 P 2 1 P 1 1 ( A 1 A n ) Q 1 1 Q 2 1 Q t 1

( A B C D ) 1 = Q 1 Q 2 Q t ( A 1 A n ) 1 P 1 P 2 P s

利用上面的分块矩阵的初等变换的步骤很容易算出, ( A B O D ) , ( A 0 C D ) , ( O B C D ) , ( A B C O ) ,这四类分

块矩阵的逆矩阵了。

例6 已知 A = ( 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ) ,求 A 1

A = ( 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ) = ( B B B B ) ,其中 B = ( 1 1 1 1 ) ,则有 B 1 = 1 2 ( 1 1 1 1 ) = 1 2 B

利用方法一:

( B B B B | E E ) r 1 + r 2 ( 2 B 0 B B | E E E ) r 1 × 1 2 B 1 ( E 0 B B | 1 2 B 1 1 2 B 1 E ) r 2 × B 1 ( E 0 E E | 1 2 B 1 1 2 B 1 B 1 ) r 2 r 1 ( E 0 0 E | 1 2 B 1 1 2 B 1 1 2 B 1 1 2 B 1 ) r 2 × ( 1 ) ( E 0 0 E | 1 2 B 1 1 2 B 1 1 2 B 1 1 2 B 1 ) = ( E | A 1 ) A 1 = 1 4 ( B B B B ) = 1 4 A

利用方法二:

( B 1 0 0 B 1 ) ( E 0 0 1 2 E ) ( E 1 2 E 0 E ) ( E 0 E E ) ( B B B B ) = ( E 0 0 E )

( B B B B ) 1 = ( B 1 0 0 B 1 ) ( E 0 0 1 2 E ) ( E 1 2 E 0 E ) ( E 0 E E ) 1 4 ( B B B B ) = 1 4 A

6. 利用初等循环矩阵进行求逆

A为 n × n 阶循环矩阵,如果存在多项式, f ( x ) = i = 1 i = n a i x i 满足 f ( P ) = i = 1 i = n a i P i ,称多项式 f ( x ) 为A的关联多项式,这里的系数 a i ( i = 1 , 2 , , n ) 正好是第一行元素。

初等循环矩阵定义:A的行和为 s = i = 1 i = n a i ,则

(i) A E n = s E n

(ii) 当A可逆,且 s 0 ,则有 A 1 E n = 1 s E n

证明: A E n = ( i = 1 i = n a i P i ) E n = i = 1 i = n ( a i P i E n ) = s E n ,两边取逆 E n 1 A 1 = 1 s E n 1

定理 [5]:设A为 n × n 阶循环矩阵且可逆,则可以证明下列结论成立

(i) ( A + c E n ) 1 = A 1 c s ( s + c n ) E n

(ii) 证明: ( A + c E n ) 1 ( A + c E n ) = I ,所以有 ( A + c E n ) 1 A = I c ( A + c E n ) 1 ,A的逆右乘两端得:

( A + c E n ) 1 = A 1 c s c ( A + c E n ) 1 E n ,化简移项得: ( A + c E n ) 1 ( A E n + c E n 2 ) E n = E n ,所以有 ( A + c E n ) 1 E n = A 1 c s ( s + c n ) E n 成立。

结论:通过(i)我们知道,如果可用 A + c E n 来表示一个初等循环矩阵B,其中A是一个初等循环矩阵并且A的逆矩阵易求。

例7设 B = ( 7 11 7 7 7 7 7 11 11 7 7 7 7 7 11 7 ) ,求 B 1

解: B = 4 ( 0 1 0 0 0 0 0 1 1 0 0 0 0 0 1 0 ) + 7 ( 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ) ,设 A = 4 ( 0 1 0 0 0 0 0 1 1 0 0 0 0 0 1 0 ) s = 4 c = 7 n = 4

则有 B 1 = A 1 c s ( s + c n ) E 4 = 1 128 ( 7 7 25 7 25 7 7 7 7 7 7 25 7 25 7 7 )

推论:通过构造初等循环矩阵,利用一个公式就能轻松地求出一个较为复杂的矩阵的逆。同时我们还可以与上面的分块矩阵求逆相结合,把一个求高阶的循环矩阵变成求低阶的循环矩阵。

7. 总结

本篇文章利用矩阵的初等变换,分块矩阵的初等变换以及哈密顿–凯莱定理对矩阵求逆的方法,进行了总结分析。以上求逆矩阵的方法只是一些初步认识,通过归纳总结可以看出,在求逆矩阵的过程中,必须熟练掌握运用相关的可逆矩阵定理和性质,同时运用恰当的技巧可使计算更加简便。但是文章最后我们用了初等循环矩阵求逆,这是针对一类特殊类型的矩阵求逆,所以我们能看到矩阵求逆的方法多种多样,对于其他类型的矩阵求逆,则需要继续深入研究。

NOTES

1李扬高等代数强化讲义[Z]. 2021。

参考文献

[1] 徐仲. 高等代数考研教案[M]. 西安: 西北工业大学出版社, 2006.
[2] 王萼芳, 石生明. 高等代数[M]. 第5版. 北京: 高等教育出版社, 2019.
[3] 张家宝. 浅谈求逆矩阵的几种方法[J]. 数学学习与研究, 2020(10): 4-5.
[4] 钱吉林. 高等代数题解精粹[M]. 修订版. 北京: 中央民族大学出版社, 2005.
[5] 朱亚培, 薛娇. 关于循环矩阵求逆方法的研究[J]. 现代经济信息, 2017(12): 388+391.