huise.m

基于混沌扰动算法来解决pso算法早熟问题

function [GM f]=huise(data,H)
T=length(data); 
X0=data;

for i=2:T                      
XR1(i)=(X0(1)-H(2)/H(1))*exp(-1*H(1)*(i-1))+H(2)/H(1);
end
K(1)=data(1);
XR1(1)=data(1);

for i=2:T                     %还原计算出预测值
        K(i)=XR1(i)-XR1(i-1);
end

GM=minf(K,data);
f=K;

%N为预测的年数