N = 15;
M = 4;
O_g = [2, 3, 4, 5];
O_b = [3, 5, 11, 13];

%input = randi([0,1],[1,N]);
output_g = zeros(1,N);
output_b = zeros(1,N);

output_g(2:end) = input(1:end-1);
output_b(2:end) = input(1:end-1);

for n=1:1:2^N-1
    output_g(2:end) = output_g(1:end-1);
    output_b(2:end) = output_b(1:end-1);
    for i=1:1:M
        output_g(1) = bitxor(output_g(1), output_g(O_g(i)));
        output_b(1) = bitxor(output_b(1), output_b(O_b(i)));
    end
end

figure;

subplot(3, 2, 1);
stairs(input);
title('Wejscie')
ylim([-0.5 1.5]);

autocorr_g = xcorr(output_g);
subplot(3, 2, 3);
stairs(output_g);
title('Dobry')
ylim([-0.5 1.5]);
subplot(3, 2, 4);
plot(autocorr_g);
title('Dobry - autokorelacja')

autocorr_b = xcorr(output_b);
subplot(3, 2, 5);
stairs(output_b);
title('Słaby')
ylim([-0.5 1.5]);
subplot(3, 2, 6);
plot(autocorr_b);
title('Słaby - autokorelacja')