矩阵的特征值和征向量、二次型.docVIP

数学工具软件教案 授课题目：矩阵的特征值和特征向量、二次型 授课时间：2012年2月28日、3月1日 教学目的与要求：学会用MATLAB软件求矩阵的特征值和特征向量；学会用MATLAB软件将二次型化为标准型；学会用MATLAB软件编程来判断二次型的正定性 教学重点与难点：用MATLAB软件求矩阵的特征值和特征向量；用MATLAB软件编程来判断二次型的正定性 教学方法：讲授法 教学手段：多媒体教学 讲授内容： 一、特征值与特征向量 矩阵与向量相乘，即表示矩阵对向量的变换(Transformation)。一般说来，向量在变换的作用下将发生旋转(Rotation)、反射(Reflection)和放大缩小。但对于任何一个矩阵来说，总存在那么一些特殊的向量，在对其变换的作用下，向量的方向不变，而仅长短发生变化。这种向量就是所谓的特征向量。 定义：设是阶方阵，是一个数。如果存在非零的列向量，使得 成立，则称数为方阵的特征值(Eigenvalue)，非零列向量称为方阵的属于特征值的特征向量(Eigenvector)，该方程称为特征方程(Eigenvalue Equation)。 的全体特征值的和称为矩阵的迹(Trace)。它等于的主对角元素的和。 用Matlab计算特征值和特征向量的命令如下： d=eig(A) 仅计算A的特征值(以向量形式d存放) [V,D]=eig(A) 其中：D为由特征值构成的对角阵，V为由特征向量作为列向量构成的矩阵。且使 AV=VD 成立 trace(A) 计算矩阵A的迹 例1：求方阵的特征值、特征向量和迹 解： A=[2 2 -2;2 5 -4;-2 -4 5]; [V D]=eig(A) trace(A) V = -0.2981 0.8944 0.3333 -0.5963 -0.4472 0.6667 -0.7454 0 -0.6667 D = 1.0000 0 0 0 1.0000 0 0 0 10.0000 trace(A) ans = 12 答：特征值为：（二重），。对应于特征值的全部特征向量为：其中不能同时为零。对应于特征值的全部特征向量为：其中不能为零。矩阵的迹为： 二、矩阵的相似对角化 设，都是阶方阵，若存在 阶可逆矩阵，使：，则称矩阵，是相似的。 设是阶方阵，若与对角矩阵相似，则称可对角化。 定理1：阶方阵可对角化的充分必要条件是有个线性无关的特征向量。 例2：判断下列方阵是否可对角化。若可对角化，求出可逆阵P，使为对角阵。 ； 解(1)： A=[4 6 0;-3 -5 0;-3 -6 1]; [V D]=eig(A) V = 0 0.5774 -0.8944 0 -0.5774 0.4472 1.0000 -0.5774 0 D = 1 0 0 0 -2 0 0 0 1 rank(V) ans = 3 答：A可对角化，且 (2)： A=[0 1 0;-1 2 0;-1 1 1]; [V D]=eig(A) V = 0 0.6325 0.4511 0 0.6325 0.4511 1.0000 0.4472 0.7701 D = 1 0 0 0 1 0 0 0 1 rank(V) ans = 2 答：A不可对角化。 定理2：方阵可对角化的充分必要条件是它的几何重数等于代数重数。的特征值的几何重数为方程组的解空间的维数；的特征值的代数重数为作为特征根的重数。 下述函数可用来判断矩阵是否可对角化，若可对角化返回1，否则返回0。 function y=trigle(A) %可对角化返回1，否则返回0。 y=1;c=size(A); if c(1)~=c(2) y=0; return; end e=eig(A);n=length(A); while 1 if isempty(e) return; end d=e(1); f=sum(abs(e-d)10*eps); g=n-rank(A-d*eye(n)); if f~=g y=0; return;