data = csvread('D:\article\demo\seasonal\stl.csv');
sample = size(data,1);
timespan = 6;
imfindex = 3;
for i = 1:(sample-timespan)
x(i,:) = data(i:i+timespan-1,imfindex)';
end
y = data(timespan+1:end,imfindex);
predictspan = 7;
dividindex = size(x,1) - predictspan + 1;
x_train = x(1:dividindex-1,:);
y_train = y(1:dividindex-1,:);
x_test = x(dividindex:end,:);
y_test = y(dividindex:end,:);
[train_input,inputps] = mapminmax(x_train');
[train_output,outputps] = mapminmax(y_train');
inputnum = size(x,2);
hiddennum = 12;
len = inputnum*hiddennum+hiddennum;
%PSO参数
c1 = 2;
c2 = 2;
sizepop = 20;
k = 100;
w = 0.9
param = rand(sizepop,len)