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

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

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

%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(f(1:100), wavemag(1:100,5)); axis([0 100 0 20]); xlabel ('Frequency (Hz)'); ylabel ('Magnitude');

figure
for n=2:8
subplot(4,2,n); plot(f(1:100), wavemag(1:100, n)); axis([0 100 0 20]); %axis off;
end