traitement lecture fichier et plot figure

Version imprimable

Bonjour,
Pour simplifier mon code et tracer toutes mes données, voilà mon code :

Code:

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
clear all
close all
clc
load gsom.mat;
idxLignesZeros = find(gsom(:,1) == 0);
sous_matrices = mat2cell(gsom, diff([idxLignesZeros ; size(gsom,1)+1]), size(gsom,2));
TimeStep=0.1;
for i=1:60
      eval(['v' num2str(i) '= sous_matrices(i,1)']);%v1=sous_matrices{1,1};
      eval(['T' num2str(i) '=[]']);%T1=[];
end

for t = v1{1,1}(:,4):TimeStep:(length(v1{1,1})-1+(v1{1,1}(:,4)/TimeStep))*TimeStep;
    T1=[T1 t];
end

for t = v2{1,1}(:,4):TimeStep:(length(v2{1,1})-1+(v2{1,1}(:,4)/TimeStep))*TimeStep;
    T2=[T2 t];
end
for t = v3{1,1}(:,4):TimeStep:(length(v3{1,1})-1+(v3{1,1}(:,4)/TimeStep))*TimeStep;
    T3=[T3 t];
end
for t = v4{1,1}(:,4):TimeStep:(length(v4{1,1})-1+(v4{1,1}(:,4)/TimeStep))*TimeStep;
    T4=[T4 t];
end
for t = v5{1,1}(:,4):TimeStep:(length(v5{1,1})-1+(v5{1,1}(:,4)/TimeStep))*TimeStep;
    T5=[T5 t];
end
for t = v6{1,1}(:,4):TimeStep:(length(v6{1,1})-1+(v6{1,1}(:,4)/TimeStep))*TimeStep;
    T6=[T6 t];
end
for t = v7{1,1}(:,4):TimeStep:(length(v7{1,1})-1+(v7{1,1}(:,4)/TimeStep))*TimeStep;
    T7=[T7 t];
end
for t = v8{1,1}(:,4):TimeStep:(length(v8{1,1})-1+(v8{1,1}(:,4)/TimeStep))*TimeStep;
    T8=[T8 t];
end
for t = v9{1,1}(:,4):TimeStep:(length(v9{1,1})-1+(v9{1,1}(:,4)/TimeStep))*TimeStep;
    T9=[T9 t];
end
for t = v10{1,1}(:,4):TimeStep:(length(v10{1,1})-1+(v10{1,1}(:,4)/TimeStep))*TimeStep;
    T10=[T10 t];
end
for t = v11{1,1}(:,4):TimeStep:(length(v11{1,1})-1+(v11{1,1}(:,4)/TimeStep))*TimeStep;
    T11=[T11 t];
end
for t = v12{1,1}(:,4):TimeStep:(length(v12{1,1})-1+(v12{1,1}(:,4)/TimeStep))*TimeStep;
    T12=[T12 t];
end
for t = v13{1,1}(:,4):TimeStep:(length(v13{1,1})-1+(v13{1,1}(:,4)/TimeStep))*TimeStep;
    T13=[T13 t];
end
for t = v14{1,1}(:,4):TimeStep:(length(v14{1,1})-1+(v14{1,1}(:,4)/TimeStep))*TimeStep;
    T14=[T14 t];
end
for t = v15{1,1}(:,4):TimeStep:(length(v15{1,1})-1+(v15{1,1}(:,4)/TimeStep))*TimeStep;
    T15=[T15 t];
end
for t = v16{1,1}(:,4):TimeStep:(length(v16{1,1})-1+(v16{1,1}(:,4)/TimeStep))*TimeStep;
    T16=[T16 t];
end
for t = v17{1,1}(:,4):TimeStep:(length(v17{1,1})-1+(v17{1,1}(:,4)/TimeStep))*TimeStep;
    T17=[T17 t];
end
for t = v18{1,1}(:,4):TimeStep:(length(v18{1,1})-1+(v18{1,1}(:,4)/TimeStep))*TimeStep;
    T18=[T18 t];
