新建 文本文档.txt

指数平滑法预测数据

/************************************************************************ 
* 本算法用指数平滑法预测数据 
* 输入: 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; 
} 
}