数值分析实验报告(七).doc

《数值分析实验报告(七).doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、数值分析实验报告（七）%1.实验名称：曲线拟合与函数逼近。%1.实验目的：最小二乘法解决曲线拟合问题0用最小二乘法求拟合下面的函数表X・1.00-0.75-0.50-0.2500.25y-0.22090.32950.88261.43922.0032.5645四.程序：#include#includeusingnamespacestd;intArray(double***Arr,intn){double**p;inti;p=(double**)malloc(n*sizeof(double*));if

2、(!p)return0;for(i=0;i

3、输入节点函数值：”vvendl;double*Y;Y=newdouble[n];for(i=0;i

4、Al存分配失败";elsefor(i=0;i<3;i++){for(j=();j<3;j++){A[i]U]=0;for(i=0;i

5、*X[i];A[2][0]=A[0][2];A[2][1]=A[1][2];A⑵[2]+=X[i]*X[i]*X[i]*X[i];1}for(i=0;i<3;i++){B[i]=O;}for(i=0;i

6、A[k]U];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];}coutvv”拟合曲线为:”vvendl;cout«ny=n«S[0]«n+n«S[l]«Hx+n«S[2]«HxA2u«endl«endlvv”误差为:n«endl;fbr(i=0;i

7、]+S[2]*X[i]*X[i]-Y[i]«endl;}五•运行结果：wE:caihanlinS值计算Debug7.exe11共有N个节点：6请输入节点值:-1.00-0.75-0.50-0.2500.25请输入节点函数值:-0.22090.32950.88261.43922.0032.5645拟合曲线为:y=2.00221+2.22904x+0.0280857x*2误差为：YO-Y(tO)=0.248986Yl-Y(tl)=0.0167249Y2-Y(t2)=l.11259Y3-Y(t3)=0.564762Y4-Y(t4

8、)=-0.000793333Y5-Y(t5)=-O.560538Pressanykeytocontinue六.体会：通过这次报小二乘法的上机实验，我对最小二乘法有了更深刻的认识，同时在编写程序的时候，遇到不少的困难，报后在同学的帮助下，顺利解决问题。信计（1）班010108蔡汉林