Analyse spectrale

Chargement et affichage du signal :

x = np.loadtxt("../_static/consommation.csv",
               delimiter=",", skiprows=1, usecols=2)

Spectre du signal…

La phase est, comme la plupart des signaux réels, difficile à analyser… Mais le module est informatif.

En zoomant sur les fréquences basses du spectre (matplotlib.pyplot.xlim), on peut identifier les fréquences principales. La fonction scipy.signal.find_peaks peut aider à trouver les pics.

peaks, _ = find_peaks(module, height=1e8)

On peut alors afficher les fréquences des raies principales, et la période correspondante.

0.006 h^-1 :	phénomène qui se répète 0.006 fois par heure soit toutes les 168.143 h (ou 7.006 jours)
0.042 h^-1 :	phénomène qui se répète 0.042 fois par heure soit toutes les  24.000 h (ou 1.000 jours)
0.083 h^-1 :	phénomène qui se répète 0.083 fois par heure soit toutes les  12.000 h (ou 0.500 jours)

On retrouve là les périodicités hebdomadaire, quotidienne et même sur 12 heures de la consommation électrique.