end
for t = v19{1,1}(:,4):TimeStep:(length(v19{1,1})-1+(v19{1,1}(:,4)/TimeStep))*TimeStep;
    T19=[T19 t];
end
for t = v20{1,1}(:,4):TimeStep:(length(v20{1,1})-1+(v20{1,1}(:,4)/TimeStep))*TimeStep;
    T20=[T20 t];
end
for t = v21{1,1}(:,4):TimeStep:(length(v21{1,1})-1+(v21{1,1}(:,4)/TimeStep))*TimeStep;
    T21=[T21 t];
end
for t = v22{1,1}(:,4):TimeStep:(length(v22{1,1})-1+(v22{1,1}(:,4)/TimeStep))*TimeStep;
    T22=[T22 t];
end
for t = v23{1,1}(:,4):TimeStep:(length(v23{1,1})-1+(v23{1,1}(:,4)/TimeStep))*TimeStep;
    T23=[T23 t];
end
for t = v24{1,1}(:,4):TimeStep:(length(v24{1,1})-1+(v24{1,1}(:,4)/TimeStep))*TimeStep;
    T24=[T24 t];
end
for t = v25{1,1}(:,4):TimeStep:(length(v25{1,1})-1+(v25{1,1}(:,4)/TimeStep))*TimeStep;
    T25=[T25 t];
end
for t = v26{1,1}(:,4):TimeStep:(length(v26{1,1})-1+(v26{1,1}(:,4)/TimeStep))*TimeStep;
    T26=[T26 t];
end
for t = v27{1,1}(:,4):TimeStep:(length(v27{1,1})-1+(v27{1,1}(:,4)/TimeStep))*TimeStep;
    T27=[T27 t];
end
for t = v28{1,1}(:,4):TimeStep:(length(v28{1,1})-1+(v28{1,1}(:,4)/TimeStep))*TimeStep;
    T28=[T28 t];
end
for t = v29{1,1}(:,4):TimeStep:(length(v29{1,1})-1+(v29{1,1}(:,4)/TimeStep))*TimeStep;
    T29=[T29 t];
end
for t = v30{1,1}(:,4):TimeStep:(length(v30{1,1})-1+(v30{1,1}(:,4)/TimeStep))*TimeStep;
    T30=[T30 t];
end
for t = v31{1,1}(:,4):TimeStep:(length(v31{1,1})-1+(v31{1,1}(:,4)/TimeStep))*TimeStep;
    T31=[T31 t];
end
for t = v32{1,1}(:,4):TimeStep:(length(v32{1,1})-1+(v32{1,1}(:,4)/TimeStep))*TimeStep;
    T32=[T32 t];
end
for t = v33{1,1}(:,4):TimeStep:(length(v33{1,1})-1+(v33{1,1}(:,4)/TimeStep))*TimeStep;
    T33=[T33 t];
end
for t = v34{1,1}(:,4):TimeStep:(length(v34{1,1})-1+(v34{1,1}(:,4)/TimeStep))*TimeStep;
    T34=[T34 t];
end
for t = v35{1,1}(:,4):TimeStep:(length(v35{1,1})-1+(v35{1,1}(:,4)/TimeStep))*TimeStep;
    T35=[T35 t];
end
for t = v36{1,1}(:,4):TimeStep:(length(v36{1,1})-1+(v36{1,1}(:,4)/TimeStep))*TimeStep;
    T36=[T36 t];
end
for t = v37{1,1}(:,4):TimeStep:(length(v37{1,1})-1+(v37{1,1}(:,4)/TimeStep))*TimeStep;
    T37=[T37 t];
end
for t = v38{1,1}(:,4):TimeStep:(length(v38{1,1})-1+(v38{1,1}(:,4)/TimeStep))*TimeStep;
    T38=[T38 t];
end
for t = v39{1,1}(:,4):TimeStep:(length(v39{1,1})-1+(v39{1,1}(:,4)/TimeStep))*TimeStep;
    T39=[T39 t];
