值分析第四版 颜庆津 计算实习题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