[Débutant] filtre passe bande

[metall1986]

Bonjour,

j'ai un souci, je vous explique.
J'ai un signal.mat (son nom est signal_RMN1.mat)
pour afficher ce signal j'ai fait:

Code :

1
2
3
4
sig = signal_RMN1
load signal_RMN1.mat
plot(sig)
size(sig)

jusque là tout va bien. Mais là où je bloque c'est que je dois filtrer ce signal en haute et basse fréquence. Donc faire un filtre passe bande.

On me dit que:

Le signal RMN est la somme de toutes les contributions de chaque noyau, chacun émettant un signal :

S(t)=rα . (e^(t/Trα)). cos (2.pi.fα.t)

rα la concentration de noyaux de type α dans l’échantillon
Trα son temps de relaxation 
fα sa fréquence de résonance.
Le signal RMN est dans la gamme radiofréquence, mais il est démodulé pour abaisser sa fréquence avant numérisation. 
Ainsi fα est basse fréquence et peut même être négative selon le type de noyau et la fréquence centrale de démodulation.

Je sais que ma fréquence d'échantillonnage est fe = 10kHz.
Le signal est une mesure sur un échantillon qui contient un seul type de noyaux résonnant à 1000Hz après démodulation. Ce signal est fortement bruité en Haute et basse fréquence.

Je dois filtrer ce signal et renvoyer le signal filtrer. mais le problème c'est qu'en cours j'ai déjà fait ce genre d'exercice, mais j'avais la fréquence de coupure, j'avais rα....
Là je n'ai pas plus de donnée. Si quelqu'un sais comment faire ça serait sympa de m'aider.

Merci par avance pour ceux qui prendront de leur temps pour m'aider.
Cordialement