欢迎来到天天文库
浏览记录
ID:59826134
大小:86.91 KB
页数:3页
时间:2020-11-24
《曲线拟合最小二乘法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(do
2、uble*));if(!p)return0;for(i=0;i>n;cout<<"请?输º?入¨?节¨²点Ì?值¦Ì(ꡧXi)ê?:êo"<3、e*X;X=newdouble[n];for(i=0;i>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、){A[0][0]+=1;//cout<5、[i]*Y[i];}for(k=0;k<2;k++){for(i=k+1;i<3;i++){doubleT=A[i][k]/A[k][k];B[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<<"拟6、a合?曲¨²线?为a:êo"<
3、e*X;X=newdouble[n];for(i=0;i>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、){A[0][0]+=1;//cout<5、[i]*Y[i];}for(k=0;k<2;k++){for(i=k+1;i<3;i++){doubleT=A[i][k]/A[k][k];B[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<<"拟6、a合?曲¨²线?为a:êo"<
4、){A[0][0]+=1;//cout<5、[i]*Y[i];}for(k=0;k<2;k++){for(i=k+1;i<3;i++){doubleT=A[i][k]/A[k][k];B[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<<"拟6、a合?曲¨²线?为a:êo"<
5、[i]*Y[i];}for(k=0;k<2;k++){for(i=k+1;i<3;i++){doubleT=A[i][k]/A[k][k];B[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<<"拟
6、a合?曲¨²线?为a:êo"<
此文档下载收益归作者所有