def main(): cfs = np.logspace(np.log10(125), np.log10(16e3), 32) ths = calc_thresholds_rate( model=cochlea.run_zilany2014, cfs=cfs, model_pars={'species': 'human'}, asr_filter=True ) # ths = calc_thresholds_rate( # model=cochlea.external.run_matlab_auditory_periphery, # cfs=cfs, # model_pars={}, # asr_filter=True # ) human_ths = calc_human_hearing_thresholds(cfs) fig, ax = plt.subplots() human_ths.plot(ax=ax, logx=True, style='--', linewidth=5) ths.plot(ax=ax, logx=True) plt.show()
def main(): ths = calc_thresholds_rate( model=cochlea.run_zilany2014, model_pars={'species': 'cat'} ) print(ths) ths.plot(logx=True) plt.show()
def main(): cfs = np.logspace(np.log10(125), np.log10(16e3), 32) ths = calc_thresholds_rate(model=cochlea.run_zilany2014, cfs=cfs, model_pars={'species': 'human'}, asr_filter=True) # ths = calc_thresholds_rate( # model=cochlea.external.run_matlab_auditory_periphery, # cfs=cfs, # model_pars={}, # asr_filter=True # ) human_ths = calc_human_hearing_thresholds(cfs) fig, ax = plt.subplots() human_ths.plot(ax=ax, logx=True, style='--', linewidth=5) ths.plot(ax=ax, logx=True) plt.show()