www.gusucode.com > MATLAB2008应用程序接口编程技术源码程序 > MATLAB2008应用程序接口编程技术源码程序/code/第7章/7.5/示例2/示例2.m

    [m,n]=size(a)
alf=0.2;   % 设定平滑系数的数值
for j=1:2
s(1,j)=a(1,1);
end
for i=2:m
for j=1:2
if j==1
s(i,j)=alf*a(i,1)+(1-alf)*s(i-1,j);
else
s(i,j)=alf*s(i,j-1)+(1-alf)*s(i-1,j);
              end
end
end
% 计算预测方程的系数 a(t),b(t)
temp=alf/(1-alf);  ％将循环式中每次运算的常数提取出来计算，以减少计算量
for i=1:m
at(i,1)=2*s(i,1)-s(i,2);
bt(i,1)=temp*(s(i,1)-s(i,2));
end
%得出的预测模型为 Y＝at(m,1)+bt(m,1)*T
%求初步预测值
for i=1:4
y1(i)=at(m,1)+bt(m,1)*i;
end
y1