计算机的思想.docVIP

数值分析上机题 姓名：徐玉珉 学号：093216 第一章 20．（上机题）舍入误差与有效数 设，其精确值为。 （1）编制按从大到小的顺序，计算的通用程序。 （2）编制按从小到大的顺序，计算的通用程序。 （3）按两种顺序分别计算，，，并指出有效位数。（编制程序时用单精度） （4）通过本上机题，你明白了什么？ 解：（1）从大到小的顺序 #include stdio.h Void main () { float j,s,sum=0,N=0; printf(input N\n); scanf(%d,N); for(j=2;j=N;j++) { s=1/(j*j-1); sum+=s; } printf(%d\n, sum); } （2）按从小到大的顺序 #include stdio.h Void main () { float j,s,sum=0,N=0; printf(input N\n); scanf(%d,N); for(j=N;j=2;j--) {s=1/(j*j-1); sum+=s; } printf(%d\n, sum); } （3） 从大到小的顺序的值 从小到大的顺序的值 精确值 有效位数 从大到小 从小到大 0.740049 0.74005 0.740049 6 5 0.749852 0.7499 0.7499 4 4 0.749852 0.749999 0.749999 3 6 （4）通过本题，可以看出按两种不同的顺序计算的结果是不相同的，按从大到小的顺序计算的值与精确值有较大的误差，而按从小到大的顺序计算的值与精确值较吻合。从大到小的顺序计算得到的结果的有效位数少。计算机在进行数值计算时会出现“大数吃小数”的现象，导致计算结果的精度有所降低，我们在计算机中进行同号数的加法时，采用绝对值较小者先加的算法，其结果的相对误差较小。 第二章 20．（上机题）Newton迭代法 （1）给定初值及容许误差，编制Newton法解方程根的通用程序。 （2）给定方程，易知其有三个根，，。 1．由Newton方法的局部收敛性可知存在，当时，Newton迭代序列收敛于根。试确定尽可能大的。 2．试取若干初始值，观察当，，，，时Newton序列是否收敛以及收敛于哪一个根。 （3）通过本上机题，你明白了什么？ 解：（1）编制的通用程序： #include stdio.h #includeiostream.h #includemath.h #define error 0.000001 //给定容许误差 float f_x(float x) //定义函数f(x) { float a; a=x*x*x/3-x; //f(x)的表达式; return(a); } float d_f(float x) //定义函数d_f(x)，计算f(x)的导函数 { float a; a=x*x-1; //f(x)导函数的表达式; return (a); } void main(void) { float x0,x1,a; int k=0; cout请输入初值x0:; cinx0; do { a=-f_x(x0)/d_f(x0); x1=x0+a; k++; x0=x1; } while(fabs(a)error); coutk\tx0; //输出迭代的次数和根值 } （2）计算迭代序列收敛于根的尽可能大的的函数为： #includeiostream.h #includemath.h void delay(int n) //定义延时函数 {for(n=10000;n0;n--);} #define error 0.000001 float f(float x) //定义函数f(x) { float f; f=x*x*x/3-x; //f(x)的表达式; return(f); } float d_f(float x) //定义函数df(x)，计算f(x)的导函数 { float df; df=x*x-1; //f(x)导函数的表达式; 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