计算方法实验报告

《计算方法实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、计算方法实验报告班级：学号：姓名：成绩:1舍入误差及稳定性一、实验目的（1）通过上机编程，复AI巩岡以前所学程序设汁语言及上机操作指令;（2）通过上机计算，了解舍入误差所引起的数值不稳定性二、实验内容1、用两种不M的顺序计算分析其误差的变化/?=!2、已知连分数/=%+,利用下而的算法计算4b}+a2l(b2+%/(•••+“"lb,)、Md：(/=n—1，"—2，"”0)/=d()写一程序,读入…，义，计算并打印/3、给!I!一个有效的算法和一个无效的算法计算积分dx(n=0，l”“，10)^72-1已知其精确值为1~N（1）编制按从大到小的顺序计算的程序（2）编制按从小到人的顺序计算S

2、/V的税序（3）按两种顺序分别计算S眶），Snxw，53（麵，并指出冇效位数三、实验步骤、程序设计、实验结果及分析1、用两种不同的顺序计算分析其误差的变化/:=!（1）实验步骤：分別从1〜10000和从10000〜1两种顺序进行计算，应包含的尖文件有stdio.h和math.h（2）程序设计：#includcvoidmain(){floatn=0.0;inti;floatm=0.0;intj;for(i=1;i<=l0000;i++)//顺序计算{n=n+1.0/(i*i);}printf("%-100f*，n);printf(,,u);for(j=10000;j>=1

3、;j-)//逆序计算{m=m+1.0/(j*j);}printf(n%-7fm);printf("");}(3)实验结果:程序运行结果：‘C:Docu>entsandSettingsLX205001Docu>entsDebug.6447251.644834I*ressanykeytocontinue2、已知连分数/=么+7———_-b}+a2/[b2+a3/(...+“"/bt(1)实验步骤：利用4=么，4=夂+~(Z=n-1，AZ—2”.”0),f=d0,计算/(2)程序设计#includevoidmain(){floatf，b[100]，a[100]，d[

4、100];inti,n;printf("n=");scanf("%d"，&n);for(i=0;i=0;i—){d[n]=b[n];d[i]=b[i]+a[i+l]/d[i+l】;f=d[O];}printf("%r，f);}(3)实验结果程序运行结果：c:、*C:Docu>entsandSettingsLX205001MyDocu>entsDebugn=4b[0]=2a[0]=45b[l]=16a

5、tl]=8b[2]=45at2]=lb[3]=6a[3]=92.499326Pressanykeytocontinue3、给出一个有效的算法和一个无效的算法计算积分4x+l(1)实验步骤利用C语言编写程序，分别使用数值稳定的和数值不稳定的计算公式所建立的递推公式进行计算。(2)程序设计%无效的算法％#includevoidmain()doubley_0=(l/ll+l/55)/2，y_l;intn=10;printf("y[10]=%-20广，y_0);while(1){y_l=4/n-16*y_0;printf(”y[%d]=%-20r，n-l，y_l);if(n<=l)

6、printf(n");y_o=y_i;n-；if(n%3==0)printf(HH);}}程序运行结果：CA*C:Docu>entsandSettingsLX205001lyDocu>entsDebugVii"算方法.exe*ly[10]=0-000000y[91=0.000000yLS]=0-000000y[71=0-000000y[61=0.000000p[5]=0-000000V[41=0-000000y[31=1.000000V[21=-15.000000y[l1=242.000000ST[0]=-3868.000000Pressanykeytocontinue%有

7、效的算法％#include#includevoidmain(){doubley_0=log(5),y_l;intn=l;printf(,,y[0]=%-20fn,y_0);while(1){y_l=(l/4)*(l/n-(l/4)*y_0);printf(,,y[%d]=%-20f,,n,y_l)；if(n>=10)break;y_o=y_i;n++;if(n%2==0)printf(H