function [LPBand, LMBand, HMBand, HPBand] = audiofilt(input, HLowPass, HLowMid, HHighMid, HHighPass)

LP = filter(HLowPass, input); %Filtering original signal into different freq bands
LP1 = flipud(LP);
LP2 = filter(HLowPass, LP1);
LPBand = flipud(LP2); %Flipping values and running through filter backwards to reduce phase when summed

LM = filter(HLowMid, input);%Filtering original signal into different freq bands
LM1 = flipud(LM);
LM2 = filter(HLowMid, LM1);
LMBand = flipud(LM2);%Flipping values and running through filter backwards to reduce phase when summed

HM = filter(HHighMid, input); %Filtering original signal into different freq bands
HM1 = flipud(HM);
HM2 = filter(HHighMid, HM1);
HMBand = flipud(HM2); %Flipping values and running through filter backwards to reduce phase when summed

HP = filter(HHighPass, input); %Filtering original signal into different freq bands
HP1 = flipud(HP);
HP2 = filter(HHighPass, HP1);
HPBand = flipud(HP2); %Flipping values and running through filter backwards to reduce phase when summed