矩阵分析二章_第1页
矩阵分析二章_第2页
矩阵分析二章_第3页
矩阵分析二章_第4页
矩阵分析二章_第5页
已阅读5页,还剩106页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、 第二章第二章 矩阵与矩阵的矩阵与矩阵的jordan标准形标准形 矩阵的基本概念矩阵的基本概念定义定义:设:设为数域为数域 上的多项式,则称上的多项式,则称 ( )(1,2,;1,2, )ijaim jnf111212122212( )( )( )( )( )( )( )( )( )nnmmmnaaaaaaaaaa为多项式矩阵或为多项式矩阵或 矩阵。矩阵。定义定义 如果 矩阵 中有一个 阶 子式不为零,而所有 阶子式(如果有的话)全为零,则称 的秩为 ,记为零矩阵的秩为0。定义定义 一个 阶 矩阵称为可逆的,如果有一个 阶 矩阵 ,满足这里 是 阶单位矩阵。 称为 矩阵的逆矩阵,记为 。( )

2、ar(1)r 1r ( )arrank ( )arnn( )b( ) ( )( ) ( )abbaeen( )b( )a1( )a定理定理 一个 阶 矩阵 可逆的充要必要是 一个非零的常数。定义 下列各种类型的变换,叫做 矩阵的初等变换:(1) 矩阵的任二行(列)互换位置;(2) 非零常数 乘矩阵的某一行(列);(3) 矩阵的某一行(列)的 倍加到另一行(列)上去,其中 是 的一个多项式。 对单位矩阵施行上述三种类型的初等变换便得相应得三种 矩阵得初等矩阵 n( )adet( )ac( ) ( , ), ( ( ), ( , ( )p i j p i cp i j( ) 定理定理 对一个 的

3、矩阵 的行作初等行变换,相当于用相应的 阶初等矩阵左乘 。对 的列作初等列变换,相当于用相应的 阶初等矩阵右乘 。定义定义 如果 经过有限次的初等变换之后变成 ,则称 与 等价,记之为mn( )am( )a( )a( )an( )a( )a( )b( )( )ab( )b定理定理 与 等价的充要条件是存在两个可逆矩阵 与 ,使得( )a( )b( )p( )q( )( ) ( ) ( )bpaq 矩阵矩阵smith标准形的存在性标准形的存在性 定定 理理 任意一个非零的 型的 矩阵都等价于一个对角矩阵对角矩阵,即mn12( )( )( )( )00rddad 其中 是首项系数为1的多项式且称这

4、种形式的 矩阵为 的smith标准形。 称为 的不变因子。1,( )ird1( )( )(1,2,1)iiddir( )a12( ),( ),( )rddd( )a例 122221( )1a将其化成smith标准形。2222222223232432432110( )11101000000a解:32322432100100000001000000(1) 例 22(1)( )(1)a 将其化成smith标准形。2(1)( )(1)a 解:2(1)(1)(1)(2)1 32222(1)2021(1)(1)11(1)(1) 例 3将其化为smith标准形。222223232123( )43532344

5、21a22222421( )32321234353234a解:22222222421437333443532344214373334210222222210421437333412024133437342222221200104313410001043134222232100010344311000100132210001000110001000(1)(1)将其化为smith标准形。例 411( )1aaaaa100010( )001000aaaaa解:221000()1000100010000()10001000aaaaaaa223100001()0001000100001()000()10

6、00aaaaaa331000010000()100010000100001()000aaaa4100001000010000()a矩阵标准形的唯一性定定 义义: 为一个为一个 矩阵且矩阵且 对对于任意的正整数于任意的正整数 , , 必有非零的必有非零的 阶子式,阶子式, 的全部的全部 阶子式的最大公因式阶子式的最大公因式 称为称为 的的 阶阶行列式因子行列式因子。 ( )a( ( )rank ark1kr( )ak( )ak( )kd( )ak显然,如果 ,则行列式因子一共有 个。例 1 求的各阶行列式因子。解:( ( )rank arr22221( )1a由于 ,所以 。显然 而且其余的7各

7、2 阶子式也都包含 作为公因子,所以另外(1, )1 1( )1d2222321(1)( )1(1)( )1fg( ( ), ( )fg32323( )( )ad 2( )d注意 :观察 三者之间的关系。定理: 等价(相抵) 矩阵有相同的各阶行列式因子从而有相同的秩。 设 矩阵 的smith标准形为123( ),( ),( )ddd( )a12( )( )( )( )00rddad容易计算上面的标准形的各阶行列式因子为1121212( )( )( )( )( )( )( )( )( )rrddddddddd显然有:112211( )( )( )( )( )( )( )( )rrrddddddd

8、d由于 与上面的smith标准形具有相同的各阶行列式因子,所以 的各阶行列式因子为而又是由这些行列式因子唯一确定的,于是我们得到定 理: 的smith标准形是唯一的。例 1 求下列 矩阵的smith标准形。( )a( )a12( ),( ),( )rddd12( ),( ),( )rddd( )a2222121000000(1)0(1)00000(2)nacacaca00120120(3)12002000解 :(1)容易计算出12224434( )1,( )(1)( )(1) ,( )(1)dddd 122234( )1,( )(1),( )(1),( )(1)dddd 221(1)( )(1

9、)(1)a (2)首先观察此矩阵的元素排列规律,显然下面看 阶行列式因子。有一个阶子式要注意,即( )()nnda1n 1n 121 211nncacc ccac容易计算出 从而 121121( )( )( )1( )1,( )1,( )1,( )()nnnnddddddda1( )1nd111()na(3) 4111(2)定理定理 矩阵 与 等价的充要条件是对于任何的 ,它们的 阶行列式因子相同。定理定理 矩阵 与 等价的充要条件是 与 有相同的不变因子。( )a( )bkk( )a( )b( )b( )a与一般的数字矩阵一样,我们有下面的推论:推论 矩阵 可逆的充要条件为 与单位矩阵等价。

10、推论 矩阵 可逆的充要条件为 可以表示成一系列初等矩阵的乘积。 ( )a( )a( )a( )a初等因子和矩阵的相似初等因子和矩阵的相似设 矩阵 的不变因子为在复数域内将它们分解成一次因式的幂的乘积:( )a12( ),( ),( )rddd 11112221221211221212() ()()()()()()()()ssrsrreeeseeeseeersddd其中 是互异的复数, 是非负整数。因为 ,所以满足如下关系1,sije1|( )(1,1)iiddir112111222212000rrssrseeeeeeeee定义 在上式中,所以指数大于零的因子称为 矩阵 的初等因子( )a()

11、,0,1, ,1,ijejijeir js例 如果 矩阵 的不变因子为( )a 1222323341(1)(1) (1)(1) (1) (2)dddd 则 的初等因子为( )a2, ,1, 2323(1) ,(1) ,(1) ,(1) ,(2)例 如果 矩阵 的秩为4,其初等因子为( )a5 62233, ,1,(1) ,(1) ,() i 3() i,求 的smith标准形。( )a解:首先求出 的不变因子( )a 233342321(1) () ()(1)(1)1diiddd 从而 的smith标准形为定理 阶 矩阵 与 等价的充要条件是它们有相同的秩且有相同的初等因子。( )a22323

12、1000000(1)0000( )00(1)000000(1)(1)00000000a n( )a( )b定理 设 矩阵为准对角形矩阵,则 与 的初等因子的全体是 的全部初等因子。此定理也可推广成如下形式:( )( )( )bac( )b( )c( )a定理 若 矩阵则 各个初等因子的全体就是 的全部初等因子。12( )( )( )( )taaaa12( ),( ),( )taaa( )a例 1 求 矩阵的初等因子,不变因子与标准形。解:记22000000( )00(1)10022a21223( ),( ),(1)1( )22aaa那么对于 ,其初等因子为 由上面的定理可知 的初等因子为因为

13、的秩为4,故 的不变因子为123( )00( )0( )000( )aaaa3( )a,1,1 ( )a, , ,1,1,1 ( )a( )a 4321(1)(1),(1),1dddd 因此 的smith标准形为( )a1000000( )00(1)0000(1)(1)a 例 2 判断下面两个 矩阵是否等价?22223141( )1122a22122( )2232211b例 3 求下面 矩阵不变因子1000100015432例 4 求下列 矩阵的行列式 因子与不变因子121001000001nnaaaa数字矩阵的相似与 矩阵的等价定理: 设 是两个 阶的数字矩阵,那么 与 相似的充分必要条件为

14、它们的特征矩阵 与等价。定义: 对于数字矩阵 ,我们称 的不变因子为 的不变因子,称 的初等因子为 的初等因子。 ,a bnabiaibaiaaiaa 对于任何一个数字矩阵 所以 ,于是可得下面两个定理定理: 两个同阶的方阵 相似的充分必要条件是它们有相同的初等因子。定理:两个同阶的方阵 相似的充分必要条件是它们有相同的行列式因子(或不变因子)。例 设 ,证明:,a b,a b,aian()rankian0(1) 阶矩阵与11aaaanaaba相似;(2) 阶矩阵与n11aaaa不相似。 矩阵的jordan标准形定义: 称 阶矩阵11aaban为jordan块。设 为jordan块,称准对角形

15、矩阵111iiiiiinnaaja 12,sj jj为jordan标准形矩阵。由前面的例题和定理可知jordan块的初等因子为,从而jordan标准形矩阵的初等因子为12sjjjj()inia1212() ,() ,()snnnsaaa于是可以得到下面的定理定理: 设 的初等因子为则,这里,n naca1212() ,() ,()snnnsaaaaj12sjjjj其中 我们称 是矩阵 的jordan标准形。特别地,我们有定理: 可以对角化的充分必要条件是111iiiiiinnaaja ,(1,2, )isjaaa的初等因子都是一次因式。例 1 求矩阵的jordan标准形。解: 先求出 的初等因

16、子。对 运用初等变换可以得到110430102a aia所以 的初等因子为211043010211(1) (2)ia2(1) ,2a故 的标准形为或a110010002j200011001j例 2 求矩阵的jordan标准形。解: 先求出 的初等因子。对 运用初等变换可以得到112336224aaia1123362241(2)ia 所以 的初等因子为a, ,2 故 的jordan标准形为或000000002j000020000ja求jordan标准形的另一种方法:特征矩阵秩的方法.具体操作步骤:(1)先求出该矩阵的特征多项式及其特征值(2)其jordan标准形的主对角线上都是 的特征值,并且特

17、征值 在主对角线上出现的次数等于 作为特征根的重数。对于每个特征值 ,求出以它为主对角元的各级jordan 块的数目 ,首先求出 那么以 为主对角元的 jordan 块的总数是aiii()in()irank aii这里 为该矩阵的阶数,而以 为主对角元的 级 jordan 块的数目是依次先求出直至满足条件()()iinnrank ainit11( ;)()()2()tiittiin trank airank airank ai(1;),(2;),( ;)iiinnn t为止。(3)根据第二步求出的各级jordan块的数目,就可以写出 的一个jordan标准形。例 1 用矩阵秩的方法求出矩阵的j

18、ordan标准形。()(1;)(2;)( ;)iiiinnnn ta2321822143a解: 先求出 的特征多项式及其特征值。 对于特征值 ,它是 的1重根,从而 在 的 jordan 标准形的主对角线上出现一次,因此 中主对角元为1 的jordan块只有一个且它为一阶的。a2( )232182(1)(3)2143fia11( )f1aj对于特征值 ,先求 所以 从而23(3 )rank ai13213231520842146000ai(3 )2rank ai2()2321nn特征值 是 的两重根,从而 在 的jordan标准形 的主对角线上出现两次,因此 中主对角元为 3的jordan块只

19、有一个且它为二阶的。故 的标准形为或2aj( )f23ja100031003j310030001j例 2 用矩阵秩的方法求出矩阵的jordan标准形。解:首先求出其特征值,显然其特征多项式为1234012300120001a4( )(1)fia所以 为 的4重根,从而 在 的 jordan 标准形 的主对角线上出现四次,下面计算 中主对角元为1 的jordan块的数目,先计算 , 容易得到那么中主对角元为 的jordan块数是由此立即可得其jordan标准形为( )fajj( )()431nnrank ai1()rank ai()3rank ai11100011000110001a如何求相似变

20、换矩阵? 设 阶方阵 的jordan标准形为 ,则存在可逆矩阵 使得najp1p apj,称 为相似变换矩阵。对于相似变换矩阵的一般理论我们不作过多的讨论,只通过具体的例题说明求 的方法。例 1 求方阵的jordan标准形及其相似变换矩阵 。pp308316205ap解: 首先用初等变换法求其jordan标准形:230831620510001000(1)ia故 的初等因子为从而 的jordan标准形为 再求相似变换矩阵: 设所求矩阵为 ,则 ,对于 按列分块记为21,(1)aa100011001jp1p apjp123,pxxx于是有从而可得1231231231223,100,011001,a

21、pa xxxaxaxaxpjxxxxxxx 1122323,axxaxxaxxx 整理以后可得三个线性方程组前面的两个方程为同解方程组,可以求出它们的一个基础解系:可以取 ,但是不能简单地取,这是因为如果 选取不当会使得第三个非齐次线性方程组无解。由于1232()0()0()ia xia xia xx120,1,0,2,0,1tt 11x22x2x12, 的任意线性组合都是前两个方程组的解,所以应该取 使得第三个非齐次方程有解,即其系数矩阵与增广矩阵有相同地秩,容易计算出其系数矩阵的秩为1,从而应该使得增广矩阵 的秩也为1。即21122xkk2,iax22124082,306204kiaxkk

22、容易看出只需令 就会使得上述矩阵的秩为1,于是再由第三个方程解出一个特解为123,2kk 212324,3,2tx31,0,0tx ,那么所求相似变换矩阵为例 2 求方阵的jordan标准形及其相似变换矩阵 。123041,130020pxxx126103114a 解: 首先用初等变换法求其jordan标准形:21261311410001000(1)ia故 的初等因子为从而 的jordan标准形为 再求相似变换矩阵: 设所求矩阵为 ,则 ,对于 按列分块记为a21,(1)a100011001j123,pxxxp1p apj于是有从而可得1231231231223,100,011001,apa

23、xxxaxaxaxpjxxxxxxx1122323,axxaxxaxxx整理以后可得三个线性方程组前面的两个方程为同解方程组,可以求出它们的一个基础解系:可以取 ,但是不能简单地取,这是因为如果 选取不当会使得第三个非齐次线性方程组无解。由于1232()0()0()ia xia xia xx 121,1,0,3,0,1tt 11x22x2x12,的任意线性组合都是前两个方程组的解,所以应该取 使得第三个非齐次方程有解,即其系数矩阵与增广矩阵有相同地秩,容易计算出其系数矩阵的秩为1,从而应该使得增广矩阵的秩也为1。即21122xkk2,iax122122263,113113kkiaxkk容易看只

24、要 就会使得上述增广矩阵的秩为1。令 ,于是再由第三个方程解出一个特解为12kk121kk2122,1,1tx32,0,1tx ,那么所求相似变换矩阵为从而有123122,110011pxxx1100011001p ap一般地,设 ,则存在 阶可逆矩阵 使得其中 为jordan块,记这里n nacnp121tjjp apjij12,tpp ppin nipc那么有记 ,又可得121122,1,2,tttiiiap apappj p jpjappjit12,iiiiinpxxx112121iiiiiiiiiiininiinaxxaxxxaxxx注意: 是矩阵 的对应于特征值 的特征向量,特征向量

25、 的选取应该保证特征向量 可以求出,同样特征向量 的选取应该保证特征向量 可以求出,依此类推,并且使得线性无关。jordan标准形的某些应用例 1 对于方阵1 ixai1 ix2ix2ix3ix12,iiiinxxx126103114a 求 。解:首先用初等变换法求其jordan标准形:10a21261311410001000(1)ia故 的初等因子为a21,(1)从而 的jordan标准形为再求相似变换矩阵 且 ,那么 按照前面例题的方式,容易计算出 a100011001jp1p apj10101apj p122110011p从而10101122100110011001100110219206011210930113101031apj p 例 2 求解常系数线性微分方程组解: 令1132123313383825dxxxdtdxxxxdtdxxxdt 112233308316 ,205dxdtxdxdxaxxd

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论