清华大学贾仲孝老师高等数值分析第二次实验(共19页).docx

上传人:晟*** 文档编号:6926628 上传时间:2021-09-15 格式:DOCX 页数:20 大小:386.76KB
下载 相关 举报
清华大学贾仲孝老师高等数值分析第二次实验(共19页).docx_第1页
第1页 / 共20页
清华大学贾仲孝老师高等数值分析第二次实验(共19页).docx_第2页
第2页 / 共20页
清华大学贾仲孝老师高等数值分析第二次实验(共19页).docx_第3页
第3页 / 共20页
清华大学贾仲孝老师高等数值分析第二次实验(共19页).docx_第4页
第4页 / 共20页
清华大学贾仲孝老师高等数值分析第二次实验(共19页).docx_第5页
第5页 / 共20页
点击查看更多>>
资源描述

高等数值分析第二次实验作业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) 观察

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 实用文档资料库 > 公文范文

Copyright © 2018-2021 Wenke99.com All rights reserved

工信部备案号浙ICP备20026746号-2  

公安局备案号:浙公网安备33038302330469号

本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。