-
Notifications
You must be signed in to change notification settings - Fork 0
/
apply_iir_filters.m
49 lines (39 loc) · 1.55 KB
/
apply_iir_filters.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
function y = apply_iir_filters(data,Fs,gain)
% return composit data
%low pass filter with cutoff frequency of 170Hz
filter_obj = iir_low_pass(170,Fs);
dummy_data = filter(filter_obj,data);
composite_data = dummy_data * gain(1);
%band pass filter 170 - 310 Hz
filter_obj = iir_band_pass(170,310,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(2);
%band pass filter 310 - 600 Hz
filter_obj = iir_band_pass(310,600,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(3);
%band pass filter 600 - 1k Hz
filter_obj = iir_band_pass(600,1000,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(4);
%band pass filter 1k - 3k Hz
filter_obj = iir_band_pass(1000,3000,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(5);
%band pass filter 3k - 6k Hz
filter_obj = iir_band_pass(3000,6000,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(6);
%band pass filter 6k - 12k Hz
filter_obj = iir_band_pass(6000,12000,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(7);
%band pass filter 12k - 14k Hz
filter_obj = iir_band_pass(12000,14000,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(8);
%band pass filter 14k - 16k Hz
filter_obj = iir_band_pass(14000,16000,Fs);
dummy_data = filter(filter_obj,data);
composite_data = composite_data + dummy_data * gain(9);
y = composite_data;