一、冪法（乘冪法）與反冪法第八章：特征值問題的計算方法1、冪法求按模最大的特征值和特征向量clear，clcA=[210;131;014];k=0;u0=[111]';u=A*u0;[t,i]=max(abs(u));mu=u(i);u=u/mu;whilenorm(u-u0)>10^(-7)&k<10^4u0=u;k=k+1;u=A*u0;[t,i]=max(abs(u));mu=u(i);u=u/mu;endu,mu,k2、冪法的原點平移：加速收斂使得：A=[210;131;014];d=2.14;B=A-d*eye(3);k=0;u0=[111]’;u=B*u0;[M,i]=max(abs(u));mu=u(i);u=u/mu;whilenorm(u-u0,inf)>10^(-7)&k<10^4k=k+1;u0=u;u=B*u;[M,i]=max(abs(u));mu=u(i);u=u/mu;endmu=mu+d,u,kclear,clcA=[210;131;014];n=length(A);[L,U,P]=lu(A);u0=ones(n,1);

u=P*u0;

u=U\(L\u);[m,i]=max(abs(u));mu=u(i);u=u/mu;k=0;whilenorm(u-u0)>=10^(-15)&k<10^4;k=k+1;u0=u;u0=P*u0;u=U\(L\u0);[m,i]=max(abs(u));mu=u(i);u=u/mu;endu,mu=1/mu,k3

反冪法：求按模最小的特征值clear,clc,formatlongA=[210;131;014];r=1.26795;n=length(A);B=A-r*eye(n);[L,U,P]=lu(B);u0=ones(n,1);u=P*u0;u=U\(L\u);[m,i]=max(abs(u));

mu=u(i);u=u/mu;k=0;whilenorm(u-u0)>=10^(-12)&k<10^4;k=k+1;u0=u;u=P*u0;u=U\(L\u);[m,i]=max(abs(u));mu=u(i);u=u/mu;endu,mu=1/mu+r,k4

帶位移反冪法：求特征向量（高精度）Givens變換，又稱平面旋轉變換若只需將向量的某一個分量化為0時，采用Givens變換。稱下列矩陣為Givens變換矩陣：行行列列n=3時記取適當的可使得：令推廣到一般情形：令則有function[c,s]=giv(a,b)ifb==0c=1;s=0;elseifabs(b)>abs(a)tau=a/b;s=1/sqrt(1+tau^2);c=s*tau;elsetau=b/a;c=1/sqrt(1+tau^2);s=c*tau;end二、Jacobi迭代法求對稱矩陣的全部特征值和特征向量1、Givens（旋轉）變換2、用Givens變換把向量的某些分量化為零clear,clcn=5;a=rand(5,1),forp=2:n[c,s]=giv(a(1),a(p));r1=c*a(1)+s*a(p);r2=-s*a(1)+c*a(p);a(1)=r1;a(p)=r2;enda93、用Givens變換把矩陣的某些元素化為零clear,clcA=[3,2,2;4,1,2;0,2,1],n=3;r=zeros(1,n);t=zeros(1,n);forp=1:2[c,s]=giv(A(p,p),A(p+1,p));forj=1:n

r(j)=c*A(p,j)+s*A(p+1,j);t(j)=-s*A(p,j)+c*A(p+1,j);endA(p,:)=r;A(p+1,:)=t;endA10作業1.編程解8.8用帶原點