benda/nonlinearbaseline2025
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

final touches on figure 3 and captions

Browse Source
This commit is contained in:
Jan Benda
2026-01-30 10:26:29 +01:00
parent 533a4245ca
commit c825c783dc
3 changed files with 8 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
twobeats.py
View File

@@ -223,9 +223,9 @@ def plot_stimulus(ax, s, tmax, eodf, freqs, c=0.1):
ax.show_spines('')
ax.plot(1000*time, eod, **s.lsEOD)
if len(ams) > 0:
ax.plot(1000*time, am, **ams)
ax.plot(1000*time, am + 0.1, **ams)
ax.set_xlim(0, 1000*tmax)
ax.set_ylim(-1.02 - 2*c, 1.02 + 2*c)
ax.set_ylim(-1.02 - 2*c, 1.1 + 2*c)
ax.text(0.5, 1.2, label, ha='center', transform=ax.transAxes)
@@ -293,7 +293,7 @@ if __name__ == '__main__':
s = plot_style()
fig, axs = plt.subplots(5, 4, cmsize=(s.plot_width, 0.6*s.plot_width),
height_ratios=[1, 0, 2, 1.3, 3, 0.7, 5])
fig.subplots_adjust(leftm=7.5, rightm=4.5, topm=1.5, bottomm=4,
fig.subplots_adjust(leftm=7.5, rightm=4, topm=1.5, bottomm=4,
wspace=0.4, hspace=0.4)
fmax = 250
tmin = 0.106
@@ -326,9 +326,9 @@ if __name__ == '__main__':
mark_freq(axs[4, 3], freqs, powers[3], df1, '', s.psF0)
mark_freq(axs[4, 3], freqs, powers[3], df1, '', s.psF01, 130, 1.5)
mark_freq(axs[4, 3], freqs, powers[3], abs(df1) + abs(df2) - 2,
f'$\\Delta f_2 + \\Delta f_1={abs(df1) + abs(df2):.0f}$\\,Hz', s.psF012, 20, angle=40)
f'$\\Delta f_2 + \\Delta f_1$\n$={abs(df1) + abs(df2):.0f}$\\,Hz', s.psF012, 0, toffs=-1, angle=45)
mark_freq(axs[4, 3], freqs, powers[3], abs(df1) - abs(df2),
f'$\\Delta f_2 - \\Delta f_1={abs(df1) - abs(df2):.0f}$\\,Hz', s.psF01_2, 50, toffs=5, angle=40)
f'$\\Delta f_2 - \\Delta f_1$\n$={abs(df1) - abs(df2):.0f}$\\,Hz', s.psF01_2, 65, toffs=-2, angle=45)
fig.common_yticks(axs[3, :])
fig.common_yticks(axs[4, :])
axs[3, 0].xscalebar(1, 0, 20, 'ms', ha='right')