《数值分析》上机题.doc

《数值分析》上机实习题 秦九韶算法 1．源程序 #includestdio.h #includeiostream.h //输出文件指针 FILE *stream = NULL; typedef float * pFloat; //和用户交互输入多项式各系数 bool InputParam(int n,pFloat a,float x) { cout秦九韶算法\n; cout请输入最高次数n:; cinn; if(n=0) return false; a = new float[n+1]; for(int i=0;i=n;i++) { cout请输入a[i]:; cina[i]; } cout请输入x的值:; cinx; return true; } //按照格式打印多项式到文件中 void OutputTheMultinomial(int n,pFloat a,float x) { fprintf(stream,此多项式为:\n); fprintf(stream,P(x)=); for(int i=n;i=2;i--) { fprintf(stream,(%f)*x**%d+,a[i],i); } fprintf(stream,(%f)*x+(%f)\n,a[1],a[0]); fprintf(stream,x=%f,x); } //根据多项式各系数及给定的x值计算 float QinJiuShao(int n,pFloat a,float x) { float r=a[n]; for(int i=n-1;i=0;i--) { r = x*r+a[i]; } return r; } //计算第一题 void Calculate1() { fprintf(stream,秦九韶计算结果1\n); int n = 0; float *a = NULL;//指向存放多项式系数的指针 float x = 0; char next = Y; //和用户交互进行输入并计算 while(next == Y || next == y) { if(InputParam(n,a,x)) { //输出多项式 OutputTheMultinomial(n,a,x); //输出结果 fprintf(stream,\n计算结果为:P(%f)=%f\n,x,QinJiuShao(n,a,x)); if(a) { delete []a; a=NULL; } } cout是否进行下一次运算?(Y/N):; cinnext; } } //计算第二题 void Calculate2() { fprintf(stream,秦九韶计算结果2\n); //多项式系数 float a[] = {2.0f,-3.36f,2.81f,-2.74f,-3.27f,0.22f}; int n = 5; float x = 0; //按每次递增0.05计算 for(int i = 0;i = 20;i++) { x=0.8+0.05*i; //输出多项式 OutputTheMultinomial(n,a,x); //输出结果 fprintf(stream,\n计算结果为:P(%f)=%f\n,x,QinJiuShao(n,a,x)); } } void main() { //打开输出结果文件 if((stream = fopen( result.txt, w )) == NULL ) { printf( 打不开输出文件\n ); return; } //通过用户输入计算第一题 Calculate1(); //通过循环生成 x值计算第二题 Calculate2(); //关闭文件 if(stream) fclose(stream); } 2．输出结果 秦九韶计算结果1 此多项式为: P(x)=(7.540000)*x**4+(11.080000)*x**3+(3.820000)*x**2+(0.440000)*x+(-0.480000) x=3.250000 计算结果为:P(3.250000)=1262.865112 此多项式为: P(x)=(2.790000)*x**4+(9.850000)*x**3+(14.150000)*x**2+(5.380000)*x+(7.240000) x=3.250000 计算结果为:P(3.250000)=823.586670 此多项式为: P(x)=(9.3600