alexbrandmeyer@643
|
1 //
|
alexbrandmeyer@643
|
2 // agc.h
|
alexbrandmeyer@643
|
3 // CARFAC Open Source C++ Library
|
alexbrandmeyer@643
|
4 //
|
alexbrandmeyer@643
|
5 // Created by Alex Brandmeyer on 5/30/13.
|
alexbrandmeyer@643
|
6 //
|
alexbrandmeyer@643
|
7 // This C++ file is part of an implementation of Lyon's cochlear model:
|
alexbrandmeyer@643
|
8 // "Cascade of Asymmetric Resonators with Fast-Acting Compression"
|
alexbrandmeyer@643
|
9 // to supplement Lyon's upcoming book "Human and Machine Hearing"
|
alexbrandmeyer@643
|
10 //
|
alexbrandmeyer@643
|
11 // Licensed under the Apache License, Version 2.0 (the "License");
|
alexbrandmeyer@643
|
12 // you may not use this file except in compliance with the License.
|
alexbrandmeyer@643
|
13 // You may obtain a copy of the License at
|
alexbrandmeyer@643
|
14 //
|
alexbrandmeyer@643
|
15 // http://www.apache.org/licenses/LICENSE-2.0
|
alexbrandmeyer@643
|
16 //
|
alexbrandmeyer@643
|
17 // Unless required by applicable law or agreed to in writing, software
|
alexbrandmeyer@643
|
18 // distributed under the License is distributed on an "AS IS" BASIS,
|
alexbrandmeyer@643
|
19 // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
alexbrandmeyer@643
|
20 // See the License for the specific language governing permissions and
|
alexbrandmeyer@643
|
21 // limitations under the License.
|
alexbrandmeyer@643
|
22
|
alexbrandmeyer@643
|
23 #ifndef CARFAC_AGC_H
|
alexbrandmeyer@643
|
24 #define CARFAC_AGC_H
|
alexbrandmeyer@643
|
25
|
alexbrandmeyer@643
|
26 #include <vector>
|
alexbrandmeyer@643
|
27 #include "common.h"
|
alexbrandmeyer@643
|
28
|
alexbrandmeyer@643
|
29 struct AGCParams {
|
alexbrandmeyer@643
|
30 AGCParams() {
|
alexbrandmeyer@643
|
31 num_stages = 4;
|
alexbrandmeyer@643
|
32 agc_stage_gain = 2.0;
|
alexbrandmeyer@643
|
33 time_constants.resize(num_stages);
|
alexbrandmeyer@643
|
34 agc1_scales.resize(num_stages);
|
alexbrandmeyer@643
|
35 agc2_scales.resize(num_stages);
|
alexbrandmeyer@643
|
36 agc1_scales[0] = 1.0;
|
alexbrandmeyer@643
|
37 agc2_scales[0] = 1.65;
|
alexbrandmeyer@643
|
38 time_constants[0] = 0.002;
|
alexbrandmeyer@643
|
39 for (int i = 1; i < num_stages; ++i) {
|
alexbrandmeyer@643
|
40 agc1_scales[i] = agc1_scales[i - 1] * sqrt(2.0);
|
alexbrandmeyer@643
|
41 agc2_scales[i] = agc2_scales[i - 1] * sqrt(2.0);
|
alexbrandmeyer@643
|
42 time_constants[i] = time_constants[i - 1] * 4.0;
|
alexbrandmeyer@643
|
43 }
|
alexbrandmeyer@643
|
44 decimation = {8, 2, 2, 2};
|
alexbrandmeyer@643
|
45 agc_mix_coeff = 0.5;
|
alexbrandmeyer@643
|
46 }
|
alexbrandmeyer@643
|
47 int num_stages;
|
alexbrandmeyer@643
|
48 FPType agc_stage_gain;
|
alexbrandmeyer@643
|
49 FPType agc_mix_coeff;
|
alexbrandmeyer@643
|
50 std::vector<FPType> time_constants;
|
alexbrandmeyer@643
|
51 std::vector<int> decimation;
|
alexbrandmeyer@643
|
52 std::vector<FPType> agc1_scales;
|
alexbrandmeyer@643
|
53 std::vector<FPType> agc2_scales;
|
alexbrandmeyer@643
|
54 };
|
alexbrandmeyer@643
|
55
|
alexbrandmeyer@643
|
56 struct AGCCoeffs {
|
alexbrandmeyer@643
|
57 int num_agc_stages;
|
alexbrandmeyer@643
|
58 FPType agc_stage_gain;
|
alexbrandmeyer@643
|
59 FPType agc_epsilon;
|
alexbrandmeyer@643
|
60 int decimation;
|
alexbrandmeyer@643
|
61 FPType agc_pole_z1;
|
alexbrandmeyer@643
|
62 FPType agc_pole_z2;
|
alexbrandmeyer@643
|
63 int agc_spatial_iterations;
|
alexbrandmeyer@643
|
64 FPType agc_spatial_fir_left;
|
alexbrandmeyer@643
|
65 FPType agc_spatial_fir_mid;
|
alexbrandmeyer@643
|
66 FPType agc_spatial_fir_right;
|
alexbrandmeyer@643
|
67 int agc_spatial_n_taps;
|
alexbrandmeyer@643
|
68 FPType agc_mix_coeffs;
|
alexbrandmeyer@643
|
69 FPType agc_gain;
|
alexbrandmeyer@643
|
70 FPType detect_scale;
|
alexbrandmeyer@643
|
71 FPType decim;
|
alexbrandmeyer@643
|
72 };
|
alexbrandmeyer@643
|
73
|
alexbrandmeyer@643
|
74 struct AGCState {
|
alexbrandmeyer@643
|
75 ArrayX agc_memory;
|
alexbrandmeyer@643
|
76 ArrayX input_accum;
|
alexbrandmeyer@643
|
77 int decim_phase;
|
alexbrandmeyer@643
|
78 };
|
alexbrandmeyer@643
|
79
|
alexbrandmeyer@643
|
80 #endif // CARFAC_AGC_H |