矩阵对角化及应用论文
矩阵的可对角化及其应用
附件:分类号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.设A∈P椭,P是数域,若存在可逆矩阵X∈P棋”,使X.1AX为对角矩阵,则称矩阵A 可对角化. 关于矩阵对角化,有以下几个结论: 引理1[2] I'I阶方阵A可对角化的充要条件是A有行个线性无关的特征向量. 引理2r23 若行阶方阵A有行个不同的特征值,则A可对角化. 引理3E3]扎阶方阵A可对角化的充要条件是对于A的每个志重特征根.:l。,有
3
1
取X=(Xl X2 X。),则X_1=
3
1— 3
—1)T X3=(1 1 1)r.
1
1
3
3
一号o o
2 3
一丢o 1 .A=X
3
o
X一.
1
1
3
3
O
O
1
口P
一丢o o
口i∞
口∥ =X
0一百1
口≯
厶Hale Waihona Puke OOX一1 口lo)
0
口;o)
1
·+2·(一丢)‘ ·一(一号)‘ ·一(一号)‘ 口i0)
1 3
·一(一丢)‘1-{-2·--号)‘·一--专)‘
口5卜1)),口≯=丢(口{卜1)+口;H)),‘
证明li m口中:竺i::±冀:±.坚(f:1,2,3).
证明
口, 口∥ 口5D
将锑椎关系组改写为下列矩阵形式:
0..1—..—1—
2
2
O
1—2
1—2
1—2
0 1—2
口{卜1) 口扩" 口5H’
n
1
1
u 虿虿
口(o)
高等代数论文
安徽师范大学数学与应用数学系“高等代数选讲”课程论文题目:对角化的讨论及应用姓名:***学号:************安徽师范大学数学与应用数学系数学与应用数学专业2012级2013年8月25日对角化的讨论及应用摘要:本文主要讨论了线性变换的对角化以及实以对称矩阵的对角化的问题,线性变换的对,角化实质上也是矩阵的对角化,分析对角化问题,讨论矩阵是否可与对角矩阵相似,若相似,则有相同的特征值,即可用一定的初等变换将之化为对角阵,以及对角阵在解题材上下班具有比较简便的求法,化一个矩阵为对角阵,不但可以使矩阵运算简化,而且在理论上和应用上都具有十分重要的意义 .关键词: 对角化实对称矩阵特征值相似标准形式.( 一 )线性变换的对角化.1212111,dim(),(),(),()[],()()()(),, ,.:()()()()(),()()t t j j t F n N L g x h x F x g h h g x x x x x l j x x x λλλννσνσσσσλλλσνννλλλλ+-+=∈∈∈=+++=----≠=-j j j,l 正文一对角化的条件:设是数域上的线性空间又设则多项式的运算满足乘法交换律知引理1,设是的两两不同的特征值则和是直和证明g 当时g 令g 1212112,0,1,2,.()()()()00()(0)()()()()()()()()()()()(),()0,0,,(j t lj lid l t t j j j j j j t j t c λαααλανσσσλασσααασασασασαλαλαμμμσν++=∈==-===++=++==≠=j j,l j j j j j j j j 这也是一个多项式,设其中由g g 有g g g g g 因为g g 而g 所以设是的所有的两两不同的特征值.记121212)()(),1,2,,(),()dim(),()(){,,},();(),()(,,)tj j j t j j t III III j t III c n III n III III III diag n μμμμννννννηηηνσηλησλλλσσ=⊕⊕====j 若是的一个基,则将合并得到的向量组线性无关并且是的一组基:引理2:如果而含有个向量.记则是的一组基记则在下的矩阵是推论:如果有个两两不同的特征值,则可对12121()():{,,},,,,,(),dim(()),,dim(),(),t n III VI c m c n F n N L σσσααααααννσννσνμ+∈=≥=∈∈角化.证明:注意到特征子空间的维数是正整数,则此时每一个特征子空间的维数只能是1,故可对角化.引理3:若有m 个向量,m<n,则不可对角化.证明:(反证)如果有基下的矩阵是对角矩阵因此所以与已知矛盾故不可对角化.定理1:设是数域上的线性空间21,,dim()j t tj nμμμσσλ==∑是的所有的两两不同的特征值,的可对角化的充分必要条件是:1212(),,,dim():(1),,(2)(1),n t j j j js j j j j L F m m m m n m Fψμμμψμωωψμμω⨯∈=++==∈A A n A 二 对角化的计算方法:现在考虑F 上n n 矩阵A 的对角化的计算问题,注意到A 否相似于对角矩阵,也就是是否可对角化.设是(也是A 的)的所有的两两不同的特征值,齐次线性方程组(E -A)X=0的解空间记为我们有可对角化的充分必要条件是的重数如果成立将的基(即(121212():{,,},,1,2,,(),,,,j n i i i t n III i n III μξξξξλξψλλλξξξ==n A -1E -A)X=0的一个基础解系),j=1,2,t合并起来得到向量组于是A 而在下的矩阵是对角矩阵B=diag()(3)若取C=(),则C 是可逆的,并且C AC=B三 相似标准形现在假定A 可对角化,我们来研究与A 相似的对角矩阵是否在某种意义下"惟一"?也就是说,如果G,H 都是12121112,,,,()()()()()(){,}t t j t i n n x x x x x x λλλμμμλλλμμμψεεε---=---A 对角矩阵,并且G H,G 与H 有什么进一步的关系?引理4:如果G=diag(),H=diag()都是对角矩阵,并且G H,则G 与H 只相差主对角线上排列的不同:证明:先证必要性:如果G H,则G 与H 的特征多功多项式相等因此,则结论成立.再证充分性,G 是在基12,()(),,,,t VI F P C σλλλ⨯⨯∈∈=-1n n n n 下的矩阵可以适当排列这个基得到另一个基(VI)使得在下的矩阵恰是H:定义1,设A Mat 如果存在F 上可逆矩阵C,使C AC 是对角矩阵G=diag(),则称G 为A 的(相似)标准形.定理2:如果矩阵A 可对角化,则它的相似标准形在引理4下意义下惟一.四 可对角化的等价情况设P 为数域A P 当时则下列条件等价: 1,A 的每个若尔当块皆为1级的, 2,A 的最小多项式无重根, 3,A 的最后一个不变因子无重根, 4,A 的初等因子是一次的, 5,A 的特征多项式无重根.112*,,,0n n λλλλλ<>⎛⎫ ⎪⎪ ⎪⎝⎭-1-1二实对称矩阵的对角化引理1'任意n 阶复矩阵A 必相似于上三角形矩阵,即存在可逆矩阵P, 使得P AP=其中为A 的全部特征根引理2'实对称矩阵的特征根为实数.引理3'设A 与B 为n 阶实矩阵,则A 与B 在实数域上相似的充要条件是 A 与B 在复数域上相似,即有实可逆矩阵P,使得P AP=B 的充要条件 1,μμμμ=-≠≠-1-12 是有可逆复矩阵Q,使得Q AQ=B. 证明:必要性显然.下证充分性,设有可逆复矩阵Q,使得Q AQ=B,且令Q=C+iD,其中C,D 是 实矩阵,i 而AQ=QB,于是AC=QB,AD=DB,因Q 可逆,故|Q|=|C+iD|0 即|C+iD|不是零多项式,则有实数,使得|C+D|0,即实矩阵C+D 可逆 令P=C+D,则有AP=AC 1*",0n μμλλ⎛⎫ ⎪⇒ ⎪ ⎪⎝⎭-1-1T -1+AD=CB+BD=PB,于是Q AP=B,即A 与B 在R 上相似 定理1'(1) n 阶实矩阵A 的特征根都是实数的充要条件是存在正交矩阵Q,使得 Q AQ= Q AQ 为上三角形矩阵.(2)当A 的特征根都为实数且A 为正交矩阵时,则A 为对称矩阵. 证明:(1) "由引理1'知,存在可逆矩阵P,使得P AQ= 由1*,0"n λλ⎛⎫ ⎪= ⎪ ⎪⎝⎭⇐-1-1-1-1-1-1引理3'知,当在复数域上相似时,必有实域上相似,因而可令P 为实可逆矩阵,令P=QT,其中,Q 为正交矩阵,T 为上三角形矩阵,且主对角线上的元素都为正数 ,于是Q AQ=T(P AP)T 因T 为上三角形矩阵,则T 也为上三角形矩阵,可见Q AQ 为上三角形矩阵之积,所以,Q AQ 为上三角形矩阵 "设对于实矩阵A 11212*,,,0,,,n n n λλλλλλλλ⎛⎫⎪=⎪ ⎪⎝⎭=-1-1T T T -1-1T ,有正交矩阵Q,使得Q AQ 则为A 的全部特征根,且都是实数.(2)由(1)知存在正交矩阵Q,使得 Q AQ= Q AQ=B 为上三角形矩阵,因A 是正交矩阵,Q 是正交矩阵 Q 是正交矩阵,则B 作为正交矩阵之积自然也是正交矩阵,于是B B 而B 为上三角形矩阵, B 为下三角形矩阵,故)T ====T T T T T B 必为对角形矩阵,易得A (QBQ QB Q QBQ A-1T T T T T T T 定理2' n 阶实对称矩阵A 必正交相似于对角形矩阵,即有正交矩阵Q,使得 Q AQ= Q AQ=B 为对角矩阵.证明:因A 为实对称矩阵,则由引理2'知A 的特征根都为实数,又由定理1' 的(1)知,有正交矩阵Q,使得Q AQ=B 为上三角形矩阵,而A 是对称的,所以 B =(Q AQ)=Q AQ=B,但B 为下三角形矩阵,故B 必有对角形矩阵定理3'若n 阶实矩阵A 既正定又10,0110,,01I ⎛⎫⎪⎪ ⎪⎝⎭⎛⎫ ⎪= ⎪ ⎪⎝⎭T T正交,那么,A=I(单价阵)证明:因A 正定,则A 的特征根都是正实数,又A 是正交阵,则A 的特征根只能 均为1,从而,有正交矩阵Q,使得,Q AQ=所以A=Q Q(三)练习12311:,155,1,,λλλλλλλξ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦-+-===---1233T 1122例设A=212A 是否可对角化?如果A 可对角化,求可逆矩阵C,221 使得C AC 是对角矩阵.解:A 的特征多项式是|E A|=()()故A 的全部特征值是解齐次线性方程组(E A)X=0, 得到它的一个基础解系:{=(1 1 1)}, 这也是A 属于5的特征子空间的一个基,解齐次2,λξξξξξ-3T T 23123-1线性方程组(E A)X=0,得到它的两个基础解系: {=(-1 0 1),=(0 1 -1)},这是A 属于-1的特征子空间的一个基,令C=(,),则C 是可逆的, C AC=diag(5,-1,-1)⨯⨯T T T T T T T T T 例 2 设A,B 是两个n n 实对称矩阵,且B 是正定矩阵,证明存在一个n n 实 可逆矩阵T, 使得T AT 与T BT 同时为对角形. 证明: B 是正定矩阵,故存在可逆矩阵Q,使得Q BQ=EA 是实对称矩阵,故Q AQ 也是实对称矩阵,故存在正交矩阵P,使得P Q AQP=(QP)AQP 为对角矩阵 故令T=QP,则T AT 为对角矩阵. 则T BT=(QP T T T T )BQP=P EP=P P=E 则T BT 也是对角矩阵.123123,:2(6),2,6,2(1,1,0),(1,1,0),6(1,2,3)112E λλλλλλλλλ⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦--======-===-⎡⎤⎢⎥-⎢⎥⎢⎣⎦k2T 1T T 231-11 例 3 已知 A=24-2求A -3-35解可求得det(A-)=-()所以A 的特征值为 对应于有两个线性无关的特征向量P P 对应于的特征向量为P 故A 可对角化,-1 则P=10013111,(2,2,6),()5*222123*222*43*23*23*2k k k k k k k k k k diag νν+++=⎥==⎡⎤---⎢⎥--⎢⎥⎢⎥--⎣⎦-1k-1k k k -1k k k k k k k k k P AP= 所以A P P Pdiag(2,2,6)P 66+6 =+2*6+2*6666+3*6121412,lim 21007111:,,,247111(,,),247111,,,lim 0247n n n n n diag →∞→∞⎡⎤-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦==n -1n-1n 例4 已知 A=求A 的值解A 有三个互异的特征值所以存在可逆阵P使得P AP= 而A Pdiag()P 故A1231231,1,(0,1,1),1,,,0,0λλλλλλ=-=====<>=====T 123T 112312311223例 5 设三阶实对称矩阵A 的特征值为 对应于的特征向量为P 求A解:因为A 为三阶实对称矩阵,故必可对角化,又因是A 的二重 特征值,故A 的与特征值1对应的线性无关的特征向量有两个,设为P P 且都与P 正交,设所求特征向量为X=(x ,x ,x )则P X 即x +x x x 由x xx 123123,(1,0,0),(0,1,1),(0,1,1),(1,0,0),1)|||110,0,100110110εεεεεενν⎧⎪==-⎨⎪-⎩======-⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎥==⎢⎥⎥⎢⎥⎥⎢-⎢⎣=T T232T T T123123-1T-1得P P xP P P 规范化得|P |P |P 010作正交矩阵P=(,)=则P P 有A=P PP 010010011001010000100101000⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-⎡⎤⎡⎤⎢⎥⎥⎢⎥⎢⎥=-⎢⎥⎥⎢⎥⎢⎥⎢⎥⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎢-⎢⎣T P 010所以A=参考文献:1 屠伯埙 : 线性代数—方法导引 上海:复旦大学出版社,1968.2 李师正: 高等代数解题方法与技巧 高等教育出版社 20043 旋武杰: 高等代数 高等教育出版社 4陈重穆 等: 高等代数 北京:高等教育出版社 19905张禾瑞,郝炳新高等代数北京:高等教育出版社19836 北京大学数学系几何与代数教研室代数小组编北京:高等教育出版社7 扬子胥:高等代数习题解(修订版)山东科学技术出版社2001。
矩阵的正交对角化及其应用分析
矩阵的正交对角化及其应用分析矩阵是线性代数中重要的概念之一,广泛应用于各个领域。
在矩阵理论中,正交对角化是一种重要的矩阵分解方法,它在解决线性方程组、求解特征值和特征向量等问题中具有重要的应用。
正交对角化是指将一个矩阵通过相似变换转化为对角矩阵,并且这个相似变换矩阵是正交矩阵的过程。
具体而言,对于一个n阶方阵A,如果存在一个正交矩阵P,使得P^(-1)AP=D,其中D为对角矩阵,那么我们称A可被正交对角化。
为了实现正交对角化,我们需要先找到矩阵A的特征值和特征向量。
设λ为A的一个特征值,x为对应的特征向量,那么有Ax=λx。
通过求解特征方程det(A-λI)=0,我们可以得到A的所有特征值。
对于每个特征值,我们可以通过求解(A-λI)x=0得到对应的特征向量。
接下来,我们需要将特征向量组成一个矩阵P,使得P的每一列都是A的一个特征向量。
然后,我们可以计算P^(-1)AP,得到一个对角矩阵D。
这个过程中,我们需要确保特征向量的线性无关性,以及特征向量的单位化。
正交对角化的应用非常广泛。
首先,它可以用于解决线性方程组。
对于一个矩阵A,如果它可被正交对角化,那么我们可以通过正交对角化将线性方程组转化为对角方程组,从而更加方便地求解。
其次,正交对角化可以用于求解特征值和特征向量。
通过正交对角化,我们可以将一个矩阵转化为对角矩阵,对角矩阵的对角线元素就是原矩阵的特征值。
同时,对角矩阵的每一列就是原矩阵的特征向量。
此外,正交对角化还可以用于矩阵的相似变换。
相似矩阵具有相同的特征值,因此通过正交对角化,我们可以将一个矩阵转化为对角矩阵,从而更容易进行矩阵的运算和分析。
需要注意的是,并非所有的矩阵都可以被正交对角化。
一个矩阵可被正交对角化的充要条件是它是对称矩阵。
对称矩阵具有一些特殊的性质,例如它的特征值都是实数,特征向量之间正交等。
总结来说,矩阵的正交对角化是一种重要的矩阵分解方法,它在解决线性方程组、求解特征值和特征向量等问题中具有广泛的应用。
矩阵的对角化及其应用
矩阵的对角化及其应用13届分类号:单位代码:10452临沂大学理学院毕业论文(设计)矩阵的对角化及其应用2013年3月20日临沂大学理学院2013届本科毕业论文(设计)摘要矩阵可对角化问题是矩阵理论中的一个重要问题,可对角化矩阵作为一类特殊的矩阵,在理论上和应用上有着十分重要的意义.本文对可对角化矩阵做出了较全面的概括和分析,并利用高等代数和线性代数的有关理论总结出了矩阵可对角化的若干条件,同时也讨论了化矩阵为对角形的方法,最后总结出可对角化矩阵在求方阵的高次幂)利用特征值求行列式的值)由特征值和特征向量反求矩阵)判断矩阵是否相似)向量空间)线性变换等方面的应用.关键词:对角化;特征值;特征向量;相似;线性变换临沂大学理学院2013届本科毕业论文(设计)ABSTRACTMatrix 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 ofmatrix several conditions diagonolization, also discussed the matrix of the diagonal shape of 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临沂大学理学院2013届本科毕业论文(设计)目录1 引言 (1)2矩阵对角化 (1)2.1可对角化的几个条件 (1)2.2可对角化的矩阵的性质 (3)2.3 矩阵的对角化 (5)2.3.1 用矩阵初等变换将矩阵对角化的方法 (5)2.3.2 利用内积构造齐次线性方程组的方法 (7)3 矩阵对角化的应用 (10)3.1 求具有线性递推关系( 组) 的数列的通项式与极限 (10)3.2 求解行列式的值 (14)3.3对角矩阵的其他方面的应用.................................... 15 4 结论 .......................................................... 19 参考文献 ..................................................... 19 致谢 (21)临沂大学理学院2013届本科毕业论文(设计)1 引言对角化矩阵在求解一类具有递推关系式的数列的通项与极限及一类三对角线行列式、求方阵的高次幂、利用特征值求行列式的值、由特征值和特征向量反求矩阵、判断矩阵是否相似、在向量空间、线性变换等方面的应用.对角矩阵贯穿于高等代数之中,有着十分重要的作用.定义1.1 对角矩阵是一个主对角线之外的元素皆为 0 的矩阵.对角线上的a元素可以为 0 或其他值.因此行列的矩阵= 若符合以下的性质: nnAa,,ij,ij,10,, ij,1,2,,…,nij,=0,,.形如. ,,,,01,,V定义1.2 矩阵可对角化:设是维线性空间的一个线性变换,如果存在n,V的一组基,使在这组基下的矩阵为对角矩阵,则称线性变换可对角化. ,, 定义1.3 矩阵是数域上的一个维方阵,如果存在数域上的级可逆APPnn,1TAT矩阵,使为对角矩阵,则称矩阵可对角化. AT2矩阵对角化通俗地说就是经过矩阵的一系列行、列变换(初等变换)后,能得到一个只有主对角线上元素不全为零,而其他位置全为零的另一个矩阵(这个矩阵称为对角阵),这个过程就叫做矩阵的对角化,并不是所有矩阵都能对角化. 2.1可对角化的几个条件矩阵可对角化在求矩阵的高次幂中有重要应用, 矩阵的对角化有多种判别方法.本节对矩阵对角化作一点讨论,nn,22,PABB, 引理2.1 设,,且=,,.则存在可逆矩阵,ABAABBA,P使,可同时对角化. ABnn,,Pdiag,,,,,…,引理2.2 如果=有个互不相同的对角元,对某Pn,,12nnn,,P个,则当切仅当本身是对角阵. BPBBP,BE0,,r2AA,由于任意一个幂等矩阵A必相似于对角矩阵.而且每个与对角,,,,00,,n矩阵都可以进行谱分解,即=,A,其中是的特征值,为幂等阵.那么AA,A,iiiii,1任意有限个幂等阵的线性组合是否对角化,有如下结论:1临沂大学理学院2013届本科毕业论文(设计)定理1.1 以=,为个数,Akkk,,,,,,…kkk,,…,n1122nn12nij,为个幂等阵,且两两可换,即=,,则可对,,,,,…,,,,,An,,12nijji 角化.证明为个幂等阵,且两两可换.由引理1可知,存在可逆阵,,,,,…,n12n ,1,1,使可同时对角化.即,…,,,,,,,…,,,,PPP,,,PP12n1nnnn1111 ,1,1,,,…,是对角阵PkP,PkP,.==++…Akkk,,,,,,…,,,,,,1n11221122nn,1,1PkP,PkkP,,,,,…+k++.由知,,,…,是对角阵,,,,nnnn11221n 也为对角阵,故可对角化. Akk,,,,,…+k1122nn如果矩阵只有两个不同的特征值,可有如下结论:nn,,P定理2.2 设,,为其两个不同的特征值,则可对角化存在AA,,,12 ,,,,幂等矩阵,使得=+,其中为幂等阵. ,AE,,,,211,E,,11-1证明必要性:若可对角化,存在可逆矩阵,使=相似APPAP,,11E,,,22,1PP,于对角阵,则= A,,0,,,1 = ,PEP,,,,,1,,,E,,,,21,,,,0,,,1,1,, =+, PPPEP,,,,,2111E2,,0,,,1,, =+, PP,E,,,,2111E2,,000,,,,,,,1,,112,PPPP且相似于== ,PP,,,,,,EEE222,,,,,,,,,,故为幂等阵,即存在幂等阵使得=+. ,,AE,,,211,,,,充分性:若存在使=+.因为为幂等阵,故存在可逆阵,使,AE,T,,,211 00,,,,,1,1,,,得=,则=+TT ,TTAE,,,,,,,211EE22,,,,2临沂大学理学院2013届本科毕业论文(设计),,0,,,1= ,TET,,,,,,,,E,,,,21,,,,,E,,11,1= TT,,E,22,,,E,,11,1,1TAT故= TT,,E,22,,即可对角化. AAB,如果满足条件的情况,有如下结论: ABBA=nn,nn,,P,P定理2.3 假设个互不相同的特征值,对某个个ABn有,则有AB,当且仅当同时对角化. ABBA=,1TPAP,证明必要性.由有个互不相同的特征值,则可对角化.设,AAndiag,,,,…,其中=.则T,,12,n,,11,,11,1,1,1,1PAPPBPPABPPBAPPBPPAPTPBPPBPT=====.即与T,,,,,1,1PBPPBP可交换,由引理2知是对角阵,从而是可对角化矩阵. B ,1AB,充分性.可同时对角化,故存在可逆阵,使得,PAPP,,1,1其中,为对角阵,,,1. BPP,,,,22,,11,1,1,,11=====. ABBAPPPP,,PP,,PP,,PPPP,,21121221对定理,我们可得到矩阵只有两个不同的特征值时可对角化的判别方法: A22,,=,,,,AE,,,/若,则可对角化,否则不可对角化.其中. AA,,,,122.2可对角化的矩阵的性质是数域上的一个可对角化的阶矩阵,是定理2.2.1 设APA,,,,,…,n12t 阶矩阵,使AA,,…,An的互不相同的特征根,则存在12t1+AAAA,,,,,,…; ,1122tt2+=E,EAAA,,…为单位矩阵; ,12t23AA,; ,ii,140,AAij,,,0为零矩阵,其中. ATBT,,ijii1证明由上一个阶可逆矩阵,使得 APTn,可对角化,则存在3临沂大学理学院2013届本科毕业论文(设计),0,,1,,,2,1,,TATB,, ,,…,,0,t,,其中的重数为,由于 s,ii10,,,,,,,,……,,,,,,,,10 B,,,,…+,,,,1t01,,,,,,,,……,,,,,,,,01,,,,记,所以 ,,,BBB,,…+1122tt,,11ATBTTBBBT,,,,,,,…+ ,,tt1122,,,111= TBTTBTTBT,,,,,…+t t1122,,11,,TBTTBT,…+= ,,,,tt11,1记,其中 ,,,AAA,,…+ATBT,1122tti故. AAAA,,,,,,…+1122tt 2由每个为对角形幂等阵,则, BBBE,,,…+B,12ti,1,1,,,111TET=ETBBBT,,…+===AAA,,…+TBTTBTTBT,,…+,,,t12t12t12故 AAA,,…+=E12t,,11,12,,113TBTTBT由,则== ATBT,ATBTTBT,,,,,iiiiiii,121,,1==,TBTTBBTTBTiiii=Ai2故. AA,ii,,11,1,,11TBTTBTTBBT4ij,TBTTBT当时,====0;0为零矩阵 AA,,,,,ijijijij故 AAij,,0,ij4临沂大学理学院2013届本科毕业论文(设计)15115,,,,1,2,3例2.2.1 设是数域上的矩阵.是矩阵的特征根,A,,20158PA,,,,876,,, 100231,,,,,,,,,1则存在可逆矩阵,使得=,其中T,342TAT,020B,,,,,,,,112003,,,,,,652,,,,,1, T,,431,,,,111,,,100,,,,,,,,,,,,由于,记 BBB,,23B,,,02130123,,,,,,,,,,,,001,,,,,,,,11ATBTTBBBT,,,,23所以,,123,,,111TBTTBTTBT,,23= ,,,,123,1=,其中AAA,,23ATBT,123ii,,,,,121041293111,,,,,,,,,,,,,且满足: AAA,,,,,,,18156,16124,222123,,,,,,则 ,,,,,,,,,,652431222,,,,,,123AAAA,,,; ,1232AAAE,,,; ,12323AA,i,1,2,3,; ,,,ii40,AAij,,,0为零矩阵. ,,,ij通过一个具体的可对角化矩阵,鲜明地反映了上述性质是成立的.2.3 矩阵的对角化2.3.1 用矩阵初等变换将矩阵对角化的方法VV数域P上维线性空间的一个线性变换判定其是否在中能找到一组基n使它在此基下的矩阵为对角形矩阵; 当这种基存在时, 如何去寻求它是线性代数学上一个十分重要的问题,利用矩阵的初等变换法解决此问题.,1TAT若矩阵在数域上可对角化,则有上可逆矩阵使=为对角阵.APPBT5临沂大学理学院2013届本科毕业论文(设计) 于是的主对角线上的元素即为的全体特征值,并且可表示为,BATQQQ,…12s,1,,,111i,1,2,…,s其中为初等矩阵,,于是,,又也BQQQAQQQ,……QQis1112sS,,1是初等矩阵,由初等矩阵与矩阵的初等变换的关系,即知相当于对施AQAQ11行了一次初等行变换和一次初等列变换,将此种变换称为对施行了一次相似变A注:为单位矩阵E换.又由,可进行如下初等变换,则可将化A TEQQQ,…,,12S 为对角矩阵,且可求得: BTAB,,,,对施行一系列相似变换A,对只施行其中的初等列变换. E,,,,,,,,,,,,ET,,,,当不可对角化时,也可经相似变换化简后,求得其特征值,判定它可否对角AA 化.-1,,,111T类似地,可有=,做如下的初等变换则可将化为对角形矩阵AQQQE…s11s,,且可求得或由求的特征值,判定可否对角化: BBAAT对施行一系列相似变换A,1AEBT,,,,,,,,,对对只施行其中的初等行变换. E,,,,并且在施行相似变换时,不必施行一次行变换后接着施行一次列变换这样进行,可施行若干次行(或列)变换后再施行若干次相应的列(或行)变换,只要保持变换后,最后所得矩阵与相似即可. Ajk为叙述简便,这里用表示第行,表示第列,表示用数乘第行iicrkr,riiji jk后加到第行上,表示用数乘第列后加到第列上. iickc,ij注意到初等矩阵的逆矩阵,,11,1,1PijPijPikPijkPijkPijk,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,1A故用左乘A,相当于对施行了变换,右乘A,相当于对A施rkr,Pijk,,,,,ji行了变换. ckc,ji例 2.3.1 求如下矩阵的特征值,并判定它们可否对角化,若可,则将其对角化: 1111,,511,,,,,1111,,,,,,21602,,,,,,,,,1111,,,,,311,,,,1111,,,,;511,411,,,,,,,,,rr,cc,31131CC解由=,知与相似.A,,,,602,,,,402A,,,,,, ,,,,,311002,,,,C2,2,2,2EC,C易知,的特征值为的秩为,所以不可对角化,从而知的特2A6临沂大学理学院2013届本科毕业论文(设计)2,2,2,征值为且不可对角化. A11111111,2111,,,,,,,,,,,,1111,,22000200,,,,,,,,,,,,1111,,20200020,,,,,,1111,,20020002rri,,,2,3,4cci,,,2,3,4,,,,,,i11i,,,,,,,,,,,2由,,,,,,,,100010001000,,,,,,01000100,1100,,,,,,,,,,,,00100010,1010,,,,,,,,,,,,00010001,1001,,,,,,,,2000,,111,,,2,,,,0200222,,,,,,00200200,,,,,,00020020,,1,,,2,3,4rri, ,,,1i,知可对角化,的BB111400021,,,,,,,,1,,cci,,,2,3,4ii4444,,,,,,,,,,1000,,311,,,,,1,,,,,11 00444,,,,131,1010,,,,,,,1,,,,444,1001,,,,113,,,,,1,,,,444 111,,1,,444,2,,,,,,311,,2,1,,,1,,,2,2,2,2。
矩阵的可对角化及应用
矩阵的可对角化及应用矩阵的可对角化是线性代数中一个重要的概念,它在解决线性方程组、求解特征值和特征向量等问题中有着广泛的应用。
本文将从可对角化的定义、判定条件、可逆矩阵的可对角化以及应用等方面进行论述。
首先,我们来定义可对角化矩阵。
一个n阶方阵A称为可对角化的,如果存在一个可逆矩阵P,使得P^-1AP=D,其中D是一个对角矩阵。
那么如何判定一个矩阵是否可对角化呢?根据线性代数基本定理,一个n阶方阵A可对角化的充分必要条件是:A有n个线性无关的特征向量。
具体地说,对于特征值λ及其对应的特征向量v,如果A存在n个线性无关的特征向量v1,v2,...,vn,对应特征值λ1,λ2,...,λn,则A可对角化。
需要注意的是,A不一定有n个不同的特征值,但是至少要有n个线性无关的特征向量。
关于可对角化矩阵,还有一个重要的结论是:若A可对角化,且λ1,λ2,...,λn是A的n个特征值,则A的k次幂矩阵A^k可对角化,且对应的特征值是λ1^k,λ2^k,...,λn^k。
下面我们来讨论可逆矩阵的可对角化。
对于一个可逆矩阵A,它可以写成对角化的形式A=PDP^-1,其中D是对角矩阵,P是可逆矩阵。
根据可逆矩阵的性质,P^-1A=DP^-1,即A和D有相似的特征值和相同的特征多项式。
在实际应用中,矩阵的可对角化具有许多重要的应用。
首先,对于一个可对角化的矩阵A,我们可以很方便地求解线性方程组Ax=b。
我们可以先通过对A进行对角化,得到A=PDP^-1,然后将方程组转化为DP^-1x=P^-1b,令y=P^-1x,则有Dy=P^-1b,由于D是对角矩阵,所以这个方程组的解可以很容易地得到。
最后通过y=P^-1x,我们可以求得方程组Ax=b的解x。
其次,可对角化矩阵在求解特征值和特征向量的问题上也有着重要的应用。
对于一个n阶可对角化矩阵A,我们可以通过对A进行对角化,得到A=PDP^-1。
由于D是对角矩阵,它的对角线元素就是A的特征值。
矩阵对角化论文 最终版
矩阵对角化的判定条件及应用摘要:对角矩阵是矩阵中简单的一种,在高等代数中占有极其重要的位置。
本文归纳总结了矩阵对角化的若干方法,并且分情况讨论了有n个特征根的n阶矩阵的对角化方法。
然后基于定义及判定定理,引出了实对矩阵、反循环矩阵的若干重要性质,为读者对矩阵对角化中求特征值、特征向量、求可逆矩阵、使对角化问题,提供了简便、快捷的求解途径。
最后列举了几种常用矩阵的对角化问题。
关键词:矩阵,对角化,特征根,反循环矩阵Criterions of Matrix Diagonalization and GeneralizationsDiagonal matrix is a simple matrix in algebra occupies an extremely important position. This article summarizes some of the matrix diagonalization method, and the Points of discussion are n characteristic roots of the n-order matrix diagonalization. Then determined based on the definition and theorem, leads to the fact of the matrix, a number of important anti-cyclic nature of the matrix, matrix diagonalization for the readers in the eigenvalue, eigenvector, Inverse Matrix, so that diagonalization, to provide a simple, efficient way to solve. Finally, several commonly cited problem of matrix diagonalization.Key words: matrix, diagonalization, eigenvalues, anti-cyclic matrix引言可对角化矩阵是一类重要矩阵,它与有限维线性空间的线性变换的对角化问题密切相关,因而是矩阵论重点研究的内容之一。
矩阵的对角化及其应用
学院2016届本科毕业论文(设计) 矩阵的对角化及其应用学生:学号:专业:数学与应用数学指导老师:答辩时间:2016.5.22 装订时间:2016.5.25A Graduation Thesis (Project)Submitted to School of Science, Hubei University forNationalitiesIn Partial Fulfillment of the Requiring for BS DegreeIn the Year of 2016 Diagonalization of the Matrix and its ApplicationsStudent Name Student No.:Specialty: Mathematics and Applied Mathematics Supervisor:Date of Thesis Defense:2016.5.22 Date of Bookbinding: 2016.5.25摘要矩阵在大学数学中是一个重要工具,在很多方面应用矩阵能简化描述性语言,而且也更容易理解,比如说线性方程组、二次方程等. 矩阵相似是一个等价关系,利用相似可以把矩阵进行分类,其中与对角矩阵相似的一类矩阵尤为重要,这类矩阵有很好的性质,方便我们解决其它的问题. 本文从矩阵的对角化的诸多充要条件及充分条件着手,探讨数域上任意一个n阶矩阵的对角化问题,给出判定方法,研究判定方法间的相互关系,以及某些特殊矩阵的对角化,还给出如幂等矩阵、对合矩阵、幂幺矩阵对角化的应用.关键词:对角矩阵,实对称矩阵,幂等矩阵,对合矩阵,特征值,特征向量,最小多项式AbstractThe matrix is an important tool in college mathematics, and can simplify the description language based on the application of matrix in many ways. So it is easier to understand in many fields, for example, linear equations, quadratic equations. In many characteristics, the matrix similarity is an very important aspect. We know that the matrix similarity is an equivalence relation by which we can classify matrix, the diagonal matrix is very important. This kind of matrix has good properties, and it is convenient for us to solve other problems, such as the application of similar matrix in linear space. In this paper, we first discuss many necessary and sufficient conditions of diagonalization of matrix and then give some applications of special matrix diagonalization.Key words: diagonal matrix,real symmetric matrix,idempotent matrix,involutory matrix,the eigenvaule,the feature vector,minimalpolynomial目录摘要 (I)Abstract (II)绪言 (1)课题背景 (1)目的和意义 (1)国外概况 (1)预备知识 (2)相关概念 (2)矩阵的对角化 (4)特殊矩阵的对角化 (14)矩阵对角化的应用 (22)总结 (24)致谢 (25)参考文献 (26)独创声明 (28)1 绪言本课题研究与矩阵的对角化相关的问题,从对角化的判定展开论述,结合其它学术期刊的结论加上自己的体会,希望能让读者更好的理解矩阵及其对角化的妙处.1.1 课题背景在由大学数学系几何与代数教研室前代数小组编、王萼芳与石生明修订、高等教育出版的《高等代数》一书中,我们为了方便线性方程组的运算引入了矩阵的概念. 在线性方程组的讨论中我们看到,线性方程组的系数矩阵和增广矩阵反应出线性方程组的一些重要性质,并且解方程组的过程也表现为变换这些矩阵的过程.除线性方程组之外,还有大量的各种各样的问题也提出矩阵的概念,并且这些问题的研究常常反应为有关矩阵的某些方面的研究,甚至于有些性质完全不同的、表面上完全没有联系的问题,归结为矩阵问题以后却是相同的. 在二次型中我们用矩阵研究二次型的性质,引入了矩阵合同、正定、负定、半正定、半负定等概念及其判别方法.在向量空间中用矩阵研究线性变换的性质,引入矩阵相似的概念,这是一种等价关系,利用它我们把矩阵分类,其中与对角矩阵相似的矩阵引起的我们的注意,由此我们对线性变换归类,利用简单的矩阵研究复杂的,方便我们看待问题,进而又引入对角型矩阵、λ矩阵及若尔当标准型.本文主要由矩阵定义和向量空间研究矩阵的对角化,从不同角度揭示矩阵对角化的判定及其性质,还给出特殊矩阵的对角化及其相应的应用.1.2 课题研究的目的和意义课题研究的意义:(1) 研究矩阵对角化的判定定理及应用,为其它学术研究提供便捷的工具;(2) 比较全面的介绍矩阵的对角化,方便读者的整体理解和应用;1.3 国外概况实数域、复数域等数域上的矩阵的对角化研究已经很成熟,涉及特征值、最小多项式、线性变换方面的对角化证明也已完善,四元素体上矩阵的广义对角化也有小有成就,矩阵对角化与群环域的结合方面的研究也有所突破. 实对称矩阵、正交矩阵、分块儿矩阵的对角化已完善,矩阵的应用也渐渐出现在更多的学科和科研当中. 矩阵的同时对角化、同时次对角化,以及对角化与秩的恒等式等方面的研究基本完善.2 预备知识给出本文容所涉及的一些定义,方便对后面定理证明的理解. 定义1 常以n m P ⨯表示数域P 上n m ⨯矩阵的全体,用E 表示单位矩阵.定义2 n 阶方阵A 与B 是相似的,如果我们可以找到一个n 阶非奇异的方阵矩阵T n n P ⨯∈,使得AT T B 1-=或者BT T A 1-=.根据定义我们容易知道相似为矩阵间的一个等价关系:①反身性:AE E A 1-=; ②对称性:若A 相似于B ,则B 相似于A ; ③传递性:如果A 相似于B ,B 相似于C ,那么A 相似于C .定义3 n 阶方阵A 与B 是合同的,如果我们可以找到一个n 阶非奇异方阵T n n P ⨯∈,使得B =T T AT 或者BT T A T =.根据定义我们容易知道合同也为矩阵间的一个等价联系:①反身性:A =AE E T ;②对称性:由AT T B T =即有11)(--=BT T A T ;③传递性:由111AT T A T=和2122T A T A T =有)()(21212T T A T T A T =.定义4 式为⎪⎪⎪⎪⎪⎭⎫⎝⎛⋯⋯⋯m b b b 000000021 的m 阶方阵叫对角矩阵,这里i b 是数(),2,1m i ⋯⋯=.定义5 方阵A n n P ⨯∈,若BT T A 1-=,T 非奇异,B 是对角阵,则称A 可相似对角化.定义6 方阵A n n P ⨯∈,若BT T A T =,T 非奇异,B 是对角阵,则称A 可合同对角化.定义7 矩阵的初等变换:⑴互换矩阵的第i 行(列)于j 行(列); ⑵用非零数c P ∈乘以矩阵第i 行(列);⑶把矩阵第j 行的t 倍加到第i 行.定义 8 由单位矩阵经过一次初等行(列)变换所得的矩阵称为初等矩阵. 共有三种初等矩阵:①单位矩阵经过初等变换⑴得),(j i P 且),(),(1j i P j i P =-;②单位矩阵经过初等变换⑵得))((t i P 且)/1(())((1t i P t i P =-;③单位矩阵经过初等变换⑶得))(,(t j i P 且))(,())(,(1t j i P t j i P -=-.定义9 设方阵n n P B ⨯∈,若E B =2,就称B 为对合矩阵.定义10 设方阵n n P A ⨯∈,若A A m =,就称A 为幂幺矩阵.定义 11 设方阵C n n P ⨯∈,若C C =2,就称C 为幂等矩阵.定义 12 设方阵n n P A ⨯∈,P ∈λ,若存在向量,满足X Al λ=,我们就称λ是A 的特征值,X 是A 属于特征值λ的特征向量.定义13 n n P A ⨯∈,定义)(λA m 为矩阵A 的最小多项式 ,)(λA m 的一个根为A 而且比其他以A 为根的多项式的次数都低,)(λA m 首项系数是1.3 矩阵的对角化本章介绍数域P 上n 阶方阵阵的对角化问题.先给出矩阵对角化几个一般的充要、充分条件及其证明.引理1 如果1μ,…,k μ是矩阵Q 的不同的特征值,而,1i α…,i ir α是属于特征值i λ的线性无关的特征向量,2,1=i …,k ,那么,11α…,11r α,…,1k α,…,kkr α也线性无关.证明:假设++12121111ααt t …++1111r r t α…++11k k t α…k k kr kr t α+=0,P t ij ∈,令+11i i t α…+i i i i k k t α=i η,则i i i Q ηλη=(2,1=i ...k ,), 且 ++21ηη...k η+=0 (1)分别用,,,2Q Q E …1,-k Q 左乘以(1)两端,再由引理4得:i i i m Q ηλη=,(1...2,1-=k m ;t i ,...,1=),由此有⎪⎪⎪⎩⎪⎪⎪⎨⎧=++=++=++=++---.0......................................,0...,0...,0...12121112222121221121k k k k k k K k K k ηληληληληληληληληληηη 该线性方程组的系数矩阵为⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=---11211211111k k k k k D λλλλλλ,D 为德蒙行列式,又由)...2,1(k i i =λ互异有0≠D . 根据克拉默法则就有0=i η,即+11i i t α…+i i i i k k t α=0,再由i ir i αα,...,1线性无关得:)...2,1(0...21k i t t t i i i i k ===== ,故k i kr ir r αααα...,...,,...,1111线性无关.推论1 属于不同特征值的特征向量是线性无关的. 定理1 Q n n P ⨯∈与对角阵相似⇔Q 有n 个特征向量,它们是线性无关的.证明:Q 可以对角化⇔存在可逆矩阵21,(T T T =,…,)n T 使得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=-n QT T λλλ00211 ,即⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n T QT λλλ0021⇒ ,,(21QT QT …),n QT =(,,21T T λλ…n T λ,).因此Q 可以对角化⇔存在i T (2,1=i …n ,)P ∈使得i i i T QT λ=,也即Q 有n 个线性无关的特征向量.根据这个定理判定一个方阵是否可以对角化,必须从求解这个矩阵的特征多项式入手,虽然很直接,但考虑其计算量很大,加之特征值与特征向量只能分开求解,下面会介绍更简便的方法.推论2如过方阵n n P Q ⨯∈有n 个不同的特征值,那么该矩阵可对角化.证明:由Q 有n 个不同的特征值及引理1的推论有Q 有n 个线性无关的特征向量,再由定理1即有Q 可以对角化. 注意:该推论为对角化的充分条件.定理2 t μμμ,...,,21(互不相同)是B n n P ⨯∈的特征值,),...,2,1(t i P i =∈μ, B 可对角化⇔n t B E r ti i)1()(1-=-∑=μ (r 表示矩阵的秩).证明:0)(=-X B E i μ的基础解系的一组基向量的个数为:)(B E r n i --μ,我们可以得到关于i μ的线性无关的特征向量的个数是)(B E r n i --μ(),...,2,1t i =,再由引理1推出矩阵B 有))((1B E r n ti i --∑=μ个线性无关的特征向量.根据定理1就有:n 阶方阵B 可对角化⇔B 有n 个线性无关的特征向量 ⇔))((1B E r n ti i--∑=μ=n ,⇔n t B E r ti i)1()(1-=-∑=μ.定理3 n n P Q ⨯∈与对角矩阵相似的充要条件:)...,2,1(t i P i =∈λ且i i r Q E n =--)(λ (i r 表示i λ的代数重数).证明:设i λ的线性无关的特征向量为i ir i i βββ,...,,21,由引理1有:titr ir r βββββ,...,,...,,...,,111211线性无关.若n r r r t =+++...21,那么Q 就有n 个线性无关的特征向量⇔Q 可以对角化. 若Q 与对角矩阵相似,则Q 的属于不同特征值的特征向量总数一定为n . 否则根据定理1就可以推出t λλλ,...,,21线性相关,矛盾.相较于定理1,定理3的优点在于判定一个矩阵是否可以对角化着点于特征向量的重数,方便了许多,也易于计算.下面利用定理1结合矩阵的秩给出矩阵可对角化的另一判别方法.引理2 设n 阶方阵B A ,n n P ⨯∈,则有)()()(B r A r B A r +≤+.证明:先证)()(],[B rank A rank B A rank +≤……(2). 根据矩阵秩的定义有 ≤],[B A r n n 2⨯阶矩阵],[B A 的线性无关的行数≤方阵A 的线性无关的行数+方阵B 的线性无关的行数 ≤)()(B r A r +.对方阵矩阵=+A B ⎥⎦⎤⎢⎣⎡E E A B ],[,由(2)式有],[)(B A r A B r ≤+,所以)()()(B r A r B A r +≤+.引理3 对于n 阶方阵D C ,有n B r A r AB r -+≥)()()(.证明:先证⎪⎪⎭⎫⎝⎛≤⎪⎪⎭⎫ ⎝⎛=+D O T C r D O O C r D r C r )()(……(3),其中T 为任意n 阶方阵.显然当D C ,中有一个为O 时结论成立;另设q D r p C r ==)(,)(,则C 有p 阶子式01≠M ,D 有q 阶子式02≠M .于是⎪⎪⎭⎫⎝⎛D O T C 有q p +阶子式0*2121≠=M M M OM , 因此)()(D r C r q p D O T C r +=+≥⎪⎪⎭⎫ ⎝⎛.要证n B r A r AB r -+≥)()()(,只需证明:)()()(B r A r n AB r +≥+运用分块矩阵的初等变换有:⎪⎪⎭⎫⎝⎛-→⎪⎪⎭⎫ ⎝⎛-→⎪⎪⎭⎫ ⎝⎛→⎪⎪⎭⎫ ⎝⎛A O E B O A B E AB A O E AB OO E n nnn, 有初等变换不改变矩阵的秩以及式(3)有:)()()(B r A r A O E B r n AB r n +≥⎪⎪⎭⎫ ⎝⎛-=+.另证:令p A r =)(,则存在可逆矩阵D C ,使得CAD =⎪⎪⎭⎫ ⎝⎛O OO E p,若令⎪⎪⎭⎫ ⎝⎛--p n E O O O C 11-D =H ,则)(H r p n -=以及A H +=11--D C . 又因为任意矩阵左乘以与其行数相等的非奇异方阵或者右乘以与其列数相等的非奇异方阵不改变这个矩阵的秩,因此)(B r =r (B D C 11--) =)()(HB r AB r + ≤r(AB)+r(H) p n AB r -+≤)(.引理3的一般形式:(Syl 希尔维斯特不等式)设A ,B ,C n n P ⨯∈分别为t k k j j i ⨯⨯⨯,,矩阵,则)()()()(B r BC r AB r ABC r -+≥. 证明:要证)()()()(B r BC r AB r ABC r -+≥只需证明)()()()(BC r AB r B r ABC r +≥+, 因为分块矩阵的初等变换不会改变矩阵的秩,而⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛BC B O ABO E E O E C O E B O O ABC E O A E , 也即⎪⎪⎭⎫ ⎝⎛-→⎪⎪⎭⎫ ⎝⎛-→⎪⎪⎭⎫ ⎝⎛→⎪⎪⎭⎫ ⎝⎛BC B O AB B BC AB O B O AB ABC B O O ABC , 再有定理(3)就得)()(BC rank AB rank BC B O ABrank B O O ABC rank +≥⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛. 推论3设t B B B ,...,,21为数域P 上的n 阶方阵,则)...()1()(...)()(2121t t B B B r n t B r B r B r +-≤+++.定理4 设n 阶方阵Q n n P ⨯∈,21μμ≠,且0))((21=--Q E Q E μμ,则Q 可对角化.证明:由21μμ≠,0))((21=--Q E Q E μμ有矩阵Q 的特征值为1μ或2μ,根据引理2,引理3得:n Q E r Q E r =-+-)()(21μμ,从而Q 的特征向量(线性无关)共有n Q E r n Q E r n =--+--)()(21μμ个.由定理1即得矩阵Q 可对角化.定理4' 设n 阶方阵Q n n P ⨯∈,t μμμ,...,,21两两互不相等,若0))(())((121=--⋯---Q E Q E Q E Q E t t μμμμ则Q 与对角阵相似.证明:根据0))(())((121=--⋯---Q E Q E Q E Q E t t μμμμ有Q 的特征值在t μμμ,......,21中取得. 再由引理3的推论有n t Q E r Q E r Q E r t )1()(...)()(21-≤-++-+-μμμ,从而方阵Q 的线性无关的特征向量的个数为)(...)()(21Q E r n Q E r n Q E r n t --++--+--μμμ))(...)()((21Q E r Q E r Q E r tn t -++-+--=μμμn t tn )1(--≥n =.又因为n Q r ≤)(,故方阵Q 的线性无关的特征向量的个数为n ,由此矩阵Q 可对角化.推论4在定理4的前提条件下我们可以得到如下结论:n t Q E r Q E r Q E r t )1()(...)()(21-=-++-+-μμμ.定理4是判定矩阵相似与对角矩阵的充要条件,若矩阵阶数较高,计算量依然很大,特征值仍然需要计算,下面给出类似于定理4的充要条件.定理5 设t μμμ,...,,21(互不相同)是n n P Q ⨯∈的的特征值,重数分别为t s s s ,...,,21且n s s s t =+++...21,Q 可对角化⇔0)(1=-∏=Q E ti iμ.证明:先证明必要性Q 与V =⎪⎪⎪⎪⎪⎭⎫⎝⎛T μμμ21相似,则存在非奇异矩阵T 满足 ⎪⎪⎪⎪⎪⎭⎫⎝⎛==-t t E E E T TVT Q μμμ221111-T , 其中),...2,1(t i E i =为i s 阶单位矩阵,于是1)()(--=-T V E T Q E i i μμ=12211)()()(-⎪⎪⎪⎪⎪⎭⎫⎝⎛---T E E E T t t i i i μμμμμμ , 从而有∏∏=-=-=-ti it i iTV E T Q E 111)()(μμ⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---=∏∏∏i t t i iii i E E E T )()()(2211μμμμμμ1-T . 由于),...,2,1(0)(t j E ij j i ==-∏μμ,因此0)(=-∏Q E ii μ.再证充分性:对于n 阶矩阵Q ,存在可逆矩阵T ,使得1211--⎪⎪⎪⎪⎪⎭⎫⎝⎛==T J J J T TJT Q t, ),...,2,1(t i J i =是Jordan 块,若=j J j j E μ),...2,1(t j =,Q 就可以对角化,而1)()(--=-T J E T Q E i i μμ12211)()()(-⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---=T E J E J E J T t t i i i μμμ , 12211)()()()(-⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛---=-∏∏∏∏T E J E J E J T Q E i t t i iii i ii μμμμ. 所以,若0)(=-Q E i μ,则因T 可逆有),...,2,1(0)(t j J E ij i i ==-∏μ,又因为当j i ≠时,(j i μμ≠0)≠,)(j j i J E -μ可逆,所以0)(=-j j j J E μ,即),...,2,1(t j J E j j j ==μ. 引理4 nn PX ⨯∈,1∂,2∂…...m ∂是X 的关于特征值λ的特征向量,我们有imi i k ∂∑=1(m i k i ,...,2,1,=不全为0,P k i ∈)也是X 的关于λ的特征向量.证明:已知i i X ∂=∂λ,则i i i i k X k ∂=∂λ,也即i i i i k Xk ∂=∂λ,因此i mi i i m i i k k X ∂=∂∑∑==11λ,又i k 不全为0,因此01≠∂∑=i m i i k ,由特征向量的定义有i mi i k ∂∑=1是矩阵X 的属于特征值λ得特征向量.定理6 t μμμ,...,,21(互不相同)是n 阶矩阵Q 的所有特征值,它们的代数重数依次是t s s s ,...,,21,则方阵Q 与对角矩阵相似⇔),...,2,1()(t j s A r j j ==,∏≠-=ji i j Q E A )(μ.证明:先证必要性.Q 可对角化⇒存在可逆矩阵T 使得121),...,,(-=T Tdiag Q t μμμ,从而)(Q E A ji i j -=∏≠μ12211)()()(-≠≠≠⎪⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---=∏∏∏T E E E T j i t t i ji ij i i μμμμμμ11)(-≠⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛-=∏T O E O T t ji jj iμμ, 其中j O 为j s 阶0矩阵,j E 为j s 阶单位矩阵(),...,2,1(t j =. 因T 可逆,且j i μμ≠,所以有j j j j ji i j s E r E r A r ==-=∏≠)())(()(μμ),...,2,1(t j =.再证充分性:用反证法.假设方阵Q 不与对角矩阵相似,由几何重数≤代数重数得:至少存在一个整数q ,使得q q s n Q E r ->-)(μ,于是当q j ≠时,由引理3有∑∏≠≠---≥-=ji i ji i j n t Q E r Q E r s )2()())((μμ∑≠--->ji j n t s n )2()(∑≠----=ji i s n t n t )2()1(j j s s n n =--=)(.矛盾,假设不成立,故Q 与对角矩阵相似.定理7 t μμμ,...,,21(互不相同)是n 级方阵Q n n P ⨯∈的所有特征根,若对任意m +∈Z 满足)()(Q E r Q E r i m i -=-μμ,则矩阵Q 与对角矩阵相似.证明:设t μμμ,...,,21的重数分别为t s s s ,...,,21,由Hamilton Cayley -定理(高等代数第三版,高等教育)得:O Q E Q E Q E t s t s s =---)...()()(2121μμμ,再有引理3的推论就有t s t s s Q E r Q E r Q E r )(...)()(2121-++-+-μμμ))...()(()1(11t s t s Q E Q E r n t --+-≤μμn t )1(-=.对任意正整数m ,有)()(Q E r Q E r i m i -=-μμ,因此n t Q E r Q E r Q E r t )1()(...)()(21-≤-++-+-μμμ.从而有方阵Q 的线性无关的特征向量的个数为)(...)()(21Q E r n Q E r n Q E r n t --++--+--μμμ )](...)()([21Q E r Q E r Q E r tn t -+-+--=μμμn n t tn =--≥)1(.又n Q r ≤)(,从而Q 的线性无关的特征向量的个数小于或等于n ,因此Q 共有n 个线性无关的特征向量,再根据定理1就有矩阵Q 与对角矩阵相似. 接下来介绍最小多项式在矩阵对角化中的应用.定理8 n 阶方阵Q 与对角矩阵相似⇔矩阵Q 的最小多项式)(μQ m 无重根. 证明:先证必要性.Q 和对角阵相似⇒存在非奇异矩阵n n P T ⨯∈,满足1211--⎪⎪⎪⎪⎪⎭⎫⎝⎛==T T TVT Q n μμμ,从而有m m V T Q T =-1,令t μμμ,...,,21)(n t ≤是方阵Q 的互不相同的特征值,记))...()(()(21t f μμμμμμμ---= =t t t t s s s ++++--μμμ111.... 因为 T E s Q s Q s Q T T Q f T t t t t )...()(11111++++=----ET T s QT T s T Q T s T Q T t t t t 1111111...------++++==E s V s V s V t t t t ++++--111...)(V f =.又 E s V s V s V V f t t t t ++++=--111...)(⎪⎪⎪⎪⎪⎭⎫⎝⎛++⎪⎪⎪⎪⎪⎭⎫ ⎝⎛+⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=---t tt t n t t t n t t s s s s s s (111)2111121μμμμμμ⎪⎪⎪⎪⎪⎭⎫⎝⎛+++++++++=---k t n t n kt t kt t s s s s s s (1)112121111μμμμμμ 0)()()(21=⎪⎪⎪⎪⎪⎭⎫⎝⎛=n f f f μμμ. 所以0)(=Q f ,于是)()(μμf m Q ,然而)(μf 无重根,故)(μQ m 无重根.再证充分性:)(μQ m 的互不相同的根是t μμμ,...,,21,由)(μQ m 无重根就有:)(μQ m ))()...()((121t t μμμμμμμμ----=-,于是0))...()(()(21=---=Q E Q E Q E Q m t Q μμμ.令i i q Q E r =-)(μ,则i μ的特征子空间的维数为i q n -,因此Q 总共有s q n q n q n t =-++-+-)(...)()(21个线性无关的特征向量,且n s ≤. 又因为n t q q q t )1(...21-≤+++,故n q n q n q n s t ≥-++-+-=)(...)()(21.从而n s =,也即矩阵Q 有n 个线性无关的特征向量,由定理1就得Q 可以对角化.4某些特殊矩阵的对角化4.1 实对称矩阵的对角化问题实对称矩阵这种矩阵很特别,在诸多方面的到运用,如常用来研究对称变换,对线性变换进行分类.而研究对称矩阵的对角化,是进行分类的初步.引理5 ]每一个n 阶复矩阵都存在一个上三角矩阵与其相似,并且上三角矩阵主对角线上的元素为复矩阵的特征值.对任意n n C A ⨯∈,可逆矩阵T ,使得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=-n AT T λλλ211*,其中n λλλ,...,,21是矩阵A 的特征值. 引理6 实对称矩阵的特征值为实数.证明:设0λ实对称矩阵A 的一个特征值,则存在非零向量⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n x x x 21η,满足 ηλη0=A . 令⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=n x x x 21η,i x 称为i x 的共轭复数),...2,1(n i =,则ηλη0=A .观察下面式子ηηηηηηηη)()()('='=''='A A A A ,上式左边等于ηηλ'0,右边等于ηηλ'0, 故ηηλ'0=,0ηηλ'又0...11≠++='n n x x x x ηη,故00λλ=,即0λ是一个实数.引理7 设N M ,为n n ⨯实方阵,我们有如下结论:N M ,在实数域上相似⇔N M ,在复数域C 上相似.证明:必要性显然,下面证明充分性.N M ,在复数域上相似⇒∃n 级可逆复矩阵,使得NP P M 1-=.令iD A P +=,n n R D A ⨯∈, ,则ND DM NA AM iD A N M iD A ==⇒+=+,)()(.所以对任意λ属于R 都有)()(D A N D A λλ+=+ (4)记D A x h λ+=)((实数系多项式),因为0)(≠=+=P iD A i h ,所以0)(≠x h . 因此,D A λ+有有限个实数根,则存在η属于R ,使得0≠+D A η. 由(4)式得)()(1D A N D A M ηη++=-, 也即N M ,在实数域上相似. 定理9 ⑴n 级实对称矩阵A 的特称根全是实数⇔存在正交矩阵T ,满足D AT T AT T ='=-1,D 是上三角矩阵.⑵A 正交且特征值全是实数⇒A 是对称矩阵.证明:先证明必要性,根据引理5有,存在可逆矩阵P ,使得⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=-n AP P λλλ211*. 再根据引理7,矩阵如果在复数域上相似则一定在实数域上相似,因此可以令P QT =为实矩阵,Q 乃正交矩阵,T 是上三角矩阵且主对角上元素全是实数,于是就有⎪⎪⎪⎪⎪⎭⎫ ⎝⎛==---n T AP P T AQ Q λλλ21111*)( 由T 是上三角矩阵知他的逆1-T 也是上三角矩阵,再由上三角矩阵之积仍然是上三角矩阵知AQ Q 1-为上三角矩阵.再证充分性:A 为n 阶实矩阵,且存在正交矩阵Q 使得AQ Q AQ Q '=-1为上三角矩阵,即AQ Q AQ Q n '=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=-λλλ211*, 由此易知n λλλ,...,,21为实数且为A 的特征根.⑵由⑴容易得到AQ Q AQ Q '=-1D =为上三角矩阵(Q 是正交矩阵),又正交矩阵的积为正交矩阵,从而D 为正交矩阵. 因而1-='D D ,但是1-D 是上三角矩阵,而D '为下三角矩阵,故D 必为对角矩阵.从而A Q QD Q D Q Q QD A ='=''=''=')(,也即A 为对称矩阵.引理8 设A 是对称变换,1V 是-A 子空间,则1V 的正交补也是-A 子空间. 定理10 对任意n 级实对称矩阵A ,存在n 阶正交矩阵T ,使得AT T AT T 1-='为对角矩阵.证明定义A 是与A 对应的对称变换,只要证A 有一组标准正交基(n 个向量组成).下面用数学归纳法进行证明.当1=n 时结论明显成立.假设对1-n 结论成立. 对n 维欧氏向量空间n R ,1β为线性变换A 的一个特征向量,对应的特征值是1λ. 将1β单位化,并记为1α,再作1α的生成向量空间)(1αL 的正交补,记为1V ,由引理8有1V 是对称变换A 的不变子空间,他的维数为1-n ,显然A 限制在1V 上仍然是对称变换1V A ,根据假设1V A 有特征向量n ααα,...,,32做成1V 的标准正交基,从而n αα,...,1使n R 的标准正交基,又是A 的n 个特征向量.根据归纳假设定理得证. 例4.1 已知⎪⎪⎪⎪⎪⎭⎫⎝⎛----=0111101111011110A , 求正交矩阵T 使得AT T 1-为对角矩阵.解:第一步,求矩阵A 的特征值. 由μμμμμ111111111111--------=-A Eμμμμμμμμ1111100101011102---------=110101111)1(3μμ----=)3()1(3+-=μμ 由此有1(3重),-3为A 的特征值.第二步,求特征值1对应的特征向量. 将1=μ带入下式⎪⎪⎩⎪⎪⎨⎧=+--=-++-=-++-=+--.0,0,0,04321432143214321x x x x x x x x x x x x x x x x μμμμ (5) 得基础解系为)0,0,1,1(1=μ, )0,1,0,1(2=μ, )1,0,0,1(3-=μ..将基础解系正交化,得)0,0,1,1(1=β,)0,1,21,21(2-=β,)1,31,31,31(3-=β..再将上式单位化,有)0,0,21,21(1=η,)0,62,61,61(2-=η, )123,121,121,121(3-=η. .上式为属于特征值1(三重)的三个标准正交特征向量.同理可求得特征值-3的标准正交特征向量为)2/1,2/1,2/1,2/1(4--=η.特征向量4321,,,ηηηη构成4R 的一组标准交基,所求正交矩阵()4321,,,ηηηη''''=T , 此时⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=-31111AT T . 4.2幂等矩阵定理11幂等矩阵A 与对角矩阵⎪⎪⎭⎫ ⎝⎛O OO E r相似. 证明:根据A A =2有,矩阵A 的最小多项式)(λA m 整除λλ-2. 因02=-λλ无重根,由引理5 就有)(λA m 无重根,再由定理8就得矩阵A 可对角化.4.3对合矩阵定理12对合矩阵A 可对角化. 证明:E A =2⇒)(λA m 12-λ,易知12-λ=0无重根,根据引理5得)(λA m 无重根,再根据定理8,A 能够对角化.4.4幂幺矩阵引理9 λ是矩阵X 的任一特征根λ⇔是X 的最小多项式的根. 证明:用反证法假设0λ是矩阵X 的特征根而不是其最小多项式)(λX m 的根,则有1))(,(0=-λλλX m ,故存在多项式)(),(λϕλφ,使得1)()())((0=+-λλϕλλλφX m , 将X 带入上式有E X m X E X X X =+-)()())((0ϕλφ, 即有 E E X X =-))((0λφ.所以E X 0λ-可逆(即00≠-E X λ),与0λ是矩阵X 的特征根矛盾. 故假设不成立,定理得证.定理13幂幺矩阵A 与对角矩阵相似.证明:因为E A m =,所以矩阵A 的最小多项式)(λA m 整除1-m λ(m 为正整数),而1-m λ无重根,根据引理9就得)(λA m 无重根,再由定理8即得矩阵A 与对角矩阵相似.注意:幂幺矩阵A 与对角矩阵⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n λλλ21相似,其中),...,2,1(1n i mi ==λ. 4.5矩阵的逆、伴随矩阵的对角化定理14n n P A ⨯∈能够对角化⇒*1,A A -可对角化. 证明:(I)根据题设条件,存在非奇异矩阵T 满足⎪⎪⎪⎪⎪⎭⎫⎝⎛=-n AT T μμμ211 由矩阵T 可逆就有),...,2,1(0n i i =≠μ,从而⎪⎪⎪⎪⎪⎭⎫⎝⎛=-----1121111n T A T μμμ , 从而1-A 与对角矩阵相似. (II)由1*-=A A A 得⎪⎪⎪⎪⎪⎭⎫⎝⎛==---22111*1///μμμA A A T A T A T A T从而*A 也与对角矩阵相似.4.6某些正交矩阵的对角化 4.6.1二阶正交矩阵的对角化问题设⎪⎪⎭⎫⎝⎛=22211211b b b b A 是正交矩阵,根据正交矩阵的性质就有 ⎪⎩⎪⎨⎧=+=+=+01122211211222212221211b b b b b b b b , 从而二阶正交矩阵A 有两种形式⎪⎪⎭⎫ ⎝⎛-=ααααcos sin sin cos A 或⎪⎪⎭⎫⎝⎛-=ααααcos sin sin cos A . 定理15若⎪⎪⎭⎫⎝⎛-=ααααcos sin sin cos A ,则矩阵A 与对角矩阵相似.证明:A 的特征多项式为μαμαμαμαμ-----=-cos sin sin cos E A1cos 22+-=αμμ由0=-E A μ得矩阵A 的特征值为1cos cos ,1cos cos 2221--=-+=ααμααμ.当01cos 2≠-α时,容易得到21μμ≠,故正交矩阵A 有两个不同的特征值,容易看出此时正交矩阵A 有两个线性无关的特征向量,由定理1即有正交矩阵A 可对角化.当01cos 2=-α,即1cos ±=α时,⎪⎪⎭⎫ ⎝⎛=1001A 或⎪⎪⎭⎫⎝⎛--1001,此时正交矩阵A 显然与对角矩阵相似.定理16若 ⎪⎪⎭⎫⎝⎛-=ααααcos sin sin cos A ,那么A 与对角矩阵相似. 该定理的证明与定理⑴类似,在此不做赘述.4.6.2几类三阶正交矩阵的对角化定理17正交矩阵⎪⎪⎪⎭⎫ ⎝⎛=33322322110000b b b b b A 可对角化.证明:由A 是正交矩阵可得,111±=a .当111=b 时,⎪⎪⎪⎭⎫⎝⎛-=θθθθcos sin 0sin cos 0001A , )1cos 2)(1(cos sin 0sin cos 0012+--=----=-θλλλλθθθλθλλA E . i)当01cos 2≠-θ时,矩阵A 有三个不同的特征值,分别为1cos cos 21--=θθλ,1cos cos 22++=θθλ,13=λ.。
矩阵论文
上饶师范学院数学与计算机科学学院本科毕业论文论文题目:矩阵可对角化的充分必要条件专业:数学与应用数学班级:09数(3)学号:***************指导教师姓名:龚和林上饶师范学院数学与计算机科学学院2013年04 月摘要在高等代数中,方阵A可对角化当且仅当它可相似于对角矩阵,即存在一个可逆矩阵P,使APP1 是对角矩阵。
因为对角矩阵的特征值与特征向量是已知的,从而,对矩阵的对角化有积极的意义。
本文给出了矩阵四种可对角化的充分必要条件和相应的证明。
关键词方阵;特征值;特征向量;对角化2目录绪论.................................... 错误!未定义书签。
1 矩阵可对角化的概念 (5)1.1特征值、特征向量的概念 (5)1.2矩阵可对角化的概念 (6)2 矩阵可对角化的充分必要条件 (3)2.1利用特征向量判断矩阵可否对角化 (3)2.2 利用特征根的性质判断矩阵可否对角化 (4)2.3利用最小多项式判定矩阵可否对角化 (5)2.4 利用线性变换相关知识判断矩阵可否对角化 (5)3 矩阵可对角化的应用 (7)结论 (16)参考文献 (16)致谢 (17)34绪论矩阵是高等代数中的重要组成部分,是许多数学分支研究的重要工具。
而对角矩阵作为矩阵中比较特殊的一类,其形式简单,研究起来也非常方便。
研究矩阵的对角化及其理论意义也很明显,矩阵相似是一种等价关系,对角化相当于对一类矩阵在相似意义下给出了一种简单的等价形式,这对理论分析是方便的。
相似的矩阵拥有很多相同的性质,比如特征多项式、特征根、行列式……线性代数中矩阵是否可以对角化,是矩阵的一条很重要的性质。
矩阵对角化也是《高等代数》和《线性代数》中矩阵理论这一部分的主要内容。
人们对此研究得出了很多有用的结论。
诸如一些充要条件:n 阶方阵A 可以对角化的充要条件是它有n 个线性无关的特征向量;方阵A 可以对角化的充要条件是它的最小多项式没有重根;还有复方阵A 可以酉相似于对角形矩阵的充要条件是它为正规矩阵,此外,还有一些充分条件。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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。