矩阵对角化及应用论文

1、丽水学院2012届学生毕业论文矩阵对角化及应用理学院 数学082 缪仁东 指导师：陈巧云摘 要:本文是关于矩阵对角化问题的初步研究,对矩阵对角化充要条件的归纳,总结,通过对实对称矩阵,循环矩阵,特殊矩阵对角化方法的计算和研究,让读者对矩阵对角化问题中求特征值、特征向量,求可逆矩阵,使对角化,提供了简便,快捷的求解途征.关键词:对角矩阵;矩阵对角化;实对称矩阵;特征值;特征向量.矩阵对角化是矩阵论的重要组成部分,在矩阵论中占有重要的作用,研究矩阵对角化问题很有实用价值,关于矩阵对角化问题的研究,这方面的资料和理论已经很多.但是他们研究的角度和方法只是某个方面的研究,没有进行系统的分类归纳和总结.

2、因此,我就针对这方面进行系统的分类归纳和总结,对一些理论进行应用和举例,给出算法.特别给出了解题时方法的选择.1矩阵对角化概念及其判定所有非主对角线元素全等于零的n阶矩阵,称为对角矩阵或称为对角方阵.定义1.1 矩阵是数域上的一个n级方阵. 如果存在一个上的n级可逆矩阵,使 为对角矩阵,则称矩阵可对角化. 矩阵能否对角化与矩阵的特征值特征向量密切相关.定义 1.2 设是一个阶方阵,是一个数,如果方程组                 &

3、#160;                   (1)存在非零解向量,则称为的一个特征值,相应的非零解向量称为属于特征值的特征向量（1）式也可写成,                        &#

4、160;           (2)这是个未知数个方程的齐次线性方程组,它有非零解的充分必要条件是系数行列式                         ,       

5、60;                             (3) 即                     &#

6、160;           上式是以为未知数的一元次方程,称为方阵的特征方程 其左端是的次多项式,记作,称为方阵的特征多项式                 显然,的特征值就是特征方程的解特征方程在复数范围内恒有解,其个数为方程的次数（重根按重数计算）,因此,阶矩阵有个特征值设阶矩阵的特征值为,由多项式的根与系数之间的关系,不难证明（）;（）.若

7、为 的一个特征值,则一定是方程的根, 因此又称特征根,若为方程的重根,则称为的重特征根方程 的每一个非零解向量都是相应于的特征向量,于是我们可以得到求矩阵的全部特征值和特征向量的方法如下：     第一步：计算的特征多项式；     第二步：求出特征方程的全部根,即为的全部特征值；     第三步：对于的每一个特征值,求出齐次线性方程组：           

8、0;          的一个基础解系,则的属于特征值的全部特征向量是          （其中是不全为零的任意实数） 设是数域, Mn () 是上n ×n 矩阵构成的线性空间, A Mn () , 为A 的t 个互不相同的特征值,高等代数第二版（北京大学数学系几何与代数教研室编）第四版（张和瑞、郝炳新编）课程中,我们学过了矩阵可对角化的若干充要条件如：(1) A 可对角化当且仅当A 有n 个线性无关的特征

9、向量;(2) A 可对角化当且仅当特征子空间维数之和为n ;(3) A 可对角化当且仅当A 的初等因子是一次的;(4) A 可对角化当且仅当A 的最小多项式无重根 我们知道线性变换A 的特征多项式为f () ,它可分解成一次因式的乘积则V 可分解成不变子空间的直和其中 = | ;V引理 1.1:设A, B 都是n 阶矩阵, 则秩( AB) 秩( A) + 秩( B) - n.定理 1.1:设A 是实数域F上的一个n 阶矩阵, A 的特征根全在F 内, 若, ,., 是A 的全部不同的特征根, 其重数分别为, ,. , 那么 () 可对角化的充要条件是 j=1, 2,.k () 当( 1) 式成

