数值计算—拟合曲线.docVIP

实验报告 班级：计科1101 姓名：杨玲 学号：201126100130 一、实验目的 用最小二乘法求多项式拟合曲线。 二、实验方法 1.首先用Matlab绘制出拟合曲线，得到拟合函数。 2.用C++语言计算出系数，就是根据书上的公式，计算出各个值，然后再利用之前的顺序消元法求出各个系数值。 三、实验代码 clear clc x=[1,3,4,5,6,7,8,9,10] y=[10,5,4,2,1,1,2,3,4] z=polyfit(x,y,2) plot(x,y,r-*) hold on; f=poly2sym(z) ezplot(f,[0,10]) #include iostream #include iomanip using namespace std; void main() { int x[9]={1,3,4,5,6,7,8,9,10}, y[9]={10,5,4,2,1,1,2,3,4}; double sum1=0; for(int p=0;p9;p++) { sum1+=x[p]; } double sum2=0; for(int q=0;q9;q++) { sum2+=x[q]*x[q]; } double sum3=0; for(int g=0;g9;g++) { sum3+=x[g]*x[g]*x[g]; } double sum4=0; for(int w=0;w9;w++) { sum4+=x[w]*x[w]*x[w]*x[w]; } double sum5=0; for(int r=0;r9;r++) { sum5+=y[r]; } double sum6=0; for(int t=0;t9;t++) { sum6+=x[t]*y[t]; } double sum7=0; for(int f=0;f9;f++) { sum7+=x[f]*x[f]*y[f]; } //=================================================== float a[3][4]={9,sum1,sum2,sum5,sum1,sum2,sum3,sum6,sum2,sum3,sum4,sum7},XX[3]; //--------------------------------------- for (int k=0;k2;k++) { for (int i=k+1;i3;i++) { float l=a[i][k]/a[k][k]; for (int j=k+1;j4;j++) { a[i][j]=a[i][j]-l*a[k][j]; } } } //消元 for (int m=2;m=0;m--) { float sum=0; for (int j=m+1;j3;j++) { sum+=a[m][j]*XX[j]; } XX[m]=(a[m][3]-sum)/a[m][m]; } //计算 //=========================================== for (int h=0;h3;h++) { coutXX[h+1]=XX[h]endl; } couty=XX[2]x^2XX[1]x+XX[0]endl; cout最低点位置:((-1)*XX[1]/(2*XX[2]),(4*XX[0]*XX[2]-XX[1]*XX[1])/(4*XX[2]))endl; } 四、实验结果及其讨论 实验结果有很小的误差，这是计算机本身的保留的小数位引起的 五、总结 通过本次实验，我开始会使用了Matlab软件，并且产生了浓厚的兴趣！ 另外，这次实验的代码，可能有点笨拙，不是很技巧，因为是直接计算出各个值在进行计算的，循环的次数很多，效率不是很高。