Smooth de courbe

Bonjour

Je doit tracer une courbe avec une quantité astronomique de donnée!! Ayant une oscillation trés importante de cette courbe, je souhaite utiliser la fonction smooth.

J'ai utilisé les méthodes 'moving', 'rlowess' et 'rloess' qui me donne des courbes à chaque foi différente. Seule les méthodes 'rlowess' et 'rloess' donnent des courbes interessante, cependant mon PC met un temps fou à les tracer. Pouvez vous me dire comment faire pour augmenter la vitesse de tracé? (si c'est possible).

Merci.

voilà le code que j'utilise

Code:

---

muu=smooth(Xpos(:,4),150,'moving');

---

a savoir que la matrice Xpos possèdent environ 150000 lignes...

Citation:

---

Envoyé par mygwel

Je doit tracer une courbe avec une quantité astronomique de donnée!! Ayant une oscillation trés importante de cette courbe, je souhaite utiliser la fonction smooth.

---

Tu pourrais nous dire à quoi correspondent physiquement ces données et aussi nous montrer un aperçu du tracé ?

Citation:

---

Envoyé par mygwel

mon PC met un temps fou à les tracer

---

Est-ce le tracé qui prend du temps ou bien le calcul de l'approximation par les fonctions que tu cites ?

je pense que c'est le calcul du tracé qui prend du temps. Lorsque j'utilise la méthode 'moving' j'en ai pour qq secondes, alors que si je prend la méthode 'rloess' ou 'rlowess', ça prend 15min... j'ai essayer d'augmenter le span, mais ça prend toujours autant de temps...

Physiquement, c'est le tracé du coefficient de frottement entre deux matériaux en fonction du nombre de cycles...

Je ne pense pas que se soit possible sans diminuer le nombre d'échantillon du signal.

Tu ne peux pas prendre 1 cycle sur 2 ou 1 cycle sur 5 ?

au départ le signal est saisi à une fréquence de 100Hz, que j'ai déjà réduit à 20Hz. Si je prend encore un cycle sur deux, j'ai peur de perdre beaucoup d'informations.

Citation:

---

Envoyé par mygwel

j'ai peur de perdre beaucoup d'informations.

---

As-tu au moins essayé ?

Et je fais comment pour extraire une oscillation sur 2? J'ai un fichier texte qui fait 150000 lignes, et le nombre de point de chaque oscillation n'est pas identique. Je me suis déjà arranger pour ne faire les calculs que sur les parties positives et négatives (je calcul les 2 séparément)

Tu peux faire quelque chose comme ceci :

Code:

---

1 2	pas = 5; plot(Xneg(1:pas:end,3),-Xneg(1:pas:end,4))

---

Enfin je suppose que tes données sont triées... sinon il faudra utiliser SORT pour les trier d'abord.

Code:

---

1 2	[Xneg(:,4),idx] = sort(Xneg(:,4)); Xneg(:,3) = Xneg(idx,3);

---

J'arrive un peu tard sur ce fil que je viens de découvrir.

Pour savoir si tu peux sous-échantillonner ton signal, effectue une transformée de Fourier (FFT). En fréquentiel, si le premier 5ième (resp. xième) des données seulement est significatif, tu peux sous-échantillonner par 5 (resp. x).

RemiZ

Salut!
Je pense que Dut t'a posé la bonne question:

Citation:

---

Tu pourrais nous dire à quoi correspondent physiquement ces données et aussi nous montrer un aperçu du tracé ?

---

Ce serait bien que tu y répondes.
Jean-Marc Blanc