10、立时, 的列空间就是A 的属于特征根的特征子子空间.证明: () 设A 可对角化, 则存在可逆阵T, 使这里右边是分块对角矩阵, 为阶单位阵, 于是有 = = = j=1,2, .k.反之,若 i=1,2,.k, 反复用引理可得 j=1,2,.,k.这里用到了齐次线性方程组的解空间的维数不大于的重数不大于这个结论.于是又 从而 i=1,2,.k. 这样的矩阵可以对角化.()设( )式成立,则A 可对角化.故A的最小多项式为从而 即 这就是说,列空间包含在的特征子空间中,但是由(1), 的列空间的维数是n,它正是的特征子空间的维数,所以结论() 成立.推论: 设A 为实数域F上的n阶矩阵,A 的

11、特征根全为F 内,且, 是A的全部不同的特征根, 其维数分别为, , 若秩,秩,则A 可以对角化,且的列向量组的极大无关组恰是属于 的极大线性无关的特征向量组,的列向量组的极大无关组恰是属于的极大无关的特征向量组.例1: 判断A=能否对角化,并求特征向量.解: 易知A的特征根 =-2 , =1. = 和 =的秩分别为2与1,故A可对角化. 又因为可以选取和为的列空间的一个基,是属于的特征向量. 定理和推论把判断矩阵是否对角化的问题与求它的特征向量的问题联系起来,给出了一个不用解线性方程组而求得可对角化矩阵的特征向量的方法, 在矩阵的不同特征根较少时, 这个方法较方便.2实对称矩阵对角化的计算方

12、法我们知道任意实对称矩阵,总正交相似于一对角阵. 该对角阵的对角元即为实对称矩阵的特征值, 正交相似变换矩阵的各列构成相应的特征向量. 给定一实对称阵A ,如何求正交相似变换矩阵P ,使为对角阵. 理论上的解决方法为:首先利用特征方程: | I - A | = 0 求出全部特征值,针对不同特征值求出相应的完全特征向量系,合在一起构成实对称阵A 的完全特征向量系. 再利用施密特正交化法得到A 的规范化正交特征向量系. 以此作为列向量得到正交相似变换矩阵P , 为对角阵, 参见文献5 . 此方法理论可行,但在具体操作时,由于要事先求出实对称阵A 的全部特征值,操作上有如下困难: (1) 特征方程:

13、 | I- A | = 0 给出困难; (2) 特征方程求根困难(5 次以上的代数方程没有统一的求根公式) . 因此有必要寻求方法.定义2.1 (瑞雷商) 设A 为n 阶实对称阵,对于任一n维非零列向量x ,称R ( x) =( A x , x)/( x , x) 为关于向量x 的瑞雷商.引理2.1 设A 为n 阶实对称阵, . 为A 的特征值. 定义2.2设w 为n 维列向量,且= 1 ,则n 阶矩阵H = I - 2 称为Householder 阵.引理2.2Householder 矩阵具有如下性质:(1) (2) ( H 是正交阵) .引理2.3设x , y , x y , ,则存在Ho

14、useholder 矩阵H, 使Hx = y. 其中定理2.1 设A 是实对称矩阵, x ( = 1) 是A 的一个特征值和相应的特征向量,则存在P 为一个正交阵,使Px = = . 且的第一行和第一列的第一个元素为,其余元素均为零.证 设A 是实对称矩阵, .为A 的特征值. 根据引理2.1 ,利用多元函数求极值的拉格朗日乘数法,可求得 及相应的规范化特征向量 . 不妨假设 = 1 ,由引理2.3 ,存在为一个正交阵,使=.且的第一行和第一列的第一个元素为 , 其余元素均为零. 设, 为对称阵,故 也为对称阵,设 及 为最大特征值及相应的规范化特征向量,则根据引理2.3 ,存在为一个正交阵,

15、使.且 的第一行和第一列除 外其余元素均为零. 令,容易验证亦为正交阵, 满足:依此类推, 存在正交阵, , , 使得.,则=D, 其中D 为对角阵,令,则,即为将实对称阵对角化的正交相似变换矩阵.例2: 设矩阵, 为A的特征值.按上面的算法进行对角化,求出正交矩阵P及特征根和特征向量. 解: (1)利用瑞雷商和多元函数求极值的拉格朗日乘数法,可求得 = 18 ,相应的特征向量为(2) 计算正交矩 =,满足且,至此已实现对角化. 借此可求得= =9 , = - 9. 相应的特征向量分别为,.3循环矩阵对角化方法的研究在复数域C 上,形如的矩阵,称关于元素列的循环矩阵.已知n阶循环矩阵,并令 ,

