research/mysignal
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

update

Browse source
This commit is contained in:
François Boulogne 2026-01-06 13:46:32 +01:00
parent 489c255abd
commit 882f35cb5c
1 changed files with 6 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

8
mysignal/phasefreq.py
View file

@ -594,7 +594,7 @@ def get_phases(e_time, e_filtered, e_frequencies, s_filtered, n_modes, callback=
return pd.DataFrame(res)
def plot_phases(e_time, e_filtered, e_frequencies, s_filtered, n_modes, callback=analyze_signal_wavelet, output=None):
def plot_phases(e_time, e_filtered, e_frequencies, s_filtered, n_modes, df, output=None):
"""
"""
@ -603,7 +603,11 @@ def plot_phases(e_time, e_filtered, e_frequencies, s_filtered, n_modes, callback
for mod in range(n_modes):
res = callback(e_time, e_filtered[mod], s_filtered[mod], freq_rtol=0.3)
ax[mod, 0].set_title(f'Freq: {res['freq_e']:.3f}, Phase: {res['phase']:.3f}, Delay: {res['delay']:.3f}' )
freq = df[df['mode'] == mod].get('freq')
phase = df[df['mode'] == mod].get('phase')
delay = df[df['mode'] == mod].get('delay')
ax[mod, 0].set_title(f'Freq: {freq:.3f}, Phase: {phase:.3f}, Delay: {delay:.3f}' )
ax[mod, 0].plot(e_time, e_filtered[mod], label='e')
ax[mod, 0].plot(e_time, s_filtered[mod], label='s')
#ax[mod].set_xlim(left=20, right=230)