end
for t = v40{1,1}(:,4):TimeStep:(length(v40{1,1})-1+(v40{1,1}(:,4)/TimeStep))*TimeStep;
    T40=[T40 t];
end
for t = v41{1,1}(:,4):TimeStep:(length(v41{1,1})-1+(v41{1,1}(:,4)/TimeStep))*TimeStep;
    T41=[T41 t];
end
for t = v42{1,1}(:,4):TimeStep:(length(v42{1,1})-1+(v42{1,1}(:,4)/TimeStep))*TimeStep;
    T42=[T42 t];
end
for t = v43{1,1}(:,4):TimeStep:(length(v43{1,1})-1+(v43{1,1}(:,4)/TimeStep))*TimeStep;
    T43=[T43 t];
end
for t = v44{1,1}(:,4):TimeStep:(length(v44{1,1})-1+(v44{1,1}(:,4)/TimeStep))*TimeStep;
    T44=[T44 t];
end
for t = v45{1,1}(:,4):TimeStep:(length(v45{1,1})-1+(v45{1,1}(:,4)/TimeStep))*TimeStep;
    T45=[T45 t];
end
for t = v46{1,1}(:,4):TimeStep:(length(v46{1,1})-1+(v46{1,1}(:,4)/TimeStep))*TimeStep;
    T46=[T46 t];
end
for t = v47{1,1}(:,4):TimeStep:(length(v47{1,1})-1+(v47{1,1}(:,4)/TimeStep))*TimeStep;
    T47=[T47 t];
end
for t = v48{1,1}(:,4):TimeStep:(length(v48{1,1})-1+(v48{1,1}(:,4)/TimeStep))*TimeStep;
    T48=[T48 t];
end
for t = v49{1,1}(:,4):TimeStep:(length(v49{1,1})-1+(v49{1,1}(:,4)/TimeStep))*TimeStep;
    T49=[T49 t];
end
for t = v50{1,1}(:,4):TimeStep:(length(v50{1,1})-1+(v50{1,1}(:,4)/TimeStep))*TimeStep;
    T50=[T50 t];
end
for t = v51{1,1}(:,4):TimeStep:(length(v51{1,1})-1+(v51{1,1}(:,4)/TimeStep))*TimeStep;
    T51=[T51 t];
end
for t = v52{1,1}(:,4):TimeStep:(length(v52{1,1})-1+(v52{1,1}(:,4)/TimeStep))*TimeStep;
    T52=[T52 t];
end
for t = v53{1,1}(:,4):TimeStep:(length(v53{1,1})-1+(v53{1,1}(:,4)/TimeStep))*TimeStep;
    T53=[T53 t];
end
for t = v54{1,1}(:,4):TimeStep:(length(v54{1,1})-1+(v54{1,1}(:,4)/TimeStep))*TimeStep;
    T54=[T54 t];
end
for t = v55{1,1}(:,4):TimeStep:(length(v55{1,1})-1+(v55{1,1}(:,4)/TimeStep))*TimeStep;
    T55=[T55 t];
end
for t = v56{1,1}(:,4):TimeStep:(length(v56{1,1})-1+(v56{1,1}(:,4)/TimeStep))*TimeStep;
    T56=[T56 t];
end
for t = v57{1,1}(:,4):TimeStep:(length(v57{1,1})-1+(v57{1,1}(:,4)/TimeStep))*TimeStep;
    T57=[T57 t];
end
for t = v58{1,1}(:,4):TimeStep:(length(v58{1,1})-1+(v58{1,1}(:,4)/TimeStep))*TimeStep;
    T58=[T58 t];
end
for t = v59{1,1}(:,4):TimeStep:(length(v59{1,1})-1+(v59{1,1}(:,4)/TimeStep))*TimeStep;
    T59=[T59 t];
end
for t = v60{1,1}(:,4):TimeStep:(length(v60{1,1})-1+(v60{1,1}(:,4)/TimeStep))*TimeStep;
    T60=[T60 t];
end

