clear;
close all;
fname = 'zene_frost.wav';
raw_name = 'zene_raw.wav';

p_raw = audioread(raw_name);
[p_frost, fs] = audioread(fname);
n = size(p_frost,1);
t = (0 : n-1)' / fs;

figure;
plot(t, p_frost, t, p_raw(:,1), t, mean(p_raw, 2));
legend('frost', 'input', 'mean');

fc = 500;
[b,a]  = butter(3, fc/(fs/2), 'high');

p_raw = filtfilt(b, a, p_raw);
p_frost = filtfilt(b, a, p_frost);

sel = t >= 30 & t < 35;
figure;
plot(t(sel), p_frost(sel,1), t(sel), p_raw(sel,1));
legend('frost', 'input');
