function y = wavesets(sensors,sensors2,wl)

name = struct('strings', {{'COUNTER'; 'INTERPOLATED'; 'AF3'; 'F7'; 'F3'; 'FC5'; 'T7';'P7';

'O1'; 'O2'; 'P8';'T8'; 'FC6' ;'F4'; 'F8'; 'AF4'; 'RAW_CQ'; 'CQ_AF3'; 'CQ_F7'; 'CQ_F3';

'CQ_FC5'; 'CQ_T7'; 'CQ_P7'; 'CQ_O1'; 'CQ_O2'; 'CQ_P8'; 'CQ_T8'; 'CQ_FC6'; 'CQ_F4';

'CQ_F8'; 'CQ_AF4'; 'CQ_CMS'; 'CQ_DRL'; 'GYROX'; 'GYROY'}} );

N = wl;

if mod(N,2)==0

k=-N/2:N/2-1; % N even

else

k=-(N-1)/2:(N-1)/2; % N odd

end

for segmentid=1:floor(length(sensors)/wl);

segmentt=sensors(wl*(segmentid-1)+1:segmentid*wl,:);

segmentt2=sensors2(wl*(segmentid-1)+1:segmentid*wl,:);

for j=3:18%j=1:size(segmentt,2)

fourier = abs(fft(segmentt(:,j),128));

f(:,j-2) = fourier(64:128);

fourier = abs(fft(segmentt2(:,j),128));

f2(:,j-2) = fourier(64:128);

end

N = ceil(size(f,1)/2);

for i=3:18

g(i-2,1) = sum(f(8:12,i-2))/(12-8+1); %alpha

g(i-2,2) = sum(f(12:39,i-2))/(39-12+1); %beta

g(i-2,3) = sum(f(4:7,i-2))/(7-4+1); %theta

g(i-2,4) = sum(f(1:4,i-2))/(4-1+1); %delta

g2(i-2,1) = sum(f2(8:12,i-2))/(12-8+1); %alpha

g2(i-2,2) = sum(f2(12:39,i-2))/(39-12+1); %beta

g2(i-2,3) = sum(f2(4:7,i-2))/(7-4+1); %theta

g2(i-2,4) = sum(f2(1:4,i-2))/(4-1+1); %delta

end

x = [1 2 3 4];

xx= [1.5 2.5 3.5 4.5];

K = 0.5;

for i = 3: 18

j=i-2;

subplot(4,4,j)

title(name.strings(i));

bar1=bar(x, g(j,:)', 'FaceColor', 'b', 'EdgeColor', 'b');

set(bar1,'BarWidth',K/2);

%bar(g(j,:)'); figure(gcf)

hold on;

bar2=bar(xx, g2(j,:)', 'FaceColor', 'r', 'EdgeColor', 'r');

set(bar2,'BarWidth',K/2);

%bar(g2(j,:)', 'r' ); figure(gcf)

hold all;

pause(0.25);

end

pause(1);

hold off;

end