clear
close all

measurement = open('Fisher_H_IR.mat');

fs = measurement.fs_H;
IR = measurement.IR_H;
IR0 = IR(1,:)';
t = (0 : size(IR0,1)-1)'/fs;
f = (0 : size(IR0,1)-1)'/size(IR0,1)*fs;

figure
subplot(1,2,1);
plot(t, IR0);
xlabel('t -> [s]');
title('Mért impulzusválasz')

subplot(1,2,2);
semilogx(f, 20*log10(abs(fft(IR0))));
xlim([20, 20e3])
xlabel('f -> [Hz]');
ylabel('|H(f)| -> [dB]');
title('Mért átviteli karakterisztika');