季节性差分处理matlab.txt

用matlab编写季节性差分处理程序！！ ！！ ！！ ！！ ！1

test=datamain(1);
len = test(2);%原始数据个数
i=[1:1:len]';
X=[(i.^0),i,(i.^2)];
Y=[];
for tm=1:1:len
    test=datamain(tm);
    Y=[Y;test(1)];%读取原始数据
end
YOrigin=Y;%保留原始数据
subplot(3,1,1)
plot(i,YOrigin(i))
title('实际数据走势图')

%%%%%%%%周期性一阶差分处理（保留后面的）%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
s = 7; %周期数据
XS1 = [] ;
YS1 = [] ;
for j= s+1:1:len
    xj = i(j);
    XS1 = [XS1 xj];   %处理后的横坐标
    yj = Y(j) - Y(j-s) ;
    YS1 = [YS1 yj];   %处理后的纵坐标
end 
subplot(3,1,2)
plot(XS1,YS1)
title('周期性一阶差分处理后走势图')

%%%%%%%%周期性2阶差分处理（保留后面的）%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
s = 7; %周期数据
XS2 = [] ;
YS2 = [] ;
for j= s+1:1:length(XS1)
    xj = XS1(j);
    XS2 = [XS2 xj];   %处理后的横坐标
    yj = YS1(j) - YS1(j-s) ;
    YS2 = [YS2 yj];   %处理后的纵坐标
end 
subplot(3,1,3)
plot(XS2,YS2)
title('周期性2阶差分处理后走势图')