%filters data and plots power spectrum for 64 electrode array
%requires raw data and sampling frequency

function [data] = PlotRaw4( wave, S );

waveFilt = BudapestbandPass(wave(:,2:65),1000,5,200);

Lrotate=ones(8,8);
for m=1:8
    for n=1:8
        Lrotate(m,n)=n*8-m+1;
    end
end
Lrotate=reshape(Lrotate,1,64);

Rrotate=ones(8,8);
for n=1:8
    for m=1:8
        Rrotate(m,n)=8*(8-m)+n
    end
end
Rrotate=reshape(Rrotate,1,64);

figure
plot(wave(:,1), waveFilt(:,27)); axis([1200 1700 -0.1 0.1]); 
figure
for n=1:64
subplot(8,8,Rrotate(n)); plot(wave(:,1), waveFilt(:,n)); 
axis([1200 1700 -0.1 0.1]); axis off;
end
