duoxiangshi.rar_二次多项式

#include<iostream.h> #include<fstream.h> #include<math.h> #define num 100 float x[num],y[num],error[num],sumX[num],sumY[num],A[num],xi[num][num],x_0,y_0; int n,cishu; float f1(float a,float b,float c,float x); float power(float &a,int n) { float b=1; for(int i=0;i<n;i++) { b*=a; } return b; } void in(); void out(); void qiujie(); void main() { in(); int i,j,cishu; sumX[1]=sumY[1]=0; for(i=1;i<=n;i++) { sumX[1]+=x[i-1]; sumY[1]+=y[i-1]; } cishu=2; sumX[0]=(float)n; for(i=2;i<=2*cishu;i++) { for(j=1;j<=n;j++) { sumX[i]+=power(x[j-1],i); } } for(i=2;i<=cishu+1;i++) { for(j=1;j<=n;j++) { sumY[i]+=power(x[j-1],i-1)*y[j-1]; } } //给未知量为ai的方程组的系数赋值 for (i=0;i<cishu+1;i++) { for(j=0;j<cishu+1;j++) {xi[i][j]=sumX[j+i];} xi[i][j]=sumY[i+1]; } //以下是用列主元素法求得函数的系数ai int k,c,r,p; float d,s,l[200][200]; for(i=0;i<cishu+1;i++) { //选取列主元素 for(j=i+1;j<cishu+1;j++) { c=(fabs(xi[j][i])>fabs(xi[i][i]))?j:i; for(p=0;p<cishu+2;p++) {s=xi[i][p]; xi[i][p]=xi[c][p]; xi[c][p]=s;} } // 高斯法消去，化为上三角 for(k=i+1;k<cishu+1;k++) { l[k][i]=xi[k][i]/xi[i][i]; for(r=i;r<cishu+2;r++) xi[k][r]=xi[k][r]-l[k][i]*xi[i][r]; } } //回代,求解方程组的解 A[cishu]=xi[cishu][cishu+1]/xi[cishu][cishu]; for(i=cishu-1;i>=0;i--) { d=0; for(j=i+1;j<cishu+1;j++) d=d+xi[i][j]*A[j]; A[i]=(xi[i][cishu+1]-d)/xi[i][i]; } //误差分析 for (i=0;i<n;i++) error[i]=fabs(f1(A[0],A[1],A[2],x[i])-y[i])/y[i]; y_0=f1(A[0],A[1],A[2],x_0); out(); } float f1(float a,float b,float c,float x) { float e; e=a+b*x+c*x*x; return e; } void in() { int i; float b; ifstream infile("in4.dat",ios::in); infile>>b; n=(int)b; for (i=0;i<n;i++) { infile>>b; x[i]=b; } for(i=0;i<n;i++) { infile>>b; y[i]=b; } infile>>b; x_0=b; infile.close (); } void out() { ofstream outfile("out4.dat",ios::out); outfile<<"输入结点如下："; outfile<<'\n'; outfile<<"x="; for(int i=0;i<n;i++) { outfile<<x[i]<<'\t'; } outfile<<'\n'<<"y="; for(i=0;i<n;i++) { outfile<<y[i]<<'\t'; } outfile<<'\n'<<"待求结点x为："; outfile<<x_0<<'\t'<<"对应的y为："<<y_0<<'\n'; outfile<<"相对误差如下："<<'\n'; for(i=0;i<n;i++) { outfile<<x[i]<<'\t'<<error[i]*100<<'%'<<'\n'; } outfile.close(); }