lanczos算法求特征值matlab程序

《lanczos算法求特征值matlab程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、利用Lanczos算法化对称阵为三对角阵T,并利用Strum算法对T求解最大特征值和最小特征值.无中断情况:functionlanczos(A)[m,n]=size(A);q1=ones(m,1)/norm(ones(m,1),2);a1=q1'*A*q1;r1=A*q1-a1.*q1;b1=norm(r1,2);q2=r1./b1;a(1)=a1;b(1)=b1;fori=2:ma2=q2'*A*q2;a(i)=a2;r1=A*q2-a2.*q2-b1*q1;ifr1==0disp('ÖжÏÁË!');break;end

2、b1=norm(r1,2);q1=q2;q2=r1./b1;ifia1a=a1;endifb1>bb=b1;endendmax=0;min=0;a1=a;b1=b;c=(a1+b1)/2;whileabs(b

3、1-a1)>0.0005ifSs(A,c,n)>0a1=c;elseb1=c;endc=(a1+b1)/2;endmax=ccleara1b1ca1=a;b1=b;c=(a1+b1)/2;whileabs(b1-a1)>0.0005ifn>Ss(A,c,n)b1=c;elsea1=c;endc=(a1+b1)/2;endmin=c计算同号数:functions=Ss(A,c,n)p0=0;p(1)=1;p(2)=(A(1,1)-c)*p(1)-A(1,2)^2*p0;fori=3:n+1p(i)=(A(i-1,i-1)-c)

4、*p(i-1)-A(i-1,i-2)^2*p(i-2);ends=0;fori=n+1:-1:2ifp(i)==0s=s+1;elseifp(i)*p(i-1)>0s=s+1;endend例子:A=1513750800189263025170613>>lanczos(A);T=16.80005.15360005.15361.25577.26170007.2617-3.61473.58290003.58290.26570.98540000.98543.2933max=18.5853min=-10.0070