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