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?=