clear all;
close all;
clc;

load('eegData');
amountOfDataToPlot = 100;
referenceChannel1 = 60;
referenceChannel2 = 63;

% get some data to plot from before the filter
preReferenceData = eegData(1,1:amountOfDataToPlot);

% reference the data using two channels
referenceChannel = (eegData(referenceChannel1,:) + eegData(referenceChannel2,:))/2;
% subtract the reference channel from the data
for channelCounter = 1:size(eegData,1)
    eegData(channelCounter,:) = eegData(channelCounter,:) - referenceChannel;
end

% get some data to plot post filter
postReferenceData = eegData(1,1:amountOfDataToPlot);

% plot some stuff
plot([1:1:amountOfDataToPlot],preReferenceData,'linewidth',3);
hold on;
plot([1:1:amountOfDataToPlot],postReferenceData,'linewidth',3);
hold off;
title('Pre versus Post Reference Data');