实验5大作业.docVIP

不同拟合模型的比较 小组成员：李伟庆，康家磊，苏善杰，李樯，高行行， 实验内容 实验结果 方案一 Q=28.1658 方案二P(X)= 1 56.25 0 3.596 2 2.848 2.5673 1.17781 3 2.79043 2.11707 -1.67586 Q=2631.04 方案三P(X)=+ 1 118.81 0 3.596 2 2.848 2.5673 1.17781 3 2.79043 2.11707 -1.67586 4 2.21603 1.69087 -1.33819 Q=2140.91 方案四 A=89.9 b=-1456.01 Q=107232 实验结果分析 附代码 方案一 #includeiostream #includemath.h using namespace std; int main() { double x[25]={0,0.2,0.6,1.0,1.3,1.6,1.7,1.8,1.9,2.2,2.3,2.5,2.6,2.9,3.1,3.4,3.8,4.1,4.4,4.7,4.8,4.9,5.0,5.1,5.3}; double y[25]={0,-2.5,-4.0,-5.7,-3.5,-2.0,-1.0,2,3.5,4.0,7.0,7.5,9.9,10.9,11.9,13.5,13.0,11.9,9.0,6.5,4.0,1.5,0.0,-2.5,-5}; double G[25][4],GT[4][25],GTG[4][4]={0},b[4]={0}; for(int j=0;j4;j++) { for(int i=0;i25;i++) { G[i][j]=pow(x[i],j); GT[j][i]=G[i][j]; } } for(int i=0;i4;i++) { for(int j=0;j4;j++) { for(int k=0;k25;k++) { GTG[i][j]=GTG[i][j]+GT[i][k]*G[k][j]; } } for(int l=0;l25;l++) { b[i]=b[i]+GT[i][l]*y[l]; } } double u[4][4]={0},l[4][4]={0}; for(int k=0;k4;k++) { for(int j=k;j4;j++) { u[k][j]=GTG[k][j]; for(int s=0;sk;s++) { u[k][j]=u[k][j]-l[k][s]*u[s][j]; } for(int i=k+1;i4;i++) { l[i][k]=GTG[i][k]; for(int s=0;sk;s++) { l[i][k]=l[i][k]-l[i][s]*u[s][k]; } l[i][k]=l[i][k]/u[k][k]; } } } double Y[4],a[4]={0}; for(k=0;k4;k++) { Y[k]=b[k]; for(int s=0;sk;s++) Y[k]=Y[k]-l[k][s]*Y[s]; } for(k=3;k=0;k--) { double sum=0; for(int s=k+1;s4;s++) sum=sum+u[k][s]*a[s]; a[k]=(Y[k]-sum)/u[k][k]; } for(i=0;i4;i++) couta[i]\t; //计算Q1 double q1=0; for(i=0;i25;i++) q1=q1+pow(a[0]+a[1]*x[i]+a[2]*pow(x[i],2)+a[3]*pow(x[i],3)-y[i],2); cout\nq1 Q1结果; } 方案二（变量“jie”控制拟合多项式的次数，jie=3为3次，jie=4为4次） #includeiostream #includemath.h using namespace std; const int jie=4; int main() { double x[25]={0,0.2,0.6,1.0,1.3,1.6,1.7,1.8,1.9,2.2,2.3,2.5,2.6,2.9,3.1,3.4,3.8,4.1,4.4,4.7,4.8,4.9,5.0,5.1,5.3}; doub