计算方法实验(c语言_六项试验).docVIP

计算方法实验报告 实验一 舍入误差与数值稳定性 目的与要求： 通过上机编程，复习巩固以前所学程序设计语言； 通过上机计算，了解舍入误差所引起的数值不稳定性。 通过上机计算，了解运算次序对计算结果的影响，从而尽量避免大数吃小数的现象。 实验内容： 通过正反两个实例的计算，了解利用计算机进行数值计算中舍入误差所引起的数值不稳定性，深入理解初始小的舍入误差可能造成误差积累从而对计算结果的巨大影响。 通过实际编程，了解运算次序对计算结果的影响，了解实数运算符合的结合律和分配律在计算机里不一定成立。 1 对 n = 0,1,2,…,20 计算定积分 = 算法 1 利用递推公式 = - 5 n = 1,2,…,20 取 ln6- ln5 0.182 322 算法 2 利用递推公式 n = 20,19,…,1 注意到 取 算法一： #includestdio.h #includeconio.h #includemath.h #includestdlib.h void main(void) { float y_0=log(6.0)-log(5.0),y_1; int n=1; system(cls);/*清屏*/ printf(y[0]=%-20f,y_0); while(1) { y_1=1.0/n-5*y_0; printf(y[%d]=%-20f,n,y_1);/*输出*/ if(n=20)break; y_0=y_1; n++; if(n%3==0)printf(\n); } getch();/*保持用户屏幕*/ } 运行结果： Y[0]=0.182322 Y[1]=0.088392 Y[2]=0.058039 Y[3]=0.043128 Y[4]=0.034310 Y[5]=0.028448 Y[6]=0.024428 Y[7]=0.020719 Y[8]=0.021407 Y[9]=0.004076 Y[10]=0.079618 Y[11]=-0.307181 Y[12]=1.619237 Y[13]=-8.019263 Y[14]=40.167744 Y[15]=-200.772069 Y[16]=1003.922729 Y[17]=-5019.554688 Y[18]=25097.828125 Y[19]=-125489.085938 Y[20]=627445.500000 算法二： #includestdio.h #includeconio.h #includemath.h #includestdlib.h main() { float y_0=(1/105.0+1/126.0)/2,y_1; int n=20; system(cls);/*清屏*/ printf(y[20]=%-20f,y_0); while(1) { y_1=1/(5.0*n)-y_0/5.0; printf(y[%d]=%-20f,n-1,y_1); if(n=1)break; y_0=y_1; n--; if(n%3==0)printf(\n); } getch(); } 运行结果： Y[20]=0.008730 Y[19]=0.008254 Y[18]=0.008876 Y[17]=0.009336 Y[16]=0.009898 Y[15]=0.010520 Y[14]=0.011229 Y[13]=0.012040 Y[12]=0.012977 Y[11]=0.014071 Y[10]=0.015368 Y[9]=0.016926 Y[8]=0.018837 Y[7]=0.021233 Y[6]=0.024325 Y[5]=0.028468 Y[4]=0.034306 Y[3]=0.043139 Y[2]=0.058039 Y[1]=0.088392 Y[0]=0.182322 说明：从计算结果可以看出，算法1是不稳定的，而算法2是稳定的。 计算方法实验报告 实验二 拉格朗日插值与牛顿插值 一、目的与要求： 熟悉拉格朗日插值多项式和牛顿插值多项式，注意其不同特点； 实验内容： 通过拉格朗日插值和牛顿插值多项式的两个实例的计算，了解两种求解方法，分析他们的优缺点。 程序与实例 算法 输入x,y(i=0,1,2,(,n),令L(x)=0; 对=0,1,2,(,n计算 ll(x)= L( L+l(x)y 程序与实例 已知函数表 x 0.56160 0.56280 0.56401 0.56521 y 0.82741 0.82659 0.82577 0.82495 用