clc
clear
y=[];%自己輸入
Data=y;%共？個(gè)數(shù)據(jù)
SourceData=Data（1:6001）;?%前600個(gè)訓(xùn)練集
step=120;%后120個(gè)測試
TempData=SourceData;
TempData=detrend（TempData）;%去趨勢線
TrendData=SourceData-TempData;%趨勢函數(shù)

H=adftest（TempData）;
difftime=0;
SaveDiffData=[];
while?~H
SaveDiffData=[SaveDiffDataTempData（11）];
TempData=diff（TempData）;%差分，平穩(wěn)化時(shí)間序列
difftime=difftime+1;%差分次數(shù)
H=adftest（TempData）;%adf檢驗(yàn)，判斷時(shí)間序列是否平穩(wěn)化
end

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