计算方法实验报告.docxVIP

计算方法实验报告班级：学号：姓名：成绩：1 舍入误差及稳定性一、实验目的（1）通过上机编程，复习巩固以前所学程序设计语言及上机操作指令；（2）通过上机计算，了解舍入误差所引起的数值不稳定性二、实验内容1、用两种不同的顺序计算，分析其误差的变化2、已知连分数，利用下面的算法计算：写一程序，读入计算并打印3、给出一个有效的算法和一个无效的算法计算积分4、设，已知其精确值为（1）编制按从大到小的顺序计算的程序（2）编制按从小到大的顺序计算的程序（3）按两种顺序分别计算并指出有效位数三、实验步骤、程序设计、实验结果及分析1、用两种不同的顺序计算，分析其误差的变化（1）实验步骤：分别从1~10000和从10000~1两种顺序进行计算，应包含的头文件有stdio.h和math.h（2）程序设计：#includestdio.h#includemath.hvoid main(){ float n=0.0;inti; float m=0.0;int j; for(i=1;i=10000;i++) //顺序计算 { n=n+1.0/(i*i); }printf(%-100f,n);printf(\n); for(j=10000;j=1;j--) //逆序计算 { m=m+1.0/(j*j); }printf(%-7f,m);printf(\n);}（3）实验结果：程序运行结果：2、已知连分数，计算：（1）实验步骤：利用，，计算（2）程序设计#includestdio.hvoid main(){float f,b[100],a[100],d[100];inti,n;printf(n=);scanf(%d,n);for(i=0;in;i++){printf(b[%d]=,i);scanf(%f,b[i]);printf(a[%d]=,i);scanf(%f,a[i]);}for(i=n-1;i=0;i--){d[n]=b[n];d[i]=b[i]+a[i+1]/d[i+1];f=d[0];}printf(%f,f);}（3）实验结果程序运行结果：3、给出一个有效的算法和一个无效的算法计算积分（1）实验步骤利用C语言编写程序，分别使用数值稳定的和数值不稳定的计算公式所建立的递推公式进行计算。（2）程序设计%无效的算法%#includestdio.h#includemath.hvoid main(){double y_0=(1/11+1/55)/2,y_1;int n=10;printf(y[10]=%-20f,y_0);while (1){y_1=4/n-16*y_0;printf(y[%d]=%-20f,n-1,y_1);if (n=1) printf(\n );y_0=y_1;n--;if (n%3==0) printf(\n);}}程序运行结果：%有效的算法%#includestdio.h#includemath.hvoid main(){double y_0=log(5),y_1;int n=1;printf(y[0]=%-20f,y_0);while (1){y_1=(1/4)*(1/n-(1/4)*y_0);printf(y[%d]=%-20f,n,y_1);if (n=10) break;y_0=y_1;n++;if (n%2==0) printf(\n);}getch();}程序运行结果：4、设，已知其精确值为（1）实验步骤先编程按从大到小的顺序计算的程序，再编程按从小到大的顺序计算的程序，然后按两种顺序分别计算。（2）程序设计1.从大到小#includestdio.h#includemath.h#includeiostreamusing namespace std;void main(){float n=0.0;inti;int N;coutplease input Nendl;cinN;for(i=N;i1;i--){n=n+1.0/(i*i-1);N=N-1;}printf(%-100f,n);printf(\n);}程序运行结果：N=20时N=5时N=1530时2.从小到大#includestdio.h#includemath.h#includeiostreamusing namespace std;void main(){float n=0.0;inti;int N;coutplease input Nendl;cinN;for(i=2;i=N;i++){n=n+1.0/(i*i-1);}printf(%-100f,n);printf(\n);}程序运行结果：N=20时N=5时N=1530时（3）①用从小到大的C语言算法计算S1000，S10000，S30000，结果如下：N=1000时：N=10000时：N=30000