曲线拟合 给定10个点 拟合曲线 复制即用 c++.docVIP

#include iostream #include fstream #include string #includeiomanip.h #define F1(t) (1+(t)*(t)*(2*(t)-3)) #define F2(t) ((t)*(t)*(-2*(t)+3)) #define F3(t) ((t)*(1+(t)*((t)-2))) #define F4(t) ((t)*(t)*((t)-1)) using namespace std; void tangent(double *y,double *dy,int n) //计算各个节点的切矢，采用自由端点条件 //追赶法求解三对角方程组 { int i; double *a,*b,*c,*d; a=new double[n]; b=new double[n]; c=new double[n]; d=new double[n]; //追赶法计算系数 for(i=1;i n;i++) a[i]=1; b[0]=2; b[n-1]=2; for(i=1;i n-1;i++) b[i]=4; for(i=0;i n-1;i++) c[i]=1; d[0]=3*(y[1]-y[0]); d[n-1]=3*(y[n-1]-y[n-2]); for(i=1;i n-1;i++) d[i]=3*(y[i+1]-y[i-1]); //追赶法回代 double *L,*U; L=new double[n]; U=new double[n]; U[0]=b[0]; for(i=1;i n;i++) { L[i]=a[i]/U[i-1]; U[i]=b[i]-L[i]*c[i-1]; } double *yy; yy=new double[n]; yy[0]=d[0]; for(i=1;i n;i++) yy[i]=d[i]-L[i]*yy[i-1]; dy[n-1]=yy[n-1]/U[n-1]; for(i=n-2;i=0;i--) dy[i]=(yy[i]-c[i]*dy[i+1])/U[i]; //delete a,b,c,d,L,U,yy; } double spline(double *x,double *y,double *dy,double xa,int n) //三次Hermite插值 //数据点x[n],y[n],切矢dy[n] //返回xa处的插值数值 //要求x严格递增，x[0] =xa =x[n-1] { int i; for(i=0;i n-1;i++) if(xa =x[i+1]) break; double t = (xa-x[i])/(x[i+1]-x[i]); return F1(t)*y[i]+F2(t)*y[i+1]+F3(t)*dy[i]+F4(t)*dy[i+1]; } int main() { double s_ua,xa_ua;//A类型循环变量 double s_ub,xa_ub;//B类型循环变量 double s_uc,xa_uc; //C类型循环变量 double s_a_result,xa_a_result;// //double s_b_result,xa_b_result;// double average_s,average_xa; double dy_ua[4]; double dy_ub[4]; double dy_uc[4]; double x_ua[4]={0.000,0,0,0}; double y_ua[4]={0.0000,0,0,0}; double x_ub[4]={0.000,0,0,0}; double y_ub[4]={0.0000,0,0,0}; double x_uc[4]={0.000,0,0,0}; double y_uc[4]={0.0000,0,0,0}; //DNA //double X[1