关于一个特殊矩阵的讨论
Discussion on a Special Matrix
DOI: 10.12677/pm.2024.146237, PDF, HTML, XML, 下载: 9  浏览: 38 
作者: 赵志辉, 王亚林, 方晓峰:火箭军工程大学基础部,陕西 西安
关键词: 矩阵向量矩阵的秩特征值与特征向量Matrix Vector Rank of Matrix Eigenvalues and Eigenvectors
摘要: 矩阵是线性代数中的重要表示工具和主要研究对象,不同结构的矩阵有着不同的价值和意义。本文主要从一个n维列向量和n维行向量的乘积获得的矩阵出发,通过典型例题剖析和重要理论分析相结合的方式,研究线性代数中重要理论之间的关联关系,为学生有效学习线性代数知识和教师设计线性代数的课堂教学提供帮助。
Abstract: Matrix is an important representation tool and a main research object in linear algebra, and matrices with different structures have different values and meanings. This paper mainly starts from the matrix of the product of multiplication of an n-dimensional column vector and an n-dimensional row vector. The paper studies the correlation between the important theories in linear algebra to provide help for students to effectively learn the knowledge of linear algebra and teachers to design the classroom teaching of linear algebra through the combination of typical example analysis and important theoretical analysis.
文章引用:赵志辉, 王亚林, 方晓峰. 关于一个特殊矩阵的讨论[J]. 理论数学, 2024, 14(6): 162-169. https://doi.org/10.12677/pm.2024.146237

1. 引言

线性代数是工科院校的一门重要基础课程,是研究多元线性问题理论工具,是工程数学中重要的组成部分,也是学生学习后续课程的数学基础。线性代数的特点是高度抽象、符号多、定义定理多,对初学者来说感觉有点繁杂。然而线性代数的知识点之间存在着密切的关系。在教学的过程中如果能够充分挖掘知识点之间的关联关系,往往可以达到事半功倍的效果。在矩阵构成的集合中,有一个特殊结构的矩阵在知识体系构建和解决问题等方面中具有重要价值。这个特殊的矩阵就是

( a 1 b 1 a 1 b 2 a 1 b n a 2 b 1 a 2 b 2 a 2 b n a n b 1 a n b 2 a n b n ) ,

其中 a i , b i ( i=1,2,,n ) ,这里特别要求 n2 ,因为 n=1 时,该矩阵退化为一个数,后面的讨论中将不再重复强调这一点。该矩阵可以建立与矩阵的运算、矩阵的秩、线性方程组、矩阵的特征值与特征向量、二次型五个重要知识之间的关联关系,为后续的教学提供有利的支撑。

2. 矩阵运算

2.1. 矩阵乘法及交换律

α= ( a 1 , a 2 ,, a n ) T β= ( b 1 , b 2 ,, b n ) T ,则

α β T =( a 1 a 2 a n )( b 1 , b 2 ,, b n )=( a 1 b 1 a 1 b 2 a 1 b n a 2 b 1 a 2 b 2 a 2 b n a n b 1 a n b 2 a n b n ) .

上述乘法运算既可以从矩阵乘法定义的角度体会矩阵乘法运算的规则,也可以从分块矩阵的视角去理解分块矩阵的乘法规则[1]。同时,作如下计算

β T α=( b 1 , b 2 ,, b n )( a 1 a 2 a n )= k=1 n a k b k =k ,

这一结果可看作向量数量积的矩阵乘法描述。结合上面两种方式的运算,不难发现下面的结论:

命题1 α= ( a 1 , a 2 ,, a n ) T β= ( b 1 , b 2 ,, b n ) T ,则

(1) α β T =O 的充要条件是 α=0 β=0

(2) β T α0 的必要条件是 α β T O

(3) α β T β T α

命题1证明比较简单,这里不再证明。注意结论(2)只是必要条件,而结论(3)说明矩阵乘法不满足交换律,这是矩阵运算中要特别注意到地方。

在矩阵乘法的教学中,应用此矩阵简单直观,可以很好的诠释矩阵运算的规则,以及不满足交换律的性质。

