Discussion :

[niepoc]

bonjour tout le monde

j'ai deux matrices que je souhaiterais comparer par ligne pour pouvoir garder que celles qui m'interessent, elles sont de tailles differentes en nombre de lignes et de colonnes

voilà le debut de mes deux matrices

M =



1.0e+003 *



2.0050 0.0010 0.0010 0.0100 0.0107

2.0050 0.0010 0.0020 0.0100 0.0125

2.0050 0.0010 0.0030 0.0110 0.0121

2.0050 0.0010 0.0040 0.0150 0.0156

2.0050 0.0010 0.0050 0.0130 0.0126

2.0050 0.0010 0.0060 0.0300 0.0288

2.0050 0.0010 0.0070 0.0150 0.0152

2.0050 0.0010 0.0080 0.0110 0.0144

2.0050 0.0010 0.0090 0.0160 0.0162

2.0050 0.0010 0.0100 0.0170 0.0187

2.0050 0.0010 0.0110 NaN 0.0138

2.0050 0.0010 0.0120 0.0220 0.0224

2.0050 0.0010 0.0130 0.0270 0.0244

2.0050 0.0010 0.0140 0.0250 0.0237

2.0050 0.0010 0.0150 0.0290 0.0287

2.0050 0.0010 0.0160 0.0120 0.0141

2.0050 0.0010 0.0170 0.0270 0.0238

2.0050 0.0010 0.0180 NaN 0.0168

2.0050 0.0010 0.0190 NaN 0.0120

2.0050 0.0010 0.0200 0.0110 0.0111

OT440 =



1.0e+003 *



2.0050 0.0010 0.0010 0.0002

2.0050 0.0010 0.0040 0.0007

2.0050 0.0010 0.0050 0.0009

2.0050 0.0010 0.0060 0.0026

2.0050 0.0010 0.0070 0.0010

2.0050 0.0010 0.0080 0.0008

2.0050 0.0010 0.0090 0.0013

2.0050 0.0010 0.0120 0.0028

2.0050 0.0010 0.0130 0.0009

2.0050 0.0010 0.0150 0.0005

2.0050 0.0010 0.0180 0.0005

2.0050 0.0010 0.0190 0.0002

2.0050 0.0010 0.0210 0.0002

2.0050 0.0010 0.0220 0.0016

2.0050 0.0010 0.0230 0.0010

2.0050 0.0010 0.0240 0.0006

2.0050 0.0010 0.0250 0.0009

2.0050 0.0010 0.0260 0.0009

2.0050 0.0010 0.0300 0.0010

2.0050 0.0010 0.0310 0.0006

l'idée est de comparer les trois premieres colonnes de chaque matrice et de ne garder que les lignes de M ou l'on trouve les trois premieres valeurs dans AOT440, (en fait ces trois premieres valeurs sont des dates)

j'ai essayé avec intersect, setxor, celà me donne les lignes avec les valeurs et non la position de la ligne.

j'ai du mal avec l'indexation aussi

est ce que quelqu'un a une idéé

voilà le debut de mon code

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
clear all;
 
close all;
 
 
 
%% CHARGEMENT DES DEUX MATRICES A COMPARER
 
 
 
load('matrice.mat')%données agoufou 2005 pour AOT440
 
AOT440;
 
load ('365 fichiers lus, dernier au 20051231.mat')
 
MAT;
 
TABFINAL;
 
M=MAT(:,1:5) %données Agoufou OMI, valeur cimel et moyenne pavé 3
 
[x y]=size(MAT);
 
[n m]=size(AOT440);
 
%% COMPARAISON DES DEUX MATRICES
 
 
 
a=M(:,1:3);
 
b=AOT440(:,1:3);
 
c=intersect(a,b)
 
 
 
for n=1:length(b)
 
    for m=1:length(a)
 
    if b(n,:) == a(m,:)

aprés je bloque

merci d'avance

[Caro-Line]

Regarde un peu mieux l'aide de la fonction INTERSECT : tu peux avoir plus que d'un paramètre de sortie.

[c, ia, ib] = intersect(a, b) also returns column index vectors ia and ib such that c = a(ia) and c = b(ib) (or c = a(ia,: ) and c = b(ib,: )).

En combinant avec l'option 'rows' tu devrait t'approcher du résultat.

[niepoc]

merci j'ai reussi à trouver mes numeros de lignes avec ce code suivant

Code :

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

1
2
3
4
5
a=M(:,1:3);
 
b=AOT440(:,1:3);
 
[c, ia, ib] = intersect(a, b,'rows');

maintenant je dois affecter dans une nouvelle matrice les lignes que je garde et là ça pose probleme

voilà le code

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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
clear all;
 
close all;
 
 
 
%% CHARGEMENT DES DEUX MATRICES A COMPARER
 
 
 
load('matrice.mat')%données agoufou 2005 pour AOT440
 
AOT440;
 
load ('365 fichiers lus, dernier au 20051231.mat')
 
MAT;
 
TABFINAL;
 
M=MAT(:,1:5) %données Agoufou OMI, valeur cimel et moyenne pavé 3..........mat=[M(:,1:3), M(:,5)] sans la valeur cimel
 
[x y]=size(MAT);
 
[n m]=size(AOT440);
 
%% COMPARAISON DES DEUX MATRICES
 
 
 
%mettre les matrices avec le même nombre de colonnes pour trouver les
 
%lignes à conserver, valeur donnée par 'ia'
 
a=M(:,1:3);
 
b=AOT440(:,1:3);
 
[c, ia, ib] = intersect(a, b,'rows');
 
 
 
 
 
ligne=ia'; %ia devient un vecteur
 
L=length(ligne) %taille du vecteur
 
agou05=zeros(L,6) %affectation de la nouvelle matrice
 
 
 
for n=1:L;
 
    agou05(n,:)=M(ligne(n),:);
 
    agou05(n,6)=AOT440(n,4);
 
end

et le message d'erreur est le suivant

??? Error using ==> unknown

Subscripted assignment dimension mismatch.



Error in ==> comparaison at 28

agou05(n,:)=M(ligne(n),:);

je comprends pas je cherhce dans tous les sens je vois pas où est l'erreur

merci

[niepoc]

c bon j'ai trouvé en insistant

voilà la fin du code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
ligne=ia'; %ia devient un vecteur
 
L=length(ligne) %taille du vecteur
 
 
 
for n=1:L;
 
    agou05(n,:)=M(ligne(n),:);
 
end
 
agou05=[agou05, AOT440(:,4)]

[Jerome Briot]

Ou plus simplement :

Code :

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

1
2
[c, ia, ib] = intersect(M(:,1:3), AOT440(:,1:3),'rows');
agou05 = M(ia,:);

Ou encore :

Code :

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

1
2
idx = ismember(M(:,1:3),AOT440(:,1:3),'rows');
agou05 = M(idx,:)