Discussion :

[membreComplexe12]

Salut tous,

j'aurais besoin une fois de plus de votre aide si ça ne vous dérange pas

j'ai une courbe d'hystéresis (cf. PJ) qui n'est pas aussi belle que celle que j'ai mis en image. En fait lors de la mesure il y a de petites perturbation (très légères) qui me gêne pour le traitement de la courbe.

Ce que je voudrais faire c'est filtrer ma courbe afin de ne pas avoir de petites oscillations parasites et donc avoir un signal parfaitement lisse tout comme l'image en PJ.

le soucis est que je n'y connait absolument rien en filtrage de données :
- quel type de filtre ?
- pourquoi ?
- principe de fonctionnement ?
- déjà implémenté sous matlab, lequel choisir ?

j'espère que vous pourrez m'éclairer tous ceci...

merci d'avance

EDIT :

j'ai trouvé sur le forum ce bout de code qui est très bien :

Code :

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

1
2
3
4
5
6
7
m=50;
b=ones(1,m);
a=m;
x=2*sin(2*pi*.013*(1:500))+randn(1,500);
y=filter(b,a,x);
subplot(2,1,1);plot(x);
subplot(2,1,2);plot(y);

mais comme j'ai compris ce code filtre les donnés d'un vecteur X alors que moi je veux filtrer un couple (X,Y) de données (c'est pareil ?)

Le soucis :

1°) le premier soucis est que je ne comprends pas vraiment le principe (par exemple pourquoi il y a une matrice unité qui intervient ?)
2°) ensuite le deuxieme soucis :
-> là la première courbe est très bruité et la deuxieme ne l'ai presque plus
-> or dans mon cas j'ai une courbe qui est peu bruité (mais c'est quand même dérangeant), ma courbe ressemble en fait à la deuxieme courbe de cet exemple et je voudrais savoir si il est possible de filtrer encore mieux que cela ?

merci d'avance
A+

[FR119492]

Salut!
En fait, les cycles d'hystérèse d'un matériau ferromagnétique sont très mal définis: si tu prélèves des échantillons sur des livraisons successives de tôles de "même" type chez le même fournisseur, tu as de fortes chances d'avoir des cycles d'hystérèse assez différents. Considérer un certain cycle lissé comme une "vérité" est donc une illusion. Je te conseille donc une approche beaucoup plus simple: tu imagines une fonction qui ressemble le plus possible à celle représentée sur ta pièce jointe et tu la programmes.
Jean-Marc Blanc

[membreComplexe12]

merci de cette réponse.

L'idée est pas mal mais je ne sais pas quelle courbe prendre pour fiter ce type de chose (d'autant plus qu'ici ce n'est qu'un cas simplifié et en réalité la hauteur des boucles peu évoluer et il y a plusieurs cycles)

en fait je pense que le plus simple est bien de filtrer la courbe de ce bruit de mesure

ps : (oscillations dues à mon banc de mesure (un peu spécial ) qui vibre légèrement. Si je le bloque pendant un petit instant le bruit disparait)

EDIT :

en fouillant sur le net j'ai trouvé ceci pour filtre passe bas :

Code :

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

1
2
 
Ynew=filter(a,[1 a-1],Y)

j'ai testé mais je ne sais pas comment régler "a" de façon optimale, c'est au "taton" que ça se fait ?

=> le filtrage induit un décalage de ma courbe c'est normal ?

[le fab]

en fouillant sur le net j'ai trouvé ceci pour filtre passe bas :

Code :

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

1
2
 
Ynew=filter(a,[1 a-1],Y)

j'ai testé mais je ne sais pas comment régler "a" de façon optimale, c'est au "taton" que ça se fait ?

=> le filtrage induit un décalage de ma courbe c'est normal ?

pour régler a, teste plusieurs valeurs entre 0 et 1 et choisi " la plus belle" courbe filtrée
à a=1 tu récupère le signal de départ (pas de filtrage)
à a=0 tu récupère un vecteur nul (trop de filtrage)

le décalage est normal, pour l'enlever, il faut filtrer dans l'autre sens :

Code :

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

1
2
3
4
5
 
Ynew=filter(a,[1 a-1],Y);
ynew = ynew(end:-1:1);
Ynew=filter(a,[1 a-1],Ynew);
ynew = ynew(end:-1:1);

par contre prendre un a plus proche de 1 quand tu filtres 2 fois que lorsque tu filtres une fois

[membreComplexe12]

pour régler a, teste plusieurs valeurs entre 0 et 1 et choisi " la plus belle" courbe filtrée
à a=1 tu récupère le signal de départ (pas de filtrage)
à a=0 tu récupère un vecteur nul (trop de filtrage)

le décalage est normal, pour l'enlever, il faut filtrer dans l'autre sens :

Code :

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

1
2
3
4
5
 
Ynew=filter(a,[1 a-1],Y);
ynew = ynew(end:-1:1);
Ynew=filter(a,[1 a-1],Ynew);
ynew = ynew(end:-1:1);

par contre prendre un a plus proche de 1 quand tu filtres 2 fois que lorsque tu filtres une fois

merci beaucoup pour ces explications