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

upd

Browse source
This commit is contained in:
François Boulogne 2025-09-01 13:45:37 +02:00
parent 241abc1657
commit 0c423acead
2 changed files with 13 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

2
mysignal/__init__.py
View file

@ -1 +1 @@
__version__ = '0.1.9'
__version__ = '0.1.10'

19
mysignal/phasefreq.py
View file

@ -574,19 +574,24 @@ def analyze_signal_wavelet(time, e_signal, s_signal, freq_rtol=0.01):
def plot_phases(e_time, e_filtered, s_filtered, n_modes, callback=analyze_signal_wavelet):
fig, ax = plt.subplots(nrows=n_modes, figsize=(12, 6))
fig, ax = plt.subplots(nrows=n_modes, ncols=2, figsize=(12, 6), gridspec_kw={'width_ratios': [3, 1]})
for mod in range(n_modes):
res = callback(e_time, e_filtered[mod], s_filtered[mod], freq_rtol=0.3)
ax[mod].set_title(f'Freq: {res['freq_e']:.3f}, Phase: {res['phase']:.3f}, Delay: {res['delay']:.3f}' )
ax[mod].plot(e_time, e_filtered[mod], label='e')
ax[mod].plot(e_time, s_filtered[mod], label='s')
ax[mod].set_xlim(left=20, right=230)
for a in ax:
ax[mod, 0].set_title(f'Freq: {res['freq_e']:.3f}, Phase: {res['phase']:.3f}, Delay: {res['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)
ax[mod, 1].plot(e_time, e_filtered[mod], label='e')
ax[mod, 1].plot(e_time, s_filtered[mod], label='s')
ax[i+1, 1].set_xlim(left=0.5 * e_time.mean(), right=0.5*e_time.mean() + 10/e_frequencies[i])
for a in ax[:, 0]:
a.legend()
plt.tight_layout();