欢迎来到天天文库
浏览记录
ID:56731225
大小:62.50 KB
页数:3页
时间:2020-07-06
《曲线拟合最小二乘法 C++程序.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、课题八曲线拟合的最小二乘法实验目标:在某冶炼过程中,通过实验检测得到含碳量与时间关系的数据如下,试求含碳量y与时间t内在关系的拟合曲线。t05101520253035404550Y01.272.162.863.443.874.154.374.514.584.02最小二乘法:#include#includeusingnamespacestd;intArray(double***Arr,intn){double**p;inti;p=(double**)malloc(n*sizeof(double*));if(!p)retu
2、rn0;for(i=0;i>n;cout<<"请?输º?入¨?节¨²点Ì?值¦Ì(ꡧXi)ê?:êo"<3、cin>>X[i];}cout<<"请?输º?入¨?节¨²点Ì?函¡¥数ºy值¦Ì(ꡧYi)ê?:êo"<>Y[i];}if(!Array(&A,3))cout<<"内¨²存ä?分¤?配?失º¡ì败㨹!ê?";else{for(i=0;i<3;i++){for(j=0;j<3;j++){A[i][j]=0;}}for(i=0;i4、i]*X[i];A[1][0]=A[0][1];A[1][1]=A[0][2];A[1][2]+=X[i]*X[i]*X[i];A[2][0]=A[0][2];A[2][1]=A[1][2];A[2][2]+=X[i]*X[i]*X[i]*X[i];}}for(i=0;i<3;i++){B[i]=0;}for(i=0;i5、i]=B[i]-T*B[k];for(j=k+1;j<3;j++){A[i][j]=A[i][j]-T*A[k][j];}}}S[2]=B[2]/A[2][2];for(i=1;i>=0;i--){doubleTemp=0;for(intj=i+1;j<2;j++)Temp=Temp+A[i][j]*S[j];S[i]=(B[i]-Temp)/A[i][i];}cout<<"拟a合?曲¨²线?为a:êo"<6、êo"<
3、cin>>X[i];}cout<<"请?输º?入¨?节¨²点Ì?函¡¥数ºy值¦Ì(ꡧYi)ê?:êo"<>Y[i];}if(!Array(&A,3))cout<<"内¨²存ä?分¤?配?失º¡ì败㨹!ê?";else{for(i=0;i<3;i++){for(j=0;j<3;j++){A[i][j]=0;}}for(i=0;i4、i]*X[i];A[1][0]=A[0][1];A[1][1]=A[0][2];A[1][2]+=X[i]*X[i]*X[i];A[2][0]=A[0][2];A[2][1]=A[1][2];A[2][2]+=X[i]*X[i]*X[i]*X[i];}}for(i=0;i<3;i++){B[i]=0;}for(i=0;i5、i]=B[i]-T*B[k];for(j=k+1;j<3;j++){A[i][j]=A[i][j]-T*A[k][j];}}}S[2]=B[2]/A[2][2];for(i=1;i>=0;i--){doubleTemp=0;for(intj=i+1;j<2;j++)Temp=Temp+A[i][j]*S[j];S[i]=(B[i]-Temp)/A[i][i];}cout<<"拟a合?曲¨²线?为a:êo"<6、êo"<
4、i]*X[i];A[1][0]=A[0][1];A[1][1]=A[0][2];A[1][2]+=X[i]*X[i]*X[i];A[2][0]=A[0][2];A[2][1]=A[1][2];A[2][2]+=X[i]*X[i]*X[i]*X[i];}}for(i=0;i<3;i++){B[i]=0;}for(i=0;i5、i]=B[i]-T*B[k];for(j=k+1;j<3;j++){A[i][j]=A[i][j]-T*A[k][j];}}}S[2]=B[2]/A[2][2];for(i=1;i>=0;i--){doubleTemp=0;for(intj=i+1;j<2;j++)Temp=Temp+A[i][j]*S[j];S[i]=(B[i]-Temp)/A[i][i];}cout<<"拟a合?曲¨²线?为a:êo"<6、êo"<
5、i]=B[i]-T*B[k];for(j=k+1;j<3;j++){A[i][j]=A[i][j]-T*A[k][j];}}}S[2]=B[2]/A[2][2];for(i=1;i>=0;i--){doubleTemp=0;for(intj=i+1;j<2;j++)Temp=Temp+A[i][j]*S[j];S[i]=(B[i]-Temp)/A[i][i];}cout<<"拟a合?曲¨²线?为a:êo"<6、êo"<
6、êo"<
此文档下载收益归作者所有