Discussion :

[Ptinéwik]

Ola, bonjour à tous!!!

Alors voila je vous explique ce que je veux faire....
J'ai un fichier de référence M1=12*4 avec en colonnes:
1-Longitude
2-Latitude
3-Distance à la terre
4-Bathymétrie

Et donc pour 12 stations

J'ai une autre matrice de M2=2994*6
Je veux faire un programme qui permette de réattribuer distance à la cote et bathymétrie (que j'ai dans M1) à M2

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
clear all
close all
 
load M1
load M2
 
[n,p]=size(M1);
[n2,p2]=size(M2);
 
A=ones(n2,10);
[n3,p3]=size(A);
 
 
for i=1:n
 
    clear f
    f=find(M2(:,1)==M1(i,2));
 
 
     cp=0;
     for i2=1:n2  
 
     cp=cp+1;  
     [nf,pf]=size(f); 
     A(i2,:)=[AmmoniumSOMLITbis(f,:) repmat(REF(i,:),nf,1)];
 
     end
end

Mais j'ai le message

Code :

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

??? Subscripted assignment dimension mismatch.

Donc là, je trouve dans M2 toutes les lignes où j'ai la valeur correspondante de M1.

Ce que je veux lui faire faire, c'est créer repmat qui est la répétition des données de M1 autant de fois que trouvé dans M2, puis de concaténer M1 et M2.
Mais le pb qui se pose est que cela ne fonctionne que pour une la dernière valeur de i...je suis dedans depuis hier et je commence à avoir des lignes de code qui sortent de mes oreilles....

[Ptinéwik]

Bon je vais faire un message plus simple

Donc voila 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
 
clear all
close all
 
load M1
load M2
 
[n,p]=size(M1);
[n2,p2]=size(M2);
 
cp=0;
for i=1:n
    cp=cp+1;
    clear f
    f=find(M2(:,1)==M1(i,2));
    [nf,pf]=size(f);
    A=repmat(M1(cp,:),nf,1);
 
    B=                      ;
 
end

Je veux que B soit une matrice donnant la concaténation horizontale de A pour chaque valeur de i...

Au secours

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
cp=0;
B=[];
for i=1:n
    cp=cp+1;
    clear f
    f=find(M2(:,1)==M1(i,2));
    [nf,pf]=size(f);
    A=repmat(M1(cp,:),nf,1);

    B= [B A];

end