[2018年最新整理]二次样条插值及其C语言的实现.docVIP

1、实验要求：用C语言编程实现三次样条插值函数及曲线拟合的最小二乘法的求解。 2、实验目的： (、通过编程加深对三次样条插值及曲线拟合的最小二乘法的理解； (、学习用计算机解决工程问题，主要包括数据处理与分析。 3、 程序说明： (、三次样条插值函数的求法： 分析：在解线性方程组确定M值时，因为矩阵是三对角阵，所以我用克拉默法则=进行求解，其中D为三对角阵的行列式，值为2^12。 程序如下： #includestdio.h void main() { double x[12],y[12],u[11],v[11],g[12]; //x[12]用于存放点的横坐标，y[12]用于存放点的纵坐标，u[11]用于存放μ1~μ10 //v[11]用于存放λ1~λ10，g[12]用于存放g0~g11 int i,j=0; int AllD=1; //三角对阵的行列式 double Mv=1,Mu=1; //Mv表示λ1*λ2*…*λ10，Mu表示μ1*μ2*…*μ10 double M[12],X3[11],X2[11],X1[11],X0[11]; //M[12]用于存放M0~M11，X3[11]用于存放x^3的系数，X2[11]用于存放x^2的系数, //X1[11]用于存放x的系数,X0[11]用于存放常数项的系数 printf(Please input x(i):\n); for(i=0;i12;i++) x[i]=0; for(i=0;i12;i++) scanf(%lf,x[i]); printf(Please input y(i):\n); for(i=0;i12;i++) scanf(%lf,y[i]); for(i=1;i11;i++) { u[i]=(x[i]-x[i-1])/(x[i+1]-x[i-1]); //μ1~μ10的求解公式 v[i]=1-u[i]; //λ1~λ10的求解公式 g[i]=(6/(x[i+1]-x[i-1]))*((y[i+1]-y[i])/(x[i+1]-x[i])-(y[i]-y[i-1])/(x[i]-x[i-1])); //g1~g10的求解公式 } g[0]=(6/(x[1]-x[0]))*((y[1]-y[0])/(x[1]-x[0])-0.75); //g0的求解公式 g[11]=(6/(x[11]-x[10]))*(18-(y[11]-y[10])/(x[11]-x[10])); //g11的求解公式 printf(u(1) to u(10):\n); for(i=1;i11;i++) { printf( %f , u[i]); j++; if(j%4==0) printf(\n); } printf(\n); j=0; printf(v(1) to v(10):\n); for(i=1;i11;i++) { printf( %f , v[i]); j++; if(j%4==0) printf(\n); } printf(\n); j=0; printf(g(0) to g(11):\n); for(i=0;i12;i++) { printf( %f , g[i]); j++; if(j%4==0) printf(\n); } for(i=1;i=12;i++) AllD=AllD*2; for(i=1;i=10;i++) { Mv=Mv*v[i]; Mu=Mu*u[i]; } //采用克拉默法则进行线性方程组的求解x[i]=D[i]/D double D[12]; D[0]=g[0]*AllD/2+Mv*g[11]; D[11]=g[11]*AllD/2+Mu*g[0]; for(i=1;i=10;i++) D[i]=g[i]*AllD/2; for(i=0;i=11;i++) M[i]=D[i]/AllD; for(i=1;i=11;i++) { X3[i]=(M[i]-M[i-1])/(6*(x[i]-x[i-1])); //计算想x^3系数的公式 X2[i]=(x[i]*M[i-1]-x[i-1]*M[i])/(2*(x[i]-x[i-1])); //计算想x^2系数的公式 X1[i]=(-x[i]*