Discussion :

[j_ere_my]

Bonjour,
petite question ...

j'ai un vecteur A, j'utilise la fonction "find" pour définir un nouveau vecteur B suivant les conditions désirées sur d'autres vecteurs.
maintenant, je veux le vecteur C qui serait les autres valeurs de A

sachant, que mes conditions sont diverses et nombreuses, je ne veux pas faire un "find" pour mon C.
comment puis-je obtenir C, qui est donc le vecteur A sans les lignes de B ??


merci

[Loïc B.]

Vu ici :

http://www.developpez.net/forums/d68...ments-vecteur/

La fonction à utiliser serait setdiff.

Bonne journée

[j_ere_my]

Bon, je vais préciser ma question

je ne cherche pas à enlever les valeurs identiques mais les lignes

A ==> B ==> C
a1 a2 a1
a2 a4 a3
a3 an a5
a4
a5
an

sachant que les valeurs de a1 peuvent etre egales a a2, par exemple

[Caro-Line]

Peux-tu montrer exactement comment tu as trouvé B ?
Normalement avec find tu récupères les indices.
Donc les indices des lignes restantes sont facilement récupérables avec SETDIFF.

Exemple :

Code :

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

1
2
3
4
A=rand(1,10)
idx=find(A>0.5)
B=A(idx)
C=A(setdiff(1:length(A),idx))

Sachant que cela peut s'écrire plus "joliment" avec l'indexage logique :

Code :

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

1
2
3
4
A=rand(1,10)
idx=A>0.5
B=A(idx)
C=A(~idx)

[Loïc B.]

Bon, je vais préciser ma question

je ne cherche pas à enlever les valeurs identiques mais les lignes

A ==> B ==> C
a1 a2 a1
a2 a4 a3
a3 an a5
a4
a5
an

sachant que les valeurs de a1 peuvent etre egales a a2, par exemple

??? Tu dis ne pas vouloir enlever les valeurs mais les lignes. (a1 peut être égal à a2). Donc si je reprends ton exemple avec des chiffres tu veux :
A : 1 1 2 3 4 5
B : 1 2 3
C : 1 4 5

Je ne vois pas comment cela est possible, vu que avec les éléments que tu nous donne, on a aucun moyen de distinguer le premier "1" du second "1".

[Caro-Line]

Je ne vois pas comment cela est possible, vu que avec les éléments que tu nous donne, on a aucun moyen de distinguer le premier "1" du second "1".

Si par les indices.
La solution que j'ai proposée fonctionnera très bien.

[j_ere_my]

ok ca a marche, merci beaucoup

en fait, il me faut 3 vecteurs dans A

Code :

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

