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

updating figure

Browse Source
This commit is contained in:
saschuta 2024-06-14 11:20:17 +02:00
parent 83ab41c4ac
commit 3d65e87d67
11 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
data/nonlin_regime_afe_0.01_nfft_400000_len_500_psd.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.01_nfft_400000_len_500_psdf.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.04_nfft_400000_len_500_psd.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.04_nfft_400000_len_500_psdf.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.1_nfft_400000_len_500_psd.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.1_nfft_400000_len_500_psdf.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.2_nfft_400000_len_500_psd.npy Normal file
View File

Binary file not shown.

BIN
data/nonlin_regime_afe_0.2_nfft_400000_len_500_psdf.npy Normal file
View File

Binary file not shown.

BIN
nonlin_regime.pdf
View File

Binary file not shown.

BIN
nonlin_regime.png
View File

Binary file not shown.
Side by Side Swipe Overlay

Before

Width:  |  Height:  |  Size: 57 KiB

After

Width:  |  Height:  |  Size: 58 KiB

5
nonlin_regime.py
View File

@ -231,6 +231,7 @@ def nonlin_regime(yposs=[450, 450, 450], printing=False, beat='',
ylim = ax_u1.get_ylim() ylim = ax_u1.get_ylim()
ax_u1.set_ylim(0, ylim[-1]) ax_u1.set_ylim(0, ylim[-1])
ax_u1.set_xlabel('Contrast [$\%$]') ax_u1.set_xlabel('Contrast [$\%$]')
ax_u1.set_ylabel('Amplitude [Hz]')
start = 200 # 1000 start = 200 # 1000
mults_period = 3 mults_period = 3
@ -259,7 +260,7 @@ def nonlin_regime(yposs=[450, 450, 450], printing=False, beat='',
trials_nr = 1 trials_nr = 1
nfft_here = values_nfft_full_model() nfft_here = values_nfft_full_model()
stimulus_length_here = 500 # values_stimuluslength_model_full() stimulus_length_here = 500#500 # values_stimuluslength_model_full()
nfft_here = 20 * 20000 nfft_here = 20 * 20000
a_f2s = [c_nrs_orig[c_nn]] a_f2s = [c_nrs_orig[c_nn]]
a_f2s = [c_nrs_orig[c_nn]] a_f2s = [c_nrs_orig[c_nn]]
@ -378,7 +379,7 @@ def nonlin_regime(yposs=[450, 450, 450], printing=False, beat='',
#if a == 2: # if (a_f1s[0] != 0) & (a_f2s[0] != 0): #if a == 2: # if (a_f1s[0] != 0) & (a_f2s[0] != 0):
#fish = '3 fish: $' + label_f_eod_name_core_rm() + '$\,\&\,' + f_vary_name() + '\,\&\,' + f_stable_name() # + '$'#' $\Delta '$\Delta$ #fish = '3 fish: $' + label_f_eod_name_core_rm() + '$\,\&\,' + f_vary_name() + '\,\&\,' + f_stable_name() # + '$'#' $\Delta '$\Delta$
beat_here = 'Contrast$\,=%s$' % (int(np.round(c_nrs[c_nn]))) + '$\%$' # +'$' beat_here = r'$\rm{Contrast}=%s$' % (int(np.round(c_nrs[c_nn]))) + '$\%$' # +'$'
#'$c_{1}=%s$' % (int(np.round(c_nrs[c_nn]))) + '$\%$' + ',\,$c_{2}=%s$' % (int(np.round(a_f2_cm[0]))) + '$\%$' # +'$' #'$c_{1}=%s$' % (int(np.round(c_nrs[c_nn]))) + '$\%$' + ',\,$c_{2}=%s$' % (int(np.round(a_f2_cm[0]))) + '$\%$' # +'$'
#plt.suptitle(f_vary_name(freq=int(freq1), delta=True)+', '+f_stable_name(freq=int(freq2), delta=True)) #plt.suptitle(f_vary_name(freq=int(freq1), delta=True)+', '+f_stable_name(freq=int(freq2), delta=True))