矩阵对角化及应用论文

合集下载

矩阵的可对角化及其应用

矩阵的可对角化及其应用

附件:分类号O15商洛学院学士学位论文矩阵的可对角化及其应用作者单位数学与计算科学系指导老师刘晓民作者姓名陈毕专业﹑班级数学与应用数学专业07级1班提交时间二0一一年五月矩阵的可对角化及其应用陈毕(数学与计算科学系2007级1班)指导老师刘晓民摘要:矩阵可对角化问题是矩阵理论中的一个重要问题,可对角化矩阵作为一类特殊的矩阵,在理论上和应用上有着十分重要的意义。

本文对可对角化矩阵做出了全面的概括和分析,并利用高等代数和线性代数的有关理论给出了矩阵可对角化的若干条件,同时也讨论了化矩阵为对角形的求解方法,最后总结出可对角化矩阵在求方阵的高次幂﹑利用特征值求行列式的值﹑由特征值和特征向量反求矩阵﹑判断矩阵是否相似﹑向量空间﹑线性变换等方面的应用.关键词:对角化;特征值;特征向量;相似;线性变换Matrix diagonolization and its applicationChen Bi(Class 1,Grade 2007,The Depart of Math and Calculation Science)Advisor:Lecturer Liu Xiao MinAbstract: Matrix diagonolization problem is an important problem in matrix theory diagonolization matrix, as a kind of special matrix, in theory and application has the extremely vital significance. This paper has made diagonolization matrix analysis and generalization, and using higher algebra and linear algebra are given the relevant theory of matrix several conditions diagonolization, also discussed the matrix of the diagonal shape of solving method, and finally summarized; diagonolization matrix in high power, the policy of using eigenvalue beg determinant by characteristic value and value, feature vector reverse matrix, judgment matrix is similar, vector Spaces, the application of linear transformation, etc.Key words: The diagonalization; Eigenvalue; Feature vector; Similar; Linear transformation引言所谓矩阵可对角化指的是矩阵与对角阵相似,而说线性变换是可对角化的指的是这个线性变换在某一组基下是对角阵(或者说线性变换在一组基下的矩阵是可对角化的),同样可以把问题归到矩阵是否可对角化。

浅谈矩阵的对角化问题(浓缩稿)

浅谈矩阵的对角化问题(浓缩稿)

浅谈矩阵的对角化问题(浓缩版)学号:0807402069 学生姓名:马莉莹 指导老师:朱广俊数学科学学院,2008级,数学与应用数学(师范)摘要:矩阵的对角化是矩阵理论中的一个重要问题,本文利用高等代数的有关理论给出了矩阵可对角化的若干条件;从初等变换、线性方程组、特征子空间等不同角度探究了将一般矩阵和实对称矩阵对角化的若干方法;最后,分析了一些特殊矩阵的对角化问题,如幂等矩阵、幂零矩阵、实对称矩阵和Hermite 矩阵等. 关键词:对角化,特征值,特征向量,相似变换,线性变换.Abstract: Diagonalization of Matrix is an important problem in the matrix theory. We give several conditions of matrix diagonalization by the use of higher algebra related theory. We give some methods of diagonalization of general matrix and real symmetric matrix from different aspects, such as elementary transformation, system of linear equations and characteristic subspace. In the end, we analysis the diagonalization of some special matrix, such as idempotent matrix, nilpotent matrix ,real symmetric matrix and hermite matrix. Keywords : diagonalization ,eigenvalue ,eigenvectors ,similarity transformation ,linear transformation.一.矩阵相似对角化的条件由于矩阵的类型和所在数域的不同,其对角化的条件也不同. 1.任意数域上矩阵相似对角化的条件 充要条件设1,,m λλ 为n 阶方阵A 的m 个互异的特征值,且它们的重数分别为1,,m s s ,1,2,,i m = .A 可对角化⇔A 有n 个线性无关的特征向量⇔对于A 的每个特征值i λ,其代数重数等于其几何重数 ⇔()i i r n s λ-=-I A ⇔A的最小多项式无重根⇔1()mii λ=-=∏I A 0⇔对于A 的每个特征值i λ,都有2()()r r λλ-=-I A I A⇔A 的初等因子都是1次的 ⇔A与某个循环矩阵相似充分条件A 有n 个不同特征值⇒A可对角化A的零化多项式无重根⇒A可对角化2.复数域上Hermite 矩阵必可酉相似于对角矩阵.3.实数域上对称矩阵必可正交相似于对角矩阵.二.矩阵对角化的若干方法(一)一般矩阵对角化的方法特征向量法是将矩阵对角化的常规方法,用该方法解决问题时需要求解齐次线性方程组,过程繁琐.下面介绍其它四种将矩阵对角化的方法. 1.矩阵乘积运算法设12,,,s λλλ 是A在数域F 上全部互异的特征值.其重数分别为12,,,s n n n ,且1sii nn ==∑,记i V λ为A 的属于i λ()1,2,,i s = 的特征子空间. 对()i λ-=I A X 0,有:(1)若A 可对角化,则对A 的每一特征值i λ,都有i n 个与之对应的线性无关的特征向量. (2)A 可对角化的充要条件是对于A 的每个特征值i λ,()ii dim V n λ=.采用类比推测,可得定理1.定理1:设12,,,s λλλ 是A 在数域F 上全部互异的特征值,其重数分别为12,,,s n n n ,且1sii nn ==∑,记i W =()1sj j j iλ=≠-∏I A ()1,2,,i s = . 对()()()12s λλλ---= I A I A I A 0,有:(1)若A 可对角化,则矩阵i W 的列向量组中有对应于i λ的i n 个线性无关的特征向量. (2)A 可对角化的充要条件是()i i rank n =W ()1,2,,i s = .定理1表明,要构造可对角化矩阵A 的相似变换矩阵P ,只需对每一特征值i λ,从矩阵乘积()1sj j j i λ=≠-∏I A 中找出i n 个与之对应的线性无关的特征向量,以这样所得的in n=∑个特征向量为列作一个n 阶矩阵即可.例1:设12202120221001⎛⎫ ⎪ ⎪= ⎪ ⎪⎝⎭A ,求可逆矩阵P ,使得1-P A P 为对角矩阵. 解:由2(1)(5)(1)0λλλλ=-+--=I A ,得 11λ=-(二重),25λ=,31λ= ()()()()()123()50λλλ---=----=因为 I A I A I A I A I AI A ,所以A 可对角化.当11λ=-(二重)时:()()()()123584404840448000λλ--⎛⎫ ⎪-=--=-⎪= ⎪-- ⎪⎝⎭--W I A I A I A I A 取1W 中两个线性无关的特征向量()()12844,04,8,4,0TT=--=--,,,αα. 当25λ=时:()()()()21388808880888000λλ=--⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭=---W I A I A I A I A 取2W 中的特征向量()38,8,8,0T=α当31λ=时:()()()()312000000000050008λλ=--=--⎛⎫ ⎪⎪= ⎪ ⎪--⎝⎭W I A I A I A I A 取3W 中的特征向量()40,0,0,8T=-α.令()1234=,,,P αααα,则1(1,1,5,1)diag -=--P A P2.Jordan 标准形法由于复数域C 上任意n 阶矩阵A 都相似于一个Jordan 矩阵J ,所以存在可逆矩阵P ,使得1-=P A P J .如果J 为对角矩阵,则A 可对角化,否则,A 不可对角化.由于矩阵P 可逆,所以存在一系列的初等矩阵12,,,t P P P ,使得12t = P P P P .于是有:1112112t t ---= P P P A P P P J .可对A 先施行一次初等行变换后,接着施行一次相应的初等列变换,我们称此种初等变换为对A 施行了一次相似变换.显然,可对A 施行一系列的相似变换,将A 化为Jordan 形矩阵J .例2:设460350361⎛⎫⎪-- ⎪ ⎪--⎝⎭=A ,求可逆矩阵P ,使得1-P A P 为对角矩阵. 解:将A 化为Jordan 标准形3121121346026026011350010010(1)(1)361361001r r r r c c c c --⎛⎫⎛⎫⎛⎫+⨯+⨯⎪⎪⎪--−−−−−−→−−−−−−→ ⎪⎪⎪+⨯-+⨯- ⎪ ⎪ ⎪---⎝⎭⎝=⎝⎭⎭A1221200(2)0102001r r c c -⎛⎫+⨯- ⎪−−−−−−→ ⎪+⨯ ⎪⎝⎭由A 的Jordan 标准形知,矩阵A 可对角化且它的特征值为-2,1,1.上述过程对A 共施行了三次相似变换,且三次初等列变换对应的矩阵分别为:123100100120110,010,010001101001⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪=-== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪-⎝⎭⎝⎭⎝⎭P P P所以123120110121⎛⎫⎪==-- ⎪ ⎪--⎝⎭P P P P ,且1211--⎛⎫⎪= ⎪ ⎪⎝⎭P A P .3.λ矩阵标准形法引理1:设A 是n 阶方阵,则必能用初等变换将λ-I A 变为对角矩阵:12()()()()n t t t λλλλ⎛⎫⎪⎪= ⎪ ⎪⎝⎭T 并且多项式 ()(1,2,,)i t i n λ= 的所有根恰好是A 的所有特征值.定理2:设A 是n 阶方阵,{}12()(),(),()n diag t t t λλλλ= T 是对角形λ矩阵,()λP ,()λQ 是可逆的λ矩阵,且满足()()()()λλλλ-=P I A Q T .如果()()()((),)((),())()TTTTTTλλλλλλλ--−−−−−−−−−−→Q I A P I A I T Q Q I.即对()T λ-I A 作初等行变换和初等列变换,使其变为对角矩阵()λT .I 随着()T λ-I A 行的变化而变为()T λQ .则(1) 若12(),(),()n t t t λλλ 的所有根12,,s λλλ 都在F 内,则12,,s λλλ 就是A 的所有特征值.(2) 对于A 的特征值12,,s λλλ ,设第12,,,m ik k k 行是()i λT 的全部为零的行,则()T i λQ 的第12,,,m ik k k 行即构成iV λ的基.其中iV λ为特征值i λ的特征子空间.(3)A 可对角化⇔,(1,2,)i i i r m i s λ∀== ,此处i r 是i λ的重数.根据定理2即可得到λ矩阵标准形法: (1) 作初等变换:()()()((),)((),())()TTTTTTλλλλλλλ--−−−−−−−−−−→Q I A P I A I T Q Q I设{}12()(),(),,()n diag t t t λλλλ= T ,求出12(),(),,()0n t t t λλλ= 的所有解. (2) 若12(),(),,()0n t t t λλλ= 的解都在F 内,并且对每个解i λ都有()i λT 中零行的数目 等于i λ的重数,则A 可对角化,转(3);否则A 不可对角化,结束.(3) 对于A 的任一特征值i λ,若()i λT 的第12,,,m i k k k 行都为零,则取出()T i λQ 的第 1k ,2k , ,m ik 行构作:1111((),,(),,(),,())m s m sT TTTk kk s k s λλλλ= T Q Q Q Q则12112(,,,)sm m s m diag λλλ-= T AT I I I .例3:设132132264⎛⎫⎪=--- ⎪ ⎪⎝⎭A ,求可逆矩阵T ,使得1-T A T 为对角矩阵. 解:作初等变换:()2112100100100,33601002011222410021T λλλλλλλλ--⎛⎫⎛⎫⎪ ⎪-=-+-→-+-+- ⎪ ⎪ ⎪ ⎪---⎝⎭⎝⎭I A I 按上述方法:(1)记2100002()00λλλλ⎛⎫⎪= ⎪ +⎪⎝⎭-T ,100()112201T λλ⎛⎫⎪=-+- ⎪ ⎪-⎝⎭Q 则1230,2λλλ===(2)当120λλ==时,(0)T 中零行的数目0=的重数2=-当32λ=时,(2)T 中零行的数目2=的重数1=-.所以A 可对角化.(3)当120λλ==时,()()()1001000,00001120021T ⎛⎫ ⎪=- ⎪ ⎪-⎝⎭T Q 取(0)T Q 中与(0)T 中零行所对应的特征向量()11,1,2T=-α,()22,0,1T=-α 当32λ=时,()()()1001002,200011200221T ⎛⎫ ⎪=-- ⎪ ⎪-⎝⎭T Q 取(2)T Q 中与(2)T 中零行所对应的特征向量()31,1,2T=--α.令()123121,,101212--⎛⎫ ⎪== ⎪ ⎪--⎝⎭T ααα,则1002-⎛⎫⎪⎪ ⎪⎝⎭T A T =4. 数字矩阵对角形法若矩阵A 在数域F 上可对角化,则存在F 上的可逆矩阵T ,使得1-=T AT B 为对角矩阵,且B 的主对角线上的元素为A 的全体特征值.由于矩阵T 可逆,所以存在一系列的初等矩阵12,,,s T T T ,使得12s = T T T T .于是:11111112s s s ----- B =TA T =T T T A T T T ,做初等变换:⎛⎫⎛⎫→⎪ ⎪⎝⎭⎝⎭A B I T . 即对A 施行一系列的初等行变换和初等列变换,使其变为对角矩阵B ,对I 只施行相应的初等列变换变为T .在施行初等变换时,可施行若干次行(或列)变换后再施行若干次相应的列(或行)变换,只要保持变换后所得矩阵与A 相似即可.例4:若1111111111111111⎛⎫ ⎪-- ⎪= ⎪-- ⎪--⎝⎭A ,求可逆矩阵T ,使得1-T A T 为对角矩阵.解:作初等变换:200002001111002011110002111111111111444100031110100444001013114440011131444-⎛⎫ ⎪ ⎪⎛⎫⎪ ⎪-- ⎪ ⎪⎪ ⎪-- ⎪⎪⎪--⎛⎫ ⎪=→ ⎪ ⎪ ⎪⎝⎭⎪ ⎪--- ⎪ ⎪⎪ ⎪ ⎪ ⎪--- ⎪ ⎪⎝⎭⎪ ⎪--- ⎪⎝⎭A I 所以A 可对角化.令1111444311144413114441131444⎛⎫ ⎪ ⎪ ⎪---⎪=⎪ ⎪--- ⎪ ⎪ ---⎪⎝⎭T ,则有120000200002002--⎛⎫ ⎪⎪= ⎪ ⎪⎝⎭T A T .利用初等变换将矩阵对角化时,我们可以从变换后的最终矩阵中直接读出相似变换矩阵和对角矩阵,大大简化了求解过程.(二)实对称矩阵对角化的方法Schmidt 正交法是将实对称矩阵对角化的基本方法,使用该方法时需要牢记公式且计算量较大.下面我们介绍另外两种方法. 1.直接正交法该方法从向量正交的基本定义出发,直接从特征子空间中求出正交向量,易于理解和掌握,且在特征值出现重根的情况下,计算量也大为减少.例5:设 1333313333133331---⎛⎫ ⎪--- ⎪= ⎪--- ⎪---⎝⎭A ,求正交矩阵P , 使得1-P A P 为对角矩阵. 解:由3(4)(8)0λλλ-=+-=I A ,得14λ=-(三重),28λ=. 设41234(,,,)T x x x x R =∈X当14λ=-时,解齐次线性方程组(4)--=I A X 0,得1243x x x x =+-.先取一个特征向量1(1,1,0,0)T =α. 设特征向量22222(,,,)T a b c d =α.因2α与1α正交,从而有220a b +=.又因为2222a b d c =+-,所以可得2222a d c =-. 取211(,,0,1)22T =-α.再设特征向量33333(,,,)T a b c d =α.因3α与1α和2α都正交,从而有330a b +=,33311022a b d -+=.又因为3333a b d c =+-,所以可得333a c =-.取3(2,2,6,2)T =---α. 现将1α,2α,3α都单位化:122,,0,022T⎛⎫= ⎪⎪⎝⎭β,2666,,0,663T ⎛⎫=- ⎪ ⎪⎝⎭β,33333,,,6626T⎛⎫=--- ⎪ ⎪⎝⎭β. 当28λ=时,可求得单位特征向量:41111,,,2222T⎛⎫=-- ⎪⎝⎭β.令1234(,,,)=P ββββ,则()14,4,4,8T diag ----P AP =P AP =.2.度量矩阵法对于n 维欧氏空间V ,令1,,n αα是它的一个基,它的度量矩阵()()()()1111,,,,n n n n ⎛⎫⎪= ⎪⎪⎝⎭A αααααααα是正定矩阵,于是A 合同于单位矩阵I ,即可求得n 阶可逆矩阵U ,使得T =U AU I .利用U 和V 的基1,,n αα作一个新基:121(,,,)(,,)n n = βββααU .那么,新基的度量矩阵即为:()()()()1111,,,,n Tn n n ⎛⎫⎪= ⎪ ⎪⎝⎭=U A U Iββββββββ.所以12,,,n βββ是欧式空间V 的标准正交基.例6:设0111101111011110-⎛⎫ ⎪- ⎪= ⎪- ⎪-⎝⎭A ,求正交矩阵P , 使得1-P A P 为对角矩阵. 解:由3(1)(3)0λλλ-=-+=I A ,得11λ=(三重),23λ=-. 当11λ=时,解齐次线性方程组()-=I A X 0,得基础解系 1(1,1,0,0)T =α,2(1,0,1,0)T =α,3(1,0,0,1)T =-α当23λ=-时,解齐次线性方程组(3)--=I A X 0,得基础解系4(1,1,1,1)T =--α 则 1234,,,αααα是4R 一组基.记其度量矩阵为B ,那么21101210112004-⎛⎫ ⎪-⎪= ⎪-- ⎪⎝⎭B 对矩阵⎛⎫ ⎪⎝⎭B I 作合同变换:⎛⎫ ⎪⎝⎭B I =2110121011200004100001000010001-⎛⎫ ⎪- ⎪ ⎪--⎪ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭→1000010000100001263026663003630002102⎛⎫ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪⎝⎭.取263026663003630002102⎛⎫-⎪ ⎪ ⎪⎪⎪= ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭U ,则有1111T ⎛⎫⎪⎪= ⎪ ⎪⎝⎭U B U . 利用U 和基1234,,,αααα作新基:12341234(,,,)(,,,)=ββββααααU .则: 122,,0,022T⎛⎫= ⎪⎪⎝⎭β, 2666,,,0663T⎛⎫=- ⎪ ⎪⎝⎭β. 33333,,,6662T⎛⎫=- ⎪ ⎪⎝⎭β, 41111,,,2222T⎛⎫=-- ⎪⎝⎭β.由于1234,,,ββββ的度量矩阵T =U B U I ,故1234,,,ββββ是4R 的标准正交基.令1234(,,,)=P ββββ,则P 是正交矩阵且1T -P AP =P AP .三.特殊矩阵的对角化 1.幂等矩阵定理3:n 阶幂等矩阵A一定可以对角化,并且A的相似标准形是 0r⎛⎫⎪⎝⎭I ,其中()r rank =A ,r I 是r阶单位矩阵.证明: 因为2=A A ,所以A 有零化多项式2()(1)g λλλλλ=-=-,因为()g λ无重根,所以A可对角化.而A 的特征值只有0和1,所以A 的相似标准形是0r⎛⎫⎪⎝⎭I ,其中()r rank =A .由该定理可以推出幂等矩阵的若干性质: 性质1:幂等矩阵A 的迹等于A 的秩.证明:设A 是数域F 上的一个n 阶幂等矩阵,()r rank =A .如果0r =,则()0()rank tr ==A A .如果r n =,则=A I .从而()()rank n tr ==A A .下面设0r n <<.由A 的相似标准形0r⎛⎫⎪⎝⎭I 得: ()((,0))()r rank r tr diag tr ===A I A .性质2:任意n 阶矩阵A 都可以表示成为一个可逆矩阵与一个幂等矩阵的乘积. 证明:设n 阶方阵A 的秩为r ,则存在n 阶可逆矩阵,P Q 使得: 000r ⎛⎫=⎪⎝⎭I PA Q 所以1111100()()0000r r -----⎛⎫⎛⎫== ⎪ ⎪⎝⎭⎝⎭I I A PQ P Q Q Q . 令11--=B P Q ,1000r -⎛⎫=⎪⎝⎭I C Q Q .易知B 为可逆矩阵.因为2=C C ,所以C 为幂等矩阵.即任意n 阶矩阵A 都可以表示成为一个可逆矩阵与一个幂等矩阵的乘积.2.幂零矩阵引理2:若()f λ 为A 的特征多项式,()m λ为A 的最小多项式,则()()f m ==A A 0. 引理3:设12,,,n λλλ 为n 阶矩阵A 的特征值,则对任意的多项式()f x 有()f A 的特征值为12(),(),,()n f f f λλλ .幂零矩阵具有下列性质:性质3:A 为幂零矩阵的充分必要条件是A 的特征值全为0.证明:(必要性) 若A 为幂零矩阵,则存在正整数k ,使得k =A 0.令0λ为A 的任意一个特征值,则存在≠α0,使得0λ=A αα.由引理3知0k λ为k A 的特征值. 所以存在 ≠β0,使得 0k k λ=A ββ,从而有00k λ=即有00λ=.又由k =A 0,知00kk ==⇒=A A A ,所以 0(1)(1)00k k ⨯-=-=-=-⋅=I A A A . 所以00λ=为A 的特征值.由0λ的任意性知A 的特征值全为0.(充分性)因为A 的特征值全为0, 所以A 的特征多项式为()n f λλλ=-=I A ,由引理2知()n f ==A A 0,所以A 为幂零矩阵.性质4:若A 为幂零矩阵且≠A 0,则A 不可对角化.证明:若A 可对角化,则存在可逆矩阵P ,使得1-=A P DP ,此处D 是n 阶对角形.若A 为 幂零矩阵,则存在正整数k ,使得k =A 0,即: 11()k k k --===A P DP P D P 0,因为1110kk k k k ---=====P D P P D P P P D D D ,所以有: 10,,-====D D 0A P DP 0, 与题设矛盾.3.幂幺矩阵性质5:幂幺矩阵在复数域上可对角化.证明:若A 为幂幺矩阵,则存在正整数k ,使得k =A I ,所以A 有零化多项式()1k g λλ=-. 因为在复数域上,()g λ的根都是k 次单位根,故()g λ无重根,所以A 可对角化.注意:A 在实数域上不一定可对角化! 例如0110-⎛⎫=⎪⎝⎭A ,满足4=A I ,即A 为幂幺矩阵,但是2()1f λλλ=-=+I A 在实数域上无根,所以A 在实数域上不可对角化.4.实对称矩阵性质6:实对称矩阵的不同特征值的特征向量相互正交.性质7:设λ是实对称矩阵的k 重特征值,则对应于特征值λ,矩阵有k 个线性无关的特征向量. 定理4:设A是一个n n ⨯实对称矩阵.则存在一个正交矩阵P,使得()112,,,Tn diag λλλ-== P AP PAP ,并且i λ是实数,1,2,,i n = .证明:设A的互不相等的特征值为12,,,()s s n λλλ≤ ,并且它们的重数依次为1212,,,()s s r r r r r r n +++= .则对于特征值(1,2,,)i i s λ= ,恰有i r 个线性无关的实特征向量.把它们正交化并单位化,即得i r 个单位正交的特征向量.由12s r r r n +++= 知,这样的特征向量共可得n 个.由于不同特征值的特征向量正交,故这n 个单位特征向量两两正交,以它们为列向量作成正交矩阵P ,则:1T -=P AP P AP 为一个实对称矩阵111,,,,,,s s sdiag r r λλλλ⎛⎫⎪ ⎪⎝⎭.5.Hermite 矩阵欧氏空间实质上是实数域上的一个内积空间.类似地考虑复数域上的内积空间—酉空间和酉空间上的线性变换.与正交变换和实对称矩阵类似,酉空间中有酉变换与Hermite 矩阵.性质8:设n n C ⨯∈A 是Hermite 矩阵,则A 的特征值均为实数.证明:设λ为A 的特征值,α为其对应的特征向量,即λ=A αα,那么: (,)(,)(,)(,)(,)(,)λλλλ=====ααααααααααααA A 但(,)0>αα,所以λλ=,即λ为实数.性质9:设n n C ⨯∈A 是Hermite 矩阵,则对应于A 的不同特征值的特征向量必正交. 证明:设,λμ是A的两个不同的特征值,,αβ分别是它们所对应的特征向量,则有λ=A αα,μ=A ββ.(,)(,)(,)(,)(,)(,)λλμμ=====αβαβαβαβαβαβA A ,即()(,)0λμ-=αβ.由于A 的特征值为实数,也即()(,)0λμ-=αβ.又因为λμ≠,所以(,)0=αβ,即,αβ正交.引理4:设n n C ⨯∈A ,则存在一个酉矩阵P ,使得1-P A P 是一个上三角形矩阵.定理5:设n n C ⨯∈A ,并且A是Hermite 矩阵,则存在一个酉矩阵P , 使得()112,,,Hn diag λλλ-== P AP PAP ,并且i λ是实数,1,2,,i n = .证明:由引理4知存在一个酉矩阵P ,使得 ()1H ij n n g -⨯===G P AP P AP 是一个上三角形矩阵.又P 是一个酉矩阵,故G 也是Hermite 矩阵.于是,对任意,,1i j i j n ≤<≤,都有ij ji g g =,这迫使当1,2,,,1,2,,,i n j n i j ==≠ 时,有0ij g =;并且i ii g λ=是实数,1,2,,i n = .因此,Hermite 矩阵必定可以对角化,且它的特征多项式的复数根都是实数.。

