Discussion :

[membreComplexe12]

Bonjour tous,

j'ai utilisé il y a un petit moment une machine qui doit me donner un signal qui ressemble à un signal triangulaire. A présent, je veux visualiser les résultats avec matlab qui sont dans le fichier texte ci joint (fichier 'courbe.txt' que j'ai compressé).

J'ai donc chargé mon fichier texte avec matlab et j'obtiens la courbe que j'ai mis en pièce jointe qui parait bien ressembler à un signal périodique quasi triangulaire (je changement d'amplitude à mi chemin est normal).

mais en fait, si on zoom (image n°2 en PJ) alors on se rend compte que ma machine à bien sauvegardée la courbe que je souhaite mais elle a sauvegardé aussi une deuxième courbe

Si on regarde la pièce jointe n°2 on voit en bleu les morceaux de courbes qui ne sont pas normaux :

Code :

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

1
2
3
4
5
 
figure (1); hold on;
data1=load('courbe.txt');
plot(data1(:,1),data1(:,2),'-','linewidth',2)
plot(data1(:,1),data1(:,2),'.r')

voyez vous comment faire pour virer dans mon fichier texte les points qui apparaissent en bleu sur ma courbe ?

moi je ne vois pas du tout ....

[Invité]

Bonsoir,

Les données semblent avoir les même abscisses, avec la fonction unique donc :

Code :

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

1
2
[uData1 ia ic] = unique(data1(:,1));
plot(uData1, data1(ia,2))

[membreComplexe12]

Bonsoir,

Bonsoir Jerome et merci beaucoup pour ton aide !

Les données semblent avoir les même abscisses, avec la fonction unique donc :

Code :

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

1
2
[uData1 ia ic] = unique(data1(:,1));
plot(uData1, data1(ia,2))

ça marche parfaitement
par contre, ça m'interesse de savoir comment tu as remarqué ceci car même en bien le fichier je n'ai pas l'impression qu'il y a exactement les mêmes abscisses

Code :

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

1
2
[uData1 ia ic] = unique(data1(:,1));
plot(uData1, data1(ia,2))

je n'ai pas totalement compris la fonction "unique" tel que tu l'utilise donc je me suis fais un petit exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
A=[1 2;
   2 2;
   3 2;
   4 2;
   4 8;
   5 2]
 
[uData1 ia ic] = unique(A(:,1));
Xnew=uData1
Ynew=A(ia,2)

je comprends que uData1 est le vecteur A(:,1) auquel on a enlever les premiers indices (pourquoi les premiers d'ailleurs et pas le second ? il y a moyen de changer ceci ?)

et d'apres ce que j'ai compris de l'aide, si je fais :

Code :

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

Xdepart=Xnew(ic)

alors Xdepart est bien mon vecteur de depart.

mais ce qui me semble bizarre c'est que ic à une taille plus grande que Xnew donc pourquoi je n'ai pas d'erreur qui dit que la taille de mon vecteur est trop petite ?

[Invité]

par contre, ça m'interesse de savoir comment tu as remarqué ceci

Visuellement, en zoomant et utilisant le datacursor.

je comprends que uData1 est le vecteur A(:,1) auquel on a enlever les premiers indices (pourquoi les premiers d'ailleurs et pas le second ? il y a moyen de changer ceci ?)

Paramètre 'first' / 'last' (défaut)

mais ce qui me semble bizarre c'est que ic à une taille plus grande que Xnew donc pourquoi je n'ai pas d'erreur qui dit que la taille de mon vecteur est trop petite ?

Petit passage de la doc

[C,ia,ic] = unique(A) also returns index vectors ia and ic, such that C = A(ia) and A = C(ic).

C'est parce qu'il y a des indices en double, ceux que font des valeurs non-uniques
Petit exemple :

Code :

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

1
2
3
unik = [1 2 3];
ic = randi(3,1,10);
unik(ic)

[membreComplexe12]

merci beaucoup Jerome pour tout
je crois avoir bien saisi à présent