figure
plot(T1,v1{1,1}(:,1),'b');hold on ; plot(T2,v2{1,1}(:,1),'b');plot(T3,v3{1,1}(:,1),'b');plot(T4,v4{1,1}(:,1),'b');plot(T5,v5{1,1}(:,1),'b');plot(T6,v6{1,1}(:,1),'b');plot(T7,v7{1,1}(:,1),'b');plot(T8,v8{1,1}(:,1),'b');plot(T9,v9{1,1}(:,1),'b');plot(T10,v10{1,1}(:,1),'b');plot(T11,v11{1,1}(:,1),'b');plot(T12,v12{1,1}(:,1),'b');plot(T13,v13{1,1}(:,1),'b');plot(T14,v14{1,1}(:,1),'b');plot(T15,v15{1,1}(:,1),'b');plot(T16,v16{1,1}(:,1),'b');plot(T17,v17{1,1}(:,1),'b');plot(T18,v18{1,1}(:,1),'b');plot(T19,v19{1,1}(:,1),'b');plot(T20,v20{1,1}(:,1),'b');plot(T21,v21{1,1}(:,1),'b');plot(T22,v22{1,1}(:,1),'b');plot(T23,v23{1,1}(:,1),'b');plot(T24,v24{1,1}(:,1),'b');plot(T25,v25{1,1}(:,1),'b');plot(T26,v26{1,1}(:,1),'b');plot(T27,v27{1,1}(:,1),'b');plot(T28,v28{1,1}(:,1),'b');plot(T29,v29{1,1}(:,1),'b');plot(T30,v30{1,1}(:,1),'b');plot(T31,v31{1,1}(:,1),'b');plot(T32,v32{1,1}(:,1),'b');plot(T33,v33{1,1}(:,1),'b');plot(T34,v34{1,1}(:,1),'b');plot(T35,v35{1,1}(:,1),'b');plot(T36,v36{1,1}(:,1),'b');plot(T37,v37{1,1}(:,1),'b');plot(T38,v38{1,1}(:,1),'b');plot(T39,v39{1,1}(:,1),'b');plot(T40,v40{1,1}(:,1),'b');plot(T41,v41{1,1}(:,1),'b');plot(T42,v42{1,1}(:,1),'b');plot(T43,v43{1,1}(:,1),'b');plot(T44,v44{1,1}(:,1),'b');plot(T45,v45{1,1}(:,1),'b');plot(T46,v46{1,1}(:,1),'b');plot(T47,v47{1,1}(:,1),'b');plot(T48,v48{1,1}(:,1),'b');plot(T49,v49{1,1}(:,1),'b');plot(T50,v50{1,1}(:,1),'b');plot(T51,v51{1,1}(:,1),'b');plot(T52,v52{1,1}(:,1),'b');plot(T53,v53{1,1}(:,1),'b');plot(T54,v54{1,1}(:,1),'b');plot(T55,v55{1,1}(:,1),'b');plot(T56,v56{1,1}(:,1),'b');plot(T57,v57{1,1}(:,1),'b');plot(T58,v58{1,1}(:,1),'b');plot(T59,v59{1,1}(:,1),'b');plot(T60,v60{1,1}(:,1),'b');grid
title('Distribution des trajectoires de 60 paquets de véhicules sorties selon le modèle GSOM') 
xlabel('temps de simulation (s)')
ylabel('position (m)')

Mon but est de rassembler tout ce qui est marqué en rouge dans une seule boucle générale au lieu de faire un traitement pour chaque sous_matrice (et pareil pour plot pour avoir ma figure) pour avoir un code général comportant toutes les données présentes dans "gsom.mat".
Merci d'avance

07/06/2014, 18h42
Invité

Bonjour,

On ne donnera jamais assez ce lien... ni ne répètera jamais assez que générer N variables au lieu d'avoir un tableau de N éléments est une mauvaise méthode.
Comment créer dynamiquement des variables nommées A1, A2, A3, ..., AN ?
07/06/2014, 20h45
ghada97
Bonsoir Winjerome,
J'ai utilisé déjà la fonction "eval" pour :
Code:

