插值和拟合算法的例子.doc

插值和拟合算法的例子：//---------------------------------------------------------------------------//用最小二乘法拟合曲线y=a0+a1*x+a2*x^2+a3*x^3+...+an*x^n;//x,y为自变量和因变量数组;n为为系数个数;m为数据的组数;sum用于存和的数组;//l为sum数组边界;r为系数数组,也用于存x^n*y的和;a用于存放正则方程组的系数方阵;//返回值:0,正常;-1,正则方程系数方阵异常;-2,数据太少,无法拟合;int LeastSquarFit(double x[],double y[],int n,int m,double sum[],? ? ? ? ? ? ? ? double an[],double dMatrix[][5]){?int i,j,k,k1,rtn=0;?double w,dx,dy,xx,eps;? ? ? ? if(mn) ? ? ? ? //数据的组数小于所求系数的个数,无法拟合曲线!? ? ? ? ? ? ? ? return -2;? ? ? ? sum[0]=m;? ? ? ? for(j=0;jn;j++)? ? ? ? ? ? ? ? an[j]=0.0;? ? ? ? for(j=1;j2*n-1;j++)? ? ? ? ? ? ? ? sum[j]=0.0;? ? ? ? for(i=0;im;i++)? ? ? ? {? ? ? ? ? ? ? ? dx=1.0;? ? ? ? ? ? ? ? dy=y[i];? ? ? ? ? ? ? ? xx=x[i];? ? ? ? ? ? ? ? for(j=1;j2*n-1;j++)? ? ? ? ? ? ? ? {? ? ? ? ? ? ? ? ? ? ? ? dx*=xx;? ? ? ? ? ? ? ? ? ? ? ? sum[j]+=dx; ? ? //x的j次方的和? ? ? ? ? ? ? ? }? ? ? ? ? ? ? ? for(j=0;jn;j++)? ? ? ? ? ? ? ? {? ? ? ? ? ? ? ? ? ? ? ? an[j]+=dy; ? ? ? //y与x的j次方的乘积的和? ? ? ? ? ? ? ? ? ? ? ? dy*=xx;? ? ? ? ? ? ? ? }? ? ? ? }? ? ? ? k1=0;? ? ? ? for(j=0;jn;j++)? ? ? ? {? ? ? ? ? ? ? ? k=k1;? ? ? ? ? ? ? ? for(i=0;in;i++)? ? ? ? ? ? ? ? {? ? ? ? ? ? ? ? ? ? ? ? dMatrix[i][j]=sum[k]; //正则方程的系数为各次x的和? ? ? ? ? ? ? ? ? ? ? ? k++;? ? ? ? ? ? ? ? }? ? ? ? ? ? ? ? k1++;? ? ? ? }? ? ? ? int nans[5];? ? ? ? rtn=SolvLinEqua(dMatrix,n,an,nans,1e-10);? ? ? ? return rtn;}//---------------------------------------------------------------------------//用高斯消除法解n阶一次线性方程组,先化方阵为对角阵,然后求解//Matrix为系数方阵；n为方阵的阶数；an作为输入参数时为方程组右端常数项，//an作为返回值时为方程组的解；min为误差限，最小为1e-8；nans为是计算中要用的中间数组//函数返回值：0为正常；-1，系数矩阵异常int SolvLinEqua(double dMatrix[][5],int n,double an[],int nans[],double min){? ? ? ? int i,k,j,m,nstop;? ? ? ? double epss,bigst,pivot,v,w;? ? ? ? nstop=0; ? ? ? ?//返回值预置为0-正常? ? ? ? if(min1e-8) ? ?//误差限至少为1e-8? ? ? ? ? ? ? ? min=1e-5;? ? ? ? epss=min*min;? ? ? ? for(i=0;in;i++)? ? ? ? ? ? ? ? nans[i]=-1; ? ? //nans用于记录列标，预置为-1? ? ? ? for(k=0;kn;k+