浅谈矩阵对角化及其应用(米亚兄) - 天津商业大学商学院【优秀资料】

浅谈矩阵对角化及其应用(米亚兄) - 天津商业大学商学院【优秀资料】

浅谈矩阵对角化及其应用(米亚兄)- 天津商业大学商学院【优秀资料】(可以直接使用,可编辑完整版实用资料,欢迎下载)浅谈矩阵对角化及其应用写在前面:结识高等代数已经快一年了,我们从最初的认识行列式,一直到到现在的欧几里得空间,逐一学习了线性方程组、矩阵、多项式、二次型、线性空间、线性变换,现在就浅谈一下自己对矩阵对角化及其应用的认识。

众所周知:n维向量空间V中的线性变换δ可否对角化的问题是高等代数中十分重要的内容,而δ可对角化的充要条件是δ关于V的矩阵A可对角化。

内容摘要:文章综述了矩阵可以对角化的条件,讨论了可对角化矩阵的基本性质和结论,给出了矩阵(特殊矩阵如是对称阵)对角化的基本方法,以及对应特征多项式的性质,最后讨论其在特征值、特征向量方面的应用。

关键词:矩阵对角化特征多项式特征值特征向量导言:文章由矩阵可对角化出发,说明矩阵可对角化的条件、讨论了可对角化矩阵的基本性质和结论,给出了矩阵(特殊矩阵如是对称阵)对角化的基本方法,以及对应特征多项式的性质,最后讨论其在特征值、特征向量方面的应用。

具体内容:1、矩阵可对角化的条件:1)设δ是n维线性空间的一个线性变换,δ的矩阵可以在某一组基下维对角矩阵的充分必要条件是δ有n 个线性无关的特征向量。

2)方块矩阵A被称为可对角化的,如果它相似于对角矩阵,就是说,如果存在一个可逆矩阵P使得P−1AP是对角矩阵。

3)设A 是数域F上的n阶矩阵,如果存在F上n阶可逆矩阵T,使得T1-AT=∧,那么,就说矩阵A 是可以对角化的。

可对角化矩阵的基本性质和结论:1)数域F上n阶矩阵A可以对角化的充要条件是A有n个线性无关的特征向量。

2)数域F上n阶矩阵A在F内有n个不同的特征根,那么A可以对角化。

3) 属于不同特阵值的特征特真向量是线性无关的。

4)如果在n 维空间V 中,线性变换δ的特征多项式在数域P 中有n 个不同的根,即δ有n 个不同的特征值,那么在某组基下的矩阵是对角形的。

对角化矩阵的应用本科

对角化矩阵的应用本科

