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

function [data] = PlotPowerSpec1( 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(1:100), wavemag(1:100, n+1)); axis([0 100 0 40]); 
end