第四章特征值与特征向量工程技术中的振动问题与稳定性问题,数学中矩阵的对角化与微分方程组的求解问题,还有其他一些实际问题,都可以归结为求矩阵的特征值与特征向量本章先介绍矩阵的特征值与特征向量的概念,再引入相似矩阵的概念,并讨论矩阵的相似对角化,最后介绍几维向量空间的正交性及实对称矩阵的相似对角化。4.0引例R第三章提到的高光谱图像由搭载在不同空间平台上的成像光谱仪,以数十至数百个连续且细分的光谱波段对目标区域同时成像得到,其中每个波段对应一幅灰度图像,也称为高光谱图像的一个通道,在处理高光谱图像时,同时挖拥多个通道中的信息比单独前沿视角处理每个通道更有效。通常,成像区域的一些特征可同时出现在多特征值与特征尚量个通道中,这意味着高光谱图像中包含了大量几余信息主成分分析方法可将原始多个通道的高光谱图像中的大部分信息用少数几幅彼此不相关的合成图像来表示.比如,图4.1(a),(b),(c)分别表示高光谱图像中的三个光谱波(a)波段7(b)波段30(c)波段120(d)第一主成分(e)第二主成分()第三主成分图4.1
第四章特征值与特征向量142段图像,主成分分析方法将这三个波段的图像重新线性组合,得到了三幅主成分图像,如图4.1(d),(e),(f)所示,第一主成分图像(d)的方差占比率为97.08%.这意味着,第一主成分图像包含了原始三个波段图像97.08%的信息,如此,信息得以浓缩和简化.那么如何得到这些主成分图像呢?这个间题的答案与本章内容密切相关,特征值与特征向量的概念与计算4.1在实际问题中,常常碰到这样的问题,即对于一个给定的n阶方阵A,是否存在非零的n维向量α,使得Aα与α平行,即存在常数入,使得Aα=入α成立,在数学上,这就是特征值与特征向量的问题定义设A是n阶方阵,如果存在数入和n维非零向量Q,使(4.1)Aa=a,则称入为方阵A的一个特征值,α为方阵A对应于特征值入的一个特征向量典型例随精讲例1设开征售包里的极有Ao202AaAB=7由定义可知,1与2就是A的两个特征值,α1与α2就是A分别对应于特征值1与2的特征向量;而β则不是A的特征向量从几何上看,矩阵A分别乘向量α1α2与β的结果如图4.2所示,Aα2相当于将向量α2增大一倍.这说明,如果α是A的特征向量,那么Aα相当于对α作一次“伸缩变换VNA
4.1特征值与特征向量的概念与计算143X2Aa.AaB-C1xAp图4.2例2设方阵A满足A2=A.试证A的特征值只有0或1.证设入是A的特征值,α是A对应于入的特征向量,则e典型例题精讲A=Qα0).于是特征值、特征α=Aα=AQ=A(Aα)=A(Aα)=(Aα)=Q,向量的性质所以(12-)Q= 0.因为α0所以2=(1)=0,即=0或=1.对于n阶方阵A,如果α是A对应于特征值入的一个特征向量,那么对任意的k≠0,A(ka)=k(Aα)=k(a)=入(ka)所以,kα也是A对应于特征值入的特征向量设i,α2,.,α都是A对应于特征值入的特征向量,且kiQ1+k2α2++k,α+0,则A(kiai+k2a2+...+krar)=k(Aai)+k2(AQ2)+.**+k(Aar)= k (Aai)+ k2 (Aa2)+...+ kr(Aar)=入(kia1+k22+...+kar)所以kiα1+k2Q2++hα也是A对应于特征值入的特征向量设V是n阶方阵A对应于特征值入的所有特征向量以及零向量所组成的集合,即V=(|Aa=AQ,AeC,aEC".由以上分析可知V对向量的加法和数乘封闭,故V构成子空间.我们称V为A的特征子空间在例1中,Aα1=α1,AQ2=2α2,不难证明A对应于特征值入1=1的所有特征向量都可以由α1线性表出,对应于入2=2的所有特征向量都可以由α2线性表出从几何上看,V,对应于过原点与点(1,1)的直线1,V,对应于过原点与点(2,1)的直线12(图4.3),即V,=11上的所有向量},V=[12上的所有向量
144第四章特征值与特征向量X2AaoX图4.3下面讨论对于给定的方阵A=(a)nxn,怎样求A的特征值与特征向量设是方阵A对应于特征值入的特征向量,即Aα=入α(α半O)(I-A)α=0于是,α是齐次线性方程组(入I-A)X=0即(-au)a1-a12r2-...-ainan=0,-a21t1+(-a22)a2-.-a2mn=0.an121-an222-+(-anm)en=0的非零解,因此,(ΛI-A)X=0的解空间就是A的特征子空间,它的维数为dimV=n-R(I-A)(I-A)X=0的基础解系为 V的基,因为齐次线性方程组有非零解的充要条件是系数行列式等于零,所以det (I - A) = 0.我们将方程det(入I-A)=0称为方阵A的特征方程,特征方程的根就是特征值,故有时又将特征值称为特征根,若入是单根,则称入为A的单特征根:若入是重根,则称入为A的重特征根由以上分析可得求方阵A的特征值与特征向量的计算步骤如下:1°求特征方程det(入I-A)=0的全部相异根入1,2,入(k≤n):分别求(入,I-A)X=0(i=1.2..,k)的基础解系αi,Qi2,,αir,则20kiai+k2ai2+...+krQi(ki,k2,".,hr,不全为零)就是A对应于特征值入,的全部特征向量3的特征值与特征向量例3求方阵A012>-32-3入+2=(Λ-1)(入-2),A的特征值为解det(I-A)-1入=12=2
一4.1特征值与特征向量的概念与计算145对于入1=1,齐次线性方程组(iI-A)X=0的系数矩阵为相应简化的齐次线性方程组为21-22=0,其基础解系为α1=(1,1),故对应于>1=1的A的全部特征向量为hα1(k1≠0).对于>2=2,齐次线性方程组(>2I-A)X=0的系数矩阵为相应简化的齐次线性方程组为21-22=0,其基础解系为α2=(2,1),故对应于>2=2的A的全部特征向量为k2α2(k2≠0).2-21例4求A=-224的特征值与特征向量,2-24解入-12-2-12Y-2det (-A)=21+2-42入+2-4-2-4入+201-2入-2入14-2142入+6-4=(入-2)2入+600入2= (入 2)2(入 + 7),A的特征值为入1=2(2重),入2=-7.对于>1=2,齐次线性方程组(A1I-A)X=0的系数矩阵为21-2224002-44000相应简化的齐次线性方程组为31=-2r2+2r3,其基础解系为α1=-2,1,0),Q2=(2,0,1)T,A对应于入1=2的全部特征向量为kα+k2Q2(h1,kz不全为零)对于入2=-7,齐次线性方程组(>2I-A)X=0的系数矩阵为2N2-5相应简化的齐次线性方程组为2其基础解系为 Q3=(1,2,-2)},A对应To:13于入2=-7的全部特征向量为kQ3(k3≠0)