Discussion :

[diego45]

Je suis totalement d'accord.
Et le code ne tourne que partiellement, je comprend pas pourquoi cela écrase les valeurs et reste bloqué à 6.


Donc pour en revenir à ton message précedent, j'ai lancé le code, et j'obtiens bien effectivement l'ensemble de mes chemins dans une structure 1x1445 cell.

Après je vais pouvoir les appeler plus facilement comme ceci :

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
 
 
for k = 1:numel(chemin_spread2)
 
 
    o = 0;
 
 
    fid = fopen(chemin_spread2{k,:}, 'r');
    M = textscan(fid, '%s %*f', 'TreatAsEmpty', '.', 'HeaderLines', 1);
    fclose(fid);
    formatIn = 'dd/mm/yy';
    J = datenum(M{1},formatIn);
 
    udate = unique([udate ; J]);
 
    o = o+1;
 
    C = NaN(numel(udate), o+1);
    C(:,1) = udate;
 
    o = 1;
 
    fid = fopen(chemin_spread2{k,:}, 'r');
    M = textscan(fid, '%s %f', 'TreatAsEmpty', '.', 'HeaderLines', 1);
    fclose(fid);
 
    formatIn = 'dd/mm/yy';
    J = datenum(M{1},formatIn);
 
    idx = ismember(udate, J);
    o = o+1;
    C(idx,o) = M{1,2};
 
end

[Jerome Briot]

OK donc maintenant deuxième partie, récupérer les dates uniques :

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
udate = [];
formatIn = 'dd/mm/yy';
 
for k = 1:numel(chemin_spread2)      
 
    fid = fopen(chemin_spread2{k}, 'r');
    M = textscan(fid, '%s %*f', 'TreatAsEmpty', '.', 'HeaderLines', 1);
    fclose(fid);
 
    J = datenum(M{1}, formatIn);
 
    udate = [udate ; J];
 
end
 
udate = unique(udate);

Le vecteur udate contient-il toutes les dates uniques ?

[diego45]

effectivement, il contient l'ensemble des dates uniques.

Et donc normalement, en faisant le code suivant, je devrais avoir ce que je cherche :

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
 
 
C = NaN(numel(udate),numel(chemin_spread2));
o = 1;
C(:,1) = udate;
for k = 1:numel(chemin_spread2)
    fid = fopen(chemin_spread2{k,:}, 'r');
    M = textscan(fid, '%s %f', 'TreatAsEmpty', '.', 'HeaderLines', 1);
    fclose(fid);
    o=o+1;
    formatIn = 'dd/mm/yy';
    J = datenum(M{1},formatIn);
 
    idx = ismember(udate, J);
 
    C(idx,o) = M{1,2};
 
end

Si pour toi le code est bon, chez moi ça à l'air de tourner ça prend à peine 2 minutes VS des heures avant !!! Un énorme merci

[Jerome Briot]

Tu as compris. Pour le reste, tu es le seul à pouvoir vérifier que le résultat final est le bon.

Pour finir, tu peux améliorer la lisibilité du code en mettant par exemple l'incrémentation de la variable o juste avant de l'utiliser :

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
C = NaN(numel(udate),numel(chemin_spread2));
C(:,1) = udate;

o = 1;
for k = 1:numel(chemin_spread2)
    
    fid = fopen(chemin_spread2{k}, 'r');
    M = textscan(fid, '%s %f', 'TreatAsEmpty', '.', 'HeaderLines', 1);
    fclose(fid);

    formatIn = 'dd/mm/yy';
    J = datenum(M{1},formatIn);
    
    idx = ismember(udate, J);
    
    o = o+1;
    C(idx,o) = M{2};
    
end