大規模矩阵问题的Krylov.pptVIP

大规模矩阵问题的Krylov 子空间方法综述;Krylov子空间方法综述;背景;投影方法;解方程组的投影法的矩阵表示 　　设n×m阶矩阵V=[v(1), v(2), …v(m)]与W=[w(1), w(2), …w(m)]的列分别构成K与L的一组基 .记z=x(1)-x(0),z=Vy,有 　 当　　　非奇异时（通常情况成立，见定理1.1），有 从而得到迭代公式 ;投影方法的最优性 １. (误差投影)设A为对称正定矩阵, x(0)为初始近似解,且K=L,则x(1)为采用投影方法得到的新近似解的充要条件是 其中, 且x为(1.1)的精确解.;;;Kyrlov子空间;Krylov子空间的标准正交基;;;;;Krylov子空间方法解线性方程组;解方程组的投影法的矩阵表示 　　设n×m阶矩阵V=[v(1), v(2), …v(m)]与W=[w(1), w(2), …w(m)]的列分别构成K与L的一组基 .记z=x(1)-x(0),z=Vy,有 　 当　　　非奇异时（通常情况成立，见定理1.1），有 从而得到迭代公式 ;;;;Krylov子??间方法解线性方程组;;;;;Krylov子空间方法解线性方程组;回顾;;;Krylov子空间方法解线性方程组;;;Krylov子空间方法解线性方程组;;回顾;;回顾;;;;Krylov子空间方法解线性方程组;;;;;Krylov子空间方法解矩阵特征值问题;;;;;;;;;改进技术;;;;;尚待解决的问题;CG法解方程组Ax=b;;;function x=cg(A,b,tol) t=cputime; n = numel(b); x0=rand(n,1) r=b-A*x0; s=b-A*x0; x=x0; j=1; while (norm(r)tol) count(j)=j-1; res(j)=norm(r); alpha(j)=dot(r,r)/dot(s,(A*r)); b1=dot(r,r); (续);x=x+(alpha(j))*s; r=r-(alpha(j))*(A*s); b2=dot(r,r); beta(j+1)=b2/b1; s=r+(beta(j+1))*s; j=j+1; end t=cputime-t plot(count,res),xlabel(迭代步),ylabel(残量),title(残量随迭代步的变化),set(findobj(gca,Type,line,Color,[0 0 1]),... Color,red,... LineWidth,2)