ita_filterRir.m 400 Bytes
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
function state = ita_filterRir(params,state)
% Noam Shabtai
% Institution of Technical Acoustics
% RWTH Aachen
% nsh@akustik.rwth-aachen.de
% 24.7.2014

if ~params.stages.rir.filter return; end

h = state.rir_signal;
N = params.filter.order;
band = params.filter.normalized_band;

[numerator, denominator] = butter(N, band);
hf = filter(numerator, denominator, h);

state.rir_signal_filtered = hf;