function [filterBank, centralFreqs] = createFilterBank(sr, numBands, filterOrder)

% Define low and high frequencies
freqLow = 100;
freqHigh = 7999;
logDifference = (log10(freqHigh) - log10(freqLow))/numBands;

for i = 1:numBands
    %logDifference
    %freqLow
    freqHigh = 10^(log10(freqLow) + logDifference);
    centralFreqs(i) = 10^((log10(freqHigh) + log10(freqLow))/2);
    fInit = fdesign.bandpass('N,F3dB1,F3dB2', filterOrder, freqLow/(sr/2), freqHigh/(sr/2));
    filter = design(fInit, 'butter');
    filterBank(i) = filter;
    freqLow = freqHigh;
end