1 2 3 eval(['v' num2str(i) '= sous_matrices(i,1)']);%v1=sous_matrices{1,1}; eval(['T' num2str(i) '=[]']);%T1=[];
Mais ce qui m'embête ce sont les mêmes boucles :
Code:

1 2 3 4 5 6 7 for t = v1{1,1}(:,4):TimeStep:(length(v1{1,1})-1+(v1{1,1}(:,4)/TimeStep))*TimeStep; T1=[T1 t]; end . . .
qui se répètent 60 fois .
Je veux bien les résumer en une seule simple boucle en utilisant la fonction "eval", mais je n'y arrive pas, c'est pour ça j'ai demandé votre aide :?
07/06/2014, 21h55
Invité

Citation:

Envoyé par Winjerome

générer N variables au lieu d'avoir un tableau de N éléments est une mauvaise méthode.

As-tu compris cette partie ? Ainsi que le message que voulais faire passer le sujet de FAQ ?
Tu te retrouves avec N variables que tu dois gérer indépendamment (au lieu d'une seule, tableau de N éléments facilement manipulable en boucle) et comme tu as pu le constater la syntaxe avec eval est atroce !

De plus je ne vois absolument aucun intérêt à passer par ces variables vi au lieu d'écrire directement sous_matrices(i,1).
Et attention : v1{1,1}(:,4):TimeStep:(length(v1{1,1})-1+(v1{1,1}(:,4)/TimeStep))*TimeStep; vouloir manipuler un vecteur avec l'opérateur colon « : » n'est pas correct, il faut des scalaires.

Je suis débutante en Matlab :(:?
Je l'utilise seulement pour lire mon fichier de données "gsom.txt" et faire des représentations graphiques ..
L'intérêt des Vi est de remplir mes vecteurs Ti dans les boucles pour résultat final arriver à faire ça :

Code:

1
2
 
plot(T1,v1{1,1}(:,1),'b');hold on ; plot(T2,v2{1,1}(:,1),'b');plot(T3,v3{1,1}(:,1),'b');plot(T4,v4{1,1}(:,1),'b');plot(T5,v5{1,1}(:,1),'b');plot(T6,v6{1,1}(:,1),'b');plot(T7,v7{1,1}(:,1),'b');plot(T8,v8{1,1}(:,1),'b');plot(T9,v9{1,1}(:,1),'b');plot(T10,v10{1,1}(:,1),'b');plot(T11,v11{1,1}(:,1),'b');plot(T12,v12{1,1}(:,1),'b');plot(T13,v13{1,1}(:,1),'b');plot(T14,v14{1,1}(:,1),'b');plot(T15,v15{1,1}(:,1),'b');plot(T16,v16{1,1}(:,1),'b');plot(T17,v17{1,1}(:,1),'b');plot(T18,v18{1,1}(:,1),'b');plot(T19,v19{1,1}(:,1),'b');plot(T20,v20{1,1}(:,1),'b');plot(T21,v21{1,1}(:,1),'b');plot(T22,v22{1,1}(:,1),'b');plot(T23,v23{1,1}(:,1),'b');plot(T24,v24{1,1}(:,1),'b');plot(T25,v25{1,1}(:,1),'b');plot(T26,v26{1,1}(:,1),'b');plot(T27,v27{1,1}(:,1),'b');plot(T28,v28{1,1}(:,1),'b');plot(T29,v29{1,1}(:,1),'b');plot(T30,v30{1,1}(:,1),'b');plot(T31,v31{1,1}(:,1),'b');plot(T32,v32{1,1}(:,1),'b');plot(T33,v33{1,1}(:,1),'b');plot(T34,v34{1,1}(:,1),'b');plot(T35,v35{1,1}(:,1),'b');plot(T36,v36{1,1}(:,1),'b');plot(T37,v37{1,1}(:,1),'b');plot(T38,v38{1,1}(:,1),'b');plot(T39,v39{1,1}(:,1),'b');plot(T40,v40{1,1}(:,1),'b');plot(T41,v41{1,1}(:,1),'b');plot(T42,v42{1,1}(:,1),'b');plot(T43,v43{1,1}(:,1),'b');plot(T44,v44{1,1}(:,1),'b');plot(T45,v45{1,1}(:,1),'b');plot(T46,v46{1,1}(:,1),'b');plot(T47,v47{1,1}(:,1),'b');plot(T48,v48{1,1}(:,1),'b');plot(T49,v49{1,1}(:,1),'b');plot(T50,v50{1,1}(:,1),'b');plot(T51,v51{1,1}(:,1),'b');plot(T52,v52{1,1}(:,1),'b');plot(T53,v53{1,1}(:,1),'b');plot(T54,v54{1,1}(:,1),'b');plot(T55,v55{1,1}(:,1),'b');plot(T56,v56{1,1}(:,1),'b');plot(T57,v57{1,1}(:,1),'b');plot(T58,v58{1,1}(:,1),'b');plot(T59,v59{1,1}(:,1),'b');plot(T60,v60{1,1}(:,1),'b');grid

De 1) si c'est juste pour effectuer le tracé et rien d'autre, il est absolument inutile de vouloir créer des variables pour stocker tout ceci.
De 2) je t'ai parlé de l'intérêt de créer ces atroces variables vi par rapport à utiliser directement sous_matrices(i,1) : c'est à dire aucun.

