高等数值分析第二次实验作业T1. 构造例子特征值全部在右半平面时, 观察基本的Arnoldi 方法和GMRES 方法的数值性态, 和相应重新启动算法的收敛性.Answer:(1) 构造特征值均在右半平面的矩阵A:根据实Schur分解,构造对角矩阵D由n个块形成,每个对角块具有如下形式,对应一对特征值这样D=diag(S1,S2,S3Sn)矩阵的特征值均分布在右半平面。生成矩阵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;(2) 观察