gaosinihe.rar_函数拟合预测_拟合 预测_高斯函数 拟合_高斯拟合程序

#include<iostream> #include<math.h> #include<iomanip> using namespace std; double max(double array[3]); double gt(double x[11],double y[11]); double GS(double a[3][3],double b[3],double e[],double s[]); void main() { double a[11]={1,1,1,1,1,1,1,1,1,1,1},b[11],c[11],e[11],f[11],g[11]; cout<<"请输入变量t值:"; for(int i=0;i<11;i++) { cin>>b[i]; c[i]=b[i]*b[i]; } cout<<"请输入y值:"; for(int i=0;i<11;i++) cin>>e[i]; f[0]=gt(a,a),f[1]=gt(a,b),f[2]=gt(a,c),f[3]=gt(b,b),f[4]=gt(b,c),f[5]=gt(c,c); g[0]=gt(a,e),g[1]=gt(b,e),g[2]=gt(c,e); double array[3][3]={{f[0],f[1],f[2]},{f[1],f[3],f[4]},{f[2],f[4],f[5]}}; GS(array,g,e,b); } double gt(double x[11],double y[11]) { double sum=0; for(int i=0;i<11;i++) sum=x[i]*y[i]+sum; return sum; } double GS(double a[3][3],double b[3],double e[11],double s[11]) { double y[11]={0}; double c[3]={0}; double x0[3]={0}; int i,j,k; double x[3]; double r,sum=0; double w; for(k=1;;k++) { for(i=0;i<3;i++) { for(j=0;j<3;j++) { sum=a[i][j]*x0[j]+sum; } x[i]=x0[i]+(b[i]-sum)/a[i][i]; c[i]=fabs(x[i]-x0[i]); x0[i]=x[i]; sum=0; } r=max(c); if(r<0.00000001) { cout<<"解得的二次拟合曲线系数为:\n"; for(i=0;i<3;i++) cout<<"A"<<i<<"="<<x[i]<<endl; break; } } for(i=0;i<11;i++) { y[i]=x[0]+x[1]*s[i]+x[2]*s[i]*s[i]; w=fabs(y[i]-e[i]); cout<<y[i]<<" "<<e[i]<<endl; cout<<"Y"<<i<<"的误差为"<<w<<setprecision(15)<<endl; } return 0; } double max(double array[3]) { double a=array[0]; int i; for(i=0;i<3;i++) { if(a<array[i]) a=array[i]; } return a; }