用matlab实现矩阵对角.pdfVIP

用 matlab 实现矩阵的对角 对角阵在实际上的应用特别广泛， 对角阵解决现实问题上很方便， 通过对角 矩阵可以最简单地处理物力问题， 也可以解出线性方程组的解； 最普遍的是可以 直接知道相似矩阵的行列式值， 秩，特征值等， 所以可以说研究对角化问题是特 别重要。对角化的最快，最方便的方法是利用 matlab 软件。 1 2 3 一般 形式的矩阵为对角矩阵 n （空白处为零）。 在相似变换下， 方阵 A 的许多重要性质 （如 行列式，秩，特征值等）保持不便， 因此我们可以通过相似变换将矩阵 A 化简，并利用化简后的矩阵来研究与矩阵有 关的问题。 -1 我们讨论的主要问题是：对于 n 阶方阵 A ，是否存在可逆方阵 P，使 P AP 为 对角矩阵。这就是矩阵的对角化问题。 -1 1 0 2 0 0 矩阵 A = -4 3 0 通过相似变换下简化为对角矩阵 B= 0 1 0 1 0 2 0 0 1 通过矩阵 B 可以直接知道矩阵的行列式，秩，特征值，对应方程组的解等 重要性质。 但任何矩阵不一定可以对角化。一个矩阵是否可以对角化有如下的判断方法： 1）判断 A 是否实对称矩阵，茹是一定可对角化，因为 A 是实对称矩阵，则 有 （1）A 的全部特征值是实数。 （2 ）A 的不同特征值对应的特征向量正交。 （3 ） -1 T A 一定相似于对角矩阵，且存在正交矩阵 T ，使得 T AT=T AT= 。 的对角 元素是 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