%%%%%%%%%%%%%%%%%% BASE
filename = 's14.xls';
A = xlsread(filename);
B_protected = xlsread(filename);
%%%%%%%%%%%%%%%%% REMOVE NaN
x=0;
for i = 2996:-1:1
    z = 0;
%%%%%%%% CHANGE J HERE(columns you wanna work on)
    for j = 15:1:22
       if(isnan(A(i,j)))
            z = z+1;
       end
    end
    if(z>=1)
        A(i,:)=[];
        x=x+1;
    end
end
%%%%%%%%%%%%%%%%% FILL NaN
%{
for i = 2996:-1:1
    for j = 1:1:23
       if(isnan(A(i,j)))
           if((i - 1 > 0))
               if(~isnan(A(i - 1,j)))
                A(i,j) = (A(i-1, j) + A(i+1, j))/2;
               end
           end
       end
    end
end
%}
%%%%%%%%%%%%%%%%% SELECT COLUMNS FIRST DATA SET
range = size(A,1);
%%%%%%%%%%%COLUMN FROM EXCEL -2
firstRows = 13;
for i = 1:1:range
    for j = 1 : 1 : 4
    dataSet1(i,j) = A(i,firstRows + j);
    end
end
%%%%%%%%%%%%%%%%% SELECT COLUMNS SECOND DATA SET
%%%%%%%%%%%COLUMN FROM EXCEL -2
secondRows = 17;
for i = 1:1:range
    for j = 1 : 1 : 4
    dataSet2(i,j) = A(i,secondRows + j);
    end
end
%%%%%%%%%%%%%%%%%
for i = 1:1:range
    outputData(i, 1) = (dataSet1(i,1) + dataSet1(i,2) + dataSet1(i,3) + dataSet1(i,4)) / 4;
    outputData(i, 2) = (dataSet2(i,1) + dataSet2(i,2) + dataSet2(i,3) + dataSet2(i,4)) / 4;
end
sortedData = sortrows(outputData, 1);
plot(sortedData);
grid on;
%%%%%%%%%%%%%%%%%
title('Zestawienie częstotliwości skoków od napięcia elektrofiltru');
xlabel('Liczba próbek');
ylabel('Wartość napięcia');
legend({'Napięcie elektrofiltru [kV]','Częstotliwość przeskoków'},'Location','northeast')
print(gcf,'zad1.png','-dpng','-r1024');