Discussion :

[Newenda]

Bonjour

Alors voilà,n je sollicite votre aide pour un petit problème que je vais essayé de simplifier :

J'ai une matrice DATA de valeurs, disons que chaque colonne de la matrice correspond a plusieurs valeurs d'un jour (vecteur Jour).

Code :

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

1
2
3
4
Jour = 1        2         3          5         9
DATA = 5        8          4          7        5
      1        4          7          5        5
      7        8          4          4        6

Bon j'essaye de ploter ça avec PCOLOR seulement il me garde la même valeur entre le jour 3 et le jour 5, c'est à dire qu'il me met la valeur du jour 3 dans jour 4 alors que là justement il n'y a pas de valeur.

Donc voici ma question, comment faire pour refaire une matrice DATA de valeur ou il y aura des 0 dans les colonnes ou les jours n'existeront pas.

Déjà recréer le vecteur jour afin d'avoir un vecteur avec TOUS les jours, même ceux n'ayant pas de mesure, par :

Code :

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

jour_tous = jour(1):jour(end)

Mais ensuite je bloque pour recréer une matrice DATA avec des colonnes de 0 là où les jours ont été ajoutés. Peut etre en créant une matrice DATA de 0 de la bonne dimension et ajouter les valeurs là où il faut mais j'ai besoin d'aide là.

Merci

[Jerome Briot]

Mais ensuite je bloque pour recréer une matrice DATA avec des colonnes de 0 là où les jours ont été ajoutés.

Simplement comme ceci :

Code :

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

1
2
3
4
5
6
7
8
Jour = [1        2         3          5         9]
 
DATA = [5        8          4          7        5
        1        4          7          5        5
        7        8          4          4        6]
 
 
A(:,Jour) = DATA

[christophe_halgand]

Prends les numéros des jours comme des indices de colonnes c'est à dire un truc du style :

Code :

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

1
2
3
4
5
for a=1:max(jour)
    if jour(a)~=a
        DATA(:,a)=zeros(nb_line,nb_col)
    end
end

C'est l'idée que tu dois voir... lol

Christophe

PS : trop rapide Dut...

[Newenda]

Simplement comme ceci :

Code :

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

1
2
3
4
5
6
7
8
Jour = [1        2         3          5         9]
 
DATA = [5        8          4          7        5
        1        4          7          5        5
        7        8          4          4        6]
 
 
A(:,Jour) = DATA

Merci bien mais c'est là qu'on remarque que trop simplifier le problème n'est pas toujours opportun... En réalité mes jours ne sont pas 1,2,3 etc.. mais plutôt sous cette forme (AnsMoisJour):
20090201
20090202
20090203
20090205
20090209

Donc là il va chercher une matrice à plus de 20 millions de valeurs...

Je peux aussi avoir du temps en Heure, voire en minute mais toujours sous cette forme ou le pas de temps correspondra à +1
(exemple en heure :
2009020100
2009020104
2009020106
2009020203
2009020209)

Prends les numéros des jours comme des indices de colonnes c'est à dire un truc du style :

Code :

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

1
2
3
4
5
for a=1:max(jour)
    if jour(a)~=a
        DATA(:,a)=zeros(nb_line,nb_col)
    end
end

C'est l'idée que tu dois voir... lol

Christophe

PS : trop rapide Dut...

Je vais étudier ça, mais j'aimerais éviter les boucles autant que faire ce peut.

En tout cas merci à vous.

[Jerome Briot]

En réalité mes jours ne sont pas 1,2,3 etc.. mais plutôt sous cette forme (AnsMoisJour):
20090201
20090202
20090203
20090205
20090209

Donc là il va chercher une matrice à plus de 20 millions de valeurs...

Tu te rappelles de la fonction DATENUM ? Mais si, tu t'en rappelles !

Code :

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

1
2
3
4
5
6
7
8
9
D = ['20090201'
     '20090202'
     '20090203'
     '20090205'
     '20090209']
 
Jour = datenum(D,'yyyymmdd');
 
Jour = Jour-Jour(1)+1

Encore une fois, c'est l'idée... à toi de l'adapter à ton problème

[Newenda]

Jour = datenum(D,'yyyymmdd');

Oui ça e le fais juste après, car je plot en temps réel donc j'ai besoin de ce temps matlab.

Jour = Jour-Jour(1)+1[/CODE]

Oui tient pas con ça !

Mais en faite je me prends la tête parce qu'à la base, j'ai fait un méthode de calcul de moyennage de mes valeurs en fonction du temps peu orthodoxe:

******************************************************

En effet, mes données brutes sont de cette forme

Code :

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

1
2
3
4
5
6
 
temps               bf1    bf2     bf3     bf4     bf5
733581.93072453700 ; 5     8     6     4     5 
733581.93087268516 ; 4     7     5     6    4
733581.93132986105 ; 4     7     5     6     1
[...]

Avec en 1ere colonne le temps et ensuite les différentes valeurs (correspond à la valeur d'un bande fréquentiel: bf). Moi je transforme le temps pour avoir une valeur identique par pas que je souhaite (Jour, Heure, Minute). Si je veux une moyenne par heure par exemple :

Code :

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

1
2
3
4
5
temps_pas = datestr(temps,'yyyymmddHH')
temps_pas = str2num(temps_pas); % en valeur numérique

(%% pour un pas de temps en minute il suffirait d'écrire:
temps_pas = datestr(temps,'yyyymmddHHMM'))

Ensuite j'utilise la fonction UNIQUE pour moyenner toutes les valeurs comprises dans une heure

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
temps_heure_uni = unique(temps_pas) ;
 
 
for i = 1:length(temps_heure_uni);
 
    bf1_heure(i) = mean(bf1(temps_pas==temps_heure_uni(i)));
    bf2_heure(i) = mean(bf2(temps_pas==temps_heure_uni(i)));
    bf3_heure(i) = mean(bf3(temps_pas==temps_heure_uni(i)));
    bf4_heure(ii) = mean(bf3(temps_pas==temps_heure_uni(i)));
    bf5_heure(i) = mean(bf3(temps_pas==temps_heure_uni(i)));
 
end

pour finir je rechange le format de mon temps pour le plottage :

Code :

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

temps_heure_uni_matlab = datenum(num2str(temps_pas_uni),'yyyymmdd');

et je plot


******************************************************

Ce qui serait beaucoup plus simple c'est que je moyenne à TOUTES les heures directement et pas seulement à celles qui ont des valeurs comme ici.

Peut etre est ce simple et cela résoudrait tout mon problème d'une manière beaucoup plus efficace et simple.

Une sorte de boucle come cela :

Code :

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

1
2
3
for i = début_temps:1_heure:fin_temps
   moyenne_heure(i)=mean(data(temps>i & temps<i+1))
end

Merci

[Newenda]

Bon j'ai reprit l'idée de DUT, ça marche mais finalement je vais essayer de moyenner par pas:

Code :

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

1
2
3
4
5
6
7
 
 
j=1 ;
for i = début_temps : pas_temps : fin_temps
   moyenne_heure(j)=mean(data(temps>i & temps<i+1))
j=j+1
end

ça marche apparemment et ça me met des NaN dans les pas de temps sans valeurs

Merci