chaotic_prediction_testing.m

来自「This program is prepared as an one-step 」· M 代码 · 共 25 行

%% Testing the network
ERROR = [];
input_test = x;
a = 0;
b = 0;
M = ones(1,3);
y = 0;
for r = 1200 : 1499
    input = input_test(r-2:r);
    for j = 1:3
        for i = 1:3
            M(j) = M(j) * exp ( -((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
figure
plot(input_test(1200:1500)); hold on
plot(x(1200:1500),'r');
% plot(ERROR, 'DisplayName', 'ERROR', 'YDataSource', 'ERROR'); figure(gcf)
%