Discussion :

[Cypette]

Salut,

J'ai un problème avec une matrice A: Je dois faire la moyenne pour les cellules de chaque chaque colonne en ignorant les valeurs égales à 0.
Mon idée était d'utiliser nanmean, c'est à dire de convertir les 0 en NaN puis d'utiliser nanmean.
Mon problème est que j'utilise un tableau de cellule en argument, donc mon code:

Code :

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

1
2
A(A==0)=Nan;
nanmean(A);

ne marche pas. Comment puis-je faire?

Merci beaucoup

Cypette

[Invité]

Bonjour,

J'ai un problème avec une matrice A

j'utilise un tableau de cellule en argument

Qu'utilises-tu exactement ? Aurais-tu un exemple ?

[Cypette]

Je vais essayer d'être clair !

Sur mon fichier, j'ai des valeurs toutes les 3 min pour un certains nombre de jour. Je dois maintenant pour les heures entre 6h et 7h; 7h et 8h , 9h et 10h ainsi que pour 15h et 16h/16h et 17h/17h et 18h/18h et 19h/19h et 20h faire la moyenne des ces valeurs sur l'heure entière. Mon problème est que je ne veux pas que dans la moyenne les 0 entrent en compte.

Sur mon code j'ai d'abord créée une matrice B: B=[nombre de lignes à prendre en compte pour la moyenne, num de ligne du début du bloc à prendre en compte, date début, heure]

Ensuite je fais la moyenne pour chaque bloc de 1 h.

C'était pas très compliqué puisqu'il y a toujours 480 minutes dans une heure.

Voici mon fichier de base radar511PLSDt.txt

Voici 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
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
function [ ] = traitementbisbisPL4( )
tic
fmt = '%d\t%d\t%d\t%d\t%d\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\n';
 
beta = 'DonneesAgregees/radar511PLSD.txt';
 
fid = fopen(beta);
A = textscan(fid, fmt);
fclose(fid);
 
fidWrite = fopen('DonneesAgregees/radar511PLSDAGR.txt', 'w+');
 
total_size = length(A{1});
 
n=121;
i=1;
count=0;
 
while (n<=total_size)
 
    index=n;
 
    while (index <= n+19)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=6;
    count=0;
    i=i+1;
 
    while (index <= n+39)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=7;
    count=0;
    i=i+1;
 
    while (index <= n+59)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=8;
    count=0;
    i=i+1;
 
    while (index <= n+79)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=9;
    count=0;
    i=i+1;
 
    index=index+100;
 
     while (index <= n+199)
 
        count=count+1;
        index=index+1;
 
     end
 
     B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=15;
    count=0;
    i=i+1;
 
    while (index <= n+219)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=16;
    count=0;
    i=i+1;
 
    while (index <= n+239)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=17;
    count=0;
    i=i+1;
 
    while (index <= n+259)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=18;
    count=0;
    i=i+1;
 
    while (index <= n+279)
 
        count=count+1;
        index=index+1;
 
    end
 
    B{i,1}=count-1;
    B{i,2}=index-20;
    B{i,3}=A{2}(index-20);
    B{i,4}=19;
    count=0;
    i=i+1;
 
    n=n+480;
 
end
 
n=length(B);
 
