% generate a perfectly reconstructive filterbank
fs=48000;
octFiltBank = octaveFilterBank('1/3 octave','FrequencyRange',[25,20000],'FilterOrder',6,'SampleRate',fs)
% fvt = fvtool(octFiltBank,"NFFT",2^16);
% set(fvt,FrequencyScale="log")
% zoom(fvt,[.01 24 -20 1])
% x=zeros(65536,1);
% x(4096,1)=1;
[x fs]=audioread('Eq-Filter_1_3.wav');
y = octFiltBank(x);
audiowrite('Filterbank_1_3Oct_eq.wav',y,fs,'BitsPerSample',32);