jgrewe/gp_neurobio
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

buh

Browse Source
This commit is contained in:
efish 2018-11-29 16:45:59 +01:00
parent 28242c59d4
commit 3db401a3ea
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
code/eod_chirp_beat.py
View File

@ -47,7 +47,7 @@ def plot_chirp(eodf, eodf1, phase, axis):
y = chirp_eod * 0.4 + eod y = chirp_eod * 0.4 + eod
p, t = pd.detect_peaks(y, 0.1) p, t = pd.detect_peaks(y, 0.1)
axis.plot(time*1000, y) axis.plot(time*1000, y, color = 'royalblue')
axis.plot(time[p]*1000, (y)[p], lw=2, color='k') axis.plot(time[p]*1000, (y)[p], lw=2, color='k')
axis.plot(time[t]*1000, (y)[t], lw=2, color='k') axis.plot(time[t]*1000, (y)[t], lw=2, color='k')
axis.spines["top"].set_visible(False) axis.spines["top"].set_visible(False)
@ -63,18 +63,18 @@ ax2 = fig.add_subplot(222)
ax3 = fig.add_subplot(223) ax3 = fig.add_subplot(223)
ax4 = fig.add_subplot(224) ax4 = fig.add_subplot(224)
plot_chirp(600, 650, 0, ax1) plot_chirp(600, 650, 0, ax2)
plot_chirp(600, 650, np.pi, ax3) plot_chirp(600, 650, np.pi, ax4)
plot_chirp(600, 620, 0, ax2) plot_chirp(600, 620, 0, ax1)
plot_chirp(600, 620, np.pi, ax4) plot_chirp(600, 620, np.pi, ax3)
ax1.set_ylabel('EOD [mV]', fontsize=22) ax1.set_ylabel('EOD [mV]', fontsize=22)
ax1.set_title('$\Delta$f = 50 Hz', fontsize = 18) ax1.set_title('$\Delta$f = 20 Hz', fontsize = 18)
ax1.yaxis.set_tick_params(labelsize=18) ax1.yaxis.set_tick_params(labelsize=18)
ax1.set_xticklabels([]) ax1.set_xticklabels([])
ax2.set_title('$\Delta$f = 20 Hz', fontsize = 18) ax2.set_title('$\Delta$f = 50 Hz', fontsize = 18)
ax2.set_xticklabels([]) ax2.set_xticklabels([])
ax2.set_yticklabels([]) ax2.set_yticklabels([])
ax3.set_ylabel('EOD [mV]', fontsize=22) ax3.set_ylabel('EOD [mV]', fontsize=22)