Problème de boucle

**Raphh** · 29/11/2008, 10h50

Bonjour,

je suis débutant à matlab et je dois réaliser une simulation de diffusion atomique. Mon programme est primitif : des boucles for et if.

Le principe : je dispose d'une matrice 30x30 composée d'atomes 1 et 2 (rouge et bleu) et de 0: les lacunes. Je fais un balayage de la matrice pour déplacer les lacunes. Pour ne pas reprendre les mêmes lacunes je balaye un coup les indices pairs et un coup les indices impairs.
Je balaye les bords séparément.

Mon problème : le programme fonctionne bien pour les premières itérations, mais il se met à tourner dans le vide après un certain temps (assez court), les matrices ne bougeant plus... je ne comprends pas!

j'ai mis une grande boucle en t pour faire tourner les autres le nombre de fois que je veux.

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
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
 
mat(1:30,1:15)=1;
mat(1:30,16:30)=2;
for p=1:15
mat(floor(rand(1)*30),floor(12+6*rand(1)))=0;
end
for q=1:7
mat(floor(rand(1)*30),floor(1+10*rand(1)))=0;
mat(floor(rand(1)*30),floor(20+10*rand(1)))=0;
end
A=mat
 
 
 
hold on
for i=1:30 
    for j=1:30
        if A(i,j)==1
        plot(j/30,1-(i-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
         elseif A(i,j)==2
        plot(j/30,1- (i-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
    end
    end
end
 
for t=1:1000000
 
for m=2*(1:14)
    for n=2*(1:14)
        if A(m,n)==0
        a=floor(rand(1)*4)
              if a==0 && A(m-1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==0 && A(m-1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m+1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m+1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==2 && A(m,n+1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==2 && A(m,n+1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==3 && A(m,n-1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)             
              elseif a==3 && A(m,n-1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              end
                    if a==0 
                        A(m-1,n)=0
plot(n/30,1-(m-2)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==1 
                        A(m+1,n)=0 
plot(n/30,1-(m)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    elseif a==2 
                        A(m,n+1)=0
plot((n+1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    elseif a==3
                        A(m,n-1)=0
plot((n-1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    end
                    else continue 
         end
    end
end
 
 
 
for m1=2*(1:14)+1
    for n1=2*(1:14)+1
        if A(m1,n1)==0
        a=floor(rand(1)*4)
              if a==0 && A(m-1,n1)==1
                 A(m1,n1)=1
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==0 && A(m1-1,n1)==2
                 A(m1,n1)=2
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m1+1,n1)==1
                 A(m1,n1)=1
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m1+1,n1)==2
                 A(m1,n1)=2
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)             
              elseif a==2 && A(m1,n1+1)==1
                 A(m1,n1)=1
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==2 && A(m1,n1+1)==2
                 A(m1,n1)=2
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==3 && A(m1,n1-1)==1
                 A(m1,n1)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==3 && A(m1,n1-1)==2
                 A(m1,n1)=2
plot(n1/30,1-(m1-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              end
                    if a==0 
                        A(m1-1,n1)=0
plot(n1/30,1-(m1-2)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    elseif a==1 
                        A(m1+1,n1)=0
plot(n1/30,1-(m1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                  
                    elseif a==2 
                        A(m1,n1+1)=0
plot((n1+1)/30,1-(m1-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    elseif a==3
                        A(m1,n1-1)=0
plot((n1-1)/30,1-(m1-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    end
                    else continue 
         end
    end
end
 
 
 
for m=1 
    for n=2:29
        if A(m,n)==0
        a=floor(rand(1)*3)
 
              if a==0 && A(m+1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==0 && A(m+1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m,n+1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m,n+1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)            
              elseif a==2 && A(m,n-1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)            
              elseif a==2 && A(m,n-1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              end
 
                    if a==0 
                        A(m+1,n)=0
plot(n/30,1-(m)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==1 
                        A(m,n+1)=0
plot((n+1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==2
                        A(m,n-1)=0
plot((n-1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    end
                    else continue 
        end       
    end
end
 
 
 
 
for m=30 
    for n=2:29
        if A(m,n)==0
        a=floor(rand(1)*3)
              if a==0 && A(m-1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)
              elseif a==0 && A(m-1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m,n+1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)             
              elseif a==1 && A(m,n+1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)               
              elseif a==2 && A(m,n-1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)             
              elseif a==2 && A(m,n-1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              end
                    if a==0 
                        A(m-1,n)=0
plot(n/30,1-(m-2)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                     
                    elseif a==1 
                        A(m,n+1)=0
plot((n+1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==2
                        A(m,n-1)=0
plot((n-1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    end
                    else continue 
        end       
    end
end
 
 
 
for n=1 
    for m=2:29
        if A(m,n)==0
        a=floor(rand(1)*3)
              if a==0 && A(m-1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)               
              elseif a==0 && A(m-1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)               
              elseif a==1 && A(m+1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m+1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)             
              elseif a==2 && A(m,n+1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)            
              elseif a==2 && A(m,n+1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)             
              end
                    if a==0 
                        A(m-1,n)=0
plot(n/30,1-(m-2)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==1 
                        A(m+1,n)=0
plot(n/30,1-(m)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==2 
                        A(m,n+1)=0
plot((n+1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    end
                    else continue 
        end       
    end
end
 
 
 
 
for n=30 
    for m=2:29
        if A(m,n)==0
        a=floor(rand(1)*3)
              if a==0 && A(m-1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==0 && A(m-1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m+1,n)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
pause(0.05)              
              elseif a==1 && A(m+1,n)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)             
              elseif a==2 && A(m,n-1)==1
                 A(m,n)=1
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
 
 
pause(0.05)             
              elseif a==2 && A(m,n-1)==2
                 A(m,n)=2
plot(n/30,1-(m-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
pause(0.05)              
              end
                    if a==0 
                        A(m-1,n)=0
plot(n/30,1-(m)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==1 
                        A(m+1,n)=0
plot(n/30,1-(m)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                    
                    elseif a==2
                        A(m,n-1)=0
plot((n-1)/30,1-(m-1)/30,'o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8)
pause(0.05)                   
                    end
                    else continue 
        end       
    end
end
 
continue
end

En clair je ne comprends pas pourquoi les boucles de marchent plus après avoir été parcourues, mais vraiment pas!!

si quelqu'un a le courage de regarder, cela me serait d'une grande aide!
merci d'avance

**Caro-Line** · 29/11/2008, 11h32

Peut-on savoir quel est l'objectif final ?

Est-il nécessaire de visualiser (c'est joli mais j'ai planté mon MATLAB

) les modifications graphiquement ?

Ou veux-tu tout simplement avoir comme résultat une matrice où tous les zéros seraient en bout de ligne ?

Dans le dernier cas on peut trouver plusieurs solutions très rapides.
Exemple (utilisation de l'indexage logique ) :

Code :

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

1
2
3
4
5
6
%Boucle sur les lignes
for aLig=1:size(A,1)
    idx=(A(aLig,:)==0);%indices logiques des emplacements des 0
    A(aLig,:) = [A(aLig,~idx),A(aLig,idx)];%reagencement de la ligne
end
A

**Raphh** · 29/11/2008, 11h43

euh... les zéros vont se retrouver en bout de ligne?? c'est pas le but...
Le but est plutôt de montrer graphiquement que les atomes vont se mélanger. J'ai l'intention de faire une capture vidéo après...

Dans un premier temps je lui avais fait choisir les indices aléatoirement : avec randi mais il changeait toujours les mêmes atomes, du coup j'ai pensé au pair/impair.
La il tourne une fois puis, plus moyen de relancer les boucles...

**Caro-Line** · 29/11/2008, 11h50

Ok donc j'ai bien rien compris

Tu pourrais essayer d'expliquer algorithmiquement ce que tu veux car c'est pas évident de comprendre à partir du programme.
Le début ça va (enfin je ne suis sure de rien) : tu as une matrice avec des atomes rouges d'un côté, des bleus de l'autre, et dans tout ça il y a des trous que tu places aléatoirement.

C'est quoi la loi qui fait déplacer les lacunes ?

As-tu une idée d'à partir de quelle valeur de t tu estimes que plus rien ne bouge ?

**Raphh** · 29/11/2008, 11h59

ok, désolé. je vais détailler le principe pour faire bouger les lacunes.

D'abord je cherche les indices m et n tels que A(m,n)=0.
Puis je mets un rand pour que cette lacune soit remplacée par un des 4 atomes qui l'entourent.
Par exemple pour qu'elle soit remplacée par l'atome de la ligne précédente je vais mettre A(m-1)=0 (la lacune monte) et A(m,n)=1 si l'atome en A(m-1,n) était un atome1

Je fais ca pour les indices pairs puis impairs. Puis je précise l'opération sur les bords, par exemple sur la première ligne m=1, il faut que j'enlève le cas A(m-1,n).

Je viens de tester les boucles sans le t, le programme tourne une fois correctement en traitant les indices pairs, impairs, et les bords sans prendre deux fois le même atome : parfait.
Je refais ctrl+ entrée et la il ne change que les bords.
une fois que les lacunes ne sont plus au bord : crt+entrée ne fait plus d'effet.

voilà voilà...

**Caro-Line** · 29/11/2008, 12h32

Yep c'est normal (même si ce n'est pas toujours à t=4), au bout d'un moment toutes les lacunes se trouvent à des indices pair+impair ou impair+pair...

Tu peux le constater à la fin de ton programme en faisant :

Code :

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

1
2
[i,j]=find(A==0);
[i,j]

Les couples (i,j) correspondent aux indices des lacunes.

Donc avec ton programme tu mets de côté presque la moitié des emplacements.

Par contre si j'ai bien compris, à chaque t tu ne veux pas forcément bouger toutes les lacunes mais seulement certaines ?

Peut-être pourrais-tu prendre le truc différemment ie à chaque tour :
- Récupérer tous les indices tels que A = 0
- Prendre aléatoirement quelques indices parmi ceux-ci (s'inspirer de Comment obtenir un tirage aléatoire sur un nombre donné d'élément d'une matrice ? )
- Bouger ceux choisis

Je pense aussi qu'il y aurait moyen de faire plus ... simple, ou court surtout

Par exemple j'ai modifié le début de ton programme ainsi :

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
%"Capitalisation"
redmarker = struct('Marker','o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8);
bluemarker = struct('Marker','o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8);
whitemarker = struct('Marker','o','MarkerEdgeColor','w','MarkerFaceColor','w','MarkerSize',8);
 
%Matrice initiale
hold on
for i=1:30
    for j=1:30
        if A(i,j)==1
            plot(j/30,1-(i-1)/30,redmarker)
        elseif A(i,j)==2
            plot(j/30,1-(i-1)/30,bluemarker)
        end
    end
end
 
for t=1:1000000
    %les pairs "intérieurs"
    for m=2*(1:14)
        for n=2*(1:14)
            %A(m,n) est une lacune
            if A(m,n)==0
                %nombre alétoire entre 0 et 3
                a=floor(rand(1)*4)
                switch a
                    case 0
                        if A(m-1,n)==1
                            plot(n/30,1-(m-1)/30,redmarker)
                        elseif A(m-1,n)==2
                            plot(n/30,1-(m-1)/30,bluemarker)
                        end
                        pause(0.05)
                        A(m,n)=A(m-1,n);
                        A(m-1,n)=0
                        plot(n/30,1-(m-2)/30,whitemarker)
                        pause(0.05)
                    case 1
                        if A(m+1,n)==1
                            plot(n/30,1-(m-1)/30,redmarker)
                        elseif A(m+1,n)==2
                            plot(n/30,1-(m-1)/30,bluemarker)
                        end
                        pause(0.05)
                        A(m,n)=A(m+1,n);
                        A(m+1,n)=0
                        plot(n/30,1-(m)/30,whitemarker)
                        pause(0.05)
                    case 2
                        if A(m,n+1)==1
                            plot(n/30,1-(m-1)/30,redmarker)
                        elseif A(m,n+1)==2
                            plot(n/30,1-(m-1)/30,bluemarker)
                        end
                        pause(0.05)
                        A(m,n)=A(m,n+1);
                        A(m,n+1)=0
                        plot((n+1)/30,1-(m-1)/30,whitemarker)
                        pause(0.05)
                    case 3
                        if A(m,n-1)==1
                            plot(n/30,1-(m-1)/30,redmarker)
                        elseif A(m,n-1)==2
                            plot(n/30,1-(m-1)/30,bluemarker)
                        end
                        pause(0.05)
                        A(m,n)=A(m,n-1);
                        A(m,n-1)=0
                        plot((n-1)/30,1-(m-1)/30,whitemarker)
                        pause(0.05)
                end
            end
        end
    end

Nota : pour l'initialisation de la matrice, avec le rand, je ne dois pas avoir de chance mais 2 fois sur 3 le floor(rand) me donne un indice 0 et ça plante. Tu pourrais surement l'améliorer

**Raphh** · 29/11/2008, 13h58

MERCI!!!
le problème est réglé, et le code est bien simplifié! Ca fait plaisir!

Par contre, existe il une fonction qui permettrait de représenter par exemple le nombre d'atomes 1 suivant l'axe x? Pour visualiser la répartition, et la progression (pénétration) de ces atomes dans les autres.

**Caro-Line** · 29/11/2008, 14h11

Envoyé par Raphh

Par contre, existe il une fonction qui permettrait de représenter par exemple le nombre d'atomes 1 suivant l'axe x? Pour visualiser la répartition, et la progression (pénétration) de ces atomes dans les autres.

Peux-tu préciser ? (faire un dessin peut-être).
Tu peux modifier l'affichage des labels avec la propriété 'xticklabel' des axes.

Un truc du style (à améliorer c'est décalé), et surtout à remodifier à chaque modification :

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
%Matrice initiale
hold on
a=axes;
for i=1:30
    for j=1:30
        if A(i,j)==1
            plot(j/30,1-(i-1)/30,redmarker)
        elseif A(i,j)==2
            plot(j/30,1-(i-1)/30,bluemarker)
        end
    end
end
%Nombre de 1 dans chaque colonne
for j=1:30
    num(j)=sum(A(:,j)==1);
end
set(a,'xtick',0:1/29:1,'xticklabel',num)

**Jerome Briot** · 29/11/2008, 18h38

Code :

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

1
2
3
4
5
6
7
8
9
10
hold on
for i=1:30
    for j=1:30
        if A(i,j)==1
            plot(j/30,1-(i-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8)
        elseif A(i,j)==2
            plot(j/30,1- (i-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8)
        end
    end
end

peut être simplifié comme ceci :

Code :

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

1
2
3
4
5
[i,j] = find(A==1);
plot(j/30,1-(i-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8);
hold on
[i,j] = find(A==2);
plot(j/30,1-(i-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8);

Ou encore :

Code :

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

1
2
3
4
5
6
7
 
[j,i] = meshgrid(1:size(A,1));
idx = (A==1);
plot(j(idx)/30,1-(i(idx)-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','r','MarkerSize',8);
hold on
idx = (A==2);
plot(j(idx)/30,1-(i(idx)-1)/30,'o','MarkerEdgeColor','k','MarkerFaceColor','b','MarkerSize',8);

**Raphh** · 30/11/2008, 18h18

Bonjour,

en fait j'aimerais montrer comment évoluer la population d'atomes rouges par exemple, au début ils seront tous à gauche puis il vont s'étaler et se mélanger avec les bleus ce qui va donner ces allures de courbes :

j'espère etre clair...