矩阵初等变换及应用

矩阵初等变换及应用王法辉摘要：矩阵初等变换是高等代数的重要组成部分。本文对初等变换进行了研究探讨，详细介绍了与矩阵初等变换有关的基础知识。在阐述矩阵初等变换方法及应用原理的基础上，首先重点讨论该方法在解决高等代数相关计算问题上的应用，如求多项式的最大公因式、求逆矩阵解矩阵方程、求解线性方程组、判定向量的线性相关性、化二次型为标准型、求空间的基等。尤其是利用矩阵初等变换法求空间的基(解空间、特征子空间、核、值域等)的问题的计算,以具体实例生动的展示出问题的内在关系,最后给出了该方法在解决实际问题中的应用。本文理论分析与实际相结合，凸现了矩阵初等变换法直接、便利、有效的威力与作用。关键词：矩阵初等变换；最大公因式；线性相关性；二次型；空间的基1导言在线性方程组的讨论中我们看到，线性方程组的一些重要性质反映在它的系数矩阵和增广矩阵的性质上，并且解方程组的过程也表现为变换这些矩阵的过程。在数学的学习和应用中，矩阵理论是高等代数的重要组成部分，矩阵初等变换方法更是贯穿高等代数理论的始终。应用初等变换证明命题过程容易被接受，同时也是解决高等代数相关计算问题最直接、便利、有效的方法。此外，还有大量的各种各样的，表面上看完全没有联系的问题的解决,都可以通过相同的方法实现：矩阵的初等变换。因此，对矩阵初等变换方法及应用进行探讨，无疑是十分必要和重要的。目前，有许多文献涉及到对矩阵初等变换方法该的讨论，但比较零散。在研读文献的基础上，对矩阵初等变换的内涵进一步挖掘，使矩阵初等变换方法的威力作用得以充分展示是重要也是必要的。2矩阵及其初等变换2.1矩阵由nm个数)j,,,2,1(miaij（i=1,2,,j=1,2,n，）排成m行n列的数表mnmmnnaaaaaaaaaA212222111211称为m行n列的矩阵，简称nm矩阵。2.2矩阵的初等变换及初等矩阵矩阵有行列之分，因此有如下定义定义1矩阵的初等行（列）变换是指如下三种变换(1)交换矩阵某两行（列）的位置，记为jirr)(jicc；(2)把某一行（列）的k倍加到另一行（列）上，记为jikrr)(jikcc；(3)用一个非零常数k乘以某一行（列），记为ikr)(ikc，k0；矩阵的初等行变换及初等列变换统称为矩阵的初等变换。定义2由单位矩阵E经过一次初等变换得到的方阵称为初等矩阵。有以下3种形式(1)互换矩阵E的i行和j行的位置，得1101111011),(jiP；(2)用数域P种非零数c乘E的i行，得111))((cciP；(3)把矩阵E的j行的k倍加到i行，有1111))(,kkjiP（。定义3如果B可以由A经过一系列初等变换得到，矩阵A与B称为等价的。2.3矩阵初等变换的若干性质矩阵的初等变换改变了矩阵的元素,但矩阵初等变换具有以下性质（1）对矩阵A施行初等行（列）变换，其列（行）向量组之间的线性关系保持不变。（2）对矩阵A施行初等行变换相当于左乘相应的初等矩阵，施行初等列变换相当于右乘相应的初等矩阵。（3）可逆矩阵可以表示成一系列初等矩阵的乘积。初等矩阵的逆矩阵也是初等矩阵。（4）初等变换不改变矩阵的秩。3矩阵初等变换在高等代数计算问题中的应用矩阵初等变换与线性方程组的求解密不可分，不仅给解线性方程组带来了极大方便，同时也发展和完善了矩阵理论本身，极大丰富了矩阵理论的应用领域。矩阵的初等变换方法更是贯穿高等代数理论的始终，在高等代数有关理论的证明及相关计算问题中更是起着巨大的作用。3.1求多项式的最大公因式3.1.1基本概念以][xP表示数域P上的一元多项式环。定义1（最大公因式）设)()(xgxf，是][xP中两个多项式，][xP中多项式)(xd称为)()(xgxf，的一个最大公因式，如果它满足(1))(xd是)()(xgxf，的公因式；(2))()(xgxf，的公因式全是)(xd的因式。定义2以][xP中的一元多项式为元素的矩阵称为多项式矩阵。定义3以下3种变换称为多项式矩阵的初等行变换(1)交换多项式矩阵的某两行；(2)用零次多项式(P中不等于零的数)乘以多项式矩阵的某一行；(3)用一个多项式乘以多项式矩阵的某一行再加到另一行。且分别称以上三种变换为第1类，第2类，第3类多项式矩阵的初等行变换。所说的初等行变换总是指多项式矩阵的行初等变换，所说的矩阵总是指多项式矩阵。3.1.2主要结果在高等代数中，求数域P上两个多项式的最大公因式通常是利用辗转相除法，当多项式的次数较高时，辗转相除法计算较繁琐。由于多项式辗转相除法主要表现为系数间的运算，因此通常利用分离系数法，使运算相对简化。同样地，为了简化求多项式最大公因式的运算，考虑将要求最大公因式的两个多项式的系数与二行矩阵表示式对应起来。考虑][xP中的多项式)0()()0()(01110111mmmmmnnnnnbbxbxbxbxgaaxaxaxaxf其中iajbP(0,1,2,;0,1,2,)injm,引入如下记号当mn时,（)(xf,)(xg）011011bbbbaaaannnn；当mn时，（)(xf,)(xg）010111000bbbaaaaaammmnn。由于多项式的最大公因式具有以下基本性质（1）（)(xf,)(xg）=（)(xg,)(xf）；（2）若（)(xf,)(xh）=1，则（)(xf,)(xg）=（)(xf,)()(xhxg）；（3）（)(xf,)(xg）=（)()(xkgxf,)(xg）,Pk；因此，如上引入的二行矩阵反映了以下事实（1）交换二行矩阵两行的位置，得到的矩阵仍然对应这两个多项式的最大公因式；（2）二行矩阵某一行的k倍加于另一行得到的矩阵仍然对应这两个多项式的最大公因式。上述事实意味着数域P上多项式的最大公因式（)(xf,)(xg）可以利用二行矩阵进行初等行变换求得。具体实施步骤为（1）根据多项式的系数作出（)(xf,)(xg）对应的二行矩阵；（2）利用第1、2类初等行变换使得二行矩阵中的行出现端首(左端或右端)为0；（3）向左(或向右)平移二行矩阵中某行，使得这一行端首的0去掉。这表明（)(xf,)(xg）的次数在降低。反复利用（1）、（2）、（3）直到出现二行矩阵的两行元素对应成比例为止。3.1.3计算举例例1已知数域P上的一元多项式7787)(346xxxxxf，7373)(235xxxxg求))(),((xgxf。解构造二行矩阵A并实施初等行变换070370373140731400070370377087012131rrA12149070370300731407314rr首项不为零将第一行元素轮换使其14147070033990707022将第二行元素轮换使其首项不为零1228271414777070000000332727999970700707002222rr212431477000002727997070022rr将第一行元素轮换使其不为零770000027270700700将第二行元素轮换使其首项不为零1212777000000000000272770070007007000rr700000001001000第二行除以（）第二行元素共轮换过3次，所以最大公因式为1)(3xxd。例2求多项式3442)(234xxxxxf，3452)(23xxxxg，6116)(23xxxxh的最大公因式。解构造三行矩阵A并进行初等行变换12r12443109000254302543161160161160rA01090025430161160对第二行进行轮换，使其首项不为21312,1090010900051430514300062060620600rrrr轮换001031000054251400090100131010053514100901131221,61,51rrrrrr0003-10006-2009-01003100062000901103,7532轮换rr1212039003900026000260001300013000rr轮换所以3))(),(),((xxhxgxf。3.2求逆矩阵解矩阵方程3.2.1可逆矩阵定义若对n级矩阵A有n级矩阵B使EBAAB则称A是可逆的，B称为A的可逆矩阵。其中E为n级单位矩阵。3.2.2初等变换求逆的原理和步骤由于可逆矩阵A可表示为一系列初等矩阵的乘积，故由EAA1有EAPPEAPPSs11因此有如下求逆步骤（1）构造nn2的矩阵EA｜；（2）对上述矩阵实行初等行变换，当用初等行变换把A化为单位阵，则E的位置变成A的逆矩阵，即EA｜1AE｜需要指出的是在此过程中只能用初等行变换。如果用列变换，则需把E置于A的下方变成2nn矩阵且只能使用列变换把A化为单位矩阵，同时E化为A的逆矩阵，即EA初等列变换1AE利用与求逆矩阵相同的原理，矩阵初等变换可用于解矩阵方程。3.2.2计算举例例1求A=032203120的逆矩阵。解构造矩阵，由3204900102030011201000320102030011202323rrEA｜649100001120010203320490001120010203232192rrrr13)2(264910032401043600164910064802012918003123132AErrrrrr｜得1A=649324436例2设A=031221312，B=520211，求X使得BAX解构造矩阵EAC｜并实施初等行变换C=520310222111312131221,2,rrrrrr50050131300222123323,5rrrrr23100100100222132122rrr231001001024001=BAE1｜得X=1AB=2310243.3求解线性方程组3.3.1有关概念与结论考虑n元线性方程组nnnnnnnnnbxaxabxaxabxaxa