Bonjour,

J'aimerais réaliser un filtre coupe bande (il me semble que c'est la traduction pour un notch filter) sur une serie temporelle (Y) contenant une echeance mensuelle sur 26 ans (soit 312 points).
J'ai désaisonnalisé ma série (en y enlevant la climatologie sur 12 mois) et j'aimerais enlever definitivement de ce signal la bande centrée sur 12 mois soit la bande [11mois 13mois].
Comme conseiller dans la faq, je comptais utiliser la fonction fir1
Code : Sélectionner tout - Visualiser dans une fenêtre à part
B = FIR1(N,Wn,'stop')
avec :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
 
fs=1;
Wn=2*[7e-2 9e-2]/fs;
Par contre, je n'est aucune idée de l'ordre (N) que je dois prendre et apres comment appliquer le filtre à ma série pour réaliser la fft et tracer la psd en fonction de la periode.

Si vous aviez une idée

Merci