function MAPE=Predict(data,m,tau,P,lmd,MaxStep)
%單步預測

% disp('--------最大預測步數---------')
%  MaxStep=round(1/lmd)
% % MaxStep=5;
% % MaxStep=MaxStep+5;%避免預測太長控制在10步
deltaT=1;



for i=1:MaxStep
    i
    newdata=data(1:length(data)-MaxStep+i-1);
    pre_value(i) = FunctionChaosPredict(newdata,length(newdata),P,deltaT,tau,m,1);%調用AOLMM進行多步預報
end

%figure(1)
 plot(length(data)-MaxStep+1:length(data),data(length(data)-MaxStep+1:length(data)),'b-',length(data)-MaxStep+1:length(data),pre_value,'r-*');
legend('實測值','局域法預測值');
% %-----------求標準差-----------------
% delt=std(data(length(newdata)+1:length(data))-newdata1(length(newdata)+1:length(data)))

disp('-----------求RMSE----------------------')
old=data(length(data)-MaxStep+1:length(data))';
q=mean(old);
RMSE=norm(pre_value-old,2)/norm(old-q,2)
disp('-------------中誤差-----------------')
v=abs(pre_value-old);
MSE=sqrt(sum(v.^2)/length(old))

disp('----------MAPE平均絕對百分比誤差---------------')
for i=1:length(old)
    a2(i)=pre_value(i)-old(i);
    a(i)=abs((old(i)-pre_value(i))/old(i));
end

MAPE=sum(a)/length(old)*100

figure
plot(a2)
b=max(a2)-min(a2);
axis([1 length(a2) min(a2)-b max(a2)+b])
grid;xlabel('預測步數');ylabel('誤差')

xuhao=(length(data)-MaxStep+1:length(data));
jieguo=[xuhao;pre_value;old;a2;a]';