回声状态网络ESN

dataset = csvread('D:\matlab2016a\data.csv');
x = dataset(:,1:3);
y = dataset(:,4);

x_train = x(1:145,:);
y_train = y(1:145);
x_test = x(146:end,:);
y_test = y(146:end);

[samples,inputnum] = size(x_train);
hiddennum = 10;
outputnum = 1;

w_in = rand(hiddennum,inputnum);
w_state = rand(hiddennum,hiddennum);
w_back = rand(hiddennum,outputnum);

s(1,:) = zeros(1,hiddennum);
for i = 1:samples
    s(i+1,:) = tanh(w_in*x_train(i,:)'+w_state*s(i,:)'+w_back*y_train(i));
end
s = s(2:end,:);
w_out = (inv(s'*s)*s'*y_train)';

for i = 1:samples
    predict_train(i) = w_out*s(i,:)';
end

train_RMSE = sqrt((sum(predict_train-