东南大学数值分析上机题(上).doc

数值分析上机报告 姓名： 学号： 专业： 2013年10月27日 第一章 舍入误差与有效数 设，其精确值为。 （1）编制按从大到小的顺序，计算的通用程序。 （2）编制按从小到大的顺序，计算的通用程序。 （3）按两种顺序分别计算，，，并指出有效位数。（编制程序时用单精度） （4）通过本上机题，你明白了什么？ 解： （1） #includestdio.h void main() { float n,i,s; printf(please input n=); scanf(%f,n); for(i=2,s=0;i=n;s+=1/(i*i-1),i++); printf(s=%f\n,s); } （2） #includestdio.h void main() { float n,i,s; printf(please input n=); scanf(%f,n); for(i=n,s=0;i=2;s+=1/(i*i-1),i--); printf(s=%f\n,s); } 按从大到小顺序：=0.740049 有效位数6位 =0.749852 有效位数3位 =0.749852 有效位数3位 按从小到大顺序：=0.740050 有效位数5位 =0.749900 有效位数6位 =0.749999 有效位数6位 通过上述实验数据可以看出此次算法使用从小到大的顺序进行得到的数据相对而言更精确，可以得到这样的启示：在计算数值时，要先分析不同算法对结果的影响，避免大数吃小数的现象，找出能得到更精确的结果的算法。 第二章 （上机题）Newton迭代法 （1）给定初值及容许误差，编制Newton法解方程根的通用程序。 （2）给定方程，易知其有三个根，，。 1．由Newton方法的局部收敛性可知存在，当时，Newton迭代序列收敛于根。试确定尽可能大的。 2．试取若干初始值，观察当，，，，时Newton序列是否收敛以及收敛于哪一个根。 （3）通过本上机题，你明白了什么？ 解： （1）#includeiostream.h #includemath.h #define eps 0.000001 float f(float x) { float f; f=x*x*x/3-x; return(f); } float df(float x) { float df; df=x*x-1; return (df); } void main(void) { float x0,x1,a; int k=0; cout请输入初值x0:; cinx0; do { a=-f(x0)/df(x0); x1=x0+a; k++; x0=x1; } while(fabs(a)eps); coutk\tx0; } （2）1 #includeiostream.h #includemath.h void delay(int n) {for(n=10000;n0;n--);} #define eps 0.01 float f(float x) { float f; f=x*x*x/3-x; return(f); } float df(float x) { float df; df=x*x-1; return (df); } int judgement(float z) { int count=5; float x0,x1,type,type1; x0=z; while(count--0) { x1=x0-f(x0)/df(x0); type=fabs(x1); type1=fabs(x1-x0); coutcount=count\ttype=type\ttype1=type1\n; if(fabs(x1-x0)eps) return 1; x0=x1; delay(30000); } return 0; } void main(void) { float delta=0; int flag=1; while(flag==1) { cout方程的根为:\n; delta+=eps; flag=judgem