predict.m

来自「本文是通过文中方法来求最小嵌入维程序.子程序,相重构程序.」· M 代码 · 共 48 行

 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]';