数值计算方法第四章 特征值与特征向量.pptVIP

* ? 工程中有许多振动问题，如桥梁和建筑的振动，机械机件的振动飞机积机翼的颤动等。这些问题的求解均归结为求的矩阵特征值问题。一些稳定性分析问题及相关问题可以转化为求矩阵特征值与特征向量的问题。 第四章 矩阵特征值与特征向量的计算 主要内容：1、幂法和反幂法 2、JACOBI方法 3、Q-R算法 定义：按下（4-1，4-4）计算矩阵取模最小的特征值与相应的特征向量的方法称为反幂法。 1、???幂法和反幂法 定义：按下（4-1，4-4）计算矩阵A取模最大的特征值与相应的特征向量的方法称为幂法。 方法： 其中 为 相应的特征向量的近似。 ? 由最后可得特征值表示 注：在计算过程中为防止溢出，需归一化。 幂法实际使用公式为： 例1 用幂法求矩阵 的取模最大的特征值与相应的特征向量。 取x(0)=(0,0,1)T 取由求误差不超过10-3 解：由算法知， 事实上A的特征值为 与相应的 特征向量为（1，-1，1）T。 继续有k=0，1，2，---，9。 X1（k）=0,0,0.5,---, 2.8436517 . =1,2,2.5,---,2.9996973. X2（k）=0,-1,-2---,-2.9993946 . X3（k）=1,2,2.5,---,2.9996973. =1,2,2.5,---,2.9996973. y1（k）=0,0,0.2,---, 0.9219772 . y2（k）=0,-0.5,-0.8,---,-0.9996973 . y3（k）=1,1,1,---,1. 因2.9996973-2.999092410-3 所以 ,相应的特征向量 因2.9996973-2.999092410-3 所以 ,相应的特征向量 A=[2 -1 0;0 2 -1;0 -1 2]; u0=[0;0;1]; u0max=max(u0);u0min=min(u0); if abs(u0max)=abs(u0min) a=u0max; else a=u0min; end y0=sign(a)*u0/abs(a); u1=A*y0; for i=1:3 if(a==u0(i)) r=i; end end k=0; t1=sign(a)*u1(r); lmda=t1; fprintf( k u(k)_1 u(k)_2 u(k)_3 lamda\n) fprintf(%2d %12.9f %12.9f %12.9f %12.9f\n,k,y0,lmda) for k=1:10 u0=u1; u0max=max(u0);u0min=min(u0); if abs(u0max)=abs(u0min) a=u0max; else a=u0min; end y0=sign(a)*u0/abs(a); u1=A*y0; for i=1:3 if(a==u0(i)) r=i; end end u1=A*y0; t1=sign(a)*u0(r); lmda=t1; x=y0; fprintf(%2d %12.9f %12.9f %12.9f %12.9f\n,k,y0,lmda) end k u(k)_1 u(k)_2 u(k)_3 lamda 0 0.000000000 0.000000000 1.000000000 2.000000000 1 0.000000000 -0.500000000 1.000000000 2.000000000 2 0.200000000 -0.800000000 1.000000000 2.500000000 3 0.4285714