16、称为循环矩阵基本列(其中= 为单位矩阵).循环矩阵基本列有如下特点:都是循环矩阵; ,即;n 阶循环矩阵K有n 个特征根: 关于元素列的n 阶循环矩阵A 可用循环矩阵基本列表示为,反之,能用循环矩阵基本列线性表示的矩阵,则一定是循环矩阵.循环矩阵的性质 性质1 同阶循环矩阵的和矩阵为循环矩阵. 性质2 同阶循环矩阵的乘积满足交换律. 性质3 同阶循环矩阵的乘积为循环矩阵. 性质4 循环矩阵的逆矩阵为循环矩阵.n阶矩阵A 关于多项式函数f (x) 生成的矩阵为f (A) ,A 的特征根与f (A)的特征根有下面的结论:命题3.1 设f (x) 是一个n - 1 次多项式函数,若是矩阵A 的特征根

17、,则f () 是矩阵f (A) 的特征根.命题3.2 设f (x) 是一个n - 1 次多项式函数,若矩阵A 相似于矩阵B ,则矩f (A) 相似于矩阵f (B) .考察n 阶循环矩阵K,K的特征多项式为:如果n 阶循环矩阵A 记为不难求得K中与特征值相应的特征向量,记:, 则由命题3.1 得,可以验证.将这n 个两两正交的向量单位化,可得标准正交基,令矩阵则命题3.3任意n 阶循环矩阵 在复数域C 上都可对角化,即 =推论n 阶循环矩阵A 可逆的充要条件是(i=0,1,.,n-1).例3:求四阶循环矩阵的特征根,并对角化.解: 令 得 ,由于, 所以A的特征根分别为:=10 , =-2-2i

18、, =-2, =-2+2i, 4特殊矩阵特殊对角化的研究前面对实对称矩阵循环矩阵的对角化问题作了研究,本部分主要讨论,当矩阵只有两个特征根时的对角化问题,方法简捷. 对于数域F 上的n 阶矩阵A ,若仅有的两个特征根都在F 内,并且可以对角化,不通过解线性方程组求特征向量,而用初等变换求出可逆矩阵T,使为对角形矩阵.定理4.1 设数域F 上的n 阶矩阵A 可以对角化,其特征根为,如果 P,B 为列满秩矩阵,那么(i) A 的属于 的线性无关的特征向量为 的个列向量;A 的属于的线性无关的特征向量为B 的个列向量.(ii) 令T = ( P ,B) ,则T 可逆,且有其中 有个, 有个.证 因为

19、初等矩阵不改变矩阵的秩,且为列满秩,则的重数.（i）根据矩阵的初等变换和分块矩阵的运算性质,可得 从而因为列满秩矩阵,则的个列向量为齐次线性方程组的基础解系,亦即的 个列向量为的属于的线性无关的特征向量.又可以对角化,且的重数为 ,则有可逆矩阵Q,使得, 令,则有 = = 由于 的列向量为 的列空间的基,则的 个列向量为齐次线性方程组的基础解系, 的 个列向量为 的属于的线性无关的特征向量.(ii) 因矩阵的属于不同特征根的特征向量线性无关,且特征向量的个数之和等于的阶数,于是, 令 即有例4:令矩阵,求可逆矩阵T,使得为对角形式.解: 方法一,先求A的特征根= 则 = 1 (二重) , =

20、- 1. 可见,此例为定理所述的情况.对矩阵作初等列变换,即所以,由定理4.1 知,A 的属于 = - 1 的线性无关的特征向量为;A 的属于 = 1 的线性无关的特征向量为 , 令,则有. 这与1 的结果一致.方法二 在矩阵中,亦可取,这时则A 的属于=1 的线性无关的特征向量为 , ;A 的属于=- 1 的线性无关的特征向量为令,则有.5常规矩阵对角化方法的新探 众所周知,对数域上一个阶矩阵是否存在一个可逆矩阵,使得为对角形矩阵,当这种矩阵存在时,如何去寻求它.一般有关教材中都是先计算一个行列式,求出的特征值,再利用线性方程组和特征向量的有关理论及求法解决此问题的.在这里利用矩阵的初等变换