XXX学校毕业论文(设计)对角化矩阵的应用学生姓名学院专业班级学号指导教师2015年 4 月 25 日毕业论文(设计)承诺书本人郑重承诺:1、本论文(设计)是在指导教师的指导下,查阅相关文献,进行分析研究,独立撰写而成的.2、本论文(设计)中,所有实验、数据和有关材料均是真实的.3、本论文(设计)中除引文和致谢的内容外,不包含其他人或机构已经撰写发表过的研究成果.4、本论文(设计)如有剽窃他人研究成果的情况,一切后果自负.学生(签名):2015 年4月25日对角化矩阵的应用摘要矩阵对角化问题是矩阵理论中一个关键性问题.本文借助矩阵可对角化条件,可对角化矩阵性质和矩阵对角化方法来研究可对角化矩阵一些应用,包括求方阵的高次幂,反求矩阵,判断矩阵是否相似,求特殊矩阵的特征值,在向量空间中证明矩阵相似于对角矩阵,运用线性变换把矩阵变为对角矩阵,求数列通项公式与极限,求行列式的值.【关键词】对角化;特征值;特征向量;矩阵相似;线性变换Application of diagonalization matrixAbstractMatrix diagonalization problem is the key issue in the matrix theory. In this paper, by using matrix diagonalization conditions, diagonalization matrix properties and matrix diagonalization method we study some applications of diagonalization matrix, including for high-order exponent of matrix, finding the inverse matrix, matrix to determine whether it is similar, the eigenvalue of special matrix, in the vector space that matrix similar to a diagonal matrix, using linear transformation matrix is a diagonal matrix, for the series of general term formula and limit, the determinant of value.[Key words] The diagonalization; Eigenvalue; Feature vector; Similar; Linear transformation目录引言 (1)1矩阵对角化 (1)1.1矩阵对角化的几个条件 (1)1.2对角化矩阵的性质 (3)1.3 矩阵对角化的方法 (5)2对角化矩阵的应用 (5)2.1求方阵的高次幂 (5)2.2反求矩阵 (6)2.3判断矩阵是否相似 (7)2.4求特殊矩阵的特征值 (7)2.5在向量空间中应用 (7)2.6在线性变换中应用 (7)2.7求数列通项公式与极限 (8)2.8求行列式的值 (11)2.9对角化矩阵在其他方面的应用 (12)参考文献 (14)致谢 (15)引 言现如今,我们所提到的矩阵对角化其实质指的就是矩阵和对角阵存在相似的地方,其中我们学过的线性变换也是可对角化的,其原理是指在某一组基的作用下这个线性变换可以变为对角阵(或者可以说是在某一组基的作用下这个线性变换的矩阵是可对角化的),当然刚刚提到的这个问题其实我们可以把它归类到矩阵是否可对角化的问题中去,因为其两者本身就是相辅相成的.当然本篇文章我们主要是研究和探索判定矩阵可对角化的诸多条件,以及我们如何去运用矩阵对角化的有关性质,来把将矩阵化为对角形的问题进行解决.与此同时,我们也在研究和探索中发现了它在其他方面一些重要的运用.1矩阵对角化我们所涉及的矩阵都是可以对角化的,其原理是指通过矩阵的一系列初等变换(指:行、列变换)后,就能够得到一个特殊的矩阵,其特殊性在于只有在其主对角线的数上不全为零,然而其他位置的数则是全部为零(那么这个特殊的矩阵就可以被我们称为对角阵),这一整个的变换过程就被我们称为矩阵的对角化.当然值得我们注意的是,我们所学过的矩阵并非都能对角化的,这个是有条件限制的.1.1矩阵对角化的几个条件引理]1[1 设n n P B A ⨯∈,,且,2A A =,2B B =BA AB =,则存在可逆矩阵P ,使B A ,可同时对角化.引理]2[2 如果n n n P diag P ⨯∈=),,,(21λλλ 的n 个对角元互不相同,矩阵n n P B ⨯∈,那么BP PB =当且仅当B 本身就是对角阵.因为任何一个幂等矩阵)(2A A A =一定相似于一个对角矩阵⎥⎦⎤⎢⎣⎡000rE ,所以任何一个对角矩阵都是能够进行谱分解的,即∑==n i i i A A 1λ,其中i λ是矩阵A 的特征值,矩阵i A 为幂等矩阵,那么是否任意有限个幂等矩阵的线性组合都可以对角化呢?有如下结论:定理]3[1 若,2211n n k k k A ∆++∆+∆=n k k k ,,,21 是n 个数,n ∆∆∆,,, 21是n 个幂矩阵,并且他们两两可替换,)(,j i i j j i ≠∆∆=∆∆,则矩阵A 可对角化.证明 若n ∆∆∆,,, 21是n 个幂矩阵,并且两两可换,则一定有一个可逆矩阵1P ,使得n ∆∆∆,,, 21,可同时对角化.n n n n P D P P D P 111111--=∆=∆,, )(1是对角矩阵,,n D D , P D k D k D k P P D k P P D k P P D k P k k k A n n n n n n )()()()(2211112211112211+++++++=∆++∆+∆=---- ,由是对角矩阵,,n D D 1知n n D k D k D k +++ 2211同样是对角矩阵,即矩阵A 为对角化的矩阵.定理]4[2 如果n n P A ⨯∈,21λλ,是它两个不相同的特征值,那么矩阵A 可对角化⇔一定有幂等矩阵∆,满足∆-+=)(121λλλE A .证明 必要性:如果A 是一个对角化的矩阵,那么就一定会有一个可逆的矩阵P ,满足∆=⎥⎦⎤⎢⎣⎡=-2211111E E AP P λλ 是一个对角阵.()()()121211121211111211000-----⎥⎦⎤⎢⎣⎡-+=⎥⎦⎤⎢⎣⎡-+=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡-+==P E P E P E P P E P P E P PAP A λλλλλλλλλ, 并且∆相似于2121212000∆=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡---P E P P E P P E P ,若∆为幂矩阵,则一定有一个幂矩阵∆满足∆-+=)(121λλλE A .充分性:若存在∆使得∆-+=)(121λλλE A ,因为∆是幂矩阵,所以一定会有一个T ,满足T E T ⎥⎦⎤⎢⎣⎡=∆-210, ()()T E E T T E E T T E T E A ⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-+=⎥⎦⎤⎢⎣⎡-+=---2211112121121)0(0λλλλλλλλ, 因此,T E E T AT T ⎥⎦⎤⎢⎣⎡=--221111λλ, 即矩阵A 为可对角化的.定理]5[3 设矩阵n n P A ⨯∈存在n 个不同的特征值,则对于矩阵n n P B ⨯∈,BA AB =,当且仅当矩阵B A ,同时可以对角化.证明 必要性 若矩阵A 存在n 个特征值,且这些特征值是互不相同的数,则矩阵A 为对角化的矩阵.设AP P T 1-=,其中),,,(21n diag T λλλ =,则ABP P BP APP P BP P T 1111)(----==T BP P AP BPP P )(111---==,即T 与BP P 1-是可以进行交换的,因此得知BP P 1-是对角矩阵,且矩阵B 也是为对角化的矩阵.充分性 如果矩阵B A ,可以同时进行对角化,那么一定存在一个可逆阵P ,使得P D P A 11-=,P D P B 21-=(其中为21D D ,对阵),BA P D PP D P P D D P P D D P P D PP D P AB =====------11211212112111,因此我们可以通过上述的一系列条件,来求出A 的特征值,且这是两个相互不同的数.从而我们得出了矩阵对角化的成立的条件:如果∆=∆2这个条件成立,那么就认为矩阵A 可对角化,否则就认为矩阵A 不能可对角化,其中)(/)(21λλλ--=∆E A .1.2对角化矩阵的性质定理]6[4 设A 为数域P 上的一个n 阶的矩阵,且它为可对角化的,t λλλ,,,21 是A 的相互不同的特征根,则一定会有n 阶的t A A A ,,,21 满足(1)t t A A A A λλλ+++= 2211;(2)E E A A A t ,21=+++ 是单位矩阵;(3)i i A A =2;(4)j i A A j i ≠=,0,其中1-=T TB A i i . 证明 (1)如果A 可对角化,那么在数域P 上一定会存在一个可逆矩阵T ,并且它的阶数为n 阶,满足B AT T t =⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=-λλλ00211 , 其中i λ的重数为i s ,由于矩阵⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡++⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=110000111 t B λλ, 将它记为t t B B B λλλ+++ 2211,因此,)()()(1111122111----++=+++==T TB T TB T B B B T TBT A t t t t λλλλλ ,将其记为t t A A A λλλ+++ 2211,其中1-=T TB A i ,所以t t A A A A λλλ+++= 2211.(2)如果每个i B 为对角形的幂矩阵,那么E B B B t =+++ 21,E TET T TB T TB T TB A A A t t ==+++=+++----11121121 ,故E A A A t =+++ 21.(3)如果1-=T TB A i i ,那么i i i i i i i i i i A T TB T TB T B TB T TB T TB T TB T TB A ======-------112111112))((,故i i A A =2.(4)当j i ≠时, 0))((11111====-----T B TB T TB T TB T TB T TB A A j i j i j i j i ,0为零矩阵,故j i A A j i ≠=,0.例1 在数域P 上,若已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=6788152051115A 的三个特征根分别是3,2,1,则一定会有一个⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=211243132T ,满足B AT T =⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=-3000200011,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=-1111342561T ,将矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=10030102001B , 记32132B B B ++,则,3211321132)32(A A A T B B B T TBT A ++=++==-- 其中1-=T TB A i i ,于是⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------=222222111,134412163912,2566151841012321A A A , 并且满足:(1)32132A A A A ++=;(2)E A A A =++321;(3))3,2,1(2==i A A i i ; (4)j i A A j i ≠=,0.可以通过一个比较具体的可对角化矩阵,很直观地反映上述所说的性质是成立的.1.3 矩阵对角化的方法1.3.1 运用矩阵初等变换的方法在数域P 上,一个n 维空间V ,研究和探讨它能否可以找到一组基,并且在此基的作用下,所有的矩阵都是对角化的矩阵;发现这种基存在时, 如何去探索它是一个线性代数学上相当重要的问题,可以利用矩阵的初等变换的方法来解决此问题.当发现矩阵A 不能够实现对角化的时候,同样可以经过相近的一系列变换后,化简出矩阵A ,并且能够判定它是否可以对角化.类似地,可有矩阵E Q Q Q T s s 111111-----= ,做如下的初等变换,则可以将矩阵A 化简为对角形矩阵B ,并且可以求得T 或由B 求A 的一系列特征值.1.3.2 求解齐次方程组的方法设矩阵A 是实对称矩阵,则求证交矩阵T 使得),,,(211n diag AT T λλλ =-的问题,一般的解法为:(1)求其特征值; (2)求其对应的特征向量;(3)写出矩阵T 及),,,(211n diag AT T λλλ =-.从而可以求出正交矩阵T ,可以避免了商的繁琐运算.定理]7[5 设A 是实对称矩阵,则有)1(21重,-n λλ,n αααβ,,,, 321对应于21λλ,,记)(1βL 由1β生成的一个空间,且)(32n L ααα,,, 由n ααα,,, 32生成的空间.2对角化矩阵的应用2.1求方阵的高次幂例2 设在数域P 上,有一个二维的线性空间V ,21ξξ,是这个线性空间V 的一组基,那么线性变换σ在21ξξ,这组基的作用下的矩阵⎥⎦⎤⎢⎣⎡-=0112A ,试通过上述给出的条件计算出矩阵k A .解 通过分析上述的条件,我们应该先计算线性变换σ在线性空间V 的另一组基21ηη,作用下的矩阵,令[][]⎥⎦⎤⎢⎣⎡--=2111,,2121ξξηη, 则⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡---10112111011211122111011221111, 易知⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡1011011k k, 再运用上面得出的几个关系⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡---10112111011221111, 即⎥⎦⎤⎢⎣⎡+--+=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡--⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-11111210121112111101121-1-101-121kk k k k k k.2.2反求矩阵例3 设有一个实对称矩阵A ,且它的阶数为3阶,已知11321==-=λλλ,,1λ对应于T P )1,1,0(1=,求解A .解 根据矩阵A 是3阶实对称矩阵的条件,我们可以推出矩阵A 可以对角化的结论,即得出矩阵A 是由三个线性无关的特征向量组成的结论,并且132==λλ对应于T X X X P ),,(321=,因为它和1P 正交,即003211=++=⋅X X X P P ,所以可以求出T T P P )1,1,0()0,0,1(32-==,,它们分别对应132==λλ.取 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==1000100011-01101010),,(321B P P P P ,, 则B AP P =-1,于是⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-==-010********21000121211000100011011010101PBP A . 2.3判断矩阵是否相似例4 请判断下述三个矩阵是否会相似⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=300020102,300120012,300020002321A A A . 解 我们可以很容易的得出三个矩阵321,,A A A 的特征值分别都是21=λ(二重),32=λ,其中矩阵1A 已经是对角阵,所以我们只需要进一步判断两个矩阵32,A A 是否都可以对角化.通过21=λ,0)2(2=-X A E ,可以推出T )0,0,1(1=α,因为21=λ,是一个二重的特征值,但是却只有一个特征向量与之所对应,那么我们可以推出矩阵2A 与矩阵1A 不相似的结论.通过21=λ,0)2(3=-X A E ,得出T T )0,1,0(,)0,0,1(21==ηη,通过32=λ,0)3(3=-X A E ,得出T )1,0,1(3=η,通过上述所推出的结论,我们可知矩阵3A 有三个线性无关的特征向量,即矩阵3A 与矩阵1A 这两个矩阵相似. 2.4求特殊矩阵的特征值例]8[5 设有一个实对称矩阵A ,并且它的阶数为n 阶,满足A A 22=,n r A r <=)(,求出A 的全部特征值.解 假设λ为矩阵A 的一个特征值,而我们令ξ为矩阵A 的特征向量,它对应于特征值λ,因为λξξ=A ,所以ξλλξξ22==A A ,又因为A A 22=,所以λξξξ222==A A ,即λλ22=,由此我们可以推出02或=λ,根据矩阵A 是实对称矩阵的这个条件,我们可以断定矩阵A 一定能够进行对角化,即⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡=0022~ B A ,与r A r =)(,所以A 的秩数就是2的个数,以及A 有r 个2和)(r n -个0的特征值. 2.5在向量空间中应用例]9[6在n 维的V 空间中,有一个复矩阵,并且它的阶数为n 阶,还有一个复数α, 令{}{}0)(,)(21=-∈=∈-=βαβββαA E V W V A E W ,则矩阵A 相似于对角阵,并且{}021=⋂W W .证明 因为对于任意一个210W W X ⋂∈,则有βα)(0A E X -=和0)(0=-X A E α,所以0)(2=-βαA E .又因为发现矩阵A 相似于对角阵,所以我们可以推出0)(0=-X A E α与0)(2=-βαA E 两个的解空间是完全相同的,即{}021=⋂W W . 2.6在线性变换中应用例]10[7 设()1][>n X P n 是数域P 上的一个全体,且它是一个次数小于n 的多项式与零多项式,则请通过所学的进一步判断在n X P ][的任一组基下,矩阵通过微分变换τ能否变为对角形矩阵.证明 如果取()!1!211--n X X X n ,,,, , 那么矩阵可以表示为⎥⎦⎤⎢⎣⎡0001-n E ,所以有nA E λλ=-. 如果在某一组基的作用下,微分变换τ的矩阵B 为对角矩阵,由已知的矩阵B A ~可推出矩阵A 可对角化,那么就会存在一个可逆矩阵T 能够使得B AT T =-1,所以1-=TBT A .通过已知的微分变换τ的全为零,可以推出0=B ,0=A 这是不可能的,所以在n X P ][的任何一组基的作用下,微分变换τ的矩阵都不可能成为对角阵.2.7求数列通项公式与极限例]11[8 设两个数列{}{}n n q p ,都满足条件1,,21111==+=+=++q p q p q q p p n n n n n n ,则请求解nn n q p∞→lim .解 把已知条件中的几个递推关系组n n n n n n q p q q p p +=+=++11,2,通过化简改写成下面的列矩阵的形式:⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡==⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡++111111211121q p q p q p nn n n n ,由⎥⎦⎤⎢⎣⎡=1121A 和0=-A E λ,可以求出A 的21,2121-=+=λλ,并且21λλ,分别对应T T X X )1,2(,)1,2(21-==.取),(21X X X =,则⎥⎦⎤⎢⎣⎡-=-21212211X ,1210021-⎥⎦⎤⎢⎣⎡-+=X X A , 从而⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--+-++=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡-+=⎥⎦⎤⎢⎣⎡++++-++2)21()21(2)21()21(112100211111111n n n n nn n X X q p , 因此2)21()21(nn n p -++=,2)21()21(n n n q --+=, 并且2)21()21()21(2)21(2lim lim =--+-++=∞→∞→n n nn n nn n q p . 例9 已知),2,1(2,2),(,11111 =+=+=>==+++n ba b b a a b a n n n n n n βαβα这四个条件,请证明n n n n b a ∞→∞→lim lim 及存在并且相等,给出证明过程,同时请求出这两个的极限值. 证明 把已知条件中的递推关系组作进一步简化推出434,2211n n n n n n b a b b a a +=+=++,然后再改写为另一种矩阵的形式:⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡==⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡++11114341212143412121b a b a b a nn n n n ,由⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=43412121A 和0=-A E λ,可以求出A 的14121==λλ,,并且21λλ,分别对应()()TTX X 11,1221,,=-=,取()⎥⎦⎤⎢⎣⎡-==1112,21X X X ,则 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=-323131311X ,110041-⎥⎥⎦⎤⎢⎢⎣⎡=X X A , 因为⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡++1004111X b a n n ,⎥⎦⎤⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⋅+⋅-+⋅-+⋅=⎥⎦⎤⎢⎣⎡-βα324131314131324231314231111n nn n b a X ,所以βα⋅⎪⎭⎫ ⎝⎛+⋅-+⋅⎪⎭⎫ ⎝⎛+⋅-=+3242313142311n n n a ,βα⋅⎪⎭⎫⎝⎛+⋅+⋅⎪⎭⎫ ⎝⎛+⋅-=+3242313142311n n n b ,即n n n n b a ∞→∞→=+=lim 3231lim βα. 例10 设有10=x ,e x =1,)1(11≥⋅=-+n x x x n n n 这三个条件,请求出n n x ∞→lim .解 从已知的三个条件可以推出),2,1(0 =>n x n ,以及)ln (ln 21ln 11-++=n n n x x x ,令n n x a ln =,则00=a ,11=a ,)1()(2111≥+=-+n a a a n n n ,所以 ⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡==⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡-+0111012121012121a a a a a a nn nn n , 由⎥⎥⎦⎤⎢⎢⎣⎡=012121A 和0=-A E λ,求得A 的21121-==λλ,,并且21λλ,分别对应TT X X )121(,)11(21,,-==.取),(21X X X =,令⎥⎥⎦⎤⎢⎢⎣⎡-=-11211321X,121001-⎥⎥⎦⎤⎢⎢⎣⎡-=X X A , 则⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=⎥⎦⎤⎢⎣⎡⎥⎥⎦⎤⎢⎢⎣⎡-=⎥⎦⎤⎢⎣⎡+-+n n nn n X X a a )21(1)21(1320121001111, 从而推出:))21(1(32nn a --=,即))21(1(32n e x n --=,32lim e x n n =∞→.例11 设11=x ,nn x x +=+111,求n n x ∞→lim .解 令1+=n n n a a x ,根据条件nn x x +=+111,将其简化为n n n a a a +=++12,然后再写成矩阵)2(0111011112111≥⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡==⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡--+n a a a a a a n n n n n , 由⎥⎦⎤⎢⎣⎡=0111A 和0=-A E λ,求出A 的βλαλ=-==+=25125121,,且21λλ,分别对应的是T T X X )1(,)1(21,,βα==,取⎥⎦⎤⎢⎣⎡==11),(21βαX X X ,则100-⎥⎦⎤⎢⎣⎡=X X A βα, ⎥⎦⎤⎢⎣⎡--=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡++++-+112211511100n n n n nn n X X a a βαβαβα, 即2151)()(1lim lim limlim 1122111-==--=--==++∞→++++∞→+∞→∞→ααββααββαβαn n n n n n n n n n n n n a a x . 2.8求行列式的值例]12[12 设有一个n 阶的行列式,化简并求出它的值.)0(sin cos 21001cos 2100000001cos 21000001cos 21000001cos 2≠=ααααααn D ,解 按照第一列展开的21cos 2---=n n n D D D ,可以写成矩阵的另外一种形式⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡-=⎥⎦⎤⎢⎣⎡---211011cos 2n n n n D D D D α, 记矩阵⎥⎦⎤⎢⎣⎡-=011cos 2αA ,则 )2(122211≥⎥⎦⎤⎢⎣⎡==⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡----n D D A D D A D D n n n n n , 通过0=-A E λ,我们可以计算出矩阵A 的ia ia e e -==21λλ,,且21λλ,分别对应T ia T ia e X e X )1(,)1(21,,-==,取⎥⎦⎤⎢⎣⎡==-11),(21ia ia e e X X X ,则100--⎥⎦⎤⎢⎣⎡=X e e X A ia ia, 推出()()⎥⎦⎤⎢⎣⎡-⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡-----ααcos 21cos 40021221X e e X D D n ia n ia n n , 即)0(sin sin )1sin(≠+=αααn D n .例13 设有一个实对称矩阵A ,并且它的阶数是n 阶,满足条件A A =2,且r 为矩阵A 的秩,通过上述条件求出行列式A E -2的值.解 因为A A =2,X X A AX X 22λλ===,所以有0)-(2=X λλ.因为0≠X ,所以0)1-(=λλ,10或=λ.因为矩阵A 是一个n 阶的实对称矩阵,所以它相似于对角矩阵,又因为矩阵A 的秩为r ,所以一定会存在一个可逆矩阵P ,可以使得B E AP P r =⎥⎦⎤⎢⎣⎡=-0001,其中矩阵r E 表示的是r 阶单位矩阵,所以可以推出)(220022211r n E E B E PBP PP A E rn r -==-=-=----.2.9对角化矩阵在其他方面的应用例14 在某个城市的就业数据中显示,一共有30万人从事着不同的三种行业,分别是农业、工业、经商,假设在几年之间这个从业总人数都会保持不变,而且经过整个社会的普查显示:(1)在这个城市的30万人中,投身于农业的有15万人,工业的有9万人,经商的有6万人;(2)在投身于农业的人中,每年大概有%10的人转行去经商,%20的人转行去做工业;(3)在投身于工业的人中,每年大概有%20的人转行去干农业,%10的人转行去经商;(4)在投身于经商的人中,每年大概有%10的人转行去做工业,%10的人转行去干农业.现在请大概预测一下,在未来的一、二年以后,从事这三个行业的人数,以及经历多年以后,从事这三个行业的人员总数会有什么样的一个发展趋势.解 第i 年后还从事这三种行业的人员总数,我们会用一个3维的向量i X 去表示它,则T X )6,9,15(0=.如果想要求21X X ,,并且能够很精确地考察在∞→n 时,n X 的一个发展趋势,那么我们必须要引用一个3阶矩阵)(ij a A =,它的作用是用来体现从事这三种职业人员之间的转移情况.那我们就能够得出矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=8.01.01.01.07.02.01.02.07.0A ,通过矩阵的乘法法则,我们可以得出⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡===-2.79.99.12001AX X A X T ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡===04.823.1073.110212X A AX X , 所以01X A AX X n n n ==-,如果要继续进一步精确地分析n X ,那么必须要事先计算矩阵A 的n 次幂n A ,所以我们先可以将矩阵A 进行对角化,)5.0()7.0()1(8.01.01.01.07.02.01.02.07.0λλλλλλλ---=---=-E A ,所以能够得出特征值5.0,7.0,1321===λλλ,三个特征值分别代表其求出的所对应的三个特征向量321,,q q q ,于是令),,(321q q q Q =,则就会有矩阵1-=QBQ A ,从而推出1-=Q QB A n n ,0X A X n n =,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=5.07.01B ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=n nn B 5.07.01, 当∞→n 时,矩阵n B 将趋向于⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡001,从而推出矩阵n A 将趋向于1001-⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡Q Q , 因为矩阵n X 跟我们已经确定下来的常量*X 非常接近,所以可以得出1-n X 亦必趋于*X ,再通过1-=n n AX X 的转化,就能够准确得知*X 必需要满足条件**AX X =,进而可以推断出*X 是矩阵A 属于特征值11=λ的一个特征向量T T t t t t X ),,()111(*==,,,,303==++t t t t 10=t ,按照上面所讲述的规律转移,经过许多年以后,那么这三种职业的从业人数一定会趋于相等, 三者平均下来为10万人.参考文献[1] 北京大学教学系几何与代数教研室.高等代数(第二版)[M].北京:高等教育出版社,1988.[2] 胡显佑主编.线性代数挚习指导[M].天津:南开大学出版社,1997.[3] 刘九兰,张乃一,曲问薄主编.线性代数考研[M].天津:天津大学出版社,2000.5.[4] 谢国瑞主编.线性代数及应用[M].北京:高等教育出版社.1999.[5] 张学元主编.线性代数能力试题题解[M].武汉:华中理工大学出版社,2000.[6] 徐仲主编.线性代数典型题分析解集[M].西北工业大学出版社,1998,6.[7] 樊辉,钱吉林主编,代数学辞典[M].武汉:华中师范大学出艋社.1994,12.[8] 曹锡皓.高等代数[M].北京:北京师范大学出版社,1987.[9] 张远达.线性代数原理[M].上海:上海科学出版社,1981.[10] Kline Morris. Mathematical Thought from Ancient to Modern Times[M]. New York: OxfordUniversity Press, 1972.[11]Rebollo-Neira L,Fernandez Rubio J.On the Inverse Windowed Fourier transform[M].IEEET rankson Information Theory,1999.[12] Babaie-Zadeh,M. Jutten, C.,Mohimani, H. On the Error of Estimating the Sparsest Solution ofUnderdetermined Linear Systems[M].2011.致谢在开始准备着手写论文到最后定稿的整个过程中,指导教师XXX老师都是非常耐心和细心的引导我和帮助我,在此我向王老师表示由衷的感谢.王老师的严谨治学态度让我受益匪浅.在毕业论文写作的这段时间里,他时时刻刻关心着我的毕业论文的完成情况,并且经常给我指出毕业论文中的缺点与需要改正的地方,最后才能使得我可以顺利完成毕业论文.与此同时,我很感谢所有给过我帮助的老师、同学以及一起努力奋斗过的好朋友.第16 页共16页。

矩阵的对角化及其在高等数学中的应用

矩阵的对角化及其在高等数学中的应用

矩阵的对角化及其在高等数学中的应用矩阵是高等数学中的基础概念之一,它在解决线性方程组和矩阵变换问题中具有重要作用。

在实际问题中,矩阵常常需要进行对角化处理,以便更方便地求解问题。

本文将介绍矩阵的对角化及其在高等数学中的应用。

一、什么是矩阵的对角化对角化是指将一个矩阵变换为对角形式的过程,使得矩阵的主对角线上为非零元素,而其余元素均为零。

举个例子,一个2×2的矩阵A可以进行对角化,其对角化后的形式可以写成:> P^-1 * A * P = D其中P是一个可逆矩阵,D为对角矩阵。

对角矩阵只有主对角线上有非零元素,其他位置都为零。

通过对角化,矩阵变得更加简单,容易处理。

二、如何进行矩阵的对角化对于一个n×n的矩阵A,要进行对角化处理,需要满足以下条件:1.矩阵A必须有n个线性无关的特征向量,这些特征向量组成的矩阵可以写成P=[v1,v2,···,vn]。

2.对于对角矩阵D,其主对角线上的元素必须是矩阵A的n个特征值。

基于这些条件,可以得到矩阵A的对角化公式:> P^-1 * A * P = D其中P=[v1,v2,···,vn],D=[λ1,λ2,···,λn]为对角矩阵。

λ1、λ2···λn为A的特征值,v1、v2···vn为对应的特征向量。

三、高等数学中的应用在高等数学中,矩阵的对角化在求解一些实际问题中具有重要作用。

1. 矩阵的对角化在求解差分方程中的应用线性差分方程是数学中的一种经典问题。

对于一个n阶线性差分方程,其解法是先对其进行离散化处理,变成一个线性方程组。

接着,对该线性方程组进行矩阵形式的表示,就可以得到一个n×n矩阵。

通过矩阵的对角化,可以将线性方程组解放到主对角线上,从而得到差分方程的通解。

2. 矩阵的对角化在离散傅里叶变换中的应用离散傅里叶变换是一种将时域上信号变换为频域上信号的重要算法。

矩阵对角化问题 数学毕业论文

矩阵对角化问题  数学毕业论文

矩阵对角化问题高等代数中,在讲到线性空间和线性变换时,一个主要内容是讨论矩阵对角化,即在什么条件下矩阵与对角矩阵相似.而矩阵对角化的原始问题是:设V 是有限维复线性空间,A 是V 上的线性变换,能否在V 中找到一个基,使得A 在这个基下的矩阵比较简单.作为纯粹的几何问题就是V 能否分解成一些不变子空间的直和.讨论这个几何问题的证明对于了解线性空间有很大好处.本文将对V 分解成所谓根子空间的直和给出一种较为初等的证明,并由根子空间分解定理推出线性变换(或n 阶方阵)可对角化的充要条件.把这些充要条件与其他线性变换(或n 阶方阵)可对角化的充要条件进行汇总比较,从而得到线性变换的矩阵对角化的方法的优劣,便于学习和研究根据具体情况选用.1.预备知识1.1有关定义定义 1.1.1[]1 线性空间V 一个变换A 称为线性变换,如果对于V 中任意的元素αβ和数域P 中任意数K 都有A (α+β)=A (α)+A (β)A (k α)=k A(α) 定义1.1.2[]1 设A 是数域P 上的线性空间V 的线性变换,W 是V 的子空间,如果W 中的向量在A 下像仍在W 中,换句话说,对于W 中任一向量ξ,有A ξ∈W ,我们就称W 是的A 不变子空间,简称A -空间.定义 1.1.3[]1设1V ,2V 线性空间V 的子空间,如果和1V +2V 中每个向α=1α+2α,1122,V V αα∈∈是唯一的,这个和就称为直和.定义1.1.4[]1如果数域P 上的n 阶矩阵A 相似于对角阵,则A 可对角化定义1.1.5[]1设A 是数域P 上的n 阶矩阵,如果数域P 上的多项式()f x 使得()f A = 0,则称()f A 以A 为根.在以A 为根的多项式中次数最低且首相系数为1的多项式称为A 的最小多项式.定义1.1.6 设A 是数域P 上的n 维线性空间V 的线性变换,如果存在非零向量V ξ∈,数λ∈P ,m ∈N,使得()0m A λεξ-=,那么称ξ为属于λ的根向量.线性变换A 的属于特征根λ的根向量的全体,再添上零向量所组成的V 的子集是V 的一个子空间,称V 的这个子空间为A 的属于特征值λ的根子空间.Sylvester 不等式 设,A B 均为n 阶矩阵,秩(A )+秩(B )≤n +秩(AB )1.2 线性空间根子空分解定理引理 设A 是n 维复线性空间V 的线性变换, 12,...s λλλ是A 的所有不同的特征值,且12...s V V V V λλλ=++其中12,,...,s V V V λλλ是V 的全部根子空间,则i A λε-在i V λ上为幂零线性变换,而在1211......i i s V V V V V λλλλλ-+++++++上为可逆线性变换.证明 不失一般性,只证明1A λε-在1V λ上为幂零线性变换,而在23...s V V V λλλ++上为可逆线性变换.在1V λ中取一个基 12,...t γγγ, 则有正整数12,...t p p p ,使1()0i p i A λεγ-= , i = 1,2,…, t ,取p = max {}12,...t p p p , 有()10pi A λεγ-=, i = 1 ,2…t ,于是对任意γ∈1V λ,令1ti i i k γγ==∑,则1()pA λεγ- =1()pA λε-(1ti ii k γ=∑ )=11()0tPii i k A λεγ=-=∑ ,即在1V λ上,1()p A λε- =ϑ (ϑ为零变换) ,所以1A λε-在1V λ上为幂零线性变换.令W =2...s V V λλ++,若1()A W λε-不可逆,则1()A W λε-一定有一个特征根是0 ,因而1A λε-在W 上有属于特征根0 的特征向量0ξ (0ξ∈W) ,即有10()A W λεξ-=1()A λε-0ξ=0, 亦即010()A ξλξ=(0ξ≠0). 又因0ξ∈W = 2...sV V λλ++ ,所以有0ξ=23...s ξξξ++,其中ii V λξ∈ ( i = 2 ,…,s ) 于是有正整数i m ,使()0im i i A λεξ-= , i = 2 ,…,s ,令()()22...s m m s A A τλελε=--,则τ(i ξ) = ()()22...s m ms A A λελε--i ξ= 0 , i = 2 ,…, s ,从而τ(0ξ) = τ(2ξ) + … + τ(ξs) = 0 , 另一方面, 因为()010A ξλξ=,又τ(0ξ)=21()...()s m m s A A λελε--0ξ=()()2121...0s m ms λλλλ--≠这就导致了矛盾.所以1A λε-在2...s V V λλ++ 上为可逆线性变换.定理1.2.1 (根子空间分解定理) 设A 是n 维复线性空间V 的线性变换, 12,...s λλλ是A 的所有不同的特征值,i V λ是属于i λ 的根子空间, i = 1 ,2 ,…, s ,则12...S V V V V λλλ=⊕⊕⊕.证明 设A 的特征多项式为1212()()()...()s s f x x x x γγγλλλ=--- 令()()()ii i f x g x x γλ=- i = 1 ,2 ,…, s , 则12(),(),...,()s g x g x g x 互素, 于是有多项式12(),(),...,()s u x u x u x , 使1()()1si i i g x u x ==∑, 将A 代入上式, 得 1()()si i i g A u A ε==∑,(ε为单位变换), 任给ξ ∈ V ,有ξ =ε(ξ) =()()1s i i i g A u A =⎛⎫⎪⎝⎭∑ξ=1(()())siii g A u A ξ=∑, 记()()i i i g A u A ξξ=, i = 1 ,2 ,…, s ,于是12...s ξξξξ=+++. 下面证明i i V ξ∈ , i = 1 ,2 ,…,s因为()()()i i i f x x g x γλ=-,由哈密尔顿- 凯莱定理()()()i i i A g A f A γλεϑ-== (ϑ为零变换),于是有()i i i A γλεξ-=()()()i i i i A g A u A γλεξϑ-=(ϑ为零变换)即i i V λξ∈, i = 1 ,2,… , s ,所以12...S V V V V λλλ⊂+++,又显然12...S V V V V λλλ⊃+++ ,故12...S V V V V λλλ=+++.再证明上面的和是直和,设12...0,i s i V λαααα++=∈, i = 1 ,2 ,…,s 由引理知i A λε-在i V λ上为幂零变换,所以存在正整数i n ,使得在i V λ上()i n i A λεϑ-=(ϑ为零变换),又由引理 ,i A λε-在111.......i i s V V V V λλλλ-++++++上为可逆变换,所以()i n i A λε- 在111.......i i s V V V V λλλλ-++++++上也是可逆变换,于是0 =()(0)i n i A λε-=()i n i A λε-(12...s ααα++)= ()i n i A λε-i α+()i n i A λε-(1211...i i s ααααα-++++++)=()i n i A λε-(1211...i i s ααααα-++++++)从而1211...i i s ααααα-++++++=0 ,于是()1211......0i i i s αααααα-+=-+++++= , i = 1 ,2 ,… s,由零向量的表法唯一知12...S V V V V λλλ=⊕⊕⊕ 根子空间分解定理全部证完.运用根子空间分解定理可以推出一些矩阵对角化的充要条件.对角矩阵可以认你为是矩阵中最简单的一种,一些复杂的矩阵可以通过适当的方法化为对角阵.通过相应对角阵的研究学习,可以推知这些复杂矩阵的性质,促进对复杂矩阵的了解,简化很多复杂工作,给学习和研究带来很大方便.下面就矩阵对角化的充要条件作一详细论述.2. 矩阵可对角化的一些充要条件及矩阵对角化方法2.1 特征向量法定理2.1.1 设A 是n 维线性空间V 的一个线性变换, A 的矩阵可以在某一组基下为对角阵充要条件是, A 有n 个线性无关的特征向量.证明 设A 在基12,...n εεε下具有对角阵1...n λλ⎛⎫⎪ ⎪ ⎪⎝⎭.即i i i A ελε= i=1,2…n 因此, 12,,...,n εεε就是A 的n 个线性无关的特征向量.反过来,如果A 有n 个线性无关的特征向量,那么就取12,,...,n εεε为基.显然, A 在这组基下的矩阵是对角阵. 证 毕.例1. 设线性变换A 在基12,,...,n εεε下的矩阵是(1)122212221A ⎛⎫ ⎪= ⎪ ⎪⎝⎭, (2)310410482A ⎛⎫⎪=-- ⎪ ⎪--⎝⎭,问A 是否可以对角化? 解 (1)因为特征多项式为122212221E A λλλλ----=------=()()215λλ+-所以A 的特征值是-1(二重)和5把特征值-1代入齐次方程组得()()()123123123122021202210x x x x x x x x x λλλ---=⎧⎪-+--=⎨⎪--+-=⎩ (1)解得基础解系是101⎡⎤⎢⎥⎢⎥⎢⎥-⎣⎦和011⎡⎤⎢⎥⎢⎥⎢⎥-⎣⎦因此属于-1的两个线性无关的特征向量是112223,ξεεξεε=-=-把特征值5代入(1)得基础解系111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,所以属于5的全部特征向量为3123ξεεε=++ 则A 在基123,,ξξξ下的矩阵为B=100010005-⎛⎫ ⎪- ⎪ ⎪⎝⎭(2) E A λ-=310410482λλ--+-=()()212λλ-+,所特征值为1(二重)和-2. 对应特征值1的特征向量为11233620ξεεε=-+ 对应特征值-2的特征向量为23ξε=由此知A 有两个线性无关的特征向量,由定理1知A 不能对角化.运用此定理判定一个线性变换的矩阵是否可以对角化的方法简单易懂,但是过程比较繁琐.先计算一个行列式求出A 的特征值,再利用方程组和特征向量的有关理论及求法计算出A 是否有n 个线性无关的特征向量.计算过程容易出错.下面利用最小多项式给出一个线性变换的矩阵可角化的充要条件.此定理比定理2.1.1简洁实用2.2 最小多项式法引理 设A 是一个对角阵A=12A A ⎛⎫⎪⎝⎭,并设1A ,2A 的最小多项式为12(),()g x g x ,那么A 的最小多项式为12(),()g x g x 的最小公倍数[]12(),()g x g x .证明 ()g x =[]12(),()g x g x ,首先12()()()g A g A g A ⎛⎫= ⎪⎝⎭=0.因此()g x 能被A 的最小多项式整除.其次()0h A =.那么12()()()h A h A h A ⎛⎫= ⎪⎝⎭=0, 1()h A =0,2()h A =0,因而11()()g x h x ,22()()g x h x .并由此得()()g x h x .这样就证明了()g x 是A 的最小多项式. 这个结论可以推广到A 为若干矩阵组成的准对角阵的情形.即如果A=1 (00)S A A ⎛⎫ ⎪⎪ ⎪⎝⎭,i A 的最小多项式为()i g x ,i=1,2,…,s.那么A 的最小多项式为[]12(),(),...,()s g x g x g x .定理2.2.1 数域P 上n 级矩阵A 与对角阵相似的充要条件为A 的最小多项式是P 上互素的一次因式的乘积.证明 根据引理的推广形式,条件的必要性是显然的. 下面证明充分性.根据矩阵和线性变换之间的关系,我们可以定义任意线性变换A 的最小多项式,它等于其对应矩阵A 的最小多项式.所以只需证明,若数域P 上某线性空间V 的线性变换A 的最小多项式()g x 是P 上互素的一次因式的乘积1()()li i g x x a ==-∏,则A 有一组特征向量做成V的基.实际上,由于()0g A V =.由定理 1.2.1同样的步骤可证12...l V V V V =⊕⊕⊕,其中{}()0,i i V A a V ξεξξ=-=∈,把12,...l V V V 各自的基合起来就是V 的基,而每个基向量都属于某个i V ,因而是A 的特征向量. 证毕.推论 复数矩阵A 与对角阵相似的充要条件是A 的最小多项式无重根. 不利用定理2.2.1,该推论也可证明.下面给出令一种证明.证明 必要性设A 相似diag 12(,...)n λλλ,所以存在可逆矩阵T 使1T AT -=∧,(∧为对角阵),从而1i i T A T -=∧,不妨12,...k λλλ是A 的互不相同的特征根()k n ≤ 记()()()11211()......k k k k k g a a a λλλλλλλλλλ--=---=+++ 因而()11111(...)k k k k T g A T T A a A a A a E T----=+++=1111111...k k k k T A T a T A T a T AT a T ET ------+++=11...k k k a a E -∧+∧++=()g ∧ 而()11...k k k g a a E -∧=∧+∧++=1111211121(,...)(,...)...(,...)k k k k k k n n k k k diag diag a a a diag a a a λλλλλλ---++= 11111.........k k k k k n n k a a a a λλλλ-⎛⎫++ ⎪⎪ ⎪++⎝⎭=diag ()()()12(),...n g g g λλλ=0所以()g A =0.于是()()A m g λλ,但是()g λ没有重根,因而()A m λ没有重根.充分性 设12,...n λλλ为最小多项式()A m λ的互不相同的根,则由()A m λ无重根()A m λ=()()()12...k λλλλλλ---,于是()A m A =()()()12...k A E A E A E λλλ---=0 令rank ()i A E λ-=i γ,则dim I V λ=n -i γ,所以A 共有()()()12...k n n n s γγγ---=个线性无关的特征向量并且显然s n ≤.另一方面()12...1k k n γγγ+++≤-.因而又有()()()12...k s n n n n γγγ=---≥,故s n =.这就说明了A 有n 个线性无关的特征向量由定理2.1.1知A 可对角化. 证毕.例2. 判下列矩阵是否可以对角化.(1)001010110⎛⎫ ⎪ ⎪ ⎪⎝⎭(2)3131131331311313--⎛⎫⎪-- ⎪⎪-- ⎪--⎝⎭解(1)可求的A 的特征多项式为()()2010101110E A λλλλλλ--=-=-+-由于A 的最小多项式为()()211λλ-+的因式,计算得0A E -≠,0A E +≠.而()A E -()A E +=0.因此A 的最小多项式为()()11λλ-+.显然A 的最小多项式是实数域上互素的一次因式的乘积,从而由定理2.2.1知A 可对角化.(2)可求得A 的最小多项式为E A λ-=3131131331311313λλλλ-----+---+=4λ由于的最小多项式为4λ的因式,计算得A 0≠, 2A =0.因此A 的最小多项式为2λ.从而由定理2.2.1知A 不可对角化.例3 k A =E,则A 与对角阵相似.(k=1,2…)证明 由k A E =知A 为多项式()1k f x λ=-的零点,即()f A =0.因A 的最小多项式()()A m f λλ,而()f λ无重根,所以()A m λ无重根,故由推论知A 与对角阵相似.对于单纯的判断一个线性变换的矩阵能否对角化运用定理 2.2.1及其推论是很简洁方便的,它部避免了运用定理2.1.1的繁琐过程.但是对于既要判定某个数域上的线性变换的矩阵是否可对角化,对于可对角化的矩阵又要求出相似变换矩阵及矩阵特征值的题目来说运用定理2.2.1及推论是达不到要求的.而运用定理2.1.1虽然能达到要求但方法却很繁琐.下面给出的方法仅需利用矩阵的乘法运算便可判定一个矩阵是否相似与对角阵,并且在判定的过程中简洁的构造出相似变换矩阵完全不需解性方程组.2.3 矩阵的乘法运算法定理 2.3.1 设12,,...,s λλλ为n 阶矩阵A 的全部相异特征值,其重数分别为12,,...,s n n n ,1sii nn ==∑,则A 与对角阵相似的充要条件是1()si i E A λ=-∏=0.(i=1,2,…,s)证明 必要性若A 相似于阵对角阵∧,则存在可逆矩阵P 使得A =P 1...s E E λλ⎛⎫⎪ ⎪ ⎪⎝⎭1P -,其中i E 为in 阶单位矩阵(i=1,2,…,s )于是()i E A λ-=()1i P E P λ--∧=()()111...i i s s E P P E λλλλ--⎛⎫⎪ ⎪ ⎪-⎝⎭,于是1()s i i E A λ=-∏=()11s i i P E P λ-=-∧∏= P ()()1111...s i i si s s i E E λλλλ==⎛⎫- ⎪ ⎪⎪ ⎪ ⎪- ⎪⎝⎭∏∏1P - 由于()i j j iE λλ-∏=0(j=1,2,…,s ).所以1()si i E A λ=-∏=0.充分性 因为对于任何n 阶矩阵A 都存在可逆矩阵P ,使得A= P 12...S J J J ⎛⎫ ⎪⎪ ⎪ ⎪⎝⎭1P -,其中jJ 为jordan 块(j=1,2,...,s).因此要证A 可对角化,只要证j J =j j E λ(j=1,2,…,s ),由于()i E A λ-=()i P E J λ-1P -=P ()()()1122...i i i s s E J E J E j λλλ-⎛⎫ ⎪-⎪⎪ ⎪⎪-⎝⎭1P - ()()()()111221...i i i i i i s s i E J E J E A P P E J λλλλ-⎛⎫- ⎪ ⎪- ⎪-= ⎪ ⎪⎪-⎪⎝⎭∏∏∏∏所以若()0i iE A λ-=∏.则因P 可逆而有()0i j j iE J λ-=∏(j=1,2,…,s ).又当i j ≠时()0ijλλ-≠,()i jj EJ λ-可逆,所以()i j j E J λ-0≠,即j j j J E λ=(j=1,2,…,s)定理2.3.2 设12,,...,s λλλ时n 阶矩阵的全部相异特征根,其重数分别为12,...s n n n ,则A 于对角阵相似的充要条件是()j i i jW E A λ≠=-∏的秩为()j j R W n =(j=1,2,…,s ).证明 必要性()()()111...i i j j i i ji s S i E W E A P P E λλλλλ≠-≠⎛⎫- ⎪⎪=-=⎪ ⎪- ⎪⎝⎭∏∏∏=()110...0ijji js P P Eλλ-≠⎛⎫⎪ ⎪⎪- ⎪ ⎪ ⎪⎝⎭∏ 其中0,j j E 分别是j n 阶的零矩阵和单位矩阵(j=1,2,…,s ).由于P 满秩且i j λλ≠.所以()j R W =()i j j i j R E λλ≠⎛⎫- ⎪⎝⎭∏=()j j R E n =.充分性 用反证法假设()i j j i j R E λλ≠⎛⎫- ⎪⎝⎭∏不可对角化,则因几何重数≤代数重数[]5,必至少存在一整数k 使得()k R E A λ->()j R E []3,于是j k ≠时.由sylvester 不等式知j n =()k i j R E A λ≠⎛⎫- ⎪⎝⎭∏≥()()2i i jR E A s n λ≠---∑>()()2j i j n n s n ≠---∑=()()()12i j j i js n n s n n n n n ≠-=--=--=∑矛盾.所以A 可对角化.推论 1 设12,,...,s λλλ为n 阶矩阵A 的相异的特征根,其重数为12,,...,s n n n ,则矩阵j W =()k i jE A λ≠-∏的列向量中由对应于j λ的j n 个线性无关的特征向量.证明 因A 可对角化,由定理2.3.1得()i i jE A λ≠-∏=0,()jE A λ-()i i jE A λ≠-∏=()j jE A Wλ-=0.由此,j W 中每一列非零向量都是方程组()i E A λ-X=0解向量,即j λ的特征向量.又有定理2.3.2知()j j R W n =,所以j W 的列向量组中有恰好对应于j λ的j n 个线性无关的特征向量.上述的结论表明,要构造可对角化矩阵 A 的相似变换矩阵P ,完全可以不像传统的方法那样解方程组()k E A λ-X=0,而只需对每一特征值j λ(j=1,2,…,s )从矩阵乘积()ki jE A λ≠-∏中直接找出jn个与j λ对应的线性无关的特征向量,这样所得的j n n =∑个特征向量为列作一n 阶矩阵即可.推论2 若n 阶可对角化矩阵A 只有两个相异特征值1λ(k 重)和2λ(n k -重),则矩阵()1E A λ-(或()2E A λ-的n k - (或k )个线性无关列向量就是对应2λ(或1λ)的特征向量的极大无关组.这一结论进一步表明,在可对角化矩阵A 只有2个相异特征值的情况下,不仅不需要解方程组,而且不需要计算矩阵的乘积就可以把对应于不同特征值的特征向量立即求出.例4 求下列矩阵A 相似变换矩阵.(1)A =741471444-⎛⎫ ⎪ ⎪ ⎪-⎝⎭ (2)A =1220212022100001⎛⎫⎪⎪ ⎪ ⎪⎝⎭解 (1)A 的特征值1λ=12,2λ=3(二重)21541451448W E A λ-⎛⎫⎪=-=- ⎪ ⎪⎝⎭,124441W E A λ-⎛⎫⎪=-= ⎪ ⎪-⎝⎭-41-4-41由于()()120E A E A λλ--=,所以A 可对角化,有推论2知1λ的一个特征向量()11,1,1α=-(取1W 的第3列)2λ的2个线性无关的特征向量()()234,5,4,1,1,8αα=-故相似变换矩阵P =()123,,ααα=141151148-⎛⎫⎪⎪ ⎪-⎝⎭,1(12,3,3)PAP diag -=(2)A 的特征值1λ=-1(二重),2λ=5,3λ=1,而()()123W E A E A λλ=--=8448*4400-⎛⎫⎪-⎪ ⎪-- ⎪⎝⎭ ()()1300*08E A E A λλ⎛⎫ ⎪ ⎪--= ⎪ ⎪⎝⎭2W=,3W=88*80⎛⎫⎪ ⎪ ⎪⎪⎝⎭由推论2可得1λ的特征向量()()1284404αα''=--=-,,,,,8,-4,0. 23λλ,的特征向量分别为()()3400088880αα==,,,,,,,于是相似变换矩阵为P=()1234αααα,,,==84008440000-⎛⎫⎪⎪⎪--⎪⎝⎭8-488-8 P A 1P -=diag(-1,-1,5,-1).上文讨论了矩阵是否可对角化的判定及矩阵对角化方法问题,给出了简便易行的判定和求法.区别于传统的方法,定理2.3.1定理2.3.2及推论把矩阵对角化问题归结为矩阵的乘法运算,不需要解方程组就可以得到特征向量及相似变换矩阵,但是上述方法都没有达到特征值,特征向量,相似变换矩阵同步求解的效果.下面引入λ-矩阵,改进在一般情形下矩阵对角化的方法,使判定和求解一步到位并得到矩阵对角化十分简单的方法,主要依据下面两个定理.2.4 引入λ-矩阵推出数字矩阵可对角化的充要条件定理2.4.1 设A 是数域P 上的n 阶方阵,()()A E A λλ=-为其特征矩阵E 为n 阶单位阵.如果()A λ经过初等变换化为对角阵()D λ,则A 的特征值为()D λ的对角线上元素的乘积的多项式的根. (证明略)定理2.4.2 在定理2.4.1 的假设下,如果()()(),T A D λλ经初等变换化为()()(),D P λλ,且()D λ为对角阵,则(1) 对于A 的每个特征值i λ,()i P λ中与()i D λ的零行对应的行向量生成属于i λ的特征子空间.(2) 若A 的特征值都在P 内,设12,,...,s λλλ为A 的全部不同的特征值,其重数分别为12,,...,s γγγ,则A 可以对角化的充要条件是()i D λ中零行的数目=i λ的重数i γ(i=1,2,…,s )证明 (1)因为()D λ与()T A λ的秩为n ,则总有可逆的λ-矩阵()P λ,()Q λ,使()()()()()()()12(,,...,)T n P A Q diag d d d D λλλλλλλ==.即对()T A λ施行()P λ对应的一些行初等变换和()Q λ对应的一些列初等变换可使()T A λ化为对角阵()D λ,有()()()(),T P A Q E λλλ→()()(),D P λλ (1) 这里相当于初等列变换的()Q λ右乘作用在()T A λ而不作用于E.因为()()()T P A Q λλλ=()D λ,所以()()()()T T T Q A P D λλλλ==()D λ.于是对A 的每个特征值i λ有()()()T T i i i Q A P λλλ=diag(()()()12,,...,i i n i d d d λλλ)设()i D λ中有i m 个零行,相应的i m 个为的对角元记为()()()12...0i i i i i im i d d d λλλ====()1i m n ≤≤,取()T i P λ中对应的列向量1,2,...,i i i im P P P ,则()()T i i Q E A λλ- ()1,2,...,i i i im P P P =0.因为()T i Q λ可逆,所以()i E A λ- ()1,2,...,i i i im P P P =0 (2)由于()T i P λ可逆,故()12,,...,i T T T i i im P P P 列满秩,从而由(2)知12,,...,i T T Ti i im P P P 正是A 属于iλ的i m 个线性无关的特征向量,再从(1)式,注意到()i D λ中n -i m 个非零行是行满秩的.由[]7中定理1知A 属于i λ的线性无关的特征向量就是()i P λ中与()i D λ的零行对应的行向量,他们生成i λ对应的特征子空间.(2) A 可对角化⇔秩()i E A λ-=i n γ-=i n m -,即i m =i γ(i=1,2,…,s ) 证毕. 基于以上讨论我们不难得到矩阵对角化的简单方法,其步骤如下: (1)对(),TiE A E λ-作初等变换化为()()(),D P λλ,其中()()()()12(,...,)n D diag d d d λλλλ=,,则A 的特征值恰是()()()12...n d d d λλλ=0的根. (2) 如果A 的特征向量全在P 内,且对每个i λ有()i D λ中零行数目=i λ的重数,则A 可以对角化,否则不可对角化.(3) 对于每个i λ,在()i P λ中取出与()i D λ中零行对应的行向量12,,...,i i i im P P P 得A 属于i λ线性无关的特征向量.(4) 若A 可以对角化,作可逆矩阵()1121,,...,,...,,...,si i im s sm T P P P P P =,则11122(,,...,)s s T AT diag E E E λλλ-=,i E 为i γ阶矩阵.例5 判定下列矩阵可否对角化,若可以求可逆矩阵T ,使1T AT -为对角阵.(1) A =011111011⎛⎫ ⎪- ⎪ ⎪⎝⎭ (2) A =321222361-⎛⎫ ⎪-- ⎪ ⎪-⎝⎭解 ()()10100,111010011100T A E λλλλ-⎛⎫⎪=-- ⎪ ⎪--⎝⎭→2011002011111001λλλλλλλ⎛⎫-- ⎪--- ⎪ ⎪--⎝⎭ →21110010111102011λλλλλ--⎛⎫⎪-+ ⎪ ⎪---⎝⎭→210000101111020011λλλ-⎛⎫ ⎪-+ ⎪ ⎪--⎝⎭→()22100001010111001231λλλλλλλ⎛⎫-⎪-+ ⎪ ⎪⎪---++-++⎝⎭故P 的特征值是120,1λλ==(二重),因()1D 中的零行数目2λ≠的重数,故P 不可对角化.(2)()()2323100121001,22601002240121210010242103T A E λλλλλλλλλλλ⎛⎫----+⎛⎫⎪ ⎪=-+-→-- ⎪ ⎪ ⎪ ⎪-+-+-+⎝⎭⎝⎭()()2100001100001022*******012024241030242103λλλλλλλλλλλ⎛⎫⎛⎫ ⎪ ⎪→--→- ⎪ ⎪ ⎪ ⎪---+-+--+-+⎝⎭⎝⎭()()()1000010200120024121λλλλ⎛⎫⎪→- ⎪ ⎪--+--+⎝⎭故A 的特征值为12λ=(2重根), 24λ=-.又()2D 中零行数=2=1λ的重数;()4D -的零行数=1=2λ的重数,故P 可对角化,且由()()()2,2D P =100001000012000123⎛⎫ ⎪⎪ ⎪--⎝⎭可得出()()αβ''=0,1,2,=1,-2,-3是A 属于2的线性无关特征向量由()()()4D -,P-4=100001060012000123⎛⎫ ⎪- ⎪ ⎪-⎝⎭得()1,2,3γ'=-是A 属于-4的线性无关的特征向量.令T=011122233⎛⎫ ⎪-- ⎪ ⎪-⎝⎭,则1224T AT -⎛⎫⎪= ⎪ ⎪-⎝⎭参考文献[1] 北京大学数学系.高等代数.北京:高等教育出版社,第88版,1988.[2] 许以超.代数学引论[]M.上海:社会科学技术出版社,1966[3] 钱吉林.矩阵及其广义矩阵[]M.武汉:华中师范大学出版社.[4] 王心介.高等代数与解析几何[]M.北京:科学出版社,2002.[5] 张远达.线性代数原理.上海:上海教育出版社,1980.[6] 彭海明.对“矩阵特征值与特征向量同步求解方法探讨”的改进意见[]J.数学通报,1993(2):45-47.[7] 刘国洪.王宝智.利用矩阵的初等行变换对矩阵的特征值和特征向量同步求解,数学通报,1996,2.。

空间解析几何的对角化对角化与相似矩阵的计算与应用

空间解析几何的对角化对角化与相似矩阵的计算与应用

空间解析几何的对角化对角化与相似矩阵的计算与应用在线性代数中,对角化是一个重要的概念,它在几何学和矩阵计算中具有广泛的应用。

本文将探讨空间解析几何中的对角化概念及其与相似矩阵的计算和实际应用。

一、对角化的概念对角化是指将一个矩阵通过相似变换变为对角阵的过程。

在空间解析几何中,对角化可以帮助我们更好地理解和描述几何物体的性质和运动规律。

对于一个n阶矩阵A,如果存在一个可逆矩阵P,使得P^{-1}AP=D,其中D为对角阵,那么我们称矩阵A可对角化,矩阵D的主对角线上的元素即为矩阵A的特征值。

二、对角化的计算方法对角化的计算可以通过求解矩阵的特征值和特征向量来实现。

具体步骤如下:1.计算特征值:解方程|A-\lambda I|=0,其中A为给定矩阵,\lambda为标量,I为单位矩阵。

解该方程可得到矩阵A的特征值。

2.计算特征向量:将每个特征值代入方程(A-\lambdaI)\mathbf{X}=\mathbf{0},其中\mathbf{X}为特征向量。

解该齐次线性方程组即可得到特征向量。

3.构造P矩阵:将特征向量按列排成一个矩阵P,即P=[\mathbf{X_1},\mathbf{X_2},\cdots,\mathbf{X_n}],其中\mathbf{X_i}为第i个特征向量。

4.计算相似矩阵:利用矩阵P和对角阵D,可以得到相似矩阵A=PDP^{-1}。

三、对角化的应用对角化在空间解析几何中有许多应用,以下列举几个常见的应用。

1.求解线性方程组:对角化可以将一个线性方程组转化为简化的形式。

利用相似矩阵的性质,对角阵的求解相对更加简单。

通过对角化,可以更快速地求解线性方程组。

2.描述几何变换:对角化可以帮助我们描述和研究几何物体的旋转、缩放和平移等变换。

通过对角化后得到的相似矩阵,可以直观地了解几何物体的变换规律。

3.优化问题:在优化问题中,对角化可以将问题转化为更简单的形式。

例如,对角化可以将一个二次函数实现坐标的旋转和缩放,从而更便于求解最优解。

【论文】矩阵对角化方法及相关应用开题报告

【论文】矩阵对角化方法及相关应用开题报告

【关键字】论文
毕业论文开题报告
数理系数学与应用数学专业 2012 级 1 班
课题名称:矩阵对角化方法及相关应用
毕业论文起止时间:
年月日~月日(共周)
学生姓名:丁潞泷学号:
指导教师:黄斌
报告日期: 2012年6月25日
说明:
1.本报告必须由承担毕业论文课题任务的学生在接到“毕业论文任务书”、正式开始做毕业论文的第2
周或第3周末之前独立撰写完成,并交指导教师审阅。

2.每个毕业论文课题撰写本报告一份,作为指导教师、教研室主任审查学生能否承担该毕业论文课题任
务的依据,并接受学校的抽查。

此文档是由网络收集并进行重新排版整理.word可编辑版本!。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

矩阵对角化及应用理学院 数学082 缪仁东 指导师:陈巧云摘 要:本文是关于矩阵对角化问题的初步研究,对矩阵对角化充要条件的归纳,总结,通过对实对称矩阵,循环矩阵,特殊矩阵对角化方法的计算和研究,让读者对矩阵对角化问题中求特征值、特征向量,求可逆矩阵,使对角化,提供了简便,快捷的求解途征.关键词:对角矩阵;矩阵对角化;实对称矩阵;特征值;特征向量.矩阵对角化是矩阵论的重要组成部分,在矩阵论中占有重要的作用,研究矩阵对角化问题很有实用价值,关于矩阵对角化问题的研究,这方面的资料和理论已经很多.但是他们研究的角度和方法只是某个方面的研究,没有进行系统的分类归纳和总结.因此,我就针对这方面进行系统的分类归纳和总结,对一些理论进行应用和举例,给出算法.特别给出了解题时方法的选择.1.矩阵对角化概念及其判定所有非主对角线元素全等于零的n 阶矩阵,称为对角矩阵或称为对角方阵.定义1.1 矩阵A 是数域P 上的一个n 级方阵. 如果存在一个P 上的n 级可逆矩阵X ,使1X AX - 为对角矩阵,则称矩阵A 可对角化.矩阵能否对角化与矩阵的特征值特征向量密切相关.定义 1.2 设A 是一个n 阶方阵,λ是一个数,如果方程组AX X λ= (1)存在非零解向量,则称λ为的A 一个特征值,相应的非零解向量X 称为属于特征值λ的特征向量.(1)式也可写成,()0E A X λ-= (2)这是n 个未知数n 个方程的齐次线性方程组,它有非零解的充分必要条件是系数行列式=0E A λ-, (3)即1112121222120n nn n nna a a a a a a a a λλλ------=---上式是以λ为未知数的一元n 次方程,称为方阵A 的特征方程. 其左端A E λ-是λ的n 次多项式,记作()f λ,称为方阵的特征多项式.111212122212()||n nA n n nna a a a a a f E A a a a λλλλλ------=-=---111n n n n a a a λλλ--=++++显然,A 的特征值就是特征方程的解.特征方程在复数范围内恒有解,其个数为方程的次数(重根按重数计算),因此,n 阶矩阵A 有n 个特征值.设n 阶矩阵()ij A a =的特征值为12,,n λλλ,由多项式的根与系数之间的关系,不难证明(ⅰ)121122n nn a a a λλλ+++=+++;(ⅱ)12n A λλλ=.若λ为A 的一个特征值,则λ一定是方程=0A E λ-的根, 因此又称特征根,若λ为方程=0A E λ-的i n 重根,则λ称为A 的i n 重特征根.方程 ()0A E X λ-=的每一个非零解向量都是相应于λ的特征向量,于是我们可以得到求矩阵A 的全部特征值和特征向量的方法如下: 第一步:计算A 的特征多项式E A λ-;第二步:求出特征方程=0E A λ-的全部根,即为A 的全部特征值;第三步:对于的每一个特征值λ,求出齐次线性方程组:()0E A X λ-= 的一个基础解系12,,,s ξξξ,则A 的属于特征值λ的全部特征向量是 1122s s k k k ξξξ+++(其中12,,,s k k k 是不全为零的任意实数).设P 是数域, Mn (P ) 是P 上n ×n 矩阵构成的线性空间, A ∈Mn (P ) , 1,2t ,,λλλ 为A 的t 个互不相同的特征值,高等代数第二版(北京大学数学系几何与代数教研室编)第四版(张和瑞、郝炳新编)课程中,我们学过了矩阵可对角化的若干充要条件如: (1) A 可对角化当且仅当A 有n 个线性无关的特征向量; (2) A 可对角化当且仅当特征子空间维数之和为n ; (3) A 可对角化当且仅当A 的初等因子是一次的; (4) A 可对角化当且仅当A 的最小多项式无重根我们知道线性变换A 的特征多项式为f (λ) ,它可分解成一次因式的乘积1212()()()()i r r r i f λλλλλλλ=---则V 可分解成不变子空间的直和其中i V = {ξ|iri 12-==s V V V V λ⊕⊕⊕(A E );ξ∈V}引理 1.1:设A, B 都是n 阶矩阵, 则秩( AB) ≥秩( A) + 秩( B) - n.定理 1.1:设A 是实数域F 上的一个n 阶矩阵, A 的特征根全在F 内, 若1λ, 2λ,...,K λ 是A 的全部不同的特征根, 其重数分别为1r , 2r ,... k r , 那么 (Ⅰ) 可对角化的充要条件是()i j i jE A r λ≠⎛⎫-= ⎪⎝⎭∏秩 j=1, 2,.......k(Ⅱ) 当( 1) 式成立时,()ii jE A λ≠-∏ 的列空间就是A 的属于特征根iλ的特征子子空间.证明: (Ⅰ) 设A 可对角化, 则存在可逆阵T, 使{}11122,,...,k K T AT diag E E E λλλ-=这里右边是分块对角矩阵, j E 为i r 阶单位阵, 于是有()()()11i i i i j i j i j E A T E A T E T AT λλλ--≠≠≠⎛⎫⎛⎫⎛⎫⎛⎫-=-=- ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭∏∏∏秩秩秩={}()122,,...,i K K i j E diag E E E λλλλ≠⎛⎫-⎪⎝⎭∏秩=()()(){}12,,...,,i j i j i j Ki j diag E E E λλλλλλ≠⎛⎫---⎪⎝⎭∏秩 =()0,0,...0,,0,0,...,0i j j j i jdiag E r λλ≠⎛⎫⎧⎫-= ⎪⎨⎬ ⎪⎩⎭⎝⎭∏秩 j=1,2, ......k.反之,若()()ijE A r λ-=∏秩i=1,2,.....k, 反复用引理可得()()()()()22i j i i i ji jE A E A K n n r k n λλ≠≠-≥---≥---∑∑∏秩r 秩 i j i jn r r ≠=-=∑ j=1,2,...,k.这里用到了齐次线性方程组()0i E A X λ-=的解空间的维数不大于i λ的重数不大于j r 这个结论.于是又()()iii j i jE A n r λ≠≠-=-∑∑秩从而()i iA n r λ-=-秩 i=1,2,......k. 这样的矩阵可以对角化.(Ⅱ)设( Ⅰ)式成立,则A 可对角化.故A 的最小多项式为()1kii x λ=-∏从而()10kii E A λ=-=∏ 即 ()()0i ii jE A E A λλ≠--=∏这就是说,列空间包含在i λ的特征子空间中,但是由(1), ()ii jE A λ≠-∏的列空间的维数是n,它正是j r 的特征子空间的维数,所以结论(Ⅱ) 成立.推论: 设A 为实数域F 上的n 阶矩阵,A 的特征根全为F 内,且1λ, 2λ 是A 的全部不同的特征根, 其维数分别为1r , 2r , 若秩()12E A r λ-=,秩()21E A r λ-=,则A 可以对角化,且()E A λ-的列向量组的极大无关组恰是属于2λ 的极大线性无关的特征向量组,2E A λ-的列向量组的极大无关组恰是属于1λ的极大无关的特征向量组.例1: 判断A=460350361⎛⎫⎪-- ⎪ ⎪-⎝⎭能否对角化,并求特征向量.解: 易知A 的特征根1λ =-2 , 2λ =1.1E A λ- =660350363--⎛⎫ ⎪ ⎪ ⎪--⎝⎭ 和2E A λ- =360360360--⎛⎫⎪⎪ ⎪⎝⎭的秩分别为2与1,故A 可对角化. 又因为可以选取001⎛⎫ ⎪ ⎪⎪⎝⎭和210-⎛⎫⎪⎪ ⎪⎝⎭为的列空间的一个基,111⎛⎫ ⎪- ⎪ ⎪-⎝⎭是属于1λ的特征向量.定理和推论把判断矩阵是否对角化的问题与求它的特征向量的问题联系起来,给出了一个不用解线性方程组而求得可对角化矩阵的特征向量的方法, 在矩阵的不同特征根较少时, 这个方法较方便.2.实对称矩阵对角化的计算方法我们知道任意实对称矩阵,总正交相似于一对角阵. 该对角阵的对角元即为实对称矩阵的特征值, 正交相似变换矩阵的各列构成相应的特征向量. 给定一实对称阵A ,如何求正交相似变换矩阵P ,使1T P AP PAP -=为对角阵. 理论上的解决方法为:首先利用特征方程: | λI - A | = 0 求出全部特征值,针对不同特征值求出相应的完全特征向量系,合在一起构成实对称阵A 的完全特征向量系. 再利用施密特正交化法得到 A 的规范化正交特征向量系. 以此作为列向量得到正交相似变换矩阵P , 1T PAP PAP -=为对角阵, 参见文献[5 ]. 此方法理论可行,但在具体操作时,由于要事先求出实对称阵A 的全部特征值,操作上有如下困难: (1) 特征方程: | λI- A | = 0 给出困难; (2) 特征方程求根困难(5 次以上的代数方程没有统一的求根公式) . 因此有必要寻求方法.定义2.1 (瑞雷商) 设A 为n 阶实对称阵,对于任一n 维非零列向量x ,称R ( x) =( A x , x)/( x , x) 为关于向量x 的瑞雷商.引理2.1 设A 为n 阶实对称阵, 1λ≥2λ≥......≥n λ 为A 的特征值.()()()()11/{0}/{0},,max ,min,,nnx R x R Ax x Ax x x x x x λλ∈∈== 定义2.2 设w 为n 维列向量,且T w w = 1 ,则n 阶矩阵H = I - 2Tww 称为Householder 阵.引理2.2 Householder 矩阵具有如下性质: (1) TH H =(2) T TH H HH I == ( H 是正交阵) .引理2.3 设x , y ∈nR , x ≠y , X Y =,则存在Householder 矩阵H, 使Hx = y. 其中()()22/TH I x y x y x y =----定理2.1 设A 是实对称矩阵,λ, x (2X= 1) 是A 的一个特征值和相应的特征向量,则存在P 为一个正交阵,使Px =1e = ()1,0,0 0. 且TPAP 的第一行和第一列的第一个元素为λ,其余元素均为零.证 设A 是实对称矩阵, 1λ≥ 2λ≥ ...≥ n λ为A 的特征值. 根据引理2.1 ,利用多元函数求极值的拉格朗日乘数法,可求得1λ 及相应的规范化特征向量1X . 不妨假设‖1X ‖ = 1 ,由引理2.3 ,存在1P 为一个正交阵,使11P X =1e =()1,0,0, 0.且TPAP 的第一行和第一列的第一个元素为1λ , 其余元素均为零. 设111100TP AP A λ⎛⎫=⎪⎝⎭, 为对称阵,故1A 也为对称阵,设2λ 及2X 为1A 最大特征值及相应的规范化特征向量,则根据引理2.3 ,存在2Q 为一个正交阵,使()2211,0,0, 0Q x e ==.且212T Q A Q 的第一行和第一列除2λ 外其余元素均为零. 令22100P Q ⎛⎫= ⎪⎝⎭,容易验证2P 亦为正交阵, 满足:1121122212200000000T TT P P AP P Q AQ A λλλ⎛⎫⎛⎫⎪==⎪ ⎪⎝⎭ ⎪⎝⎭依此类推, 存在正交阵1p ,2p , ⋯,1n p -, 使得1n p -...2p 1p 121...T T Tn Ap p p D -=,则T PAP =D,其中 D 为对角阵,令121P P P P n -=,则TPAP D =,P 即为将实对称阵对角化的正交相似变换矩阵.例2: 设矩阵210210582811A ⎛⎫⎪=- ⎪ ⎪-⎝⎭, 1λ≥2λ≥3λ为A 的特征值.按上面的算法进行对角化,求出正交矩阵P 及特征根和特征向量.解: (1)利用瑞雷商和多元函数求极值的拉格朗日乘数法,可求得1λ = 18 ,相应的特征向量为1122,,333Tx ⎛⎫=- ⎪⎝⎭(2) 计算正交矩1p =()()211112/Tp I x e x e x e =----=122333221333212333⎛⎫- ⎪ ⎪ ⎪-- ⎪ ⎪ ⎪--- ⎪⎝⎭,满足()1111,0,0T p x e ==且111800090009TP AP ⎛⎫⎪=- ⎪ ⎪⎝⎭,至此已实现对角化. 借此可求得= 2λ=9 , 3λ = - 9. 相应的特征向量分别为2212,,333Tx ⎛⎫=--- ⎪⎝⎭,3221,,333Tx ⎛⎫=-- ⎪⎝⎭.3.循环矩阵对角化方法的研究在复数域 C 上,形如012110121230........................n n n a a a a a a a a A a a a a ---⎛⎫⎪⎪= ⎪⎪⎝⎭的矩阵,称关于元素列011,,...,n a a a -的循环矩阵.已知n 阶循环矩阵010 (00)01...0 (1)00...0K ⎛⎫⎪⎪= ⎪⎪⎝⎭,并令ii K K = (1,2,,)i n =,称121,,,....,n E K K K -为循环矩阵基本列(其中E = n K 为单位矩阵).循环矩阵基本列有如下特点: ①121,,,...,n E K K K -都是循环矩阵;②n i i K K += ,即n i iK K +=;③n 阶循环矩阵K 有n 个特征根: cossinm mx mxi n nλ=+ (0,1,,1)m n =-④关于元素列0121,,,...,n a a a a -的n 阶循环矩阵 A 可用循环矩阵基本列表示为210121...n n A a E a K a K a K --=++++,反之,能用循环矩阵基本列线性表示的矩阵,则一定是循环矩阵. 循环矩阵的性质性质1 同阶循环矩阵的和矩阵为循环矩阵. 性质2 同阶循环矩阵的乘积满足交换律.性质3 同阶循环矩阵的乘积为循环矩阵. 性质4 循环矩阵的逆矩阵为循环矩阵.n 阶矩阵A 关于多项式函数f (x) 生成的矩阵为f (A) ,A 的特征根与f (A) 的特征根有下面的结论:命题3.1 设f (x) 是一个n - 1 次多项式函数,若λ是矩阵A 的特征根,则 f (λ) 是矩阵f (A) 的特征根.命题3.2 设f (x) 是一个n - 1 次多项式函数,若矩阵A 相似于矩阵B , 则矩f (A) 相似于矩阵f (B) .考察n 阶循环矩阵K,K 的特征多项式为:()211,(n i njjnj E K ei πλλληη-=-=-=-==∏如果n 阶循环矩阵A 记为()210121...n A n A f K a E a K a K a K --==++++不难求得K 中与特征值j η相应的特征向量,记:()11...j j n x ηη-⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦, ()()22......11j j j j j j j j kx x ηηηηηη⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦则由命题3.1得()()()()()jjj j A A Ax f K x f x η==,可以验证()()()()1111000,1,.11,1n n m kmkk k m xxm k mηη---==≠⎧==-=⎨=⎩∑∑.将这n 个两两正交的向量()j x 单位化,可得标准正交基()()()011,,...,n x x x -⎫⎬⎭,令矩阵()()()21011242(1)(1)2(1)(1)(1)111 (1)1...,,...,1..................1...n n n n n n n T x x ηηηηηηηηη-------⎛⎫ ⎪⎪⎫⎪==⎬⎪⎭⎪ ⎪⎝⎭则()()())0111',...n TT x x x --==命题 3.3 任意n 阶循环矩阵()A A f K = 在复数域 C 上都可对角化,即1T AT -=11[(0)(),...,()]n A A A diag f f f ηη-推论 n 阶循环矩阵A 可逆的充要条件是()0iA f η≠(i=0,1,...,n-1).例3:求四阶循环矩阵1234412334122341A ⎛⎫⎪⎪= ⎪⎪⎝⎭的特征根,并对角化.解: 令23()1234f x x x x =+++ 得 ()()A A f K =,0100001000011000K ⎛⎫⎪⎪= ⎪ ⎪⎝⎭由于2i nei πη==, 所以A 的特征根分别为:()()0A f η=10 , ()()1A f η=-2-2i, ()()2A f η=-2, ()()3A f η=-2+2i11111111111211i i T i i ⎛⎫ ⎪--- ⎪= ⎪-- ⎪---⎝⎭, 111111*********i i T i i -⎛⎫ ⎪--⎪= ⎪-- ⎪--⎝⎭4.特殊矩阵特殊对角化的研究前面对实对称矩阵循环矩阵的对角化问题作了研究,本部分主要讨论,当矩阵只有两个特征根时的对角化问题,方法简捷. 对于数域F 上的n 阶矩阵A ,若仅有的两个特征根都在F 内,并且可以对角化,不通过解线性方程组求特征向量,而用初等变换求出可逆矩阵T,使1T AT -为对角形矩阵.定理4.1 设数域F 上的n 阶矩阵A 可以对角化,其特征根为1λ,2λ,如果()10n s n n s B I A p I λ⨯⨯-⎛⎫-⎛⎫−−−−→ ⎪ ⎪ ⎪*⎝⎭⎝⎭初等变换P,B 为列满秩矩阵,那么(i) A 的属于1λ 的线性无关的特征向量为P 的n s -个列向量;A 的属于2λ的线性无关的特征向量为B 的s 个列向量.(ii) 令T = ( P ,B) ,则T 可逆,且有11122......T AT λλλλ-⎛⎫ ⎪⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭其中1λ 有n s -个, 2λ有s 个.证 因为初等矩阵不改变矩阵的秩,且B 为列满秩,则()12s B I A λλ==-=秩秩的重数. (i )根据矩阵的初等变换和分块矩阵的运算性质,可得()())()(1,0n n s I A P B λ⨯--*=,从而()10I A P λ-= 因P 为列满秩矩阵,则P 的n s -个列向量为齐次线性方程组()10I A X λ-= 的基础解系,亦即P 的n s - 个列向量为A 的属于1λ的线性无关的特征向量. 又A 可以对角化,且2λ的重数为s ,则有可逆矩阵Q,使得11122......A Q Q λλλλ-⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭, 令1122......D λλλλ⎛⎫⎪⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭,则有()()()()111212I A I A I Q DQ I Q DQ λλλλ----=--=()()1112QI D QQ I D Q λλ----=()()112Q I D I D Q λλ--- = 10Q OQ -=由于B 的列向量为1I A λ- 的列空间的基,则B 的s 个列向量为齐次线性方程组()10I A X λ-=的基础解系, B 的s 个列向量为A 的属于2λ的线性无关的特征向量.(ii) 因矩阵A 的属于不同特征根的特征向量线性无关,且特征向量的个数之和等于A 的阶数n ,于是, 令 )(,T P B = 即有1T AT D -=例4:令矩阵001010100A ⎛⎫ ⎪= ⎪ ⎪⎝⎭,求可逆矩阵T,使得1T AT -为对角形式.解: 方法一,先求A 的特征根()0101010A f λλλλ-⎛⎫ ⎪=- ⎪ ⎪-⎝⎭= ()()211λλ-+则1λ = 1 (二重) , 2λ = - 1. 可见,此例为定理所述的情况.对矩阵1I A I λ-⎛⎫⎪⎝⎭作初等列变换,即11011000000001011000100101010010001001I A B I P λ-⎛⎫⎛⎫⎪ ⎪ ⎪ ⎪⎪ ⎪---⎛⎫⎛⎫=→= ⎪ ⎪ ⎪ ⎪*⎝⎭⎝⎭ ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭所以,由定理4.1 知,A 的属于2λ = - 1 的线性无关的特征向量为()11,0,1Ta =-;A 的属于1λ = 1 的线性无关的特征向量为()20,1,0Ta = , ()31,0,1Ta =令011100011T ⎛⎫ ⎪= ⎪ ⎪-⎝⎭,则有1111T AT -⎛⎫⎪= ⎪ ⎪-⎝⎭. 这与[1 ]的结果一致.方法二 在矩阵()I A λ-中,亦可取21λ=-,这时1011000200201011000100101010010001001I A B I P ---⎛⎫⎛⎫⎪ ⎪-- ⎪ ⎪⎪ ⎪-----⎛⎫⎛⎫=→= ⎪ ⎪ ⎪ ⎪-*⎝⎭⎝⎭ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭则A 的属于1λ=1 的线性无关的特征向量为()11,0,1Ta =-- , ()20,2,0Ta =- ;A 的属于2λ=- 1 的线性无关的特征向量为()21,0,1Ta =-令101020101T --⎛⎫ ⎪=- ⎪ ⎪-⎝⎭,则有1111T AT -⎛⎫ ⎪= ⎪ ⎪-⎝⎭.5.常规矩阵对角化方法的新探众所周知,对数域P 上一个n 阶矩阵A 是否存在一个可逆矩阵T ,使得1T AT -为对角形矩阵,当这种矩阵存在时,如何去寻求它.一般有关教材中都是先计算一个行列式,求出A 的特征值,再利用线性方程组和特征向量的有关理论及求法解决此问题的.在这里利用矩阵的初等变换解决此问题的,它比教材中的常规方法简单一些,因为不必解若干的齐次线性方程组,有时也不必计算行列式.5.1理论依据为说话方便,我们规定如果数域P 上,对n 阶矩阵存在一个可逆矩T ,使得1T AT -为对角形矩阵, 则称矩阵在数域P 上可对角化.当可对角化时, 我们说将A 对角化,即指求矩阵T ,使1T AT -为对角形矩阵.若矩阵n 在数域P 上可对角化, 则有P 上可逆矩阵T ,使得1T AT B-=为对角形矩阵.于是B 的主对角线上的元素,即为A 的全体特征值, 并且可表示:12,...S T Q Q Q = 其中i Q 为初等矩阵,i=1,2,...,s,于是,1111112......SS S B QQ Q AQ Q Q ----=,又1i Q -也是初等矩阵, 由初等矩阵与矩阵的初等变换的关系, 即知11Q AQ - , 相当于对A 施行了一次初等行变换与一次初等列变换.这里, 我们称此种初等变换为对A 施行了一次相似变换.显见, 可对A 施行一系列的相似变换化为B .又由, 12...S T EQ Q Q =(E 此处表单位矩阵)可如下进行初等变换, 则可将A 化为对角形矩阵B , 且可求得T :A AB E T ⎛⎫⎛⎫−−−−−−−→ ⎪ ⎪⎝⎭⎝⎭对施行一系列相似变换,对E 只施行其中的初等列变换. 当A 不可对角化时, 也可经相似变换化简A 后, 求得其特征值, 判定它可否对角化. 类似地, 可由111111...S S TQ Q Q E -----=,做如下初等变换则可将A 化为对角形矩阵B,且可求得T 或由B 求A 的特征值, 判定可否对角化:()()A AE B T −−−−−−−→对施行一系列相似变换,对E 只施行其中的初等行变换.并且在施行相似变换时, 不必施行一次行变换后接着施行一次列变换这样进行, 可施行若干次行或列变换后再施行若干次相应的列或行变换, 只要保持变换后, 最后所得矩阵与A相似即可.5.2 应用举例为叙述简便,这里用i r 表示i 第行,i c 表示第i 列,i j r kr +表示用数k 乘第j 行后再加到第i 行上,i j c kc +表示用数k 乘第j 列后再加到第i 列上.例5 求如下矩阵的特征值, 并判定它们可否对角化,若可则将其对角化:(1)511602311A -⎛⎫ ⎪= ⎪ ⎪-⎝⎭, (2)1111111111111111B ⎛⎫⎪-- ⎪= ⎪-- ⎪--⎝⎭. 解:(1)由31511`602202r r A +-⎛⎫ ⎪−−−→ ⎪ ⎪⎝⎭ 13411402002c c C --⎛⎫⎪−−−→= ⎪ ⎪⎝⎭,知A 与C 相似. 易得,C 的特征值为2,2,2,且2E-C 的秩为2,所以C 不能对角化,从而知A 的特征值为2,2,2且A 不可以对角化.(2)由1,2,3,41111111111112200111120201111200210001000010001000010001000010000i r r i +=⎛⎫⎛⎫⎪ ⎪-- ⎪ ⎪ ⎪ ⎪-- ⎪ ⎪--⎪ ⎪−−−−→ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭1,2,3,4i c c i -=−−−−−→ 1111,2,3,4,2,3,4441112111222202000200002000200002000210001000110011001010101010011001i i r r i c c i -=+=⎛⎫--⎛⎫ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪−−−−−→−−−−−→ ⎪ ⎪ ⎪ ⎪- ⎪ ⎪- ⎪- ⎪- ⎪ ⎪ ⎪- ⎪⎝⎭-⎝⎭20000200002000021111444311144413114441131444-⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪--- ⎪ ⎪ ⎪--- ⎪ ⎪ ⎪--- ⎪⎝⎭, 知B 可以对角化,B 的特征值为-2,2,2,2.令1111444311144413114441131444T ⎛⎫ ⎪⎪ ⎪--- ⎪=⎪ ⎪--- ⎪ ⎪ ---⎪⎝⎭, 则12000020000200002T AT --⎛⎫⎪⎪= ⎪⎪⎝⎭.当不易直接用相似变换化简判定时, 可先求出特征值, 再用相似变换.例6判定1200320000230043A -⎛⎫ ⎪-⎪= ⎪- ⎪-⎝⎭可否对角化,若可,则将其对角化. 解法1(教材中的方法)由120032000023043x x xE A x x ---=-- ()()()2461x x x =--+,知A 的特征值为4,6,-1,-1.解 齐次线性方程组()40E A X -=得一基础解系23100⎛⎫- ⎪ ⎪⎪ ⎪ ⎪ ⎪⎝⎭解 齐次线性方程组()60E A X -=得一基础解系00341⎛⎫ ⎪ ⎪⎪- ⎪ ⎪ ⎪⎝⎭解 齐次线性方程组()0E A X --=得一基础解系1100⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭,0011⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭于是可,A 可对角化,且取201031010*******01T ⎛⎫- ⎪⎪ ⎪=⎪- ⎪ ⎪ ⎪⎝⎭,则140060000100001T AT -⎛⎫⎪⎪= ⎪-⎪-⎝⎭.解法2由12003200002300431000010000100001-⎛⎫ ⎪- ⎪ ⎪- ⎪- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭ 2143,r r r r --−−−−→ 12,3412004400002300661000010000100001c c c c ++-⎛⎫ ⎪- ⎪ ⎪- ⎪- ⎪−−−−→ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭12000400001300061000110000100011--⎛⎫ ⎪ ⎪ ⎪-- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭123423,57r r r r --−−−−−→2100504003001700061000110000100011⎛⎫-- ⎪⎪⎪⎪-- ⎪⎪ ⎪⎪⎪ ⎪⎪⎪ ⎪⎝⎭214323,57c c c c --−−−−−→100004000010000621005310053001740017-⎛⎫ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪ ⎪⎝⎭知,A 可对角化,且取.21005310053001740017T ⎛⎫- ⎪ ⎪⎪ ⎪= ⎪ ⎪- ⎪ ⎪ ⎪⎝⎭,11000040000100006T AT --⎛⎫⎪⎪= ⎪- ⎪⎝⎭两法比较, 法2比法1简便, 因不必计算行列式和解几个线性方程组.上述内容为本人对各类基本常见的矩阵类型的对角化计算方法,计算技巧的一些探讨,比较传统的计算方法、计算技巧,有一些优越性.计算简便,步骤简单具体,有较强的实用性.参考文献:[1] 张禾瑞 赫炳新 高等代数[M] 第四版 北京 :高等教育出版社 1998.166-410[3] 毛纲源 线性代数[M] 解题方法与技巧归纳 第二版 华中科技大学出版社 1997,7.213-241. [4] 丘维声 抽象代数[M] 北京 :高等教育出版社 2003.160-190.[5] 王萼芳 石生明 高等代数[M] 北京 :高等教育出版社 1987.176-254. [6] 王萼芳 高等代数教程[M] 北京清华大学 1996.91-184.[7] 张爱萍 循环矩阵的性质及其对角化[J] 广西师范自然科学报,2000,12.No.8.168-170. [8] 高吉全 矩阵特征根与特征向量的同步求解方法探讨[J] 数学通报,1991.12.No.7.23-26. [9] 郭亚梅.最小多项式与矩阵的对角化[J]河南机电高等专科学校学报.2006.No.4.106-108. [10]张正成 可对角化矩阵的应用[J] 科技资讯.2007.No.24.252-253.[11]张学元 线性代数能力试题解题[M] 武汉:华中理工大学出版社, 2000.34-37 [12]向人晶 矩阵可对角化的简单判定[J] 数学通报,2003,3.No.12.13-15.[13]靳廷昌有两个特征根矩阵对角化[J] 数学通报,1997,11.No.23.53-57.[14]李世余代数学的发展和展望[J] 广西大学学报.1985.No.1.146-148.[15]周立仁矩阵同时对角化的条件讨论[J] 湖南理工学院学报.2007.Vol.20.No.1.8-10.致谢本论文是在指导师陈巧云老师细心指导下完成.陈老师认真、负责、真诚的做人态度和作为教师对学生不倦教诲的精神,令我很受触动.同时,在论文的选题、修改、定稿都凝聚了陈老师的大量心血.陈老师尽心的指导与严格的监督,促使我最终完成了论文.值此论文完成之际,我谨向陈老师致以深深的敬意和感谢!On the martix diagonatization and application College of science Mathematics 082 Miao Rendong Director:Chen QiaoyunAbstract:This paper initially studied about matrit diagonatization concluding and summarizing about the necessary condition of matrix diagonalization,Through caclulation and research on read synmetrices matrices,cycle matrix,and special matrix diagonalizational ways it proride simple and fast ways of solution on the question of matrix diagonalization in the characteristic root,charateristic rector,and reversible matrix.Key words:diagonal matrix; matrix diagonalizationv; real symmetric matrix;eigenvalue; eigenvectors。

相关文档
最新文档