矩阵与数值分析.docVIP

2013级工科硕士研究生 《矩阵与数值分析》课程数值实验题目 一、设，分别编制从小到大和从大到小的顺序程序分别计算 并指出两种方法计算结果的有效位数。 Matlab程序如下： function [si,sd]=S(N) format long; si=0;sd=0; for j=N:-1:2 si=1.0e6/(j^2-1)+si; end for j=2:N sd=1.0e6/(j^2-1)+sd; end end 在matlab命令窗口中输入：[si,sd]=S(10000) 运行结果：si =7.499000049995000e+005 sd =7.499000049994994e+005 在matlab命令窗口中输入：[si,sd]=S(1000000) 运行结果：si =7.499990000005000e+005 sd =7.499990000005200e+0051 结果分析：si为从大到小的顺序求和的值，sd为从小到大的顺序求和的值。当N分别为10000和1000000时，si分别为7.499000049995000e+005和7.499990000005000e+005，可以看出这两个数的有效值均为13位；而sd分别为7.499000049994994e+005和7.499990000005200e+005，这两个数的有效值均为16位。这就出现了我们在矩阵理