用matlab实现阵的对用matlab实现阵的对角.doc

用matlab实现矩阵的对角 对角阵在实际上的应用特别广泛，对角阵解决现实问题上很方便，通过对角矩阵可以最简单地处理物力问题，也可以解出线性方程组的解；最普遍的是可以直接知道相似矩阵的行列式值，秩，特征值等，所以可以说研究对角化问题是特别重要。对角化的最快，最方便的方法是利用matlab软件。 1 2 3 形式的矩阵为对角矩阵 一般 n （空白处为零）。 在相似变换下，方阵A的许多重要性质（如 行列式，秩，特征值等）保持不便，因此我们可以通过相似变换将矩阵A化简，并利用化简后的矩阵来研究与矩阵有关的问题。 我们讨论的主要问题是：对于n阶方阵A，是否存在可逆方阵P，使P-1AP为对角矩阵。这就是矩阵的对角化问题。 2 0 0 -1 1 0 通过相似变换下简化为对角矩阵B= 0 1 0 -4 3 0 A= 矩阵 1 0 0 1 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对角 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.3927 D= -3.7266 0 0 0 0 0.9416 0 0 0 0 1.9420 0 0 0 0 4.8430 程序说明： D 对角线上的元素为A的特征值，V为相对应的特征向量所构成的矩阵。可以看出特征值都是实数，每个特征值对应于一个特征向量，对应的特征向量是正交的，特征向量所构成的矩阵是正交阵。如果不用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； 结果一样，A有3个互不相同的特征值，故可对角化；这两种方法相互比较，用Matlab程序好计算这种问题。 3）求A的特征向量，若有n个线性无关的特征向量，则A可以对角化，否则不可以对角化。 0 1 1 -1 -1 A= 1 0 1 因为 |I-A| = -1 -1 =（-2）（+1）2 1 1 0 -1 -1 于是