【数值分析】相关数值分析多种算法代码
整理一些相关的数值分析的代码，共享给急切需要同行们！希望能在您能获多获少都会有所收获.>_<呵呵.
//***********************************************************************// 本算法用指数平滑法预测数据// 输入: k--平滑周期// n--原始数据个数// m--预测步数// alfa--加权系数// x--指向原始数据数组指针// 输出: s1--返回值为指向一次平滑结果数组指针// s2--返回值为指向二次指数平滑结果数组指针// s3--返回值为指向三次指数平滑结果数组指针// xx--返回值为指向预测结果数组指针//***********************************************************************void phyc(int k,int n,int m,double alfa,double x[N_MAX], double s1[N_MAX],double s2[N_MAX],double s3[N_MAX],double xx[N_MAX]){ double a,b,c,beta; int i; s1[k-1]=0; for(i=0;i<k;k++) s1[k-1]+=x; s1[k-1]/=k; for(i=k;i<=n;i++) s1=alfa*x+(1-alfa)*s1[i-1]; s2[2*k-2]=0; for(i=k-1;i<2*k-1;i++) s2[2*k-2]+=s1; s2[2*k-2]/=k; for(i=2*k-1;i<=n;i++) s2=alfa*s1+(1-alfa)*s2[i-1]; s3[3*k-3]=0; for(i=2*k-2;i<3*k-2;i++) s3[3*k-3]+=s2; s3[3*k-3]/=k; for(i=3*k-2;i<=n;i++) s3=alfa*s2+(1-alfa)*s3[i-1]; beta=alfa/(2*(1-alfa)*(1-alfa)); for(i=3*k-3;i<=n;i++){ a=3*s1-3*s2+s3; b=beta*((6-5*alfa)*s1-2*(5-4*alfa)*s2+(4-3*alfa)*s3); c=beta*alfa*(s1-2*s2+s3); xx=a+b*m+c*m*m; }}======================================================================转载请注明出处：http://blog.csdn.net/songzitea/article/details/13627193======================================================================