2.2. 矩阵的行列式与Binet-Cauchy定理

根据行列式的性质,容易计算 | α β T |=0 ,从这一特殊的结论出发Binet与Cauchy发现了乘积矩阵行列式的重要结论,即Binet-Cauchy定理[2]

Binet-Cauchy定理 A M m×n ( F ) B M n×m ( F ) ,则

| AB |={ 0,m>n, | A || B |,m=n, 1 k 1 < k 2 << k m m A( 1,2,,m k 1 , k 2 ,, k m )B( k 1 , k 2 ,, k m 1,2,,m ),m>n,

运用该结论(当 m>n 时, | AB |=0 )计算一些行列式更加快捷高效。

1 记矩阵

C=( sin2 α 1 sin( α 1 + α 2 ) sin( α 1 + α 3 ) sin( α 1 + α n ) sin( α 2 + α 1 ) sin2 α 2 sin( α 2 + α 3 ) sin( α 2 + α n ) sin( α 3 + α 1 ) sin( α 3 + α 2 ) sin2 α 3 sin( α 3 + α n ) sin( α n + α 1 ) sin( α n + α 2 ) sin( α n + α 3 ) sin2 α n ),n3 ,

计算 | C |

A=( sin α 1 cos α 1 sin α 2 cos α 2 sin α n cos α n ) B=( cos α 1 cos α 2 cos α n sin α 1 sin α 2 sin α n ) ,由三角函数和差角公式,有

sin α i cos α j +cos α i sin α j ={ sin2 α i ,i=j sin( α i + α j ),ij ,i,j=1,2,,n ,

于是 C=AB 。根据Binet-Cauchy定理有 | C |=| AB |=0

例1是北京科技大学2005年高等代数研究生入学考试题,此方法是计算该行列式最高效的方法。

在行列式的计算的教学中,引入这一结论,拓展计算某些特殊行列式计算方法,进一步提升学生的思维能力。

2.3. 乘法的结合律与方阵的高次幂

矩阵乘法虽然不满足交换律和消去律,但是矩阵乘法满足结合律,同时可以利用结合律计算矩阵的高次幂,而矩阵 α β T 就是最典型的实例。

命题2 α= ( a 1 , a 2 ,, a n ) T β= ( b 1 , b 2 ,, b n ) T A=α β T k= β T α ,则

(1) m + ,有 A m = k m1 A

(2) A 2 =O 的充要条件条件是 k=0

证明 (1) A m = ( α β T ) m =α ( β T α )( β T α ) ( m1 ) β T α β= ( β T α ) m1 ( αβ )= k m1 A

(2) “充分性” 当 k=0 时,由(1)的结果,有 A 2 =kA=O

“必要性” 当 A 2 =O 时,由(1)的结果,有 A 2 =kA=O 。假设 k0 ,则必有 A=O 。根据命题1的结论(1),知 α=0 β=0 ,于是 k= β T α=0 ,矛盾,从而 k=0

综合上述, A 2 =O 的充要条件条件是 k=0

命题2中结论(1)的证明过程充分体现了矩阵乘法的结合律,上述结论可进一步拓展到矩阵多项式。

在命题2条件的基础上,设关于x的多项式为 f( x )= c m x m + c m1 x m1 ++ c 1 x+ c 0 ,当

k0 时,有

f( A )= i=1 m c i k i1 A+ c 0 E .

2 已知 A=( 1 2 3 2 3 6 3 6 9 ) ,求 A n

α= ( 1,2,3 ) T β= ( 1,2,3 ) T ,则 β T α=14

由于 A=α β T ,则 A n = 14 99 ( 1 2 3 2 4 6 3 6 9 )

高次幂的计算方法很多,可以利用数学归纳法、二项展开式以及相似对角化等等。然而当一个矩阵可以分解为一个列向量和一个行向量的乘积时,利用矩阵乘法的结合律能够简化矩阵高次幂的计算,在矩阵乘法的教学,引入该矩阵,给学生扩展思维,提供计算矩阵高次幂的一种方法,更好的理解矩阵乘法不满足交换律,而满足结合律这一特性。

3. 矩阵的秩

由矩阵 α β T 的构成不难发现, α β T 的任一 k( 1kn ) 阶子矩阵都能分解成一个k维列向量与k维行向量的乘积,结合Binet-Cauchy定理容易获得 α β T 的所有2阶子式均为零,结合矩阵的秩的定义,即有 R( α β T )={ 1,α β T O 0,α β T =O 。对此结果作进一步拓展可得如下定理[3]

定理1 An阶非零矩阵,则 R( A )=1 的充要条件是存在n维非零列向量 α n维非零行向量 β T ,使得 A=α β T

由于非零列向量和非零行向量的秩均为1,因此把 A=α β T 也称为秩1分解。矩阵的低秩分解在数据压缩等方面有着重要应用。结合矩阵的初等变换还可得如下结论:

命题3 R( α β T )=1 时,有

(1) α β T r ( b 1 b 2 b n 0 0 0 0 0 0 ) ;(2) α β T ( 1 0 0 0 0 0 0 0 0 )

命题3的证明十分容易,但是它进一步揭示了矩阵的秩的含义、矩阵的行阶梯形及矩阵的标准形,对研究矩阵分类以及线性方程组求解有重要意义。

4. 线性方程组

根据命题3,可以得出关于线性方程组的相关结论:

命题4 R( α β T )=1 ,则齐次线性方程组 α β T x=0 β T x=0 同解[4]

命题4对研究 α β T x=0 的解有很大帮助。进一步研究发现,当 R( α β T )=1 时,必有 β T 0 。不妨设 b 1 0 ,可得 β T x=0 的一个基础解系为

ξ 1 =( b 2 b 1 0 0 ), ξ 2 =( b 3 0 b 1 0 ),, ξ n1 =( b n 0 0 b 1 ) ,

从而通解为 x= k 1 ξ 1 + k 2 ξ 2 ++ k n1 ξ n1 ,其中 k 1 , k 2 ,, k n1 为任意常数,即 α β T x=0 的解空间由它的一个基础解系张成。从几何角度看:当 n=2 时, α β T x=0 的解空间为过原点的直线 β T x=0 ;当 n=3 时, α β T x=0 的解空间为过原点的平面 β T x=0 ;当 n4 时, α β T x=0 的解空间为过原点的 ( n1 ) 维超平面 β T x=0 。解空间的正交补空间则由向量 β 张成,且它们的维数之和恰为n,即 ξ 1 , ξ 2 ,, ξ n1 ,β 线性无关,且可张成n维向量空间 n 。下面结合具体的示例说明这一点。

3 已知 α= ( 2,3,1 ) T β= ( 1,2,1 ) T ,求齐次线性方程组 α β T x=0 的通解。

根据上面的讨论, α β T x=0 的一个基础解系为

ξ 1 =( 2 1 0 ), ξ 2 =( 1 0 1 ) ,

通解为 x= k 1 ξ 1 + k 2 ξ 2 ,其中 k 1 , k 2 为任意常数。

Figure 1. The solution space of α β T x=0

1. α β T x=0 的解空间

图1中的平面 π 即为 α β T x=0 的解空间,以 β 为方向向量的直线为其正交子空间,同时3维向量空间 3 可由 ξ 1 , ξ 2 ,β 张成。进一步研究发现,向量 α α β T x=0 的一个基础解系有着如下的线性相关性。

命题5 R( α β T )=1 ,且 ξ 1 , ξ 2 ,, ξ n1 为齐次线性方程组 α β T x=0 的一个基础解系,则

(1) 当 β T α=0 时, ξ 1 , ξ 2 ,, ξ n1 ,α 线性相关;

(2) 当 β T α0 时, ξ 1 , ξ 2 ,, ξ n1 ,α 线性无关。

证明 (1) 当 β T α=0 时, α α β T x=0 的一个非零解,可由 ξ 1 , ξ 2 ,, ξ n1 线性表示,于是 ξ 1 , ξ 2 ,, ξ n1 ,α 线性相关;

(2) 当 β T α0 时,假设 ξ 1 , ξ 2 ,, ξ n1 ,α 线性相关。由 ξ 1 , ξ 2 ,, ξ n1 线性无关, α 可由 ξ 1 , ξ 2 ,, ξ n1 线性表示,即 α α β T x=0 的一个非零解,于是 α β T α=( β T α )α=0 。已知 β T α0 ,那么 α=0 ,从而 α β T =O ,与 R( α β T )=1 矛盾,故 ξ 1 , ξ 2 ,, ξ n1 ,α 线性无关。

命题5的直观理解可结合图1实现,这里不在赘述,这为后面相似对角化问题奠定直观的认知基础。

5. 矩阵的特值与特征向量

特征值与特征向量是线性代数中的重要理论,它是研究矩阵相似对角化等问题的理论基础[5]

5.1. 矩阵 α β T 的特征值、特征向量及特征多项式的降阶定理

任何矩阵的特征值及特征向量都与系数矩阵有着一定的联系,然而矩阵 α β T 的特征值及特征向量与向量 α,β 的联系更加密切,下面的命题可说明这一点。

命题6 λ 1 , λ 2 ,, λ n 为矩阵 α β T n个特征值,则

(1) 当 β T α=0 时,有 λ 1 = λ 2 == λ n =0 ;且 α β T O 时,对应着 ( n1 ) 个线性无关的特征向量 ξ 1 , ξ 2 ,, ξ n1 ,其中 ξ 1 , ξ 2 ,, ξ n1 α β T x=0 的一个基础解系;

(2) 当 β T α0 时,有 λ 1 = β T α λ 2 == λ n =0 ;且 α λ 1 对应的特征向量, ξ 1 , ξ 2 ,, ξ n1 λ 2 , λ 3 ,, λ n 对应 ( n1 ) 个线性无关的特征向量,其中 ξ 1 , ξ 2 ,, ξ n1 α β T x=0 的一个基础解系。

证明 (1) 当 β T α=0 时,由命题2可知 ( α β T ) 2 =O ,则 ( α β T ) 2 的特征值均为零,于是 α β T 的特征值也全为零;

(2) 当 β T α0 时,由 α β T α=α( β T α )=( β T α )α 及特征值与特征向量定义可知, β T α α β T 的非零特征值,且 α 为对应的特征向量。再由 | α β T |=0 ,可知0也是 α β T 的特征值。又 R( α β T )=1 ,所以特征值0的几何重数是 n1 。根据特征值的几何重数小于等于代数重数,可得0是 α β T n1 重的特征值,于是 β T α α β T 的单特征值,从而 λ 1 = β T α λ 2 == λ n =0

(1)、(2)中关于特征值0对应的特征向量可通过线性方程组部分的结论即可证明。

命题6的结论(2)表明,矩阵 α β T β T α 有相同的非零特征值,结论(1)和(2)均表明 α β T β T α 有相同的迹。沿着此结论作深入思考,将向量 α,β 用一般矩阵AB替代,可得特征多项式的降阶定理。

特征多项式降阶定理AB分别为 m×n n×m 矩阵, mn ,则

| ABλ E m |= ( 1 ) mn λ mn | BAλ E m | .

特征多项式降阶定理的证明这里不再赘述,定理表明ABBA有相同的非零特征值,亦有相同的迹。

5.2. 矩阵 α β T 的相似对角化

矩阵的相似对角化取决于矩阵的特征值与特征向量,根据上面讨论易得如下命题[6]

命题7 R( α β T )=1 ,则

(1) 当 β T α=0 时,矩阵 α β T 不能相似对角化;

(2) 当 β T α0 时,存在可逆矩阵P,使得 P 1 AP=Λ ,其中

P=( a 1 b 2 b n1 b n a 2 b 1 0 0 a 3 0 0 0 a n1 0 b 1 0 a n 0 0 b 1 ),Λ=( β T α 0 0 0 ) .

命题7的证明通过命题6即可完成。从命题7中亦可看出矩阵的秩与非零特征值的个数(重数计入个数)并不相同,但矩阵若能相似对角化,其非零特征值个数必等于秩。

特别地,当 β=α0 时,矩阵 α α T 是实对称矩阵,则必可以对角化。同时 α α T 可作为某n元二次型的矩阵。

5.3. 矩阵 α α T 与二次型

矩阵 α α T 对应的二次型为 f( x 1 , x 2 ,, x n )= X T ( α α T )X ,存在正交变换 X=QY ,使得

f( y 1 , y 2 ,, y n )=( α T α ) y 1 2 ,

其中正交矩阵Q可以通过命题7中可逆矩阵P正交单位化得到。

4 α β 都是第一个分量非零的 n( n>2 ) 维列向量,且 α T β=3 ,求矩阵 B=α β T 的特征值与特征向量。

x=( a 1 , a 2 ,, a n ) , y=( b 1 , b 2 ,, b n ) ,且 a 1 0 b 1 0 x T y= i=1 n a i b i x y T =( a 1 b 1 a 1 b 2 a 1 b n a 2 b 1 a 2 b 2 a 2 b n a n b 1 a n b 2 a n b n )

λ 1 , λ 2 ,, λ n B=x y T 的特征值,则 i=1 n λ i = x T y=3 。因为 x,y 都是第一个分量非零,所以 R( x y T )=1 ,由命题6可得, λ 1 =3 λ 2 == λ n =0

因为 ( x y T )x=x( y T x )=3x ,所以 k 1 x λ 1 =3 时的特征向量, ( x y T 0E )=( x 1 y 1 x 1 y 2 x 1 y n x 2 y 1 x 2 y 2 x 2 y n x n y 1 x n y 2 x n y n )( 1 y 2 y 1 y n y 1 0 0 0 0 0 0 )

所以 λ 2 == λ n =0 对应的特征向量为 ξ 2 = k 2 ( y 2 y 1 ,1,0,,0 ),, ξ n = k n ( y n y 1 ,0,0,,1 ) 其中 k i 0( i=1,2,,n )

这道题的解法的关键,首先证明B是秩1矩阵,其次要证明 x T y=3 是矩阵B唯一的非零解最后要

证明矩阵B可以相似对角化。

6. 结束语

这个特殊矩阵在线性代数的学习和教学中占有重要的地位,从这一特殊矩阵出发,将其分解为一个列向量与行向量的乘积,由其运算可以深刻理解矩阵的乘法运算,验证矩阵乘法不满足交换律这一特性;由其行列式的值,可以得到了著名的Binet-Cauchy定理;由其秩揭示了矩阵的秩的含义、矩阵的行阶梯形及矩阵的标准形,进一步可以探讨线性方程组的求解、解的结构;由其特征值和特征向量,深刻理解相似对角化理论和二次型理论。这一特殊结构的矩阵将线性代数的重要理论都关联起来,不仅形成脉络清晰的知识网络,而且为许多问题的解决提供了一种简便的思路。通过这一特殊的矩阵,便于学生将所学的知识点串联起来,形成整体的概念,内化于心,达到融会贯通的目的。

参考文献

[1] 同济大学数学系. 线性代数(第六版) [M]. 北京: 高等教育出版社, 2014.
[2] 徐数方, 等. 数值线性代数(第二版) [M]. 北京: 北京大学出版社, 2013.
[3] 斯蒂芬∙博伊德, 利芬∙范登伯格. 应用线性代数: 向量、矩阵及二乘[M]. 北京: 机械工业出版社, 2020.
[4] 曾诚, 迟楠, 吕井明. 新工科时代应用类本科院校“线性代数”课程的体系重构[J]. 遵义师范学院学报, 2022, 24(5): 115-118.
[5] 杨威, 高淑萍, 李兵斌. 论矩阵等式在线性代数中的重要性[J]. 大学教育, 2021(3): 116-119.
[6] 项巧敏. 线性代数中矩阵的秩与向量组的秩的教学方法探讨[J]. 新教育时代电子杂志(教师版), 2021(19): 190-191.