清华大学贾仲孝老师高等数值分析第二次实验.docxVIP

T1.构造例子特征值全部在右半平面时, 观察基本的Arnoldi 方法和GMRES 方法的数值性态, 和相应重新启动算法的收敛性.Answer:构造特征值均在右半平面的矩阵A：根据实Schur分解，构造对角矩阵D由n个块形成，每个对角块具有如下形式，对应一对特征值这样D=diag(S1,S2,S3……Sn)矩阵的特征值均分布在右半平面。生成矩阵A=UTAU，其中U为正交阵，则A矩阵的特征值也均在右半平面。不妨构造A如下所示：由于选择初值与右端项：x0=zeros(2*N,1);b=ones(2*N,1);则生成矩阵A的过程代码如下所示：N=500 %生成A为2N阶A=zeros(2*N);for a=1:N A(2*a-1,2*a-1)=a; A(2*a-1,2*a)=-a; A(2*a,2*a-1)=a; A(2*a,2*a)=a;endU = orth(rand(2*N,2*N));A1 = U*A*U;观察基本的Arnoldi和GMRES方法编写基本的Arnoldi函数与基本GMRES函数，具体代码见附录。function [x,rm,flag]=Arnoldi(A,b,x0,tol,m)function [x,rm,flag]=GMRES(A,b,x0,tol,m)输入:A为方程组系数矩阵，b为右端项，x0为初值，tol为停机准则，m为人为限制的最大