1
2
3
4
5
strat=find(cl1pv12>3 & cl2pv12>3 & cl3pv12>3 & cl4 .....
xst=x(strat);
 
tropo=find(cl1pv12<2 & cl2pv12<2 & cl3pv12<2 & cl4pv12<2 .....
xtr=x(tropo);

je veux donc toutes les lignes n'appartenant ni a strat ni a tropo, appelé suspe2

j'ai fais :

Code :

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

1
2
suspe=x(setdiff(1:length(x),strat));  !!!! j'enleve le vecteur strat
    suspe2=suspe(setdiff(1:length(suspe),tropo));    !!!!! j'enleve le vecteur tropo

mais j'ai pas lenght x = lenght tropo + lenght strat + lenght suspe2
j'ai un leger excedent dans le total des termes de droites

si qq'un peut m'aider ....



>> size (x)

ans =

28217 1

>> size (strat)

ans =

6364 1

>> size (suspe)

ans =

21853 1

>> size (tropo)

ans =

3536 1

>> size (suspe2)

ans =

18949 1

[Caro-Line]

Ah ben non ça ne peut pas marcher !

Dans la variable tropo tu as les indices par rapport à x pas par rapport à suspe.

Pourquoi ne donnes-tu pas sa chance à l'indexage logique ?

Code :

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

1
2
3
4
5
6
7
strat=cl1pv12>3 & cl2pv12>3 & cl3pv12>3 & cl4 .....
xst=x(strat);
 
tropo=cl1pv12<2 & cl2pv12<2 & cl3pv12<2 & cl4pv12<2 .....
xtr=x(tropo);
 
suspe2=x(~strat & ~tropo)

[j_ere_my]

??? Error using ==> and
Matrix dimensions must agree.

Error in ==> clustplot at 359
suspe=x(~strat & ~tropo);

[j_ere_my]

bon ca a marché

je te remercie caroline, j'ai laissé sa chance à l'indexage logique

Bonne soiree

[Caro-Line]

Ah oui le problème c'est que tes conditions ne sont pas sur le vecteur x...
Par curiosité c'est quoi les variables cl1pv12 ... ?
Est-il possible d'avoir les lignes entiètes des conditions ?

Parce que si tu utilises FIND il faut quand même bien que tu trouves des indices par rapport à x (sinon ce n'est pas logique, et tu risques de tomber sur des indices qui ne sont pas dans 1:length(x)).

[Caro-Line]

bon ca a marché

je te remercie caroline, j'ai laissé sa chance à l'indexage logique

Bonne soiree

Bon nos messages se sont croisés.

J'espère surtout que tu as compris comment cela fonctionne afin de pouvoir le réutiliser par la suite

[j_ere_my]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
 
cl1pv12=x1(:,13); .....
cl1pv12=cl1pv12(~isnan(x.*y))....
 
strat=(cl1pv12>3 & cl2pv12>3 & cl3pv12>3 & cl4pv12>3 & cl5pv12>3 & cl1pv24>3 & cl2pv24>3 & cl3pv24>3 & cl4pv24>3 & cl5pv24>3 & cl1pv36>3 & cl2pv36>3 & cl3pv36>3 & cl4pv36>3 & cl5pv36>3 & cl1pv48>3 & cl2pv48>3 & cl3pv48>3 & cl4pv48>3 & cl5pv48>3 & cl1pv60>3 & cl2pv60>3 & cl3pv60>3 & cl4pv60>3 & cl5pv60>3 & cl1pv72>3 & cl2pv72>3 & cl3pv72>3 & cl4pv72>3 & cl5pv72>3);
 
xst=x(strat);
yst=y(strat);
 
tropo=(cl1pv12<2 & cl2pv12<2 & cl3pv12<2 & cl4pv12<2 & cl5pv12<2 & cl1pv24<2 & cl2pv24<2 & cl3pv24<2 & cl4pv24<2 & cl5pv24<2 & cl1pv36<2 & cl2pv36<2 & cl3pv36<2 & cl4pv36<2 & cl5pv36<2 & cl1pv48<2 & cl2pv48<2 & cl3pv48<2 & cl4pv48<2 & cl5pv48<2 & cl1pv60<2 & cl2pv60<2 & cl3pv60<2 & cl4pv60<2 & cl5pv60<2 & cl1pv72<2 & cl2pv72<2 & cl3pv72<2 & cl4pv72<2 & cl5pv72<2);
% 
xtr=x(tropo);
ytr=y(tropo);
 
 
 
    suspe=x(~strat & ~tropo);
    xsus=x(suspe);
    ysus=y(suspe);
 
 
fprintf(fid,'%f %f....\n',x1(suspe,1),x1(suspe,2)....);    
fclose(fid);
 
figure (1)
    %hold on;
 
   scatter(xsus,ysus,3,'k')

j'ai bien la bonne taille pour suspe, xsus, ysus
mais qd je trace mon scater, g plus du tt le meme nombre de points, et mon ficher fid pareil

tu as raison, je dois m'embrouiller dans les indexations

pour répondre à ta curiosité, les conditions clnpvn correspondent à des informations lagrangiennes de rétrotrajectoires dans le temps, notamment ici, le tourbillon potentiel , traceur de masses d'aires.

[Caro-Line]

Donc je suppose que tous les vecteurs utilisés pour tes conditions ont exactement la même taille que tes vecteur x et y (?).

Là :

Code :

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

1
2
3
    suspe=x(~strat & ~tropo);
    xsus=x(suspe);
    ysus=y(suspe);

suspe devient un sous-ensemble de x.
Si tu veux que suspe correspondent aux indices, il faut écrire :

Code :

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

1
2
3
    suspe=~strat & ~tropo;
    xsus=x(suspe);
    ysus=y(suspe);

[j_ere_my]

ok

derniere question ....
je veux donc écrire dans un fichier les colonnes, au nombre de 100, qui correspondent au sous-ensemble "suspe"
quand je fais mon fprintf , cela me met tous dans le désordre, en rassemblant les mêmes valeurs

Code :

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

1
2
3
fprintf(fid,'%f %f %f %f %f  ....\n,x1((suspe),1),x1((suspe),2),x1((suspe),3),x1((suspe),4),x1((suspe),5),x1((su.......
fclose(fid);
!!!! x1 étant la matrice de départ (n,100)

et je pense que j'en aurais fini avec mes questions

en tout cas, merci, je vois un peu mieux comment marche MATLAB si l'on opère en vecteur et pas en boucle

[Jerome Briot]

Jette un oeil à la fonction DLMWRITE

Sinon, dans ton cas, il aurait fallut faire :

Code :

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

1
2
3
4
fmt = [repmat('%f ',1,99) '%f\n'];
fid = fopen('data.txt','wt');
   fprintf(fid, fmt, x1(suspe,:).');
fclose(fid);

[Caro-Line]

Ne me dit pas que pour ton fprintf tu as écrit tous tes x1(suspe,i) ??

Un conseil :

Code :

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

1
2
3
4
%Création du format
n=100;%on veut 100 colonnes de float
str=repmat('%f ',1,100);
fprintf(fid,[str,'\n'],x1(suspe,:));

devrait fonctionner (à moins que x1 ait plus de 100 colonnes ?).

Par contre je ne comprends pas ce que tu entends par "dans le désordre, en rassemblant les mêmes valeurs".
MATLAB écrit ce que tu lui demandes d'écrire, donc peut-être que tes x1(suspe,i) ne correspondent pas à ce que tu attendais.

[j_ere_my]

bon c'est génial, ca marche

euh non j'ai pas écrit toutes les colonnes

merci beaucoup

[Jerome Briot]

Par contre je ne comprends pas ce que tu entends par "dans le désordre, en rassemblant les mêmes valeurs".

Code :

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

1
2
3
4
5
6
7
X = [1 2 3
     1 2 3
     1 2 3]

fprintf('%d %d %d\n',X)
fprintf('\n')
fprintf('%d %d %d\n',X.')

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
X =
 
     1     2     3
     1     2     3
     1     2     3
 
1 1 1
2 2 2
3 3 3
 
1 2 3
1 2 3
1 2 3