flatmax@592: // Copyright 2013 Matt R. Flax <flatmax\@> All Rights Reserved.
flatmax@593: // Author Matt Flax <flatmax@>
flatmax@592: //
flatmax@592: // This C++ file is part of an implementation of Lyon's cochlear model:
flatmax@592: // "Cascade of Asymmetric Resonators with Fast-Acting Compression"
flatmax@592: // to supplement Lyon's upcoming book "Human and Machine Hearing"
flatmax@592: //
flatmax@592: // Licensed under the Apache License, Version 2.0 (the "License");
flatmax@592: // you may not use this file except in compliance with the License.
flatmax@592: // You may obtain a copy of the License at
flatmax@592: //
flatmax@592: //     http://www.apache.org/licenses/LICENSE-2.0
flatmax@592: //
flatmax@592: // Unless required by applicable law or agreed to in writing, software
flatmax@592: // distributed under the License is distributed on an "AS IS" BASIS,
flatmax@592: // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
flatmax@592: // See the License for the specific language governing permissions and
flatmax@592: // limitations under the License.
flatmax@592: /**
flatmax@592:     \author {Matt Flax <flatmax\@>}
flatmax@592:     \date 2013.02.08
flatmax@592: */
flatmax@592: 
flatmax@592: #include "Ear.H"
flatmax@592: 
flatmax@593: Ear::Ear(FP_TYPE fs_) {
flatmax@593:     fs=fs_; // set the specified sample rate
flatmax@593:     design();
flatmax@592: }
flatmax@592: 
flatmax@593: Ear::Ear(void) {
flatmax@593:     fs=DEFAULT_SAMPLERATE; // Use the default sample rate
flatmax@593:     design();
flatmax@592: }
flatmax@593: 
flatmax@593: Ear::~Ear(void) {
flatmax@593: }
flatmax@593: 
flatmax@593: void Ear::design(void) {
flatmax@593: 
flatmax@593:     // first figure out how many filter stages (PZFC/car.AC channels):
flatmax@593:     FP_TYPE pole_Hz = car.param.first_pole_theta * fs / (2.*M_PI);
flatmax@593:     n_ch = 0;
flatmax@593:     while (pole_Hz > car.param.min_pole_Hz) {
flatmax@593:         n_ch = n_ch + 1;
flatmax@593:         pole_Hz = pole_Hz - car.param.ERB_per_step *
flatmax@593:                   PsychoAcoustics::Hz2ERB(pole_Hz, car.param.ERB_break_freq, car.param.ERB_Q);
flatmax@593:     }
flatmax@593:     // Now we have n_ch, the number of channels, so can make the array
flatmax@593:     // and compute all the frequencies again to put into it:
flatmax@593:     car.pole_freqs.resize(n_ch, NoChange);
flatmax@593:     pole_Hz = car.param.first_pole_theta * fs / (2.*M_PI);
flatmax@593:     for (int ch = 0; ch<n_ch; ch++) {
flatmax@593:         car.pole_freqs[ch] = pole_Hz;
flatmax@593:         pole_Hz = pole_Hz - car.param.ERB_per_step *
flatmax@593:                   PsychoAcoustics::Hz2ERB(pole_Hz, car.param.ERB_break_freq, car.param.ERB_Q);
flatmax@593:     }
flatmax@593:     // now we have n_ch, the number of channels, and pole_freqs array
flatmax@593: 
flatmax@593:     max_channels_per_octave = (FP_TYPE)(log(2.) / log(car.pole_freqs[0]/car.pole_freqs[1]));
flatmax@593: 
flatmax@593:     // convert to include an ear_array, each w coeffs and state...
flatmax@601:     car.designFilters(fs, n_ch);
flatmax@601:     AGC.designAGC(fs, n_ch);
flatmax@593:     //IHC.designIHC(CF_IHC_params, fs, n_ch);
flatmax@593: }