《用matlab实现矩阵的对角.doc》由会员分享，可在线阅读，更多相关《用matlab实现矩阵的对角.doc（11页珍藏版）》请在金锄头文库上搜索。
1、用matlab实现矩阵的对角 对角阵在实际上的应用特别广泛，对角阵解决现实问题上很方便，通过对角矩阵可以最简单地处理物力问题，也可以解出线性方程组的解；最普遍的是可以直接知道相似矩阵的行列式值，秩，特征值等，所以可以说研究对角化问题是特别重要。对角化的最快，最方便的方法是利用matlab软件。1 2 3 形式的矩阵为对角矩阵 一般 n （空白处为零）。 在相似变换下，方阵A的许多重要性质（如 行列式，秩，特征值等）保持不便，因此我们可以通过相似变换将矩阵A化简，并利用化简后的矩阵来研究与矩阵有关的问题。 我们讨论的主要问题是：对于n阶方阵A，是否存在可逆方阵P，使P-1AP为对角矩阵。这就是矩
2、阵的对角化问题。 2 0 0 -1 1 0 通过相似变换下简化为对角矩阵B= 0 1 0 -4 3 0A=矩阵1 0 01 0 2通过矩阵B可以直接知道矩阵的行列式，秩，特征值，对应方程组的解等重要性质。 但任何矩阵不一定可以对角化。一个矩阵是否可以对角化有如下的判断方法： 1）判断A是否实对称矩阵，茹是一定可对角化，因为A是实对称矩阵，则有（1）A的全部特征值是实数。（2）A的不同特征值对应的特征向量正交。（3）A一定相似于对角矩阵，且存在正交矩阵T，使得T-1AT=TTAT=。的对角元素是A的特征值。 1 -1 2 -1 实对称矩阵A= -1 1 3 -2 是否可对角化？茹是将矩阵A对角
3、2 3 1 0 -1 -2 0 1化。 运用matlab程序来实现这个问题： 程序如下： A=1 -1 2 -1；-1 1 3 -2；2 3 1 0；-1 -2 0 1； V，D=eig（A） 运行结果： V= 0.4412 -0.2042 -0.8328 0.2647 0.6012 0.1266 0.4853 0.6221 -0.5683 0.4886 -0.2227 0.6234 0.3477 0.8388 -0.1462 -0.3927D= -3.7266 0 0 0 0 0.9416 0 00 0 1.9420 00 0 0 4.8430程序说明：D 对角线上的元素为A的特征值，V为相
4、对应的特征向量所构成的矩阵。可以看出特征值都是实数，每个特征值对应于一个特征向量，对应的特征向量是正交的，特征向量所构成的矩阵是正交阵。如果不用matlab软件的话实现这种（3阶以上的矩阵）问题是特别复杂及麻烦。2）求A的特征值，若n个特征值互异，则A一定可对角化。 0 1 0 - 1 0 A= 0 0 1因为 |A-E|= 0 - 1= -（1+） -6 0 -6 -6 -11 -6-（2+）（3+）=0。所以特征值分别为1=-1，2=-2，3=-3； 用Matlab程序来实现： A=0 1 0；0 0 1；-6 -11 -6； D=eig（A） 运行结果： D=-1，-2，-3； 结果一样
5、，A有3个互不相同的特征值，故可对角化；这两种方法相互比较，用Matlab程序好计算这种问题。3）求A的特征向量，若有n个线性无关的特征向量，则A可以对角化，否则不可以对角化。 0 1 1 -1 -1A= 1 0 1因为 |I-A| = -1 -1 =（-2）（+1）2 1 1 0 -1 -1 于是A的特征值为：1=2，2=-1 （二重），每个特征值对应的特征向量为X1=（1 ，1，1）T， X2=（-1， 1，0）T， X3=（-1，0，1）T，这三个向量相互线性无关，故矩阵A可对角化。 程序如下：A=0 1 1；1 0 1；1 1 0；V，D=eig（A） V= 1.0000 -1.000
6、0 -1.0000 1.0000 1.0000 0.0000 1.0000 0.0000 1.0000D= 2.0000 0 0 0 -1.0000 0 0 0 -1.00004）方阵A可对角化的充要条件：A的每个重特征值对应的线性无关的特征向量的个数等于该特征值的重数。 3 2 -1 x-3 -2 1 A= -2 -2 2 的特征多项式是 2 x+2 -2 = x3-12x 3 6 -1 -3 -6 x+1+16=（x-2）2（x+4）=0。特征根是2，2，-4。 特征根-4对应的特征向量为（1/3，-2/3，1）。 特征根2对应的特征向量为（-2，1，0）（1，0，1）。 可以看出特征向量
7、的个数等于该特征值的重数。程序如下：A=3 -2 -1；-2 -2 2；3 6 -1；V，D=eig（A） V= 0.8890 0.2673 0.1654 -0.2540 -0.5345 0.3737 0.3810 0.8018 0.9127D= 2 0 0 0 -4 0 0 0 2从上面可以看出这两种方法来解得的结果一样，很容易知道用Matlab程序来处理即方便又速度快。下面讲将矩阵对角化的方法：1） 令A是数域F上的一个n阶矩阵，如果A的特征多项式fA（x）在F内有n个单根，那么存在一个n阶可逆矩阵T，使12。T-1AT=3 n 3 2 -1 x-3 -2 1设A= -2 -2 2的特征多
8、项式是 2 x+2 -2 = 3 6 -1 -3 -6 x+1 X3-12x+16=（x-2）2（x+4）；特征根是2，2，-4；对于特征根-4，求出齐次线性方程组 -7 -2 1 x1 0 2 -2 -2 x2 = 0 的一个基础解系（1/3，-2/3，1）。-3 -6 -3 x3 0 对于特征根2，求出一个基础解系（-2，1，0）（1，0，1）。由于基础解系所含解向量的个数都等于对应的特征根的重数，所以A可以对角化，取 1/3 -2 1 -4 0 0 T= -2/3 1 0 那么T-1AT= 0 2 0 。 1 0 1 0 0 2 Matlab程序如下：A=3 2 -1；-2 -2 2；3
9、 6 -1；V，D=eig（A） V= 0.8890 0.2673 0.1654 -0.2540 -0.5345 0.3737 0.3810 0.8018 0.9127 D= 2 0 0 0 -4 0 0 0 2T=inv（V）*A*V T= 2.0000 0 0 0.0000 -4.0000 0.0000 -0.0000 -0.0000 2.0000可以看出解向量的个数等于特征跟的重数。2） 设A是一个n阶实对称矩阵，那么存在一个n阶正交矩阵U，使得UTAU是对角形。 0 1 1 -1设A= 1 0 -1 1 ，求正交矩阵Q，使QTAQ是对角矩阵。 1 -1 0 1 -1 1 1 0解：求A的特征值。 -1 -1 1解特征值方程 |E-A| = -1 1 -1 =（+3）（-1）3 = 0 。 -1 1 -1 1 -1 -1 得到A的特征值是1=-3。2=3=4=1。求A的特征向量。 当1=-3时，特征向量方程为（1E-A）x=0，对特征矩阵（1E-A）作初等行变换： -3 -1 -1 1 1 1 1 1 -1 -3 1 -1 0 -2 2 01E-A= -1 1 -3 -1 0 2 -2 0 1 -1 -1 -3 0 -2 0 -2 1 1 1 1 1 0
