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

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

waveFilt = bandFilter(wave,3,60,S/1000);
wavefft = fft(waveFilt);
wavemag = abs(wavefft);
f = (0:length(wavefft)-1)'*S/length(wavefft);

for n=1:64
subplot(8,8,n); plot(f(5:0.2*S), wavemag(5:0.2*S, n+1))
end