Au final tu peux partir sur quelque chose comme :
Code:

1 2 3 4 5 6 7 8 9 10 11 hold on for i = 1:60 % j'aurais tendance à remplacer 60 par numel(sous_matrices) debut = ... % sous_matrices{i,1}(:,4) (à voir pour obtenir un scalaire et non un vecteur) fin = ... % (size(sous_matrices{i,1},1)-1+(sous_matrices{i,1}(:,4)/TimeStep))*TimeStep (idem : scalaire et non un vecteur) T = debut : TimeStep : fin; V = sous_matrices{i,1}(:,1); plot(T, V, 'b') end title('Distribution des trajectoires de 60 paquets de véhicules sorties selon le modèle GSOM') xlabel('temps de simulation (s)') ylabel('position (m)')
Et si tu dois vraiment stocker ces T et V, fais-le dans deux uniques tableaux de cellules en remplaçant T par T{i] et V par V{i}.
Tu conviendras que c'est plus simple non ? Aucune prise de tête avec eval que je te conseille de vite oublier ;)

Merci infiniment pour tes précieuses remarques; j'apprends toujours des choses :rose2::merci:
Sinon, pour:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
clear all
close all
clc
load gsom.mat;
idxLignesZeros = find(gsom(:,1) == 0);
sous_matrices = mat2cell(gsom, diff([idxLignesZeros ; size(gsom,1)+1]), size(gsom,2));
TimeStep=0.1;
 
hold on
z=numel(sous_matrices);
for i = 1:z % j'aurais tendance à remplacer 60 par numel(sous_matrices)
    debut = sous_matrices{i,1}(:,4); % sous_matrices{i,1}(:,4) (à voir pour obtenir un scalaire et non un vecteur)
    fin = size(sous_matrices{i,1},1)-1+(sous_matrices{i,1}(:,4)/TimeStep)*TimeStep;
    T = debut : TimeStep : fin;
    V =  sous_matrices{i,1}(:,1);
    plot(T, V, 'b')
end
title('Distribution des trajectoires de paquets de véhicules sorties selon le modèle GSOM') 
xlabel('temps de simulation (s)')
ylabel('position (m)')

J'ai le message d'erreur suivant :

Code:

1
2
3
4
5
6
 
??? Error using ==> plot
Vectors must be the same lengths.
 
Error in ==> read_data at 16
    plot(T, V, 'b')

Citation:

Envoyé par Winjerome

Code:

1
2
    debut = ... % sous_matrices{i,1}(:,4) (à voir pour obtenir un scalaire et non un vecteur)
    fin = ... % (size(sous_matrices{i,1},1)-1+(sous_matrices{i,1}(:,4)/TimeStep))*TimeStep (idem : scalaire et non un vecteur)

Là, toi seule peut décider des bonnes valeurs à donner :)