計算方法与实习实验报告舍入误差与数值稳定性.docxVIP

实验一 舍入误差与数值稳定性 用两种不同的顺序计算≈1.644834，分析其误差的变化；程序：#includestdio.h#includemath.hvoid main(){int i;float sum_0=0,sum_1=0;for(i=1;i10001;i++)sum_0+=1.0/(i*i);printf(%f\n,sum_0);for(i=10000;i0;i--)sum_1+=1.0/(i*i);printf(%f\n,sum_1);}输出：sum_0=1.644725sum_1=1.644834分析讨论：在计算机内作加法运算时，首先要对加数作对阶处理。由于加数和被加数相比是一个小量，加之计算机字长有限，加数被作为机器零处理了，与被加数相加实际上是被被加数“吃掉”了。因此多个数相加，应按绝对值从小到大的顺序依次相加，可以保留更多的有效位数。实验二 方程求根——二分法 一、目的和要求1）通过对二分法的编程练习，掌握方程求根的二分法的算法；2）通过对二分法的上机运算，进一步体会二分法的特点。二、实习内容1）二分法的编程实现。2）进行有根区间和误差限的比较和讨论。三、算法流程图：计算：fa=f(a), fb=f(b), fab=fa*fb计算：c=(a+b)/2, fc=f(c)|fc|DELTA|| b-aEPS?(1)fa*fc0 or(2)fc*fb0?b=cfb=f