时间序列arima matlab程序

clc
clear
y=[];%自己输入
Data=y;%共？个数据
SourceData=Data（1:6001）; %前600个训练集
step=120;%后120个测试
TempData=SourceData;
TempData=detrend（TempData）;%去趋势线
TrendData=SourceData-TempData;%趋势函数

H=adftest（TempData）;
difftime=0;
SaveDiffData=[];
while ~H
SaveDiffData=[SaveDiffDataTempData（11）];
TempData=diff（TempData）;%差分，平稳化时间序列
difftime=difftime+1;%差分次数
H=adftest（TempData）;%adf检验，判断时间序列是否平稳化
end

u = iddata（TempData）;
test = [];
for p =