Discussion :

[vanhelsing69]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

jusqu'a là tout va bien.
mais là ou 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)

Avec rα la concentration de noyaux de type α dans l’échantillon, Trα son temps de relaxation et 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 probleme c'est quand cours j'ai déja fait ce genre d'exercice, mais j'avais la fréquence de coupure, j'avais rα....
la je n'ai pas plus de donnée. si quelqu'un sais comment faire ça serai sympa de m'aider.

merci par avance pour ceux qui prendront de leur temps pour m'aider.
cordialement

[phryte]

Bonjour.

Ce signal est fortement bruité en Haute et basse fréquence.

Si tu connais les plages de fréquences qui sont bruitées tu peux utiliser le filtre de Butterworth :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
[numz,denz]=butter(n,[f1,f2]/(fe/2)) ;
 s=filter(numz,denz,signal)

[fregate]

si le but est juste de garder une bande de fréquence, a toi de choisir au mieux la fréquence de coupure.

puisque tu as la fréquence de résonance, tu suffit de couper un peu avant et un peu après la fréquence de résonance.

si tu es à 1000 hz, tu peux couper par exemple a 794 hz et 1258 hz (10^2,9 et 10 ^3,1) - pour avoir des frequence de coupure symétrique en log

[vanhelsing69]

merci pour vos réponses,

voila j'ai réussi a continuer:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
fe=10000; %fréquence échatillonnage fe=10kHz
N=10000;
f=[-fe/2:fe/N:fe*(N*1)/(2*N)]
y=fft(sig)
y1=fftshift(y)
plot(f,abs(y1))

j'ai fait la transformée de fourier afin d'avoir mon signal en fréquentiel, mais j'aimerai savoir comment faire pour que mon signal soit a zero de 0 à 900 hz et pareil de 1100 à 10 000hz. je sais que c'est brutal mais déjà si j'arrive a faire comme ça après j'ajusterai.
car une fois cela fait il me manquera plus qu'a refaire une transformée inverse. pour avoir mon signal filtré.

merci par avance pour ceux qui prendront de leur temps pour m'aider.
cordialement

[phryte]

Bonjour.

mais j'aimerai savoir comment faire pour que mon signal soit a zero de 0 à 900 hz et pareil de 1100 à 10 000hz. je sais que c'est brutal mais déjà si j'arrive a faire comme ça après j'ajusterai.

Tu as un exemple là :
http://www.developpez.net/forums/d54...es-signal-wav/

[vanhelsing69]

oui je veux bien, faire comme ça, mais j'ai 10 000 valeurs et ce sont des vecteurs donc une valeur réelle et une valeur imaginaire.
donc je vais pas les taper un par un.
c'est justement ça le probleme

[phryte]

Bonjour.

donc je vais pas les taper un par un.

Il suffit de faire le gabarit de 0 à 900 Hz et de 1100 à 10 000 Hz !

[vanhelsing69]

Le problème c'est que je n'arrive pas a faire le gabarit. je suis tout nouveau sur matlab. Et j'ai fait deux exercices dessus mais j'étais guidé. en gros c'était une recette a suivre, j'avais les valeurs et les fonctions a utiliser. et la j'ai un signal, et je doit le filtrer et le renvoyer.

http://dl.free.fr/iQUElGcb3

voici le lien du signal que je dois filtrer. c'est un fichier zipé, dedans il y a deux autres fichiés zipé, un pour matlab version7 et un pour matlab version inférieur a la version 7. ( le fichier a traité s'appelle signal_RMN1.mat)

je vous remet le code que j'ai fait pour le moment:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
sig = signal_RMN1
load signal_RMN1.mat
plot(sig)
size(sig)
fe=10000; %fréquence échatillonnage fe=10kHz
N=10000;
f=[-fe/2:fe/N:fe*(N*1)/(2*N)]
y=fft(sig)
y1=fftshift(y)
plot(f,abs(y1))

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.
si vous pouvez m'aider, car j'ai vraiment besoin d'aide

[vanhelsing69]

merci pour votre aide.
Le code pour filtrer ce signal sera mis la semaine prochaine pour ce que ça intéresse.
cordialement