clear
close all

in = load('Untitled_002.mat');

fs = in.Result.AcqSetup.Fs;
resp = in.Result.ImpulseResponse(:,2);
[b,a] = butter(5,2e3/(fs/2));

resp_out = filter(b,a,resp);

figure
freq = (0:length(resp_out)-1)/length(resp_out)*fs;
semilogx(freq*2,20*log10(abs(fft(resp_out))));
grid on
xlim([1,5e3])