数值分析 反幂法求矩阵按模最小特征值 MATLAB实现

值分析第四版 颜庆津 计算实习题P238

结果：

%function [lam]=jingmi(ep)

%反幂法求矩阵按模最小特征值

%2015.11.8 密密编写  (*^__^*) 嘻嘻……

 

 

function [lam]=jingmi(ep)%lam是A的按模最小的特征值

%初始化

   n=501;

   a=1:n;%存储A的主对角线元素

   for i=1:n

       a(i)=(1.64-0.024*i)*sin(0.2*i)-0.64*exp(0.1/i);

   end

   b=0.16;

   c=-0.064;

   u0=ones(n,1);%迭代初始向量

   u1=ones(n,1);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%第一次迭代   

        m=sqrt(u0'*u0);

   y=u0/m;

%以下代码功能：三角分解法，从线性方程组A*u1=y求解u1

b0=y;

   L=zeros(3,n);

   U=zeros(3,n);

   

   %L中的主对角线元素均为1

   for i=1:n

        L(1,i)=1;

   end

%%%%%%%%%%%%%%%%%%%%%%doolittle分解,第一部分A=LU

   for k=1:n

      

       for j=k:min(k+2,n) %计算 U(k,j)

           sum=0; %临时和

           M=[1,k-2,j-2];%存储t的初始量备选值，t取M中的最大值

           for t=max(M):k-1

                sum=sum+L(k-t+1,t)*U(t-j+3,j);

           end