《矩阵分析与数值分析》实验报告

《《矩阵分析与数值分析》实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《矩阵分析与数值分析》实验报告院系：姓名学号：所在班号:任课老师:一.设错误！未找到引用源。，分别编制从小到大和从大到小的顺序程序计算错误！未找到引用源。并指出有效位数。程序如下：functionsum3j二input（'请输入求和个数"j":;A=0;B二0;doubleB;doubleA;forn=2:jm=rT2T;t=l./m;A二A+t;enddispC从小到大:'）s二Aforn=j：-l:2m=n2-1;t=l./m;B二B+t;enddispC从大到小:’）s二B运行结果：>>sum3请输入求和个数〃j〃：100从小到大：s=

2、0.740049504950495从大到小：s=0.740049504950495>>sum3请输入求和个数V：10000从小到大：s=0.749900004999506从大到小：s=0.749900004999500>>sum3请输入求和个数〃j〃：1000000从小到大：s=0.749999000000522从大到小：s=0.749999000000500二、解线性方程组1.分别Jacobi迭代法和Gauss-Seide1迭代法求解线性方程组。<-2100>(一1)1-210001-21兀30001-2xA0丿/47/迭代法计算停止的

3、条件为：max丫仏+i)<10"Ol

4、isp('x=');disp(x(:,k+1));ifnorm(x(:,k+l)-x(:,k))

5、

6、x(k+1)-x(k)

7、

8、

9、

10、x"(k+l)-x"(k)

11、

12、

13、00.799998799067310.599998427958700.399998056850090.19999902842505(2)Gauss-Seidel迭代法程序代码:functiongauss_seidcl(A,b,N)clc;clear;A=[-2100;l-210;01-2l;001-2];b=[-l000]';N二100;n二size(A,1);D=diag(diag(A));L=tril(-A,-1);U=triu(-A,1);Tg=inv(D-L)*U;eg=inv(D-L)*b;tol=le-06;k=1;x=zeros

14、(n,1);whilek<=Nx(:»k+1)=Tg*x(:,k)+eg;disp(k);disp(，x=');disp(x(:,k+1));ifnorm(x(:,k+1)-x(:,k))〈toldisp(5TheprocedurewassuccessfuT)disp(JCondition

15、

16、(k+1)-x"(k)

17、

18、

19、ion

20、(k+1)-x"(k)

21、

22、(l)Gauss歹ij主元消去法程序代码：functionx=Gaussmain(A,b)clc;clear;format1ongA二[1212;253-2;-2一235;1323];b二[47-10

23、]';N=length(A);x=zeros(N,1);y=zeros(N,1);c二0;d二0;A(:,N+l)=b;fork=l:N-lfori=k:4ifc<