实验5矩阵的特征值特征向量相似对角形.docVIP

实验5 矩阵的特征值特征向量Matlab命令求方阵的特征值和特征向量以及矩阵的相似对角化，并解决相关问题. 实验原理 （一）预备知识? 线性代数中的、知识； ? 本实验所用Matlab命令提示：eig(A)返回A的特征值列向量，如果A是实对称的，特征值为实数．特征值也可能为复数. 2. 命令[P,X]=eig(A)输出的是对角线上的元素A的特征值对角矩阵也可以先求出方阵的特征多项式，然后根据特征多项式确定特征值（）例求3阶方阵A的特征clear (清理屏幕) A=[1 2 3;2 1 3;3 3 6]; v=eig(A) [P,X]=eig(A) 输出为： v= -1.0000 -0.0000 9.0000 P= 0.5774 0.4082 -0.7071 0.5774 0.4082 0 -0.5774 0.8165 X= -1.0000 0 0 0 -0.0000 0 0 0 9.0000 其中v为特征向量，P为以特征向量为列的矩阵，X为以特征值为对角线元素的对角矩阵. 当然，我们也可以先求出方阵的特征多项式，然后根据特征多项式确定特征值求3阶方阵 的特征A=[3 2 4;2 0 2;4 3 2]; PA=poly(A) %A的特征多项式x =roots(PA) 矩阵的相似对角化 若n阶方阵A有n个线性无关的特征向量，则A与对角形矩阵相似，即存在 可逆矩阵P，使得为对角矩阵.实对称矩阵一定与对角矩阵相似，而且存在正交矩阵P，使得为对角矩阵. 例3 设方阵,求一个可逆矩阵,使得为对角矩阵. 解 (用命令[P,X]=eig(A)) 输入： clear; A=[4,1,1;2,2,2;2,2,2]; A=sym(A); [P,X]=eig(A) (输出的特征向量没有对角化) 输出为： P= [1,0,-1] [1,-1,1] [1,1,1] X= [6,0,0] [0,0,0] [0,0,2] 从输出结果可知，A的特征值6，0，2，特征向量是，且线性无关，所以根据线性代数中相似对角化的知识，可得就是要求的可逆矩阵，且 . 例4 方阵是否与对角矩阵相似? 解 只要检查A的线性无关的特征向量的数目是否为2即可.所以 输入: clear; A=[1,0;2,1]; [P,X]=eig(A) 输出为: P= 0 0.0000 1.0000 -1.0000 X= 1 0 0 1 输出结果显示,1是A的二重特征值,但是两个特征向量线性相关,因此矩阵A不与对角矩阵相似. 例5 对于实对称矩阵,求一个正交矩阵,使得为对角矩阵. 解 输入: clear; A=[0,1,1,0;1,0,1,0;1,1,0,0;0,0,0,2]; [P,X]=eig(A) 输出为: P= 0 0.5774 -0.7152 0.3938 0 0.5774 0.0166 -0.8163 0 0.5774 0.6987 0.4225 1.0000 0 0 0 X= 2.0000 0 0 0 0 2.0000 0 0 -1.0000 0 0 0 0 -1.0000 从显示结果可知,P就是所求的正交矩阵. 如果要验证P是正交矩阵,而且是对角矩阵,输入: transpose(P)*P x1=inv(P)*A*P x2= transpose(P)*A*P 输出为： 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000 x1= 2.0000 0 0 0 0 2.0000 -0.0000 0.0000 0 -0.0000 -1.0000 -0.0000 0 -0.0000 0.0000 -1.0000 x2=