raab/GP2023_chirp_detection
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

checking the instantaneos_frequency

Browse Source
This commit is contained in:
wendtalexander 2023-01-16 20:07:38 +01:00
parent 133aff26a1
commit 990e858d02
1 changed files with 6 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
code/chirp_sim.py
View File

@ -26,19 +26,20 @@ plt.close()
# trying thunderfish fakefish chirp simulation --------------------------------- # trying thunderfish fakefish chirp simulation ---------------------------------
samplerate = 44100 samplerate = 44100
freq, ampl = fakefish.chirps(eodf=500, chirp_contrast=0.2) freq, ampl = fakefish.chirps(eodf=500, chirp_contrast=0.2)
data = fakefish.wavefish_eods(fish='Alepto', frequency=freq, fphase0=3) data = fakefish.wavefish_eods(fish='Alepto', frequency=freq, phase0=3)
# filter signal with bandpass_filter # filter signal with bandpass_filter
data_filterd = bandpass_filter(data*ampl, samplerate, 495, 505) data_filterd = bandpass_filter(data*ampl+1, samplerate, 0.01, 1.99)
data_freq_time, data_freq = instantaneos_frequency(data_filterd, samplerate) data_freq_time, data_freq = instantaneos_frequency(data*ampl, samplerate)
fig, ax = plt.subplots(4, 1, figsize=(20 / 2.54, 12 / 2.54), sharex=True) fig, ax = plt.subplots(4, 1, figsize=(20 / 2.54, 12 / 2.54), sharex=True)
ax[0].plot(np.arange(len(data))/samplerate, data*ampl) ax[0].plot(np.arange(len(data))/samplerate, data*ampl+1)
ax[1].plot(np.arange(len(data_filterd))/samplerate, data_filterd) ax[1].plot(np.arange(len(data_filterd))/samplerate, data_filterd)
ax[2].plot(np.arange(len(freq))/samplerate, freq) ax[2].plot(np.arange(len(freq))/samplerate, freq)
ax[3].plot(data_freq_time[1:], data_freq[1:]) ax[3].plot(data_freq_time[1:], data_freq[1:])
plt.show() plt.show()
embed()