数值逼近上机练习答案..docxVIP

下载本文档

11
0
约2.45万字
约 32页
2017-01-13 发布于重庆
举报
版权申诉

数值逼近上机练习答案..docx

1、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

数值逼近上机练习答案.

实验一分别用循环型Arnoldi算法和循环型GMRES算法求解线性方程组。取，初始值，迭代终止条件为。要求输出数值近似解和迭代步数。M-文件Arnoldi算法function X=Arnoldi(X0,m,A,b,ep)nk=0;while(max(abs(A*X0-b))=ep)r0=b-A*X0;n=length(b);v=zeros(n,m+1);h=zeros(m,m);p=sqrt(sum(r0.^2));v(:,1)=r0./p;for k=2:m; tmp=0; for i=1:k-1 h(i,k-1)=(A*v(:,k-1))*v(:,i); tmp=tmp+h(i,k-1)*v(:,i); end vk=A*v(:,k-1)-tmp; h(k,k-1)=sqrt(sum(vk.^2)); v(:,k)=vk/h(k,k-1);endtmp=0;for i=1:m h(i,m)=(A*v(:,m))*v(:,i); tmp=tmp+h(i,m)*v(:,i);endvm1=A*v(:,m)-tmp;hm1=sqrt(sum(vm1.^2));v(:,m+1)=vm1/hm1;e1=zeros(m,1);e1(1)=1;y=inv(h)*p*e1;z=v(:,1:m)*y;X0=X0+z;nk=nk+1;endnkX=X0;Command windows：A=[1 0 0 0;1 1 0 0;1 1 1 0;1 1 1 1]b=[1 1 1 1]m=2x0=[0 0 0 0]ep=1.0e-8结果：Arnoldi(X0,m,A,b,ep)nk = 23ans = 1.0000 0.0000 -0.0000 -0.0000GMRES算法function X=GMRES(X0,m,A,b,ep)nk=0;while(max(abs(A*X0-b))=ep)r0=b-A*X0;n=length(b);v=zeros(n,m+1);h=zeros(m,m);p=sqrt(sum(r0.^2));v(:,1)=r0./p;for k=2:m; tmp=0; for i=1:k-1 h(i,k-1)=(A*v(:,k-1))*v(:,i); tmp=tmp+h(i,k-1)*v(:,i); end vk=A*v(:,k-1)-tmp; h(k,k-1)=sqrt(sum(vk.^2)); v(:,k)=vk/h(k,k-1);endtmp=0;for i=1:m h(i,m)=(A*v(:,m))*v(:,i); tmp=tmp+h(i,m)*v(:,i);endvm1=A*v(:,m)-tmp;hm1=sqrt(sum(vm1.^2));v(:,m+1)=vm1/hm1;H=[h;[zeros(1,m-1),hm1]];e1=zeros(m+1,1);e1(1)=1;y=pinv(H)*p*e1;z=v(:,1:m)*y;X0=X0+z;nk=nk+1;endnkX=X0;Command windowsA=[1 0 0 0;1 1 0 0;1 1 1 0;1 1 1 1]b=[1 1 1 1]m=2x0=[0 0 0 0]ep=1.0e-8结果：GMRES(X0,m,A,b,ep)nk = 12ans = 1.0000 0.0000 -0.0000 -0.0000用追赶法、线性插值法和双参数法求解阶三对角方程组。其中，要求分别输出当时的计算误差和计算时间。M文件双参数法function [X,deta]=sdjscs(a,b,c,n,f)if(length(f)=n length(f)=n) s=zeros(n,1); t=zeros(n,1); tic s(1)=0; s(2)=f(1)/c(1); t(1)=1; t(2)=-b(1)/c(1); for i=3:n s(i)=(f(i-1)-b(i-1)*s(i-1)-a(i-1)*s(i-2))/c(i-1); t(i)=-(b(i-1)*t(i-1)+a(i-1)*t(i-2))/c(i-1); end x=zeros(n,1); x(1)=(f(n)-a(n)*s(n-1)-b(n)*s(n))/(a(n)*t(n-1)+b(n)*t(n)); for i=2:n x(i)=t(i)*x(1)+s(i); end X=x; A=[zeros(1,n);[diag(a(2:end)) zeros(n-1,1)]]+diag(b)+[[zeros(n-1,1) diag(c(1:end-1))];zeros(1,n)]; deta=norm(f-A*X,2); tocelse X=you are wrong!;end 窗口： n=1000;a