clear all
clc

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;%保留原始數據


th=arx(YOrigin,7)
fd=predict(YOrigin,th)

preErrArraylist=[];%誤差列表
preErrSum=0;%誤差和
for i=1:1:1063
    preErr=abs((YOrigin(i)-fd(i))/YOrigin(i));
    preErrArraylist=[preErrArraylist,preErr];
    preErrSum = preErrSum + preErr;
end
preAverageError=preErrSum/length(preErrArraylist) %平均預測誤差
i=1:1:1063;
plot(i,preErrArraylist(i))