21、解决此问题的,它比教材中的常规方法简单一些,因为不必解若干的齐次线性方程组,有时也不必计算行列式.5.1理论依据为说话方便,我们规定如果数域上,对阶矩阵存在一个可逆矩,使得为对角形矩阵, 则称矩阵在数域上可对角化.当可对角化时, 我们说将A对角化,即指求矩阵,使为对角形矩阵.若矩阵在数域上可对角化, 则有上可逆矩阵,使得为对角形矩阵.于是的主对角线上的元素,即为的全体特征值, 并且可表示： 其中为初等矩阵,i=1,2,.,s,于是,又也是初等矩阵, 由初等矩阵与矩阵的初等变换的关系, 即知 , 相当于对施行了一次初等行变换与一次初等列变换.这里, 我们称此种初等变换为对施行了一次相似变换.显见

22、, 可对施行一系列的相似变换化为.又由, （E此处表单位矩阵）可如下进行初等变换, 则可将化为对角形矩阵, 且可求得：,对E只施行其中的初等列变换.当A不可对角化时, 也可经相似变换化简A后, 求得其特征值, 判定它可否对角化.类似地, 可由,做如下初等变换则可将A化为对角形矩阵B,且可求得T或由B求A的特征值, 判定可否对角化：,对E只施行其中的初等行变换.并且在施行相似变换时, 不必施行一次行变换后接着施行一次列变换这样进行, 可施行若干次行或列变换后再施行若干次相应的列或行变换, 只要保持变换后, 最后所得矩阵与相似即可.5.2 应用举例为叙述简便,这里用表示i第行,表示第i列,表示用数

23、k乘第j行后再加到第i行上,表示用数k乘第j列后再加到第i列上.例5 求如下矩阵的特征值, 并判定它们可否对角化,若可则将其对角化:(1), (2). 解：（1）由 ,知A与C相似. 易得,C的特征值为2,2,2,且2E-C的秩为2,所以C不能对角化,从而知A的特征值为2,2,2且A不可以对角化. （2）由 ,知B可以对角化,B的特征值为-2,2,2,2.令, 则.当不易直接用相似变换化简判定时, 可先求出特征值, 再用相似变换.例6判定可否对角化,若可,则将其对角化.解法1（教材中的方法）由 ,知A的特征值为4,6,-1,-1.解 齐次线性方程组得一基础解系解 齐次线性方程组得一基础解系解

24、齐次线性方程组得一基础解系,于是可,A可对角化,且取,则.解法2由 知,A可对角化,且取.,两法比较, 法2比法1简便, 因不必计算行列式和解几个线性方程组.上述内容为本人对各类基本常见的矩阵类型的对角化计算方法,计算技巧的一些探讨,比较传统的计算方法、计算技巧,有一些优越性.计算简便,步骤简单具体,有较强的实用性.参考文献：1 张禾瑞 赫炳新 高等代数M 第四版 北京 ：高等教育出版社 1998.1664103 毛纲源 线性代数M 解题方法与技巧归纳 第二版 华中科技大学出版社 1997,7.213241.4 丘维声 抽象代数M 北京 ：高等教育出版社 2003.160190.5 王萼芳 石

25、生明 高等代数M 北京 ：高等教育出版社 1987.176254.6 王萼芳 高等代数教程M 北京清华大学 1996.91184.7 张爱萍 循环矩阵的性质及其对角化J 广西师范自然科学报,2000,12.No.8.168170.8 高吉全 矩阵特征根与特征向量的同步求解方法探讨J 数学通报,1991.12.No.7.2326.9 郭亚梅最小多项式与矩阵的对角化J河南机电高等专科学校学报2006No.410610810张正成 可对角化矩阵的应用J 科技资讯.2007.No.24.252253.11张学元 线性代数能力试题解题M 武汉：华中理工大学出版社, 2000.343712向人晶 矩阵可对角化的简单判定J 数学通报,2003,3.No.12.1315.13靳廷昌 有两个特征根矩阵对角化J 数学通报,1997,11.No.23.5357.14李世余 代数学的发展和展望J 广西大学