Discussion :

[Ptinéwik]

Bonjour, alors voila qq lignes du programme que j'a écrit.
Le but est de renvoyer les résultats de la boucle Data2 dans la matrice Datatout.
En fait, il s'agit du calcul des moyennes mensuelles de 13 paramètres sur 12 stations, en sachant que chaque station est prise au début de la boucle.

Le problème est que la concaténation ne se fait pas comme il faut, et que je me retrouve avec juste ma 13ème station de mémorisée...

Quelqu'un pourrait-il m'aider?

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
cp0=-12;
for k=1:12
    cp0=cp0+13;
    S=load(char(Site(k,3)));
    Data=S.(char(Site(k,3)));
    Data=[Data(:,6:7) Data(:,13:end)];
    clear S
 
    %if TEMP==1
 
    Data2=[];
    Datatout=[];
 
    cp=0;
    for i=1:12
        clear f
        cp=cp+1;
        f=find(Data(:,2)==i);
        if ~isempty(f)
            Data2(cp,:)=nanmean(Data(f,3:end));
 
        end
    end
 
    Datatout(:,(cp0:cp0+12))=Data2;
end

[Caro-Line]

Déjà dès tes 1ères lignes tu écrases ta variable Data :

Code :

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

Data=[Data(:,6:7) Data(:,13:end)];

Tu ne gardes que les 6,7 et 13 jusqu'à la fin colonnes du Data que tu viens de calculer, à la prochaine itération cette variable sera écrasée : est-ce ce que tu souhaites ?
Après j'avoue que sans données pour tester c'est un peu difficile de savoir ce que tu as en entrée et ce que tu veux en sortie.

[EDIT]Ah oui et surtout tu réinitialises Datatout à vide [] à chaque itération donc forcément...Cette initialisation doit être faite avant la boucle for principale.

[Ptinéwik]

Et oui, c'était cet écrasement qui venait tout vider au fur et à mesure...trop bete, mais à force d'etre dessus, je n'ai pas vu cette petite erreur....merci beaucoup de ton aide...comme d'habitude d'ailleurs....