clc;

%figure, plot(sig(1:1000,1)) % wykres sygnalu zaszumionego

%figure, plot(sig(1:1000,2)) % wykres sygnalu bodzca

r=find(sig(:,2)); % wyciagniecie wszystkich impulsow z sygnalu bodzca

p=r(2:end)-r(1:end-1);

hist(p); % rozklad sygnalu bodźca

dlugosc_usredniania = min(p);

dane = zeros(length(r),dlugosc_usredniania);


j=1;
for i = r.'
    dane(j,:) = sig(i:i+dlugosc_usredniania-1,1).';
    j=j+1;
end
%
%
%srednia = mean(dane(801:900,:));
%odchylenie_std = std(dane(801:900,:))/sqrt(length(dane(801:900,:)));
%
%figure, plot(1:dlugosc_usredniania, srednia, 'r', 1:dlugosc_usredniania, srednia+odchylenie_std, 'b', 1:dlugosc_usredniania, srednia-odchylenie_std, 'b')
%figure, plot(odchylenie_std)
%


srednia_1 = mean(dane(1:100,:));
odchylenie_std_1 = std(dane(1:100,:))/sqrt(length(dane(1:100,:)));


srednia_2 = mean(dane(501:600,:));
odchylenie_std_2 = std(dane(501:600,:))/sqrt(length(dane(501:600,:)));

srednia_3 = mean(dane(801:900,:));
odchylenie_std_3 = std(dane(801:900,:))/sqrt(length(dane(801:900,:)));

figure, plot(1:dlugosc_usredniania, srednia_1, '-r', 1:dlugosc_usredniania, srednia_1+odchylenie_std_1, '.r', 1:dlugosc_usredniania, srednia_1-odchylenie_std_1, '.r',1:dlugosc_usredniania, srednia_2, '-b', 1:dlugosc_usredniania, srednia_2+odchylenie_std_2, '.b', 1:dlugosc_usredniania, srednia_2-odchylenie_std_2, '.b',1:dlugosc_usredniania, srednia_3, '-g', 1:dlugosc_usredniania, srednia_3+odchylenie_std_3, '.g', 1:dlugosc_usredniania, srednia_3-odchylenie_std_3, '.g')