affichage progressif

Version imprimable

Bonjour,
j'ai fabriqué le code suivant sous Scilab et qui affiche un triangle de Sierpinski dans un triangle équilatéral (partant de (0;0), on prend de manière répétitive le milieu du point précédent et d'un des sommets A, B ou C choisis au hasard). Le tout est présenté sous forme matricielle.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 nbpoints=10000 P=zeros(2,nbpoints); for i=2:nbpoints tirage=rand(); if tirage<1/3 then A=[[0.5,0];[0,0.5]];V=[0;0]; else if tirage<2/3 then A=[[0.5,0];[0,0.5]];V=[1/2;0]; else A=[[0.5,0];[0,0.5]];V=[1/4;sqrt(3)/4]; end end P(:,i)=A*P(:,i-1)+V; plot(P(1,i),P(2,i),"or") end
Obtenir un affichage de la figure une fois terminée est facile: il suffit de sortir l'avant-dernière ligne de la boucle et on obtient sans problème la figure attendue.
Le fait de placer cette ligne à cet endroit devrait me permettre d'observer un affichage progressif (point après point) de la figure. Or, cet affichage n'est pas lisible car on a un clignotement incessant très désagréable de l'image.
Quelqu'un connaît-il la solution à ce problème qui dépasse d'ailleurs le cadre stricte de cet exemple particulier, c'est à dire obtenir un affichage progressif lisible d'une construction.
Pour info, je travaille avec la version 5.4.0 de Scilab sur un Macbook Pro.
Merci d'avance pour vos lumières.

Essaie ceci :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
nbpoints=10000;
P=zeros(2,nbpoints);
 
f = figure();
f.pixmap = "on";
f.background = color("white");
 
for i=2:nbpoints 
    tirage=rand();
    if tirage<1/3 then
        A=[0.5,0;0,0.5];V=[0;0];
    else
        if tirage<2/3 then 
            A=[0.5,0;0,0.5];V=[1/2;0]; 
        else 
            A=[0.5,0;0,0.5];V=[1/4;sqrt(3)/4];
        end
    end    
    P(:,i)=A*P(:,i-1)+V;
    plot(P(1,i),P(2,i),"or")
    drawnow()
end

12/03/2013, 11h30
arkhimedes

Bonjour Dut,
tout d'abord merci pour cette réponse très rapide.
J'ai testé votre proposition: cela semble réduire le clignotement qui malheureusement reste pourtant bien présent.
Je profite de cette réponse pour demander un complément d'info: je n'utilise pas Scilab depuis très longtemps mais je me demande s'il est possible de faire apparaître ce genre de figure progressivement à l'aide par exemple d'un curseur (style ce qu'on peut avoir facilement dans un logiciel tel que GeoGebra), c'est à dire interagir de manière un peu dynamique avec cette figure ?
Scilab est un excellent programme mais malgré de nombreuses recherches sur internet, j'ai l'impression que cette question d'affichage progressif ou d'interactivité avec la figure ne se résoud pas facilement...
Merci en tout cas de vous intéresser à mon problème :D

Citation:

Envoyé par arkhimedes

J'ai testé votre proposition: cela semble réduire le clignotement qui malheureusement reste pourtant bien présent.

Voici une autre technique :
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 nbpoints=10000; P=zeros(2,nbpoints); f = figure(); f.pixmap = "on"; f.background = color("white"); plot(%nan,%nan,"ro"); a = gca(); a.tight_limits = "on"; a.auto_scale = "off"; a.data_bounds = [0,0,1,1]; for i=2:10000//nbpoints tirage=rand(); if tirage<1/3 then A=[0.5,0;0,0.5];V=[0;0]; else if tirage<2/3 then A=[0.5,0;0,0.5];V=[1/2;0]; else A=[0.5,0;0,0.5];V=[1/4;sqrt(3)/4]; end end P(:,i)=A*P(:,i-1)+V; a.children.children.data = P(:,1:i).'; drawnow() end
Par contre, je n'arrive pas à fixer les limites des axes...

Citation:

Envoyé par arkhimedes

Je profite de cette réponse pour demander un complément d'info: je n'utilise pas Scilab depuis très longtemps mais je me demande s'il est possible de faire apparaître ce genre de figure progressivement à l'aide par exemple d'un curseur (style ce qu'on peut avoir facilement dans un logiciel tel que GeoGebra), c'est à dire interagir de manière un peu dynamique avec cette figure ?
Scilab est un excellent programme mais malgré de nombreuses recherches sur internet, j'ai l'impression que cette question d'affichage progressif ou d'interactivité avec la figure ne se résoud pas facilement...

Le plus simple serait de démarrer une nouvelle discussion dans le forum Scilab en nous donnant plus d'informations sur ce que tu cherches à faire exactement.