原创 指数平滑法预测数据(转贴)

2008-4-13 22:37 2000 2 2 分类: 工业电子
/************************************************************************
* 本算法用指数平滑法预测数据
* 输入: 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;
}
}

文章评论0条评论)

登录后参与讨论
相关推荐阅读
用户907182 2013-11-17 16:56
基于Matlab的IIR数字滤波器设计方法比较及应用
  0 引言   数字滤波器(Digital Filter)是指输入、输出都是离散时间信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器在数字信号处理中起着...
用户907182 2012-06-09 08:04
评论:@信号链杂谈's Blog 博客中提到的“信号链基础知识# 63:利用 SPICE 设计 TEC 温度环路 ...”
TEC温度环路PID控制...
用户907182 2011-02-01 23:17
用DSP实现增量式光电编码器的细分
用DSP实现增量式光电编码器的细分 技术分类: 微处理器与DSP  | 2011-01-04 21ic   1 引言   目前,各类伺服驱动器及其应用中广泛采用光栅装置作为速度测量、位置测量的敏感元...
用户907182 2011-02-01 23:16
基于MSP430F149单片机的光电编码器位置检测系统设计
基于MSP430F149单片机的光电编码器位置检测系统设计 技术分类: 测试与测量  | 2011-01-31 微计算机信息 李红果   在实际的工业位置控制领域中,为了提高控制精度,准确地对控制对...
用户907182 2011-02-01 23:11
业界首款可直接安装在电机上的定制化电机控制模块
业界首款可直接安装在电机上的定制化电机控制模块 技术分类: 工业电子  | 2010-12-03 : EDN China    智能功率器件及解决方案供应商意法半导体发布业界首款整合工业级以太网接口...
用户907182 2011-02-01 23:05
MEMS惯性传感器优势解析THELMA制程和低成本封装方法
MEMS惯性传感器优势解析:THELMA制程和低成本封装方法 技术分类: 模拟设计  | 2010-11-29 意法半导体 Benedetto Vigna   意法半导体公司推出一系列惯性传感器,极...
我要评论
0
2
1
2
3
4
5
6
7
8
9
0
关闭 站长推荐上一条 /3 下一条