数值分析实验报告二.doc

下载文档

10
0
约 8页
2017-10-07 发布于重庆
举报
版权申诉
保障服务

数值分析实验报告二.doc

1、本文档共8页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

数值分析实验报告二

关于sor算法的一些实验一程序源代码（使用二范数控制精度） function [x,NN]=sor(a,b,w,esp,N) if nargin4 esp=10e-10;N=200; end if nargin==4 N=200; end NN=0; n=size(a,1);s1=0; x1=ones(1,n);x2=zeros(1,n); while(((sum((x1-x2).^2)).^0.5=esp)(NNN)) x1=x2; NN=NN+1; for i=1:n if i==1 s2=0; for j=1:n s2=s2+a(i,j)*x1(j); end x2(i)=x1(i)+w*(b(i)-s2)/a(i,i); elseif i~=1 s1=0; for j=1:(i-1) s1=s1+a(i,j)*x2(j); end s2=0; for j=i:n s2=s2+a(i,j)*x1(j); end x2(i)=x1(i)+w*(b(i)-s1-s2)/a(i,i); end end end x=x2; 二检验书上结论 1 迭代的速度与松弛因子的选择有关 a=[-4 1 1 1;1 -4 1 1;1 1 -4 1;1 1 1 -4]; b=[1 1 1 1]; for i=1:10 w=0.9+0.1*i; [x,NN]=sor(a,b,w,1e-5); vpa(x,10),NN,w end 运行结果： ans = [ -.9999940999, -.9999948735, -.9999955457, -.9999961298] NN = 22 w = 1 ans = [ -.9999946977, -.9999955972, -.9999963441, -.9999969643] NN = 17 w = 1.1000 ans = [ -.9999984706, -.9999988509, -.9999991309, -.9999993323] NN = 13 w = 1.2000 ans = [ -1.000001519, -.9999992183, -1.000000116, -1.000000520] NN = 12 w = 1.3000 ans = [ -.9999982042, -1.000000770, -1.000000348, -.9999989677] NN = 15 w = 1.4000 ans = [ -.9999978969, -.9999995670, -1.000001281, -.9999982763] NN = 19 w = 1.5000 ans = [ -1.000000916, -.9999976894, -1.000002316, -.9999988757] NN = 25 w = 1.6000 ans = [ -.9999981218, -1.000000921, -1.000001526, -.9999980029] NN = 36 w = 1.7000 ans = [ -1.000002795, -1.000000621, -.9999977532, -1.000002784] NN = 56 w = 1.8000 ans = [ -1.000001441, -.9999962143, -1.000001956, -.9999999063] NN = 118 w = 1.9000 可见w=1.3时迭代次数最少为12次。但是结果与书上给出的略有不同。可能是程序的细节编写不同所导致。 2 当w=2的情况 w=2.0;[x,NN]=sor(a,b,w,1e-5);NN,vpa(x,10) NN = 200 ans = [ -1.57384563