matlab写二次指数平滑

function y=Expon(X,a,b)
    n=length(X);
    y=zeros(n);
    t=zeros(n);
    y(1)=X(1);%初始化初值
    t(1)=X(1);
    for i=2:n
        y(i)=a*X(i)+(1-a)*(y(i-1)+t(i-1));
        t(i)=b*(y(i)-y(i-1))+(1-b)*t(i-1);
    end

end