吴文刚的论文11.docVIP

线性方程组迭代方法的总结及程序实现 线性方程组线性方程组迭代方法程序实现迭代方法 (1) 将式(1)用矩阵表示为Ax=b,并记 , 于是A=D-E-F,按前面假定有det(D)≠0.这样, 公式可表示为 x(k+1)=D-1(E+F)x(k)+D-1b, (2) x(k+1)=(D-E)-1Fx(k)+(D-E)-1b, (3) x(k+1)=(D-ωE)-1[((1-ω)D+ωF)x(k)+ωb] (4) 显然,公式(2), (3), (4)具有统一的形式 x(k+1)=Bx(k)+f (5) 2.1 设x为方程(1)的解,即有x=Bx+f,将式(5)与之相 减有,递推得到x(k+1)-x=Bk(x(0)-x), 两边取极限可知 2.2 设λ为B的特征值,ξ为相应的特征向量,则Bkξ= λkξ,可见 由此可得迭代收敛定理: 2.3 此程序在Matlab中的实现 L=zeros(n);U=zeros(n); for i=1: n 　　for j=i+1: n 　　L(,j i)=B(,j i);U(,i i)=B(,i i);U(