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()
Beispiel #2
0
def main():

    ths = calc_thresholds_rate(
        model=cochlea.run_zilany2014,
        model_pars={'species': 'cat'}
    )
    print(ths)

    ths.plot(logx=True)
    plt.show()
Beispiel #3
0
def main():

    ths = calc_thresholds_rate(
        model=cochlea.run_zilany2014,
        model_pars={'species': 'cat'}
    )
    print(ths)

    ths.plot(logx=True)
    plt.show()
Beispiel #4
0
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()