数值上机作业重点分析.docx

《数值分析》上机作业 （第一二三章） 学院：电 气 工 程 学 院 班级： 电气13级硕士2班 教师： 石佩虎 老师 姓名： 高蕾 学号： 132046 第一章 实验1 舍入误差与有效数 设，其精确值为。 (1) 编制按从大到小的顺序 QUOTE SN=122-1+132-1+…+1N2-1 ，计算 QUOTE SN 的通用程序； (2) 编制按从小到大的顺序 QUOTE SN=1N2-1+1(N-1)2-1+…+122-1 ，计算的通用程序； (3) 按两种顺序分别计算 QUOTE S102 、、，并指出有效位数（编制程序时用单精度）； (4) 通过本上机题你明白了什么？ 解答如下： (1). 按从大到小的顺序计算 QUOTE SN 的通用程序如下所示： n=input(Please Input an N (N1):); y=0; accurate=1/2*(3/2-1/n-1/(n+1)); %精确值 for i=2:1:n %从大到小的顺序 x=1/(i^2-1); x=single(x); y=y+x; end error= accurate-y; format long; disp(____________________________________________________); disp(The value of Sn from large to small is:); disp(y); disp(The value of error is:); disp(error); (2) 编制按从小到大的顺序计算的通用程序如下所示： n=input(Please Input an N (N1):); y=0; accurate=1/2*(3/2-1/n-1/(n+1)); for i=n:-1:2 x=1/(i^2-1); x=single(x); y=y+x; end error= accurate-y; format long; disp(____________________________________________________); disp(The value of Sn from large to small is:); disp(y); disp(The value of error is:); disp(error); (3) 计算结果： 按从大到小的顺序计算得： N误差值有效数字位数0.740049570.749852140.74985213 按从小到大的顺序计算得： N误差值有效数字位数0.740049570.749900070.74999907 (4) 总结： 当我们采用不同的计算顺序，对于同一个计算式，会得出不同的结果。本例中，分别采用从大到小和从小到大两种顺序，结果有较大差异，有效位数最少的情况只有3位。 当按从大到小的顺序进行时，当越大， QUOTE 1N2-1 就越小，与已计算出的和值相比数量及特别小，很可能被舍去为0，使得计算结果变小，有效位数减少，结果不精确。这种现象随着越大，结果的差异就越大，误差就越大。而从小加到大时，越往后小数之间数量级相差不大，前面的和值与大数的数量级相差不大，之前的情况不会出现，误差都是符合要求的，有效数字也有7位。 所以在进行数值计算时，若忽略运算次序，可能会造成大数吃小数的问题，影响到结果的准确性。所以我们要根据计算式的性质，设计合理有效的算法，进行计算。 第二章 实验2 Newton迭代法 (1)给定初值 QUOTE x0 及容许误差，编制Newton法解方程 QUOTE fx=0 根的通用程序。 (2)给定方程 QUOTE fx=x33-x=0 ，易知其有三个根。  = 1 \* GB3 ①由Newton方法的局部收敛性可知存在 QUOTE δ0 ，当 QUOTE x0∈(-δ，δ) 时Newton迭代序列收敛于根x2*，试确定尽可能大的；  = 2 \* GB3 ②试取若干初始值，观察当时Newton