Fonction Scipy.signal.Welch avec trous dans signal
Bonjour,
Je réalise une analyse PSD d'un signal. Cette analyse se fait sur 17 overlapping de 200 secondes donnant un temps d'analyse de 30 minutes.
Normalement, sans trous dans mon signal, j'ai 512 lignes toutes les 200 secondes.
J'utilise alors le code suivant :
Code:
1
2
3
4
5
6
7
| T=200
fs=2.56
ts=1/fs
N=int(T/ts)
df_temp=pd.concat([df.loc[(df.Time>=df.loc[0,'Time']+m*60*30) & (df.Time<(1+m)*60*30+df.loc[0,'Time'])],pd.DataFrame()],ignore_index=True) # m is just to select the block of 30 minutes we want
signal=df_temp.Elevation.values #Convert in numpy format
f , pxx=sp.signal.welch(signal,fs=fs,window='hann',nperseg=N) |
avec df la base de donnée initiale et df_temp la base de donnée contenant 30 minutes de signal.
Ma question est de savoir si la fonction Welch est toujours performante si j'ai un trou de 70 secondes par exemple (car du coup entre les deux lignes avant et après le trou, la fréquence d'acquisition n'est plus de 2.56 Hz et le nombre de ligne varie.
Merci de votre aide
