chaotic_prediction2(test).m

来自「This program is to one-step EEG predicti」· M 代码 · 共 26 行

%% Testing the network
load W
load C
load f
m  = 7;
ERROR = [];
input_test = test_set;
a = 0;
b = 0;
M = ones(1,7);
y = 0;
for r = 1 : 299
    input = input_test(r-6:r);
    for j = 1:7
        for i = 1:7
            M(j) = M(j) * exp ( -0.5*((input(i) - C(i,j)) / W(i,j))^2 );
        end
        a = a + f(j) * M(j);
        b = b + M(j);
    end
    y = a / b;
    input_test(r+1) = y;
    ERROR = [ERROR ,0.5*(input_test (r+1)-x(r+1)).^2];
end