Discussion :

[Bonois]

Bonjour;

J'ai une longue matrice, dans la première colone des dates a des pas de temps variable sous le format yyyymmdd HH:MM et dans la deuxiemme colonne des valeurs correspondantes.

je voudrais crié une nouvelle matrice avec pour chaque journnée la valeurs max. c'est a dire prendre les valeurs maximuim journalière mais je ne sais pas comment m'y prendre. autre chose quand il y a deux fois la même valeur max , prendre qu'une seule.

est que vous pouvez me donner une piste.

merci
exemple:

la matrice initiale

19801017 00:00 0.407
19801017 01:24 0.577
19801017 02:16 0.917
19801017 02:26 1.01
19801017 04:55 0.917
19801017 09:01 0.728
19801017 12:05 0.623
19801017 18:59 0.577
19801019 00:43 0.49
19801019 04:49 0.577
19801019 18:51 0.49
19801020 14:36 0.452
19801020 16:06 0.413
19801021 16:27 0.413
19801021 18:43 0.375
19801023 08:25 0.344
19801023 10:16 0.375
19801023 14:16 0.375
19801023 20:42 0.344
19801024 09:53 0.314
19801024 11:22 0.344
19801024 13:14 0.344
19801024 13:42 0.375
19801024 16:01 0.577
19801024 17:27 0.917

la matrice que je voudrais avoir:

19801017 02:26 1.01
19801019 04:49 0.577
19801020 14:36 0.452
19801021 16:27 0.413
19801023 14:16 0.375
19801024 17:27 0.917

[Invité]

Bonjour,

Je pense que tu peux t'en sortir en utilisant la fonction

Code :

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

doc unique

sur le vecteur contenant les dates et utilisant le 'm' de la documentation Matlab

[Bonois]

Merci , ca marche voici le programme que j'ai fait

Code :

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

1
2
3
4
5
6
7
8
9
10
[A,C]=textread('QCrappone_80_08.txt','%15c%f');
B= datenum(A,'yyyymmdd');
D = unique (B)
 
for i=1:length(D) % nombre de jours dans la période 1980-2008
    rr = D(i)
    idx=find(B == rr);
    G(i,: ) = max (C(idx));
end
Qmax=[D,G]

merci encore

[Jerome Briot]

Une solution sans boucle :

Code :

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

1
2
3
4
5
6
7
8
9
[A,B,C] = textread('QCrappone_80_08.txt','%s %s %f');
 
An  = datenum(A,'yyyymmdd');
 
[AC,idx] = sortrows([An C]);
 
[pasbesoin,idx2,pasbesoin] = unique(AC(:,1),'last');
 
strcat(A(idx(idx2)),{' '},B(idx(idx2)),{' '},num2str(C(idx(idx2)),'%-.3f'))

et une autre sur le même principe :

Code :

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

1
2
3
4
5
6
7
[A,B,C] = textread('QCrappone_80_08.txt','%s %s %s');
 
[AC,idx] = sortrows(strcat(A,C));
 
[pasbesoin,idx2,pasbesoin] = unique(A,'last');
 
strcat(A(idx(idx2)),{' '},B(idx(idx2)),{' '},C(idx(idx2)))

[Bonois]

C'est un coup de maître