for j=1:n
 
    C{j,1}=sum(A{5}(B{j,2}:B{j,1}+B{j,2}));
    C{j,2}=mean(A{6}(B{j,2}:B{j,1}+B{j,2}));
    C{j,3}=mean(A{7}(B{j,2}:B{j,1}+B{j,2}));
    C{j,4}=mean(A{8}(B{j,2}:B{j,1}+B{j,2}));
    C{j,5}=mean(A{9}(B{j,2}:B{j,1}+B{j,2}));
    C{j,6}=mean(A{11}(B{j,2}:B{j,1}+B{j,2}));
    C{j,7}=mean(A{12}(B{j,2}:B{j,1}+B{j,2}));
    C{j,8}=mean(A{13}(B{j,2}:B{j,1}+B{j,2}));
    C{j,9}=mean(A{14}(B{j,2}:B{j,1}+B{j,2}));
    C{j,10}=mean(A{15}(B{j,2}:B{j,1}+B{j,2}));
    C{j,11}=mean(A{16}(B{j,2}:B{j,1}+B{j,2}));
    C{j,12}=mean(A{17}(B{j,2}:B{j,1}+B{j,2}));
    C{j,13}=mean(A{18}(B{j,2}:B{j,1}+B{j,2}));
    C{j,14}=mean(A{19}(B{j,2}:B{j,1}+B{j,2}));
    C{j,15}=mean(A{20}(B{j,2}:B{j,1}+B{j,2}));
    C{j,16}=mean(A{21}(B{j,2}:B{j,1}+B{j,2}));
    C{j,17}=mean(A{22}(B{j,2}:B{j,1}+B{j,2}));
    C{j,18}=mean(A{23}(B{j,2}:B{j,1}+B{j,2}));
    C{j,19}=sum(A{24}(B{j,2}:B{j,1}+B{j,2}));
    C{j,20}=mean(A{25}(B{j,2}:B{j,1}+B{j,2}));
    C{j,21}=mean(A{26}(B{j,2}:B{j,1}+B{j,2}));
    C{j,22}=mean(A{27}(B{j,2}:B{j,1}+B{j,2}));
    C{j,23}=mean(A{28}(B{j,2}:B{j,1}+B{j,2}));
    C{j,24}=mean(A{30}(B{j,2}:B{j,1}+B{j,2}));
    C{j,25}=mean(A{31}(B{j,2}:B{j,1}+B{j,2}));
    C{j,26}=mean(A{32}(B{j,2}:B{j,1}+B{j,2}));
    C{j,27}=mean(A{33}(B{j,2}:B{j,1}+B{j,2}));
    C{j,28}=mean(A{34}(B{j,2}:B{j,1}+B{j,2}));
    C{j,29}=mean(A{35}(B{j,2}:B{j,1}+B{j,2}));
    C{j,30}=mean(A{36}(B{j,2}:B{j,1}+B{j,2}));
    C{j,31}=mean(A{37}(B{j,2}:B{j,1}+B{j,2}));
    C{j,32}=mean(A{38}(B{j,2}:B{j,1}+B{j,2}));
    C{j,33}=mean(A{39}(B{j,2}:B{j,1}+B{j,2}));
    C{j,34}=mean(A{40}(B{j,2}:B{j,1}+B{j,2}));
    C{j,35}=mean(A{41}(B{j,2}:B{j,1}+B{j,2}));
    C{j,36}=mean(A{42}(B{j,2}:B{j,1}+B{j,2}));
 
    fprintf(fidWrite,'%d\t%d\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\n',B{j,3},B{j,4},C{j,1},C{j,2},C{j,3},C{j,4},C{j,5},C{j,6},C{j,7},C{j,8},C{j,9},C{j,10},C{j,11},C{j,12},C{j,13},C{j,14},C{j,15},C{j,16},C{j,17},C{j,18},C{j,19},C{j,20},C{j,21},C{j,22},C{j,23},C{j,24},C{j,25},C{j,26},C{j,27},C{j,28},C{j,29},C{j,30},C{j,31},C{j,32},C{j,33},C{j,34},C{j,35},C{j,36});
 
end
 
toc

[Invité]

Code :

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

fmt = '%d\t%d\t%d\t%d\t%d\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f';

Tu peux aussi spécifier textscan(..., 'delimiter', '\t').
Tes boucles while pourraient se simplifier en une simple différence.

Pour le reste, le plus simple serait de te ramener à un tableau numérique à la place d'un tableau de cellules.
Rajouter l'option 'CollectOutput', true à textscan sera un bon début.

[Cypette]

Excuse moi je n'ai pas très bien saisi ton message, notamment comment fais-t-on exactement pour employer CollectOutPut dans mon code?

Merci bcp

[Invité]

Idem que délimiter :

Code :

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

textscan(..., 'CollectOutput', true)

[Cypette]

Je ne comprends pas bien parce que ça me découpe mon tableau en plusieurs sous-tableaux (cad A{1} est une matrice avec les 5 premières colonnes A{2} avec les 4 prochaines colonnes etc...). De plus pour A{2} cela me divise mes valeurs par 10^4...

[Invité]

Cela rassemble les colonnes par formats adjacents communs :

Code :

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

fmt = '%d\t%d\t%d\t%d\t%d\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%d\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f\t%f';

Le mieux serait selon moi de n'utiliser qu'un seul et même format, cela te faciliterait grandement les choses

[Cypette]

Pour tester le code lorsque je mets :

Code :

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

A(A{1}(:,5)==0)=NaN

Cela m'indique le même message d'erreur qu'avant:

>> traitementbisbisPL4
Conversion to cell from double is not possible.

Error in traitementbisbisPL4 (line 13)
A(A{1}(:,5)==0)=NaN

[Invité]

Code :

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

A{1}((A{1}(:,5)==0),5)=NaN

[Cypette]

Merci bcp et une dernière question: pour remplacer tt les 0 par NaN sauf les 0 de la 6 eme et 9eme colonne cmt puis-je faire?

[Invité]

Tu peux sélectionner les colonnes comme ceci : [1:5 7 8 10:end]