计算方法上机实习题大作业(实验)探究.docVIP

计算方法实习题大作业 ，分析其误差的变化 2、已知连分数，利用下面的算法计算： 写一程序，读入计算并打印 3、给出一个有效的算法和一个无效的算法计算积分 4、设，已知其精确值为 （1）编制按从大到小的顺序计算的程序 （2）编制按从小到大的顺序计算的程序 （3）按两种顺序分别计算并指出有效位数 三、实验步骤、程序设计、实验结果及分析 1、用两种不同的顺序计算，分析其误差的变化 （1）流程图 （2）程序设计： a.顺序计算 #includestdio.h #includemath.h //顺序 void main() { double sum=0; int n=10000,i; for(i=1;i=n;i++) { sum=sum+1/pow(i,2); if(i%1000==0) printf(sum[%d]=%-30f\n,i,sum); } }b.逆序计算 #includestdio.h #includemath.h //逆序 void main() { double sum=0; int n=1,i; for(i=10000;i=n;i--) { sum=sum+1/pow(i,2); if(i%1000==0) printf(sum[%d]=%-30f\n,i,sum); } printf(sum[%d]=%-30f\n,i,sum); } 结果分析：两种不同顺序计算结果是一样的，顺序计算误差从一开始就很小，而逆序计算误差最开始十分大，后来结果正确。 2、已知连分数，计算： （1）流程图 （2）实验步骤： 利用 ，，计算 （3）程序设计 #includestdio.h #includemath.h void main() { float a[1024],b[1024],d[1024]; int i,n; printf(请输入n值\n); scanf(%d,n); printf(请依次输入a[i]的值\n); for(i=1;i=n;i++) { printf(a[%d]=,i); scanf(%f,a[i]); } printf(请依次输入b[i]的值\n); for(i=0;i=n;i++) { printf(b[%d]=,i); scanf(%f,b[i]); } d[n]=b[n]; for(i=n-1;i=0;i--) d[i]=b[i]+a[i+1]/d[i+1]; printf(f=%f,d[0]); } （3）实验结果 程序运行结果： 3、给出一个有效的算法和一个无效的算法计算积分 （1）实验步骤 利用C语言编写程序，分别使用数值稳定的和数值不稳定的计算公式所建立的递推公式进行计算。 （2）流程图 （3）程序设计 #includestdio.h #includemath.h void main() { int i; double y[11]; y[0] = 0.25*log(5.0); printf(有效算法输出结果：\n); printf(y[0]=%-20f, y[0]); for (i = 1; i = 10; i++) { y[i] = (double)(1 / (4.0*i) - y[i - 1] / 4.0); printf(y[%d]=%-20f, i, y[i]); if (i % 3 == 2) printf(\n); } printf(\n无效算法输出结果：\n); y[10] = (double)((1 / 11.0 + 1 / 55.0) / 2.0); printf(y[10]=%-20f, y[10]); for (i = 10; i 0; i--) { y[i - 1] = (double)(1.0 / i - 4.0 * y[i]); printf(y[%d]=%-20f, i - 1, y[i - 1]); if (i % 3 == 0) printf(\n); } } （4）实验结果及分析 程序运行结果： 结果分析：无效算法数值不稳定，误差造成的影响大 4、设，已知其精确值为 （1）实验步骤 先编程按从大到小的顺序计算的程序，再编程按从小到大的顺序计算的程序，然后按两种顺序分别计算。 （2）流程图 （3）程序设计 #includestdio.h void main() { int N; double SN[30000]; SN[30000] = (3.0 / 2.0 - 1.0 / 30000.0 - 1 / 30001.0) / 2.0; for (N = 30000; N = 2; N--)