特征值与特征向量.doc_第1页
特征值与特征向量.doc_第2页
特征值与特征向量.doc_第3页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

特征值与特征向量一、基本理论设A是n阶方阵. 如果存在非零向量x使Ax=x对某个常数C成立, 则称是A的特征值 , x是属于特征值的特征向量。()=|IA|是n次多项式, 称为A的特征多项式。特征多项式的n个根就是A的全体特征值。特征值的性质: n个特征值之和等于A的迹,即对角元之和;n个特征值之积等于A的行列式。上三角矩阵的对角元就是全体特征值. 实对称方阵A的全体特征值为实数,且有n个线性无关的特征向量(x1,xn),将其作为列向量排成矩阵P,则AP=PD, 其中D=Diag(1, n)为对角阵,对角元为A的特征值.P-1AP=D, 实对称矩阵A相似到对角阵D,称对角化.求矩阵n阶方阵A的特征值与特征向量的步骤:1. 求出A的特征多项式2. 求出特征多项式的n个根1, n3. 对每个特征值i,解齐次线性方程组 (iI-A)x=0,解空间的任意一组基即是A的属于特征值i的特征向量.二、Matlab实现实现一:1. 用Matlab命令 p=poly(A) 得到A的特征多项式的系数p(p是n+1维向量) 2. 用命令 roots(p) 特征多项式的根: roots(p)3. 对每个特征值i ,解齐次线性方程组 (iI-A)x=0,可利用Matlab提供的null(iI-A),或先前编写的函数nulbasis(iI-A) (见齐次线性方程一节内容)实现二: 直接使用 Matlab命令 P D = eig(A), 同时返回特征向量组 P 和全体特征值 D.三、例子例. 求的特征值与特征向量.A=3, 2, 1; 0, -1, -2; 0, 0, 3 A = 3 2 1 0 -1 -2 0 0 3 特征多项式p=poly(A) p = 1 -5 3 9 roots(p) ans = 3.0000 3.0000 -1.0000 特征值为3,-1.X1=nulbasis(3*eye(3)-A) X1 = 1.0000 0 0 -0.5000 0 1.0000 属于特征值3的特征向量:,X2=nulbasis(-eye(3)-A) X2 = -0.5000 1.0000 0 属于特征值-1的特征向量:P=X1 X2 P = 1.0000 0 -0.5000 0 -0.5000 1.0000 0 1.0000 0 inv(P)*A*P ans = 3 0 0 0 3 0 0 0 -1 解二:Matlab的eig(A)P D=eig(A) P = 1.0000 -0.4472 0 0 0.8944 -0.4472 0 0 0.8944D = 3 0 0 0 -1 0 0 0 3 inv(P)*A*P ans

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论