def main():
fs = 100e3
### Make sound
t = np.arange(0, 0.1, 1 / fs)
s = dsp.chirp(t, 80, t[-1], 20000)
s = cochlea.set_dbspl(s, 50)
pad = np.zeros(10e-3 * fs)
sound = np.concatenate((s, pad))
### Run model
anf = cochlea.run_zilany2014(
sound, fs, anf_num=(100, 0, 0), cf=(125, 20000, 100), seed=0, powerlaw="approximate", species="human"
)
### Accumulate spike trains
anf_acc = th.accumulate(anf, keep=["cf", "duration"])
anf_acc.sort("cf", ascending=False, inplace=True)
### Plot auditory nerve response
fig, ax = plt.subplots(2, 1)
th.plot_signal(signal=sound, fs=fs, ax=ax[0])
th.plot_neurogram(anf_acc, fs, ax=ax[1])
plt.show()
def main():
fs = 48e3
### Make sound
tmax = 0.03
t = np.arange(0, tmax, 1 / fs)
s = dsp.chirp(t, 80, t[-1], 16000)
sound = cochlea.set_dbspl(s, 50)
### Run model
anf = run_matlab_auditory_periphery(
sound,
fs,
anf_num=(100, 50, 20),
cf=(125, 16000, 80),
seed=0,
)
### Accumulate spike trains
anf_acc = th.accumulate(anf, keep=['cf', 'duration'])
anf_acc.sort('cf', ascending=False, inplace=True)
### Plot auditory nerve response
fig, ax = plt.subplots(2, 1, sharex=True)
th.plot_signal(signal=sound, fs=fs, ax=ax[0])
th.plot_neurogram(anf_acc, fs, ax=ax[1])
plt.show()
def main():
fs = 48e3
### Make sound
t = np.arange(0, 0.1, 1 / fs)
s = dsp.chirp(t, 80, t[-1], 20000)
s = cochlea.set_dbspl(s, 50)
s = np.concatenate((s, np.zeros(int(10e-3 * fs))))
### Run model
anf_trains = cochlea.run_holmberg2007(
s,
fs,
anf_num=(100, 0, 0),
seed=0,
)
### Plot auditory nerve response
anf_acc = th.accumulate(anf_trains, keep=['cf', 'duration'])
anf_acc.sort('cf', ascending=False, inplace=True)
fig, ax = plt.subplots(2, 1)
th.plot_signal(signal=s, fs=fs, ax=ax[0])
th.plot_neurogram(anf_acc, fs, ax=ax[1])
plt.show()
def main():
fs = 100e3
### Make sound
t = np.arange(0, 0.1, 1/fs)
s = dsp.chirp(t, 80, t[-1], 20000)
s = cochlea.set_dbspl(s, 50)
s = np.concatenate( (s, np.zeros(10e-3 * fs)) )
### Run model
anf = cochlea.run_zilany2009(
s,
fs,
anf_num=(100,0,0),
cf=(80, 20000, 100),
seed=0,
powerlaw='approximate'
)
### Plot auditory nerve response
anf_acc = th.accumulate(anf, keep=['cf', 'duration'])
anf_acc.sort('cf', ascending=False, inplace=True)
cfs = anf.cf.unique()
fig, ax = plt.subplots(2,1, sharex=True)
th.plot_neurogram(
anf_acc,
fs,
ax=ax[0]
)
th.plot_raster(
anf[anf.cf==cfs[30]],
ax=ax[1]
)
ax[1].set_title("CF = {}".format(cfs[30]))
plt.show()
def main():
fs = 48e3
### Make sound
tmax = 0.03
t = np.arange(0, tmax, 1/fs)
s = dsp.chirp(t, 80, t[-1], 16000)
sound = cochlea.set_dbspl(s, 50)
### Run model
anf = run_matlab_auditory_periphery(
sound,
fs,
anf_num=(100,50,20),
cf=(125, 16000, 80),
seed=0,
)
### Accumulate spike trains
anf_acc = th.accumulate(anf, keep=['cf', 'duration'])
anf_acc.sort('cf', ascending=False, inplace=True)
### Plot auditory nerve response
fig, ax = plt.subplots(2, 1, sharex=True)
th.plot_signal(
signal=sound,
fs=fs,
ax=ax[0]
)
th.plot_neurogram(
anf_acc,
fs,
ax=ax[1]
)
plt.show()
def main():
fs = 100e3
# Make sound
t = np.arange(0, 0.1, 1/fs)
s = dsp.chirp(t, 80, t[-1], 20000)
s = cochlea.set_dbspl(s, 50)
s = np.concatenate((s, np.zeros(int(10e-3*fs))))
# Run model
anf = cochlea.run_zilany2009(
s,
fs,
anf_num=(100,0,0),
cf=(80, 20000, 100),
seed=0,
powerlaw='approximate'
)
# Plot auditory nerve response
anf_acc = th.accumulate(anf, keep=['cf', 'duration'])
anf_acc.sort_values('cf', ascending=False, inplace=True)
cfs = anf.cf.unique()
fig, ax = plt.subplots(2,1, sharex=True)
th.plot_neurogram(
anf_acc,
fs,
ax=ax[0]
)
th.plot_raster(
anf[anf.cf==cfs[30]],
ax=ax[1]
)
ax[1].set_title("CF = {}".format(cfs[30]))
plt.show()
def main():
fs = 100e3
# Make sound
t = np.arange(0, 0.1, 1/fs)
s = dsp.chirp(t, 80, t[-1], 20000)
s = cochlea.set_dbspl(s, 50)
pad = np.zeros(int(10e-3 * fs))
sound = np.concatenate( (s, pad) )
# Run model
anf = cochlea.run_zilany2014(
sound,
fs,
anf_num=(100,0,0),
cf=(125, 20000, 100),
seed=0,
powerlaw='approximate',
species='human',
)
# Accumulate spike trains
anf_acc = th.accumulate(anf, keep=['cf', 'duration'])
anf_acc.sort_values('cf', ascending=False, inplace=True)
# Plot auditory nerve response
fig, ax = plt.subplots(2,1)
th.plot_signal(
signal=sound,
fs=fs,
ax=ax[0]
)
th.plot_neurogram(
anf_acc,
fs,
ax=ax[1]
)
plt.show()
def main():
fs = 48e3
### Make sound
t = np.arange(0, 0.1, 1/fs)
s = dsp.chirp(t, 80, t[-1], 20000)
s = cochlea.set_dbspl(s, 50)
s = np.concatenate( (s, np.zeros(10e-3 * fs)) )
### Run model
anf_trains = cochlea.run_holmberg2007(
s,
fs,
anf_num=(100,0,0),
seed=0,
)
### Plot auditory nerve response
anf_acc = th.accumulate(anf_trains, keep=['cf', 'duration'])
anf_acc.sort('cf', ascending=False, inplace=True)
fig, ax = plt.subplots(2,1)
th.plot_signal(
signal=s,
fs=fs,
ax=ax[0]
)
th.plot_neurogram(
anf_acc,
fs,
ax=ax[1]
)
plt.show()
