Discussion :

[gmachi10]

Bonjour. Je ne sais pas si je trouverai une solution sur ce forum, mais je pense que nombreux "collègues" qui développe des applications sous matlab ont du être confrontés à ce problème.

Voilà je réalise plusieurs acquisition d'un signal (10,15,...). en fonction du temps.

Si on se place dan le cadre théorque, les signaux que je possède doivent être rigoureusement égaux. Cependant les signaux sont assujetti à un bruit, en général blanc (ne pose à priori pas de problème), mais peuvent être assujetti à des valeurs aberrantes qui perturbent un signal (vibration, courant mal stabilisé....).

J'aimerais savoir comment faire pour détecter à partir d'une acquisition de N signaux s'il y a 1 ou plusieurs signal aberrant auquel avec un intervalle de confiance, on puisse rejeter le signal.

J'ai déjà essayé quelques méthodes mais sans résultat fiable.

[phryte]

Salut.

1 ou plusieurs signal aberrant auquel avec un intervalle de confiance, on puisse rejeter le signal.

Je calculerai la moyenne et l'écart-type sur un horizon (pondéré ou pas) compatible avec la dynamique de tes signaux. Au-delà d'un seuil (ou d'un %) le signal est rejeté. (Si tu ne t'intéresse qu'à l'amplitude).

[gmachi10]

Voilà ce que j'ai déjà essayé.
boucle for sur nombre d'acquisitions
j'enlève le i-ème enregistrement, suivi du calcul sur chacune des lignes de l'écart type à la fin des lignes je calcul l'écart type moyen que je stocke. Et cela donc sur mon nombre total d'acqusition.
Ensuite je compare l'écart type moyen avec chacun des ecart type, mais comme mes signaux dépendent de mon échantillon, je ne peux définir de réel critère.

Avez vous d'autres soliutions?

[phryte]

Salut.
Quelle est la nature de ton signal (constant, évolutif) comment calcules-tu la moyenne...

[gmachi10]

Pour répondre à ta question, je vais donner une image simple.

Je réalise sur une mesure sur un échantillon sur une durée d'acquisition après excitation de 10s.
Sur cette période là, ma première courbe dirons nous est un sinus d'amplitude A et de fréquence f.

L'amplitude en fonction du temps n'est donc pas la même. cependant si je réalise une mesure 2s après, je doit retrouver cette même courbe.

Je souhaite donc que si au lien de trouver le signal Asin(2pif t) , je trouve des aberrations sur le signal saut d'amplitude, déphasage important, ça me les détecte.

La moyenne est une moyenne normale. somme:/nbacquisition

[phryte]

somme:/nbacquisition

Cela est la moyenne globale.
Il faut faire une moyene glissante et un écart-type glissant pour détecter les perturbations. Quelle est la période d'échantillonnage ?

[gmachi10]

Pour éviter les problèmes de vocabulaire, le mot échantillon correspond pour moi à un matériau à caractériser.

Mes mesures s'effectuent sur 1024 ,2048 ou 4096 pts.

[phryte]

Salut.

des aberrations sur le signal

Il faudrait caractériser le signal ou avoir une référence.
Matlab possède la fonction "coherence(s1,s2,nfft,fs)".

[gmachi10]

En fait le signal dépend de l'échantillon dont on ne connait pas forcément sa nature donc on ne peut pas le cractériser. La seule chose est que quelque soit le matériau, un mesure de sa réponse doit être reproductible. Ce que je cherche à éliminer c'est un signal de mesure qui serait aberrant par rapport aux autres. Dans aberrant, j'attends contient au moins une valeur qui s'écarte totalement des autres valeurs et qui erronne donc le signal;

[phryte]

Salut.
Peux-tu nous faire voir un ou deux signaux.

[gmachi10]

Voici un signal .
Normalement sur un même échantillon, tous les signaux sont identiques en tout point , au bruit blanc près.

UNe aberration sur le signal peut se traduire par un déphasage, problème au niveau de l'amplitude.....

[noliv]

As tu essayé de voir si tu pouvais t'en sortir en fitant les courbes obtenues sur des échantillons de même nature? une fois que tu as une loi "théorique", tu devrais pouvoir réjecter les valeurs abberantes. Forcément, ce n'est pas automatique... A moins que tu puisses fiter tous les échantillons avec des polynome d'ordres égaux.